Coker A.K. Fortran Programs for Chemical Process Design, Analysis, and Simulation
Подождите немного. Документ загружается.
C
C
C
22
C
C
C
33
C
C
C
44
55
18
12
19
13
C
66
C
C
77
C
C
88
READ TANK DIAMETER, D2: in.
READ LIQUID DEPTH, H2: in.
READ TANK LENGTH, L2: in.
READ (3, *, ERR=I9) D2, H2, L2
GO TO 77
READ TANK DIAMETER, D3: in.
READ LIQUID DEPTH, H3: in.
READ TANK LENGTH, L3: in.
READ (3, *, ERR=I9) D3, H3, L3
GO TO 88
READ TANK DIAMETER, D4: in.
READ LIQUID DEPTH, H4: in.
READ TANK LENGTH, L4: in.
READ (3, *, ERR=I9) D4, H4, L4
GO TO 99
READ (3, *, ERR=I9) D5, H5, L5
GO TO Iii
WRITE (*, 12)
FORMAT (6X, "DATA FILE DOES NOT EXIST')
GO TO 999
WRITE (*', 13)
FORMAT(6X, 'ERROR MESSAGE IN THE DATA VALUE')
GO TO 999
CALCULATE THE VOLUME OF THE TANK WITH FLAT HEADS
CALL FLAT
PRINT THE RESULTS OF THE TANK WITH FLAT HEADS
CALL OFLAT
GO TO 5
CALCULATE THE VOLUME OF THE TANK WITH DISHED HEADS
CALL DISHED
PRINT THE RESULTS OF THE TANK WITH DISHED HEADS
CALL ODISH
GO TO 5
CALCULATE THE VOLUME OF THE TANK WITH ELLIPTICAL HEADS
CALL ELLIP
307
C
99
C
C
iii
C
C
999
PRINT THE RESULTS OF THE TANK WITH ELLIPTICAL HEADS.
CALL OELLIP
GO TO 5
CALCULATE THE VOLUME OF THE TANK WITH HEMISPHERICAL HEADS
CALL HEMISP
PRINT THE RESULTS OF THE TANK WITH HEMISPHERICAL HEADS.
CALL OHEMIS
GO TO 5
CALCULATE THE PARTIAL VOLUME OF A VERTICAL TANK
CALL VERT
PRINT THE RESULTS OF THE PARTIAL VOLUME OF A VERTICAL
VESSEL
CALL OVERT
GO TO 999
CLOSE (3, STATUS='KEEP')
CLOSE (i)
STOP
END
THIS PROGRAM CALCULATES THE PARTIAL VOLUME FOR A HORIZONTAL
CYLINDRICAL VESSEL WITH FLAT HEADS
SUBROUTINE FLAT
REAL L1
COMMON PI
COMMON~DATA1~
DI, HI, L1
COMMON/PROGI/
VHC, VHCCFT, VHCGAL, VHCBBL
R1 = DI/2.
CLI = 2.*SQRT(DI*HI-(HI**2))
VALI = (RI-HI)/RI
VAL2 = (RI-HI)*SQRT(2.*RI*HI-(HI**2))
ACI = ((RI**2)*ACOS(VALI))-VAL2
AREAl = (0.25*PI*DI**2)
VHC = ACI*LI
CONVERT THE RESULTS TO FT^3, GALLONS AND BARRELS
CALL CONVT (VHC, VHCCFT, VHCGAL, VHCBBL)
RETURN
308
END
***************************************************************
THIS PROGRAM CALCULATES THE PARTIAL VOLUME FOR A HORIZONTAL
CYLINDRICAL VESSEL WITH DISHED HEADS
**************************************************************
SUBROUTINE DISHED
REAL L2
COMMON PI
COMMON/DATA2/ D2, H2, L2
COMMON/PROG2/
VT, VTCFT, VTGAL, VTBBL
R2 = D2/2.
CL2 = 2.*SQRT(D2*H2-(H2**2))
VAL3 = (R2-H2)/R2
VAL4 = (R2-H2)*SQRT(2.*R2*H2-(H2**2))
AC2 = ((R2**2)*ACOS(VAL3))-VAL4
VHC2 = AC2*L2
VDH = 0.215483*H2*(I.5*D2-H2)
VT = 2.*VDH+VHC2
CONVERT THE RESULTS TO CUBIC FEET, GALLONS AND BARRELS.
CALL CONVT (VT, VTCFT, VTGAL, VTBBL)
RETURN
END
***************************************************************
THIS PROGRAM CALCULATES THE PARTIAL VOLUME FOR A HORIZONTAL
CYLINDRICAL VESSEL WITH ELLIPTICAL HEADS
***************************************************************
SUBROUTINE ELLIP
REAL L3
COMMON PI
COMMON/DATA3/ D3, H3, L3
COMMON/PROG3/VTELL,VTECFT,VTEGAL,VTEBBL
R3 = D3/2.
CL3 = 2.*SQRT(D3*H3-(H3**2))
VAL5 = (R3-H3)/R3
VAL6 = (R3-H3)*SQRT(2.*R3*H3-(H3**2))
AC3 = ((R3**2)*ACOS(VAL5))-VAL6
VHC3 = AC3*L3
VELL = 0.5236*(H3**2)*(I.5*D3-H3)
VTELL = 2.*VELL+VHC3
CONVERT THE RESULTS TO CUBIC FEET, GALLONS AND BARRELS.
CALL CONVT (VTELL, VTECFT, VTEGAL, VTEBBL)
RETURN
END
*************************************************************
309
THIS PROGRAM CALCULATES THE PARTIAL VOLUME FOR A HORIZONTAL
CYLINDRICAL VESSEL WITH HEMISPHERICAL HEADS
SUBROUTINE HEMISP
REAL L4
COMMON PI
COMMON/DATA4/ D4, H4, L4
COMMON/PROS4/ VTHS, VTHCFT, VTHGAL, VTHBBL
R4 = D4/2.
CL4 = 2.*SQRT(D4*H4-(H4**2))
VAL7 = (R4-H4)/R4
VAL8 = (R4-H4)*SQRT(2.*R4*H4-(H4**2))
AC4 = ((R4**2)*ACOS(VAL7))-VAL8
VHC4 = AC4*L4
YELL4 = 0.5236*(H4**2)*(I.5*D4-H4)
VHS = 2.*VELL4
VTHS = 2.*VHS+VHC4
CONVERT THE RESULTS TO CUBIC FEET, GALLONS AND BARRELS.
CALL CONVT (VTHS, VTHCFT, VTHGAL, VTHBBL)
RETURN
END
THIS PROGRAM CALCULATES THE PARTIAL VOLUME FOR A VERTICAL
CYLINDRICAL VESSEL
SUBROUTINE VERT
COMMON PI
COMMON/DATA5/ D5, H5, L5
COMMON/PROG5/ VVT, VVCFT, VVGAL, VVBBL
VVT = 0.25*PI*(D5**2)*H5
CONVERT THE RESULTS TO CUBIC FEET, GALLONS AND BARRELS.
CALL CONVT (VVT, VVCFT, VVGAL, VVBBL)
RETURN
END
THIS PROGRAM CONVERTS THE RESULTS TO CUBIC FEET, GALLONS, AND
BARRELS.
SUBROUTINE CONVT (VOL, VLFT3, VLGAL, VLBAR)
VLFT3 = VOL/1728.0
310
i00
ii0
1
2
120
1
2
3
130
VLGAL = VOL*O.004329
VLBAR = VOL*0.00010307
RETURN
END
***************************************************************
THIS PROGRAM PRINTS THE RESULTS OF THE TANK VOLUME WITH FLAT
HEADS
***************************************************************
SUBROUTINE OFLAT
REAL L1
COMMON/DATA1/ DI, HI, L1
COMMON/PROGI/ VHC, VHCCFT, VHCGAL, VHCBBL
WRITE (i, i00)
FORMAT (//,25X, 'TANK VOLUME WITH FLAT HEADS',/,78(IH*))
WRITE (i, ii0) DI, HI, L1
FORMAT (5X,'TANK DIAMETER, in.:', T60, F9.3,/,
5X,'LIQUID DEPTH, in.:', T60, F9.3,/,
5X,'TANK LENGTH, in.:', T60, F9.3)
WRITE (i, 120) VHC, VHCCFT, VHCGAL, VHCBBL
FORMAT (5X,'TANK VOLUME, in.^3: ' , T60, E14.6,/,
5X,'TANK VOLUME, ft.^3: ' , T60, E14.6,/,
5X,'TANK VOLUME, gal:', T60, E14.6,/,
5X,'TANK VOLUME, bbl:', T60, E14.6)
WRITE (i, 130)
FORMAT (78(IH-))
140
150
RETURN
END
***************************************************************
THIS PROGRAM PRINTS THE RESULTS OF THE TANK VOLUME WITH DISHED
HEADS.
***************************************************************
SUBROUTINE ODISH
REAL L2
COMMON/DATA2/ D2, H2, L2
COMMON/PROG2/ VT, VTCFT, VTGAL, VTBBL
WRITE(I, 140)
FORMAT (//,25X,'TANK VOLUME WITH DISHED HEADS',/,78(IH*))
WRITE (I, 150) D2, H2, L2
FORMAT (5X,'TANK DIAMETER, in.:', T60, F9.3,/,
5X,'LIQUID DEPTH, in.:', T60, F9.3,/,
5X,'TANK LENGTH, in.:', T60, F9.3)
311
160
170
180
190
200
210
220
WRITE (i, 160) VT, VTCFT, VTGAL, VTBBL
FORMAT (5X,'TANK VOLUME, in^3: ' , T60, E14.6,/,
1 5X,'TANK VOLUME, ft.^3: ' , T60, E14.6,/,
2 5X,'TANK VOLUME, gal:', T60, E14.6,/,
3
5X,'TANK VOLUME, bbl:', T60, E14.6)
WRITE (I, 170)
FORMAT (78(IH-))
RETURN
END
***************************************************************
THIS PROGRAM PRINTS THE RESULTS OF THE TANK VOLUME WITH
ELLIPTICAL HEADS
***************************************************************
SUBROUTINE OELLIP
REAL L3
COMMON/DATA3/ D3, H3, L3
COMMON/PROG3/
VTELL, VTECFT, VTEGAL, VTEBBL
WRITE (i, 180)
FORMAT (//,25X, 'TANK VOLUME WITH ELLIPTICAL HEADS',/,78(IH*))
WRITE (i, 190) D3, H3, L3
FORMAT (5X,'TANK DIAMETER, in.:', T60, F9.3,/,
5X,'LIQUID DEPTH, in.:', T60, F9.3,/,
5X,'TANK LENGTH, in.:', T60, F9.3)
WRITE (i, 200) VTELL, VTECFT, VTEGAL, VTEBBL
FORMAT (5X,'TANK VOLUME, in.^3: ' , T60, E14.6,/,
5X,'TANK VOLUME, ft.^3: ' , T60, E14.6,/,
5X,'TANK VOLUME, gal:', T60, E14.6,/,
5X,'TANK VOLUME, bbl:', T60, E14.6)
WRITE (i, 210)
FORMAT (78(IH-))
RETURN
END
***************************************************************
THIS PROGRAM PRINTS THE RESULTS OF THE TANK VOLUME FOR
HEMISPHERICAL HEADS.
***************************************************************
SUBROUTINE OHEMIS
REAL L4
COMMON~DATA4~
D4, H4, L4
COMMON/PROG4/ VTHS, VTHCFT, VTHGAL, VTHBBL
WRITE (i, 220)
FORMAT
(//,25X,'T/~IK VOLUME WITH HEMISPHERICAL HEADS',/,78(1H*))
312
230
240
250
260
270
280
290
WRITE (i, 230) D4, H4, L4
FORMAT (SX,'TANK DIAMETER, in.:', T60, F9.3,/,
5X,'LIQUID DEPTH, in.:', T60, F9.3,/,
5X,'TANK
LENGTH, in.:', T60, F9.3)
WRITE (i, 240) VTHS, VTHCFT, VTHGAL, VTHBBL
FORMAT (5X,'TANK VOLUME, in.^3: ' , T60, E14.6,/,
5X,'TANK VOLUME, ft.^3: 9 , T60, E14.6,/,
5X,'TANK
VOLUME, gal:', T60, E14.6,/,
5X,'TANK
VOLUME, bbl:', T60, E14.6)
WRITE (i, 250)
FORMAT (78(IH-))
RETURN
END
***************************************************************
THIS PROGRAM PRINTS THE RESULTS OF THE PARTIAL VOLUME FOR A
VERTICAL CYLINDRICAL VESSEL
***************************************************************
SUBROUTINE OVERT
REAL L5
COMMON~DATA5~
D5, H5, L5
COMMON/PROG5/ VVT, WCFT, WGAL, VVBBL
WRITE (i, 260)
FORMAT (//,25X,'VERTICAL TANK VOLUME', /, 78(IH*))
WRITE (r, 270) DS, HS, L5
FORMAT (SX,'TANK DIAMETER, in.:', T60, F9.3,/,
5X,'LIQUID DEPTH, in.:', T60, F9.3,/,
5X,'TANK
LENGTH, in.:', T60, F9.3)
WRITE (I, 280) VVT, VVCFT, VVGAL, WBBL
FORMAT (5X,'TANK VOLUME, in.^3: ' , T60, E14.6,/,
5X,'TANK
VOLUME, ft.^3: ' , T60, E14.6,/,
5X,'TANK
VOLUME, gal:', T60, E14.6,/,
5X,'TANK VOLUME, bbl:', T60, E14.6)
WRITE (i, 290)
FORMAT (78(IH-))
FORM FEED THE PRINTING PAPER TO THE TOP OF THE NEXT PAGE.
WRITE (i, *) CHAR(12)
RETURN
END
313
PROGRAM PROG43
*****************************************************************
A PROGRAM FOR PRELIMINARY VESSEL DESIGN.
THIS PROGRAM WILL ESTIMATE SHELL THICKNESS AND WEIGHT,
HEAD THICKNESS AND WEIGHT FOR CYLINDRICAL PRESSURE VESSELS.
IN ADDITION, HEAD THICKNESS AND WEIGHT ARE CALCULATED FOR
HEMISPHERICAL, ELLIPSOIDAL AND TORISPHERICAL HEADS. TOTAL VESSEL
WEIGHT IS CALCULATED (SHELL AND HEADS) EXCLUSIVE OF NOZZLES,
ATTACHMENTS AND VESSEL INTERNALS.
*****************************************************************
TKS = SHELL THICKNESS, in.
TKHH = HEMISPHERICAL HEAD THICKNESS, in.
TKELL = ELLIPSOIDAL HEAD THICKNESS, in.
TKTOR = TORISPHERICAL HEAD THICKNESS, in.
P = INTERNAL DESIGN PRESSURE, psig.
R = INSIDE RADIUS, in.
S = ALLOWABLE STRESS, ib/in^2.
E = JOINT EFFICIENCY
C = CORROSION ALLOWANCE, in.
ID = INSIDE DIAMETER OF SHELL OR HEAD, in.
RC = CROWN RADIUS OF TORISPHERICAL HEAD, in.
LN = SHELL LENGTH, ft.
DENV = DENSITY OF VESSEL MATERIAL, ib/ft^3.
OD = OUTSIDE DIAMETER OF HEAD, in.
HFAC = HEAD BLANK DIAMETER FACTOR
SF = STRAIGHT FLANGE LENGTH, in.
BD = HEAD BLANK DIAMETER, in.
WTH = WEIGHT OF HEAD, ib
WTS = WEIGHT OF VESSEL SHELL, lb.
TOTWT = TOTAL VESSEL WEIGHT, lb.
***************************************************************
REAL LNI, LN2, LN3, IDI, ID2, ID3
COMMON PI
COMMON/DATA1/ IDI, LNI, PI, SFI, Cl, DENS1, ESI, EHI, S1
COMMON/RESI/ TKSI, TKHI, WTSI, WTHI, TOTWTI
COMMON/DATA2/ ID2, LN2, P2, SF2, C2, DENS2, ES2, EH2, $2
COMMON/RES2/ TKS2, TKHS2, WTS2, WTH2, TOTWT2
COMMON~DATA3~
ID3, LN3, P3, SF3, C3, DENS3, ES3, EH3, $3, RC
COMMON/RES3/
TKS3, TKHS3, WTS3, WTH3, TOTWT3
OPEN (UNIT=3, FILE='DATA43.DAT', STATUS='OLD', ERR=IS)
OPEN (UNIT=l, FILE='PRN')
PI = 3.1415927
READ (3, *, ERR=I9) I
GO TO (ii, 22, 33), I
READ INTERNAL DIAMETER, IDI: in.
READ SHELL LENGTH, LNI: in.
READ INTERNAL DESIGN PRESSURE, PI: psig.
READ STRAIGHT FLANGE LENGTH, SFI: in.
314
ii
READ CORROSION ALLOWANCE, CI: in.
READ DENSITY OF VESSEL MATERIAL, DENS1: ib/ft^3.
READ JOINT EFFICIENCY OF SHELL, ESI
READ JOINT EFFICIENCY OF HEAD, EHI
READ ALLOWABLE STRESS, SI: Ib/in^2.
READ (3, *
READ (3, *
READ (3, *
, ERR=I9) IDI, LNI, P1
, ERR=I9) SFI, CI, DENSl
, ERR=I9) ESl, EHI, Sl
GO TO 44
READ INTERNAL DIAMETER, ID2: in.
READ SHELL LENGTH, LN2: in.
READ INTERNAL DESIGN PRESSURE, P2: psig.
READ STRAIGHT FLANGE LENGTH, SF2: in.
READ CORROSION ALLOWANCE, C2: in.
READ DENSITY OF VESSEL MATERIAL, DENS2: Ib/ft^3.
READ JOINT EFFICIENCY OF SHELL, ES2:
READ JOINT EFFICIENCY OF HEAD, EH2:
READ ALLOWABLE STRESS, $2: ib/in^2.
22
33
18
13
19
14
44
READ (3, *
READ (3, *
READ (3, *
, ERR=I9) ID2, LN2, P2
, ERR=I9) SF2, C2, DENS2
, ERR=I9) ES2, EH2, $2
GO TO 55
READ INTERNAL DIAMETER, ID3: in.
READ SHELL LENGTH, LN3: in.
READ INTERNAL DESIGN PRESSURE, P3: psig.
READ STRAIGHT FLANGE LENGTH, SF3: in.
READ CORROSION ALLOWANCE, C3: in.
READ DENSITY OF VESSEL MATERIAL, DENS3: ib/ft^3.
READ JOINT EFFICIENCY OF SHELL, ES3:
READ JOINT EFFICIENCY OF HEAD, EH3:
READ CROWN RADIUS OF TORISPHERICAL HEAD, RC: in.
READ (3, *, ERR=I9) ID3, LN3, P3
READ (3, *, ERR=I9) SF3, C3, DENS3
READ (3, *, ERR=I9) ES3, EH3, S3, RC
GO TO 66
WRITE (*, 13)
FORMAT (6X, 'DATA FILE DOES NOT EXIST')
GO TO 999
WRITE (*, 14)
FORMAT (6X, 'ERROR MESSAGE IN THE DATA VALUE')
GO TO 999
CALCULATE THE SHELL THICKNESS, WEIGHT AND HEAD THICKNESS FOR
HEMISPHERICAL HEADS.
CALL HEMISP
315
C
C
55
C
66
999
PRINT THE RESULTS OF THE VESSEL WITH HEMISPHERICAL HEADS.
CALL OHEMISP
GO TO 5
CALCULATE THE SHELL THICKNESS, WEIGHT AND HEAD THICKNESS FOR
ELLIPSOIDAL HEADS.
CALL ELLIP
PRINT THE RESULTS OF THE VESSEL WITH ELLIPSOIDAL HEADS.
CALL OELLIP
GO TO 5
CALCULATE THE SHELL THICKNESS, WEIGHT AND HEAD THICKNESS FOR
TORISPHERICAL HEADS.
CALL TORIS
PRINT THE RESULTS OF THE VESSEL WITH TORISPHERICAL HEADS.
CALL OTORIS
GO TO 999
CLOSE (3, STATUS='KEEP')
CLOSE (i)
STOP
END
**************************************************************
VESSEL DESIGN WITH HEMISPHERICAL HEADS
**************************************************************
SUBROUTINE HEMISP
REAL LNI,IDI,ODI
COMMON PI
COMMON~DATA1/
IDI, LNI, PI, SFI, CI, DENS1, ESl, EHI, S1
COMMON/RESl/ TKSI, TKHI, WTSI, WTHI, TOTWTI
CALCULATE SHELL THICKNESS, in
R1 = IDI/2.
TKSl = (PI*RI)/(SI*ESI-0.6*PI)+CI
ODI = IDI + 2.0*TKSl
CALCULATE WEIGHT OF VESSEL SHELL, ib
WTSl = ((ODI+TKSl)*PI*LNI*TKSI*DENSl)/144.
CALCULATE WEIGHT OF HEAD, ib
TKHI = (PI*RI)/(2.*SI*EHI-0.2*PI)
316