Nonlinear Isotropic Lamé Moduli

• Nonlinear isotropic Lamé moduli are a generalization of classical Lamé parameters, defined through quadratic expansions of strain-energy density to capture finite strain effects.
• They quantify both volumetric and deviatoric resistance using tangent and secant incremental moduli, ensuring material stability under arbitrary loading conditions.
• Practical identification methods rely on spectral Hessians and Taylor expansions to match small-strain limits while accurately modeling large deformations and pressure effects.

Nonlinear isotropic Lamé moduli generalize the classical Lamé parameters, $\lambda$ and $\mu$, to finite strains and nonlinear constitutive models in hyperelasticity. These moduli, defined either as scalar-valued functions or as parameters obtained from quadratic expansion of the strain-energy density, encapsulate the volumetric and deviatoric resistance of isotropic materials under arbitrary loading, which becomes essential when modeling materials subject to large deformations or strong pressure regimes. The nonlinear context requires reconciling small-strain linearized moduli with strain-dependent "tangent" or incremental moduli, ensuring stability and physical realism over the entire deformation range.

1. Formal Definition and Generalization of Lamé Moduli

For an arbitrary isotropic hyperelastic energy density $\psi(F)$—where $F$ is the deformation gradient—nonlinear Lamé moduli are precisely defined through quadratic expansion about the undeformed state, using spectral decomposition $$F = U\,\mathrm{diag}(\lambda_1,\lambda_2,\lambda_3)V^T$$ and principal stretches $\lambda_i$ (Chen et al., 2024). The first Piola-Kirchhoff stress principal values are $p_i = \frac{\partial\psi}{\partial\lambda_i}$.

PK1-linearization yields: $$\lambda_L = \left.\frac{\partial^2\psi}{\partial\lambda_1 \partial\lambda_2}\right|_{(1,1,1)}, \qquad \mu_L = \frac{1}{2}\left[\left.\frac{\partial^2\psi}{\partial\lambda_1^2}\right|_{(1,1,1)} - \left.\frac{\partial^2\psi}{\partial\lambda_1 \partial\lambda_2}\right|_{(1,1,1)}\right].$$ These generalized moduli characterize the quadratic (second-order) expansion and correspond to the classical $\lambda$, $\mu$ in the small-strain limit.

For moduli associated to the tangent response under finite strain, incremental moduli such as tangent bulk modulus $K_t$ and shear modulus are computed using derivatives of stress with respect to the current state, e.g., for dilatation $F=\lambda I$, $$K_t(\lambda) = 3\,\frac{\partial\sigma_{11}}{\partial\lambda}|_{\lambda_i=\lambda}$$ (Scott, 2020).

2. Small-Strain Limit and Connection to Classical Elasticity

Under small strains, the Green strain $E = \frac{1}{2}(F^TF - I) \approx \varepsilon$ reduces the energy density to

$$\psi(E) \approx \frac{\lambda}{2} (\mathrm{tr}\,\varepsilon)^2 + \mu\,\mathrm{tr}(\varepsilon^2),$$

yielding canonical relations: $$E = \mu_L\,\frac{3\lambda_L + 2\mu_L}{\lambda_L + \mu_L},\qquad \nu = \frac{\lambda_L}{2(\lambda_L + \mu_L)}.$$ The nonlinear definitions exactly recover the classical linear isotropic Lamé parameters in the infinitesimal-strain limit (Chen et al., 2024), providing backward compatibility with linear theory.

3. Strain-Dependent Incremental Moduli: Tangent and Secant Values

Nonlinear elasticity requires explicit consideration of strain-dependent moduli:

• Tangent (incremental) moduli: Quantify the response to infinitesimal perturbations at a finite deformation state. For example:

$$K_t(\lambda) = 3\,\frac{\partial\sigma_{11}}{\partial\lambda}|_{\lambda_i=\lambda}$$

$$E_t(\lambda_1) = \lambda_1\,\frac{\partial\sigma_{11}}{\partial\lambda_1}$$

$$\nu_t(\lambda_1) = -\frac{\lambda_1}{\lambda_2}\,\frac{d\lambda_2}{d\lambda_1}$$

(Scott, 2020)

• Secant (average) moduli: Defined via finite differences between the current stress and the reference configuration. They involve only first derivatives of the energy (Scott, 2020).

A critical observation is that the incremental moduli reflect stability and monotonicity. Their positivity for all strains provides a stringent criterion (pointwise convexity), more restrictive than ground-state conditions.

4. Practical Methodologies for Identification and Tuning

Determination of nonlinear isotropic Lamé moduli leverages Taylor expansions, spectral Hessians, and boundary measurement strategies. For generic isotropic energy forms:

• Compute $\lambda_L$, $\mu_L$ at rest via second derivatives of $\psi(\lambda_1, \lambda_2, \lambda_3)$.
• Prescribe target small-strain properties (Young modulus $E$, Poisson ratio $\nu$), solve for $\lambda_L$, $\mu_L$ using:

$$\lambda_L = \frac{E\nu}{(1+\nu)(1-2\nu)},\qquad \mu_L = \frac{E}{2(1+\nu)}$$

• Material "normalization": Re-parameterize different material models so that their PK1-linearized responses yield identical $\lambda_L$, $\mu_L$; higher-order nonlinearities are then encoded via further parameters or transforms (e.g., the $\alpha$-nonlinearity scaling) (Chen et al., 2024).

Inverse recovery techniques have been studied in quasilinear systems. Given boundary stress measurements under affine displacement, one can uniquely and stably recover the modulus functions $\Lambda(u, \epsilon)$ and $\mu(u, \epsilon)$ (Johansson et al., 22 Jan 2026).

5. Canonical Examples and Strain-Energy Model Comparisons

Various classical and modern hyperelastic families admit closed-form expressions for their generalized Lamé moduli: | Model | Strain Energy Function ($\psi$) | $\lambda_L$ | $\mu_L$ | |----------------------|----------------------------------------------------------------------------------|----------------------------|-----------------------------------------| | Linear Corotational | $$\mu\sum(\lambda_i-1)^2 + \frac{\lambda}{2}(\sum\lambda_i-3)^2$$ | $\lambda$ | $\mu$ | | Neo-Hookean (stable) | $$\frac{\mu}{2}(-3+\sum\lambda_i^2)-\mu(J-1)+\frac{\lambda}{2}(J-1)^2$$ | $\lambda-\mu$ | $\mu$ | | Ogden (N-term) | $$\sum_p\frac{\mu_p}{\alpha_p}(-3+\sum\lambda_i^{\alpha_p})$$ | $0$ | $\frac{1}{2}\sum_p\mu_p(\alpha_p-1)$ | | Classic Neo-Hookean | $$\frac{\mu}{2}(-3+\sum\lambda_i^2)-\mu\log J + \frac{\lambda}{2}(\log J)^2$$ | $\lambda$ | $\mu$ |

Under large strains, the models diverge: LC remains quadratic, SNH softens, and Ogden can be tuned for stiffening or softening. These behaviors can be modulated independently from small-strain properties via scaling transforms (Chen et al., 2024).

6. Stability, Physical Realism, and Positivity Criteria

Incremental moduli serve as sensitive criteria for physically reasonable response. The positivity conditions for $K_t(\lambda)$ and $E_t(\lambda_1)$—which involve second derivatives of the strain energy—are global and often more restrictive than those for ground-state parameters. For the compressible neo-Hookean model:

• $K_t(\lambda) > 0$ for all $\lambda$ if and only if $\lambda_0 > 29\mu_0$.
• $E_t(\lambda_1) > 0$ for all $\lambda_1$ if and only if $\lambda_0 < 12\mu_0$.
• No overlap region exists satisfying both everywhere (Scott, 2020).

A significant phenomenon is that, although the ground-state Poisson ratio $\nu_0$ is positive, the incremental Poisson ratio $\nu_t$ can become negative for sufficiently large axial extensions.

7. Extensions: Pressure Dependence and Extreme Regimes

Under moderate to strong compression, as encountered in geophysical settings, nonlinear Lamé moduli exhibit pressure dependence. The total strain-energy is often decomposed into volumetric and shear-modulated parts: $W(F) = \Phi(J) + \Psi(J)\,W_{\rm shear}(E^*)$ with $\Phi(J)$ capturing bulk behavior (often via Birch–Murnaghan EOS), $\Psi(J)$ modulating the deviatoric term (Kennett, 2021).

Pressure-dependent moduli are then derived: $$\mu(J) = \frac{1}{2J}\Psi(J)\sum_q q a_q, \qquad \lambda(J) = J\Phi''(J) - \frac{1}{3J}\Psi(J)\sum_q q a_q$$ Experimental evidence in earth materials such as MgO indicates $\mu$ rises by $\sim 50\%$ at $100$ GPa, $\lambda$ doubles (Kennett, 2021).

8. Unification and Model Simplicity

All nonlinear isotropic families can be unified by extraction of their generalized Lamé moduli via the spectral Hessian at the rest shape, and augmentation with higher-order terms or nonlinearity scalings to achieve prescribed large-strain behaviors. The linear corotational material is uniquely the simplest nonlinear isotropic model, as it contains only quadratic terms in principal stretches and is the universal PK1-linearization (Chen et al., 2024).

This approach yields a three-parameter family—$(E, \nu, \alpha)$—enabling fully decoupled control over material stiffness, volume preservation, and nonlinearity, facilitating both engineering design and intuitive computational tuning in graphics and simulation contexts.