Cohesive Membranes under determinant constraints

1. Problématique et Contexte

L'article s'inscrit dans le domaine de la mécanique des milieux continus et du calcul des variations, plus précisément dans l'étude de la réduction de dimension pour les membranes élastiques.

• Objectif : Dériver rigoureusement un modèle réduit (2D) à partir d'un modèle tridimensionnel (3D) d'une membrane fine d'épaisseur $\rho$, lorsque $\rho \to 0$.
• Spécificités du modèle :
  • Fracture Cohésive : Le modèle inclut une énergie de surface dépendante de l'amplitude du saut de déplacement $[u]$ (modèle de type Barenblatt généralisé), permettant de capturer des phénomènes de fissuration avec seuil d'activation, contrairement aux modèles de Griffith classiques où l'énergie de surface est constante.
  • Contraintes Déterminantales : Le modèle impose des contraintes physiques fortes sur le déterminant du gradient de déformation $\nabla u$ :
    1. Régime non-pénétration : $\det \nabla u > 0$ (préservation de l'orientation).
    2. Régime incompressible : $\det \nabla u = 1$ (préservation du volume).
• Défi principal : L'absence de croissance uniforme $p$ de l'énergie volumique $W$ sous ces contraintes (notamment $\det F > 0$ ou $=1$) empêche l'application directe des théorèmes de représentation intégrale classiques. De plus, la dépendance de l'énergie de surface par rapport au saut $[u]$ complique la construction des suites de récupération (recovery sequences) tout en respectant les contraintes déterminantales.

2. Méthodologie

Les auteurs utilisent la convergence $\Gamma$ (Gamma-convergence) pour caractériser la limite asymptotique des fonctionnelles d'énergie.

• Fonctionnelle de départ :
  L'énergie totale $F_\rho(v)$ sur le domaine mince $\Sigma_\rho = \Sigma \times (-\rho/2, \rho/2)$ est donnée par :
  $$F_\rho(v) = \int_{\Sigma_\rho} W(\nabla v) dy + \int_{J_v} \psi([v], \nu_v) d\mathcal{H}^2$$
  où $W$ est l'énergie stockée et $\psi$ la densité d'énergie de surface.
• Changement d'échelle : Une transformation de variables ramène le problème sur un domaine fixe $\Sigma_1 = \Sigma \times (-1/2, 1/2)$, conduisant à une énergie redimensionnée $\rho^{-1} G_\rho(u)$.
• Stratégie de preuve :
  1. Inégalité $\Gamma$-liminf : Adaptation des arguments de Braides et Fonseca [15]. La compacité assurée par les hypothèses sur $\psi$ permet de montrer que la limite $u$ est indépendante de la coordonnée transversale $x_3$ et que la composante normale de la discontinuité $(\nu_u)_3$ s'annule presque partout.
  2. Inégalité $\Gamma$-limsup (Construction de suites de récupération) : C'est la contribution technique majeure. Il faut construire une suite $u_{\rho}$ qui converge vers une fonction limite $u$ tout en minimisant l'énergie et en respectant strictement les contraintes $\det \nabla u_\rho > 0$ ou $=1$.

3. Contributions Techniques Clés

Pour surmonter les obstacles liés aux contraintes déterminantales et à la régularité des fonctions, les auteurs développent deux outils mathématiques nouveaux :

A. Approximation des fonctions GSBV avec contraintes

Dans le cas incompressible, les approximations par fonctions affines par morceaux (usuelles en $\Gamma$-convergence) ne sont pas assez régulières pour contrôler la contrainte $\det \nabla u = 1$.

• Résultat (Proposition 3.6) : Construction d'une suite d'approximations $u_j \in W^{k,\infty}$ (lisses) hors de l'ensemble de saut, avec un déterminant contrôlé.
• Méthode : Utilisation de difféomorphismes lisses construits via des convolutions à noyaux variables (Lemma 3.4). Cela permet de mapper le domaine fissuré sur un domaine Lipschitz tout en préservant la régularité et en contrôlant le déterminant du gradient.

B. Difféomorphismes $C^\infty$ optimisant l'énergie de surface

Pour minimiser l'énergie de surface réduite, il faut orienter la normale au saut de manière optimale.

• Résultat (Lemmes 4.8 et Corollaire 5.3) : Construction d'une famille de difféomorphismes $C^\infty$, $f_{j,\rho}$, convergeant vers l'identité dans $W^{k,\infty}$.
  • Ces applications agissent comme des isométries locales près des segments de la fissure.
  • Elles rotent la normale du saut pour minimiser la composante verticale de l'énergie de surface (via l'infimum définissant $\psi_0$).
  • Cas incompressible : Une correction basée sur une équation différentielle ordinaire (ODE) est ajoutée pour garantir que $\det \nabla f_{j,\rho} = 1$ (Lemma 5.2).

4. Résultats Principaux

Sous des hypothèses de croissance et de coercivité raisonnables sur $W$ et $\psi$ (notamment (A1)-(A4) pour le cas non-pénétration et ($\tilde{A}$1)-($\tilde{A}$2) pour l'incompressibilité), les auteurs établissent que la suite de fonctionnelles $\rho^{-1} G_\rho$ $\Gamma$-converge vers une fonctionnelle limite $G_0$ définie sur $\Sigma$ :

$$G_0(u) = \begin{cases} \int_{\Sigma} QW_0(\nabla_\alpha u) dx + \int_{J_u} B\psi_0([u], \nu_u) d\mathcal{H}^1 & \text{si } u \in \mathcal{A} \\ +\infty & \text{sinon} \end{cases}$$

Où :

• $\mathcal{A}$ est l'espace des fonctions indépendantes de $x_3$.
• $W_0$ et $\psi_0$ sont les énergies réduites obtenues par minimisation sur la composante verticale.
• $QW_0$ est l'enveloppe quasi-convexe de $W_0$.
• $B\psi_0$ est l'enveloppe BV-elliptique de $\psi_0$.

Points forts des résultats :

1. Généralité : C'est la première dérivation de dimension réduite pour des membranes avec une énergie de surface dépendante du saut ($[u]$) sous contraintes déterminantales.
2. Unicité du modèle : Le modèle limite est identique pour les régimes non-pénétration et incompressible, bien que les techniques de preuve diffèrent significativement.
3. Rigueur : La preuve gère la non-linéarité forte et la perte de croissance uniforme de l'énergie volumique.

5. Signification et Impact

• Avancée théorique : Ce travail comble un vide important dans la littérature sur la mécanique des fractures et la réduction de dimension. Il résout le problème ouvert de la coexistence entre contraintes de déterminant (physiquement essentielles pour éviter les inter-pénétrations ou garantir l'incompressibilité) et énergies de fracture cohésives dépendantes du saut.
• Outils méthodologiques : Les techniques développées, notamment l'approximation $W^{k,\infty}$ sous contrainte de déterminant via des noyaux variables et la construction de difféomorphismes optimisant la normale, sont susceptibles d'être appliquées à d'autres problèmes de dimension réduite (coques, poutres) ou avec des hypothèses de croissance plus faibles.
• Applications potentielles : Ces modèles réduits sont cruciaux pour la simulation numérique efficace de membranes élastiques (comme les polymères, les tissus biologiques ou les films minces) soumises à de grandes déformations et à la rupture, où les effets d'incompressibilité et de non-pénétration sont dominants.

En résumé, cet article fournit une fondation mathématique rigoureuse pour l'analyse asymptotique des membranes élastiques fracturées, en surmontant les difficultés techniques majeures liées aux contraintes physiques non-linéaires.