Efficient strain-space hyperreduction in large-deformation solid mechanics

Technische Samenvatting: Efficiënte Strain-Space Hyperreductie in Groot-Vervormings Solide Mechanica

Probleemstelling
Model Order Reduction (MOR) technieken die opereren in de rekruimte (strain space) hebben recentelijk superieure prestaties aangetoond in computationele homogenisatie, met name voor Representatieve Volume-elementen (RVE's). De toepassing hiervan is echter grotendeels beperkt gebleven tot deze specifieke context. Het uitbreiden van deze technieken naar algemene problemen in de groot-vervormings solide mechanica brengt aanzienlijke uitdagingen met zich mee. In tegenstelling tot verplaatsingsgebaseerde formuleringen kunnen strain-space methoden Dirichlet-randvoorwaarden niet direct afdwingen, omdat de rekken van integratiepunten op de rand niet expliciet bekend zijn en randwaarden doorgaans worden voorgeschreven als verplaatsingen, niet als rekken. Bovendien vertrouwen standaard strain-space benaderingen vaak op periodieke randvoorwaarden, die niet toepasbaar zijn op macro-schaal problemen met willekeurige randwaarden.

Methodologie
De auteurs stellen een gegeneraliseerd framework voor om strain-space MOR toe te passen op problemen met willekeurige, geparametriseerde Dirichlet-randvoorwaarden. De kernmethodologie bestaat uit drie componenten:

1. Randafhandeling via Lifting: Om niet-nul Dirichlet-randvoorwaarden te voldoen, gebruiken de auteurs een lifting-strategie. Een rand-consistente referentieveld wordt offline berekend voor eenheidswaarden van verplaatsing bij elke onafhankelijke rand. Het werkelijke resultaat wordt vervolgens ontleed in dit referentieveld (geschaald door de werkelijke randwaarden) en een "fluctuatieveld". Het fluctuatieveld verdwijnt bij alle Dirichlet-randen door constructie, waardoor de gereduceerde basis (afgeleid van Proper Orthogonal Decomposition, POD) de randvoorwaarden a priori kan voldoen. Deze aanpak transformeert inhomogene randvoorwaarden naar homogene randvoorwaarden voor de gereduceerde solver.

2. Strain-Space Formuleringen: Het artikel generaliseert drie specifieke hyperreductietechnieken naar deze nieuwe setting:

   • Strain-Space ECM (Empirical Cubature Method): Hergeformuleerd om directe gereduceerde integratie op Gauss-punten uit te voeren, waarbij operaties op elementniveau worden omzeild.
   • Strain-Space E3C (Empirically Corrected Cluster Cubature): Gegeneraliseerd naar willekeurige mechanische problemen. E3C optimaliseert integratiepunten en gewichten in de rekruimte, waardoor deze verschillend kunnen zijn van de oorspronkelijke mesh Gauss-punten.
   • Strain-Space EMSL (Empirical Material Sampling and Linearisation): Uitgebreid voorbij RVE's. EMSL lineariseert het materiaalgedrag rond cluster-gewijze verwachte gemiddelde rekken. In tegenstelling tot andere methoden, die een exacte materiaalwet benaderen, integreert EMSL exact een cluster-gewijze lineaire benadering, wat resulteert in een affien probleem per belastingsstap dat geen iteratieve Newton-Raphson solver vereist tijdens de online fase.

3. Vergelijkingsbaseline: Deze strain-space methoden worden afgezet tegen de standaard verplaatsings-ruimte Energy Conserving Sampling and Weighting (ECSW) methode, die dient als een robuuste, structuur-behoudende baseline.

Belangrijkste Bijdragen

• Generalisatie van Strain-Space MOR: Het werk slaagt erin strain-space MOR technieken (ECM, E3C, EMSL) te generaliseren van computationele homogenisatie naar algemene groot-vervormings solide mechanica met willekeurige Dirichlet-randvoorwaarden.
• Nieuwige Randbehandeling: De introductie van een lifting-gebaseerde fluctuatieveld-formulering maakt het mogelijk om willekeurige randvoorwaarden op te leggen in de rekruimte zonder dat rand-rekdata nodig zijn.
• Uitgebreide Benchmarking: De auteurs bieden een rigoureuze vergelijking tussen strain-space versus verplaatsings-ruimte hyperreductie op twee hyperelastische voorbeelden betreffende inverse parametermestimatie, variërende materiaalparameters en complexe belastingspaden (trek, afschuiving en gemengde belasting).

Resultaten
De methoden werden getest op twee hyperelastische plaatproblemen met elliptische gaten:

1. Voorbeeld 1 (Matige Vervorming): Een plaat onder trekbelasting met tot 20% rek.
2. Voorbeeld 2 (Grote Vervorming & Complexe Belasting): Dezelfde geometrie onderworpen aan trek (tot 50%), afschuiving (tot 25%) en gemengde belasting, met een fijnere mesh en hogere computationele kosten.

Prestatiebevindingen:

• Nauwkeurigheid vs. Runtime: In beide voorbeelden presteerden de strain-space methoden consistent beter dan de verplaatsings-ruimte ECSW in de afweging tussen nauwkeurigheid en runtime.
• Versnellingen:
  • E3C: Bereikte versnellingen van ongeveer 10.000-voudig (Voorbeeld 1) en 30.000-voudig (Voorbeeld 2) ten opzichte van het full-order model terwijl een hoge nauwkeurigheid behouden bleef (bijv. <1% fout). E3C bood over het algemeen de hoogste nauwkeurigheid voor een gegeven computationeel budget onder de strain-space methoden wanneer voldoende integratiepunten werden gebruikt.
  • EMSL: Demonstreerde uitzonderlijke runtime-efficiëntie, met versnellingen tot 100.000-voudig (Voorbeeld 1) en 300.000-voudig (Voorbeeld 2). EMSL blonk uit wanneer de online en offline budgetten zeer beperkt waren, en presteerde vaak beter dan E3C qua snelheid, maar soms ten koste van een iets lagere nauwkeurigheid in zeer niet-lineaire regimes vergeleken met E3C met veel integratiepunten.
  • ECM: Presteerde beter dan ECSW, maar bleef over het algemeen achter bij E3C en EMSL in de specifiek geteste metrieken.
• Offline Kosten: Hoewel E3C een duurdere offline trainingsfase heeft door iteratieve optimalisatie, merken de auteurs op dat dit beheersbaar is en aanzienlijk minder is dan de tijd die nodig is om trainingsdata te generen. EMSL heeft de laagste offline kosten.

Betekenis en Claims
Het artikel claimt dat de voorgestelde methodologie verdere research naar strain-space MOR technieken buiten computationele homogenisatie rechtvaardigt. De auteurs benadrukken dat deze technieken het volgende bieden:

• Superieure Efficiëntie: Dramatische reducties in online runtime (ordes van grootte) terwijl een hoge nauwkeurigheid behouden blijft.
• Niet-Intrusiviteit: De methoden vereisen alleen verplaatsings-, rek- en spanningdata (gemakkelijk te extraheren uit black-box solvers) en een materiaalroutine. Ze omzeilen element-niveau operaties en vereisen geen kennis van de discretisatie details (bijv. elementtypen) of residu-vectoren.
• Toepasbaarheid: Het succesvol afhandelen van willekeurige Dirichlet-randvoorwaarden verwijdert een belangrijke barrière voor het toepassen van deze efficiënte technieken op macro-schaal structurele problemen.

De auteurs concluderen dat hoewel er uitdagingen blijven voor zeer niet-lineaire problemen (zoals elastoplastisiteit) en de schatting van cluster-gewijze rekken in EMSL, de gedemonstreerde prestaties suggereren dat deze methoden nieuwe engineering-toepassingen kunnen mogelijk maken in controle, optimalisatie en materiaalkarakterisering waar de computationele kosten momenteel prohibitief zijn.