Текст подпрограммы конструкторского расчета теплообменных аппаратов heat.

PROGRAM KR COMMON /A/ TKR, PKR, ROKR, R COMMON /B/ H00, H000, s00, s000, T0 COMMON /C/ A(10), C(10) COMMON /D/ F(10) COMMON /E/ B(8,10) COMMON /M/ AM(7), BM(15), IM(15), JM(15) COMMON /SS/ SN COMMON /L/ AL(9), BL(15), IL(15), JL(15), CL(7) COMMON /T5/ FI, DE, DPR, E1, E2 COMMON /T6/ DR, DELTAR, TAUR, FIT open (5, file=’heat.dat’) READ(5, 8) NV, NK1,NK2 8 FORMAT (3I5) READ (5,9) D1, D2, TRO, TZO, DC, QT 9 FORMAT (6F10.5) READ (5, 10) T1, T2, T3, T4 ,P1, P2 10 FORMAT (6F10.5) READ (5, 11) W1, W2, GM1, GM2 , DTCR 11 FORMAT (5F10.5) READ (5, 12) FI, DE, DPR, E1, E2 12 FORMAT (5F10.5) READ (5, 13) DR, DELTAR, TAUR, FIT 13 FORMAT (4F10.5) CALL KRIAG3 CALL RTEP (D1, D2, TRO, TZO, DC, T1, T2, T3, T4 ,P1, P2, QT, W1, W2, GM1, GM2 , DTCR, NV, NK1, NK2) STOP END

SUBROUTINE RTEP(, D2, TRO, TZO, DC, T1, T2, T3, T4 ,P1, P2, QT, W1, W2, GM1, GM2 , DTCR, NV, NK1, NK2) REAL MU1, MU2, LAMDA1, LAMDA2, KQ, LO, MPS COMMON/T5/ FI, DE, DPR, E1, E2 COMMON /T6/ DR, DELTAR, TAUR, FIT COMMON /A/ TKR, PKR, ROKR, R COMMON /B/ H00, H000, s00, s000, T0 COMMON /C/ A(10), C(10) COMMON /D/ F(10) COMMON /E/ B(8,10) COMMON /M/ AM(7), BM(15), IM(15), JM(15) COMMON /SS/ SN COMMON /L/ AL(9), BL(15), IL(15), JL(15), CL(7) DIMENSION D(100), JZN(100), BS(100), SL(100), NX(25) WRITE (20, 8) 8 FORMAT (/25X,’ISHODNII PARAMERTI’/) WRITE (20,9) DTCR, D1, D2, TRO, TZO, DC FORMAT (2x,’DT=’, F6.2, 2x,’D1=’, F6.3, 2x, ’D2=’, F6.3, 2x, ’TRO=’, F6.3, 2x, ’TZO=’, F6.3, 2x, ’DC=’, F6.3) WRITE (20, 19) T1, T2, T3, T4 ,P1, P2, Q, FI, DE, DPR, E1, E2 19 FORMAT (2x, ’T1=’, F6.2, 2x, ’T2=’, F6.2, 2x, ’P1=’, F6.3, 2x, ’P2=’, F6.3, 2x, ’T3=’, F6.2, 2x, ’T4=’, F6.2, 2x, ’Q=’, F6.2, 2x, ’FI=’, F6.3, 2x, ’DE=’, F10.5, 2x, ’DPR=’, F10.5, 2x, ’E1=’, F10.5, 2x, ’ 2E=’, F10.5) WRITE (20, 119) GM1, GM2, W1, W2, NV, NK1,NK2 119 FORMAT (2x, ’GM1=’, F6.2, 2x, ’GM2=’, F6.2, 2x, ’W1=’, F6.3, 2x, ’W2=’, F6.2, 2x, ’NV =’, I3, 2x, ’NK1=’, I3, 2x, ’NK2=’, I3) IF (N.EQ. 0) CALL AER IF (N.EQ. 1) CALL N2 IF (N.EQ. 3) CALL O2 IF (N.EQ. 4) CALL HE1 IF (N.EQ. 5) CALL HE2 IF (N.EQ. 6) CALL CH4 WRITE (20, 71) 71 FORMAT (25x, ‘ RASTCHTET RTEPLO’) T= (T1+T2)/2. TT= (T3+T4)/2. RO1= FRO (P1, T) MU1= FMU (RO1, T) LAMDA1= FLA (RO1, T) CP1= CRCP (RO1, T)*1000. WRITE (6, 555) RO1, MU1, LAMDA1, CP1 555 FORMAT (2x, 4E15.7) TR= TRO*D2 TZ= TZO*D2 IF (NV.EQ.5) TZ= D2+1.5*DPR IF (NV.EQ.5) TR= .866*(D2+1.5*DPR) DELTA= TR-D2 IF (NV.EQ.6) TR= DR*TRO IF (NV.EQ.6) TZ= DR*TZO IF (NV.EQ.6) DELTA= TR-TZ XN= .75 DO 65 I=1, 25 65 NX(I)= 1+((2*I-1)**2-1)*NX NT= GM1/(.785*W1*RO1*D1**2)+.5 IF (NT.LE.1) NT=1 W1=GM1/(.785*NT*RO1*D1**2) CALL ALFAT (W1, D1, DC, MU1, RO1, CP1, LAMDA1, ALFA1, DZITA1, RE1) N=NK2 IF (N.EQ. 0) CALL AER IF (N.EQ. 1) CALL N2 IF (N.EQ. 3) CALL O2 IF (N.EQ. 4) CALL HE1 IF (N.EQ. 5) CALL HE2 IF (N.EQ. 6) CALL CH4 RO2= FRO (P2, TT) MU2= FMU (RO2, TT) LAMDA2= FLA (RO2, TT) CP2= CRCP (RO2, TT)*1000. 100 IF (NV.NE. 5) CALL ALFAM (NV, W2, TR, TZ, D2, MU2, LAMDA2, RO2, ALFA2, DZITA2, RE2) DTL= ((T1-T4)-(T2-T3))/ALOG((T1-T4)-(T2-T3)) DTL=DTCR KQ=1./(D2/(ALFA1*D1)+1./ALFA2) IF (NV.EQ.5) KQ=1./(FI/ALFA1+1./ALFA2) IF (NV.EQ.6) KQ=1./(FIT/ALFA1+1./ALFA2) F1= 11.3*Q*1000./(KQ*DTL) FYD=1.-.785*(D2**2.)/(TR*TZ) IF (NV.EQ.6) FYD= 1.-.785*(D2**2.+(DR**2.-D2**2.)*DELTAR/TAUR)/(TR*TZ) FP=GM2/(W2*RO2) IF (NV.EQ.5) FYD= E1 IF (NV.EQ.6) S=FP/ E1+.785*DC**2 S= FP/ АНВ+.785*DC**2. DH=SQRT (S/ 0.785) EO= (DH-DC)/(2.*TR) NKK= E0+ .5 IF (NV.EQ.5) NKK= .5+((DH-DC)/2.-DPR)/TR DHAP= DC+2. *TR*NKK SF= .785*DHAP**2.-DC**2.) FHS= 0.785* DHAP**2. FPY= FHS*FYD-.785*FYD*DC**2. IF (NV.EQ.5) FPY= .785*(DHAP**2.-DC**2.)*E1 W2Y= GM2/ (FPY*RO2) EPS= ABS (W2-W2Y) W2=W2Y IF (NV.EQ.5) H= F1/ (FP*E2) LO=F1/ (3.14*D2*NT) IF (NV.EQ.5) LO= F1/(3.14*D1*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=’, F8.1, 2x, ‘NU1=’, F6.2, 5x, ‘PR1=’, F6.3) RETURN END SUBROUTINE ALFAM (NV, W2, TR, TZ, D2, MU2, LAMDA2, RO2, ALFA2, DZITA2,RE2) REAL MU2, LAMDA2, NU2 RE2= W2*D2*RO2/ M2 TRO= TR/ D2 TZO=TZ/ D2 IF (NV.EQ.6.AND.RE2.LE.400.OR.RE2.GT.5000) WRITE (6, 257) IF (NV.EQ.6) GO TO 300 IF (NV.EQ.1) GO TO 30 IF (NV.EQ.2) GO TO 31 IF (NV.EQ.3) GO TO 41 30 A= .185 AM= .95 B= 8.1 BM= .2 IF (RE2.LE. 2000. OR. RE2.GT. 10000) WRITE (6, 257) GO TO 300 31 A= .083 AM= .85 B= 13. BM= .21 IF (RE2.LE. 1000. OR. RE2. GT.7000) WRITE (6, 257) NTR= 100.*TRO NTZ= 100.*TZO IF (NTR.EQ. 110. AND. NTZ. EQ.120) GO TO 431 IF (NTR.EQ. 120. AND. NTZ. EQ.120) GO TO 432 IF (NTR.EQ. 115. AND. NTZ. EQ.130) GO TO 433 431 IF (RE2.LE. 1000. OR. RE2. GT.8000) WRITE (6, 257) B= 33.8 GO TO 300 433 IF (RE2.LE. 1500. OR. RE2. GT.4000) WRITE (6, 257) B= 15.9 GO TO 300 432 IF (RE2.LE. 1000. OR. RE2. GT.26000) WRITE (6, 257) B= 13. GO TO 300 41 A= .1 AM= .88 B= 19.2 BM= .1 IF (TZO.LE. 1.8. AND. TZO. GT.1.6) GO TO 334 IF (TZO.LE. 1.6. AND. TZO. GT.1.4) GO TO 333 IF (TZO.LE. 1.8. AND. TZO. GT.1.6) GO TO 334 IF (TZO.LE. 1.2) GO TO 331 334 IF (RE2.LE. 1000. OR. RE2. GT.7000) WRITE (6, 257) A= .195 AM= .8 B= 13.7 BM= .1 GO TO 300 333 IF (RE2.LE. 1000. OR. RE2. GT.7000) WRITE (6, 257) A= .1 AM= .83 B= 17.1 BM= .1 GO TO 300 331 IF (RE2.LE. 800. OR. RE2. GT.44000) WRITE (6, 257) A= .009 AM= 1.1 B= 48.5 BM= .21 GO TO 300 332 IF (RE2.LE.10. OR. RE2. GT.7000) WRITE (6, 257) A= .1 AM= .88 B= 19.2 BM= .1 GO TO 300 300 NU2= A*RE2**AM ALFA2= NU2*LAMDA2/ D2 DZITA2= B*RE2**(-BM) 257 FORMAT (11x, ‘VNIMANIE RE2 NE SOOTVETSTVUET B’) WRITE (20, 72) D2, W2, RE2, NU2, A, B 72 FORMAT (5x, ‘D2=’, F6.3, 2x, ‘W2=’, F6.3, 2x, ‘RE2=’, F8.1, 2x, ‘NU2=’, F6.2, 2x, ‘A=’, F6.3, 2x, ‘B=’, F7.4) RETURN END SUBROUTINE ALFAR (W2, DE, MU2, LAMDA2, RO2, CP2, ALFA2, F2,RE2) REAL MU2, LANDA2 RE2=W2*RO2*DE/ MU2 PR2= CP2*MU2/ LAMDA2 ST2= .168*RE2**(-.3)*PR2**(-.666) ALFA2= ST2*W2*CP2*RO2 IF (RE2.GE.30. AND.RE2. LT.100) F2= 12.8/RE**(.64) IF (RE2.GE.100. AND.RE2. LT.5000) F2= 2.65/RE**(.3) WRITE (20, 72) DE, W2, RE2, ST2, PR2, F2 72 FORMAT (5x, ‘DE=’, F6.3, 2x, ‘W2=’, F6.3, 2x, ‘RE2=’, F8.1, 2x, ‘ST2=’, F6.2, 2x, ‘PR2=’, F6.3, 2x, ‘DZAR2=’, F10.8) RETURN END

Приложение 4.

ИСХОДНЫЕ ПАРАМЕТРЫ:

dT среднее, К = 6.5900002

d трубки внутренний, м = 8.00000038E-03

d трубки внешний, м = 9.99999978E-03

относит.шаги навивки TRO= 1.8000000 TZO= 1.0000000

DC= 0.20000000

T1= 93.199997 T2= 86.199997 P1= 0.60000002 P2= 0.13000000

T3= 80.000000 T4= 86.199997

Q= 4.5599999 FI= 0.0000000 DE= 0.0000000 DPR= 0.0000000

E1= 0.0000000 E2= 0.0000000

GM1= 0.42300001 GM2= 0.86299998

W1= 0.60000002 W2= 1.1200000

NV= 3 NK1= 1 NK2= 1

РЕЗУЛЬТАТ:

D1= 8.00000038E-03 W1= 0.59394640 RE1= 34672.715 NU1= 151.25233 PR1= 1.8764647

D2= 9.99999978E-03 W2= 1.1200000 RE2= 10306.291 NU2= 233.69812 A= 8.99999961E-03 B= 48.500000

ALFA1= 2168.7649 ALFA2= 191.19182 KQ= 172.21443 F1= 5.2233992 DZAR2= 3.21266731E-39

DP1= 3.6733856 DP2= 0.58464295 H= 0.12100913 LO= 8.7552786 NT= 19

ТАБЛИЦА НАВИВКИ

___________________________________________

Номер Диаметр Колич. Колич. Длина

слоя слоя витков заходов трубки

1 0.21799999 12.100914 1 8.2833176

2 0.25400001 12.100914 1 9.6512060

3 0.28999999 12.100914 1 11.019093

4 0.32600001 12.100914 1 12.386980

5 0.36199999 6.0504570 2 6.8774333

6 0.39800000 6.0504570 2 7.5613775

7 0.43399999 6.0504570 2 8.2453203

8 0.47000000 6.0504570 2 8.9292650

9 0.50599998 6.0504570 2 9.6132078

10 0.54200000 6.0504570 2 10.297152

11 0.57800001 4.0336380 3 7.3207307

Приложение 5.

ИСХОДНЫЕ ПАРАМЕТРЫ:

dT среднее, К = 5.6900001

d трубки внутренний, м = 8.00000038E-03

d трубки внешний, м = 9.99999978E-03

относит.шаги навивки TRO= 1.8000000 TZO= 1.0000000

DC= 0.20000000

T1= 96.800003 T2= 88.800003 P1= 0.60000002 P2= 0.13000000

T3= 86.199997 T4= 91.800003

Q= 9.9499998 FI= 0.0000000 DE= 0.0000000 DPR= 0.0000000

E1= 0.0000000 E2= 0.0000000

GM1= 0.64099997 GM2= 0.86299998

W1= 0.60000002 W2= 1.1700000

NV= 3 NK1= 0 NK2= 1

РЕЗУЛЬТАТ:

D1= 8.00000038E-03 W1= 0.60747921 RE1= 34250.098 NU1= 150.06085 PR1= 1.8840746

D2= 9.99999978E-03 W2= 1.1700000 RE2= 10173.227 NU2= 230.38127 A= 8.99999961E-03 B= 48.500000

ALFA1= 2222.8120 ALFA2= 193.49686 KQ= 174.50809 F1= 13.026825 DZAR2= 3.21266731E-39

DP1= 7.6058140 DP2= 1.5233810 H= 0.30178908 LO= 15.956425 NT= 26

ТАБЛИЦА НАВИВКИ

___________________________________________

Номер Диаметр Колич. Колич. Длина

слоя слоя витков заходов трубки

1 0.21799999 30.178907 1 20.658066

2 0.25400001 15.089454 2 12.034745

3 0.28999999 15.089454 2 13.740457

4 0.32600001 15.089454 2 15.446169

5 0.36199999 15.089454 2 17.151880

6 0.39800000 15.089454 2 18.857594

7 0.43399999 10.059636 3 13.708870

8 0.47000000 10.059636 3 14.846011

9 0.50599998 10.059636 3 15.983152

10 0.54200000 10.059636 3 17.120295

11 0.57800001 10.059636 3 18.257435

Приложение 6.

ИСХОДНЫЕ ПАРАМЕТРЫ:

dT среднее, К = 12.980000

d трубки внутренний, м = 8.00000038E-03

d трубки внешний, м = 9.99999978E-03

относит.шаги навивки TRO= 1.2000000 TZO= 1.2000000

DC= 0.16000000

T1= 295.89999 T2= 278.00000 P1= 20.000000 P2= 0.11000000

T3= 258.29999 T4= 287.89999

Q= 24.969999 FI= 0.0000000 DE= 0.0000000 DPR= 0.0000000

E1= 0.0000000 E2= 0.0000000

GM1= 1.0640000 GM2= 0.86299998

W1= 2.0000000 W2= 7.8099999

NV= 2 NK1= 0 NK2= 1

РЕЗУЛЬТАТ:

D1= 8.00000038E-03 W1= 1.9781258 RE1= 157706.05 NU1= 361.61508 PR1= 0.85265493

D2= 9.99999978E-03 W2= 7.8099999 RE2= 6359.2119 NU2= 141.89627 A= 8.29999968E-02 B= 13.000000

ALFA1= 1658.6157 ALFA2= 344.44308 KQ= 273.45734 F1= 9.1452923 DZAR2= 3.21266731E-39

DP1= 6.5771208 DP2= 3.5749595 H= 0.24239165 LO= 6.4722519 NT= 45

ТАБЛИЦА НАВИВКИ

___________________________________________

Номер Диаметр Колич. Колич. Длина

слоя слоя витков заходов трубки

1 0.17199999 10.099652 2 5.4546199

2 0.19599999 10.099652 2 6.2157302

3 0.22000000 10.099652 2 6.9768400

4 0.24400000 10.099652 2 7.7379498

5 0.26800001 6.7331014 3 5.6660399

6 0.29200000 6.7331014 3 6.1734462

7 0.31599998 6.7331014 3 6.6808524

8 0.34000000 6.7331014 3 7.1882591

9 0.36399999 5.0498261 4 5.7717495

10 0.38800001 5.0498261 4 6.1523046

11 0.41200000 5.0498261 4 6.5328593

12 0.43599999 5.0498261 4 6.9134140

13 0.46000001 4.0398607 5 5.8351750

14 0.48400000 5.0498261 4 7.6745238

Приложение 7.