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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

72
{используя соотношение (7) определяем неизвестное поле температуры
на промежуточном временном слое}
for j:= Ny-1 downto 1 do
T[i,j]:=alfa[j]*T[i,j+1]+beta[j];
end;
end;
{цикл с предусловием окончен}
{выводим результат в файл}
Assign(f,'res.txt');
Rewrite(f);
Writeln(f,'Длина пластины L = ',L:6:4);
Writeln(f,'Толщина пластины H = ',H:6:4);
Writeln(f,'Число узлов по пространственной координате x в пластине
Nx = ',Nx);
Writeln(f,'Число узлов по пространственной координате y в пластине
Ny = ',Ny);
Writeln(f,'Коэффициент теплопроводности материала пластины lamda =
',lamda:6:4);
Writeln(f,'Плотность материала пластины ro = ',ro:6:4);
Writeln(f,'Теплоемкость материала пластины с = ',c:6:4);
Writeln(f,'Начальная температура T0 = ',T0:6:4);
Writeln(f,'Температура на границе x = 0 области решения Th = ',Th:6:4);
Writeln(f,'Температура на границе x = 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);
Assign(g,'tempr.txt');
Rewrite(g);
for i:=1 to Nx do
for j:=1 to Ny do
writeln(g,' ',hx*(i-1):10:8,' ',hy*(j-1):10:8,' ',T[i,j]:8:5);
close(g);
end.

Страницы