Основы автоматизированного проектирования химических производств. Миронов В.М - 62 стр.

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

62
' Начало метода прогонки по алгоритму обратного хода
j = 1
' Определение шага
oprh:
h = (xn - x0) / n
'Вычисление коэффициентов p, q, f дифференциального уравнения
FOR i = 0 TO n
x(i, j) = x0 + i * h
SELECT CASE T 'вычисление толщины z и ее производной zs=dz/dx
CASE 1 ' плоский диск
z(i, j) = S1
zs(i, j) = 0
CASE 2 ' конический диск
z(i, j) = S1 - (S1 - S2) * (x(i, j) - x0) / (xn - x0)
zs(i, j) = (S2 - S1) / (xn - x0)
CASE 3 ' гиперболический диск
v = LOG(S1 / S2) / LOG(xn / x0)
z(i, j) = S1 * (x0 / x(i, j)) ^ v
zs(i, j) = -v * S1 * x0 ^ v / x(i, j) ^ (v + 1)
END SELECT
p(i, j) = 1 / x(i, j) + zs(i, j) / z(i, j)
q(i, j) = mi * zs(i, j) / z(i, j) / x(i, j) - 1 / x(i, j) ^ 2
f(i, j) = -x(i, j) * ro * w ^ 2 / a1
NEXT 'конец цикла "Вычисление коэффициентов p, q, f уравнения"
'Вычисление коэффициентов m(i,j), k(i,j)
FOR i = 0 TO n - 2
m(i, j) = -2 + h * p(i, j)
k(i, j) = 1 - h * p(i, j) + h * h * q(i, j)
NEXT 'конец цикла "Вычисление коэффициентов m(i,j), k(i,j)"
'Вычисление коэффициентов c(0,j), d(0,j)
c(0, j) = (a1 - a0 * h) / (m(0, j) * (a1 - a0 * h) + k(0, j) * a1)
d(0, j) = k(0, j) * G1 * h / (a1 - a0 * h) + h * h * f(0, j)
FOR i = 1 TO n - 2 'Вычисление коэффициентов c(i,j), d(i,j)
c(i, j) = 1 / (m(i, j) - k(i, j) * c(i - 1, j))
d(i, j) = h * h * f(i, j) - k(i, j) * c(i - 1, j) * d(i - 1, j)
NEXT 'конец цикла "Вычисление коэффициентов c(i,j), d(i,j)"
'Вычисление y(n,j)
y(n, j) = (b1 * c(n - 2, j) * d(n - 2, j) + G2 * h)
y(n, j) = y(n, j) / (b1 * (1 + c(n - 2, j)) + b0 * h)

Страницы