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

UptoLike

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

Рубрика:

3.4. Листинг программы решения рассматриваемой

задачи (язык программирования Fortran Powerstation)

Program Tempreture_Dim2

Integer Nx,Ny,K

Real*8 Lx, Ly, lg, Ros, Cps, Lams, Alfae, Alfag, T0,

1 Tg, Te, Time

Parameter (Lx=0.1, Ly=0.15, lg=0.05, Nx=100,

Ny=100, K=100, Ros=1800., Cps=840.,

Lams=1.2, Alfae=30., Alfag=2500.,

T0=300., Tg=3000., Te=300., Time=25.)

Integer Nxlg

Real*8 hx, hy, tau

Real*8 a0, b0, d0, an, cn, dn

Real*8 TXY(0:Nx, 0:Ny)

Real*8 TX(0:Nx), TY1(1:Nx–1)

Real*8 TY(0:Ny), TX1(1:Ny–1)

Real*8 Axi(0:Nx), Bxi(0:Nx), Cxi(0:Nx), Dxi(0:Nx)

Real*8 Ayj(0:Ny), Byj(0:Ny), Cyj(0:Ny), Dyj(0:Ny)

Open(1,file='T(x,y).txt')

Open(2,file='T(0,y).txt')

hx = Nx !шаг по оси OX

hx = Lx / hx !шаг по оси OX