💻 computer science

A variational multi-phase model for elastoplastic materials with microstructure evolution

Die Arbeit stellt ein variationsbasiertes Mehrphasenmodell zur Beschreibung der Mikrostrukturentwicklung bei elastoplastischen Materialien mit linearer kinematischer Verfestigung vor, das Dissipationsabstände und Young-Maße zur zeitkontinuierlichen Phasenumwandlung kombiniert und durch einen FEM-Benchmark validiert wird.

Ursprüngliche Autoren: Sarah Dinkelacker-Steinhoff, Klaus Hackl

Veröffentlicht 2026-02-20

📖 4 Min. Lesezeit☕ Kaffeepausen-Lektüre

CC BY 4.0

Ursprüngliche Autoren: Sarah Dinkelacker-Steinhoff, Klaus Hackl

Originalarbeit lizenziert unter CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Dies ist eine KI-generierte Erklärung des untenstehenden Papers. Sie wurde nicht von den Autoren verfasst oder gebilligt. Für technische Genauigkeit konsultieren Sie das Originalpaper. Vollständigen Haftungsausschluss lesen

Das große Ganze: Ein Material, das sich verwandeln kann

Stellen Sie sich vor, Sie halten einen Klotz aus einem besonderen Material in der Hand. Wenn Sie ihn drücken oder ziehen, passiert etwas Magisches: Er ist nicht starr wie ein Stein, sondern verändert sich innerlich. Es ist, als würde er aus verschiedenen "Personen" bestehen, die je nach Situation die Führung übernehmen.

Die Autoren Sarah Dinkelacker-Steinhoff und Klaus Hackl haben ein mathematisches Modell entwickelt, um genau dieses Verhalten vorherzusagen. Es geht um Materialien, die elastisch (wie ein Gummiband, das zurückfedert) und plastisch (wie Knete, die ihre Form behält) sind, aber zusätzlich eine innere Struktur haben, die sich während der Belastung verändert.

Die Hauptakteure: Drei Metaphern

Um das komplexe Mathematik-Modell zu verstehen, nutzen wir drei einfache Bilder:

1. Das Team aus drei Spezialisten (Die Phasen)

Stellen Sie sich das Material nicht als einen homogenen Brocken vor, sondern als ein Team aus drei verschiedenen Spezialisten (wir nennen sie Phase 1, 2 und 3).

Jeder Spezialist hat seine eigene Stärke (Steifigkeit) und seine eigene Schwäche (wann er nachgibt).
Zu Beginn ist vielleicht nur Spezialist 3 im Raum (100 % des Materials).
Wenn Sie das Material stark belasten (z. B. drücken), merkt Spezialist 3: "Das ist zu viel für mich!"
Dann ruft er die anderen beiden Spezialisten an. Plötzlich arbeiten alle drei zusammen. Das Material besteht nun aus einer Mischung aus allen drei.

2. Der Übergangs-Manager (Die Wahrscheinlichkeitsmatrix)

Wie weiß das Material, wann es die Spezialisten wechseln muss? Hier kommt die "Wahrscheinlichkeitsmatrix" ins Spiel.

Stellen Sie sich das wie einen Türsteher oder einen Manager vor, der entscheidet: "Wer darf jetzt arbeiten?"
Wenn der Druck zu hoch wird, öffnet der Manager die Tür für den nächsten Spezialisten.
Das Besondere an diesem Modell ist, dass dieser Wechsel nicht langsam passiert, sondern sofort (wie ein Blitz), sobald eine bestimmte Schwelle erreicht ist. Das Modell berechnet genau, wie schnell und in welchem Verhältnis diese Wechsel passieren.

3. Der "Young-Maß"-Schwarm (Die feine Verteilung)

Das ist der kniffligste Teil, den die Autoren mit einem cleveren Trick lösen.

Wenn sich das Material verändert, bilden sich winzige Bereiche, in denen nur Spezialist 1 arbeitet, und daneben Bereiche mit nur Spezialist 2. Diese Bereiche sind so klein, dass man sie mit bloßem Auge nicht sieht.
Anstatt jedes einzelne winzige Teilchen zu berechnen (was unmöglich wäre), nutzen die Autoren ein Konzept namens Young-Maß.
Die Analogie: Stellen Sie sich einen großen Schwarm Bienen vor. Sie können nicht jede einzelne Biene verfolgen. Aber Sie können sagen: "In diesem Bereich sind 40 % rote Bienen und 60 % blaue Bienen." Das Modell berechnet also nicht die einzelnen Bienen, sondern nur den Anteil der verschiedenen Farben in einem bestimmten Bereich. So wird aus einem chaotischen Durcheinander eine glatte, berechenbare Entwicklung.

Was passiert im Computer? (Die Simulation)

Die Autoren haben dieses Modell in einen Computer programmiert und ein Test durchgeführt:

Der Test: Sie haben eine quadratische Platte mit einem Loch in der Mitte simuliert und sie von der Seite zusammengedrückt.
Das Ergebnis: Das Material verhält sich nicht wie ein einfacher Metallblock.
- Zuerst federt es zurück (elastisch).
- Dann fängt es an, sich zu verformen (plastisch).
- Gleichzeitig beginnt das "Team" zu wechseln. Neue Spezialisten (Phasen) tauchen auf, besonders dort, wo der Druck am höchsten ist (am Rand des Lochs).
- Am Ende hat sich das Material so weit verändert, dass es in einem neuen Gleichgewicht ist. Es wird nicht mehr weiter umgebaut, sondern verhält sich wieder stabil (dies nennen sie "Shakedown" – wie ein Schiff, das nach einer Welle ruhig weitersegelt).

Warum ist das wichtig?

Früher haben Ingenieure oft angenommen, dass Materialien einfach nur "nachgeben", wenn sie zu stark belastet werden. Dieses Modell zeigt aber, dass Materialien intelligent auf Belastung reagieren, indem sie ihre innere Struktur umorganisieren.

Das ist besonders wichtig für:

Sicherheitsanalysen: Um zu verstehen, warum Brücken oder Flugzeugteile reißen.
Neue Materialien: Um Materialien zu entwickeln, die sich selbst reparieren oder extremen Drücken standhalten (wie Glas, das unter Druck härter wird, oder Metalle, die sich verformen, ohne zu brechen).

Zusammenfassung in einem Satz

Die Autoren haben eine mathematische "Rezeptur" entwickelt, die beschreibt, wie ein Material unter Druck nicht einfach nur verbiegt, sondern wie es sich in eine Mischung aus verschiedenen inneren Zuständen verwandelt, um die Belastung zu überstehen – alles berechnet mit Hilfe von Wahrscheinlichkeiten und cleveren Durchschnittswerten.

Technische Zusammenfassung: Ein variationsbasiertes Mehrphasen-Modell für elastoplastische Materialien mit Mikrostrukturevolution

Autoren: Sarah Dinkelacker-Steinhoff und Klaus Hackl (Ruhr-Universität Bochum)

1. Problemstellung und Motivation

Das Ziel der Arbeit ist die Entwicklung eines allgemeinen Modells zur Beschreibung der Mikrostrukturentwicklung in elastoplastischen Materialien, die Phasenumwandlungen unterliegen. Solche Phänomene treten in verschiedenen Materialien und Skalen auf, z. B. bei:

Scherbändern in granularen Stoffen (Dilatanz).
Phasentrennungen und Umwandlungen in amorphen Metallen (Bulk Metallic Glasses) und glasartiger Silica unter hohem Druck.
Änderungen der Koordinationszahl und der Ringverteilung in amorphen Strukturen.

Herausfordernd ist die Modellierung des gekoppelten Verhaltens von plastischer Verformung und der Entstehung neuer Phasen, wobei die zugrundeliegenden Kräfte oft unterschiedlich, das makroskopische plastische Verhalten jedoch ähnlich ist (z. B. durch Mohr-Coulomb- oder Drucker-Prager-ähnliche Modelle beschreibbar). Bisherige Ansätze betrachten oft nur einzelne Aspekte; dieses Werk verbindet stochastische Wahrscheinlichkeitsmatrizen mit variationsbasierten Prinzipien, um eine zeitkontinuierliche Evolution der Mikrostruktur zu erfassen.

2. Methodik und theoretischer Rahmen

Das Modell basiert auf der Variationsmechanik für inelastische Materialien unter Verwendung der infinitesimalen Dehnungstheorie.

Thermodynamisches Fundament:
Das System wird durch zwei Potentiale beschrieben: die Helmholtz-Freie-Energie-Dichte $\Psi$ und das Dissipationspotential $\Delta$ . Die Gleichgewichtsbedingungen leiten sich aus dem Prinzip des Minimums des Dissipationspotentials ab (basierend auf dem zweiten Hauptsatz der Thermodynamik).
Einphasen-Beschreibung:
Für einen einzelnen Phasenzustand wird ein Material mit linearer kinematischer Verfestigung angenommen. Die Energie setzt sich aus einem elastischen Term und einem Verfestigungsterm zusammen. Die Dissipation wird durch eine Reibungsgrenze (Dry Friction) modelliert, was zu einer ratenunabhängigen Plastizität führt.
Verallgemeinerte Mehrphasen-Ansätze:
Für Materialien, die während der Verformung mehrere Phasen ausbilden können ( $k$ Phasen), wird ein stochastischer Ansatz gewählt:
- Volumenanteile ( $\lambda_i$ ): Beschreiben den Anteil der einzelnen Phasen ( $\sum \lambda_i = 1$ ).
- Übergangsmatrix ( $g_{ij}$ ): Eine Wahrscheinlichkeitsmatrix steuert die kinetischen Übergänge von Phase $i$ zu Phase $j$ . Die Änderungsrate der Volumenanteile wird durch $\dot{\lambda}_i = \sum (g_{ji} - g_{ij})$ bestimmt.
- Dissipationspotentiale: Es werden drei Arten von Dissipation definiert:
  1. $\Delta_{plast}$ : Plastische Verformung innerhalb der Phasen.
  2. $\Delta_{trans}$ : Dissipationsdistanz für den Phasenübergang (sofortiger Übergang modelliert durch eine Distanzfunktion $D$ ).
  3. $\Delta_{reg}$ : Eine viskose Regularisierung, um die mathematische Behandlung zu vereinfachen und die numerische Robustheit zu erhöhen.
Relaxierte Energie und Young-Maße:
Um eine kontinuierliche Evolution der Mikrostruktur zu beschreiben, werden generalisierte Young-Maße verwendet. Die effektive (relaxierte) freie Energie $\Psi_{rel}$ wird durch Minimierung der Gesamtenergie unter Berücksichtigung der Volumenmittelwerte der Dehnungen bestimmt. Dies führt zu effektiven Tensoren (z. B. effektive Steifigkeit $C_{eff}$ nach Reuß), die eine homogenisierte Betrachtung auf makroskopischer Ebene ermöglichen.
Evolutionäre Gleichungen:
Durch Minimierung des Lagrange-Funktionals bezüglich der kinetischen Variablen ( $\dot{\varepsilon}^p_i$ und $g_{ij}$ ) werden die Evolutionsgleichungen hergeleitet. Diese führen zu Fließregeln für die plastische Dehnung und Bedingungen für den Phasenübergang, die durch Yield-Funktionen ( $\phi_i$ für Plastizität, $\phi_{ij}$ für Phasenwechsel) charakterisiert sind. Neue Phasen werden initialisiert, sobald die thermodynamische Bedingung für einen Übergang erfüllt ist.

3. Numerische Ergebnisse

Die Validierung des Modells erfolgte durch Finite-Elemente-Methoden (FEM) in zwei Dimensionen unter Verwendung der Software Gmsh und Julia (Ferrite).

Materialpunkt-Verhalten (1D):
Ein zyklischer Zugtest mit drei Phasen wurde simuliert.
- Sobald die Streckgrenze erreicht ist, entstehen sofort neue Phasen ( $\lambda_1, \lambda_2$ ), während die ursprüngliche Phase ( $\lambda_3$ ) abnimmt.
- Das Modell zeigt ein "Shakedown"-Verhalten: Nach einer anfänglichen Phasenumwandlung stabilisiert sich das System auf einem konstanten Volumenanteil, und das Material verhält sich danach rein elastisch.
2D-Benchmark-Test:
Eine quadratische Platte mit einem zentralen Loch wurde unter Druckbelastung in x-Richtung analysiert (Symmetrie: Viertelmodell).
- Phasenentwicklung: Die Phasenumwandlung beginnt an der Lochrandgrenze (Spannungskonzentration) und breitet sich durch das gesamte Bauteil aus.
- Konvergenz: Es wurde eine Netzverfeinerungsstudie durchgeführt (grob, fein, superfein). Die Ergebnisse zeigen eine gute Konvergenz der Spannungen, der von-Mises-Spannungen und der Phasenvolumenanteile über die Zeit.
- Ergebnis: Nach 700 Zeitschritten stabilisiert sich der Volumenanteil der ursprünglichen Phase $\lambda_3$ bei ca. 0,54, was bedeutet, dass etwa die Hälfte des Volumens umgewandelt wurde.

4. Wichtige Beiträge

Kopplung von Stochastik und Variationsprinzipien: Die Arbeit verbindet erstmals die Beschreibung von Phasenübergängen über eine Wahrscheinlichkeitsmatrix mit dem Konzept der Dissipationsdistanz und Young-Maßen in einem variationsbasierten Rahmen.
Zeitkontinuierliche Mikrostrukturevolution: Durch die Verwendung von Young-Maßen und der Relaxation der Energie wird eine glatte, zeitkontinuierliche Entwicklung der Mikrostruktur ermöglicht, anstatt diskreter Sprünge.
Robuste Numerik: Die Einführung einer viskosen Regularisierung ( $\Delta_{reg}$ ) ermöglicht eine stabile numerische Implementierung und Konvergenz in FEM-Simulationen.
Validierung: Das Modell wurde erfolgreich an einem Materialpunkt und einem komplexen 2D-Benchmark-Problem (Platte mit Loch) validiert, wobei die physikalisch erwarteten Phänomene (lokale Umwandlung, Shakedown) korrekt abgebildet wurden.

5. Bedeutung und Ausblick

Das vorgestellte Modell bietet einen allgemeinen Rahmen, um die Kopplung von plastischer Verformung und Phasenumwandlung in komplexen Materialien (wie amorphen Metallen oder granularen Stoffen) zu verstehen. Es zeigt, dass die Phasenentwicklung über eine zugrundeliegende Übergangsmatrix steuerbar ist.

Zukünftige Arbeiten sollen das Modell erweitern, indem:

Die Übergangsraten als exponentielle Funktion von Materialparametern modelliert werden, um interne Wechselwirkungen besser abzubilden (ähnlich wie in [16, 17]).
Nicht-isotherme Situationen betrachtet werden.
Komplexere Mikrostrukturen (z. B. Lamellenstrukturen) direkt in das Modell integriert werden.

Zusammenfassend stellt die Arbeit einen signifikanten Schritt dar, um die mikroskopische Evolution von Materialien unter Last durch ein mathematisch konsistentes und numerisch handhabbares variationsbasiertes Modell zu beschreiben.