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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

88
{на нижней границе включения 2}
progonxIV(ny1+ny2+ny3+1,nx1,nx2,lamda1,0.5*(lamda1+lamda2),ro1,0.5*(
ro1+ro2),c1,0.5*(c1+c2),T);
{во включении 2}
for j:=ny1+ny2+ny3+2 to ny1+ny2+ny3+ny4 do
progonxIV(j,nx1,nx2,lamda1,lamda2,ro1,ro2,c1,c2,T);
{на верхней границе включения 2}
progonxIV(ny1+ny2+ny3+ny4+1,nx1,nx2,lamda1,0.5*(lamda1+lamda2),ro1,
0.5*(ro1+ro2),c1,0.5*(c1+c2),T);
{между включением 2 и границей y = H}
for j:=ny1+ny2+ny3+ny4+2 to Ny do
progonx(j,lamda1,ro1,c1,T);
{ СЛАУ в направлении оси Оy
в промежутке между границей x = 0 и включением 2}
for i:=2 to nx1 do
progony(i,lamda1,ro1,c1,T);
{на левой границе включения 2}
progonyIV(nx1+1,ny1+ny2+ny3,ny4,lamda1,0.5*(lamda1+lamda2),ro1,0.5*(
ro1+ro2),c1,0.5*(c1+c2),T);
{во включении 2}
for i:=nx1+2 to nx1+nx2 do
progonyIV(i,ny1+ny2+ny3,ny4,lamda1,lamda2,ro1,ro2,c1,c2,T);
{на правой границе включения 2}
progonyIV(nx1+nx2+1,ny1+ny2+ny3,ny4,lamda1,0.5*(lamda1+lamda2),ro1,
0.5*(ro1+ro2),c1,0.5*(c1+c2),T);
{в промежутке между вкл
ючением 2 и включением 3}
for i:=nx1+nx2+2 to nx1+nx2+nx3 do
progony(i,lamda1,ro1,c1,T);
{на левой границе включения 3}
progonyIV(nx1+nx2+nx3+1,ny1,ny2,lamda1,0.5*(lamda1+lamda3),ro1,0.5*(
ro1+ro3),c1,0.5*(c1+c3),T);
{во включении 3}
for i:=nx1+nx2+nx3+2 to nx1+nx2+nx3+nx4 do
progonyIV(i,ny1,ny2,lamda1,lamda3,ro1,ro3,c1,c3,T);
{на правой границе включения 3}
progonyIV(nx1+nx2+nx3+nx4+1,ny1,ny2,lamda1,0.5*(lamda1+lamda3),ro1,
0.5*(ro1+ro3),c1,0.5*(c1+c3),T);
{в промежутке между включением 3 и границей х = L}
for i:=nx1+nx2+nx3+nx4+2 to Nx-1 do
progony(i,lamda1,ro1,c1,T);
end; {цикл с предусловием окончен}

Страницы