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

UptoLike

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

63
'Вычисление y(i,j) обратным ходом
FOR i = n - 1 TO 1 STEP -1
y(i, j) = c(i - 1, j) * (d(i - 1, j) - y(i + 1, j))
NEXT
'Вычисление y(0,j)
y(0, j) = (a1 * y(1, j) - G1 * h) / (a1 - a0 * h)
'Проверка требуемой точности вычислений
IF j = 1 THEN
j = j + 1: n = 2 * n: GOTO oprh
ELSE
FOR i = 0 TO n / 2'STEP 2 ^ (j - 1)
IF ABS(y(i, j - 1) - y(i / 2 ^ (j - 1), j)) < 3 * eps THEN
'next
ELSE
j = j + 1: n = 2 * n: GOTO oprh
END IF
NEXT
END IF
'Конец метода прогонки по алгоритму обратного хода
'Вычисление значений радиальных gr(i) и окружных gt(i) напряжений
FOR i = 0 TO n
YS(i) = (y(i + 1, j) - y(i, j)) / h:
IF i = 0 THEN
YS(i) = (y(1, j) - y(0, j)) / h
ELSEIF i = n THEN
YS(i) = (y(n, j) - y(n - 1, j)) / h
END IF
gr(i) = a1 * (YS(i) + mi * y(i, j) / x(i, j))
gt(i) = a1 * (y(i, j) / x(i, j) + mi * YS(i))
NEXT i
'Вывод результатов в виде таблицы
tabl:
COLOR 7, 0: CLS
LOCATE 1, 10
PRINT "Радиальные и тангенциальные напряжения в диске"
LOCATE 3, 5: PRINT " Номер"
LOCATE 4, 5: PRINT " сечения"
LOCATE 3, 15: PRINT " Радиус"
LOCATE 4, 16: PRINT " сечения,мм"
LOCATE 3, 35: PRINT "Gr,МПа"