/MAT/LAW42 (OGDEN)

Block Format Keyword This keyword defines a hyperelastic, viscous, and incompressible material specified using the Ogden, 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.

Format


(1)	(3)	(5)	(6)	(7)	(9)	(10)
/MAT/LAW42/`mat_ID`/`unit_ID` or /MAT/OGDEN/`mat_ID`/`unit_ID`
`mat_title`
$ρ_{i}$
$ν$	$σ_{cut}$		`fct_ID`_blk	`Fscale`_blk	`M`	`I`_form
$μ_{1}$	$μ_{2}$	$μ_{3}$		$μ_{4}$	$μ_{5}$
$μ_{6}$	$μ_{7}$	$μ_{8}$		$μ_{9}$	$μ_{10}$
$α_{1}$	$α_{2}$	$α_{3}$		$α_{4}$	$α_{5}$
$α_{6}$	$α_{7}$	$α_{8}$		$α_{9}$	$α_{10}$

If M > 0


(1)	(3)	(5)	(7)	(9)
`G`₁	`G`₂	`G`₃	`G`₄	`G`₅
. . . `M` values of G (five per line)
$τ_{1}$	$τ_{2}$	$τ_{3}$	$τ_{4}$	$τ_{5}$
. . . `M` values of $τ$ (five per line)

Definition


Field	Contents	SI Unit Example
`mat_ID`	Material identifier (Integer, maximum 10 digits)
`unit_ID`	Unit identifier (Integer, maximum 10 digits)
`mat_title`	Material title (Character, maximum 100 characters)
$ρ_{i}$	Initial density (Real)	$[\frac{kg}{m^{3}}]$
$ν$	Poisson's ratio. 3 Default = 0.495 (Real)
$σ_{cut}$	Cut-off stress in tension. Default = 10³⁰ (Real)	$[Pa]$
`fct_ID`_blk	Function identifier which scales the bulk coefficient as a function of the relative volume. 6 (Integer)
`Fscale`_blk	Bulk function scale factor. Default = 1.0 (Real)
`M`	Number of viscous terms in the Prony series (order of the Maxwell model). (Integer)
`I`_form	Strain energy density formulation flag. 1 = 1 (Default) Standard strain energy density. = 2 Modified strain energy density to better model incompressibility.
$μ_{p}$	`p`^th parameter of the shear hyperelastic modulus (up to 10) (Real)	$[Pa]$
$α_{p}$	`p`^th material exponent (up to 10). (Real)
`G`_i	`i`_th multiplier of the Prony viscous term. 7 (Real)	$[Pa]$
$τ_{i}$	`i`_th time relaxation of the Prony viscous term. (Real)	$[s]$

Example (Rubber)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                  kg                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/OGDEN/1/1
rubber
#              RHO_I
                1E-6
#                 Nu           sigma_cut           funIDbulk         Fscale_bulk         M     Iform
                .495                   0                   0                   0         0         0
#               Mu_1                Mu_2                Mu_3                Mu_4                Mu_5
                2e-3               -1e-3                   0                   0                   0
#               Mu_6                Mu_7                Mu_8                Mu_9               Mu_10
                   0                   0                   0                   0                   0
#            alpha_1             alpha_2             alpha_3             alpha_4             alpha_5
                   2                  -2                   0                   0                   0
#            alpha_6             alpha_7             alpha_8             alpha_9            alpha_10
                   0                   0                   0                   0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Comments

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.LAW42 uses the following strain energy density representation of the Ogden material model.
- I_form=1 (default):
  $W (λ_{1}, λ_{2}, λ_{3}) = \sum_{p = 1}^{10} \frac{μ_{p}}{α_{p}} ({\bar{λ}}_{1}^{α_{p}} + {\bar{λ}}_{2}^{α_{p}} + {\bar{λ}}_{3}^{α_{p}} - 3) + \frac{K}{2} {(J - 1)}^{2}$
- I_form=2:
  $W (λ_{1}, λ_{2}, λ_{3}) = \sum_{p = 1}^{10} \frac{μ_{p}}{α_{p}} ({\bar{λ}}_{1}^{α_{p}} + {\bar{λ}}_{2}^{α_{p}} + {\bar{λ}}_{3}^{α_{p}} - 3) + K {(J - 1 - \ln J)}^{}$
Where,

$W$

Strain energy density

$λ_{i}$

i^th principal engineering stretch

$J$

Relative volume defined as: $J = λ_{1} \cdot λ_{2} \cdot λ_{3} = \frac{ρ_{0}}{ρ}$

${\bar{λ}}_{i} = J^{- \frac{1}{3}} λ_{i}$

Deviatoric stretch

$α_{p}$ and $μ_{p}$

Material constants coefficient pairs.

Up to 10 material constant pairs can be defined.

The initial shear modulus $μ$ and bulk modulus ( $K$ ) are given by:
$μ = \frac{\sum_{p = 1}^{10} μ_{p} \cdot α_{p}}{2}$
and
$K = μ \cdot \frac{2 (1 + ν)}{3 (1 - 2 ν)}$
Where, $ν$ is the Poisson's ratio and is only used for computing the bulk modulus.
Parameters $α_{p}$ and $μ_{p}$ must be chosen so that initial shear modulus is $μ > 0$ .
For material stability, it is required that each material constant pair $μ_{p} \cdot α_{p} > 0$ .
Poisson's ratio $ν$ is used only for computing the bulk modulus, Equation 4.
For pure incompressible materials, $ν = 0.5$ . This value leads to a finite bulk modulus ( $K$ ). Therefore, the recommended maximum Poisson's ratio for incompressible materials is $ν = 0.495$ .
Higher values of the Poisson's ratio may lead to a small time step value or divergence in case of implicit and explicit simulations.
A particular case of the Ogden material model is the incompressible Mooney-Rivlin model, which can be represented using the following equation for the strain energy density function:
$W = C_{10} (I_{1} - 3) + C_{01} (I_{2} - 3)$
Where, $I_{1}$ and $I_{2}$ are the first and second invariants of the right Cauchy-Green Tensor.
This representation can be derived from the Ogden strain energy density function when:
$μ_{1} = 2 \cdot C_{10}$
$μ_{2} = - 2 \cdot C_{01}$
$α_{1} = 2$
$α_{2} = - 2$
A simple case of the Ogden material model is the Neo-Hookean model represented using the following equation for the strain energy density function:
$W = C_{10} (I_{1} - 3)$
Where,

$I_{1}$

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:
$μ_{1} = 2 \cdot C_{10}$
$α_{1} = 2$
and
$μ_{2} = α_{2} = 0$
In cases when the bulk modulus of a material is not large enough to prevent compression, LAW42 allows the input of a function (fct_ID_blk) which scale the bulk modulus as a function of the relative volume so that the bulk modulus can be increased to maintain incompressibility.
Viscous (rate) effects are modeled in LAW42 using a Maxwell model, which can be described in a simplified manner as a system of n springs with stiffness' $G_{i}$ and dampers $η_{i}$ :
Figure 1. Maxwell Model

The Maxwell model is represented using Prony series inputs ( $G_{i}, τ_{i}$ ). The hyperelastic initial shear modulus $μ$ is the same as the long-term shear modulus $G_{\infty}$ in the Maxwell model, and $τ_{i}$ is the relaxation time:
$τ_{i} = \frac{η_{i}}{G_{i}}$
The $G_{i}$ and $τ_{i}$ values must be positive.
/VISC/PRONY can be used with this material law to include viscous effects.