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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

83
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8) во включении}
for j:= n1+2 to n1+n2 do
begin
{ai, bi, ci, fi – коэффициенты канонического представления СЛАУ с
трехдиагональной матрицей}
ai:=lamdan/sqr(hy);
bi:=2.0*lamdan/sqr(hy)+ron*cn/tau;
ci:=lamdan/sqr(hy);
fi:=-ron*cn*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+n2+1]:=2.0*an*a*tau*lamda/(2.0*an*a*tau*(lamda+lamdan
*(1-alfa[n1+n2]))+sqr(hy)*(an*lamda+a*lamdan));
beta[n1+n2+1]:=(2.0*an*a1*tau*lamdan*beta[n1+n2]+sqr(hy)*(a*lamdan
+an*lamda)*W[i,n1+n2+1])/(2.0*an*a*tau*(lamda+lamdan
*(1-alfa[n1+n2]))+sqr(hy)*(an*lamda+a*lamdan));
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8) в основном материале}
for j:= n1+n2+2 to Ny-1 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;
{определяем значение температуры на верхней границе на основе
верхнего граничного условия, в данном случае используем соотношение
(21)}
W[i,Ny]:=(2.0*a*tau*lamda*beta[Ny-1]+lamda*sqr(hy)*W[i,Ny])
/(lamda*sqr(hy)+2.0*a*tau*lamda*(1.0-alfa[Ny-1]));

Страницы