In search of constitutive conditions in isotropic hyperelasticity: polyconvexity versus true-stress-true-strain monotonicity

Cet article démontre que ni la polyconvexité ni la monotonie contrainte-vraie/déformation-vraie ne garantissent à elles seules un comportement physiquement raisonnable en hyperélasticité isotrope, suggérant que bien que leur combinaison constitue une solution prometteuse au Hauptproblem de Truesdell, aucune fonction d'énergie de déformation globale satisfaisant ces deux conditions n'a encore été identifiée.

L'essentiel

Imaginez que vous êtes un maître architecte concevant un bâtiment fait d'un matériau spécial et extensible. Votre objectif est d'écrire un ensemble de règles (une « loi de comportement ») qui prédit exactement comment ce matériau se comportera lorsqu'on le tire, le pousse ou le tord. Vous voulez être certain que vos règles ne prédisent jamais quelque chose d'impossible, comme un matériau qui se rétracterait soudainement en étant tiré plus fort, ou qui se comporterait de manière totalement différente simplement parce que vous avez fait pivoter le bâtiment.

Dans le monde de la physique, c'est le « Hauptproblem » (le Problème Principal) : comment écrire ces règles pour qu'elles soient mathématiquement saines et physiquement réalistes ?

Cet article explore deux ensembles de règles célèbres que les scientifiques ont proposés pour résoudre ce problème. Les auteurs, Wollner, Holzapfel et Neff, agissent comme des détectives testant ces règles les unes contre les autres. Ils se demandent : « Si un matériau suit la Règle A, suit-il automatiquement la Règle B ? »

Voici la décomposition de leur enquête en utilisant des analogies simples.

Les deux prétendants

1. La Polyconvexité (Le « Filet de sécurité mathématique »)
Considérez la Polyconvexité comme un filet de sécurité mathématique strict. C'est une règle qui garantit que le bâtiment ne s'effondrera pas dans un trou noir mathématique (où les solutions n'existent plus). Elle est très populaire dans les simulations informatiques car elle est facile à vérifier.

• La Promesse : Si vous utilisez cette règle, les mathématiques fonctionnent, et le matériau ne fera pas de choses étranges ou impossibles dans les équations.
• Le Piège : Les auteurs ont découvert que le simple fait qu'un matériau passe ce test du « filet de sécurité » ne signifie pas qu'il se comportera comme un matériau réel et sensé dans toutes les situations.

2. Le TSTS-M++ (Le « Sens commun de la monotonie »)
Considérez le TSTS-M++ (Monotonie de la contrainte vraie - déformation vraie) comme une règle de « Sens commun ». Elle stipule : « Si vous tirez plus fort sur le matériau, la force nécessaire pour le tirer doit continuer d'augmenter. Si vous le tordez davantage, la résistance doit continuer d'augmenter. » C'est comme étirer un élastique ; il devrait devenir plus difficile à étirer à mesure que l'on progresse, et non devenir soudainement plus facile.

• La Promesse : Cette règle garantit que le matériau se comporte de manière prévisible dans des tests spécifiques, comme l'étirement direct ou la torsion.
• Le Piège : Cette règle n'est pas non plus un remède miracle. Un matériau peut suivre cette règle et pourtant se comporter bizarrement d'autres manières.

L'enquête : Tester les règles

Les auteurs ont mis en place deux défis spécifiques pour voir si une règle pouvait remplacer l'autre.

Défi 1 : Le test d'étirement (Extension uniaxiale)

• Le Scénario : Imaginez que vous tirez sur un bloc de matériau pour l'étirer, comme du taffy (caramel mou).
• La Question : Si un matériau suit le « Filet de sécurité mathématique » (la Polyconvexité), deviendra-t-il toujours plus difficile de le tirer à mesure que vous l'étirez ?
• Le Résultat : Non. Les auteurs ont construit un modèle mathématique spécifique (un « faux matériau ») qui passait parfaitement le test de la Polyconvexité. Cependant, lorsqu'ils ont simulé l'étirement, la force nécessaire pour l'étirer augmentait, puis diminuait soudainement avant de remonter.
• L'Analogie : C'est comme une voiture qui est mathématiquement garantie d'être sûre, mais quand vous appuyez sur l'accélérateur, elle accélère, puis ralentit soudainement d'elle-même, puis accélère à nouveau. Ce n'est pas ainsi qu'une vraie voiture (ou un vrai matériau) devrait se comporter.
• Conclusion : La Polyconvexité seule ne suffit pas à garantir un comportement de « Sens commun » lors de l'étirement.

Défi 2 : Le test de torsion (Cisaillement simple)

• Le Scénario : Imaginez faire glisser le haut d'un jeu de cartes sur le côté tout en maintenant le bas immobile. C'est le « cisaillement ».
• La Question : Si un matériau suit la règle du « Sens commun » (TSTS-M++), sera-t-il toujours plus difficile de le tordre à mesure qu'on le tord davantage ?
• Le Résultat : Non. Les auteurs ont construit un autre « faux matériau » qui suivait parfaitement la règle du Sens Commun. Mais lorsqu'ils ont simulé la torsion, la résistance augmentait, chutait, puis remontait.
• L'Analogie : Imaginez une charnière de porte qui devient plus dure à pousser, puis devient soudainement lâche et facile à pousser, puis redevient dure. Cela viole le « Filet de sécurité mathématique » (plus précisément une condition appelée ellipticité de Legendre-Hadamard, qui assure la stabilité).
• Conclusion : Le Sens Commun (TSTS-M++) seul ne suffit pas à garantir la stabilité mathématique requise pour la torsion.

La vue d'ensemble : Le maillon manquant

Les auteurs concluent qu'aucune des deux règles n'est assez forte par elle-même.

• Vous avez besoin de la Polyconvexité pour garantir que les mathématiques sont stables (pas d'oscillations sauvages lors de la torsion).
• Vous avez besoin du TSTS-M++ pour garantir que le matériau se comporte de manière sensée lors de l'étirement (la force augmente toujours avec l'étirement).

L'Objectif Ultime : Le « Saint Graal » de ce domaine est de trouver un ensemble unique de règles qui satisfait les deux conditions en même temps pour toutes les déformations possibles.

• État Actuel : Les auteurs ont cherché très intensément à trouver ce « matériau parfait », mais ils n'ont pas trouvé de modèle qui fonctionne globalement (pour tous les étirements et torsions).
• Succès Partiel : Ils ont toutefois trouvé des solutions « limitées par une chaîne ». Considérez cela comme des matériaux qui se comportent parfaitement, mais seulement jusqu'à une certaine limite (comme un élastique qui fonctionne très bien jusqu'à une longueur spécifique, point auquel les règles se brisent).

Résumé pour le grand public

Cet article est un rappel à la réalité pour les scientifiques qui conçoivent des matériaux. Il dit : « Ne vous reposez pas uniquement sur un tour mathématique pour garantir que votre modèle de matériau est bon. »

• Si vous vérifiez seulement la Sécurité Mathématique (Polyconvexité), votre matériau pourrait agir bizarrement lors de l'étirement.
• Si vous vérifiez seulement le Sens Commun (TSTS-M++), votre matériau pourrait être instable lors de la torsion.

Pour véritablement résoudre le problème de la modélisation des matériaux élastiques idéaux, nous avons probablement besoin d'une combinaison de ces deux règles. Cependant, trouver une formule unique qui satisfait les deux parfaitement pour chaque situation possible reste un mystère non résolu, bien que les auteurs aient fourni de nouveaux outils et des réponses partielles pour aider les futurs chercheurs à percer le code.

Résumé technique

Résumé technique : À la recherche des conditions constitutives en hyperélasticité isotrope : polyconvexité versus monotonie de la contrainte réelle-déformation réelle

Énoncé du problème
L'article traite du « Hauptproblem » de l'élasticité finie : identifier les contraintes constitutives appropriées pour la fonction d'énergie de déformation $W$ afin de garantir un comportement matériel physiquement raisonnable dans le cas de l'hyperélasticité isotrope idéale. Si les mérites mathématiques de la polyconvexité sont bien établis (garantissant l'existence de minimiseurs et la quasi-convexité), ses conséquences mécaniques sont moins claires. À l'inverse, la condition de monotonie de la contrainte réelle-déformation réelle (TSTS-M++), récemment redécouverte, garantit des trajectoires contrainte-déformation monotones dans des modes de déformation spécifiques, mais manque des garanties de stabilité mathématique large de la polyconvexité. La question centrale est de savoir si l'une ou l'autre de ces conditions est suffisante à elle seule pour garantir un comportement physiquement raisonnable, ou si une combinaison est requise, et si une telle combinaison peut être réalisée globalement.

Méthodologie
Les auteurs emploient une approche analytique rigoureuse basée sur la théorie de l'hyperélasticité isotrope. Les éléments méthodologiques clés incluent :

• Représentation par invariants : La fonction d'énergie de déformation est paramétrée à l'aide d'un ensemble spécifique d'invariants $K_i$ (racines carrées des invariants principaux du tenseur de Cauchy-Green gauche $\mathbf{B}$), ce qui simplifie la représentation des inégalités constitutives par rapport aux élongations principales.
• Dérivation de conditions suffisantes : L'article dérive de nouvelles conditions suffisantes explicites pour la polyconvexité et la TSTS-M++ en termes des dérivées de la fonction d'énergie par rapport aux invariants $K_1, K_2, K_3$.
• Contre-exemples explicites : Pour tester la suffisance de ces conditions, les auteurs construisent des familles spécifiques de fonctions d'énergie de déformation. Ces fonctions sont conçues pour satisfaire une condition (par exemple, la polyconvexité) tout en violant explicitement le comportement physique attendu associé à l'autre (par exemple, la monotonie en extension uniaxiale).
• Preuves analytiques : L'étude repose sur la dérivation symbolique, la décomposition spectrale et l'analyse d'équations différentielles ordinaires (spécifiquement le Lemme 5.15) pour prouver les propriétés des fonctions d'énergie construites sans recourir à des calculs numériques de grande ampleur, excepté pour la visualisation.

Contributions clés et résultats
L'article apporte des réponses définitives à quatre « questions de défi » posées par Neff et al. (2024) concernant l'interaction entre la polyconvexité et la TSTS-M++ :

1. Insuffisance de la polyconvexité en extension uniaxiale (Défi ii & iii) :

   • Cas compressible : Les auteurs construisent une fonction d'énergie de déformation polyconvexe (Éq. 5.37) qui produit une réponse de la contrainte de Cauchy $\sigma_{11}$ non monotone en extension uniaxiale non contrainte. Plus précisément, la contrainte $\sigma_{11}$ atteint un maximum puis diminue à mesure que l'élongation $\lambda_1$ augmente, bien que la fonction soit polyconvexe et satisfasse la condition d'ellipticité de Legendre-Hadamard.
   • Cas incompressible : Bien qu'aucun contre-exemple n'ait été trouvé, les auteurs prouvent que toute fonction d'énergie de déformation incompressible satisfaisant les conditions suffisantes pour la polyconvexité proposées par Ball (1976) résulte automatiquement d'une réponse de contrainte réelle monotone. Cela réduit considérablement l'espace de recherche pour un contre-exemple potentiel, mais ne prouve pas l'impossibilité.

2. Insuffisance de la TSTS-M++ en cisaillement simple (Défi iv) :

   • Les auteurs construisent une fonction d'énergie de déformation (Éq. 5.69) qui satisfait la TSTS-M++ globalement mais présente une réponse de la contrainte de Cauchy non monotone en cisaillement simple.
   • Ce résultat implique que la fonction n'est pas à convexité de rang un (et donc pas polyconvexe), car la convexité de rang un est prouvée comme étant une condition nécessaire pour une contrainte de cisaillement monotone en cisaillement simple (Proposition 5.13).

3. Solutions limitées par des chaînes (Défi i) :

   • Les auteurs tentent de trouver une fonction unique satisfaisant à la fois la polyconvexité et la TSTS-M++ globalement. Ils ne parviennent pas à identifier une telle fonction définie sur l'ensemble du domaine $GL^+(3)$.
   • Cependant, ils construisent avec succès trois familles de fonctions d'énergie de déformation « limitées par des chaînes » (Propositions 5.2, 5.4, 5.5) qui satisfont les deux conditions dans des domaines de déformation restreints (par exemple, limités par des moyennes d'éléments de ligne, de surface ou de volume).

4. Nouvelles conditions suffisantes :

   • L'article établit des conditions suffisantes explicites pour la polyconvexité (Théorème 3.1) et la TSTS-M++ (Théorème 4.4) en termes des invariants $K_i$. Ces conditions soulignent les différences structurelles entre les deux contraintes, particulièrement en ce qui concerne les exigences de convexité de la fonction d'énergie et de ses dérivées.

Signification et affirmations
L'article conclut que ni la polyconvexité ni la TSTS-M++ ne sont suffisantes à elles seules pour garantir un comportement matériel physiquement raisonnable pour l'élasticité isotrope idéale.

• La polyconvexité assure la stabilité (ellipticité de Legendre-Hadamard) et la monotonie de la contrainte de cisaillement, mais échoue à garantir la monotonie de la contrainte en extension uniaxiale.
• La TSTS-M++ assure la monotonie de la contrainte en extension uniaxiale, mais échoue à garantir la monotonité de la contrainte de cisaillement (et donc la convexité de rang un).

Les auteurs suggèrent qu'une solution viable au Hauptproblem de Truesdell nécessite probablement une combinaison de ces deux conditions. Cependant, ils déclarent explicitement qu'aucune fonction d'énergie de déformation isotrope n'a été identifiée à ce jour qui satisfasse les deux contraintes globalement. La construction d'une telle fonction demeure un problème ouvert. Ils postulent que les conditions suffisantes dérivées pour les deux propriétés peuvent être mutuellement exclusives dans le cadre global, bien qu'elles puissent être réconciliées dans des domaines restreints (limités par des chaînes). Ce travail sert de mise en garde pour la modélisation constitutive, particulièrement pour les approches basées sur les données (par exemple, les réseaux de neurones) où la polyconvexité est souvent utilisée comme unique contrainte de stabilité, risquant ainsi de négliger l'exigence physique de trajectoires contrainte-déformation monotones en extension.