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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

127
bi:=2.0*lamda/sqr(h)+ro*c/tau;
ci:=lamda/sqr(h);
fi:=-ro*c*Tn[i]/tau-q*k0*ro*exp(-E/(R*T[i]));
{alfa[i], beta[i] – прогоночные коэффициенты}
alfa[i]:=ai/(bi-ci*alfa[i-1]);
beta[i]:=(ci*beta[i-1]-fi)/(bi-ci*alfa[i-1]);
end;
{определяем значение температуры на правой границе на основе
правого граничного условия, используя соотношение (25)}
T[N]:=(ro*c*sqr(h)*Tn[N]+2.0*tau*(lamda*beta[N-1]+kapa*h*Te))
/(ro*c*sqr(h)+2.0*tau*(lamda*(1-alfa[N-1])+kapa*h));
{используя соотношение (7) определяем неизвестное поле
температуры}
for i:= N-1 downto 1 do
T[i]:=alfa[i]*T[i+1]+beta[i];
{определяем максимум модуля разности температур на данной и
предыдущей итерации}
max:=abs(T[1]-Ts[1]);
for i:=2 to N do
if max < abs(T[i]-Ts[i]) then max:=abs(T[i]-Ts[i]);
until max<=eps; {поле температуры на данном временном слое
определили}
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,'Коэффициент теплообмена kapa = ',kapa:6:4);
Writeln(f,'Температура внешней среды Te = ',Te:6:4);
Writeln(f,'Тепловой эффект химической реакции q = ',q:6:4);
Writeln(f,'Предэкспонент k0 = ',k0:6:4);
Writeln(f,'Энергия активации химической реакц
ии E = ',E:6:4);
Writeln(f,'Результат получен с шагом по координате h = ',h:6:4);
Writeln(f,'Результат получен с шагом по времени tau = ',tau:6:4);
Writeln(f,'Температурное поле в момент времени t = ',t_end:6:4);

Страницы