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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

77
var {раздел описания переменных, которые мы будем использовать в
программе}
i, j, Nx, Ny : integer;
nx1, nx2, nx3, nx4 : integer;
ny1, ny2, ny3, ny4 : integer;
T : vector2;
a1, lamda1, ro1, c1 : real;
a2, lamda2, ro2, c2 : real;
a3, lamda3, ro3, c3 : real;
hx, hy, tau, t_end, time : real;
T0, L, H, Th, Tc : real;
f, g : text;
procedure progonx(j: integer; lamda,ro,c:real; var W:vector2);
{процедура, разрешающая СЛАУ с трехдиагональной матрицей
методом прогонки, в направлении оси Ох
j – номер слоя по оси у вдоль которого происходит решение СЛАУ;
lamda – коэффициент теплопроводности;
ro – плотность;
c – коэффициент теплоемкости;
W – двумерное поле температуры}
var {раздел описания локальных переменных}
i : integ
er;
alfa, beta : vector1;
ai, bi, ci, fi : real;
begin
{определяются начальные прогоночные коэффициенты на основе левого
граничного условия}
alfa[1]:=0.0;
beta[1]:=Th;
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8)}
for i:= 2 to Nx-1 do
begin
{ai, bi, ci, fi – коэффициенты канонического представления СЛАУ с
трехдиагональной матрицей}
ai:=lamda/sqr(hx);
bi:=2.0*lamda/sqr(hx)+ro*c/tau;
ci:=lamda/sqr(hx);
fi:=-ro*c*W[i,j]/tau;

Страницы