def getMononobeOkabeDryOverpressure(self,
H,
kv,
kh,
psi=math.radians(90),
delta_ad=0.0,
beta=0.0,
Kas=None,
g=9.81):
''' Overpressure due to seismic action according to Mononobe-Okabe
Args:
:H: height of the structure.
:kv: seismic coefficient of vertical acceleration.
:kh: seismic coefficient of horizontal acceleration.
:psi: back face inclination of the structure (< PI/2)
:beta: slope inclination of backfill.
:delta_ad: angle of friction soil - structure.
:Kas: static earth pressure coefficient
'''
gamma_soil = self.rho * g
phi_d = math.atan(math.tan(self.phi) / 1.25)
if (Kas == None):
Kas = self.Ka()
return mononobe_okabe.overpressure_dry(H, gamma_soil, kv, kh, psi,
phi_d, delta_ad, beta, Kas)
def getMononobeOkabeDryOverpressure(self,H,kv,kh,psi= math.radians(90),delta_ad= 0.0,beta= 0.0,Kas= None,g= 9.81):
''' Overpressure due to seismic action according to Mononobe-Okabe
Args:
:H: height of the structure.
:kv: seismic coefficient of vertical acceleration.
:kh: seismic coefficient of horizontal acceleration.
:psi: back face inclination of the structure (< PI/2)
:beta: slope inclination of backfill.
:delta_ad: angle of friction soil - structure.
:Kas: static earth pressure coefficient
'''
gamma_soil= self.rho*g
phi_d= math.atan(math.tan(self.phi)/1.25)
if(Kas==None):
Kas= self.Ka()
return mononobe_okabe.overpressure_dry(H, gamma_soil, kv, kh, psi, phi_d, delta_ad,beta,Kas)
kh= 0.096 # seismic coefficient of horizontal acceleration.
kv= kh/2.0 # seismic coefficient of vertical acceleration.
Kas= 1.0/3.0 #static earth pressure coefficient
gamma_soil= 20e3 #unit weight of soil (backfill)
psi= math.radians(90) #back face inclination of the structure (<= PI/2)
beta= 0 #slope inclination of backfill.
delta_ad= 0 #angle of friction soil - structure.
theta1= math.radians(5.3)
theta2= math.radians(5.8)
phi= math.radians(29.0985) #angle of internal friction of soil.
phi_d= math.atan(math.tan(phi)/1.25)
L= 12.5
K_ad1= mononobe_okabe.active_earth_pressure_coefficient_Kad(psi, phi_d, theta1, delta_ad,beta)
K_ad2= mononobe_okabe.active_earth_pressure_coefficient_Kad(psi, phi_d, theta2, delta_ad,beta)
AEd1= mononobe_okabe.overpressure_dry(H, gamma_soil, kv, kh, psi, phi_d, delta_ad,beta,Kas)
AEd2= mononobe_okabe.overpressure_dry(H, gamma_soil, -kv, kh, psi, phi_d, delta_ad,beta,Kas)
K_ad1Teor= 0.487898697512
ratio1= abs(K_ad1-K_ad1Teor)/K_ad1Teor
K_ad2Teor= 0.494882910177
ratio2= abs(K_ad2-K_ad2Teor)/K_ad2Teor
AEd1Teor= 1278.02671863e3/L
ratio3= abs(AEd1-AEd1Teor)/AEd1Teor
'''
print("theta1= ", math.degrees(theta1))
print("phi_d= ", math.degrees(phi_d))
print("K_ad1=", K_ad1)
print("K_ad2=", K_ad2)
def update(self):
self.overpressure_dry = mononobe_okabe.overpressure_dry(
self.H, self.gammaSoil, self.kv, self.kh, self.psi, self.phi,
self.delta_ad, self.beta, self.Kas)
self.max_stress = 2 * self.overpressure_dry / self.H
def update(self):
self.overpressure_dry= mononobe_okabe.overpressure_dry(self.H, self.gammaSoil, self.kv, self.kh, self.psi, self.phi, self.delta_ad, self.beta, self.Kas)
self.max_stress= 2*self.overpressure_dry/self.H