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
Algorithm 1: MAIN program.