A 2D discrete model with a phase-field regularisation applied to dynamic crack propagation

Lorsqu'un matériau fragile comme le verre ou l'acier trempé se brise, il ne se contente pas de se fendre simplement en deux. Au lieu de cela, une fissure se propage à travers le solide, se ramifiant souvent en un réseau complexe de fractures qui fragmente l'objet en de nombreux morceaux. Prédire exactement où ces fissures voyageront et à quelle vitesse elles se déplaceront est l'un des défis les plus difficiles en ingénierie. Si un modèle utilisé pour simuler cette destruction est défectueux, il pourrait suggérer qu'un pont échouera dans une direction sûre alors qu'il s'effondrera de manière catastrophique. Pendant des décennies, les scientifiques ont tenté de construire des modèles informatiques qui imitent la structure atomique des matériaux, en les traitant comme une grille de petits points connectés par des ressorts. Bien que ces modèles soient puissants, ils souffrent d'un défaut singulier : le chemin emprunté par la fissure dépend souvent de la rotation de la grille informatique, plutôt que de la physique du matériau lui-même. C'est comme si une fissure dans un véritable morceau de verre changeait soudainement de direction simplement parce que l'ingénieur a pivoté l'écran de l'ordinateur.

Pour résoudre ce problème, Matías Braun, un chercheur de l'Universidad Nebrija en Espagne, a développé une nouvelle façon de simuler la fracture dynamique. Son approche combine une grille traditionnelle de points avec une méthode de « champ de phase » (phase-field), une technique qui traite la fissure non pas comme une ligne nette sautant d'un point à l'autre, mais comme une zone de dommage lisse et diffuse qui se propage à travers le matériau. En mélangeant ces deux idées, Braun a créé un modèle qui permet aux fissures de voyager dans n'importe quelle direction, indépendamment de l'orientation de la grille informatique sous-jacente. Il a testé cette nouvelle méthode sur trois scénarios distincts : une plaque d'acier frappée par un impact à haute vitesse, une plaque entaillée tirée jusqu'à la rupture, et un cylindre métallique explosant de l'intérieur vers l'extérieur. Dans chaque cas, le nouveau modèle a produit des résultats qui correspondent aux expériences réelles et a éliminé le biais artificiel causé par la grille informatique, offrant une vue plus claire et plus fiable de la façon dont les matériaux échouent réellement.

Le cœur du problème auquel Braun a été confronté réside dans la manière dont les ordinateurs gèrent le monde invisible des atomes. Dans ses travaux précédents, il utilisait un modèle discret où le matériau était représenté par un réseau triangulaire, un motif de points reliés par des lignes. Lorsque le matériau était soumis à des contraintes, les connexions entre les points se rompaient, simulant une fissure. Cependant, parce que les points étaient disposés selon un motif fixe, la fissure était forcée de suivre les lignes de la grille. Si la grille était pivotée, la fissure changeait de trajectoire, même si le matériau lui-même n'avait pas changé. Il s'agit d'un problème important car, en réalité, une fissure doit suivre le chemin de moindre résistance déterminé par les propriétés du matériau, et non l'angle arbitraire d'un maillage informatique.

La solution de Braun a été d'introduire une régularisation par champ de phase. Au lieu de forcer la rupture de connexions spécifiques entre les points, il a assigné une valeur à chaque point représentant son degré d'endommagement. Cette valeur varie de zéro, signifiant que le matériau est parfaitement intact, à un, signifiant qu'il est complètement brisé. La fissure devient une région où cette valeur augmente progressivement, créant une transition fluide plutôt qu'une ligne dentelée. Cela permet à la fissure de courber et de se ramifier dans n'importe quelle direction, ignorant efficacement la structure rigide de la grille. Crucialement, Braun a calibré ce système en utilisant deux propriétés fondamentales du matériau : sa résistance à la traction, qui est la force requise pour amorcer une fissure, et son énergie de fracture, qui est l'énergie nécessaire pour maintenir la croissance de la fissure. En ancrant la simulation à ces constantes physiques réelles, le modèle garantit que les résultats dépendent du matériau, et non des mathématiques.

Pour prouver l'efficacité de cette nouvelle méthode, Braun a mené trois tests rigoureux. Le premier fut le test Kalthoff–Winkler, une expérience classique où une plaque d'acier présentant une entaille pré-découpée est frappée par un projectile à grande vitesse. Dans ce scénario, une fissure jaillit de l'entaille selon un angle spécifique. Les modèles précédents prédisaient souvent le mauvais angle selon la rotation de la grille. Dans les simulations de Braun, la fissure voyageait systématiquement à un angle d'environ 70 degrés, correspondant parfaitement aux observations expérimentales. Ce résultat restait vrai que la grille soit alignée avec la plaque ou pivotée de 15 ou 30 degrés. La vitesse de la fissure restait également stable, voyageant à environ 70 % de la vitesse du son dans le matériau, indépendamment de l'orientation de la grille.

Le second test impliquait une plaque à bord entamé soumise à la traction, une configuration utilisée pour étudier la croissance lente des fissures sous une traction constante. Ici, l'objectif était de voir si la fissure suivrait une ligne droite de l'entaille jusqu'au bord opposé. Les modèles de grille traditionnels forçaient souvent la fissure à zigzaguer ou à dévier selon l'alignement de la grille. Le modèle de Braun, cependant, produisait une fissure droite pour chaque orientation de la grille. La force requise pour briser la plaque était également constante, variant de moins de deux pour cent quelle que soit la rotation de la grille. Cela a démontré que le modèle pouvait prédire avec précision à la fois le chemin et la résistance du matériau sans l'influence artificielle du maillage informatique.

Le dernier test, le plus complexe, concernait la fragmentation d'un cylindre épais soumis à une explosion interne. Dans ce scénario, aucune fissure unique n'est guidée par une entaille ; au contraire, de nombreuses fissures doivent naître spontanément et foncer vers l'extérieur, pulvérisant le cylindre en des dizaines de morceaux. C'est un événement chaotique où la symétrie de la grille pourrait facilement dicter le motif des fragments. Lorsque Braun a lancé la simulation avec la grille centrée sur l'axe du cylindre, les fissures formaient un motif symétrique. Cependant, lorsqu'il a légèrement décalé la grille pour qu'elle ne soit plus parfaitement centrée, les fissures se sont distribuées de manière irrégulière et le nombre de fragments a changé. Cela a confirmé que le modèle était sensible à la configuration physique plutôt qu'à la simple symétrie de la grille. La simulation a prédit entre 12 et 19 fragments, un chiffre qui se situe dans la fourchette d'autres méthodes avancées, mais qui est inférieur à certaines d'entre elles, probablement parce que la zone de dommage lisse du modèle empêche la formation de fissures trop proches les unes des autres.

Tout au long de ces tests, les chercheurs ont soigneusement suivi l'énergie du système. Ils ont comparé l'énergie perdue lors de la rupture du matériau à l'énergie requise pour créer les nouvelles surfaces de fracture. Dans certains cas, il existait un léger décalage entre l'énergie calculée par le modèle de dommage et l'énergie réellement libérée par la rupture des liaisons. Il s'agit d'une limitation connue de l'approche hybride, où la force motrice de la fissure est calculée différemment de l'énergie des liaisons rompues. Cependant, le bilan de l'énergie mécanique restait précis à moins de 0,5 %, confirmant que la simulation était physiquement cohérente. Les résultats suggèrent que, bien que le modèle ne soit pas parfait, il représente une étape significative dans la compréhension de la fracture dynamique.

Les implications de ce travail dépassent la simple amélioration des graphismes informatiques. En supprimant la dépendance à l'orientation de la grille, les ingénieurs peuvent avoir la certitude que leurs simulations de crash-tests, de défaillances structurelles ou de tests de matériaux reflètent la véritable physique du matériau. Le modèle réussit à combler le fossé entre la nature discrète des atomes et le comportement continu des objets solides. Il démontre qu'avec les bons outils mathématiques, il est possible de simuler le processus chaotique et violent de la rupture sans laisser la grille de l'ordinateur dicter le résultat. L'étude conclut que cette approche offre un moyen robuste de prédire comment les matériaux échouent dans des conditions extrêmes, offrant une voie plus claire pour concevoir des structures plus sûres et plus résilientes.