🪨 Propagation de fissures et mécanique de la rupture

Le facteur d'intensité de contrainte K = σ√(πa) est l'amplitude de la singularité de contrainte en 1/√r en pointe de fissure — il condense la contrainte appliquée et la géométrie de la fissure en un seul nombre. La ténacité à la rupture Kc = √(2Eγ) = √(E·Gc) est une propriété du matériau, indépendante de la fissure ou de la charge. La rupture survient lorsque K atteint Kc, ce qui est mathématiquement identique à σ atteignant σc, mais c'est la forme utilisée dans toute l'ingénierie structurelle réelle car K et Kc peuvent être mesurés directement et comparés pour n'importe quelle forme de fissure ou chargement, pas seulement la fissure centrale idéalisée résolue par Griffith.

Comme σc diminue à mesure que la fissure grandit (σc ∝ 1/√a) tandis que la contrainte appliquée reste fixe, l'excès de contrainte au-delà du seuil critique augmente continuellement à mesure que la fissure s'allonge. Cette boucle de rétroaction positive explique pourquoi, une fois qu'une fissure dépasse la longueur critique à une contrainte donnée, elle ne progresse pas lentement — elle accélère vers la vitesse du son dans le matériau en quelques microsecondes, ce qui est exactement ce que reproduit le modèle de croissance de cette simulation lorsque la demi-longueur de la fissure dépasse le point de croisement de σc.

Les métaux peuvent se déformer plastiquement en pointe de fissure : les dislocations glissent et émoussent la pointe aiguë, répartissant la charge sur une plus grande région et consommant une grande quantité d'énergie sous forme de chaleur et de dommages au réseau plutôt que de simplement rompre des liaisons. Cela élève leur énergie de rupture effective (souvent modélisée comme un « γ effectif » de centaines à milliers de J/m², contre environ 1 J/m² pour l'énergie de surface réelle du verre) même si l'énergie de surface réelle du métal lui-même est d'un ordre de grandeur similaire. Les céramiques et les verres, dépourvus de mouvement de dislocation facile à température ambiante, restent fragiles et ne comptent que sur leur énergie de surface réelle.

Une tessellation de Voronoï divise une surface en régions, chacune contenant tous les points les plus proches d'un point germe particulier. Les motifs réels d'éclatement fragile — vitres, carreaux de céramique, fragmentation de roche — ressemblent étroitement à des cellules de Voronoï car les surfaces de rupture ont tendance à se propager le long de chemins qui divisent le matériau en régions à peu près convexes autour des sites de défaut ou d'impact. Cette simulation disperse des points germes (biaisés vers la ligne de fissure) et attribue chaque triangle du maillage de la plaque à son germe le plus proche, produisant un motif de fragments physiquement motivé plutôt qu'une rupture arbitraire en grille.

Les avions, les réservoirs sous pression, les pipelines et les ponts sont conçus en utilisant des limites d'inspection issues de la mécanique de la rupture : les ingénieurs mesurent ou supposent une taille maximale de défaut non détecté, calculent le facteur d'intensité de contrainte à la charge de conception, et exigent qu'il reste en toute sécurité en dessous de la ténacité à la rupture certifiée Kc du matériau. Cette philosophie de « conception tolérante aux dommages », issue directement des expériences sur fibres de verre de Griffith dans les années 1920 et des catastrophes de fissures de fatigue du Comet dans les années 1950, explique pourquoi les avions subissent des inspections de fissures planifiées plutôt que d'être conçus pour ne jamais se fissurer.

A. A. Griffith a publié la théorie du bilan énergétique de la rupture fragile en 1920, en la testant sur des fibres de verre. Sa théorie sous-estimait la ténacité des métaux ductiles jusqu'à ce que G. R. Irwin l'étende dans les années 1950 pour inclure le champ de contrainte en pointe de fissure (le K utilisé dans cette simulation) et une énergie de surface effective tenant compte du travail plastique, fondant la mécanique linéaire élastique de la rupture moderne (MLER). Les travaux d'Irwin ont directement suivi les enquêtes sur les accidents du de Havilland Comet en 1954, causés par des fissures de fatigue se développant à partir des découpes carrées des hublots.

La modélisation de rupture par champ de phase — simulant l'initiation et la propagation des fissures comme un champ d'endommagement lisse et évolutif plutôt qu'en suivant une géométrie explicite de pointe de fissure — est un domaine de recherche majeur, car elle peut naturellement gérer la ramification des fissures, leur fusion et des surfaces de rupture 3D complexes très difficiles pour la MLER classique. Elle est utilisée pour étudier tout, des réseaux de fractures rocheuses induites par la fracturation hydraulique à la fatigue dans les pièces métalliques fabriquées de manière additive (impression 3D), où la porosité interne crée des populations de défauts imprévisibles que l'analyse classique de Griffith à fissure unique ne peut pas capturer.