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

UptoLike

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

Рубрика:

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,'Результат получен с шагом по координате h = ',h: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 N do

writeln(g,' ',h*(i-1):6:3,' ',T[i]:8:5);

close(g);

end.

Результаты расчетов по приведенной программе при L = 0,1 м,

λ = 46 Вт/(м⋅ºC), ρ = 7800 кг/м

, с = 460 Дж/(кг⋅ºC), Т

=20 °С, Т

= 300 °С,

= 100 °С через 60 секунд процесса нагрева приведены на рис. 4.

Рис. 4. Распределение температуры по толщине пластины в момент

времени t = 60 с

Как отмечалось выше, рассмотренная расчетная схема является

неявной, т.е. для определения поля температуры приходится решать

систему линейных алгебраических уравнений. Но помимо