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

Wenn ein sprödes Material wie Glas oder gehärteter Stahl bricht, reißt er nicht einfach in zwei Hälften. Stattdessen rast ein Riss durch den Festkörper und verzweigt sich oft in ein komplexes Geflecht von Brüchen, das das Objekt in viele Teile zersplittern lässt. Genau vorherzusagen, wohin diese Risse verlaufen werden und wie schnell sie sich bewegen, gehört zu den schwierigsten Herausforderungen im Ingenieurwesen. Wenn ein Modell, das diese Zerstörung simulieren soll, fehlerhaft ist, könnte es suggerieren, dass eine Brücke in eine sichere Richtung versagt, während sie sich tatsächlich katastrophal zusammenstürzt. Seit Jahrzehnten versuchen Wissenschaftler, Computermodelle zu entwickeln, die die atomare Struktur von Materialien nachahmen, indem sie diese als ein Gitter aus winzigen Punkten behandeln, die durch Federn verbunden sind. Obwohl diese Modelle leistungsstark sind, leiden sie unter einem merkwürdigen Fehler: Der Pfad, den der Riss nimmt, hängt oft davon ab, wie das Computergitter rotiert ist, und nicht von der Physik des Materials selbst. Es ist, als ob ein Riss in einem echten Stück Glas plötzlich die Richtung ändern würde, nur weil der Ingenieur den Computerbildschirm gedreht hat.

Um dieses Problem zu lösen, hat Matías Braun, ein Forscher an der Universidad Nebrija in Spanien, einen neuen Weg entwickelt, um dynamische Fraktur zu simulieren. Sein Ansatz kombходит ein traditionelles Punktgitter mit einer „Phasenfeld“-Methode – einer Technik, die den Riss nicht als scharfe Linie behandelt, die von Punkt zu Punkt springt, sondern als eine glatte, diffuse Zone der Schädigung, die sich durch das Material ausbreitet. Durch die Verschmelzung dieser beiden Ideen schuf Braun ein Modell, das es Rissen ermöglicht, in jede beliebige Richtung zu verlaufen, unabhängig davon, wie das zugrunde liegende Computergitter orientiert ist. Er testete diese neue Methode in drei verschiedenen Szenarien: einer Stahlplatte, die von einem Hochgeschwindigkeitsaufprall getroffen wurde, einer genutten Platte, die auseinandergezogen wurde, bis sie brach, und einem dicken Metallzylinder, der von innen heraus explodierte. In jedem Fall lieferte das neue Modell Ergebnisse, die mit realen Experimenten übereinstimmten und den künstlichen Bias durch das Computergitter eliminierten, was eine klarere, zuverlässigere Sichtweise darauf bot, wie Materialien tatsächlich versagen.

Der Kern des Problems, das Braun adressierte, liegt darin, wie Computer die unsichtbare Welt der Atome handhaben. In seiner früheren Arbeit verwendete er ein diskretes Modell, bei dem das Material als ein dreieckiges Gitter, ein Muster aus Punkten, die durch Linien verbunden sind, dargestellt wurde. Wenn das Material unter Spannung gesetzt wurde, brachen die Verbindungen zwischen den Punkten, was einen Riss simulierte. Da die Punkte jedoch in einem festen Muster angeordnet waren, wurde der Riss gezwungen, den Linien des Gitters zu folgen. Wenn das Gitter rotiert wurde, änderte der Riss seinen Pfad, obwohl sich das Material selbst nicht verändert hatte. Dies ist ein erhebliches Problem, denn in der Realität sollte ein Riss dem Pfad des geringsten Widerstands folgen, der durch die Materialeigenschaften bestimmt wird, und nicht durch den willkürlichen Winkel eines Computergitters.

Brauns Lösung bestand darin, eine Phasenfeld-Regularisierung einzuführen. Anstatt den Riss dazu zu zwingen, spezifische Verbindungen zwischen den Punkten zu brechen, wies er jedem Punkt einen Wert zu, der repräsentiert, wie beschädigt er ist. Dieser Wert reicht von Null, was bedeutet, dass das Material perfekt intakt ist, bis zu Eins, was bedeutet, dass es vollständig gebrochen ist. Der Riss wird zu einer Region, in der dieser Wert allmählich zunimmt, was einen glatten Übergang statt einer gezackten Linie schafft. Dies ermöglicht es dem Riss, in jede beliebige Richtung zu krümmen und sich zu verzweigen, wodurch er die starre Struktur des Gitters effektiv ignoriert. Entscheidend ist, dass Braun dieses System mit zwei grundlegenden Materialeigenschaften kalibrierte: der Zugfestigkeit, welche die Kraft ist, die erforderlich ist, um einen Riss zu starten, und der Bruchenergie, welche die Energie ist, die benötigt wird, um den Riss wachsen zu lassen. Durch die Verankerung der Simulation an diesen realen physikalischen Konstanten stellt das Modell sicher, dass die Ergebnisse vom Material abhängen und nicht von der Mathematik.

Um zu beweisen, dass diese neue Methode funktionierte, führte Braun drei strenge Tests durch. Der erste war der Kalthoff–Winkler-Test, ein klassisches Experiment, bei dem eine Stahlplatte mit einer vorgefertigten Kerbe von einem Projektil mit hoher Geschwindigkeit getroffen wird. In diesem Szenario schießt ein Riss aus der Kerbe in einem bestimmten Winkel heraus. Frühere Modelle sagten oft den falschen Winkel voraus, je nachdem, wie das Gitter rotiert war. In Brauns Simulationen verlief der Riss konsistent in einem Winkel von etwa 70 Grad, was die experimentellen Beobachtungen perfekt widerspiegelte. Dieses Ergebnis blieb stabil, egal ob das Gitter mit der Platte ausgerichtet oder um 15 oder 30 Grad rotiert war. Auch die Rissgeschwindigkeit blieb stabil und bewegte sich bei etwa 70 Prozent der Schallgeschwindigkeit im Material, unabhängig von der Orientierung des Gitters.

Der zweite Test beinhaltete eine einseitig gekerbte Platte unter Zugspannung, ein Aufbau, der dazu dient, zu untersuchen, wie Risse langsam unter einer stetigen Zugkraft wachsen. Hier war das Ziel zu sehen, ob der Riss in einer geraden Linie von der Kerbe zur gegenüberliegenden Kante verlaufen würde. Traditionelle Gittermodelle zwangen den Riss oft dazu, basierend auf der Ausrichtung des Gitters zu zickzacken oder abzuweichen. Brauns Modell produzierte jedoch für jede Orientierung des Gitters einen geraden Riss. Die Kraft, die erforderlich war, um die Platte zu brechen, war ebenfalls konsistent und variierte um weniger als zwei Prozent, ungeachtet der Rotation des Gitters. Dies demonstrierte, dass das Modell sowohl den Pfad als auch die Festigkeit des Materials genau vorhersagen kann, ohne den künstlichen Einfluss des Computergitters.

Der dritte und komplexeste Test war die Fragmentierung eines dicken Zylinders, der einer internen Explosion ausgesetzt war. In diesem Szenario wird kein einzelner Riss durch eine Kerbe geleitet; stattdessen müssen viele Risse spontan entstehen und nach außen rasen, wodurch der Zylinder in Dutzende von Teilen zersplittert wird. Dies ist ein chaotisches Ereignis, bei dem die Symmetrie des Gitters leicht das Muster der Fragmente diktieren könnte. Als Braun die Simulation mit dem Zentrum des Zylinders als Mittelpunkt des Gitters durchführte, bildeten die Risse ein symmetrisches Muster. Als er das Gitter jedoch leicht verschob, sodass es nicht mehr perfekt zentriert war, verteilten sich die Risse unregelmäßig und die Anzahl der Fragmente änderte sich. Dies bestätigte, dass das Modell empfindlich auf den physischen Aufbau reagierte und nicht bloß auf die Symmetrie des Gitters. Die Simulation sagte zwischen 12 und 19 Fragmente voraus, eine Zahl, die im Bereich anderer fortgeschrittener Methoden lag, aber niedriger war als bei einigen anderen, wahrscheinlich weil die glatte Schadenszone im Modell verhinderte, dass zu nah beieinander liegende Risse entstanden.

Während dieser Tests verfolgten die Forscher sorgfältig die Energie des Systems. Sie verglichen die Energie, die beim Bruch des Materials verloren ging, mit der Energie, die zur Erzeugung der neuen Rissflächen erforderlich war. In einigen Fällen gab es eine kleine Diskrepanz zwischen der durch das Schadensmodell berechneten Energie und der tatsächlich durch das Brechen der Verbindungen freigesetzten Energie. Dies ist eine bekannte Einschränkung des Hybridansatzes, bei dem die Triebkraft des Risses anders berechnet wird als die Energie der gebrochenen Bindungen. Die mechanische Energiebilanz blieb jedoch innerhalb von weniger als einem halben Prozent genau, was bestätigte, dass die Simulation physikalisch fundiert war. Die Ergebnisse legen nahe, dass das Modell zwar nicht perfekt ist, aber einen bedeutenden Schritt nach vorn im Verständnis der dynamischen Fraktur darstellt.

Die Auswirkungen dieser Arbeit reichen weit über bessere Computergrafiken hinaus. Durch die Beseitigung der Abhängigkeit von der Gitterorientierung können Ingenieure darauf vertrauen, dass ihre Simulationen von Crashtests, strukturellen Versagen oder Materialtests die wahre Physik des Materials widerspiegeln. Das Modell überbrückt erfolgreich die Lücke zwischen der diskreten Natur von Atomen und dem kontinuierlichen Verhalten fester Objekte. Es zeigt, dass es mit den richtigen mathematischen Werkzeugen möglich ist, den chaotischen, gewaltsamen Prozess des Brechens zu simulieren, ohne dass das Computergitter das Ergebnis diktiert. Die Studie kommt zu dem Schluss, dass dieser Ansatz eine robuste Methode bietet, um vorherzusagen, wie Materialien unter extremen Bedingungen versagen, und somit einen klareren Weg für die Konstruktion sichererer und widerstandsfähigerer Strukturen ebnet.