NT * FI) IF (NV . EQ.5) LO = F1/(3.14 * D1 * NT * FIT) IF (EPS.GT. 0.02) GO TO 100 H = 2 * NT * LO ** TZ/(3.14 * (DC + DHAP) * NKK ) IF (NV.EQ.5) H = F1/(FP * E2) DGIDR = (DC + DHAP)/4. EL = DGIDR/D1 IF (NV.GT.5) DZITA2 = 0.03 DP2 = (DZITA2 * H * RO2 * W2 ** 2.)/(TZ * 10 ** 3.) IF (NV.EQ.5) DP2 = F2 * W2 ** 2. * F1 * RO2/(2. * FPY * 10 ** 3.) AA = EL-12. IF (AA.GT.0) MPS = 1.1 IF (AA.LT.0) MPS = .00866 * EL ** 2. - .1952 * EL +2.3234 DP1 = (DZITA1 * MPS * RO1 * LO * W1 ** 2.)/(2 * D1 * 10 ** 3.) WRITE (20, 40) ALFA1, ALFA2, KQ, F1, F2 40 FORMAT (5x, ALFA1 = , F8.2, 2x, ALFA2 = , F8.2, 2x , KQ = , F8.2, 2x, F1 = , F7.4, 2x, DZAR2 = , F10.2) WRITE ( 20, 43) DP1, DP2, H, LO, NT 43 FORMAT (5x, DP1 = , F10.3, 2x , DP2 = , F10.3, 2x, H = , F7.4, 2x, LO = , F7.4, 2x, NT = , I4) WRITE (20, 55) 55 FORMAT (15x, TABLICA NAVIVKI ) WRITE (20, 56) WRITE (20, 57) WRITE (20, 58) WRITE (20, 59) 56 FORMAT (5x, , /) 57 FORMAT (5x, NOMER DIAMETR KOLIC KOLIC DLINA ) 58 FORMAT (5x, SLOIA SLOIA VITKOV ZAHODOV TRUBKI ) NP = 0 DO 70 I = 1, NKK D (I) = DC + TR * (2 * I-1) HH = 2 * NT * D (I)/(DHAP + DC)/NKK NS = HH + .5 IF (I.EQ. NKK) NS = NT-NP NP = NP + NS JZN (I) = NS BS (I) = H/(TZ * JZN (I)) SL (I) = 3.14 * BS (I) * D (I) WRITE (20, 60) I, D (I), BS (I) , JZN (I), SL (I) 60 FORMAT (5x, I3, 2x, F7.3, 2x, F7.1, 3x, I3, 5x, F6.1) 70 CONTINUE RETURN END SUBROUTINE ALFAT (W1, D1, DC, MU1, RO1, CP1, LAMDA1, ALFA1, DZITA1, RE1) REAL MU1, LAMDA1, NU1 PR1 = CP1 * MU1/LAMDA1 RE1 = W1 * D1 * RO1/MU1 RCR = DC/2. REKR = 2300. +10500. * (D1/RCR) ** .3 IF (RE1.LE. 2300) GO TO 20 IF (RE1.GT. 2300. AND. RE1.LE. REKR) GO TO 21 IF (RE1.GT. REKR ) GO TO 22 22 NU1 = .23 * (1. +1.77 * D1/RCR) * RE1 ** .8 * PR1 ** .47) DZITA1 = .3164/RE1 ** .25 GO TO 23 21 NU1 = ( RE1/300.) * PR1 ** .37 * (1. +1.77 * D1/RCR) DZITA1 = .3164/RE1 ** .25 GO TO 23 20 NU1 = 3.66 DZITA1 = 64./RE1 23 ALFA1 = NU1 * LAMDA1/D1 WRITE (20, 72) D1, W1, RE1, NU1, PR1 FORMAT (5x, D1 = , F6 .3, 2x, W1 = , F6.3, 2x, RE1 =
