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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

92
function lamda(x:real):real;
{функция вычисления коэффициента теплопроводности по формуле
(31)}
begin
lamda:=5500/(560+x)+0.942*(1e-10)*x*sqr(x);
end;
begin
clrscr;
{с клавиатуры вводим все необходимые входные параметры}
Writeln('Введите количество пространственных узлов, N');
Readln(N);
Writeln('Введите окончание по времени, t_end');
Readln(t_end);
Writeln('Введите толщину пластины, L');
Readln(L);
Writeln('Введите плотность материала пластины, ro');
Readln(ro);
Writeln('Введите теплоемкость материала пластины, c');
Readln(c);
Writeln('Введите начальную температуру в К, T0');
Readln(T0);
Writeln('Введите температуру в К
на границе х=0, Th
');
Readln(Th);
Writeln('Введите температуру в К на границе х=L, Tc');
Readln(Tc);
{определяем расчетный шаг сетки по пространственной координате}
h:=L/(N-1);
{определяем расчетный шаг сетки по времени}
tau:=t_end/100.0;
{определяем поле температуры в начальный момент времени}
for i:= 1 to N do
T[i]:=T0;
{проводим интегрирование нестационарного уравнения
теплопроводности}
time:=0;
while time<t_end do {используем цикл с предусловием}
begin
{увеличиваем переменную времени на шаг τ}
time:=time+tau;

Страницы