Санкт-Петербургский Государственный Университет Низкотемпературных и
Пищевых Технологий
Лабораторная работа № 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
Вывод: