Nonlinear Incompressible Shear Wave Models in Hyperelasticity and Viscoelasticity Frameworks, with Applications to Love Waves

Cet article présente des équations générales pour les ondes de cisaillement non linéaires dans des matériaux hyperélastiques et viscoélastiques incompressibles, les applique à la propagation d'ondes de type Love à une interface, et démontre par des simulations numériques que, bien que la vitesse variable de ces ondes satisfasse généralement la condition d'existence linéaire, elle tend vers la vitesse du matériau le plus rapide à long terme.

L'essentiel

🌊 Le Secret des Vagues qui ne se cassent pas (ou presque)

Imaginez que vous êtes un géologue ou un ingénieur en matériaux. Votre travail consiste à comprendre comment les ondes de choc (comme celles d'un tremblement de terre ou d'une explosion) voyagent à travers la Terre ou à travers des matériaux complexes comme le caoutchouc, les tissus biologiques ou les mousses.

Dans la physique classique, on imagine souvent ces matériaux comme des ressorts parfaits : si vous les étirez, ils reviennent à leur place exactement comme avant. C'est ce qu'on appelle l'élasticité linéaire. Mais dans la vraie vie, les choses sont plus compliquées. Si vous tirez trop fort sur un élastique ou si vous écrasez une mousse, elle ne se comporte pas comme un ressort simple. Elle devient "capricieuse", résistante, et parfois, elle s'échauffe. C'est ce qu'on appelle le comportement non linéaire.

Ce papier de recherche, écrit par Shawn, Samuel et Alexei, s'attaque à un problème précis : comment décrire mathématiquement ces ondes de cisaillement (des vagues qui font bouger les choses de côté, comme quand on secoue une serviette) dans des matériaux qui changent de forme de manière complexe.

Voici les trois grandes idées du papier, expliquées simplement :

1. La règle du "Tunnel Interdite" (Les Ondes de Love)

Imaginez deux couches de matériaux collées l'une sur l'autre.

• La couche du haut (comme la croûte terrestre) est un peu plus "molle".
• La couche du bas (comme le manteau terrestre) est plus "rigide".

Quand une onde traverse l'interface entre ces deux couches, elle a un comportement étrange. Si elle va trop vite, elle s'échappe dans le bas. Si elle va trop lentement, elle reste bloquée en haut. Mais il existe une vitesse magique (ni trop lente, ni trop rapide) où l'onde reste piégée à l'interface, voyageant comme un train sur un rail. C'est ce qu'on appelle une Onde de Love.

Les chercheurs ont découvert que même si le matériau est très complexe et non linéaire (il se déforme beaucoup), cette règle du "tunnel" reste vraie. L'onde trouve toujours sa vitesse idéale pour rester coincée entre les deux couches, un peu comme un skieur qui trouve la vitesse parfaite pour rester sur une pente sans tomber ni s'envoler.

2. Le Modèle du "Gâteau Yeoh" (La Recette Mathématique)

Pour décrire ces matériaux bizarres, les scientifiques utilisent des "recettes" mathématiques appelées fonctions d'énergie.

• Les anciens modèles utilisaient des recettes simples (comme le modèle de Mooney-Rivlin).
• Les auteurs de ce papier utilisent une recette plus sophistiquée appelée modèle de Yeoh.

Imaginez que le matériau est un gâteau.

• Le modèle simple dit : "Plus vous appuyez, plus c'est dur, de façon régulière."
• Le modèle de Yeoh dit : "Au début, c'est mou. Ensuite, ça durcit. Et si vous appuyez encore plus fort, ça devient dur comme de la pierre, mais avec des nuances subtiles."

En utilisant cette recette plus précise (qui inclut des termes cubiques et quintiques, c'est-à-dire des puissances mathématiques), ils ont pu écrire une nouvelle équation qui prédit comment l'onde se comporte quand le matériau est vraiment étiré ou compressé. C'est comme passer d'une carte routière simplifiée à un GPS en 3D qui tient compte des nids-de-poule et des virages serrés.

3. Le Sirop et la Mousse (La Viscoélasticité)

Il y a un autre ingrédient secret : la viscosité.

• Imaginez que vous secouez un élastique sec. Il vibre longtemps.
• Maintenant, imaginez que vous secouez un élastique trempé dans du miel. Il vibre, mais le miel absorbe l'énergie et l'arrête vite.

Les matériaux réels (comme les tissus humains ou les sols humides) agissent souvent comme cet élastique dans le miel. Ils dissipent l'énergie. Les chercheurs ont ajouté ce facteur "miel" (viscosité) à leurs équations.
Le résultat ? Les ondes se propagent toujours, mais elles s'atténuent plus vite et leur forme change. C'est crucial pour comprendre pourquoi les tremblements de terre ne détruisent pas tout à l'infini : le sol "avale" une partie de l'énergie.

4. La Simulation : Le Laboratoire Virtuel

Puisque ces équations sont trop compliquées pour être résolues à la main (comme essayer de calculer la trajectoire de chaque goutte d'eau dans une tempête), les auteurs ont utilisé des superordinateurs pour simuler des explosions virtuelles.

Ils ont créé un "monde virtuel" avec deux couches de matériaux et ont lancé une petite explosion (comme un ballon qui éclate). Ils ont observé comment les ondes se propageaient :

• Résultat surprenant : Même dans ce monde virtuel très complexe, les ondes finissent par se comporter de manière très prévisible. Elles s'organisent, trouvent leur vitesse de croisière à l'interface, et finissent par ressembler aux ondes simples que l'on connaissait déjà, mais avec une touche de réalisme en plus.

En résumé

Ce papier est comme un manuel de conduite pour les ondes sismiques dans des matériaux réels et complexes.

• Avant : On utilisait des règles simples qui fonctionnaient bien pour les petits mouvements, mais qui échouaient quand les choses devenaient extrêmes.
• Maintenant : Grâce à ce nouveau modèle (basé sur le modèle de Yeoh et la viscosité), nous avons une carte plus précise. Nous savons que même si le matériau se comporte de manière bizarre et non linéaire, les ondes de Love (celles qui voyagent le long des interfaces) restent fidèles à leurs règles fondamentales : elles trouvent toujours le chemin du milieu pour voyager efficacement.

C'est une avancée importante pour mieux prédire les dégâts des tremblements de terre, pour concevoir de meilleurs matériaux de protection, ou même pour comprendre comment les ondes se propagent dans nos propres corps (tissus mous).

Résumé technique

Résumé Technique : Modèles d'ondes de cisaillement non linéaires dans les milieux hyperélastiques et hyper-viscoélastiques

1. Problématique et Contexte

L'étude des ondes de cisaillement (ondes S) se propageant dans des matériaux élastiques est fondamentale en sismologie, en géologie et en science des matériaux. Bien que la théorie linéaire de l'élasticité (loi de Hooke) soit suffisante pour les petites déformations, elle échoue à décrire le comportement de matériaux complexes (polymères, tissus biologiques, roches profondes) soumis à de grandes déformations ou présentant des effets non linéaires.

L'article se concentre sur un problème spécifique : la modélisation des ondes de Love non linéaires. Ces ondes de surface polarisées horizontalement se propagent à l'interface entre deux milieux élastiques de propriétés mécaniques différentes (par exemple, la croûte terrestre et le manteau supérieur). Le défi principal est de généraliser la théorie linéaire existante pour inclure :

• Des déformations finies (non linéaires) via des modèles hyperélastiques.
• Des effets de dissipation d'énergie (viscoélasticité).
• La propagation dans des géométries à deux couches (couche supérieure sur un demi-espace).

2. Méthodologie

Les auteurs adoptent une approche rigoureuse basée sur la mécanique des milieux continus, combinant dérivation analytique et simulation numérique.

A. Cadre Théorique et Dérivation Analytique

• Formulation Lagrangienne : Le modèle est développé dans le cadre de l'hyperélasticité incompressible. Les équations du mouvement sont dérivées à partir des lois de conservation (masse, quantité de mouvement, moment cinétique) et d'une fonction de densité d'énergie de déformation $W^H$.
• Modèles de Matériaux :
  • Les auteurs se restreignent aux matériaux de type néo-Hookeen généralisé (où l'énergie ne dépend que du premier invariant $I_1$ du tenseur de déformation de Cauchy-Green) pour éviter des conditions de compatibilité complexes liées au second invariant $I_2$.
  • Ils utilisent spécifiquement le modèle cubique de Yeoh comme fonction d'énergie potentielle, permettant de capturer des comportements non linéaires jusqu'à l'ordre cubique et quintique.
• Viscoélasticité : Un potentiel dissipatif $W^V$ est introduit pour modéliser la viscosité, menant à un cadre hyper-viscoélastique.
• Équations Résultantes :
  • Pour l'hyperélasticité pure, ils dérivent une équation d'onde non linéaire générale contenant des termes polynomiaux cubiques et quintiques en fonction du gradient de déplacement.
  • Pour le cas viscoélastique, une équation d'onde généralisée incluant des termes de dispersion mixte (dérivées croisées comme $u_{xxt}$) est obtenue.

B. Simulation Numérique

• Méthode : Les auteurs utilisent la méthode des lignes (Method of Lines) pour discrétiser les équations aux dérivées partielles (EDP) en espace, transformant le problème en un système d'équations différentielles ordinaires (EDO) couplées, résolu avec MATLAB (ode23).
• Configuration : Une simulation 2D $(2+1)$ d'une couche isotrope sur un demi-espace isotrope.
• Conditions Initiales : Une « explosion » gaussienne soudaine est utilisée comme source, placée soit dans la couche supérieure, soit dans le demi-espace inférieur.
• Analyse : Suivi des fronts d'ondes à l'interface et à la surface libre, et calcul de la vitesse de phase en fonction du temps.

C. Analyse de Symétrie

• Pour le cas 1D, une analyse de symétrie de Lie est effectuée pour trouver des solutions exactes invariantes, bien que ces solutions soient souvent non bornées, elles servent de référence pour valider le comportement qualitatif des simulations numériques.

3. Résultats Clés

A. Équations et Modélisation

• Les auteurs ont établi de nouvelles équations d'ondes de cisaillement non linéaires pour des matériaux incompressibles, incluant explicitement les termes de viscosité et les non-linéarités d'ordre supérieur (cubiques et quintiques).
• Pour le modèle de Yeoh cubique, les équations contiennent des termes de dispersion complexes impliquant des dérivées mixtes de déplacement.

B. Comportement des Ondes de Love Non Linéaires

• Condition d'Existence : Contrairement au cas linéaire où la vitesse de phase $v$ doit strictement satisfaire $c_1 < |v| < c_2$ (avec $c_1$ et $c_2$ les vitesses des ondes de cisaillement dans les deux milieux), les simulations montrent que dans le cas non linéaire complet, la vitesse de l'onde d'interface et de surface satisfait généralement cette condition initialement.
• Évolution Temporelle : À long terme, la vitesse de l'onde tend vers la plus grande des deux vitesses matérielles ($\max(c_1, c_2)$), que ce soit pour les ondes d'interface ou de surface.
• Effet de la Viscoélasticité : L'introduction de la viscosité entraîne une dispersion plus rapide des ondes et une atténuation significative de l'amplitude, rendant le comportement plus prévisible et stable.

C. Comparaison Linéaire vs Non Linéaire

• Les simulations numériques confirment que le comportement qualitatif des ondes non linéaires (réflexion, réfraction, piégeage dans la couche supérieure) reste similaire au cas linéaire, mais avec des oscillations de haute fréquence et une amplitude modifiée par les termes non linéaires.
• L'étude des solutions 1D montre que les effets non linéaires et dispersifs s'atténuent avec le temps à mesure que l'amplitude diminue, ramenant le système vers un mode linéaire.

4. Contributions Principales

1. Développement Théorique : Dérivation rigoureuse d'équations d'ondes de cisaillement non linéaires pour des matériaux hyperélastiques et hyper-viscoélastiques incompressibles, valables pour une fonction d'énergie de déformation arbitraire, puis spécialisée au modèle de Yeoh.
2. Intégration Viscoélastique : Inclusion réussie des effets de dissipation d'énergie via un potentiel de viscosité, produisant des équations de mouvement complètes incluant des termes de dispersion mixte.
3. Validation Numérique : Mise en œuvre de simulations 2D haute résolution montrant que la condition d'existence classique des ondes de Love ($c_1 < |v| < c_2$) est respectée dans le régime non linéaire, avec une convergence vers la vitesse maximale du milieu à long terme.
4. Correction et Précision : L'article inclut un addendum corrigeant une erreur initiale concernant l'utilisation des invariants de déformation, précisant la restriction aux matériaux néo-Hookéens généralisés pour assurer la définition correcte de la pression hydrostatique.

5. Signification et Perspectives

Ce travail est significatif car il comble un vide entre la théorie linéaire des ondes sismiques et la réalité des grandes déformations dans les matériaux géologiques et biologiques.

• Applications Sismologiques : Le modèle offre un cadre plus réaliste pour l'analyse des ondes de surface générées par des tremblements de terre, en particulier près du foyer où les déformations sont importantes.
• Science des Matériaux : Il fournit des outils pour caractériser les propriétés non linéaires et viscoélastiques de matériaux complexes (tissus mous, polymères) via l'analyse des ondes de cisaillement.
• Futur : Les auteurs suggèrent d'étendre ce travail aux ondes anisotropes, aux ondes de Rayleigh, de Scholte et de Stoneley, ainsi qu'à l'analyse de la rupture d'ondes (wave breaking) observée numériquement dans des configurations radiales.

En conclusion, l'article établit une base mathématique solide et vérifiée numériquement pour l'étude des ondes de Love non linéaires, démontrant que les effets non linéaires et dissipatifs modifient l'amplitude et la dispersion sans nécessairement violer les conditions fondamentales de propagation d'interface établies par la théorie linéaire.