Introduction

c o n s p i r e

User Guide

• Installation
• Mathematics
• Constitutive

Reference Guide

• Julia interface
• Python interface
• Rust interface

Installation

Julia

pkg> add Conspire

Python

pip install conspire

Rust

cargo add conspire

Mathematics

• Special functions

Special functions

• Lambert W function
• Langevin function
• Inverse Langevin function

Lambert W function

$$y=W_0(x)$$

Langevin function

$$\mathcal{L}(x)=\coth (x)-x^{-1}$$

Inverse Langevin function

$$x={\mathcal{L}}^{-1}(y)$$

Constitutive

• Solid constitutive models

Solid constitutive models

• Elastic constitutive models
• Hyperelastic constitutive models

Elastic constitutive models

Elastic constitutive models cannot be defined by a Helmholtz free energy density but still depend on only the deformation gradient. These constitutive models are therefore defined by a relation for some stress measure as a function of the deformation gradient. Consequently, the tangent stiffness associated with the first Piola-Kirchhoff stress is not symmetric for these constitutive models.

$${\mathcal{C}}_{iJkL}\ne {\mathcal{C}}_{kLiJ}$$

• Almansi-Hamel model

Almansi-Hamel model

{{#include ../../../../conspire.rs/src/constitutive/solid/elastic/almansi_hamel/model.md}}

Cauchy stress

$$\mathbf{\sigma }(\mathbf{F})=\frac{2\mu }{J}{\mathbf{e}}^{\prime }+\frac{\kappa }{J}\mathrm{tr}(\mathbf{e})1$$

Cauchy tangent stiffness

$${\mathcal{T}}_{ijkL}(\mathbf{F})=\frac{\mu }{J}\left[B_{jk}^{-1}{F}_{iL}^{-T}+B_{ik}^{-1}{F}_{jL}^{-T}-\frac{2}{3}{\delta }_{ij}{B}_{km}^{-1}{F}_{mL}^{-T}-2e_{ij}^{\prime }{F}_{kL}^{-T}\right]+\frac{\kappa }{J}\left[{\delta }_{ij}{B}_{km}^{-1}{F}_{mL}^{-T}-\mathrm{tr}(\mathbf{e}){\delta }_{ij}{F}_{kL}^{-T}\right]$$

Hyperelastic constitutive models

Hyperelastic constitutive models are completely defined by a Helmholtz free energy density function of the deformation gradient.

$$\mathbf{P}:\dot{\mathbf{F}}-\dot{a}(\mathbf{F})\ge 0$$ Satisfying the second law of thermodynamics (here, equivalent to extremized or zero dissipation) yields a relation for the stress.

$$\mathbf{P}=\frac{\partial a}{\partial \mathbf{F}}$$ Consequently, the tangent stiffness associated with the first Piola-Kirchhoff stress is symmetric for these constitutive models.

$${\mathcal{C}}_{iJkL}={\mathcal{C}}_{kLiJ}$$

• Arruda-Boyce model
• Fung model
• Gent model
• Mooney-Rivlin model
• Neo-Hookean model
• Saint Venant-Kirchhoff model
• Yeoh model

Arruda-Boyce model

Helmholtz free energy density

$$a(\mathbf{F})=\frac{3\mu N_b{\gamma }_0}{{\eta }_0}\left[\gamma \eta -{\gamma }_0{\eta }_0-\ln \left(\frac{{\eta }_0\sinh \eta }{\eta \sinh {\eta }_0}\right)\right]+\frac{\kappa }{2}\left[\frac{1}{2}\left(J^2-1\right)-\ln J\right]$$

Cauchy stress

$$\mathbf{\sigma }(\mathbf{F})=\frac{\mu {\gamma }_0\eta }{J\gamma {\eta }_0}{{\mathbf{B}}^{\ast }}^{\prime }+\frac{\kappa }{2}\left(J-\frac{1}{J}\right)1$$

Cauchy tangent stiffness

$$\begin{array}{rl}{\mathcal{T}}_{ijkL}(\mathbf{F})=& \frac{\mu {\gamma }_0\eta }{J^{5/3}\gamma {\eta }_0}\left({\delta }_{ik}{F}_{jL}+{\delta }_{jk}{F}_{iL}-\frac{2}{3}{\delta }_{ij}{F}_{kL}-\frac{5}{3}{B}_{ij}^{\prime }{F}_{kL}^{-T}\right)\\ & +\frac{\mu {\gamma }_0\eta }{3J^{7/3}{N}_b{\gamma }^2{\eta }_0}\left(\frac{1}{\eta {\mathcal{L}}^{\prime }(\eta )}-\frac{1}{\gamma }\right)B_{ij}^{\prime }{B}_{km}^{\prime }{F}_{mL}^{-T}+\frac{\kappa }{2}\left(J+\frac{1}{J}\right){\delta }_{ij}{F}_{kL}^{-T}\end{array}$$

Fung model

Helmholtz free energy density

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

Cauchy stress

$$\mathbf{\sigma }(\mathbf{F})=\frac{1}{J}\left[\mu +{\mu }_m\left(e^{c[\mathrm{tr}({\mathbf{B}}^{\ast })-3]}-1\right)\right]{{\mathbf{B}}^{\ast }}^{\prime }+\frac{\kappa }{2}\left(J-\frac{1}{J}\right)1$$

Cauchy tangent stiffness

$$\begin{array}{rl}{\mathcal{T}}_{ijkL}(\mathbf{F})=& \frac{1}{J^{5/3}}\left[\mu +{\mu }_m\left(e^{c[\mathrm{tr}({\mathbf{B}}^{\ast })-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}^{\prime }{F}_{kL}^{-T}\right)\\ & +\frac{2c{\mu }_m}{J^{7/3}}{e}^{c[\mathrm{tr}({\mathbf{B}}^{\ast })-3]}{B}_{ij}^{\prime }{B}_{km}^{\prime }{F}_{mL}^{-T}+\frac{\kappa }{2}\left(J+\frac{1}{J}\right){\delta }_{ij}{F}_{kL}^{-T}\end{array}$$

Gent model

Helmholtz free energy density

$$a(\mathbf{F})=-\frac{\mu J_m}{2}\ln \left[1-\frac{\mathrm{tr}({\mathbf{B}}^{\ast })-3}{J_m}\right]+\frac{\kappa }{2}\left[\frac{1}{2}\left(J^2-1\right)-\ln J\right]$$

Cauchy stress

$$\mathbf{\sigma }(\mathbf{F})=\frac{J^{-1}\mu J_m{{\mathbf{B}}^{\ast }}^{\prime }}{J_m-\mathrm{tr}({\mathbf{B}}^{\ast })+3}+\frac{\kappa }{2}\left(J-\frac{1}{J}\right)1$$

Cauchy tangent stiffness

$$\begin{array}{rl}{\mathcal{T}}_{ijkL}(\mathbf{F})=& \frac{J^{-5/3}\mu J_m}{J_m-\mathrm{tr}({\mathbf{B}}^{\ast })+3}[{\delta }_{ik}{F}_{jL}+{\delta }_{jk}{F}_{iL}-\frac{2}{3}{\delta }_{ij}{F}_{kL}+\frac{2{B_{ij}^{\ast }}^{\prime }{F}_{kL}}{J_m-\mathrm{tr}({\mathbf{B}}^{\ast })+3}\\ & -\left(\frac{5}{3}+\frac{2}{3}\frac{\mathrm{tr}({\mathbf{B}}^{\ast })}{J_m-\mathrm{tr}({\mathbf{B}}^{\ast })+3}\right)J^{2/3}{B_{ij}^{\ast }}^{\prime }{F}_{kL}^{-T}]+\frac{\kappa }{2}\left(J+\frac{1}{J}\right){\delta }_{ij}{F}_{kL}^{-T}\end{array}$$

Mooney-Rivlin model

Helmholtz free energy density

$$a(\mathbf{F})=\frac{\mu -{\mu }_m}{2}\left[\mathrm{tr}({\mathbf{B}}^{\ast })-3\right]+\frac{{\mu }_m}{2}\left[I_2({\mathbf{B}}^{\ast })-3\right]+\frac{\kappa }{2}\left[\frac{1}{2}\left(J^2-1\right)-\ln J\right]$$

Cauchy stress

$$\mathbf{\sigma }(\mathbf{F})=\frac{\mu -{\mu }_m}{J}{{\mathbf{B}}^{\ast }}^{\prime }-\frac{{\mu }_m}{J}{\left({\mathbf{B}}^{\ast -1}\right)}^{\prime }+\frac{\kappa }{2}\left(J-\frac{1}{J}\right)1$$

Cauchy tangent stiffness

$$\begin{array}{rl}{\mathcal{T}}_{ijkL}(\mathbf{F})=& \frac{\mu -{\mu }_m}{J^{5/3}}\left({\delta }_{ik}{F}_{jL}+{\delta }_{jk}{F}_{iL}-\frac{2}{3}{\delta }_{ij}{F}_{kL}-\frac{5}{3}{B}_{ij}^{\prime }{F}_{kL}^{-T}\right)+\frac{\kappa }{2}\left(J+\frac{1}{J}\right){\delta }_{ij}{F}_{kL}^{-T}\\ & -\frac{{\mu }_m}{J}\left[\frac{2}{3}{B}_{ij}^{\ast -1}{F}_{kL}^{-T}-B_{ik}^{\ast -1}{F}_{jL}^{-T}-B_{ik}^{\ast -1}{F}_{iL}^{-T}+\frac{2}{3}{\delta }_{ij}{\left(B_{km}^{\ast -1}\right)}^{\prime }{F}_{mL}^{-T}-{\left(B_{ij}^{\ast -1}\right)}^{\prime }{F}_{kL}^{-T}\right]\end{array}$$

Neo-Hookean model

{{#include ../../../../conspire.rs/src/constitutive/solid/hyperelastic/neo_hookean/model.md}}

Helmholtz free energy density

$$a(\mathbf{F})=\frac{\mu }{2}\left[\mathrm{tr}({\mathbf{B}}^{\ast })-3\right]+\frac{\kappa }{2}\left[\frac{1}{2}\left(J^2-1\right)-\ln J\right]$$

Cauchy stress

$$\mathbf{\sigma }(\mathbf{F})=\frac{\mu }{J}{{\mathbf{B}}^{\ast }}^{\prime }+\frac{\kappa }{2}\left(J-\frac{1}{J}\right)1$$

Cauchy tangent stiffness

$${\mathcal{T}}_{ijkL}(\mathbf{F})=\frac{\mu }{J^{5/3}}\left({\delta }_{ik}{F}_{jL}+{\delta }_{jk}{F}_{iL}-\frac{2}{3}{\delta }_{ij}{F}_{kL}-\frac{5}{3}{B}_{ij}^{\prime }{F}_{kL}^{-T}\right)+\frac{\kappa }{2}\left(J+\frac{1}{J}\right){\delta }_{ij}{F}_{kL}^{-T}$$

Saint Venant-Kirchhoff model

{{#include ../../../../conspire.rs/src/constitutive/solid/hyperelastic/saint_venant_kirchhoff/model.md}}

Helmholtz free energy density

$$a(\mathbf{F})=\mu \mathrm{tr}({\mathbf{E}}^2)+\frac{1}{2}\left(\kappa -\frac{2}{3}\mu \right)\mathrm{tr}(\mathbf{E}{)}^2$$

Second Piola-Kirchhoff stress

$$\mathbf{S}(\mathbf{F})=2\mu {\mathbf{E}}^{\prime }+\kappa \mathrm{tr}(\mathbf{E})1$$

Second Piola-Kirchhoff tangent stiffness

$${\mathcal{G}}_{IJkL}(\mathbf{F})=\mu {\delta }_{JL}{F}_{kI}+\mu {\delta }_{IL}{F}_{kJ}+\left(\kappa -\frac{2}{3}\mu \right){\delta }_{IJ}{F}_{kL}$$

Yeoh model

Helmholtz free energy density

$$a(\mathbf{F})=\sum _{n=1}^N\frac{{\mu }_n}{2}{\left[\mathrm{tr}({\mathbf{B}}^{\ast })-3\right]}^n+\frac{\kappa }{2}\left[\frac{1}{2}\left(J^2-1\right)-\ln J\right]$$

Cauchy stress

$$\mathbf{\sigma }(\mathbf{F})=\sum _{n=1}^N\frac{n{\mu }_n}{J}{\left[\mathrm{tr}({\mathbf{B}}^{\ast })-3\right]}^{n-1}{{\mathbf{B}}^{\ast }}^{\prime }+\frac{\kappa }{2}\left(J-\frac{1}{J}\right)1$$

Cauchy tangent stiffness

$$\begin{array}{rl}{\mathcal{T}}_{ijkL}(\mathbf{F})=& \sum _{n=1}^N\frac{n{\mu }_n}{J^{5/3}}{\left[\mathrm{tr}({\mathbf{B}}^{\ast })-3\right]}^{n-1}\left({\delta }_{ik}{F}_{jL}+{\delta }_{jk}{F}_{iL}-\frac{2}{3}{\delta }_{ij}{F}_{kL}-\frac{5}{3}{B}_{ij}^{\prime }{F}_{kL}^{-T}\right)\\ & +\sum _{n=2}^N\frac{2n(n-1){\mu }_n}{J^{7/3}}{\left[\mathrm{tr}({\mathbf{B}}^{\ast })-3\right]}^{n-2}{B}_{ij}^{\prime }{B}_{km}^{\prime }{F}_{mL}^{-T}+\frac{\kappa }{2}\left(J+\frac{1}{J}\right){\delta }_{ij}{F}_{kL}^{-T}\end{array}$$

Julia interface

Python interface

Rust interface

Contributors