Struct MuscleModel 

pub struct MuscleModel {
    pub f_max: f64,
    pub l_opt: f64,
    pub v_max: f64,
    pub pennation_angle: f64,
}

Muscle model using Hill’s three-element (active + passive + series) model.

Fields

f_max: f64

Maximum isometric force [N].

l_opt: f64

Optimal fiber length [m].

v_max: f64

Maximum shortening velocity [m/s].

pennation_angle: f64

Pennation angle at optimal length [radians].

Implementations

impl MuscleModel

pub fn new(f_max: f64, l_opt: f64, v_max: f64, pennation_angle: f64) -> Self

Create a new muscle model.

pub fn tibialis_anterior() -> Self

Typical human tibialis anterior muscle.

pub fn quadriceps() -> Self

Typical human quadriceps (rectus femoris).

pub fn active_force_length(&self, l_normalized: f64) -> f64

Active force-length relationship f_l(l̃) where l̃ = l/l_opt.

Gaussian: f_l = exp(-((l̃ - 1)/ω)^2), ω ≈ 0.45

pub fn passive_force_length(&self, l_normalized: f64) -> f64

Passive force-length relationship f_p(l̃).

Exponential rise for l̃ > 1: f_p = (exp(k_p*(l̃-1)/ε_0) - 1) / (exp(k_p) - 1)

pub fn force_velocity(&self, v_normalized: f64) -> f64

Force-velocity relationship f_v(ṽ) where ṽ = v/v_max (Hill’s equation).

Concentric: f_v = (1 - ṽ)/(1 + ṽ/a_v), a_v ≈ 0.25 Eccentric: f_v = (1.8 - 0.8*(1+ṽ)/(1-7.56*ṽ))

pub fn muscle_force(&self, l: f64, v: f64, activation: f64) -> f64

Compute muscle force [N] given fiber length [m], velocity [m/s], and activation a ∈ [0,1].

pub fn max_power(&self) -> f64

Maximum power output [W].

Trait Implementations

impl Clone for MuscleModel

fn clone(&self) -> MuscleModel

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for MuscleModel

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

Formats the value using the given formatter.