Direct stress imaging from shear wave propagation

Résumé technique : Imagerie directe des contraintes par propagation d'ondes de cisaillement

Énoncé du problème
L'imagerie quantitative des champs de contraintes dans les solides hétérogènes demeure un défi majeur en mécanique et en biomécanique. Contrairement au déplacement, qui est directement observable, la contrainte est une variable d'état tensorielle qui ne peut être mesurée directement. Les approches conventionnelles infèrent généralement la contrainte de manière indirecte à partir des déformations mesurées en utilisant des modèles constitutifs et des hypothèses concernant les propriétés des matériaux et les configurations de référence. Ces méthodes nécessitent souvent un étalonnage matériel extensif et peinent à s'appliquer aux matériaux mous, hétérogènes ou biologiquement pertinents. De plus, les techniques non destructives existantes basées sur l'effet acoustoélastique (où une contrainte préexistante modifie la propagation des ondes) reposent généralement sur des mesures de vitesse d'onde locales, des constantes élastiques d'ordre supérieur étalonnées ou des hypothèses simplifiées sur la géométrie et l'homogénéité. Par conséquent, la cartographie à haute résolution des champs de contraintes internes, particulièrement dans les matériaux complexes, est restée largement inaccessible.

Méthodologie : Imagerie Acoustoélastique (IAE)
Les auteurs introduisent l'Imagerie Acoustoélastique (IAE), un cadre non destructif conçu pour reconstruire les champs de contraintes directement à partir de la propagation des ondes de cisaillement, sans spécification explicite de modèle constitutif ni étalonnage de paramètres matériels.

• Base théorique : La méthode exploite l'effet acoustoélastique, où une contrainte préexistante modifie la réponse élastique incrémentale d'un matériau. L'équation d'onde directrice pour les déplacements infinitésimaux dans une configuration précontrainte contient des coefficients d'élasticité directement liés à l'état de contrainte. Plus précisément, dans le système de coordonnées principales, la différence entre les contraintes principales ($\sigma_1 - \sigma_2$) est encodée dans la différence entre des composantes spécifiques du tenseur d'élasticité eulérien ($A^0_{1'2'1'2'} - A^0_{2'1'2'1'}$).
• Formulation du problème inverse : L'IAE traite la récupération de la contrainte comme un problème inverse. Au lieu d'estimer les vitesses d'onde locales, le cadre cherche à reconstruire les coefficients spatialement variables des équations d'ondes directrices à partir de mesures de champs d'ondes complets.
• Apprentissage profond informé par la physique : Pour résoudre ce problème inverse de manière efficace, les auteurs emploient l'inversion de forme d'onde de cisaillement complète (FSWI) rendue possible par des réseaux de neurones informés par la physique (PINNs).
  • L'architecture utilise six réseaux de neurones entièrement connectés pour représenter les champs inconnus comme des fonctions continues de l'espace. Un réseau prédit le champ de déplacement (paramétré via une fonction de flux pour satisfaire l'incompressibilité) et le multiplicateur de Lagrange, tandis que les réseaux restants infèrent les directions de contraintes principales et les coefficients de l'équation d'onde.
  • La fonction de perte combine un terme piloté par les données (minimisant l'écart entre les champs de déplacement prédits et mesurés) et un terme informé par la physique (imposant les équations d'ondes acoustoélastiques directrices).
  • L'inversion met à jour de manière itérative la différence de contrainte ($\sigma_1 - \sigma_2$) et la direction principale ($\theta$) afin de minimiser les résidus.

Contributions clés et résultats
Les auteurs valident l'IAE par des simulations numériques et des mesures expérimentales d'élastographie par ondes de cisaillement par ultrasons (SWE).

1. Reconstruction de contraintes non uniformes et résiduelles :

   • Dans des simulations numériques impliquant un matériau présentant des contraintes résiduelles spatialement non uniformes (générées par une expansion thermique contrainte), l'IAE a réussi à reconstruire à la fois l'amplitude de la contrainte et les directions principales.
   • La méthode a démontré sa capacité à séparer la contrainte de la rigidité du matériau. Dans des matériaux hétérogènes contenant une inclusion rigide au sein d'une matrice plus molle, l'IAE a récupéré avec précision le champ de contrainte malgré l'hétérogénéité élastique. En revanche, l'SWE conventionnelle interprète les changements de vitesse induits par la contrainte comme des variations de rigidité, conduisant à des cartes de module lissées.
   • Il a été constaté que la reconstruction des coefficients d'élasticité individuels ($\alpha$ et $\gamma$) est non unique (différant par un décalage uniforme), mais que leur différence, qui détermine la contrainte, est récupérable de manière unique.

2. Résolution sous-longueur d'onde :

   • Une découverte notable est l'obtention d'une résolution spatiale sous-longueur d'onde (< 0,28$\lambda$). En utilisant le champ d'onde complet, incluant les composantes de champ proche et les distorsions de front d'onde, l'IAE résout des caractéristiques de contrainte nettement plus petites que la longueur d'onde de cisaillement.
   • Cette résolution s'est révélée largement indépendante de la fréquence d'excitation (comparaison entre 6 kHz et 12 kHz), contrastant avec l'imagerie conventionnelle de champ lointain limitée par la diffraction. Cependant, les auteurs notent que l'utilisation de données de champ d'onde partielles (par exemple, une seule composante de déplacement) dégrade la résolution, bien que des fréquences plus élevées puissent partiellement compenser.

3. Validation expérimentale :

   • Le cadre a été implémenté à l'aide d'une plateforme standard de SWE par ultrasons sur un échantillon d'hydrogel de PVA. Des ondes de cisaillement ont été générées par un actionneur PZT et enregistrées via une imagerie ultrasonore ultra-rapide.
   • Malgré de faibles distorsions apparentes dans le champ d'onde (augmentation de la vitesse moyenne de 4,3 m/s à 4,5 m/s), l'IAE a reconstruit des champs de contrainte et des directions principales qui concordent quantitativement avec les solutions de référence par éléments finis. La méthode s'est montrée robuste au bruit expérimental, maintenant sa précision pour des rapports signal sur bruit dépassant environ 30 dB.

Signification et affirmations
Les auteurs positionnent l'IAE comme un cadre général pour l'imagerie de contrainte à haute résolution qui déplace le paradigme de l'imagerie cinématique (mesure du déplacement) vers l'imagerie de variables d'état tensorielles (mesure de la contrainte).

• Mesure directe : Ce travail établit que la contrainte peut être traitée comme une variable d'état mécanique mesurable, encodée directement dans la propagation des ondes, indépendamment de lois de comportement spécifiques.
• Complémentarité : L'IAE et l'élastographie conventionnelle sont présentées comme des modalités complémentaires ; tandis que l'élastographie caractérise la rigidité du matériau, l'IAE caractérise l'état mécanique (la contrainte) du matériau.
• Applicabilité : La méthode offre une voie vers la cartographie non invasive des états mécaniques internes dans des matériaux complexes et des tissus biologiques.
• Limites : Les auteurs reconnaissent modestement les limites actuelles, notamment l'hypothèse de variations de contrainte relativement lisses et l'utilisation de représentations de champs d'ondes simplifiées. Ils notent que la validation expérimentale est actuellement limitée à des environnements contrôlés et que l'extension du cadre aux milieux compressibles et à des non-linéarités plus complexes est une étape future nécessaire.

En résumé, l'article démontre qu'en combinant l'inversion de forme d'onde complète et l'apprentissage informé par la physique, il est possible de reconstruire des champs de contraintes internes avec une haute résolution spatiale et sans étalonnage matériel préalable, offrant ainsi un nouvel outil pour comprendre le comportement mécanique des matériaux mous et des tissus biologiques.