Formulations

The following formulations are available in the PolyFEM list of possible formulations. The constants can be set in params. All formulations support boundary conditions. For the elasticity formulations we output:

Cauchy stress tensor

$$\sigma =\frac{1}{J}\frac{\partial \mathrm{\Psi }}{\partial F}{F}^T$$

Frist Piola Kirchhoff stress tensor (Wikipedia)

$$P=J\sigma F^{-T}$$

Second Piola Kirchhoff stress tensor (Wikipedia)

$$S=J{F}^{-1}\sigma F^{-T},$$

where $\mathrm{\Psi }$ is the energy density, $F$ the deformation gradient, and $J=\det \left(F\right)$.

Scalar

Laplacian

• Constants: none
• Description: solve for $-\mathrm{\Delta }u=f$

Bilaplacian (mixed)

• Constants: none
• Description: solve for $-{\mathrm{\Delta }}^{2}u=f$

Helmholtz

• Constants: $k$
• Description: solve for $-\mathrm{\Delta }u-k^{2}u=f$

Tensor

Linear Elasticity

• Constants: E/nu, lambda/mu
• Description: solve for $-\text{div}(\sigma [u])=f$ where

$$\sigma [u]=2\mu ϵ[u]+\lambda \text{tr}(ϵ[u])Iϵ[u]=\frac{1}{2}(\mathrm{\nabla }{u}^T+\mathrm{\nabla }u)$$

Hooke Linear Elasticity

• Constants: elasticity_tensor, E/nu, lambda/mu
• Description: solve for $-\text{div}(\sigma [u])=f$ where

$$\sigma [u]=C:ϵ[u]ϵ[u]=\frac{1}{2}(\mathrm{\nabla }{u}^T+\mathrm{\nabla }u)$$

where $C$ is the elasticity tensor

Incompressible Linear Elasticity (mixed)

• Constants: E/nu, lambda/mu
• Description: solve for

$$\begin{array}{rl}-\text{div}(2\mu ϵ[u]+pI)& =f\\ \text{div}(u)-{\lambda }^{-1}p& =0\end{array}$$

Saint Venant–Kirchoff Elasticity

• Constants: elasticity_tensor, E/nu, lambda/mu
• Description: solve for $-\text{div}(\sigma [u])=f$ where

$$\sigma [u]=C:ϵ[u]ϵ[u]=\frac{1}{2}(\mathrm{\nabla }{u}^T\mathrm{\nabla }u+\mathrm{\nabla }{u}^T+\mathrm{\nabla }u)$$

where $C$ is the elasticity tensor

NeoHookean Elasticity

• Constants: E/nu, lambda/mu
• Description: solve for $-\text{div}(\sigma [u])=f$ where

$$\sigma [u]=\mu (F[u]-F[u{]}^{-T})+\lambda \ln (\det F[u])F[u{]}^{-T}F[u]=\mathrm{\nabla }u+I$$

Mooney-Rivlin Elasticity

• Constants: c1/c2/k=1
• Description: solve for $-\text{div}(\sigma [u])=f$ where $\sigma [u]={\mathrm{\nabla }}_u\mathrm{\Psi }[u]$. The energy density $\mathrm{\Psi }$ is

$$\mathrm{\Psi }[u]=c_1(\widetilde{I_1}-d)+c_2(\widetilde{I_2}-d)+\frac{k}{2}{\ln }^2(J)$$

where $d$ is the dimension (2 or 3),

$$\begin{array}{r}F=\mathrm{\nabla }u+I,J=\det \left(F\right),\tilde{F}=\frac{1}{\sqrt[d]{J}}F,\tilde{C}=\tilde{F}{\tilde{F}}^T,\\ \widetilde{I_1}=\text{tr}\left(\tilde{C}\right),\text{and}\widetilde{I_2}=\frac{1}{2}({(\mathrm{tr}\tilde{C})}^2-\text{tr}\left({\tilde{C}}^2\right)).\end{array}$$

• Reference: FEBio documentation

Incompressible Ogden Elasticity

• Constants: c/m/k
• Description: solve for $-\text{div}(\sigma [u])=f$ where $\sigma [u]={\mathrm{\nabla }}_u\mathrm{\Psi }[u]$. The energy density $\mathrm{\Psi }$ is

$$\mathrm{\Psi }[u]=\sum _{i=1}^N\frac{c_i}{m_i^2}(\sum _{j=1}^d{\tilde{\lambda }}_j^{m_i}-d)+\frac{1}{2}K\ln {\left(J\right)}^2$$

where $N$, the number of terms, is dictated by the number of coefficients given, $d$ is the dimension (2 or 3), $J=\det \left(F\right)$ where $F=\mathrm{\nabla }u+I$, and ${\tilde{\lambda }}_j=J^{-\frac{1}{d}}{\lambda }_j$ are the eigenvalues of $\tilde{F}$ (same as in the Mooney-Rivlin Elasticity).

• Reference: FEBio documentation

Unconstrained Ogden Elasticity

• Constants: mus/alphas/Ds
• Description: solve for $-\text{div}(\sigma [u])=f$ where $\sigma [u]={\mathrm{\nabla }}_u\mathrm{\Psi }[u]$. The energy density $\mathrm{\Psi }$ is

$$\mathrm{\Psi }[u]=\sum _{i=1}^N\frac{2{\mu }_i}{{\alpha }_i^2}(\sum _{j=1}^d{\tilde{\lambda }}_j^{{\alpha }_i}-d)+\sum _{i=1}^N\frac{{(J-1)}^{2i}}{D_i}$$

where $N$, the number of terms, is dictated by the number of coefficients given, $d$ is the dimension (2 or 3), $J=\det \left(F\right)$ where $F=\mathrm{\nabla }u+I$, and ${\tilde{\lambda }}_j=J^{-\frac{1}{d}}{\lambda }_j$ where ${\lambda }_j$ are the eigenvalues of $F.$

• Reference: ABAQUS documentation

Viscous Damping

• Constants: phi/psi
• Description: an extra energy that represents dissipation, adding to the elastic energy in transient problems

$$R(F,\dot{F})=\psi \Vert \dot{E}(F,\dot{F}){\Vert }^2+\frac{\varphi }{2}{\text{tr}}^2\dot{E}(F,\dot{F})$$

where $F[u]=\mathrm{\nabla }u+I$ and $E[u]=\frac{1}{2}(F^{T}F-I)$.

The above corresponds to the viscous Piola-Kirchhoff stress

$$P=F(2\psi \dot{E}+\varphi \text{tr}(\dot{E})I)={\mathrm{\nabla }}_{2}R(F,\dot{F}).$$

Stokes (mixed)

• Constants: viscosity ($\nu$)
• Description: solve for

$$\begin{array}{rl}-\nu \mathrm{\Delta }u+\mathrm{\nabla }p& =f\\ -\text{div}(u)& =0\end{array}$$

Navier–Stokes (mixed)

• Constants: viscosity ($\nu$)
• Description: solve for

$$\begin{array}{rl}u\cdot \mathrm{\nabla }u-\nu \mathrm{\Delta }u+\mathrm{\nabla }p& =f\\ -\text{div}(u)& =0\end{array}$$

Implementing New Formulations

Todo

Describe how to implement a new formulation in C++.

---

Last update: 2023-10-03