import math
from actions.earth_pressure import earth_pressure
zGround= 0.0 #m z coordinate of the ground surface.
gammaSoil= 21e3 #N/m3 unit weight of soil (backfill)
H= 4 #m height of the structure.
kh= 0.086 #seismic coefficient of horizontal acceleration.
kv= kh/2.0 #seismic coefficient of vertical acceleration.
psi= math.radians(90) #back face inclination of the structure (<= PI/2)
phi= math.radians(30.0) #angle of internal friction of soil.
delta_ad= 0.0 #angle of friction soil - structure.
beta= 0.0 #slope inclination of backfill.
Kas= 1/3.0 #static earth pressure coefficient
mononobeOkabe= earth_pressure.MononobeOkabePressureDistribution(zGround, gammaSoil, H, kv, kh, psi, phi, delta_ad, beta, Kas)
overPressure= mononobeOkabe.overpressure_dry
sigmaMax= mononobeOkabe.getPressure(0.0)
sigma= mononobeOkabe.getPressure(-H/2.0)
sigmaMin= mononobeOkabe.getPressure(-H)
ratio1= (overPressure-11.3919636706e3)/11.3919636706e3
ratio2= (sigmaMax-5.69598183531e3)/5.69598183531e3
ratio3= (sigma-sigmaMax/2.0)/(sigmaMax/2.0)
ratio4= sigmaMin
'''
print('overPressure= ', overPressure/1e3, ' kN/m')
print('ratio1= ', ratio1 )
print('sigmaMax= ', sigmaMax/1e3, ' kN/m2')
quakeLoad = loadCaseManager.setCurrentLoadCase('earthquake')
xLimit = 24.5
structureHeightA = zGroundBackFill - 8.5518 #m height of the structure zone A.
structureHeightB = zGroundBackFill - 5.3593 #m height of the structure zone B.
# kh: seismic coefficient of horizontal acceleration.
# kv: seismic coefficient of vertical acceleration.
# psi: back face inclination of the structure (<= PI/2)
# delta_ad angle of friction soil - structure.
# beta slope inclination of backfill.
# Kas: static earth pressure coefficient
mononobeOkabeA = earth_pressure.MononobeOkabePressureDistribution(
zGround=zGroundBackFill,
gamma_soil=gSoil,
H=structureHeightA,
kv=0.11 / 2.0,
kh=0.11,
psi=math.radians(90),
phi=backFillSoilModel.phi,
delta_ad=0.0,
beta=0.0,
Kas=K0)
mononobeOkabeB = earth_pressure.MononobeOkabePressureDistribution(
zGround=zGroundBackFill,
gamma_soil=gSoil,
H=structureHeightB,
kv=0.11 / 2.0,
kh=0.11,
psi=math.radians(90),
phi=backFillSoilModel.phi,
delta_ad=0.0,
beta=0.0,