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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

82
ron – плотность материала включения;
c – коэффициент теплоемкости основного материала;
cn – коэффициент теплоемкости материала включения;
W – двумерное поле температуры}
Var {раздел описания локальных переменных}
j : integer;
alfa, beta : vector1;
ai, bi, ci, fi : real;
a, an : real;
begin
{определяются коэффициенты температуропроводности}
a:=lamda/(ro*c);
an:=lamdan/(ron*cn);
{определяются начальные прогоночные коэффициенты на основе
нижнего граничного условия, в данном случае используем соотношение
(20)}
alfa[1]:=2.0*a*tau/(2.0*a*tau+sqr(hy));
beta[1]:=sqr(hy)*W[i,1]/(2.0*a*tau+sqr(hy));
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8) до включения}
for j:= 2 to n1 do
begin
{ai, bi, ci, fi – коэффициенты канонического представления СЛ
АУ с
трехдиагональной матрицей}
ai:=lamda/sqr(hy);
bi:=2.0*lamda/sqr(hy)+ro*c/tau;
ci:=lamda/sqr(hy);
fi:=-ro*c*W[i,j]/tau;
{alfa[j], beta[j] – прогоночные коэффициенты}
alfa[j]:=ai/(bi-ci*alfa[j-1]);
beta[j]:=(ci*beta[j-1]-fi)/(bi-ci*alfa[j-1]);
end;
{определяются прогоночные коэффициенты на границе основного
материала и включения}
alfa[n1+1]:=2.0*a*an*tau*lamdan/(2.0*a*an*tau*(lamdan+lamda
*(1-alfa[n1]))+sqr(hy)*(a*lamdan+an*lamda));
beta[n1+1]:=(2.0*a*an*tau*lamda*beta[n1]+sqr(hy)*(a*lamdan
+an*lamda)*W[i,n1+1])/(2.0*a*an*tau*(lamdan+lamda
*(1-alfa[n1]))+sqr(hy)*(a*lamdan+an*lamda));

Страницы