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