russell_tensor 3.2.1

Tensor analysis, calculus, and functions for continuum mechanics
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

use russell_tensor::{StrError, Tensor2, t2_matmulx};

fn main() -> Result<(), StrError> {
    // Deformation gradient and Jacobian
    #[rustfmt::skip]
    let ff = Tensor2::<9>::from_std_matrix(&[[1.2, 0.1, 0.0],
        [0.2, 0.9, 0.0],
        [0.0, 0.0, 1.0],
    ])?;
    let jj = ff.determinant();

    // Second Piola-Kirchhoff stress
    #[rustfmt::skip]
    let ss = Tensor2::<6>::from_std_matrix(&[
        [10.0,  2.0, 0.0],
        [ 2.0, 20.0, 0.0],
        [ 0.0,  0.0, 5.0],
    ])?;

    // 1. Push-Forward: Material to Spatial
    let mut sigma = Tensor2::<6>::new();
    t2_matmulx(&mut sigma, 1.0 / jj, &ff, true, &ss)?;

    // 2. Pull-Back: Spatial to Material
    let mut ff_inv = Tensor2::<9>::new();
    ff.inverse(&mut ff_inv, 1e-10).unwrap();

    // Check
    let mut ss_back = Tensor2::<6>::new();
    t2_matmulx(&mut ss_back, jj, &ff_inv, false, &sigma)?;

    Ok(())
}