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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

22
{определяем расчетный шаг сетки по времени, руководствуясь
условием устойчивости}
tau:=0.25*sqr(h)/a;
{определяем поле температуры в начальный момент времени}
for i:= 2 to N-1 do
T[i]:=T0;
{определяем значения температуры на границе}
T[1]:=Tl;
T[N]:=Tr;
{проводим интегрирование нестационарного уравнения
теплопроводности}
time:=0;
while time<t_end do {используем цикл с предусловием}
begin
{увеличиваем переменную времени на шаг}
time:=time+tau;
{запоминаем поле температуры на предыдущем слое по времени}
for i:= 1 to N do
TT
[i]:=T[i];
{определяем неизвестное поле температуры по соотношениям (12)}
for i:= 2 to N-1 do
T[i]:=TT[i]+a*tau/sqr(h)*(TT[i+1]-2.0*TT[i]+TT[i-1]);
end;
{выводим результат в файл}
Assign(f,'res.txt');
Rewrite(f);
Writeln(f,'Толщина пластины L = ',L:6:4);
Writeln(f,'Число узлов по пространственной координате N = ',N);
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, Tl = ',Tl:6:4);
Writeln(f,'Температура на границе x = L, Tr = ',Tr:6:4);
Writeln(f,'Результат получен с
шагом по пространственно
й координате
x, h = ',h:6:4);
Writeln(f,'Результат получен с шагом по времени tau = ',tau:6:4);
Writeln(f,'Температурное поле в момент времени t = ',t_end:6:4);
close(f);
Assign(g,'tempr.txt');

Страницы