Efficient strain-space hyperreduction in large-deformation solid mechanics

Sintesi Tecnica: Iperriduzione Efficiente nello Spazio delle Deformazioni per la Meccanica del Continuo a Grandi Deformazioni

Problematica
Le tecniche di Riduzione dell'Ordine del Modello (MOR) che operano nello spazio delle deformazioni hanno recentemente dimostrato prestazioni superiori nella omogeneizzazione computazionale, in particolare per i Volumi Rappresentativi (RVE). Tuttavia, la loro applicazione è stata finora limitata a questo specifico contesto. L'estensione di queste tecniche alla meccanica del continuo a grandi deformazioni presenta sfide significative. A differenza delle formulazioni basate sullo spostamento, i metodi nello spazio delle deformazioni non possono imporre direttamente le condizioni al contorno (BC) di Dirichlet perché le deformazioni nei punti di integrazione al bordo non sono esplicitamente note, e i valori al contorno sono tipicamente prescritti come spostamenti, non come deformazioni. Inoltre, gli approoli standard nello spazio delle deformazioni spesso si basano su BC periodiche, che non sono applicabili a problemi macroscopici con valori al contorno arbitrari.

Metodologia
Gli autori propongono un framework generalizzato per applicare la MOR nello spazio delle deformazioni a problemi con BC di Dirichlet arbitrarie e parametrizzate. La metodologia centrale prevede tre componenti principali:

1. Gestione del Contorno tramite Lifting: Per soddisfare le BC di Dirichlet non nulle, gli autori impiegano una strategia di lifting. Viene computato offline un campo di riferimento coerente con il contorno per valori di spostamento unitari in ogni bordo indipendente. Il campo di soluzione reale viene quindi decomposto in questo campo di riferimento (scalato dai valori di contorno effettivi) e in un "campo di fluttuazione". Il campo di fluttuazione svanisce a tutti i bordi di Dirichlet per costruzione, permettendo alla base ridotta (derivata dalla Decomposizione Proporzionale Ortogonale, POD) di soddisfare le BC a priori. Questo approccio trasforma le BC in omogenee per il solutore ridotto.

2. Formulazioni nello Spazio delle Deformazioni: Il documento generalizza tre specifiche tecniche di iperriduzione in questo nuovo ambito:

   • ECM (Metodo di Cubatura Empirica) nello Spazio delle Deformazioni: Riformulato per eseguire l'integrazione ridotta direttamente nei punti di Gauss, bypassando le operazioni a livello di elemento.
   • E3C (Cubatura a Cluster Empiricamente Corretta) nello Spazio delle Deformazioni: Generalizzata per problemi meccanici arbitrari. L'E3C ottimizza i punti di integrazione e i pesi nello spazio delle deformazioni, permettendo loro di essere distinti dai punti di Gauss della mesh originale.
   • EMSL (Campionamento del Materiale e Linearizzazione Empirica) nello Spazio delle Deformazioni: Estesa oltre gli RVE. L'EMSL linearizza il comportamento del materiale attorno alle deformazioni medie medie per cluster. A differenza di altri metodi che integrano approssimativamente una legge materiale esatta, l'EMSL integra esattamente un'approssimazione lineare per cluster, risultando in un problema affine per ogni step di carico che non richiede un solutore iterativo di Newton-Raphson durante la fase online.

3. Baseline di Confronto: Questi metodi nello spazio delle deformazioni sono messi a confronto con il metodo standard di Campionamento e Pesatura a Conservazione dell'Energia (ECSW) nello spazio degli spostamenti, che funge da robusta baseline preservante la struttura.

Contributi Chiave

• Generalizzazione della MOR nello Spazio delle Deformazioni: Il lavoro estende con successo le tecniche di MOR nello spazio delle deformazioni (ECM, E3C, EMSL) dall'omogeneizzazione computazionale alla meccanica del continuo a grandi deformazioni con BC di Dirichlet arbitrarie.
• Nuovo Trattamento del Contorno: L'introduzione di una formulazione di campo di fluttuazione basata sul lifting consente di imporre arbitrarie condizioni al contorno nello spazio delle deformazioni senza richiedere dati di deformazione al contorno.
• Benchmarking Completo: Gli autori forniscono un confronto rigoroso tra la riduzione dello spazio degli spostamenti e quella dello spazio delle deformazioni su due esempi iperelastici di stima inversa dei parametri, con parametri materiali variabili e percorsi di carico complessi (trazione, taglio e carico misto).

Risultati
I metodi sono stati testati su due problemi di piastre iperelastiche con fori ellittici:

1. Esempio 1 (Deformazione Moderata): Una piastra sotto carico di trazione con deformazione fino al 20%.
2. Esempio 2 (Grandi Deformazioni e Carico Complesso): La stessa geometria sottoposta a trazione (fino al 50%), taglio (fino al 25%) e carico misto, con una mesh più fine e un costo computazionale più elevato.

Risultati delle Prestazioni:

• Accuratezza vs. Tempo di Esecuzione: In entrambi gli esempi, i metodi nello spazio delle deformazioni hanno superato costantemente il metodo ECSW nello spazio degli spostamenti nel compromesso tra accuratezza e tempo di esecuzione.
• Accelerazioni (Speedups):
  • E3C: Ha ottenuto accelerazioni di circa 10.000 volte (Esempio 1) e 30.000 volte (Esempio 2) rispetto al modello a ordine completo mantenendo un'alta accuratezza (es. <1% di errore). L'E3C ha generalmente fornito l'accuratezza più elevata per un dato budget computazionale tra i metodi nello spazio delle deformazioni quando venivano utilizzati sufficienti punti di integrazione.
  • EMSL: Ha dimostrato un'eccezionale efficienza temporale, raggiungendo accelerazioni fino a 100.000 volte (Esempio 1) e 300.000 volte (Esempio 2). L'EMSL eccelleva quando i budget online e offline erano molto limitati, spesso superando l'E3C in velocità ma talvolta a scapito di una leggera riduzione dell'accuratezza in regimi altamente non lineari rispetto all'E3C con molti punti di integrazione.
  • ECM: Ha performato meglio dell'ECSW ma è rimasto generalmente indietro rispetto a E3C ed EMSL nelle metriche specifiche testate.
• Costi Offline: Sebbene l'E3C abbia una fase di addestramento offline più costosa a causa dell'ottimizzazione iterativa, gli autori osservano che ciò è gestibile e significativamente inferiore al tempo richiesto per generare i dati di addestramento. L'EMSL ha il costo offline più basso.

Significato e Rivendicazioni
Il documento sostiene che la metodologia proposta giustifica ulteriori ricerche sulle tecniche di MOR nello spazio delle deformazioni oltre l'omogeneizzazione computazionale. Gli autori evidenziano che queste tecniche offrono:

• Efficienza Superiore: Drastiche riduzioni del tempo di esecuzione online (ordini di grandezza) mantenendo un'alta accuratezza.
• Non Intrusività: I metodi richiedono solo dati di spostamento, deformazione e stress (facilmente estraibili da solver black-box) e una routine del materiale. Essi bypassano le operazioni a livello di elemento e non richiedono la conoscenza dei dettagli di discretizzazione (es. tipi di elemento) o dei vettori residui.
• Applicabilità: Il trattamento riuscito di arbitrarie BC di Dirichlet rimuove una barriera principale all'applicazione di queste tecniche efficienti a problemi strutturali macroscopici.

Gli autori concludono che, sebbene rimangano sfide per problemi altamente non lineari (come l'elasto-plasticità) e per la stima delle deformazioni per cluster nell'EMSL, le prestazioni dimostrate suggeriscono che questi metodi potrebbero abilitare nuove applicazioni ingegneristiche in ambiti come il controllo, l'ottimizzazione e l'adattamento dei modelli di materiale, dove il costo computazionale è attualmente proibitivo.