n> ST2 = , F6.2, 2x, PR2 = , F6.3, 2x, < span align = "justify"> DZAR2 = , F10.8) RETURN END
Додаток 2
Текст підпрограми конструкторського розрахунку конденсатора-випарника з напиленням-оребрених труб RKINOT.
SUBROUTINE
character (80) trez! мінлива для повідомлень про помилки
character (10) prodcase (NK1) (1); prod = "Aer" (2); prod = "N2" (3); prod = "Ar" (4); prod = "O2" (5); prod = "He1" (6); prod = "He2" agent (prod)! обчислення властивостей кріопродукта = Ts (p1)! tисп у верхньому перетині
R1 = R (T1) * 1000. ! Теплотапаротворення = rov (p1, T1)! Щільність пари = rol (p1, T1)! Щільність рідини
select case (NK2) (1); prod = "Aer" (2); prod = "N2" (3); prod = "Ar" (4); prod = " ; O2 "(5); prod =" He1 "(6); prod =" He2 "agent (prod) 2 = T1 + dTv/2. ! Тконд середнє у верхньому перетині = ps (T2)! Рконд у верхньому перетині = R (T2) * 1000. ! Теплотапаротворення = rov (p2, T2)! Щільність пари = rol (p2, T2)! Щільність рідини = fmu (rol2, T2)! Динамічна в'язкість = fla (rol2, T2)/1000. ! Теплопровідність = 7000. ! Можна задати у вихідних даних або оцінити як в RKIMT
wp1 = q01 * d1 * Htr/(R1 * rov1 * (995. * st ** 2.-900. * d1 ** 2.))! привед. швидкість (3) = (SQRT (1. +5.4 * wp1 * rol1) -1.)/(2.7 * wp1 * rol1)! відносний рівень (2) = T1 * Htr * Ho * rol1/(.102 * R1 * rov1) ! температурна депресія (1) = dTv-0.5 * dlT! середній температурний напір (4)
! комплекс E у формулі (20) Z = E * dTконд
! Lambda ребра прийнято const = 89, полутолщіна ребра (t) - const = 0.001 = 8866.8! коефіцієнт у (16) = 3.14 * d1 * Htr! поверхню кипіння 1 трубки
! пов. конденсації - ребро + міжреберні пр-во на од. висоти трубки
Fr1 = dr + br +2. * sqrt (((cr-br)/2.) ** 2 + hr ** 2)
F2 = Nr * Fr1 * Htr! поверхню конденсації 1 трубки = Nr * Ar/(C * F1)! перші два дроби в (26) = .5 * d1/89. * log (d1 /d2) * c * d1 * 2./(d1 + d2)! множник в (25) = 1.5; dTa = .1; eps = .01! діапазон пошуку dTкіп і задана похибка = bisection (dTi, dTa, sdT, eps) расчет dTкіп половинним діленням = dTsr-dT1-GCF * dT1 ** 1.4655! dTконд з балансу температур = C * F1 * dT1 ** 1.4655! Qкіп 1 трубки = Nr * Ar * dT2 ** 0.875 ! Qконд 1 трубки = Q * 1000./Q1! кількість трубок = Q1/(F1 * dT1)! коефіцієнт тепловіддачі при кипінні = Q2/(F2 * dT2)! коефіцієнт тепловіддачі при конденсації (1, *) (1 , *) 'РЕЗУЛЬТАТИ КОНСТРУКТОРСЬКОГО РОЗРАХУНКУ' (1, *) (1, *) 'Температура у верхній частині випарника, К', T1...
