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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

139
fi:=-Tn[i]/tau;
{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;
{определяем значение температуры на правой границе на основе
граничного условия, используя соотношение (21)}
T[N]:=(2.0*a2*tau*beta[N-1]+sqr(h)*Tn[N])/(2.0*a2*tau
*(1.0-alfa[N-1])+sqr(h));
{используя соотношение (7) определяем неизвестное поле температуры
в талой зоне грунта}
for i:= N-1 downto k 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; {поле температуры на данном временн
ом слое
определили}
{увеличиваем переменную времени на шаг τ
n+1
}
time:=time+tau;
{формируем файл, содержащий шаги по времени τ
n+1
}
Writeln(f1,'tau',(k-1),' = ',tau:6:4);
end;
{выводим результат в файл}
Assign(f,'res.txt');
Rewrite(f);
Writeln(f,'Глубина грунта L = ',L:6:4);
Writeln(f,'Число узлов по координате N = ',N);
Writeln(f,'Коэффициент теплопроводности промерзшей зоны грунта
lamda1 = ',lamda1:6:4);
Writeln(f,'Плотность промерзшей зоны грунта ro1 = ',ro1:6:4);
Writeln(f,'Теплоемкость промерзшей зоны грунта с1 = ',c1:6:4);
Writeln(f,'Коэффициент теплопроводности талой зоны грунта lamda2 =
',lamda2:6:4);
Writeln(f,'Плотность талой зоны грунта ro2 = ',ro2:6:4);
Writeln(f,'Теплоемкость талой зоны грунта с2 = ',c2:6:4);
Writeln(f,'Начальная температура T0 = ',T0:6:4);
Writeln(f,'Температура на
границе х = 0'
,Tc:6:4);
Writeln(f,'Температура фазового перехода',Tfr:6:4);

Страницы