1 COMMON/A/AW, HO, HO0, DHO, EO, SO, TO, HU, HU0, DHU, *EU, SU, TU, WO, WU, N, | LO, LU, FO, FU, MO, MU | 2 COMMON/B/EY, ELAST | 3 COMMON/C/NN, TK, NK | 4 REAL LO, LU, MO, MU, N, NK | 5 DIMENSION TK(20), NK(20) | 6 OPEN(5,FILE = “TIEIN.TXT”) | 7 OPEN(6,FILE=“TIEOUT.TXT”) | 8 READ(5,1)LO, LU, HO, HU, DC, DW, DENS | 9 1 FORMAT(7F10.0) | 10 AC = 3.14159*DC**2/4. | 11 MO = DENS*AC*LO/2./10000. | 12 WO = −MO*10. | 13 MU = DENS*AC*LU/2./10000. | 14 WU = *10. | 15 AW = 3.14159*DW**2/4./10000. | 16 READ(5,1)SY, ELAST, SC0, SW0 | 17 EY = SY/ELAST | 18 = SC0*AC | 19 TO = SW0*AW | 20 TU = SW0*AW | 21 READ(5, 2)NN, TMAX | 22 2 FORMAT(I5, F10.0) | 23 READ(5,3)(TK, NK, = 1,NN) | 24 3 FORMAT(40F5.2) | 25 STIF1 = 2. */(LO/2.) + ELAST*AW/HO | 26 STIF2 = 2. */(LU/2.) + ELAST*AW/HO + ELAST*AW/HU | 27 STIF12= −ELAST*AW/HO | 28 = MO*MU | 29 = −(STIF1*MU + STIF2*MO) | 30 = STIF1*STIF2-STIF12**2 | 31 = SQRT(B**. *A*C) | 32 =SQRT((/(2. *)) | 33 | 34 | 35 | 36 DO 4 = 1, NN | 37 4 TK = TK* | 38 DT=/(4. * 3.14159) | 39 = 0. | 40 UO = (WO−TO)/(2. *N/(LO/2.)) | 41 UU = (WU + TO − TU)/(2. *N/(LU/2.)) | 42 VO = 0. | 43 VU = 0. | 44 EOPL = 0. | 45 EUPL = 0. | 46 EW0 = SW0/ELAST | 47 HO0 = (HO + UO − UU)/(1. + EW0) | 48 HU0 = (HU + UU)/(1. + EW0) | 49 CALL EVAL(UO, UU, EOPL, EUPL, GO, GU) | 50 5 = + DT | 51 CALL NHIST | 52 UOP = UO + VO*DT | 53 UUP = UU + VU*DT | 54 VOP = VO + GO*DT | 55 VUP = VU + GU*DT | 56 EOPLP = EOPL | 57 EUPLP = EUPL | 58 CALL EVAL (UOP, UUP, EOPLP, EUPLP, GOP, GUP) | 59 UO1 = UO + (VO + VOP)/2. *DT | 60 UU1= UU + (VU + VUP)/2. *DT | 61 VO1 = VO + (GO + GOP)/2. * DT | 62 VU1 = VU + (GU + GUP)/2. *DT | 63 EOPL1 = EOPL | 64 EUPL1 = EUPL | 65 CALL EVAL (UO1, UU1, EOPL1, EUPL1, GO1, GU1) | 66 UO =UO + (VO+VO1)/2. *DT | 67 UU = UU + (VU + VU1)/2. *DT | 68 VO = VO + (GO+GO1)/2. *DT | 69 VU = VU + (GU + GU1)/2. *DT | 70 CALL EVAL (UO, UU, EOPL, EUPL, GO, GU) | 71 WRITE (6,6)T, N, UO, UU, TO, TU | 72 WRITE (6,6)EO, SO, EU, SU | 73 6 FORMAT(1X,5(E10.4,1X)) | 74 IF(T.GT.TMAX) GO TO 7 | 75 GO TO 5 | 76 7 CLOSE | 77 CLOSE | 78 STOP | 79 END |
|