Achieving Material Robustness via Symmetric Stress Finite Element Discretizations

Résumé technique : Atteinte de la robustesse matérielle via des discrétisations par éléments finis de contraintes symétriques

Énoncé du problème
En mécanique des milieux continus, les formulations variationnelles telles que l'élasticité de Hellinger–Reissner et les formulations vitesse–contrainte pour les écoulements incompressibles impliquent des tenseurs de contraintes symétriques. La symétrie du tenseur de contrainte de Cauchy est physiquement liée à la conservation du moment cinétique. Lors de la discrétisation de ces problèmes à l'aide d'éléments finis conformes en $H(\text{div})$, les praticiens doivent choisir entre imposer la symétrie fortement (exigeant que les tenseurs discrets soient symétriques ponctuellement) ou faiblement (imposant la symétrie via un multiplicateur de Lagrange). Bien que la symétrie faible soit souvent préférée dans la littérature en raison de la complexité perçue de la construction d'éléments fortement symétriques, cet article examine l'impact de ce choix sur la précision numérique, en particulier pour les matériaux présentant des lois constitutives anisotropes.

Méthodologie
Les auteurs développent un cadre théorique et une série de benchmarks numériques pour comparer les schémas à symétrie forte et à symétrie faible.

1. Formulation mathématique : L'étude considère un problème statique linéarisé général pour les fluides et les solides, où la contrainte $\sigma$ est liée à la déformation $\varepsilon(u)$ et à un terme de champ externe $\tilde{F}$ (représentant l'anisotropie). Le problème est formulé comme un système de point selle.
2. Schémas de discrétisation :
   • Symétrie forte : Les schémas de Johnson-Mercier-Křížek (JMK) et de Hu-Zhang (HZk) sont utilisés. Ces schémas imposent $\sigma_h \in \Sigma^{\text{sym}}_h$ ponctuellement.
   • Symétrie faible : Les schémas PEERS et Arnold-Falk-Winther (AFWk) sont utilisés. Ces schémas introduisent un multiplicateur de Lagrange $\omega_h$ (représentant la partie antisymétrique du gradient) pour imposer la symétrie au sens faible.
3. Construction des benchmarks : Les auteurs construisent des « solutions manufacturées » où la contrainte exacte est identiquement nulle ($\sigma \equiv 0$) mais où les champs de déplacement et de rotation sont non triviaux. Ces solutions sont dérivées de lois constitutives inspirées par :
   • Les solides isotropes linéaires.
   • Les solides transversement isotropes (matériaux renforcés par des fibres).
   • Les fluides polaires (réseaux de polymères de cristaux liquides).
     L'amplitude du terme de force anisotrope est contrôlée par un paramètre d'échelle $\delta$.
4. Analyse théorique : Une théorie abstraite des espaces de Hilbert est développée pour les problèmes de point selle génériques. Les auteurs définissent des discrétisations « préservant la structure » et dérivent des estimations d'erreur a priori pour expliquer le comportement numérique observé.

Contributions clés

• Définition de la robustesse matérielle : L'article introduit le concept de « robustesse matérielle », défini comme la propriété selon laquelle un schéma numérique produit une approximation de contrainte nulle pour toute configuration de contrainte exacte nulle, indépendamment de la loi constitutive ou de l'amplitude de la force anisotrope.
• Identification des modes de défaillance : Les auteurs démontrent que les schémas à symétrie faible peuvent produire des approximations de contrainte arbitrairement médiocres (des erreurs échelonnant avec le paramètre de force $\delta$) même lorsque la contrainte exacte est nulle. Cela se produit spécifiquement lorsque le tenseur de rotation exact $\omega = \text{skw}(\nabla u)$ n'appartient pas à l'espace discret des multiplicateurs de Lagrange $\Xi_h$.
• Théorie unificatrice : Une théorie rigoureuse est présentée montrant que les schémas à symétrie forte satisfaisant des propriétés d'inclusion de noyau spécifiques ($\ker B_h \subseteq \ker B$) sont intrinsèquement robustes sur le plan matériel. En revanche, les schémas à symétrie faible manquent de cette invariance à moins que l'espace de rotation discret ne contienne le tenseur de rotation exact.
• Estimations d'erreur : De nouvelles estimations d'erreur a priori sont dérivées pour les discrétisations partiellement préservant la structure, montrant explicitement comment l'erreur de contrainte dépend de la qualité de l'approximation du tenseur de rotation dans le cas à symétrie faible.

Résultats
Les expériences numériques sur des maillages non structurés confirment les prédictions théoriques :

• Cas isotrope : Pour les solides isotropes linéaires où la rotation exacte est constante (et donc facilement représentable), les schémas forts et faibles se comportent de manière similaire, les erreurs étant dominées par les erreurs d'arrondi.
• Cas anisotropes : Pour les solides transversement isotropes et les fluides polaires, où le tenseur de rotation est complexe :
  • Schémas à symétrie forte (JMK, HZ3) : Produisent systématiquement des erreurs de contrainte au niveau de la précision machine (état de contrainte nulle), indépendamment du paramètre d'échelle $\delta$.
  • Schémas à symétrie faible (PEERS, AFWk) : Produisent des erreurs de contrainte qui échelonnent linéairement avec $\delta$. Même pour des schémas d'ordre très élevé (par exemple AFW3), si le tenseur de rotation exact n'est pas contenu dans l'espace discret $\Xi_h$, l'erreur de contrainte reste importante.
• Comportement transitoire : Un exemple de fluide polaire transitoire montre que, bien que les deux schémas s'accordent initialement, le schéma à symétrie faible non robuste échoue à se relâcher vers le bon état stationnaire de contrainte nulle, tandis que le schéma robuste à symétrie forte y parvient.

Signification et affirmations
L'article affirme que le choix entre l'imposition forte et faible de la symétrie n'est pas simplement une question de complexité de mise en œuvre, mais affecte fondamentalement la fidélité physique de la solution numérique. Les auteurs soutiennent que l'imposition forte de la symétrie est cruciale pour la robustesse matérielle, garantissant que les schémas numériques ne génèrent pas de contraintes spuriaires dans des configurations sans contrainte, en particulier dans les applications impliquant des invariants de rotation ou des matériaux anisotropes.

Le travail plaide pour une adoption plus large des éléments finis à symétrie forte (tels que JMK et Hu-Zhang), notant que ces éléments sont de plus en plus disponibles dans des logiciels open source comme Firedrake et NGSolve. Les auteurs concluent que pour les problèmes sensibles à la conservation du moment cinétique, les discrétisations à symétrie forte offrent une garantie nécessaire de précision que les schémas à symétrie faible ne peuvent fournir sans conditions restrictives sur les espaces discrets.