PILEDISP

PROGRAM PILEDISP
!**************************************************************
!** THIS PROGRAM IS USED TO CALCULATE THE PILE DISPLACEMENT  **
!** IN BOTH NORMAL AND EARTHQUAKE SITUATIONS                 **
!** AUTHOR : S.F.TUNG                                        **
!** PRACTICE PROGRAM (2017/10/11)                         01 **
!**************************************************************

IMPLICIT NONE
REAL :: N, D, BH, ALPHA, E0, KH0, KH, BETA, H0, EC 
REAL :: IPILE, FC, MPIlE, VPILE, DELTAC, DELTA, STA10, STA15
INTEGER :: CASEUA

OPEN (UNIT=10, FILE="PILEDISP.INP")
READ(10,*) CASEUA, N, D, BH, MPILE, VPILE, FC

IF (CASEUA == 1) THEN
  ALPHA = 1.0
ELSE 
  ALPHA = 2.0
END IF
STA10 = 1.00
STA15 = 1.50
E0 = 2.5*N
KH0 = E0*ALPHA/30.0
BH = 0.8*D
KH = KH0*(BH/30.0)**(-0.75)
EC = 15000.0*(FC)**0.5
IPILE = (3.1415926*(D)**4)/64.0
BETA = ((KH*D)/(4.0*EC*IPILE))**0.25
H0 = MPILE/VPILE
DELTA = (VPILE*(1.0+BETA*H0))/(2.0*EC*IPILE*BETA**3)

OPEN (UNIT=11, FILE="PILEDISP.OUT")
WRITE(11,*) "THIS PROGRAM IS USED TO BE CALED THE PILE DISPLACEMENT"
WRITE(11,*) "N VALUE OF THE LAYER, ", N
WRITE(11,*) "THE PILE DIAMETER, D =", D, "CM"
WRITE(11,*) "EFFECTIVE WIDTH OF PILE, BH =", BH, "CM"
WRITE(11,*) "EARTHEQUAKE FACTOR, ALPHA =", ALPHA
WRITE(11,*) "STIFFINESS OF PILE, E0 =", E0, "KGF/CM2"
WRITE(11,*) "HORIZONTAL STIFFINESS FACTOR OF GROUND, KH =", KH, "KGF/CM2"
WRITE(11,*) "BETA = ((KH*D)/(4.0*EC*IPILE))**0.25, BETA =", BETA
WRITE(11,*) "PILE EXPOSED LENGTH H0 =", H0, "CM"
WRITE(11,*) "STIFFINESS OF CONCRETE, EC =", EC, "KGF/CM2"
WRITE(11,*) "I VALUE OF PILE = (PI()*(D)**4)/64.0, IPILE =", IPILE, "CM4"
WRITE(11,*) "DISPLACEMENT OF PILE, DELTA =", IPILE, "CM4"

IF (CASEUA == 1) THEN
  IF(ABS(DELTA)-STA10 <= 0 ) THEN
    WRITE(11,*) "ACCORDING TO CODE, THE DISP OF PILE IS OK"
    ELSE
    WRITE(11,*) "ACCORDING TO CODE, THE DISP OF PILE IS NG"
  END IF
ELSE
  IF(ABS(DELTA)-STA15 <= 0 ) THEN
    WRITE(11,*) "ACCORDING TO CODE, THE DISP OF PILE IS OK"
    ELSE
    WRITE(11,*) "ACCORDING TO CODE, THE DISP OF PILE IS NG"
  END IF
END IF
ENDPROGRAM