Quasilinear Lamé Systems

• Quasilinear Lamé systems are nonlinear elasticity models where stress is defined by a spatially invariant elastic tensor that couples displacement and strain.
• The framework employs infinite-order linearization and tensor extraction techniques to ensure unique and stable recovery of material parameters from boundary measurements.
• It distinguishes between isotropic and anisotropic constitutive models, enabling efficient parameter identification even with limited boundary data.

Quasilinear Lamé systems constitute a general class of nonlinear elasticity models in which the stress response is described by a space-independent, nonlinear elastic tensor that depends on both displacement and strain. The governing equations emerge as quasilinear, with the constitutive law $\sigma(u) = C(u, \epsilon(u)) : \epsilon(u)$ for the displacement field $u$ and strain tensor $\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$. Prominent within this theory is the inverse problem: recovering nonlinear material parameters from boundary stress measurements linked to prescribed boundary displacements, an area where recent advances have delivered unique and stable parameter identification for broad classes of isotropic and anisotropic tensors—even with measurements restricted to finite or single boundary points (Johansson et al., 22 Jan 2026).

1. Mathematical Formulation and Constitutive Models

Let $\Omega \subset \mathbb{R}^n$ ($n \geq 2$), a bounded $C^2$ domain with displacement field $u: \Omega \rightarrow \mathbb{R}^n$. The infinitesimal strain is $\epsilon(u) \in S_n(\mathbb{R})$, and the constitutive law prescribes

$$\sigma(u) = C(u, \epsilon(u)) : \epsilon(u) \quad \text{in} \quad \Omega,$$

where $$C: \mathbb{R}^n \times S_n(\mathbb{R}) \rightarrow \mathbb{R}^{n \times n \times n \times n}$$ is the elastic tensor, with Frobenius contraction ":". The force equilibrium (neglecting body forces) requires

$u$0

subject to Dirichlet condition $u$1 on $u$2.

Material Law Classification:

• Isotropic case: The tensor specializes to

$u$3

where $u$4 and $u$5, $u$6 are scalar nonlinear Lamé moduli.

• Anisotropic case: More general forms include terms such as

$u$7

admitting modulation of diagonal components and full fourth-order tensor symmetries.

Under standard symmetry and strong-ellipticity conditions, the quasilinear system admits, for each constant base $u$8 and small boundary perturbation $u$9, a unique solution $\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$0 in $\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$1.

2. Boundary Data and Measurement Protocols

The identification framework utilizes controlled Dirichlet data:

$\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$2

where $\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$3 (base displacement), $\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$4 (small parameter), $\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$5 (finite set of symmetric probing matrices), and $\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$6 for $\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$7.

Measurement Sets:

• Basis and probing sets: Canonical choices are $\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$8 or enriched $\epsilon(u) = \frac{1}{2}(\nabla u + \nabla u^T)$9 for anisotropic cases.
• Boundary stress observation: At each $\Omega \subset \mathbb{R}^n$0, construct the local displacement-to-traction map:

$\Omega \subset \mathbb{R}^n$1

where $\Omega \subset \mathbb{R}^n$2 is the outward unit normal. Data are assembled as pointwise vectors $\Omega \subset \mathbb{R}^n$3 at a finite set of boundary points $\Omega \subset \mathbb{R}^n$4 and for $\Omega \subset \mathbb{R}^n$5.

In the isotropic two-modulus case, measurements at a single boundary point with suitable normal and a single matrix $\Omega \subset \mathbb{R}^n$6 are sufficient.

3. Uniqueness and Stability in Parameter Recovery

3.1 Linearization for Uniqueness

The first-order linearization in $\Omega \subset \mathbb{R}^n$7 yields the Fréchet derivative:

$\Omega \subset \mathbb{R}^n$8

By appropriate choice of $\Omega \subset \mathbb{R}^n$9 (spanning $n \geq 2$0) and evaluating at boundary points with normals $n \geq 2$1 in general position, one establishes:

• Uniqueness (Theorem 3.1): If two tensors agree on pointwise data for all $n \geq 2$2 and $n \geq 2$3, then $n \geq 2$4.
• Isotropic case (Corollary 3.2): Measurements at a single point and matrix suffice to uniquely determine $n \geq 2$5 and $n \geq 2$6.

3.2 Lipschitz Stability

The inverse recovery further benefits from unconditional, quantitative stability:

$n \geq 2$7

with $n \geq 2$8 explicit and dependent only on $n \geq 2$9 and the geometry of $C^2$0.

For isotropic two-modulus recovery:

$C^2$1

4. Nonlinear Expansions and Tensor Extraction

The map $C^2$2 is $C^2$3 in $C^2$4, supporting a Taylor expansion:

$C^2$5

with $C^2$6.

Isolation of material-parameter derivatives is effected by structural calculus:

$C^2$7

enabling explicit multilinear recovery by variation of $C^2$8 across a sufficiently rich set.

5. Reconstruction Algorithms and Explicit Formulas

Infinite-Order Linearization (algorithmic steps):

• Obtain Taylor coefficients $C^2$9 at $u: \Omega \rightarrow \mathbb{R}^n$0 up to desired order $u: \Omega \rightarrow \mathbb{R}^n$1 for each $u: \Omega \rightarrow \mathbb{R}^n$2.
• Recursively resolve parameters: for $u: \Omega \rightarrow \mathbb{R}^n$3, solve $u: \Omega \rightarrow \mathbb{R}^n$4.
• For $u: \Omega \rightarrow \mathbb{R}^n$5, higher-order coefficients isolate $u: \Omega \rightarrow \mathbb{R}^n$6.
• Application of polarization identities yields complete symmetric $u: \Omega \rightarrow \mathbb{R}^n$7-linear forms $u: \Omega \rightarrow \mathbb{R}^n$8.

Explicit formulas rely on linear system inversion for first-order coefficients, and tensorial contraction and orthogonalization for higher orders. Scalar equations are obtained via inner products with orthonormalized test vectors corresponding to the boundary normals.

6. Structural Features and Implications for Inverse Problems

• Displacement-strain coupling ($u: \Omega \rightarrow \mathbb{R}^n$9): Quasilinearity arises from $\epsilon(u) \in S_n(\mathbb{R})$0, but separation by linearization around constant displacements splits $\epsilon(u) \in S_n(\mathbb{R})$1 and $\epsilon(u) \in S_n(\mathbb{R})$2 dependencies.
• Strong ellipticity of $\epsilon(u) \in S_n(\mathbb{R})$3: Guarantees well-posedness, Sobolev regularity, and invertibility for the linearized boundary map.
• Tensor symmetries: Minor and major symmetries of $\epsilon(u) \in S_n(\mathbb{R})$4 reduce effective parameterization and ensure compatibility with probing-set cardinality.
• Isotropy vs. anisotropy: Isotropic (two-modulus) recovery is achievable with minimal data, while anisotropic cases require an increased but finite probe set.
• Infinite-order expansion: Taylor-based order-by-order extraction decomposes the nonlinear inverse problem into sequential linear algebraic inversions.

These combined features render quasilinear Lamé systems tractable for constructive, stable, and finite-dimensional recovery of nonlinear elastic parameters, even from highly restricted boundary stress information (Johansson et al., 22 Jan 2026).