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

UptoLike

ВУЗ: 

ТГАСУ | Томск

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

Рубрика: 

do j = 0, Ny
write(1, '(101e10.3)' ) Txy(0:Nx, j)
write(2, '(e10.3, A, e10.3)' ) j*hy, ' ', Txy(0, j)
end do
close(1)
close(2)
END
!процедура трехточечной прогонки
Subroutine TDMA(N, a, b, c, d, T)
integer n
real*8 A(0:n), B(0:n), C(0:n), D(0:n), T(0:n)
real*8 P(0:n), Q(0:n)
P(0) = b(0) / a(0)
Q(0) = d(0) / a(0)
do i = 1, N
P(i) = b(i) / (a(i) – c(i) * P(i–1))
Q(i) = (d(i) + c(i) * Q(i–1)) / (a(i) – c(i) * P(i–1))
end do
T(n) = Q(n)
do i = N–1, 0, –1
T(i) = P(i) * T(i+1) + Q(i)
end do
End
!прогоночные коэффициенты по оси OX
Subroutine ABCD_X(Nx, Axi, Bxi, Cxi, Dxi, Lams, Ros,
1 Cps, hx, tau, TY1, a0, b0, d0,
2 an, cn, dn)
integer Nx
real*8 Lams, Ros, Cps, hx, tau
20

Страницы