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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

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

Страницы