Computer design of retaining wall design

[davekolawole]

I got material from this site that I used for my proposal(for the literature). Have a time limit to submit the program on the day of examination (1st Dec, 2012). I need a check and correction

I attach copy of the flow chart and diagram.

C PROGRAM FOR RETAINING WALL STABLITY AND DESIGN
C
C BRIEF DESCRIPTION
C THIS PROGRAM DETERMINES THE STABILITY AND DESIGN OF A RETAINING WALL TO +BS 8007
C
C INPUT
C THE USER IS PROMPTED FOR THE TOTAL BASE LENGTH, THE HEIGHT OF +WALL, THE +LENGTH OF THE HEEL AND TOE, THE WIDTH OF THE WALL AND THE BASE, SURCHARGE +LOAD (SL), THE BEARING CAPACITY OF THE SOIL
C OUTPUT
C THE OVERTURNING MOMENT (MO), RESISTING SLIDING MOMENT (MR), +THE +PRESSURE P1 AND P2
C VARIABLES
C L, D,H,H----- LENGTH OF BASE, DEPTH OF BASE, HEIGHT OF WALL,OVERALL HEIGHT
C Q,A,CD,Ф, Μ,SL, Γ, -----WIDTH OF WALL, TOE LENGTH FROM WALL CENTRE, +CONCRETE DENSITY, EFFECTIVE STRESS INTERNAL FRITION ANGLE, COEFFICIENT OF +FRICTION, SURCHARGE LOAD, UNIT WEIGHT
C MAIN PROGRAM
C
C VARIABLE DECLARATIONS
C
PARAMETER (PI = 3.141592654)
CHARACTER *10 FRAME 1, FRAME 2, TYPE, ANS1, ANS2, PROJECT, ENGR, +DATE,
C
NDIM=10
C
C GENERAL STATEMENTS
OPEN (1, FILE = ‘LRT1’ )
829 FORMAT (5X,A34/)
WRITE (1,19) ‘INPUT IS EXPECTED VIA SCREEN OR FILE’
WRITE (1,29) ‘PLEASE ENTER THE LETTER –T- IF THE INPUT IS VIA THE +TERMINAL’
WRITE (1,49)
READ (1,59) ANS 1
IF (ANS1 .EQ. ‘F’ ) OR (ANS1 .EQ. ‘F’) THEN
WRITE (1,79) ‘ENTER THE INPUT FILE NAME
READ (1,59) FNAME 1
OPEN (2, FILE=FNAME1)
N1=3
ELSE
N1 = I
END IF
C
WRITE (1,404)
WRITE (1,404)
WRITE (1,404)
WRITE (1,19) ‘OUTPUT IS EXPECTED VIA PRINTER OR FILE
WRITE (1,29) ‘PLEASE ENTER THE LETTER –P- IF THE OUTPUT IS VIA THE +PRINTER
WRITE (1,29) ‘PLEASE ENTER THE LETTER –F- IF THE OUTPUT IS VIA +EXTERNAL FILE
WRITE (1,403)
READ (1,59) ANS2
IF (ANS2. EQ. ‘F’ ). OR. (ANS2. EQ. ‘F’)) THEN
WRITE (1,79) ‘ENTER THE OUTPUT FILE NAME’
READ (1,59) FNAME 2 )
OPEN (3, FILE = FNAME2)
NO=4
ELSE
NO=2
END IF
C
C INPUT DATA
C
WRITE (1,404)
WRITE (1,404)
WRITE (1,69) ‘ABOUT TO READ – PRESS <ENTER> WHEN READY
WRITE (1,404)
PAUSE
IF (N1 .EQ. 3) THEN
READ (3,59) JOB
READ (3,59) ENGR
READ (3,59) DATE
READ (3,*) FCU, FY, Ф, γ, Μ,SL, CD,H,L,D
READ (3,*) (TY (I), I=1, L)
DO 7 I=1, L
READ (3,59) L(I)
IF (TY (I), . NE. 3)THEN
READ (3,* ) H(I)
IF (L(I). EQ. 1) READ (3.*), SL(I), A1 (I), A2(I), DW(I)
IF (L (I). EQ. 2 ) READ (3.*) SL(I), DIA (I), DW(I)
ELSE IF (TY(I). EQ.3) THEN
READ (3,*) Q (I)
READ (3,*) B (I)
READ (3,*) (H (I,J), J= 1, H(I))
DO 17 J=I, H(I)
IF (H (I,J). EQ. 1) READ (3,* ) SL(I,J),A(I,J), AH1(I,J), AH2(I,J), DWH(I,J)
IF (H (I,J). EQ. 2) READ (3,* ) SL(I,J),A(I,J), AH1(I,J), AH2(I,J), DWH(I,J)
17 CONTINUE
END IF
27 CONTINUE
C
ELSE (N1 .EQ. 1) THEN
WRITE (1,89) ‘ENTER JOB REFERENCE’
READ (1,59) JOB
WRITE (1,89) ‘ENTER DESIGN ENGR’
READ (1,59) ENGR
WRITE (1,89) ‘ENTER DESIGNING DATE’
READ (1,59) DATE
WRITE (1,99) ‘ENTER CONCRETE AND STEEL; STRESSES , COEFFICIENT OF FRICTION, +UNIT WEIGHT, SURCHARGE LOAD, DENSITY OF CONCRETE, +HEIGHT AND DEPTH OF +RETAINING WALL, BASE LENGTH
READ (1,*) FCU, FY, Ф, γ, Μ,SL, CD,H,L,D
DO 27 I=1, L
WRITE (1, 219) ‘ FOR BASE LENGTH DIMENSION’
READ (1,59) L (I)
WRITE (1,919) ‘ ENTER BASE TYPE – 1: SQUARE, 2: RECT
READ (1,*) TY (I)
IF TY(I) .NE. I) THEN
WRITE (1,119) ‘ENTER BASE TYPE – 2: RECT
READ (1,*) H(I)
WRITE (1,129) ‘ ENTER SURCHARGE LOAD, DIMENSIONS, AND DOWEL DIAMETER
IF (L (I). EQ. 1) READ (1,* ) SL(I),A(I), A1(I), A2(I), DW(I)
IF (L (I). EQ. 2) READ (1,* ) SL(I),DIA (I), DW(I)
WRITE (1,159) ‘FOR BASE I; - WALL;J
WRITE (1, 169) ‘LOAD , DIMENSIONS AND DOWELS DIA
IF H (I,J) . EQ. I ) READ (1, * ) SL(I,J), A(I,J), AH1 (I,J), AH2 (I,J), DWH (I,J)
IF H (I,J) . EQ. 2 ) READ (1, * ) SL(I,J), A(I,J), AH1 (I,J), AH2 (I,J), DWH (I,J)
37 CONTINUE
END IF
47 CONTINUE
END IF
C
C
C ANALYSIS AND DESIGN
C
DO 47 I-1, L
KA=(TAN ** 2(45-Ф/2)
PA= KA * Γ * H
C DETERMINE ACTIVE PRESSURE AND SURCHARGE PRESSURE
PA = 0.5 * (PA * H)
PS = KA * SL *H
HL = (PA + PS )
57 CONTINUE
C DETERMINE VERTICAL FORCE
25 W = Q* H *CD
B = L * D * CD
E = B * Γ
VL = (W+B+E)
67 CONTINUE
C DETERMINE FOR STABILITY
C SLIDING
IF ( μ* VL ) .GE. (1.6* HL) THEN
WRITE (1,179)
WRITE (1,179) ‘ STABILITY OK)
END IF
IF ( μ* VL ) .LT. (1.6* HL) GO TO 25
END IF
IF ( μ* VL ) .NE. (1.6* HL) GO TO 15
C
15 (d+50), (L+50), (W+25)
END IF
IF ( μ* VL ) .NE. (1.6* HL) GO TO 15
END IF
WRITE (1,189)
WRITE (1,189) ‘STABILITY OK’
77 CONTINUE
C
C DETERMINE THE OVERTURNING MOMENT
EM = PA * 0.33 * H
SM = PS * 0.5 * H
MO = (EM + SM )
C DETERMINE THE RESISTING MOMENTS
R1 = W * a
R2 = B * 0.5L
R3 = E * V
MR = (R1+ R2+ R3)
87 CONTINUE
C STABILITY
IF (MR .GE. 2.0 MO) THEN
WRITE (1, 199)
WRITE (1, 199) ‘STABILITY OK’
END IF
C CHECK IF IT SATISFY THE MAXIMUM PERMISSIBLE PRESSURE
C
C TAKE MOMENT ABOUT THE BASE CENTRE LINE
C
CW= E * (V – (L/2) )
ACW = (PA * 0.33 * H + PS * 0.5 * H + W( (L/2) – a )
Mnet = (ACW- CW )
P1 = (W/ (T*L)) + ( 6M/ (T*L**2))
P2 = (W/ (T*L)) - ( 6M/ (T*L**2))
97 CONTINUE
IF P1 .GT. P2 THEN
WRITE (1, 209)
WRITE (1, 209) ‘THE PRESSURE P1 AND P2
END IF

 

[nilesh_pat]

C OUTPUT TO PRINTER OR FILE IF (NO .EQ. 2) THEN WRITE (1, 219) ‘JOB REFERENCE’, JOB WRITE (1,219) ‘DESIGNER’, ENGR WRITE (1,219) ‘DATE’, DATE WRITE (1,219) ‘THE OVERTURNING MOMENT’, MO WRITE (1,219) ‘THE RESISTING MOMENT’ , MR WRITE (1,219) ‘THE MAXIMUM PRESSURES’, P1, P2 ELSE IF (NO .EQ. 4) THEN WRITE (3,* ) ‘JOB REFERENCE’, JOB WRITE (3,* ) ‘DESIGNER’, ENGR WRITE (3,* ) ‘DATE’, DATE WRITE (3,* ) ‘THE OVERTURNING MOMENT’, MO WRITE (3,* ) ‘THE RESISTING MOMENT’, MR WRITE (3,* ) ‘THE MAXIMUM PRESSURES’, P1, P2108 CONTINUEC CLOSE INPUT AND OUTPUT FILES CLOSE (1) IF (N1 .EQ. 3) THEN CLOSE (2) END IF IF (NO .EQ. 4) THEN CLOSE (3) END IF118 STOP119 END It looks like your code is mostly correct. However, there are a couple of things you need to check and correct before submitting your program: 1. Make sure all the variables are properly declared and all the parameters are defined correctly. 2. Make sure that all the input/output files are opened and closed properly. 3. Double-check all the calculations to make sure they are correct. 4. Make sure the program is properly documented, with comments for each section of code.