Struct SoilModel 

pub struct SoilModel {
    pub cohesion: f64,
    pub friction_angle: f64,
    pub dilatancy_angle: f64,
    pub young_modulus: f64,
    pub poisson_ratio: f64,
}

Soil model using the Mohr-Coulomb yield criterion.

Models elastic-perfectly-plastic soil behavior with friction and cohesion.

Fields

cohesion: f64

Cohesion c [Pa].

friction_angle: f64

Internal friction angle φ [radians].

dilatancy_angle: f64

Dilatancy angle ψ [radians].

young_modulus: f64

Young’s modulus E [Pa].

poisson_ratio: f64

Poisson’s ratio ν.

Implementations

impl SoilModel

pub fn new( cohesion: f64, friction_angle_deg: f64, dilatancy_angle_deg: f64, young_modulus: f64, poisson_ratio: f64, ) -> Self

Create a new Mohr-Coulomb soil model.

pub fn medium_sand() -> Self

Typical medium-density sand.

pub fn stiff_clay() -> Self

Typical stiff clay.

pub fn yield_function(&self, sigma1: f64, sigma3: f64) -> f64

Mohr-Coulomb yield function F(σ1, σ3).

F = (σ1 - σ3) - 2ccos(φ) - (σ1 + σ3)*sin(φ) F < 0: elastic; F = 0: on yield surface.

pub fn is_yielding(&self, sigma1: f64, sigma3: f64) -> bool

Check if the stress state is yielding.

pub fn plastic_multiplier(&self, d_eps_p: f64) -> f64

Plastic multiplier Δλ for a given incremental plastic strain magnitude.

Δλ = Δεp / (1 + sin(ψ))

pub fn critical_state_line_q(&self, p: f64) -> f64

Critical state line (CSL) deviator stress q as a function of mean stress p [Pa].

q = M_cs * p where M_cs = 6*sin(φ) / (3 - sin(φ)) for compression.

pub fn bulk_modulus(&self) -> f64

Bulk modulus K = E / (3*(1-2ν)) [Pa].

pub fn shear_modulus(&self) -> f64

Shear modulus G = E / (2*(1+ν)) [Pa].

pub fn ucs(&self) -> f64

Uniaxial compressive strength (unconfined) [Pa].

qu = 2ccos(φ) / (1 - sin(φ))

pub fn k0_jaky(&self) -> f64

Ko (coefficient of earth pressure at rest) via Jaky’s formula.

Ko = 1 - sin(φ)

Trait Implementations

impl Clone for SoilModel

fn clone(&self) -> SoilModel

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for SoilModel

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.