Пример #1
0
    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)
Пример #2
0
  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)
Пример #3
0
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)
Пример #4
0
 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
Пример #5
0
 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