conspire/constitutive/solid/elastic/almansi_hamel/
mod.rs

1#[cfg(test)]
2mod test;
3
4use crate::{
5    constitutive::{
6        Constitutive, ConstitutiveError, Parameters,
7        solid::{Solid, TWO_THIRDS, elastic::Elastic},
8    },
9    math::{IDENTITY, Rank2},
10    mechanics::{CauchyStress, CauchyTangentStiffness, DeformationGradient, Scalar},
11};
12
13#[doc = include_str!("doc.md")]
14#[derive(Debug)]
15pub struct AlmansiHamel<P> {
16    parameters: P,
17}
18
19impl<P> Constitutive<P> for AlmansiHamel<P>
20where
21    P: Parameters,
22{
23    fn new(parameters: P) -> Self {
24        Self { parameters }
25    }
26}
27
28impl<P> Solid for AlmansiHamel<P>
29where
30    P: Parameters,
31{
32    fn bulk_modulus(&self) -> &Scalar {
33        self.parameters.get(0)
34    }
35    fn shear_modulus(&self) -> &Scalar {
36        self.parameters.get(1)
37    }
38}
39
40impl<P> Elastic for AlmansiHamel<P>
41where
42    P: Parameters,
43{
44    #[doc = include_str!("cauchy_stress.md")]
45    fn cauchy_stress(
46        &self,
47        deformation_gradient: &DeformationGradient,
48    ) -> Result<CauchyStress, ConstitutiveError> {
49        let jacobian = self.jacobian(deformation_gradient)?;
50        let inverse_deformation_gradient = deformation_gradient.inverse();
51        let strain = (IDENTITY
52            - inverse_deformation_gradient.transpose() * &inverse_deformation_gradient)
53            * 0.5;
54        let (deviatoric_strain, strain_trace) = strain.deviatoric_and_trace();
55        Ok(deviatoric_strain * (2.0 * self.shear_modulus() / jacobian)
56            + IDENTITY * (self.bulk_modulus() * strain_trace / jacobian))
57    }
58    #[doc = include_str!("cauchy_tangent_stiffness.md")]
59    fn cauchy_tangent_stiffness(
60        &self,
61        deformation_gradient: &DeformationGradient,
62    ) -> Result<CauchyTangentStiffness, ConstitutiveError> {
63        let jacobian = self.jacobian(deformation_gradient)?;
64        let inverse_transpose_deformation_gradient = deformation_gradient.inverse_transpose();
65        let inverse_left_cauchy_green_deformation = &inverse_transpose_deformation_gradient
66            * inverse_transpose_deformation_gradient.transpose();
67        let scaled_inverse_left_cauchy_green_deformation =
68            &inverse_left_cauchy_green_deformation * (self.shear_modulus() / jacobian);
69        let strain = (IDENTITY - &inverse_left_cauchy_green_deformation) * 0.5;
70        let (deviatoric_strain, strain_trace) = strain.deviatoric_and_trace();
71        Ok((CauchyTangentStiffness::dyad_il_jk(
72            &inverse_transpose_deformation_gradient,
73            &scaled_inverse_left_cauchy_green_deformation,
74        ) + CauchyTangentStiffness::dyad_ik_jl(
75            &scaled_inverse_left_cauchy_green_deformation,
76            &inverse_transpose_deformation_gradient,
77        )) + CauchyTangentStiffness::dyad_ij_kl(
78            &IDENTITY,
79            &(inverse_left_cauchy_green_deformation
80                * &inverse_transpose_deformation_gradient
81                * ((self.bulk_modulus() - self.shear_modulus() * TWO_THIRDS) / jacobian)),
82        ) - CauchyTangentStiffness::dyad_ij_kl(
83            &(deviatoric_strain * (2.0 * self.shear_modulus() / jacobian)
84                + IDENTITY * (self.bulk_modulus() * strain_trace / jacobian)),
85            &inverse_transpose_deformation_gradient,
86        ))
87    }
88}
conspire/constitutive/solid/elastic/almansi_hamel/mod.rs

conspire/constitutive/solid/elastic/almansi_hamel/
mod.rs
