MohrCoulombEngineering

Define a Mohr-Coulomb Engineering material.

Parameters

name : str: String naming the material.
color : Color: Material appearance. color=Color(r=0,g=170,b=127)
drainage : str: Material drainage type. drainage = 'drained_undrained', 'always_drained', 'non_porous'
cavitation_cutoff : str: cavitation_cutoff = 'no' or 'yes'
pcav : float: Cavitation pressure in kPa.
E : float: Young's Modulus, in MPa.
nu : float: Poisson's ratio.
c : float: Cohesion in kPa.
cu : float: Undrained shear strength.
phi_2d : float: Friction angle in 2D, in degrees.
phi_3d : float: Friction angle in 3D, in degrees.
dilation_cap : str: dilation_cap = 'no', 'volumetric', 'deviatoric'
ev_cr : float: Critical volumetric strain.
es_cr : float: Critical shear strain.
gamma_dry : float: Dry unit weight in kN/m^3
gamma_sat : float: Saturated unit weight in kN/m^3
tension_cutoff : str: Apply tension cutoff. tension_cutoff = 'yes' or 'no'
sigma_t : str: Tensile cutoff stress, in kPa.
softening : str: Apply softening. softening = 'yes' or 'no'
phi_residual : float: Residual shear strength.
e_cr1 : float: Critical plastic strain peak.
e_cr2 : float: Critical plastic strain residual.
K0 : float: Earth pressure coefficient at rest.
hydraulic_conductivity_option : str: hydraulic_conductivity_option = 'isotropic' or 'anisotropic'
K : float: Hydraulic conductivity, in m/day.
Kx : float: Hydraulic conductivity in x, in m/day.
Ky : float: Hydraulic conductivity in y, in m/day.
Kz : float: Hydraulic conductivity in z, in m/day.
hydraulic_model : str: hydraulic_model = 'basic' or 'van_genuchten'
e : float: Void ratio.
alpha : float: Model parameter related to the air entry pressure.
n : float: Model parameter related to the rate at which water is extracted from the soil once the air entry pressure has been exceeded.
Sr : float: Residual degree of saturation (may be slightly greater than 0).
Ss : float: Fraction of water filled pores at full saturation (may be slightly less than 1).
reducible_strength : str: Apply reducible strength. reducible_strength = 'yes' or 'no'

Examples

i = 7
model2d.add_rectangle([i,0],[i+1,2])
Face = model2d.select(p0=[i+0.5,1],types='face')
MCEMaterial = prj.MohrCoulombEngineering(name='MCEMaterial',
                                         color=Color(r=0,g=170,b=127),
                                         drainage = 'drained_undrained',
                                         cavitation_cutoff='yes',
                                         pcav=100,
                                         E=30,
                                         nu=0.2,
                                         c=5,
                                         cu=50,
                                         phi_2d=25,
                                         phi_3d=22.7,
                                         flow_rule='non_associated',
                                         psi=0,
                                         dilation_cap= 'deviatoric',
                                         ev_cr=0.01,
                                         es_cr=0.01,
                                         gamma_dry=18,
                                         gamma_sat=20,
                                         tension_cutoff='yes',
                                         sigma_t=0,
                                         softening='yes',
                                         phi_residual=20,
                                         e_cr1=0.01,
                                         e_cr2=0.012,
                                         K0=0.58,
                                         hydraulic_conductivity_option='anisotropic',
                                         K=0.001,
                                         Kx=0.001,
                                         Ky=0.001,
                                         Kz=0.001,
                                         hydraulic_model='van_genuchten',
                                         e=1,
                                         alpha=2,
                                         n=2,
                                         Sr=0,
                                         Ss=1,
                                         reducible_strength='yes',
                                         )
model2d.set_solid(shapes=Face,material=MCEMaterial)