Численное решение двумерных нестационарных уравнений теплопроводности. Маслов Е.А. - 18 стр.

UptoLike

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

Рубрика:

hy = Ny !шаг по оcи ОУ

hy = Ly / hy ! шаг по оcи ОУ

Nxlg = int (lg / hx) !кол-во точек на гр. (ВС)

tau = K !шаг по времени

tau = Time / K !шаг по времени

TXY = T0 !нач. условие T(x,y,0)=To(x,y)

DO L = 1, K !Цикл по времени

Do j = 1, Ny–1 !направление 0X

!гр. усл.(II рода) левая гр.(AB) ao To = bo T1+do

a0 = 1.

b0 = 1.

d0 = 0.

!г.у.(III рода) правая гр.(DE) an Tn = cn Tn–1+dn

an = 1. + Alfae * hx / Lams

cn=1.

dn = Te * Alfae * hx / Lams

TY1(1:Nx–1) = TXY(1:Nx–1,j) !Ti,j

call ABCD_X(Nx, Axi, Bxi, Cxi, Dxi, Lams, Ros,

1 Cps, hx, tau, TY1, a0, b0, d0, an, cn, dn)

call TDMA(Nx, Axi, Bxi, Cxi, Dxi, TX)

TXY(0:Nx, j) = TX(0:Nx)

end Do

TXY(0:Nx, 0) = TXY(0:Nx, 1)

TXY(0:Nx, Ny) = TXY(0:Nx, Ny–1)