The Mullins effect (TB,CDM) is a phenomenon typically observed in compliant filled polymers. It is characterized by a decrease in material stiffness during loading and is readily observed during cyclic loading as the material response along the unloading path differs noticeably from the response that along the loading path. Although the details about the mechanisms responsible for the Mullins effect have not yet been settled, they might include debonding of the polymer from the filler particles, cavitation, separation of particle clusters, and rearrangement of the polymer chains and particles.
In the body of literature that exists concerning this phenomenon, a number of methods have been proposed as constitutive models for the Mullins effect. The model is a maximum load modification to the nearly- and fully-incompressible hyperelastic constitutive models already available. In this model, the virgin material is modeled using one of the available hyperelastic potentials, and the Mullins effect modifications to the constitutive response are proportional to the maximum load in the material history.
The Ogden-Roxburgh [379] pseudo-elastic model (TB,CDM,,,,PSE2) of the Mullins effect is a modification of the standard thermodynamic formulation for hyperelastic materials and is given by:
 | (4–265) |
where
W0(Fij) = virgin material deviatoric strain-energy potential η = evolving scalar damage variable Φ(η) = damage function
The arbitrary limits 0 < η
1 are imposed with η = 1 defined as
the state of the material without any changes due to the Mullins effect.
Then, along with equilibrium, the damage function is defined by:
 | (4–266) |
which implicitly defines the Ogden-Roxburgh parameter η. Using Equation 4–266, deviatoric part of the second Piola-Kirchhoff stress tensor is then:
 | (4–267) |
The modified Ogden-Roxburgh damage function [380] has the following functional form of the damage variable:
 | (4–268) |
where r, m, and β are material parameters and Wm is the maximum virgin potential over the time interval
:
 | (4–269) |
The tangent stiffness tensor
for a constitutive model defined by Equation 4–265 is expressed as follows:
 | (4–270) |
The differential for η in Equation 4–268 is:
 | (4–271) |