🔢 mathematics

Polyconvexity for Cosserat nonlinear elasticity and nonlinear couple-stress theory

This paper establishes the existence of minimizers for a class of nonlinear Cosserat and couple-stress elastic energies by introducing an independent microrotation and relative stretch to achieve polyconvexity, thereby proving that the constrained model is equivalent to a rotationally regularized deformation problem involving the curvature of the polar factor.

Original authors: Diana Ciotir, Ionel-Dumitrel Ghiba, Patrizio Neff

Published 2026-08-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Diana Ciotir, Ionel-Dumitrel Ghiba, Patrizio Neff

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a piece of rubber being stretched, twisted, or squeezed. In the world of physics, predicting exactly how that material will behave under such forces is a matter of finding the most efficient shape it can take. Scientists use a mathematical tool called energy to describe this efficiency; the material naturally settles into the form that requires the least amount of energy. For decades, researchers have relied on a specific set of rules to prove that such a stable, low-energy shape actually exists and can be found. These rules work beautifully when the material is described by how its internal grid stretches and shears. However, there is a class of materials where the most natural way to describe the stretch does not fit neatly into these standard rules. When scientists try to apply the old methods to these specific materials, the mathematics breaks down, leaving them unable to prove that a stable solution exists at all. This gap is particularly important for materials that have a tiny internal structure, like a sponge or a composite, where the tiny grains inside can rotate independently of the overall stretching.

The researchers Diana Ciotir, Ionel-Dumitrel Ghiba, and Patrizio Neff have developed a new way to bridge this gap. They focused on a type of material model where the internal structure is allowed to rotate on its own, a concept known as a Cosserat continuum. In these models, the material has two distinct parts: the main shape it takes, and the orientation of its tiny internal grains. The difficulty arises because the most natural way to measure the stretch in these materials depends on the relationship between the overall shape and the internal rotation. When scientists tried to prove that a stable solution exists using the standard mathematical tools, they found that the stretch measurement did not behave well enough to guarantee a solution. The internal rotation could wiggle and shift in ways that made the energy calculation unstable, preventing the proof from working.

To solve this, the team introduced a clever mathematical trick. Instead of treating the internal rotation as a fixed part of the stretch, they treated it as a completely independent variable that could move and change freely. They then linked the two together by requiring that the stretch, when viewed from the perspective of the internal rotation, must be a specific type of symmetric, positive shape. This constraint forces the internal rotation to align perfectly with the way the material is actually deforming. By separating the rotation from the stretch in the initial setup, the researchers could use a powerful mathematical property called convexity to prove that a stable solution must exist. The key to making this work was adding a term to the energy calculation that penalizes rapid changes in the direction of the internal rotation. This term acts like a stiffener, ensuring that the rotation field does not wiggle wildly, which in turn allows the mathematics to hold together and prove that a solution exists.

The team proved that this approach works for two different types of energy models. The first type is more complex, requiring the energy to depend on the stretch, its cofactor, and its determinant. The second type is simpler, depending only on the stretch and the determinant, but it requires the stretch to grow at a specific rate to ensure stability. In both cases, the proof shows that a stable configuration exists. Crucially, once the solution is found, the independent internal rotation is no longer just a mathematical helper; it turns out to be exactly the rotation that would naturally arise from the deformation of the material. This means the new method successfully recovers the physical reality of the material while using a more flexible mathematical path to get there.

The researchers also clarified what their method does not do. They showed that while their approach controls the rotation of the material's internal structure, it does not control every possible detail of how the material bends or curves on a microscopic level. It is a selective control, focusing specifically on the rotational aspect rather than the full complexity of the material's second-order shape. Furthermore, the method guarantees that the material does not flip inside out locally, but it does not prove that the material will never pass through itself in a global sense. The solution exists, but it is a solution for a specific, well-behaved version of the problem.

This work is significant because it opens the door to proving the existence of stable solutions for a wide range of materials that were previously out of reach for standard mathematical analysis. It allows scientists to use the most physically natural descriptions of stretch for these complex materials without getting stuck in mathematical dead ends. The method provides a robust framework for understanding how materials with internal structure behave under large deformations. By establishing that a solution exists, the researchers have laid the groundwork for more accurate simulations and a deeper understanding of these complex materials, ensuring that the models used to describe them are grounded in solid mathematical reality. The result is a new, reliable path for analyzing materials where the internal rotation plays a critical role, bridging the gap between physical intuition and mathematical proof.

Technical Summary: Polyconvexity for Cosserat Nonlinear Elasticity and Nonlinear Couple–Stress Theory

Problem Statement
The paper addresses a fundamental tension in nonlinear elasticity: the conflict between mechanically natural strain measures and the structural hypotheses required by the direct method of the calculus of variations. While classical polyconvexity theory (Ball, 1976) is formulated in terms of the deformation gradient $F = D\phi$ and its minors $(F, \text{Cof } F, \det F)$ , many physically transparent energies are naturally expressed in terms of the right stretch tensor $U = \sqrt{F^T F}$ . A primary example is the Biot strain $U - \mathbb{1}$ . However, weak convergence of $F_k$ does not imply useful weak convergence of $U_k$ . Consequently, convexity or polyconvexity formulated directly in $U$ fails to provide the compactness mechanism necessary for existence proofs in first-gradient problems involving only the deformation $\phi$ .

Methodology: The Cosserat Lifting
To resolve this, the authors employ a "Cosserat lifting" strategy. They introduce an independent microrotation field $R: \Omega \to SO(3)$ and define a relative stretch tensor $U = R^T D\phi$ . This variable is not generally symmetric. The energy functional is then formulated in terms of the lifted variables $(U, \text{Cof } U, \det U)$ and a curvature variable $\alpha = R^T \text{Curl } R$ .

The methodology relies on two key analytical mechanisms:

Compactness via Curvature: For rotation fields in $W^{1,q}(\Omega; SO(3))$ with $q \ge 2$ , the row-wise $\text{Curl } R$ controls the full gradient $DR$ (via the Neff–Münch identity). This provides strong compactness for the sequence of rotations $R_k$ , which is essential for passing to the limit in the product $R_k^T D\phi_k$ .
Weak Continuity of Minors: The identities $\text{Cof}(R^T D\phi) = R^T \text{Cof } D\phi$ and $\det(R^T D\phi) = \det D\phi$ allow the minors of the lifted stretch to be identified through the weak continuity of the corresponding minors of the classical deformation gradient.

The paper analyzes two distinct growth regimes:

Cofactor-Dependent Regime: The energy density is convex in $(U, \text{Cof } U, \det U)$ and assumes separate coercivity for the cofactor term.
Reduced Determinant-Dependent Regime: The energy depends only on $(U, \det U)$ without independent cofactor coercivity. Existence here relies on the condition $p > 9/4$ (where $p$ is the growth exponent for $U$ ), which allows the combination of weak $L^{p/2}$ -convergence of cofactors with compact Sobolev convergence of deformations to identify the weak limit of the determinant via the Piola identity.

Key Contributions and Results
The paper establishes the existence of minimizers for a class of nonlinear elastic energies that are not directly covered by standard polyconvex theory in $F$ .

Existence Theorems: The authors prove existence for the auxiliary lifted Cosserat problem (Theorem 4.1) under both the cofactor-dependent and reduced determinant-dependent assumptions.
Polar Recovery and Constraint: A critical result is the proof that for the constrained model, the condition $R^T D\phi \in \text{Sym}^+(3)$ (positive semi-definite) combined with $\det D\phi > 0$ implies that the independent microrotation $R$ must coincide with the polar factor of the deformation gradient, i.e., $R = \text{polar}(D\phi)$ almost everywhere (Proposition 5.1).
Equivalence to Couple-Stress Models: This recovery establishes a bijection between the lifted Cosserat problem and a constrained deformation problem where the energy contains the curvature of the polar factor. The resulting model is identified as a "selectively rotationally regularized" nonlinear elastic model of couple-stress type.
Scope of Applicability: The framework accommodates regularized Biot-type energies, reduced Neo-Hookean models, and cofactor-dependent models (e.g., Mooney-Rivlin, Ciarlet-Geymonat). It explicitly excludes energies that lack global convexity in the lifted variable $U$ (such as the standard quadratic Hencky energy) or those that do not satisfy the necessary coercivity thresholds (e.g., purely quadratic reduced Neo-Hookean without regularization).

Significance and Claims
The paper claims to extract and systematize a variational mechanism previously implicit in the work of Neff and others, providing a complete self-contained argument for existence in this specific setting.

Selective Regularization: The authors emphasize that the result is not an existence theorem for the unregularized first-gradient Biot functional. The curvature term controls the variation of the polar rotation but does not control the variation of the stretch tensor or the full second gradient $D^2\phi$ . Thus, the model is weaker than full second-gradient regularization but stronger than standard first-gradient models regarding rotational stability.
Local vs. Global: The condition $\det D\phi > 0$ ensures local orientation preservation but does not imply global injectivity or exclude interpenetration.
Dimensional Reduction: The choice of the curvature variable $R^T \text{Curl } R$ (rather than the full gradient $DR$ ) is deliberate. It preserves the intrinsic Cosserat structure and provides an objective measure of spatial incompatibility, making the formulation particularly suitable for future dimensional reductions to Cosserat plate, shell, and rod models.
Relation to Existing Theory: The work complements the general micropolar existence theory (where $R$ and $\phi$ remain independent) by imposing a constraint that recovers the polar factor. It also extends the classical polyconvexity theory of Ball to stretch-based energies that were previously difficult to handle without the lifting device.

In summary, the paper provides a robust variational framework for proving existence for a broad class of stretch-based energies by lifting them into a Cosserat setting, utilizing curvature for compactness, and recovering the polar decomposition to link the result back to a constrained nonlinear elasticity problem.

Technical Summary: Polyconvexity for Cosserat Nonlinear Elasticity and Nonlinear Couple–Stress Theory

More like this