Аппроксимация зависимости

 

Q0= f (t0; tк) и Ne= f (t0; tк)

для агрегата  1 ступень S3-900

                       2 ступень S3-315

Произведем расчет для tк=20°C

8.2.1 Исходное уравнение для Q0= f (t0; tк)

Q0/Vh=exp*(A11tк)*(t0+90)(a+d1tк)                                            (8.8)

где: Vh=VhS3-900+VhS3-315;

    VhS3-900= 792 м3/ч;

VhS3-315= 317 м3/ч;

Vh=792+317=1109 м3

Для точек:

t0=-40°С;          Q0=170,26;

t0=-55°С;          Q0=79,477;

Записываем исходное уравнение в виде:

Q0/Vh=а+(t0+90)b                                                   (8.9)

где: а=(A11tк);

    b=(a+d1tк)

логарифмируя обе стороны получаем:

ln(170,26/1109)=ln a+b*ln (-40+90)

ln(79,477/1109)=ln a+b*ln (-55+90)

решаем систему уравнений

ì_-1,873888= l n a+b*3,912023

î - 2,635747= ln a+b*3,555349

    0,761859=b*0,356674

откуда b=0,761859/0,356674=2,136;

Подставляем значение в любое уравнение получаем а:

ln a=-10.230005

a=0.000036

При подстановке коэффициентов в уравнение получаем:

Q0/Vh=3,607*10-5(t0+90)2,136                                   (8.10)

Рассчитываем при tк=30°С

Для точек:

t0=-40°С;          Q0=164,41;

t0=-55°С;          Q0=76,266;

Записываем исходное уравнение в виде:

Q0/Vh=а+(t0+90)b                                                       (8.11)

логарифмируя уравнение получаем:

ln (Q0/Vh)= ln a+b*ln (t0+90)

Подставляем значения и вычисляем:

ln(164,41/1109)=ln a+b*ln (-40+90)

ln(76,266/1109)=ln a+b*ln (-55+90)

решаем систему уравнений

ì_-1,908852= ln a+b*3,912023

î - 2,676988= ln a+b*3,555349

0,768136=b*0,356674

откуда b=0,768136/0,356674=2,153608;

Подставляем значение в любое уравнение получаем а:

ln a=-10.109749

a=4.068*10-5

При подстановке получаем исходное уравнение:

Q0/Vh=4.068*10-5/(t0+90)2,0936;                                                  (8.12)

Аналогично проводим расчет для других температур tк и получаем значения коэффициентов а и b:

1) tк=20°С; а=3,607*10-5; b=2,136;

2) tк=25°С; а=3,424*10-5; b=2,144;

3) tк=30°С; а=3,252*10-5; b=2,154;

4) tк=35°С; а=3,092*10-5; b=2,162;

5) tк=45°С; а=2,803*10-5; b=2,179;

Произведем расчет коэффициентов с1 и d1  в уравнении:

b=c1+d1tк                                                                 (8.13)

Для значений tк:

tк=20°С;           b=2,136;

tк=25°С;           b=2,144;

подставляем значения в уравнение и вычисляем:

ì_2,136= a+d1 *20

î 2,179= a+d 1 *45

 25d1 =0,0427

Откуда d1=0,0427/25=17,07*10-4;

Подставляя в первое уравнение значение d1 получаем с1:

с1=2,10184;

Исходное уравнение будет иметь вид:

b=2,1+17,08*10-4tк

Произведем расчет коэффициентов А1 и В1 в уравнении:

а=exp(A11tк)                                                         (8.14)

Для значений tк:

 tк=20°С; а=3,607*10-5;

 tк=45°С; а=2,803*10-5;

Логарифмируя получаем:

ln a= A11tк

Подставляя значения решаем систему уравнений:

ìln 3,706*10-5 = A11 *20

î ln 2,803*10-5 = A11 *45

ì_-10,23005 = A11 *20

î - 10,482334= A11 *45

  -0.252329= В1 *25

Откуда В1= -0.252329/25= -1,00931*10-2, подставляя значения В1 в уравнение получаем:

А1=-10,028143

Исходное уравнение будет иметь вид:

а=ехр (-10,028-1,00931*10-2*tк)

Получаем значения коэффициентов:

А1=-10, 028143;      с1=2,102;

В1=-1,00931*10-2;             d1=17,08*10-4

 

Проверка:

Подставляем в первоначальное уравнение:

Q0/Vh=exp*(A11tк)*(t0+90)(a+d1tк)

значение коэффициентов и значения

tк=45°С; t0= - 45°С

получаем:

Q=1109*exp*(10,028-1,00931*10-2*45)*(-45+90)(2,102+17,08*0,0001*45)=124,36

Значение Q0=123,93 при tк=45°С; t0= - 45°С.

По вычисленным значениям коэффициентов а и b строим график рис. 8.6.

 

8.2.2 Исходное уравнение для Q0= f (t0; tк)

Ne/Vh=(a2tк+b2)*t0+(c2tк+d2)                                               (8.16)

Vh=1109 м3

Произведем расчет для tк=20°С;

Для точек:

t0=-40°С;          Ne =91,646;

t0=-55°С;          Ne =78,456;

Преобразуем исходное уравнение:

Ne =аt0+b                                                                    (8.17)

где: а=(A2tк2);

    b=(с2tк +d2)

подставляем значения и вычисляем

ì_91,646= a*-40+b

î 78,456= a*-55+b

13,190=а*15

откуда а=0,879333

91,646=0,879333*-40+b

b=126,81933

Подставляем коэффициент в уравнение:

Ne =0,879333t0+126,81933                                             (8.18)

Аналогично проводим расчет для tк = (25; 30; 35; 40; 45)°С и получаем значения коэффициентов а и b:

1) tк=20°С; а=0,879333;   b=126,81933;

2) tк=25°С; а=0,894333;   b=136,71333;

3) tк=30°С; а=0,9090666; b=146,5566;

4) tк=35°С; а=0,9233333; b=156,32333;

5) tк=40°С; а=0,9373333; b=166,05333;

5) tк=45°С; а=0,952;         b=175,77;

Произведем расчет коэффициентов А2 и В2  в уравнении:

а=(A2tк2)                                                                      (8.19)

Для значений tк:

tк=20°С;           а=0,879333;

tк=45°С;           а=0,952;

подставляем значения в уравнение и решаем систему:

ì_ 126,81933= с2 *20+d2

î 175,77= c 2 *45+ d 2

48,95067= с2 *25

Откуда с2=48,95067/25=1,9580268;

тогда: 126,81933= 1,9580268*20+d2

d2=87,6588;

Исходное уравнение будет иметь вид:

b=1,9580268tк+87,6588                                                    (8.20)

Значения коэффициентов:

А2=0,0029066;        с2=1,9580268;

В2=0,821203;           d2=87,6588;

Исходное уравнение при подстановке и вычислении:

Ne=(0,0029066tк +0,821203)t0+(1,9580268+87,6588)

Ne/Vh=(0,262*10-5tк +0,74*10-3)t0+(1,765510-3tк+0,079)

 

Проверка: при t0= -50°С; tк=35°С

Подставляем в первоначальное уравнение:

Ne/Vh=(0,262*10-5*35+0,74*10-3)*-50+(1,765510-3*35+0,079)=110,02 кВт

Neр= Ne=110,02 кВт

По полученным результатам строим график рис. 8.7.

 


 

 

 

 


9. Таблица рабочих режимов СХУ

 







Таблица рабочих режимов СХУ

 

Таблица 9.1

 

Параметры

Температура забортной воды tw °C

10 15 20 30

Трюм

Давление МПа        
- кипения 0,08455 0,04372 -0,01382 -0,00729
- конденсации 10,5481 10,9461 12,4733 13,5218
Температура кипения °C -36,779 -37,643 -38,988 -40,036
Температура на входе в КМ, °C -16,852 -17,735 -18,982 -20,033
Температура нагнетания, °C 82,190 82,988 84,322 85,799
Температура перед РВ испарительной системы, °C 30 30 30 30
Температура охл. воздуха, °C -28 -28 -28 -28
Сила тока эл. дв.КМ 63,119 64,528 67,733 71,329
Коэффициент регулирования 0,188 0,207 0,253 0,321

Морозильный агрегат

Давление МПа        
- кипения -0,4587 -0,4517 -0,4332 -0,4102
- конденсации 10,5481 10,9461 12,4733 13,5218
Температура кипения °C -50,5481 -50,006 -49,508 -48,736
Температура нагнетания, °C 82,190 82,988 84,322 85,799
Температура перед РВ испарительной системы, °C 30 30 30 30
Производительность МК 67,908 65,015 63,141 61,136
Сила тока эл. дв.КМ 109,314 109,679 110,251 111,0348

 


 


10. Выводы и рекомендации

 








Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  



double arrow
Сейчас читают про: