# MohrCoulombEngineering

Define a Mohr-Coulomb Engineering material.

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'

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)