Ogden (hyperelastic model)

The Ogden material model is a hyperelastic material model used to describe the non-linear stress-strain behaviour of complex materials such as rubbers, polymers, and biological tissue. The model was developed by Ray W. Ogden in 1972.^[1] The Ogden model, like other hyperelastic material models, assumes that the material behaviour can be described by means of a strain energy density function, from which the stress-strain relationships can be derived. These materials can generally be considered to be isotropic, incompressible and strain rate independent.

Ogden Material Model

In the Ogden material model, the strain energy density is expressed in terms of the principal stretches $\,\!\lambda_j$, $\,\!j=1,2,3$ as:

$W\left( \lambda_1,\lambda_2,\lambda_3 \right) = \sum_{p=1}^N \frac{\mu_p}{\alpha_p}\left( \lambda_1^{\alpha_p} + \lambda_2^{\alpha_p} + \lambda_3^{\alpha_p} -3 \right)$

where N, $\,\!\mu_p$ and $\,\!\alpha_p$ are material constants. Under the assumption of incompressibility one can rewrite as

$W\left( \lambda_1,\lambda_2 \right) = \sum_{p=1}^N \frac{\mu_p}{\alpha_p}\left( \lambda_1^{\alpha_p} + \lambda_2^{\alpha_p} + \lambda_1^{-\alpha_p}\lambda_2^{-\alpha_p} -3 \right)$

In general the shear modulus results from

$2\mu = \sum_{p=1}^{N} \mu_p \alpha_{p}.$

With N = 3 and by fitting the material parameters, the material behaviour of rubbers can be described very accurately. For particular values of material constants, the Ogden model will reduce to either the Neo-Hookean solid (N = 1, α = 2) or the Mooney-Rivlin material (N = 2, α₁ = 2, α₂ = − 2).

Using the Ogden material model, the three principal values of the Cauchy stresses can now be computed as

$\sigma_{\alpha} = p + \lambda_{\alpha}\frac{\partial W}{\partial \lambda_{\alpha}}$

where use is made of $\,\!\sigma_{\alpha} = \lambda_{\alpha}P_{\alpha}$.

Uniaxial tension

We now consider an incompressible material under uniaxial tension, with the stretch ratio given as $\lambda=\frac{l}{l_0}$. The principal stresses are given by

$\sigma_{\alpha} = p + \sum_{p=1}^N \mu_{p} \lambda_{\alpha}^{\alpha_p}$

The pressure p is determined from incompressibility and boundary condition σ₂ = σ₃ = 0, yielding:

$\sigma_{\alpha} = \sum_{p=1}^N\left( \mu_{p} \lambda_{\alpha}^{\alpha_p} - \mu_{p}\lambda_{\alpha}^{-\frac{1}{2}\alpha_p} \right)$

Other Hyperelastic Models

For rubber and biological materials, more sophisticated models are necessary. Such materials may exhibit a non-linear stress-strain behaviour at modest strains, or are elastic up to huge strains. These complex non-linear stress-strain behaviours need to be accommodated by specifically tailored strain-energy density functions.

The simplest of these hyperelastic models, is the Neo-Hookean solid.

$W(\mathbf{C})=\frac{\mu}{2}(I_1^C-3)$

where μ is the shear modulus, which can be determined by experiments. From experiments it is known that for rubbery materials under moderate straining up to 30-70%, the Neo-Hookean model usually fits the material behaviour with sufficient accuracy. To model rubber at high strains, the one-parametric Neo-Hookean model is replaced by more general models, such as the Mooney-Rivlin solid where the strain energy W is a linear combination of two invariants

$W(\mathbf{C})=\frac{\mu_1}{2}\left(I_1^C -3 \right) -\frac{\mu_2}{2}\left(I_2^C - 3\right)$

The Mooney-Rivlin material was originally also developed for rubber, but is today often applied to model (incompressible) biological tissue. For modeling rubbery and biological materials at even higher strains, the more sophisticated Ogden material model has been developed.

References

1. ^ Ogden, R. W., (1972). Large Deformation Isotropic Elasticity - On the Correlation of Theory and Experiment for Incompressible Rubberlike Solids, Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, Vol. 326, No. 1567 (Feb. 1, 1972), pp. 565-584.

• F. Cirak: Lecture Notes for 5R14: Non-linear solid mechanics, University of Cambridge.
• R.W. Ogden: Non-Linear Elastic Deformations, ISBN 0-486-69648-0
• K. Weinberg: Lecture Notes for Zur Methode der finiten Elemente in der Mechanik II: Nichtlineare Probleme, TU Berlin [in English].
  http://mechanik.tu-berlin.de/weinberg/Lehre/fem2/Chapter4.pdf