High-order isogeometric phase-field fracture modeling: a numerical-experimental study

Lorsque les ingénieurs conçoivent des ponts, des avions ou même des outils du quotidien, ils doivent anticiper la manière dont les matériaux vont faillir. La forme de défaillance la plus spectaculaire est la fracture, où une pièce solide se divise soudainement. Pendant des décennies, les scientifiques se sont appuyés sur deux manières principales pour étudier cela : les tests physiques en laboratoire et les simulations informatiques. Les tests physiques sont une preuve indéniable du comportement d'un matériau, mais ils sont coûteux, chronophages et destructifs. Les simulations informatiques offrent une alternative plus rapide et moins onéreuse, permettant aux concepteurs de tester des milliers de variations avant qu'une seule pièce de métal ne soit découpée. Cependant, un défi persistant a été l'écart entre les deux. De nombreux modèles informatiques sont construits pour décrire des matériaux qui se cassent comme du verre, alors que le monde réel est rempli de métaux qui se plient et s'étirent avant de rompre. Ce décalage force souvent les ingénieurs à choisir entre un modèle mathématiquement simple mais physiquement inexact, ou un modèle complexe mais difficile à exécuter.

Une équipe de chercheurs issus d'universités italiennes et allemandes a récemment exploré un juste milieu surprenant. Ils ont cherché à savoir si un modèle informatique conçu spécifiquement pour des matériaux fragiles, semblables au verre, pouvait prédire avec succès le comportement d'un métal ductile, capable de se déformer. Leur sujet était un alliage d'aluminium spécifique, couramment utilisé dans la construction et la fabrication, connu pour sa capacité à se déformer considérablement avant de rompre. Les chercheurs ont entrepris de voir s'ils pouvaient utiliser un cadre mathématique « fragile » pour simuler un événement réel « ductile », en se concentrant spéciflement sur la façon dont les fissures démarrent et se propagent dans un échantillon de métal présentant des entailles pré-découpées. Leur objectif n'était pas de recréer chaque minuscule détail de la flexion du métal, mais de voir si le modèle pouvait prédire avec précision deux choses critiques : la force maximale que le métal pourrait supporter avant de faillir, et le chemin exact que prendrait la fissure en déchirant le matériau.

L'expérience a débuté par des tests physiques sur des éprouvettes métalliques en forme d'haltère, une forme standard pour tester la résistance des matériaux. Les chercheurs ont utilisé un alliage d'aluminium connu sous le nom de EN AW-6060, qui contient du silicium et du magnésium. Ils ont préparé deux versions de l'éprouvette : une de 4 millimètres d'épaisseur et une autre de 6 millimètres. Chaque pièce présentait deux petites coupures triangulaires, ou encoches, positionnées légèrement excentrées par rapport au milieu. Ces encoches agissent comme des points faibles où une fissure est susceptible de débuter. L'équipe a placé ces éprouvettes dans une machine de test puissante qui les a tirées pour les séparer jusqu'à ce qu'elles se rompent. Ils ont observé que l'épaisseur du métal changeait radicalement le résultat. L'éprouvette plus épaisse de 6 millimètres a développé une fissure qui a voyagé diagonalement d'une encoche à l'autre, les reliant en une ligne unique et nette. En revanche, l'éprouvette plus mince de 4 millimètres a développé deux fissures distinctes qui se sont courbées l'une de l'autre et ne se sont jamais rejointes. Cette différence de comportement est un exemple classique de la façon dont les contraintes physiques de l'épaisseur d'un matériau peuvent altérer son mode de défaillance, un phénomène bien connu en ingénierie mais difficile à capturer dans des modèles informatiques simples.

Pour simuler ces événements, les chercheurs se sont tournés vers une méthode informatique sophistiquée appelée analyse isogéométrique, qui utilise des courbes mathématiques avancées pour représenter la forme de l'objet plus précisément que les méthodes traditionnelles basées sur des grilles. Ils ont employé un type spécifique de modèle appelé modèle de fracture par champ de phase. Dans cette approche, une fissure n'est pas représentée comme une ligne nette et dentelée que l'ordinateur doit suivre manuellement. Au lieu de cela, la fissure est traitée comme une zone lisse et floue où le matériau perd progressivement sa force. Cette variable de « champ de phase » agit comme un variateur d'intensité, abaissant la capacité du matériau à tenir ensemble, d'un état intact à un état complètement rompu. Les chercheurs ont choisi une version de ce modèle qui avait été conçue à l'origine pour des matériaux fragiles, ce qui signifie qu'elle suppose que le matériau se brise sans flexion ou étirement significatif. C'était un choix audacieux car l'aluminium testé est ductile, ce qui signifie qu'il s'étire et absorbe de l'énergie par déformation plastique avant de finir par se rompre.

L'équipe a exécuté ses simulations en utilisant les mêmes dimensions et propriétés de matériaux que les tests physiques. Ils ont ajusté un paramètre clé de leur modèle, appelé longueur interne, qui contrôle la largeur de la zone de dommage floue autour de la fissure. En ajustant soigneusement cette valeur, ils ont découvert qu'ils pouvaient reproduire les différents comportements observés en laboratoire. Lorsqu'ils ont fixé la longueur interne à une valeur plus grande, la simulation a produit un motif de fissure reliant les deux encoches en diagonale, correspondant à l'éprouvette épaisse de 6 millimètres. Lorsqu'ils ont réduit la longueur interne, la simulation a montré deux fissures distinctes et courbes qui ne se connectaient pas, reflétant parfaitement le comportement de l'éprouvette plus mince de 4 millimètres. Cela constituait une découverte significative car cela suggérait qu'un cadre mathématique unique et relativement simple pouvait capturer des phénomènes physiques complexes simplement en ajustant un seul paramètre, sans avoir besoin de modéliser explicitement la physique compliquée de l'étirement du métal.

Les résultats n'étaient pas seulement visuellement similaires ; ils étaient aussi numériquement précis. Le modèle informatique a prédit la force maximale que le métal pouvait supporter avant de rompre avec une précision remarquable. Pour l'éprouvette la plus épaisse, la force prédite était à moins de 1 % de la force réellement mesurée. Pour l'éprouvette la plus mince, la prédiction était à moins de 3,5 %. Ces chiffres sont exceptionnellement proches pour une simulation qui a intentionnellement ignoré la déformation plastique complexe qui se produit dans les métaux réels. Les chercheurs ont noté que, bien que le modèle ne puisse pas reproduire l'adoucissement lent et progressif du métal après avoir atteint son pic de résistance, il excellait à prédire le moment de la rupture et le chemin que prendrait la fissure. Cela fait de ce modèle un outil puissant pour les ingénieurs qui doivent connaître les limites de sécurité d'une conception avant qu'elle ne soit construite.

L'étude conclut que l'utilisation d'un modèle de fracture fragile pour simuler des matériaux ductiles peut être une stratégie de conception d'ingénierie valide et conservatrice. Par « conservatrice », les auteurs entendent que le modèle a tendance à prédire la défaillance à des charges légèrement inférieures à ce que le matériau peut réellement supporter, ce qui est une approche sûre pour les applications critiques de sécurité. Le succès de cette approche repose sur la formulation mathématique spécifique utilisée, qui permet à l'énergie d'être dissipée tôt dans le processus, imitant l'absorption d'énergie d'un matériau ductile même sans calculer explicitement l'étirement. Les chercheurs reconnaissent qu'il s'agit d'un compromis phénoménologique ; le modèle ne raconte pas toute l'histoire de ce qui se passe à l'intérieur du métal après son écoulement, mais il raconte la partie la plus importante : quand et comment il va se rompre.

Ce travail ouvre une nouvelle perspective sur la façon dont nous pouvons utiliser les simulations informatiques pour comprendre le monde physique. Il suggère que nous n'avons pas toujours besoin du modèle le plus complexe pour obtenir la bonne réponse. Parfois, un modèle plus simple, appliqué avec les bons outils mathématiques et une compréhension attentive de ses limites, peut fournir des prédictions précures et fiables. Pour l'alliage d'aluminium testé, les chercheurs ont démontré qu'un modèle fragile pouvait prédire avec succès à la fois la force requise pour rompre la pièce et le chemin complexe que suivrait la fissure, simplement en ajustant l'échelle de la zone de dommage. Cette découverte offre une voie pratique pour les ingénieurs concevant des pièces métalliques, leur permettant d'utiliser des simulations de haut niveau efficaces pour garantir la sécurité et la performance sans s'enliser dans le coût computationnel immense de la modélisation de chaque détail microscopique de la déformation plastique. L'étude témoigne de la puissance de la combinaison des tests rigoureux en laboratoire et des méthodes numériques avancées pour résoudre des défis d'ingénierie du monde réel.