Санкт-Петербургский Государственный Университет Низкотемпературных и 

Пищевых Технологий













Лабораторная работа № 4


«Разработка математической модели холодильной установки при переменной температуре конденсации»






















Выполнил: Дорогин В.А.

Проверил: Круглов А. А.











2004 г.

Цель: по исходным данным и графику предшествующей математической модели разработать математическую модель холодильной установки при условии, что температура конденсации переменна. В дополнение к подсистеме уравнений, составленной раньше, необходимо добавить соотношения:

уравнение теплового баланса.

уравнение теплового баланса холодильной установки с учётом статической модели компрессора.

уравнение теплопередачи конденсатора.

уравнение оценивающее коэффициент эффективности градирни.


1.  EMBED Equation.3  ;  где Nн – количество насосов

2.  EMBED Equation.3  

3.  EMBED Equation.3  ;

где:  EMBED Equation.3  

4.  EMBED Equation.3   ;  где: tМТ – температура мокрого термометра


Искомые величины: 

X[1] = t0 X[5] = tw2

X[2] = tв2 X[6] = QКД

X[3] = b X[7] = tk

X[4] = tw1


Текст программы:

PROGRAM IK3;

USES SKMATH2;

const tpm=1.0;toc=31.0;fn=9300;kn=0.00030;fvo1=230;kvo=0.011;

      vvo1=4.694;nvo1=24.9;nuvo=0.32;fkd=40;kkd=0.660;vbn1=0.0022;

      tmt=25.3;nugr=0.32;nubn=0.37;nbo1=24.9;nbn1=1.5;rw=2421;nvo=9;

      cv=1;rv=1.3;cw=4.19;nvn=9;

      pw=1000;tcg=10.5;t=7550;

var Qtp,Q0:real;

    x:massiv;

procedure sistnelur(x:massiv;var f:massiv);

begin

f[1]:=cv*rv*nvo*vvo1*(tpm-x[2])-Qtp;

f[2]:=kvo*fvo1*nvo*((tpm-x[2])/ln((tpm-x[1])/(x[2]-x[1])))-Qtp;

f[3]:=(722.9075+24.7827*x[1]+0.1880*x[1]*x[1]-5.4947*x[4]-0.0917*x[1]*x[4])*x[3]-Qtp;

f[4]:=(38.4166-1.7149*x[1]-0.0405*x[4]+2.0972*x[1]*x[1]+0.0672*x[1]*x[4])*x[3]+Qtp-x[6];

f[5]:=kkd*fkd*((x[5]-x[7])/ln((x[4]-x[7])/(x[4]-x[5])))-x[6];

f[6]:=cw*(x[5]-x[7])*vbn1*nvn*pw-x[6];

f[7]:=(x[5]-x[7])/(x[5]-tmt)-nugr;

end;

begin

Qtp:=kn*fn*(toc-tpm)+nvo*Nvo1*Nuvo;

nelur(7,sistnelur,x);

q0:=722.9075+24.7827*x[1]+0.1880*x[1]*x[1]-5.4947*x[4]-0.0917*x[1]*x[4];

writeln('Qtp=',Qtp:6:2, '         Q0=',Q0:6:2);

readln;

end.

Полученные результаты:

N X[1] X[2] X[3] X[4] X[5] X[6] X[7]

0 -9,0000 -1,0000 0,7000 24,7000 27,7000 85,0000 40,0000

1 -7,3170 -1,8298 0,3971 32,6198 27,4030 189,8137 25,1111

2 -7,3377 -1,8298 0,4113 36,0675 29,1629 214,9394 26,5696

3 -7,3377 -1,8298 0,4229 38,1551 31,2016 216,1373 28,5937

4 -7,3377 -1,8298 0,4326 39,8349 32,8495 217,1300 30,2295

5 -7,3377 -1,8298 0,4363 40,4403 33,4426 217,5143 30,8688

6 -7,3377 -1,8298 0,4366 40,4930 33,4942 217,5494 30,8688


Найденные корни системы:

X[1] = – 7,3377

X[2] = – 1,8298

X[3] = 0,4366

X[4] = 40,4934

X[5] = 33,4946

X[6] = 217,5497

X[7] = 30,8693

Qtr = 155,412

Qкм = 355,929


Вывод: