Instability of Two-Term Mooney-Rivlin Model

In engineering practice, rubber material models are widely used in large deformation nonlinear response simulations. One such model is the well-known Mooney-Rivlin model, described with other rubber material models in the ADINA Theory and Modeling Guide [1].

In this News we want to focus on the fact that some rubber material models are unstable, in the sense that an unstable material response for certain strain levels may be obtained. If then, in the simulation of a structure, strain levels corresponding to instability are reached, the Newton-Raphson nonlinear iterative solution may have difficulty to converge, and the calculated response may be physically unrealistic, even if no difficulties in the iterative solution have been encountered.

(a) Material 1 described by neo-Hookean model

Figure 1  Uniaxial and biaxial tension plots

For illustration, consider three incompressible materials with Young’s modulus E = 1.0 at small strains. The three materials are represented, respectively, using the neo-Hookean, two-term Mooney-Rivlin, and Ogden models. The large strain response is obtained using the following constants in the material characterizations:

Neo-Hookean material model: C₁ = 1/6

Mooney-Rivlin material model: C₁ = 1/8, C₂ = 1/24

Ogden material model: μ₁ = 0.497, α₁ = 1.3; μ₂ = 0.001, α₂ = 5.0; μ₃ = -0.0079, α₃ = -2.0.
These material constants are taken from equation (7.2.31) of Ref. [2], except that the μ₁, μ₂ and μ₃ are scaled so that the small-strain Young's modulus E = 1.0.

Using these constants, the uniaxial and biaxial tension plots for the three materials described by the models are as shown in Figure 1. Note that these curves have positive slopes over the complete strain ranges, which seems to imply that no instabilities are present and hence will be encountered in a complex nonlinear 2D or 3D simulation with these material constants.

However, when looking deeper and using the ADINA option to plot the stability curves, see Figure 2, we see that the neo-Hookean and Ogden material descriptions are stable, but the Mooney-Rivlin description is unstable (stability indicator < 0.0) in biaxial tension for strains greater than about 0.6.

(c) Material 3 described by Ogden model

Figure 2  Stability plots

Consider the following two analyses in which the instability in the Mooney-Rivlin model is triggered, and in which the subsequent response is physically unrealistic.

Plane stress sheet in biaxial tension

A rubber sheet of material 2 is modeled with one plane stress element and the Mooney-Rivlin model, and two independent forces are applied. First, equal forces are applied to stretch the sheet (point A in the graph). Then the horizontal force is increased, while the vertical force remains constant. As the movie and figure below show, the horizontal displacement decreases as the horizontal force increases. This is surely a physically unrealistic material response.

The next movie shows the same numerical experiment, except that the material 3 described by the Ogden model (with constants given above) is used. The horizontal displacement increases as the horizontal force increases, which is physically realistic.

Cylinder under axial force and internal pressure

A cylinder of material 2, modeled with 27/4 3-D elements and the Mooney-Rivlin model, is stretched by an axial force, then expanded by a deformation-dependent internal pressure. Initially, as the internal pressure is increased, the response is continuous and realistic. But when the internal pressure is increased to above a critical value, the response suddenly jumps in an unrealistic way. A movie is shown at the top of this page.

We can give some general observations:

The neo-Hookean material model (with only one constant) is always stable.

For the Mooney-Rivlin material model, the strain level at which instability occurs depends on the ratio of C₂/ C₁. If the ratio C₂/ C₁ is decreased, the biaxial strain corresponding to instability increases, but the material is always unstable after this strain level is reached.

The Ogden material model is stable for the material constants used above, but for different choices of material constants, the material model may become unstable.

Hence, for any of the rubber material models available, including the material model presented in Reference [3], it is important to look at the stability curves plotted in ADINA prior to performing a general nonlinear simulation using these models.

Keywords:
Rubber material model, Mooney-Rivlin, Ogden, neo-Hookean, hyperelastic, material instability

References

1. ADINA Theory and Modeling Guide, Report ARD 08-7.

2. R. W. Ogden, Non-linear elastic deformations, Ellis Horwood Limited, 1984.

3. T. Sussman and K.J. Bathe, “A model of incompressible hyperelastic material behavior using spline interpolations of tension-compression test data”, Communications in Numerical Methods in Engineering, 25, 53-63, 2009.