conspire 0.6.0

The Rust interface to conspire.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57

//! Solid constitutive models.

pub mod elastic;
pub mod elastic_hyperviscous;
pub mod elastic_plastic;
pub mod elastic_viscoplastic;
pub mod hyperelastic;
pub mod hyperelastic_viscoplastic;
pub mod hyperviscoelastic;
pub mod thermoelastic;
pub mod thermohyperelastic;
pub mod viscoelastic;

const TWO_THIRDS: Scalar = 2.0 / 3.0;
const FIVE_THIRDS: Scalar = 5.0 / 3.0;

use crate::{
    constitutive::{Constitutive, ConstitutiveError},
    math::{
        ContractFirstSecondIndicesWithSecondIndicesOf, ContractSecondIndexWithFirstIndexOf,
        IDENTITY, IDENTITY_00, Rank2, Tensor, TensorArray, ZERO_10,
    },
    mechanics::{
        CauchyRateTangentStiffness, CauchyStress, CauchyTangentStiffness, Deformation,
        DeformationError, DeformationGradient, DeformationGradientGeneral, DeformationGradientRate,
        DeformationGradientRates, DeformationGradients, FirstPiolaKirchhoffRateTangentStiffness,
        FirstPiolaKirchhoffStress, FirstPiolaKirchhoffTangentStiffness, Scalar,
        SecondPiolaKirchhoffRateTangentStiffness, SecondPiolaKirchhoffStress,
        SecondPiolaKirchhoffTangentStiffness, Times,
    },
};
use std::fmt::Debug;

impl<C> Constitutive for C where C: Solid {}

/// Required methods for solid constitutive models.
pub trait Solid
where
    Self: Constitutive,
{
    /// Returns the bulk modulus.
    fn bulk_modulus(&self) -> Scalar;
    /// Returns the shear modulus.
    fn shear_modulus(&self) -> Scalar;
    /// Calculates and returns the Jacobian.
    fn jacobian<const I: usize, const J: usize>(
        &self,
        deformation_gradient: &DeformationGradientGeneral<I, J>,
    ) -> Result<Scalar, ConstitutiveError> {
        match deformation_gradient.jacobian() {
            Err(DeformationError::InvalidJacobian(jacobian)) => Err(
                ConstitutiveError::InvalidJacobian(jacobian, format!("{self:?}")),
            ),
            Ok(jacobian) => Ok(jacobian),
        }
    }
}