Struct Fung 

pub struct Fung {
    pub bulk_modulus: Scalar,
    pub shear_modulus: Scalar,
    pub extra_modulus: Scalar,
    pub exponent: Scalar,
}

The Fung hyperelastic solid constitutive model.¹

Parameters

• The bulk modulus $\kappa$.
• The shear modulus $\mu$.
• The extra modulus $\mu_m$.
• The exponent $c$.

External variables

• The deformation gradient $\mathbf{F}$.

Internal variables

• None.

Notes

• The Fung model reduces to the Neo-Hookean model when $\mu_m\to 0$ or $c\to 0$.

---

1. Y.C. Fung, Am. J. Physiol. 213, 1532 (1967). ↩

Fields

bulk_modulus: Scalar

The bulk modulus $\kappa$.

shear_modulus: Scalar

The shear modulus $\mu$.

extra_modulus: Scalar

The extra modulus $\mu_m$.

exponent: Scalar

The exponent $c$.

Implementations

impl Fung

pub fn extra_modulus(&self) -> Scalar

Returns the extra modulus.

pub fn exponent(&self) -> Scalar

Returns the exponent.

Trait Implementations

impl Clone for Fung

fn clone(&self) -> Fung

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for Fung

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

Formats the value using the given formatter.

impl Elastic for Fung

fn cauchy_stress( &self, deformation_gradient: &DeformationGradient, ) -> Result<CauchyStress, ConstitutiveError>

\boldsymbol{\sigma}(\mathbf{F}) = \frac{1}{J}\left[\mu + \mu_m\left(e^{c[\mathrm{tr}(\mathbf{B}^* ) - 3]} - 1\right)\right]{\mathbf{B}^* }' + \frac{\kappa}{2}\left(J - \frac{1}{J}\right)\mathbf{1}

fn cauchy_tangent_stiffness( &self, deformation_gradient: &DeformationGradient, ) -> Result<CauchyTangentStiffness, ConstitutiveError>

\begin{aligned}
\mathcal{T}_{ijkL}(\mathbf{F}) =
\,& \frac{1}{J^{5/3}}\left[\mu + \mu_m\left(e^{c[\mathrm{tr}(\mathbf{B}^* ) - 3]} - 1\right)\right]\left(\delta_{ik}F_{jL} + \delta_{jk}F_{iL} - \frac{2}{3}\,\delta_{ij}F_{kL} - \frac{5}{3} \, B_{ij}'F_{kL}^{-T} \right)
\\
&+ \frac{2c\mu_m}{J^{7/3}}\,e^{c[\mathrm{tr}(\mathbf{B}^* ) - 3]}B_{ij}'B_{km}'F_{mL}^{-T} + \frac{\kappa}{2} \left(J + \frac{1}{J}\right)\delta_{ij}F_{kL}^{-T}
\end{aligned}

fn first_piola_kirchhoff_stress( &self, deformation_gradient: &DeformationGradient, ) -> Result<FirstPiolaKirchhoffStress, ConstitutiveError>

Calculates and returns the first Piola-Kirchhoff stress.

fn first_piola_kirchhoff_tangent_stiffness( &self, deformation_gradient: &DeformationGradient, ) -> Result<FirstPiolaKirchhoffTangentStiffness, ConstitutiveError>

Calculates and returns the tangent stiffness associated with the first Piola-Kirchhoff stress.

fn second_piola_kirchhoff_stress( &self, deformation_gradient: &DeformationGradient, ) -> Result<SecondPiolaKirchhoffStress, ConstitutiveError>

Calculates and returns the second Piola-Kirchhoff stress.

fn second_piola_kirchhoff_tangent_stiffness( &self, deformation_gradient: &DeformationGradient, ) -> Result<SecondPiolaKirchhoffTangentStiffness, ConstitutiveError>

Calculates and returns the tangent stiffness associated with the second Piola-Kirchhoff stress.

impl Hyperelastic for Fung

fn helmholtz_free_energy_density( &self, deformation_gradient: &DeformationGradient, ) -> Result<Scalar, ConstitutiveError>

a(\mathbf{F}) = \frac{\mu - \mu_m}{2}\left[\mathrm{tr}(\mathbf{B}^* ) - 3\right] + \frac{\mu_m}{2c}\left(e^{c[\mathrm{tr}(\mathbf{B}^* ) - 3]} - 1\right)

impl Solid for Fung

fn bulk_modulus(&self) -> Scalar

Returns the bulk modulus.

fn shear_modulus(&self) -> Scalar

Returns the shear modulus.

fn jacobian<const I: usize, const J: usize>( &self, deformation_gradient: &DeformationGradientGeneral<I, J>, ) -> Result<Scalar, ConstitutiveError>

Calculates and returns the Jacobian.