Efficient strain-space hyperreduction in large-deformation solid mechanics

Technische Zusammenfassung: Effiziente Strain-Space-Hyperreduktion in der großverformungsbasierten Festmechanik

Problemstellung
Modellordnungsreduktionstechniken (Model Order Reduction, MOR), die im Dehnungsraum (strain space) operieren, haben kürzlich eine überlegene Leistung in der computergestützten Homogenisierung gezeigt, insbesondere für Repräsentative Volumenelemente (RVEs). Ihre Anwendung war jedoch weitgehend auf diesen spezifischen Kontext beschränkt. Die Anwendung dieser Techniken auf allgemeine Probleme der großverformungsbasierten Festmechanik stellt erhebliche Herausforderungen dar. Im Gegensatz zu verschiebungsbasierten Formulierungen können Strain-Space-Methoden Dirichlet-Randbedingungen nicht direkt erzwingen, da die Dehnungen an den Integrationspunkten der Grenze nicht explizit bekannt sind und Randwerte typischerweise als Verschiebungen, nicht als Dehnungen, vorgegeben werden. Zudem verlassen sich Standard-Strain-Space-Ansätze oft auf periodische Randbedingungen, die jedoch für makroskopische Probleme mit beliebigen Randwerten nicht anwendbar sind.

Methodik
Die Autoren schlagen ein verallgemeinertes Framework vor, um Strain-Space-MOR auf Probleme mit beliebigen, parametrisierten Dirichlet-Randbedingungen anzuwenden. Die Kernmethodik umfasst drei Hauptkomponenten:

1. Randbehandlung mittels Lifting: Um nicht-verschwindende Dirichlet-Randbedingungen zu erfüllen, verwenden die Autoren eine Lifting-Strategie. Ein randkonsistentes Referenzfeld wird offline für Einheitsverschiebungenswerte an jeder unabhängigen Grenze berechnet. Das tatsächliche Feld wird dann in dieses Referenzfeld (skaliert mit den tatsächlichen Randwerten) und ein „Fluktuationsfeld“ zerlegt. Das Fluktuationsfeld verschwindet konstruktionsbedingt an allen Dirichlet-Rändern, wodurch die reduzierte Basis (abgeleitet aus der Proper Orthogonal Decomposition, POD) die Randbedingungen a priori erfüllt. Dieser Ansatz transformiert inhomogene Randbedingungen in homogene Randbedingungen für den reduzierten Solver.

2. Strain-Space-Formulierungen: Die Arbeit generalisiert drei spezifische Hyperreduktionstechniken für diesen neuen Kontext:

   • Strain-Space ECM (Empirical Cubature Method): Diese wurde so umformuliert, dass sie die reduzierte Integration direkt an den Gauss-Punkten durchführt und dadurch Operationen auf Elementebene umgeht.
   • Strain-Space E3C (Empirically Corrected Cluster Cubature): Generalisiert für beliebige mechanische Probleme. E3C optimiert Integrationspunkte und Gewichte im Dehnungsraum, sodass diese von den ursprünglichen Mesh-Gauss-Punkten abweichen können.
   • Strain-Space EMSL (Empirical Material Sampling and Linearisation): Über RVEs hinaus erweitert. EMSL linearisiert das Materialverhalten um clusterweise erwartete Durchschnittsdehnungen. Im Gegensatz zu anderen Methoden, die näherungsweise eine exakte Materialgesetz integrieren, integriert EMSL exakt eine clusterweise lineare Approximation, was zu einem affinen Problem pro Lastschritt führt, das während der Online-Phase keinen iterativen Newton-Raphson-Solver erfordert.

3. Vergleichsbasis: Diese Strain-Space-Methoden werden gegen die Standard-Displacement-Space-Methode ECSW (Energy Conserving Sampling and Weighting) getestet, die als robuster, struktur-erhaltender Baseline dient.

Zentrale Beiträge

• Generalisierung der Strain-Space-MOR: Die Arbeit überträgt Strain-Space-MOR-Techniken (ECM, E3C, EMSL) erfolgreich von der computergestützten Homogenisierung auf die allgemeine großverformungsbasierte Festmechanik mit beliebigen Dirichlet-Randbedingungen.
• Neuartige Randbehandlung: Die Einführung einer Lifting-basierten Fluktuationsfeld-Formulierung ermöglicht die Implementierung beliebiger Randbedingungen im Dehnungsraum, ohne Dehnungsdaten an den Rändern zu benötigen.
• Umfassendes Benchmarking: Die Autoren bieten einen rigorosen Vergleich zwischen Strain-Space- und Displacement-Space-Hyperreduktion anhand zweier hyperelastischer Beispiele zur inversen Parameterbestimmung, bei variierenden Materialparametern und komplexen Belastungspfaden (Zug, Scherung und gemischte Belastung).

Ergebnisse
Die Methoden wurden an zwei hyperelastischen Plattenproblemen mit elliptischen Löchern getestet:

1. Beispiel 1 (Moderate Deformation): Eine Platte unter Zugbelastung mit bis zu 20 % Dehnung.
2. Beispiel 2 (Große Deformation & Komplexe Belastung): Dieselbe Geometrie, die Zug (bis zu 50 %), Scherung (bis zu 25 %) und gemischte Belastung ausgesetzt ist, mit einem feineren Mesh und höherem Rechenaufwand.

Leistungsbefunde:

• Genauigkeit vs. Laufzeit: In beiden Beispielen übertrafen die Strain-Space-Methoden die Displacement-Space-ECSW-Methode konsistent im Trade-off zwischen Genauigkeit und Laufzeit.
• Beschleunigungen (Speedups):
  • E3C: Erzielte Beschleunigungen von etwa dem 10.000-fachen (Beispiel 1) und dem 30.000-fachen (Beispiel 2) im Vergleich zum Full-Order-Modell bei gleichbleibend hoher Genauigkeit (z. B. < 1 % Fehler). E3C lieferte bei ausreichender Anzahl an Integrationspunkten generell die höchste Genauigkeit für ein gegebenes Rechenbudget unter den Strain-Space-Methoden.
  • EMSL: Demonstrierte eine außergewöhnliche Zeiteffizienz mit Beschleunigungen von bis zu 100.000-fach (Beispiel 1) und 300.000-fach (Beispiel 2). EMSL zeichnete sich aus, wenn die Online- und Offline-Budgets sehr begrenzt waren, und übertraf E3C oft in der Geschwindigkeit, jedoch teilweise auf Kosten einer etwas geringeren Genauigkeit in hochgradig nichtlinearen Regimen im Vergleich zu E3C mit vielen Integrationspunkten.
  • ECM: Performte besser als ECSW, lag jedoch im Vergleich zu E3C und EMSL in den spezifisch getesteten Metriken generell zurück.
• Offline-Kosten: Obwohl E3C aufgrund der iterativen Optimierung eine teurere Offline-Trainingsphase aufweist, merken die Autoren an, dass dies handhabbar ist und signifikant weniger Zeit beansprucht als die Generierung von Trainingsdaten. EMSL weist die niedrigsten Offline-Kosten auf.

Bedeutung und Ansprüche
Das Paper behauptet, dass die vorgeschlagene Methodik die weitere Forschung in Strain-Space-MOR-Techniken über die computergestützte Homogenisierung hinaus rechtfertigt. Die Autoren heben hervor, dass diese Techniken Folgendes bieten:

• Überlegene Effizienz: Drastische Reduktionen der Online-Laufzeit (Größenordnungen) bei gleichzeitiger Beibehaltung hoher Genauigkeit.
• Nicht-Intrusivität: Die Methoden benötigen nur Verschiebungs-, Dehnungs- und Spannungsdaten (leicht aus Black-Box-Solvern extrahierbar) sowie eine Materialroutine. Sie umgehen Operationen auf Elementebene und erfordern keine Kenntnis der Diskretisierungsdetails (z. B. Elementtypen) oder Residuenvektoren.
• Anwendbarkeit: Die erfolgreiche Handhabung beliebiger Dirichlet-Randbedingungen entfernt eine wesentliche Barriere für die Anwendung dieser effizienten Techniken auf makroskopische Strukturprobleme.

Die Autoren kommen zu dem Schluss, dass trotz bestehender Herausforderungen für hochgradig nichtlineare Probleme (wie Elasto-Plastizität) und die Schätzung von clusterweisen Dehnungen in EMSL, die demonstrierte Leistungsfähigkeit darauf hindeutet, dass diese Methoden neue ingenieurtechnische Anwendungen in Bereichen wie Steuerung, Optimierung und Materialmodell-Anpassung ermöglichen könnten, in denen der Rechenaufwand derzeit prohibitiv hoch ist.