Лекция примеры моделирования

Q_р - расход тепла, равный Q_р =  ^. G₂^. C ^. T_2.
Подставляя в уравнение теплового баланса получим
Это дифференциальное уравнение решается с учетом того, что и объем V и температура Т₂ изменяются во времени, т.е.
Решая относительно температуры вещества Т₂, получим:
(3.4)

Объединив уравнения (3.1-3.4) в единую систему уравнений, получим математическую модель для процесса, протекающего в гидравлической емкости с паровой рубашкой.

Для решения дифференциальных уравнений в данной системе уравнений можно использовать метод Эйлера. Определив последовательность решений, составляем блок-схему решения задачи.


БЛОК-СХЕМЫ РЕШЕНИЯ ЗАДАЧИ ИМЕЮТСЛЕДУЮЩИЙ ВИД (см.рис.19. и рис.20.).
ПРОГРАММА РЕШЕНИЯ ЗАДАЧИ НА ЯЗЫКЕ TURBOPASCAL
program isitgich;
var g2, a1,a2,f1,f2,tk,td,t2,vd,c1,c2,dt,t1,t,g1,k,r1,s,h,rd,r,g,p:real;
i:integer;
begin
a1:=10010;
a2:=3003;
f1:=21;
tk:=132;
td:=103;
t2:=22;
vd:=0.04;
c1:=4200;
c2:=5010;
dt:=0.5;
t1:=24;
t:=0;
g1:=0.03;
k:=0.004;
r1:=1.7;
rd:=8010;
h:=3;
r:=901;
g:=11;
for i:=1 to 1000 do begin
s:=3.14*r1*r1;
f2:=2*3.14*r1*r1* h+s;
p:=r*g*h;
g2:=k*sqrt(r*g*h);
h:=h+((g1-g2)/s)*dt;
td:=td+(((a1*f1*(tk-td)-(a2*f2)*(td-t2))/(vd*rd*c1)))*dt;
t2:=t2+(((g1*t1)/(s*h))-((g2*t2)/(s*h))+((a2*f2)*(td-t2)/(c2*r*s*h))-(t2*(g1-g2)/(s*h)))*dt;
t:=t+dt;
writeln('t=' ,t:6:3,'t2=' ,t2:6:3,'td=' ,td:6:3,'g2=' ,g2:6:3,'h=' ,h:6:3);
end;
end.

Нет
Да

19-расм.

нет
да

нет
да