Разностные методы решения задач теплопроводности. Кузнецов Г.В - 87 стр.

UptoLike

ВУЗ: 

ТПУ | Томск

Составители: 

Рубрика: 

87
Writeln(f,'Начальная температура T0 = ',T0:6:4);
Writeln(f,'Температура на границе х=0 области решения Th = ',Th:6:4);
Writeln(f,'Температура на границе х=L области решения Tc = ',Tc:6:4);
Writeln(f,'Результат получен с шагом по координате x hx = ',hx:6:4);
Writeln(f,'Результат получен с шагом по координате y hy = ',hy:6:4);
Writeln(f,'Результат получен с шагом по времени tau = ',tau:6:4);
Writeln(f,'Температурное поле в момент времени t = ',t_end:6:4);
close(f);
{определяем коэффициент температуропроводности}
a1:=lamda1/(ro1*c1);
a2:=lamda2/(ro2*c2);
a3:=lamda3/(ro3*c3);
{определяем поле температуры в начальный момент вр
емени}
for i:= 1 to Nx do
for j:= 1 to Ny do
T[i,j]:=T0;
{проводим интегрирование нестационарного уравнения
теплопроводности}
time:=0;
while time<t_end do {используем цикл с предусловием}
begin
{увеличиваем переменную времени на шаг τ}
time:=time+tau;
{ СЛАУ в направлении оси Ох
в промежутке между границей y=0 и включением 3}
for j:=1 to ny1 do
progonx(j,lamda1,ro1,c1,T);
{на нижней границе включения 3}
progonxIV(ny1+1,nx1+nx2+nx3,nx4,lamda1,0.5*(lamda1+lamda3),ro1,0.5*(
ro1+ro3),c1,0.5*(c1+c3),T);
{во включении 3}
for j:=ny1+2 to ny1+ny2 do
progonxIV(j,nx1+nx2+nx3,nx4,lamda1,lamda3,ro1,ro3,c1,c3,T);
{на верхней границе вкл
ючения 3}
progonxIV(ny1+ny2+1,nx1+nx2+nx3,nx4,lamda1,0.5*(lamda1+lamda3),ro1,
0.5*(ro1+ro3),c1,0.5*(c1+c3),T);
{между включением 3 и включением 2}
for j:=ny1+ny2+2 to ny1+ny2+ny3 do
progonx(j,lamda1,ro1,c1,T);

Страницы