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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

29
T[i]:=T0;
{проводим интегрирование нестационарного уравнения
теплопроводности}
time:=0;
while time<t_end do
{используем цикл с предусловием}
begin
{увеличиваем переменную времени на шаг τ}
time:=time+tau;
{определяем начальные прогоночные коэффициенты на основе левого
граничного условия}
alfa[1]:=1.0;
beta[1]:=0.0;
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8)}
for i:= 2 to N-1 do
begin
{ai, bi, ci, fi – коэффициенты канонического представления СЛАУ с
трехдиагональной матрицей}
ai:=0.5*lamda*(2*i-1)/(sqr(h)*(i-1));
ci:=0.5*lamda*(2*i-3)/(sqr(h)*(i-1));
bi:=ai+ci+ro*c/tau;
fi:=-ro*c*T[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;
{определяем значение температуры на правой границе}
T[N]:=Th;
{используя соотношение (7) определяем неизвестное поле
температуры}
for i:= N-1 downto 1 do
T[i]:=alfa[i]*T[i+1]+beta[i];
end;
{выводим результат в файл}
Assign(f,'res.txt');
Rewrite(f);
Writeln(f,'Радиус цилиндра R = ',R:6:4);
Writeln(f,'Число узлов по координате N = ',N);
Writeln(f,'Коэффициент теплопроводности материала пластины lamda =
',lamda:6:4);
Writeln(f,'Плотность материала пластины ro = ',ro:6:4);

Страницы