Nonlinear Incompressible Shear Wave Models in Hyperelasticity and Viscoelasticity Frameworks, with Applications to Love Waves

This paper presents general nonlinear equations for shear displacements in incompressible hyper-viscoelastic materials, applies them to model nonlinear Love waves at material interfaces, and validates the findings through (2+1)-dimensional numerical simulations that demonstrate how wave speeds evolve over time while generally satisfying linear existence conditions.

The Gist

Imagine the Earth's crust not as a rigid, unyielding rock, but as a giant, complex jelly. Sometimes it's stiff like hard gelatin; other times, it's stretchy like a rubber band or gooey like honey. When an earthquake happens, it sends ripples through this "jelly" just like dropping a stone into a pond sends ripples across the water.

This paper is about understanding how those ripples—specifically a type called Love waves—behave when the "jelly" is really stretchy, really thick, or a mix of both.

Here is the breakdown of the research using simple analogies:

1. The Old Way vs. The New Way

The Old Way (Linear):
For a long time, scientists treated the Earth like a simple spring. If you pull it a little, it snaps back a little. If you pull it twice as hard, it snaps back twice as hard. This is called "Hooke's Law." It works great for tiny, gentle shakes. But if you have a massive earthquake, the ground doesn't just stretch a little; it squishes, twists, and stretches wildly. The old "spring" math breaks down here.

The New Way (Nonlinear Hyperelasticity):
The authors of this paper decided to treat the Earth like a rubber band or playdough. These materials don't just stretch; they get stiffer the more you pull them. They also have memory (viscoelasticity), meaning they don't snap back instantly; they slowly ooze back into shape, like honey dripping off a spoon.

2. What are Love Waves?

Imagine a sandwich. You have a top layer of bread (the Earth's crust) and a bottom layer of filling (the mantle).

• Love waves are like shaking the top slice of bread side-to-side while the bottom stays still.
• The wave gets "trapped" in the top layer, bouncing back and forth between the surface and the boundary where the bread meets the filling.
• In the old math, these waves only exist if the top layer is "slower" (softer) than the bottom layer. If the bottom is softer, the wave escapes and disappears.

3. The Big Discovery: The "Speed Trap"

The researchers built a new, super-complex math model to see what happens when these waves are huge and the ground is acting like a stretchy rubber band.

• The Rule: They found that even in this chaotic, stretchy world, the waves still follow the same basic rule as the simple spring model: The wave gets trapped in the top layer only if the top layer is "slower" than the bottom.
• The Twist: In the real, messy, nonlinear world, the wave speed isn't constant. It's like a car driving on a road that changes texture.
  • At first, the wave might speed up or slow down depending on how hard the ground is squished.
  • But as time goes on and the wave settles down, it eventually picks a "cruising speed" that matches the faster of the two layers. It's like a runner who starts sprinting wildly but eventually settles into a steady jog that matches the speed limit of the track.

4. The "Jelly" vs. The "Honey" (Viscoelasticity)

The team also added viscosity (thickness) to their model.

• Without Viscosity (Pure Jelly): The waves bounce around forever, getting smaller but never truly stopping.
• With Viscosity (Honey): The waves lose energy quickly. Imagine running through waist-deep water; you slow down fast. The math showed that this "thickness" prevents the waves from getting too crazy and keeps the model realistic. It acts like a shock absorber, smoothing out the jagged edges of the earthquake's energy.

5. How They Tested It (The Computer Simulation)

Since you can't easily shake a giant block of the Earth in a lab, they used a computer to simulate it.

• They created a digital "sandwich" (a top layer and a bottom layer).
• They dropped a virtual "bomb" (a Gaussian explosion) in the middle to simulate an earthquake.
• They watched the ripples spread out.
• The Result: The computer showed that even with the complex, stretchy math, the waves behaved surprisingly predictably. They stayed trapped in the top layer, and their speed eventually settled down to match the faster material, just like the simple theory predicted.

6. Why Does This Matter?

• Better Earthquake Models: By understanding how the ground acts like a stretchy rubber band rather than a stiff spring, we can predict how much damage an earthquake might do to buildings.
• Medical Applications: This same math applies to human tissues (like muscles and tendons), which are also stretchy and gooey. Understanding how waves move through "rubber" helps doctors use ultrasound to see inside the body better.
• Oil and Gas Exploration: Geologists use these waves to find oil. If they understand the nonlinear math, they can get a clearer picture of what's underground.

The Bottom Line

This paper is like upgrading the map of a city. The old map said, "The road goes straight." The new map says, "The road curves, gets bumpy, and changes speed depending on how heavy your car is."

The researchers proved that even with all these new, complicated curves and bumps, the traffic (the earthquake waves) still follows the main rules of the road: it stays in the right lane and eventually finds a steady speed. This gives scientists a much more accurate tool to understand the Earth's most violent movements.

Technical Summary

1. Problem Statement

The paper addresses the modeling of nonlinear shear waves (specifically Love-type waves) propagating through incompressible hyperelastic and hyper-viscoelastic materials. While linear elasticity theory (Hooke's Law) is standard for small deformations in seismology and material science, it fails to capture the behavior of materials undergoing large finite deformations (e.g., polymers, rubbers, biological tissues, and geological formations near earthquake foci).

The specific challenges addressed include:

• Deriving governing equations for shear displacements in incompressible media where the strain energy density function ($W^H$) is arbitrary.
• Incorporating viscoelastic effects (energy dissipation) into the nonlinear framework.
• Analyzing the propagation of Love waves (horizontally polarized transverse waves) at the interface between two layers with different mechanical properties, specifically investigating how nonlinearities affect wave speed, existence conditions, and stability compared to the linear case.

2. Methodology

A. Theoretical Framework

The authors utilize a Lagrangian (material) framework for continuum mechanics.

• Constitutive Models: The study focuses on generalized neo-Hookean materials, where the strain energy density $W^H$ depends only on the first invariant ($I_1$) of the Cauchy-Green deformation tensor, avoiding the complexities of the second invariant ($I_2$).
• Specific Model: The Cubic Yeoh model is selected as a representative constitutive law:
  $$W^H = \mu(I_1 - 3) + \frac{1}{2}M(I_1 - 3)^2 + \frac{1}{3}A(I_1 - 3)^3$$
  This introduces cubic and quintic nonlinear terms.
• Viscoelasticity: A hyper-viscoelastic approach is adopted, adding a dissipative potential $W^V$ dependent on the rate of deformation ($\dot{C}$). The specific potential used is linear in the pseudo-invariant $J_2 = \text{tr}(\dot{C}\dot{C})$.
• Derivation: Starting from the conservation of mass, linear momentum, and angular momentum, the authors derive the equations of motion. They explicitly handle the hydrostatic pressure ($p$) required to enforce incompressibility ($J=1$).

B. Governing Equations

• General Nonlinear Shear Wave Equation: For an arbitrary $W^H(I_1)$, the displacement $v$ in the $Y$-direction satisfies:
  $$\rho_0 v_{tt} = 2W'(v_X^2 + v_Z^2)(v_{XX} + v_{ZZ}) + 4W''(v_X^2 + v_Z^2)(v_X^2 v_{XX} + 2v_X v_Z v_{XZ} + v_Z^2 v_{ZZ})$$
• Viscoelastic Extension: The equation is extended to include a viscosity term $Q[v]$ involving mixed derivatives (e.g., $v_{xxt}$), representing dispersion and dissipation.
• 1D Reduction: For symmetry analysis, the model is reduced to one spatial dimension.

C. Numerical and Analytical Techniques

• Method of Lines (MOL): The 2D partial differential equations (PDEs) are discretized in space using finite differences on a non-uniform mesh (clustered near the source and interface) and solved as a system of coupled ODEs using MATLAB's ode23.
• Symmetry Analysis: Lie symmetry methods are applied to the 1D reduced equation to find exact, scaling-invariant solutions.
• Simulation Setup: Simulations involve a "Gaussian explosion" initial condition in a two-layered domain (upper layer of thickness $L$, lower half-space) to track wave propagation, reflection, and refraction at the interface and free surface.

3. Key Contributions

1. Derivation of General Nonlinear Equations: The paper presents a rigorous derivation of the nonlinear shear wave equation for incompressible hyperelastic materials with arbitrary strain energy functions, specifically highlighting the role of the cubic Yeoh model.
2. Hyper-Viscoelastic Formulation: It extends the model to include viscoelasticity via a pseudo-strain energy density, resulting in a complex PDE with mixed derivative terms that account for energy dissipation and dispersion.
3. Correction of Previous Literature: The authors identify and correct a subtle error in previous derivations regarding the substitution of invariants before expanding the Piola-Kirchhoff stress tensor, clarifying the necessity of restricting to generalized neo-Hookean materials ($W^H_2 = 0$) to ensure well-defined pressure.
4. Exact Solutions via Symmetry: The authors derive exact traveling wave solutions for the 1D hyper-viscoelastic model using Lie symmetry analysis, providing a benchmark for numerical behavior.

4. Results

A. Linear vs. Nonlinear Behavior

• Existence Condition: In the linear theory, Love waves exist only if the wave speed in the upper layer ($c_1$) is less than that in the lower half-space ($c_2$), i.e., $c_1 < |v| < c_2$.
• Nonlinear Findings: Numerical simulations of the fully nonlinear model demonstrate that variable wave speeds of interface and surface waves generally satisfy the linear existence condition $c_1 < |v| < c_2$ initially. However, as time progresses, the wave speed tends toward the larger of the two material wave speeds ($\max(c_1, c_2)$).
• Wave Trapping: When $c_1 < c_2$, waves are trapped in the upper layer (waveguide effect). When $c_1 > c_2$, energy decays into the half-space, and no stable Love wave forms.

B. Viscoelastic Effects

• Dispersion and Damping: The inclusion of the viscosity potential leads to rapid dispersion of wave energy. The amplitude of the wave decays significantly faster compared to the purely hyperelastic case.
• Stability: Viscoelasticity stabilizes the wave profile, preventing the formation of unbounded gradients that might occur in purely nonlinear elastic models. The wave envelopes become more predictable as the viscosity coefficient ($\eta$) increases.

C. Numerical Observations

• Wave Speed Evolution: Initially, surface and interface velocities differ. Over time, the outermost wave reflected from the surface dominates, causing the observed phase velocity to converge to the maximum material speed.
• Wave Breaking (Appendix): The authors report the discovery of numerically breaking wave solutions in radially symmetric configurations. While the displacement $u$ remains continuous, its spatial derivative ($u_r$) develops sharp gradients and spurious oscillations, indicating a numerical singularity or "breaking" similar to shock formation.

5. Significance and Applications

• Seismology: The model provides a more accurate framework for understanding surface waves (Love waves) near earthquake foci where finite displacements and nonlinear material responses are significant. It suggests that nonlinear effects may balance dispersion in linear models.
• Material Science: The derived equations are applicable to the characterization of soft, complex materials like biological tissues, rubbers, and foams under large deformations.
• Theoretical Advancement: By combining hyperelasticity with viscoelasticity in a rigorous continuum mechanics framework, the paper bridges the gap between linear wave theory and fully nonlinear finite deformation theory, offering new tools for analyzing wave propagation in heterogeneous media.

In conclusion, the paper successfully generalizes Love wave theory to the nonlinear, incompressible, and viscoelastic regimes, demonstrating that while nonlinearities alter wave speeds and profiles, the fundamental existence conditions derived from linear theory remain robust, with the system eventually tending toward the fastest available wave speed in the medium.