Приложение Г. Программа расчета смесителя

Программа расчета смесителя

Обозначения, принятые в программе:

N – число компонентов в потоке;

C1[1..N] – содержание компонентов в 1-м входном потоке, масс. доли;

C2[1..N] – содержание компонентов во 2-м входном потоке, масс. доли;

С3[1..N] – содержание компонентов в выходном потоке, масс. доли;

G1,G2 – расходы входных потоков, кг/ч;

G – расход выходного потока, кг/ч;

CP1,CP2 – удельные теплоемкости входных потоков;

CP – удельная теплоемкость выходного потока;

a[1..N], b[1..N], c[1..N], d1..N] – коэффициенты для расчета значений теплоемкости компонентов потока в зависимости от температуры;

Т1,Т2 – температуры входных потоков;

Т – значение температуры выходного потока;

Тx – значение температуры выходного потока на предыдущей итерации;

Е – заданная точность вычисления температуры выходного потока.

Program smesit;

var

g,cp,cp1,cp2,e,g1,g2,t1,t2,t,tx:real;

i,k:integer;

c1,c2,c3,c4,a,b,c,d:array[1..4] of real;

f1,f2:text;

begin

assign (f1,'smesit.dta');

assign(f2,'smesit.rez');

reset(f1);rewrite(f2);

for i:=1 to 4 do read(f1,a[i]);

for i:=1 to 4 do read(f1,b[i]);

for i:=1 to 4 do read(f1,c[i]);

for i:=1 to 4 do read(f1,c1[i]);

for i:=1 to 4 do read(f1,c2[i]);

read(f1,e,g1,g2,t1,t2);

g:=g1+g2;

for i:=1 to 4 do c3[i]:=(g1*c1[i]+g2*c2[i])/g;

cp1:=0; cp2:=0;

for i:=1 to 4 do begin

cp1:=cp1+(a[i]+b[i]*t1+c[i]*t1*t1+d[i]*t1*t1*t1)*c1[i];

cp2:=cp2+(a[i]+b[i]*t2+c[i]*t2*t2+d[i]*t2*t2*t2)*c2[i];

end;

t:=(t1+t2)/2;

repeat

tx:=t;cp:=0;

for k:=1 to 4 do

cp:=cp+(a[k]+b[k]*t+c[k]*t*t+ d[k]*t*t*t)*c3[k];

t:=(g1*cp1*t1+g2*cp2*t2)/g/cp;

until abs(tx-t)<e;

writeln(f2,t:8:3,cp:8:3,cp1:8:3,cp2:8:3,g:8:3);

writeln(f2,c3[1]:8:3,c3[2]:8:3,c3[3]:8:3,c3[4]:8:3);

close(f2);

end.

Приложение Д