A blended approach for evolving phase fields using peridynamics: Cyclic loading in quasi-brittle fracture

Die vorgestellte Arbeit führt eine neue Feldtheorie ein, die Peridynamik mit einem zweipunktigen, history-abhängigen Phasenfeld kombiniert, um irreversible Verformungen und schädigungsbedingte Elastizitätsmoduln in quasibrittle Materialien unter zyklischer Belastung zu modellieren und dabei thermodynamische Konsistenz sowie experimentelle Übereinstimmung bei Rissbildung und Energiedissipation zu gewährleisten.

Das Wesentliche

Der unsichtbare Kleber, der vergisst: Ein neuer Weg, Risse vorherzusagen

Stellen Sie sich vor, Sie halten einen alten, spröden Betonblock in der Hand. Wenn Sie ihn belasten, passiert etwas Interessantes: Er fängt nicht nur an zu reißen, sondern er „vergisst" auch, wie er sich vorher verhalten hat. Er wird weicher, verformt sich dauerhaft und speichert die Geschichte jeder Belastung in sich.

Die Forscher Hayden Bromley und Robert Lipton haben einen neuen, cleveren Weg entwickelt, um genau dieses Verhalten – das Brechen von Materialien wie Beton oder Stein unter wiederholtem Stress – am Computer zu simulieren. Ihr Ansatz ist wie eine Mischung aus zwei Welten: der klassischen Physik (wie Dinge sich bewegen) und einer Art „Gedächtnis-Technologie" für das Material.

Hier ist die Idee, zerlegt in einfache Bilder:

1. Das Problem: Warum alte Methoden scheitern

Stellen Sie sich vor, Sie wollen vorhersagen, wie ein Eiswürfel bricht.

• Der alte Weg: Man benutzt zwei getrennte Regler. Einer sagt dem Eis, wie es sich bewegt (Newton). Der andere sagt ihm, wann es bricht (ein separater „Bruch-Algorithmus"). Das ist wie ein Orchester, bei dem der Dirigent und die Geiger nicht miteinander sprechen. Oft führt das zu ungenauen Ergebnissen, besonders wenn das Material nicht nur bricht, sondern sich auch plastisch verformt (wie Knete, die sich nicht zurückfedert).
• Der neue Weg (die „Blended"-Methode): Die Forscher sagen: „Lass uns das alles in einen Regler packen." Sie nutzen Newtons zweites Gesetz (Kraft = Masse × Beschleunigung) als einzigen Motor. Aber sie geben diesem Motor ein Gedächtnis.

2. Die Analogie: Der unsichtbare Kleber (Peridynamik)

Stellen Sie sich das Material nicht als festen Stein vor, sondern als eine Wolke aus winzigen Punkten. Jeder Punkt ist mit seinen Nachbarn durch unsichtbare, elastische Gummibänder verbunden.

• Peridynamik: Das ist die Idee, dass Punkte nicht nur mit ihren direkten Nachbarn sprechen, sondern mit einer ganzen Gruppe um sich herum (einem „Horizont").
• Das Brechen: Wenn Sie das Material zu stark dehnen, reißen diese Gummibänder. Sobald ein Band reißt, ist die Verbindung weg. Das ist der Riss.

3. Das Herzstück: Das „Gedächtnis" der Gummibänder (Phasenfeld)

Das Besondere an dieser neuen Methode ist, dass die Gummibänder ein Gedächtnis haben.

• Zwei Arten von Gedächtnis: Ein Gummiband kann sich unter Zug (Ziehen) anders verhalten als unter Druck (Stoßen).
  • Beim Ziehen: Das Band wird schwächer, dehnt sich und vergisst, wie steif es war. Es speichert, wie sehr es schon gedehnt wurde.
  • Beim Drücken: Es bleibt oft hart und elastisch.
• Der „Unloading-Ratio" (Der Entlastungs-Faktor): Stellen Sie sich vor, Sie ziehen an einem Gummiband und lassen dann los. Es kehrt nicht ganz in die Ursprungsform zurück. Ein Teil der Dehnung bleibt als „Plastizität" (bleibende Verformung) zurück. Die Forscher haben einen Parameter eingeführt (genannt $\beta$), der misst: Wie viel von der Verformung ist dauerhaft, und wie viel ist nur elastisch? Dieser Wert wird aus echten Experimenten mit Beton entnommen.

4. Wie die Simulation funktioniert (Die Geschichte des Materials)

Statt zwei separate Gleichungen zu lösen, nutzen die Forscher nur eine: Newtons Bewegungsgesetz.

• Das Material bewegt sich.
• Während es sich bewegt, „merkt" sich jedes Gummiband, wie stark es je belastet wurde.
• Wenn die Belastung zu hoch wird, wird das Band weich (der „Phasenwert" sinkt von 1 auf 0).
• Wenn es ganz auf 0 sinkt, ist das Band gerissen.
• Das Wunder: Aus dieser einzigen Bewegungsgleichung entstehen automatisch die Risse, die Energieverluste (Dämpfung) und die bleibenden Verformungen. Man muss nichts extra „vorgeben". Die Physik macht das alles von selbst.

5. Warum ist das wichtig? (Die praktischen Ergebnisse)

Die Forscher haben ihre Methode an echten Beton-Proben getestet:

• Zyklische Belastung: Sie haben den Beton hin und her gebogen (wie bei einem Bruchtest, bei dem man immer wieder nachlässt und wieder belastet). Das Modell konnte genau die „Schleifen" (Hysterese) vorhersagen, die man im Labor sieht. Es wusste, wann das Material müde wurde.
• Größeneffekt: Ein kleines Betonstück bricht anders als ein riesiger Betonbalken. Klassische Modelle scheitern oft daran, diesen Unterschied vorherzusagen. Das neue Modell hat das jedoch perfekt getroffen, ohne dass man extra Regeln für die Größe eingeben musste. Es entstand natürlich aus der Simulation heraus.
• Rissmuster: Die Computer-Simulationen sahen fast genauso aus wie die echten Risse in den Experimenten.

Zusammenfassung in einem Satz

Die Forscher haben eine Methode entwickelt, bei der ein Computer-Modell ein Material wie einen lebenden Organismus behandelt: Es bewegt sich nach den Gesetzen der Physik, hat aber ein Gedächtnis für seine Belastungsgeschichte, wodurch es Risse, Verformungen und das Brechen von sich aus genau so vorhersagt, wie es in der realen Welt passiert – ganz ohne komplizierte Zusatzregeln.

Warum das cool ist: Es ist wie ein Simulator für Materialschicksale. Man kann damit testen, wie Brücken, Gebäude oder Flugzeugteile unter jahrelangem Stress altern und wo genau sie brechen werden, bevor es passiert. Und das Beste: Es funktioniert ohne ein festes Gitter (Mesh), was bedeutet, dass Risse überall entstehen können, ohne dass der Computer „stolpert".

Technische Zusammenfassung

1. Problemstellung

Das zentrale Problem der Kontinuumsmechanik besteht darin, eine geeignete Feldtheorie zu finden, die die Bruchentwicklung (Rissbildung) in quasi-spröden Materialien (wie Beton, Mauerwerk, Geomaterialien) unter zyklischer Belastung präzise beschreibt.

• Herausforderung: Herkömmliche Ansätze trennen oft die Verschiebungsgleichung (basierend auf Newtons zweitem Gesetz) von der Phasenfeld-Evolution (oft basierend auf Ginzburg-Landau- oder Variationsprinzipien). Dies führt zu komplexen Kopplungsproblemen und Schwierigkeiten bei der Modellierung von irreversiblen Deformationen (Plastizität) und Hysterese-Effekten.
• Ziel: Entwicklung einer einheitlichen Theorie, die elastische Verformung, irreversible (plastische) Deformation und Schädigung (Damage) in einem einzigen, wohlgestellten Anfangsrandwertproblem (IBVP) vereint, ohne separate Evolutionsgleichungen für das Phasenfeld zu benötigen.

2. Methodik

Die Autoren schlagen einen hybriden (blended) Ansatz vor, der Peridynamik (eine nicht-lokale Theorie) mit einem phasenfeldbasierten Schädigungsmodell kombiniert.

• Nicht-lokale Formulierung: Die Dynamik wird ausschließlich durch Newtons zweites Gesetz gesteuert. Die Verschiebung $u(t, x)$ entwickelt sich autonom. Es gibt keine separate Differentialgleichung für das Phasenfeld; dieses ergibt sich aus dem Verschiebungsfeld.
• Zwei-Punkt-Phasenfeld und Gedächtnis:
  • Das Materialverhalten wird durch ein zwei-Punkt-Phasenfeld $\gamma(t, y, x)$ beschrieben, das die Steifigkeit der „Bindung" (Bond) zwischen zwei Punkten $x$ und $y$ charakterisiert.
  • Es werden zwei unabhängige Phasenfelder eingeführt: $\gamma_+$ für Zugspannungen und $\gamma_-$ für Druckspannungen.
  • Ein Gedächtniseffekt wird durch die maximale historische Dehnung $S^*$ implementiert.
• Schädigungs- und Plastizitätsmodell:
  • Die Dehnung wird additiv in elastische und irreversible (plastische) Anteile zerlegt.
  • Ein Entlastungsverhältnis (Unloading Ratio) $\beta$ wird eingeführt. Dieses Parameter steuert das Verhältnis zwischen plastischer Dehnung und der Steifigkeitsdegradation (elastische Weichmachung).
  • Der Wertebereich $0 \le \beta \le 1$ erlaubt die Modellierung von rein sprödem Verhalten ($\beta=0$) bis hin zu quasi-sprödem Verhalten mit ausgeprägter Plastizität und Hysterese ($0 < \beta \le 1$).
• Konstitutives Gesetz: Die Kraft zwischen zwei Punkten hängt von der Dehnung und dem Phasenfeld ab. Bei Zugversagen wird die Steifigkeit auf Null reduziert (Riss), während das Material unter Druck weiterhin elastisch widersteht (Kontaktbedingung über Risse).
• Mathematische Fundierung: Das Problem wird als Anfangsrandwertproblem in einem Banach-Raum formuliert. Die Existenz und Eindeutigkeit der Lösung wird mittels Fixpunktsätzen nachgewiesen.

3. Wichtige Beiträge und theoretische Ergebnisse

• Thermodynamische Konsistenz: Die Theorie erfüllt automatisch den Energieerhaltungssatz. Die Dissipationsrate (Energieverlust durch Schädigung und Plastizität) ist stets nicht-negativ, was den Gesetzen der Thermodynamik entspricht. Dies ergibt sich direkt aus den Bewegungsgleichungen, ohne zusätzliche Postulate.
• Rückgewinnung der Bruchenergie: Die Energie, die zur Bildung eines flachen Risses benötigt wird, lässt sich direkt aus dem Modell ableiten. Sie entspricht dem Produkt aus Energiefreisetzungsrate ($G_c$) und Rissfläche. Dies ist unabhängig von der nicht-lokalen Längenskala $\epsilon$.
• Mesh-Free Methode: Da Peridynamik keine Ableitungen, sondern Differenzenquotienten verwendet, ist die numerische Implementierung mesh-frei (keine Elemente oder geometrische Verbindungen zwischen Knoten nötig).
• Kalibrierung: Die Modellparameter (Elastizitätsmodul, Festigkeit, Bruchzähigkeit) werden direkt mit experimentellen Materialdaten verknüpft. Wichtig ist, dass die nicht-lokale Längenskala $\epsilon$ kein Materialparameter ist, sondern ein numerischer Auflösungsparameter, der kleiner als die charakteristische Prozesszone $L$ gewählt wird.

4. Numerische Ergebnisse und Validierung

Die Autoren validierten das Modell an mehreren Benchmark-Problemen und verglichen die Ergebnisse mit experimentellen Daten:

1. Zyklische 3-Punkt-Biegeversuche an Beton:
   • Das Modell konnte die Hysterese (Last-Öffnungs-Diagramm) bei zyklischer Belastung quantitativ und qualitativ genau abbilden.
   • Die plastische Dehnungskonzentration und die Rissbildung wurden erfolgreich lokalisiert.
2. Gemischte Modus-Brüche (Mixed-Mode):
   • Simulationen an Betonbalken mit versetzten Rissen zeigten eine hohe Übereinstimmung mit experimentellen Risspfaden und Kraft-Verformungs-Kurven.
   • Auch bei L-förmigen Proben ohne Vorriß (Winkler-Test) wurde der korrekte Risspfad (Übergang von Modus I zu gemischtem Modus) vorhergesagt.
3. Größeneffekt (Size Effect):
   • Das Modell erfasst den deterministischen Größeneffekt bei Betonbalken unterschiedlicher Größe. Die Spannungen bei Versagen nehmen mit zunehmender Bauteilgröße ab, was experimentell bekannt ist, aber in klassischen lokalen Theorien oft nicht korrekt abgebildet wird.
   • Die Ergebnisse stimmen mit Daten von García-Álvarez et al. (2012) und anderen Peridynamik-Modellen überein.

5. Bedeutung und Ausblick

• Einheitlicher Rahmen: Der vorgestellte Ansatz bietet einen einheitlichen Rahmen für quasi-spröde Materialien, der sowohl monotone als auch zyklische Belastungen, von statisch bis dynamisch, abdeckt.
• Effizienz und Genauigkeit: Durch die Eliminierung separater Phasenfeld-Gleichungen und die Nutzung einer mesh-freien Diskretisierung wird die Rechenkomplexität reduziert, während die physikalische Genauigkeit erhalten bleibt.
• Vorhersagekraft: Das Modell kann Rissmuster, Schädigungszonen und plastische Hysterese vorhersagen, ohne dass spezifische Regeln für die Rissausbreitung (wie in klassischen Bruchmechanik-Ansätzen) manuell definiert werden müssen.
• Zukunft: Die Arbeit legt den Grundstein für die Erweiterung auf anisotrope Medien und Verbundwerkstoffe sowie für die Nutzung tensorwertiger Phasenfelder.

Zusammenfassend stellt diese Arbeit einen bedeutenden Fortschritt in der numerischen Bruchmechanik dar, indem sie die Vorteile der Peridynamik (natürliche Rissbildung) mit der physikalischen Genauigkeit von Schädigungs-Plastizitäts-Modellen für zyklische Belastungen verbindet.