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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

81
{определяются прогоночные коэффициенты на границе включения и
основного материала}
alfa[n1+n2+1]:=2.0*an*a*tau*lamda/(2.0*an*a*tau*(lamda+lamdan*
(1-alfa[n1+n2]))+sqr(hx)*(an*lamda+a*lamdan));
beta[n1+n2+1]:=(2.0*an*a*tau*lamdan*beta[n1+n2]+sqr(hx)*(a*lamdan
+an*lamda)*W[n1+n2+1,j])/(2.0*an*a*tau*(lamda
+lamdan*(1-alfa[n1+n2]))+sqr(hx)*(an*lamda+a*lamdan));
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8) в основном материале}
for i:= n1+n2+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;
{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;
{определяем значение температуры на правой г
ранице на основе
правого граничного условия}
W[Nx,j]:=Tc;
{используя соотношение (7) определяем неизвестное поле
температуры}
for i:= Nx-1 downto 1 do
W[i,j]:=alfa[i]*W[i+1,j]+beta[i];
end; {окончание процедуры progonхIV}
procedure progonyIV(i,n1,n2:integer; lamda,lamdan,ro,ron,c,cn:real; var
W:vector2);
{процедура, разрешающая СЛАУ с трехдиагональной матрицей
методом прогонки, в направлении оси Оy при наличии включения с
отличающимися теплофизическими характеристиками
i – номер слоя по оси x вдоль которого происходит решение СЛ
АУ;
n1 – номер узла с которого начинается включение;
n2 – номер узла которым заканчивается включение;
lamda – коэффициент теплопроводности основного материала;
lamdan – коэффициент теплопроводности материала включения;
ro – плотность основного материала;

Страницы