Compare Neo-Hookean, Mooney-Rivlin, and Ogden constitutive models side by side. Compute stress-stretch curves and strain energy density in real time for three deformation modes.
What exactly is a "hyperelastic" material? Is it just a fancy word for rubber?
🎓
Basically, yes! Hyperelasticity is the theory we use to model materials like rubber, silicone, and some biological tissues. The key idea is that the material is perfectly elastic—it returns to its original shape after you remove the force—and its stress comes from a "strain energy density" function. In practice, this means we can model huge, nonlinear deformations that metals can't handle. Try moving the "Stretch (λ)" slider above from 1 to 3; you'll see the stress rise nonlinearly for all three models.
🙋
Wait, really? So why are there three different models (Neo-Hookean, Mooney-Rivlin, Ogden) in the simulator? Which one is "correct"?
🎓
Great question! None is universally "correct"; they offer a trade-off between simplicity and accuracy. The Neo-Hookean model is the simplest, using just one material constant, $C_{10}$. It works well for moderate stretches. Mooney-Rivlin adds a second constant, $C_{01}$, which makes it more accurate for a wider range of deformations, like biaxial stretching. The Ogden model is the most flexible, using multiple terms (with parameters $\mu_i$ and $\alpha_i$) to fit complex experimental data perfectly. In the simulator, switch between them using the dropdown and see how their stress-stretch curves differ!
🙋
So if I'm simulating a car tire, which model should I pick? And how do I know what values to put in for those C10 or μ parameters?
🎓
For a tire, which undergoes complex multiaxial loading, you'd likely start with a Mooney-Rivlin or Ogden model. The parameter values come from fitting the model to experimental data from a uniaxial tension test—exactly the kind of curve you're generating in this simulator! A common case is setting $C_{10}$ and $C_{01}$ so their sum relates to the initial shear modulus. Play with the $C_{10}$ and $C_{01}$ sliders: notice how increasing $C_{10}$ raises the entire curve, while adjusting $C_{01}$ changes the curvature, especially at higher stretches.
Physical Model & Key Equations
The core of hyperelasticity is the strain energy density function, $W$, which measures the energy stored per unit volume due to deformation. For an incompressible material (like rubber), the Cauchy stress for uniaxial tension is derived from the derivative of $W$ with respect to the stretch, $\lambda$.
$$ \sigma = \lambda \frac{\partial W}{\partial \lambda}- p $$
Here, $\sigma$ is the true (Cauchy) stress, $\lambda$ is the stretch ratio (deformed length / original length), and $p$ is a hydrostatic pressure eliminated using the incompressibility condition ($\lambda_1 \lambda_2 \lambda_3 = 1$). The different models define $W$ differently.
The simulator implements the following stress-stretch relations for a uniaxial test, which are the final working equations after applying the incompressibility condition.
The physical meaning is clear when you change them: $C_{10}$ and $\mu_i$ scale the stress, while $\alpha_i$ and the $C_{01}/C_{10}$ ratio define the curve's shape.
Real-World Applications
Automotive Tires: Hyperelastic models are critical for simulating tire deformation, contact with the road, and rolling resistance. Engineers use Mooney-Rivlin or Ogden models to predict stress distributions, heat generation, and wear, ensuring safety and longevity under complex multiaxial loads.
Medical Devices (Silicone Implants & Heart Valves): Silicone breast implants and polymeric heart valves undergo large, repeated deformations. Accurate hyperelastic modeling is essential for predicting fatigue life, ensuring the material won't fail, and optimizing the design for biocompatibility and performance.
Consumer Electronics (Seals & Gaskets): The rubber seals on your smartphone or watch need to withstand compression and environmental exposure. CAE simulations using Neo-Hookean or Mooney-Rivlin models help design seals that maintain water resistance without excessive stress that could damage the housing.
Sporting Goods: From shoe soles to bicycle grips and football bladders, hyperelastic models help optimize material formulations for energy return, grip, and durability. For instance, simulating a running shoe's sole compression helps balance cushioning and responsiveness.
Common Misconceptions and Points to Note
First, understand that "material constants are not catalog values for the material." For example, when you hear that a certain rubber has a C10=0.5 MPa, you might tend to think it's a direct indicator of tensile strength or hardness. However, these constants are merely the result of fitting a specific mathematical model called a "strain energy density function." Even for the same material, the C10 value fitted for Neo-Hookean and the C10 fitted for Mooney-Rivlin are physically different values. Therefore, it's dangerous to use constants from literature or other companies' data as-is; you must always re-fit them based on your own test data.
Next, be aware of the pitfall that "not everything is determined by uniaxial tensile data alone." While uniaxial testing is simple, relying solely on it often leads to underestimating behavior in biaxial or shear modes. For instance, if you determine the C01 for Mooney-Rivlin from only uniaxial data, there are cases where you could underestimate the stress by more than 50% when an actual component undergoes biaxial deformation. In practice, ideally, you should acquire test data from compression or planar tension (close to biaxial) if possible, and adjust the parameters so the curve fits across multiple deformation modes. With this simulator, you can switch deformation modes and see at a glance how a curve fitted for uniaxial data deviates in other modes.
Finally, "a higher order N for the Ogden model is not always better." While the Ogden model can fit complex curves by increasing terms (N=1,2,3...), there's a risk of overfitting. With N=3 or higher, you might get a perfect fit within the error range of the test data, but the model can exhibit unrealistic behavior between data points (interpolation) or beyond them (extrapolation). Start with N=1 or 2, and judge the necessity based on the quality of fit. Computational cost is also a non-negligible factor in CAE.
Set stretch ratio λ (Lambda) using slLambda slider, typically ranging 1.0–3.0 for rubber materials
Adjust material coefficients: C10 (Neo-Hookean shear modulus term, ~0.1–1.0 MPa) and C01 (Mooney-Rivlin second invariant coupling, 0–0.5 MPa)
Observe real-time stress-stretch curves and strain energy density (W) outputs across Neo-Hookean, Mooney-Rivlin, and Ogden models for comparative analysis
Worked Example
Natural rubber under uniaxial tension: λ=2.5, C10=0.4 MPa, C01=0.1 MPa. Neo-Hookean model yields W=0.9 MJ/m³ and Cauchy stress σ≈1.2 MPa. Mooney-Rivlin introduces C01 coupling, reducing predicted stress to σ≈1.05 MPa by accounting for second invariant effects. Ogden model with α₁=1.3, μ₁=0.6 MPa provides nonlinear hardening at high stretches (λ>2.0), better matching experimental tear-strength data for industrial elastomer seals and gaskets.
Practical Notes
Neo-Hookean underestimates stress at high stretches (λ>2.5); use Mooney-Rivlin or Ogden for automotive seals and vibration dampers requiring accurate energy dissipation prediction
C01/C10 ratio quantifies material nonlinearity: C01/C10 > 0.2 indicates pronounced finite-strain stiffening in silicone rubbers vs. natural rubber (typically <0.1)
Validate model choice against DMA (dynamic mechanical analysis) or tensile test data before deployment in FEA for dynamic sealing applications