Mullins Effect

The Mullins effect material behavior enables you to model the stress softening of filled rubber elastomers under quasi-static cyclic loading. The model requires the specification of three Mullins coefficients.

The Mullins effect model:

  • provides an extension to model damage in hyperelastic material models;
  • is based on the theory of incompressible isotropic elasticity modified by the addition of a single variable, referred to as the damage variable;
  • assumes that only the deviatoric part of the material response is associated with damage;
  • is intended for modeling material response in situations where different parts of the model undergo different levels of damage resulting in a different material response; and
  • is applied to the long-term modulus when combined with viscoelasticity.

Mullins Effect

Hyperelastic materials are described in terms of a “strain energy potential” function U ( F ) that defines the strain energy stored in the material per unit reference volume (volume in the initial configuration). The quantity F is the deformation gradient tensor. To account for the Mullins effect, Ogden and Roxburgh propose a material description that is based on an energy function of the form U ( F , η ) , where the additional scalar variable, η , represents damage in the material. The damage variable controls the material properties in the sense that it enables the material response to be governed by an energy function on unloading and subsequent submaximal reloading different from that on the primary (initial) loading path from a virgin state. Because of the above interpretation of η , it is no longer appropriate to think of U as the stored elastic energy potential. Part of the energy is stored as strain energy, while the rest is dissipated as a result of damage.

In preparation for writing the constitutive equations for Mullins effect, it is useful to separate the deviatoric and the volumetric parts of the total strain energy density as

U = U d e v + U v o l .
U , U d e v , and U v o l are the total, deviatoric, and volumetric parts of the strain energy density, respectively. All the hyperelasticity models described in Isotropic Hyperelasticity have strain energy potential functions that are separated into deviatoric and volumetric parts. For example, the polynomial models use a strain energy potential of the form
U = i + j = 1 N C i j ( I ¯ 1 3 ) i ( I ¯ 2 3 ) j + i = 1 N 1 D i ( J e l 1 ) 2 i ,
where the symbols have the usual interpretations. The first term on the right represents the deviatoric part of the elastic strain energy density function, and the second term represents the volumetric part.

The Mullins effect is accounted for by using an augmented energy function of the form

U ( λ ¯ i , η ) = η U ˜ d e v ( λ ¯ i ) + ϕ ( η ) + U ˜ v o l ( J e l )
where U ˜ d e v ( λ ¯ i ) is the deviatoric part of the strain energy density of the primary hyperelastic behavior, defined, for example, by the first term on the right-hand-side of the polynomial strain energy function given above; U ˜ v o l ( J e l ) is the volumetric part of the strain energy density, defined, for example, by the second term on the right-hand-side of the polynomial strain energy function given above; λ ¯ i ( i = 1 , 2 ) represent the deviatoric principal stretches; and J e l represents the elastic volume ratio. The function ϕ ( η ) is a continuous function of the damage variable η and is referred to as the damage function. When the deformation state of the material is on a point on the curve that represents the primary hyperelastic behavior, η = 1 , ϕ ( η ) = 0 , U ( λ ¯ i , 1 ) = U ˜ d e v ( λ ¯ i ) + U ˜ v o l ( J e l ) and the augmented energy function reduces to the strain energy density function of the primary hyperelastic behavior. The damage variable varies continuously during the deformation and always satisfies 0 < η 1 .

The primary hyperelastic behavior is defined by using the hyperelastic material model (see Isotropic Hyperelasticity). The Mullins effect model is defined by specifying the Mullins effect parameters directly.

With the above modification to the energy function, the stresses are given by

σ ( η , λ ¯ i , J e l ) = η S ˜ ( λ ¯ i ) p ˜ ( J e l ) I ,
where S ˜ is the deviatoric stress corresponding to the primary hyperelastic behavior at the current deviatoric deformation level λ ¯ i and p ˜ is the hydrostatic pressure of the primary hyperelastic behavior at the current volumetric deformation level J e l . Thus, the deviatoric stress as a result of Mullins effect is obtained by simply scaling the deviatoric stress of the primary hyperelastic behavior with the damage variable η . The pressure stress is the same as that of the primary behavior. The model predicts loading/unloading along a single curve (that is different, in general, from the primary hyperelastic behavior) from any given strain level that passes through the origin of the stress-strain plot. It cannot capture permanent strains upon removal of load. The model also predicts that a purely volumetric deformation will not have any damage or Mullins effect associated with it.

The damage variable, η , varies with the deformation according to

η = 1 1 r e r f ( U d e v m U ˜ d e v m + β U d e v m ) ,
where U d e v m is the maximum value of U ˜ d e v at a material point during its deformation history; and r , β , and m are material parameters.

The parameters r , β , and m do not have direct physical interpretations in general. The parameter m controls whether damage occurs at low strain levels. If m = 0 , there is a significant amount of damage at low strain levels. On the other hand, a nonzero m leads to little or no damage at low strain levels.

Parameter Description
r Value of the r coefficient in the Mullins effect model. r must be greater than 1.
m Value of the m coefficient in the Mullins effect model. m must be greater than or equal to zero, and the values of m and β cannot both be zero.
β Value of the β coefficient in the Mullins effect model. β must be greater than or equal to zero, and the values of m and β cannot both be zero.
Use temperature-dependent data Specify material parameters that depend on temperature. A Temperature field appears in the data table.