# Ogden Materials

Hyperelastic materials can be used to model the isotropic, nonlinear elastic behavior of rubber, polymers, and similar materials. These materials are nearly incompressible in their behavior and can be stretched to very large strains.

In Radioss, material laws LAW42, LAW62, LAW69, LAW82, and
LAW88 utilize different strain energy density functions of the Ogden material model
^{1} to model hyperelastic materials. ^{2}

## Material Definition

- $${\epsilon}_{1}=\frac{\text{\Delta}l}{{l}_{01}}$$
- Engineering strain (also called nominal strain) in direction 1
- $${\lambda}_{1}=\frac{{l}_{1}}{{l}_{01}}$$
- Stretch in direction 1

- Uniaxial test:
${\lambda}_{1}=\lambda \text{\hspace{0.17em}}$ and ${\lambda}_{2}{}^{2}={\lambda}_{3}{}^{2}=\frac{1}{\lambda}$

- Biaxial test:
${\lambda}_{1}={\lambda}_{2}=\lambda \text{\hspace{0.17em}}$ and ${\lambda}_{3}={\lambda}^{-2}$

- Planar (shear) test:
${\lambda}_{1}=\lambda \text{\hspace{0.17em}};\text{\hspace{0.17em}}{\lambda}_{3}=1$ and ${\lambda}_{2}={\lambda}^{-1}\text{}$

## /MAT/LAW42 (Ogden)

This material model defines a hyperelastic, viscous, and incompressible material specified using the Ogden, Neo-Hookean, or Mooney-Rivlin material models. This law is generally used to model incompressible rubbers, polymers, foams, and elastomers. This material can be used with shell and solid elements.

- $$W$$
- Strain energy density
- ${\lambda}_{i}$
`i`^{th}principal engineering stretch- $$J$$
- Relative volume defined as: $$J={\lambda}_{1}\cdot {\lambda}_{2}\cdot {\lambda}_{3}=\frac{{\rho}_{0}}{\rho}$$
- ${\overline{\lambda}}_{i}={J}^{-\frac{1}{3}}{\lambda}_{i}$
- Deviatoric stretch
- ${\alpha}_{p}$ and ${\mu}_{p}$
- Material constants coefficient pairs.

Where, $\nu $ is the Poisson's ratio and is only used for computing the bulk modulus.

### Material Parameters

In general, the Ogden model can be used for strains up to 700%. The number of terms material pairs, ${\alpha}_{p}$ and ${\mu}_{p}$, needed depends on the range of experimental data that is fit and curve fitting accuracy desired. In practice, 3 material pairs fit most data. If the material pairs are not known for a particular material, then a curve fit of uniaxial test data can be done in Radioss using LAW69 or via separate fitting software.

### Neo-Hookean Model

- ${I}_{1}$
- The first invariants of the right Cauchy-Green Tensor
- ${C}_{10}$
- Material constant

This representation can be derived from the LAW42 Ogden strain energy density function when:

${\mu}_{1}=2\cdot {C}_{10}$; ${\alpha}_{1}=2$, and ${\mu}_{2}={\alpha}_{2}=0$

The Neo-Hookean model is a simple model that is typically only accurate for strains less than 20%.

### Mooney-Rivlin Model

- ${I}_{1}$ and ${I}_{2}$
- The first and second invariants of the right Cauchy-Green Tensor
- ${C}_{10}$ and ${C}_{01}$
- Material constants

This representation can be derived from the LAW42 Ogden strain energy density function when:

${\mu}_{1}=2\cdot {C}_{10}$, ${\mu}_{2}=-2\cdot {C}_{01}$, ${\alpha}_{1}=2$, and ${\alpha}_{2}=-2$

Mooney-Rivlin constants are available from a material supplier or testing company. If they are not available, then a curve fit of uniaxial test data can be done in Radioss using LAW69 or via separate fitting software. The Mooney-Rivlin material law is accurate for strains up to 100%.

### Poisson's Ratio and Material Incompressibility

The effect of different Poisson’s ratio input can be seen in Figure 2. The largest difference in the results is at higher amounts of strain. The results will match the test data better when $\nu =0.4997$ but this results in a time step that is 4 times lower than $\nu =0.495$. Thus, to balance the computation time and accuracy it is recommended to use $\nu =0.495$ for incompressible rubber material.

Higher values of the Poisson’s ratio may lead to a very small time step or divergence for explicit simulations.

- $$K$$
- Bulk modulus
- $$J=\frac{V}{{V}_{0}}=\frac{m{\rho}_{0}}{m\rho}=\frac{{\rho}_{0}}{\rho}$$
- Relative volume which simplifies to relative density if mass is constant
- $${\mathrm{f}}_{blk}\left(J\right)$$
- Bulk coefficient scale factor versus relative volume function
- $$Fscal{e}_{blk}$$
- Abscissa scale factor for function $${\mathrm{f}}_{blk}\left(J\right)$$

`fct_ID`

_{blk}input function which allows the scaling of the bulk coefficient value as a function of $$J$$. By default, there is no scaling and; thus, if the function identifier is zero and the value of the bulk scaling function is equal to 1. It is advisable to output (/ANIM/BRICK/DENS) and review the material density of LAW42 components to make sure that the density variation is small, i.e. the value of $$J$$ is close to 1 and the material is incompressible.

### Viscous (Rate) Effects

The $${G}_{i}$$ and ${\tau}_{i}$ values must be positive.

## /MAT/LAW62 (VISC_HYP)

A hyper visco-elastic material law in Radioss that can be used to model polymers and elastomers.

- $W$
- Strain energy density
- ${\lambda}_{i}$
- $${i}_{th}$$ principal stretch
- $J$
- Relative volume defined in Equation 13
- $$\beta =\frac{\nu}{\left(1-2\nu \right)}$$
- ${\alpha}_{i}$ and ${\mu}_{i}$
- Material constants coefficient pairs.

### Viscous (Rate) Effects

## /MAT/LAW69

This law, like /MAT/LAW42 (OGDEN) defines a hyperelastic and incompressible material specified using the Ogden or Mooney-Rivlin material models. Unlike LAW42, where the material parameters are input, this law computes the material parameters using test data from a uniaxial engineering stress-strain curve.

This material can be used with shell and solid elements.

`law_ID`:

`law_ID`= 1 (Ogden law):(22) $$W\left({\lambda}_{1},{\lambda}_{2},{\lambda}_{3}\right)={\displaystyle \sum _{p=1}^{5}\frac{{\mu}_{p}}{{\alpha}_{p}}\left({\overline{\lambda}}_{1}{}^{{\alpha}_{p}}+{\overline{\lambda}}_{2}{}^{{\alpha}_{p}}+{\overline{\lambda}}_{3}{}^{{\alpha}_{p}}-3\right)}+\frac{K}{2}{\left(J-1\right)}^{2}$$ `law_ID`= 2 (Mooney-Rivlin law):(23) $W={C}_{10}\left({I}_{1}-3\right)+{C}_{01}\left({I}_{2}-3\right)$

### Material Parameters

After reading the stress-strain curve (`fct_ID`_{1}),
Radioss calculates the corresponding material
parameter pairs using a nonlinear least-square fitting algorithm. For the classic
Ogden law, (`law_ID`=1), the calculated material parameter pairs
are ${\mu}_{p}$ and ${\alpha}_{p}$ where the value of `p` is defined
via the `N_pair` input. The maximum value is
`N_pair`=5 with a default value of 2. Usually no more than
`N_pair`=3 is needed for a good fit.

For the Mooney-Rivlin law (`law_ID` =2), the material parameter $${C}_{10}$$ and $${C}_{01}$$ are calculated remembering that ${\mu}_{p}$ and ${\alpha}_{p}$ for the LAW42 Ogden law can be calculated using this
conversion:

${\mu}_{1}=2\cdot {C}_{10}$, ${\mu}_{2}=-2\cdot {C}_{01}$, ${\alpha}_{1}=2$, and ${\alpha}_{2}=-2$

The minimum test data input should be a uniaxial tension engineering stress strain curve. If uniaxial compression data is available, the engineering strain should increase monotonically from a negative value in compression to a positive value in tension. In compression, the engineering strain should not be less than -1.0 since -100% strain is physically not possible.

- The experimental data curve represents a smooth monotonically increasing
function with uniform distribution of abscissa points. The number of data
points in the experimental data curve should be greater than the number of
parameter pairs (
`N_pair`). - The engineering strain is negative in compression and positive in tension. For compression test data, the engineering strain should be greater than -1.0 (100% compression maximum) but tension only stress strain data can also be used.
- If
`N_pair`≥ 3, then the test data should cover at least 100% of the tensile strain and/or 50% of the compressive strain. `N_pair`should not be set to a very large value to avoid instabilities in the fitting procedure.

This material law is stable when $${\mu}_{p}{\alpha}_{p}>0$$ (with $$p$$=1,…5) is satisfied for parameter pairs for all
loading conditions. By default, Radioss tries to fit the
curve by accounting for these conditions
(`I`_{check}=2). If a proper fit
cannot be found, then Radioss uses a weaker condition
(`I`_{check}=1: ), which ensures
that the initial shear hyperelastic modulus ($\mu $) is positive.

To determine how well the calculated material parameters represent the input test data, the Radioss Starter outputs an “averaged error of fitting” value which is recommended to not exceed 10%. For visual comparison, the stress-strain curve calculated from the strain energy density and calculated material parameters is also output by the Radioss Starter.

^{3}which are valid for incompressible materials.

- ${\epsilon}_{c}$
- Uniaxial engineering compressive stain
- ${\epsilon}_{b}$
- Equal biaxial engineering tension strain
- ${\sigma}_{c}$
- Uniaxial engineering compressive stress
- ${\sigma}_{b}$
- Equal biaxial engineering tension stress

### Material Incompressibility

Material LAW69 uses the same method to maintain incompressibility as LAW42. For additional information, refer to Poisson's Ratio and Material Incompressibility in LAW42.

### Viscous (Rate) Effects

/VISC/PRONY must be used with LAW69 to include viscous effects. Alternatively, LAW69 could be used to extract the Ogden or Mooney-Rivlin parameters and then those parameters can be used in LAW42 with viscosity added.

## /MAT/LAW82

This material model defines a hyperelastic, and incompressible material specified using the Ogden, Neo-Hookean, or Mooney-Rivlin material models. This law is generally used to model incompressible rubbers, polymers, foams, and elastomers.

- $W$
- Strain energy density
- $N$
- Number of material constants ${\alpha}_{i}$, ${\mu}_{i}$ and ${D}_{i}$
- ${\overline{\lambda}}_{i}={J}^{-\frac{1}{3}}{\lambda}_{i}$
- Deviatoric stretch
- $J$
- Relative volume as defined in Equation 2

- If $\nu =0$, ${D}_{1}$ should be entered.
- If $\nu \ne 0$, ${D}_{1}$ input is ignored and will be recalculated
and output in the Starter output using:
(28) $${D}_{1}=\frac{3(1-2v)}{\mu (1+v)}$$ - If $\nu =0$ and ${D}_{1}$=0, a default value of $\nu =0.475$ is used and ${D}_{1}$ is calculated using Equation 28

### Neo-Hookean Model

Like LAW42, LAW82 can also be simplified to a Neo-Hooken model by using:

${\mu}_{1}=2\cdot {C}_{10}$, ${\alpha}_{1}=2$ and ${\mu}_{2}={\alpha}_{2}=0$

### Mooney-Rivlin Model

Like LAW42, LAW82 can also be simplified to a Mooney-Rivlin model by using:

${\mu}_{1}=2\cdot {C}_{10}$, ${\mu}_{2}=2\cdot {C}_{01}$, ${\alpha}_{1}=2$ and ${\alpha}_{2}=-2$

### Viscous (Rate) Effects

/VISC/PRONY must be used with LAW82 to include viscous effects.

## Drücker Condition Stability Check

In LAW42 and LAW69, the Drücker stability is automatically calculated by the Radioss Starter.

Where, $i$=1,2,3 principal direction.

- $d{\tau}_{i}=J\cdot d{\sigma}_{i}$
- The change of Kirchhoff stress
- $d\tau =D:d\epsilon $
- Relationship between Kirchhoff stress and logarithmic strain

The Drücker stability criterion calculates the strain range where the material model will remain stable given a set of material parameters. This stability check cannot be made for every deformation but instead is commonly used to check material stability under uniaxial, biaxial and planar strain loading.

For example, using the following Ogden parameters:

$\begin{array}{cc}{\mu}_{1}=13.99077258830& {\alpha}_{1}=3.788192935039\\ {\mu}_{2}=-9.13454532223& {\alpha}_{2}=-7.17617341059\\ {\mu}_{3}=8.904655103235& {\alpha}_{3}=-7.27028137148\end{array}$

```
CHECK THE DRUCKER PRAGER STABILITY CONDITIONS
-----------------------------------------------
MATERIAL LAW = OGDEN (LAW42)
MATERIAL NUMBER = 1
TEST TYPE = UNIXIAL
COMPRESSION: UNSTABLE AT A NOMINAL STRAIN LESS THAN -0.3880000000000
TENSION: UNSTABLE AT A NOMINAL STRAIN LARGER THAN 0.9709999999999
TEST TYPE = BIAXIAL
COMPRESSION: UNSTABLE AT A NOMINAL STRAIN LESS THAN -0.2880000000000
TENSION: UNSTABLE AT A NOMINAL STRAIN LARGER THAN 0.2780000000000
TEST TYPE = PLANAR (SHEAR)
COMPRESSION: UNSTABLE AT A NOMINAL STRAIN LESS THAN -0.3680000000000
TENSION: UNSTABLE AT A NOMINAL STRAIN LARGER THAN 0.5829999999999
```

For a Mooney-Rivlin material, the Drücker stability should be checked since ${C}_{01}\text{or}{\mu}_{2}$ could be negative, which leads material instability.

## /MAT/LAW88

This law utilizes a tabulated uniaxial tension and compression engineering stress and strain test data at different strain rates to model incompressible materials. It is only compatible with solid elements.

^{4}

Instead, this law determines the Ogden function directly from the uniaxial engineering stress strain curve tabulated data.

Unlike other Ogden material laws, the Bulk Modulus must be input from either test data or extracted from Starter output of the LAW69 Ogden curve fit. When comparing results between LAW42 or LAW69 to LAW88, the same bulk modulus must be used.

### Unloading Behavior

Unloading can be represented using an unloading function or by providing hysteresis and shape factor inputs to a damage model based on energy.

- $${W}_{cur}$$
- Current energy
- $${W}_{\mathrm{max}}$$
- Maximum energy corresponding to the quasi-static behavior
- $$Hys$$ and $$Shape$$
- Input by user

**Tension****Loading and Unloading**- = 0
- = 1
- = -1

### Viscous (Rate) Effects

Strain rate effects can be modeled by including engineering stress strain test data at different strain rates. This can be easier than calculating viscous parameters for traditional hyperelastic material models.

## Conclusion

Make sure to use the material law that best fits the test data available.

For example, if minimal test data is available and the strains are not too large than the LAW42 Neo-Hookean Model could be used. If the loading state is known, then it is important to have test data that represents that stress state and make sure the material model fits that test data.

### References

^{1}Ogden, R. W., and Non-linear Elastic Deformations. "Ellis Horwood." New York (1984)

^{2}Miller, Kurt. "Testing Elastomers for Hyperelastic Material Models in Finite Element Analysis" Axel Products, Inc., Ann Arbor, MI (2017). Last modified April 5, 2017

http://www.axelproducts.com/downloads/TestingForHyperelastic.pdf

^{3}Axel Products, Inc. "Compression or Biaxial Extension", Ann Arbor, MI (2017). Last modified November 12, 2008

^{4}Kolling, S., P. A. Du Bois, D. J. Benson, and W. W. Feng. "A tabulated formulation of hyperelasticity with rate effects and damage." Computational Mechanics 40, no. 5 (2007): 885-899