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

Wenn Ingenieure Brücken, Flugzeuge oder sogar alltägliche Werkzeuge entwerfen, müssen sie vorhersehen, wie Materialien versagen werden. Die dramatischste Form des Versagens ist der Bruch, bei dem ein fester Körper plötzlich auseinanderreißt. Jahrzehntelang haben sich Wissenschaftler auf zwei Hauptarten verlassen, um dies zu untersuchen: physische Tests im Labor und Computersimulationen. Physische Tests sind der unbestreitbare Beweis dafür, wie sich ein Material verhält, aber sie sind teuer, zeitaufwendig und zerstörend. Computersimulationen bieten eine schnellere, kostengünstigere Alternative und ermöglichen es Designern, tausende Variationen zu testen, bevor ein einziges Stück Metall geschnitten wird. Eine hartnäckige Herausforderung war jedoch die Lücke zwischen den beiden. Viele Computermodelle sind darauf ausgelegt, Materialien zu beschreiben, die wie Glas brechen, während die reale Welt voller Metalle ist, die sich verbiegen und dehnen, bevor sie brechen. Dieses Missverhältnis zwingt Ingenieure oft dazu, sich zwischen einem Modell zu entscheiden, das mathematisch einfach, aber physikalisch ungenau ist, oder einem, das komplex, aber schwierig auszuführen ist.

Ein Forschungsteam von Universitäten aus Italien und Deutschland hat kürzlich einen überraschenden Mittelweg untersucht. Sie untersuchten, ob ein Computermodell, das speziell für spröde, glasähnliche Materialien entwickelt wurde, das Verhalten eines duktilen, verformbaren Metalls erfolgreich vorhersagen könnte. Ihr Untersuchungsobjekt war eine spezifische Aluminiumlegierung, die häufig im Bauwesen und in der Fertigung verwendet wird und für ihre Fähigkeit bekannt ist, sich vor dem Bruch erheblich zu verformen. Die Forscher wollten sehen, ob sie einen „spröden“ mathematischen Rahmen nutzen könnten, um ein „duktiles“, reales Ereignis zu simulieren, wobei sie sich speziell auf die Entstehung und Ausbreitung von Rissen in einer Metallprobe mit vorgekerbten Einkerbungen konzentrierten. Ihr Ziel war es nicht, jedes winzige Detail der Verformung des Metalls exakt nachzubilden, sondern zu prüfen, ob das Modell zwei kritische Dinge genau vorhersagen kann: die maximale Kraft, die das Metall vor dem Versagen aushalten kann, und den exakten Pfad, den der Riss nimmt, während er durch das Material reißt.

Das Experiment begann mit physischen Tests an metallischen Proben in Knochenform (Dog-Bone-Shape), einer Standardform zur Testung der Materialfestigkeit. Die Forscher verwendeten eine Aluminiumlegierung des Typs EN AW-6060, die Silizium und Magnesium enthält. Sie bereiteten zwei Versionen der Probe vor: eine mit einer Dicke von 4 Millimetern und eine weitere mit 6 Millimetern. Jedes Stück verfügte über zwei kleine, dreieckige Schnitte oder Kerben, die leicht versetzt von der Mitte positioniert waren. Diese Kerben fungieren als Schwachstellen, an denen ein Riss wahrscheinlich beginnen wird. Das Team platzierte diese Proben in einer leistungsstarken Prüfmaschine, die sie auseinanderzog, bis sie brachen. Sie beobachteten, dass die Dicke des Metalls das Ergebnis drastisch veränderte. Die dickere 6-Millimeter-Probe entwickelte einen Riss, der diagonal von einer Kerbe zur anderen verlief und sie in einer einzigen, scharfen Linie verband. Im Gegensatz dazu entwickelte die dünnere 4-Millimeter-Probe zwei separate Risse, die voneinander weg bogen und sich nie trafen. Dieser Unterschied im Verhalten ist ein klassisches Beispiel dafür, wie die physikalischen Einschränkungen der Materialdicke das Versagensverhalten verändern können – ein Phänomen, das in der Ingenieurwissenschaft wohlbekannt ist, aber in einfachen Computermodellen schwer zu erfassen ist.

Um diese Ereignisse zu simulieren, wandten die Forscher eine hochentwickelte Computermethode namens isogeometrischer Analyse an, die fortgeschrittene mathematische Kurven verwendet, um die Form des Objekts genauer darzustellen als traditionelle gitterbasierte Methoden. Sie verwendeten einen spezifischen Typ eines Modells, bekannt als Phasenfeld-Bruchmodell. In diesem Ansatz wird ein Riss nicht als scharfe, gezackte Linie dargestellt, die der Computer manuell verfolgen muss. Stattdessen wird der Riss als eine glatte, verschwommene Zone behandelt, in der das Material allmählich seine Festigkeit verliert. Diese „Phasenfeld“-Variable wirkt wie ein Dimmer, der die Fähigkeit des Materials, zusammenzuhalten, von voll intakt bis vollständig gebrochen herunterregelt. Die Forscher wählten eine Version dieses Modells, die ursprünglich für spröde Materialien entwickelt wurde, was bedeutet, dass sie davon ausgeht, dass das Material ohne nennenswertes Biegen oder Dehnen bricht. Dies war eine kühne Entscheidung, da das von ihnen getestete Aluminium duktil ist, was bedeutet, dass es sich dehnt und Energie durch plastische Verformung absorbiert, bevor es schließlich reißt.

Das Team führte die Simulationen mit denselben Abmessungen und Materialeigenschaften wie die physischen Tests durch. Sie passten einen Schlüsselparameter in ihrem Modell an, die sogenannte interne Länge, welche die Breite der verschwommenen Schadenszone um den Riss herum steuert. Durch die sorgfältige Abstimmung dieses Wertes fanden sie heraus, dass sie die im Labor beobachteten unterschiedlichen Verhaltensweisen replizieren konnten. Wenn sie die interne Länge auf einen größeren Wert setzten, erzeugte die Simulation ein Rissmuster, das die beiden Kerben diagonal verband, was der dicken 6-Millimeter-Probe entsprach. Wenn sie die interne Länge reduzierten, zeigte die Simulation zwei separate, kurvige Risse, die sich nicht trafen, was das Verhalten der dünneren 4-Millimeter-Probe perfekt widerspiegelte. Dies war ein bedeutender Befund, da es darauf hindeutete, dass ein einziger, relativ einfacher mathematischer Rahmen komplexe physikalische Phänomene erfassen kann, indem man lediglich einen Parameter anpasst, ohne die komplizierte Physik des Metallstreckens explizit modellieren zu müssen.

Die Ergebnisse waren nicht nur visuell ähnlich, sondern auch numerisch präzise. Das Computermodell sagte die maximale Kraft, die das Metall vor dem Brechen halten konnte, mit bemerkenswerter Genauigkeit voraus. Für die dickere Probe lag die vorhergesagte Kraft innerhalb von 1 Prozent der tatsächlich gemessenen Kraft. Für die dünnere Probe lag die Vorhersage innerhalb von 3,5 Prozent. Diese Zahlen sind außergewöhnlich nah an den Realwerten für eine Simulation, die die komplexe plastische Verformung, die in echten Metallen auftritt, absichtlich ignoriert hat. Die Forscher merkten an, dass das Modell zwar nicht die langsame, graduelle Erweichung des Metalls nach Erreichen der Spitzenfestigkeit reproduzieren konnte, es aber hervorragend darin war, den Moment des Versagens und den Pfad des Risses vorherzusagen. Dies macht das Modell zu einem leistungsstarken Werkzeug für Ingenieure, die wissen müssen, wo die Sicherheitsgrenzen eines Designs liegen, bevor es gebaut wird.

Die Studie kommt zu dem Schluss, dass die Verwendung eines spröden Bruchmodells zur Simulation duktiler Materialien eine gültige und konservative Strategie für das Ingenieursdesign sein kann. Mit „konservativ“ meinten die Autoren, dass das Modell dazu neigt, das Versagen bei Lasten vorherzusagen, die etwas niedriger liegen als das, was das Material tatsächlich aushalten kann, was ein sicherer Ansatz für sicherheitskritische Anwendungen ist. Der Erfolg dieses Ansatzes beruht auf der spezifischen mathematischen Formulierung, die es ermöglicht, Energie frühzeitig im Prozess abzuleiten, wodurch die Energieabsorption eines duktilen Materials nachgeahmt wird, selbst ohne die Dehnung explizit zu berechnen. Die Forscher räumen ein, dass es sich um einen phänomenologischen Kompromiss handelt; das Modell erzählt nicht die ganze Geschichte dessen, was im Inneren des Metalls nach dem Fließen geschieht, aber es erzählt den wichtigsten Teil: wann und wie es bricht.

Diese Arbeit eröffnet eine neue Perspektive darauf, wie wir Computersimulationen nutzen können, um die physische Welt zu verstehen. Sie legt nahe, dass wir nicht immer das komplexeste Modell benötigen, um das richtige Ergebnis zu erhalten. Manchmal kann ein einfacheres Modell, angewandt mit den richtigen mathematischen Werkzeugen und einem sorgfältigen Verständnis seiner Grenzen, präzise und zuverlässige Vorhersagen liefern. Für die getestete Aluminiumlegierung zeigten die Forscher, dass ein sprödes Modell erfolgreich sowohl die erforderliche Kraft zum Brechen des Teils als auch den komplizierten Pfad des Risses vorhersagen konnte, indem man einfach die Skala der Schadenszone anpasste. Dieser Befund bietet einen praktischen Weg für Ingenieure, die mit Metallen konstruieren, indem er es ermöglicht, effiziente, höherwertige Simulationen zu nutzen, um Sicherheit und Leistung zu gewährleisten, ohne sich in den immensen Rechenkosten der Modellierung jedes mikroskopischen Details der plastischen Verformung zu verlieren. Die Studie steht als Zeugnis für die Kraft der Kombination von strenger Laborprüfung mit fortgeschrittenen numerischen Methoden zur Lösung realer ingenieurtechnischer Herausforderungen.