Thermodynamically consistent phase field model for hydrogen-assisted cracking

Cet article présente un modèle de champ de phase thermodynamiquement cohérent qui simule la fissuration assistée par l'hydrogène dans les matériaux polycristallins en couplant la propagation des fissures avec la ségrégation de l'hydrogène et la réduction de l'énergie interfaciale, capturant avec succès la transition d'une rupture transgranulaire à une rupture intergranulaire sous les mécanismes de décohésion induite par l'hydrogène.

L'essentiel

Imaginez que vous possédez une structure métallique très solide, comme le châssis d'une voiture ou un pont. On s'attendrait à ce qu'elle résiste à la pression, mais parfois, des atomes d'hydrogène invisibles se glissent à l'intérieur du métal et provoquent sa rupture soudaine. Ce phénomène est appelé fragilisation par l'hydrogène. C'est comme si le métal était secrètement « empoisonné » de l'intérieur, le rendant cassant et sujet à la rupture.

Les scientifiques ont tenté de construire des modèles informatiques pour prédire exactement comment et où ce métal va se briser. Cependant, les modèles précédents présentaient une faille majeure : ils traitaient le comportement de l'hydrogène comme une règle simple et uniforme qui s'applique partout, alors que la structure interne du métal est en réalité un patchwork complexe de différents grains et de joints de grains.

Le nouveau modèle « intelligent »
Les auteurs de cet article ont créé une nouvelle simulation informatique plus sophistiquée (appelée « modèle à champ de phase ») qui agit comme une carte haute définition, thermodynamiquement cohérente. Voici comment elle fonctionne, en utilisant des analogies de la vie quotidienne :

• Le métal comme une foule : Imaginez que le métal est une pièce bondée remplie de gens (les atomes de métal). Les « joints de grains » sont les lignes invisibles séparant différents groupes de personnes. La « fissure » est un écart grandissant dans la foule.
• L'hydrogène comme un invité collant : Les atomes d'hydrogène sont comme des invités collants qui adorent se cacher dans les espaces vides entre les gens. Ils ont une préférence particulière : ils aiment encore plus se coller aux bords de la fissure et aux lignes entre les groupes (les joints de grains) qu'au milieu de la foule.
• Le problème de la « colle » : Dans un métal sain, la « colle » qui maintient les bords de la fissure ensemble est forte. Mais lorsque ces invités collants d'hydrogène se rassemblent aux bords de la fissure, ils agissent comme une huile glissante, affaiblissant la colle. Cela rend l'ouverture de la fissure beaucoup plus facile.
• L'ancienne approche vs la nouvelle approche :
  • Anciens modèles : Ils utilisaient un manuel de règles générique (l'isotherme de Langmuir-McLean) qui supposait que l'hydrogène était réparti uniformément et en équilibre parfait partout. C'est comme supposer que tout le monde dans la pièce bondée est immobile et espacé de manière égale, ce qui n'est pas le cas lors de la formation d'une fissure.
  • Nouveau modèle : Il utilise un cadre « variationnel » flexible (basé sur le formalisme de Kim-Kim-Suzuki). Au lieu d'imposer une règle rigide, il laisse l'hydrogène « migrer » naturellement là où il veut aller (les bords de la fissure et les joints de grains) en fonction des conditions locales. Il calcule précisément comment la « colle » s'affaiblit en temps réel à mesure que l'hydrogène s'accumule.

Ce qu'ils ont découvert
L'équipe a testé son nouveau modèle avec deux scénarios principaux :

1. Le test de la fissure unique : Ils ont simulé une fissure dans une seule pièce de métal. Sans hydrogène, la fissure a grandi exactement comme la physique le prédit (suivant le critère de Griffith). En ajoutant l'hydrogène, le modèle a montré que la fissure se propageait beaucoup plus facilement car l'hydrogène avait affaibli l'énergie de surface. Les résultats correspondaient parfaitement aux prédictions théoriques, prouvant ainsi l'efficacité du modèle.

2. Le test polycristallin (La grande découverte) : Ils ont simulé un métal composé de nombreux petits cristaux (grains) avec des joints entre eux.

   • Sans hydrogène : La fissure préférait traverser les grains de plein fouet (rupture transgranulaire). C'était comme une boule de démolition qui traverse les murs d'une maison parce que les murs sont plus fragiles que le mortier entre eux.
   • Avec l'hydrogène : L'hydrogène s'est fortement accumulé aux joints entre les grains, affaiblissant considérablement le « mortier » plus que les « murs ». Soudain, la fissure a changé de trajectoire. Au lieu de traverser les grains, elle a commencé à serpenter le long des joints (rupture intergranulaire). C'était comme si l'hydrogène avait transformé le mortier en sable mouillé, faisant s'effondrer la maison le long des coutures plutôt qu'à travers les briques.

Pourquoi cela importe
Ce nouveau modèle est une amélioration significative car il ne se contente pas de deviner où va l'hydrogène ; il le calcule en fonction de la thermodynamique réelle du système. Il parvient à capturer la transition d'un type de fissuration à un autre, ce qui est crucial pour comprendre pourquoi les matériaux échouent en présence d'hydrogène.

Les auteurs notent que bien que ce modèle soit une étape majeure, il se concentre actuellement sur un mécanisme spécifique (l'affaiblissement de la colle par l'hydrogène). Les travaux futurs devront ajouter d'autres facteurs complexes, comme la façon dont le métal se déforme et se tord (plasticité) et comment d'autres types de défauts interagissent avec l'hydrogène. Mais pour l'instant, ce modèle offre une manière claire, cohérente et précise de voir comment l'hydrogène transforme un métal solide en un matériau fragile.

Résumé technique

Énoncé du problème
La fragilisation par l'hydrogène (HE) représente un défi critique pour les technologies de l'hydrogène, provoquant une dégradation matérielle catastrophique et des ruptures imprévues dans les systèmes énergétiques décarbonés. Bien que des mécanismes tels que la décohésion induite par l'hydrogène (HEDE), la plasticité localisée induite par l'hydrogène (HELP) et les lacunes induites par l'hydrogène (HESIV) aient été proposés, la prédiction de la fragilisation par l'hydrogène reste difficile en raison de sa nature multiphysique s'étendant des échelles atomiques aux échelles macroscopiques. Les modèles de champ de phase (PF) existants pour la fissuration assistée par l'hydrogène couplent généralement la mécanique de rupture fragile avec le mécanisme HEDE via la diffusion de l'hydrogène. Cependant, ces modèles souffrent de deux limitations primaires : ils ne reproduisent pas explicitement la ségrégation de l'hydrogène (le phénomène responsable de la réduction de l'énergie de surface) et reposent sur l'isotherme de Langmuir-McLean pour calculer la couverture en hydrogène. Cette dernière est une relation d'équilibre macroscopique entre des phases homogènes de volume et de surface, ce qui n'est pas physiquement justifié comme une loi locale pour modéliser des microstructures hétérogènes impliquant des interactions élastiques.

Méthodologie
Les auteurs proposent un modèle de champ de phase thermodynamiquement cohérent basé sur le formalisme de Kim-Kim-Suzuki (KKS) pour simuler la fissuration assistée par l'hydrogène dans les matériaux polycristallins. Le modèle opère dans un cadre variationnel où l'énergie libre totale inclut des contributions de l'énergie de surface de la fissure, de l'énergie élastique, du couplage des joints de grains (GB) et de l'énergie libre chimique.

Les caractéristiques méthodologiques clés incluent :

• Variables de champ : Le système est décrit par un champ de dommage ($\eta$), une probabilité d'occupation de l'hydrogène ($c$) et un ensemble de paramètres d'ordre ($\phi_i$) représentant les orientations cristallographiques.
• Formalisme KKS : Au lieu d'imposer une loi de Langmuir-McLean locale, le modèle définit des phases distinctes de volume ($c_b$) et de surface de fissure ($c_{ck}$). La densité d'énergie libre chimique est construite comme une somme pondérée des contributions de volume et de fissure, contrainte par l'égalité des potentiels chimiques ($\mu_b = \mu_{ck}$). Cela permet de découpler les propriétés interfaciales des champs de concentration.
• Cohérence thermodynamique : Le modèle suppose des solutions solides idéales pour l'énergie chimique et dérive les équations d'équilibre régissantes (équilibre mécanique et équilibre du potentiel chimique) à partir du fonctionnel de l'énergie libre totale.
• Implémentation numérique : Les équations d'équilibre sont résolues à l'aide de méthodes de Fourier-spectre et d'un solveur FFT discret pour l'équilibre mécanique.

Contributions clés

1. Modélisation explicite de la ségrégation : Le modèle rend compte naturellement de la ségrégation de l'hydrogène sur les surfaces de fissures et les joints de grains sans supposer une isotherme macroscopique comme loi locale.
2. Dérivation variationnelle de la dégradation de l'énergie de surface : La dégradation effective de l'énergie de surface de la fissure induite par la ségrégation de l'hydrogène est calculée analytiquement au sein du cadre thermodynamique.
3. Capture de la transition : Le modèle démontre sa capacité à capturer la transition d'une fissuration transgranulaire vers une fissuration intergranulaire assistée par l'hydrogène via le mécanisme HEDE.

Résultats

• Validation dans les monocristaux : Les simulations numériques d'une fissure préexistante dans un monocristal sous chargement de mode I ont montré que la contrainte critique pour la croissance de la fissure ($\sigma_c$) suit le critère de Griffith. Lorsque le mécanisme HEDE était activé, $\sigma_c$ diminuait en raison d'une réduction de l'énergie de surface effective ($\Delta \gamma_{ch} \approx -4,29$ J/m$^2$), ce qui est cohérent avec les prédictions analytiques dérivées des équations du modèle.
• Application polycristalline : Les simulations dans un matériau polycristallin à base de nickel ont révélé une transition distincte dans les chemins de propagation des fissures. En l'absence d'hydrogène, la fissure se propageait de manière transgranulaire en raison d'une ténacité élevée des joints de grains. Lors de l'introduction de l'hydrogène, le modèle a prédit un passage à une fissuration intergranulaire. Cela a été attribué à une réduction drastique du rapport d'énergie de rupture joint de grain/volume (de 0,5 à 0,07), causée par des énergies de ségrégation de l'hydrogène plus élevées aux surfaces de type joint de grain par rapport aux surfaces de type volume. La carte de concentration d'hydrogène a confirmé une ségrégation significative aux joints de grains et aux surfaces de grains brisés.

Signification
L'article affirme que ce travail fournit une méthode rigoureuse et générique pour modéliser le mécanisme HEDE à l'échelle microstructurale. En utilisant le formalisme KKS, le modèle reproduit la ségrégation de l'hydrogène et la réduction de l'énergie de surface résultante sans dépendre de l'isotherme de Langmuir-McLean, qui est limitée thermodynamiquement pour les situations hétérogènes. L'étude valide que l'isotherme n'émerge naturellement que dans des limites spécifiques, homogènes et statiques, justifiant ainsi la nouvelle approche pour les microstructures complexes. Le modèle capture avec succès la transition observée expérimentalement d'une fissuration transgranulaire à une fissuration intergranulaire en utilisant des données d'entrée réalistes pour le nickel, offrant un outil robuste pour comprendre la défaillance assistée par l'hydrogène dans les métaux polycristallins. Les auteurs notent que les travaux futurs traiteront l'anisotropie de l'énergie de surface et l'incorporation d'autres mécanismes de fragilisation par l'hydrogène, tels que le mécanisme HELP, via la modélisation des dislocations.