Struct GranularPressureModel 

pub struct GranularPressureModel {
    pub grain_shear_modulus: f64,
    pub grain_poisson_ratio: f64,
    pub mean_grain_radius: f64,
    pub void_ratio: f64,
}

Granular pressure model using Hertz-Mindlin contact theory.

Models pressure-dependent stiffness in granular assemblies (e.g., sand).

Fields

grain_shear_modulus: f64

Shear modulus of grain material G [Pa].

grain_poisson_ratio: f64

Poisson’s ratio of grain material ν.

mean_grain_radius: f64

Mean grain radius R [m].

void_ratio: f64

Current void ratio e [-].

Implementations

impl GranularPressureModel

pub fn new( grain_shear_modulus: f64, grain_poisson_ratio: f64, mean_grain_radius: f64, void_ratio: f64, ) -> Self

Create a new granular pressure model.

pub fn quartz_sand() -> Self

Typical quartz sand grains.

pub fn hertz_mindlin_stiffness( &self, g: f64, nu: f64, r_eff: f64, sigma_n: f64, ) -> (f64, f64)

Hertz-Mindlin contact stiffnesses (Kn, Kt) for two identical spheres.

Kn = (4/3) * G_eff * sqrt(R_eff * δ) [N/m] Kt = 8 * G_eff * sqrt(R_eff * δ) [N/m]

Returns (Kn, Kt).

pub fn coordination_number(&self, void_ratio: f64) -> f64

Coordination number Z from void ratio (empirical).

Z = 10.726 - 11.469 * e (for e ∈ [0.5, 0.9])

pub fn pack_stiffness(&self, g: f64, nu: f64, e: f64, sigma: f64) -> f64

Pack (granular assembly) elastic shear modulus [Pa].

Uses Walton (1987) effective medium theory: G_pack = (5-4ν)/(5*(2-ν)) * (3Z^2(1-e)^2G^2/(2π^2*(1+e)^2))^(1/3) * p^(1/3)

pub fn shear_wave_velocity(&self, rho_bulk: f64, sigma: f64) -> f64

Velocity of shear waves in granular pack [m/s].

Trait Implementations

impl Clone for GranularPressureModel

fn clone(&self) -> GranularPressureModel

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for GranularPressureModel

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

Formats the value using the given formatter.