Modified augmented Lagrangian preconditioning for mixed-dimensional beam-solid coupling

Sintesi Tecnica: Precondizionamento Augmented Lagrangian Modificato per l'Accoppiamento Beam-Solid a Dimensione Mista

Definizione del Problema
Il documento affronta le sfide numeriche associate all'accoppiamento a dimensione mista di corpi solidi tridimensionali (3D) e travi di Kirchhoff–Love prive di torsione unidimensionali (1D) incorporate. Questo approccio di modellazione è essenziale per simulare materiali rinforzati con fibre, come i compositi a breve fibra o i tessuti biologici, in cui le fibre sono significativamente snelle rispetto alla matrice solida. L'accoppiamento è imposto tramite moltiplicatori di Lagrange per garantire vincoli posizionali esatti tra l'asse centrale della trave e il solido circostante.

Questa formulazione porta a un sistema lineare di equazioni con una struttura a punto di sella indefinita. Il sistema presenta difficoltà specifiche:

1. Indefinitezza: La presenza di un blocco nullo sulla diagonale (associato ai moltiplicatori di Lagrange) impedisce l'uso di solver standard definiti positivi.
2. Sottoproblemi Singolari: Il sottoproblema della trave, quando considerato isolatamente (come richiesto dai precondizionatori a blocchi), è un puro problema di Neumann a causa della mancanza di espliciti vincoli di Dirichlet sulle fibre incorporate. Ciò rende singolare la matrice di rigidezza della trave.
3. Sensibilità ai Parametri: Il condizionamento del sistema è fortemente influenzato dai parametri di modellazione, tra cui il rapporto di rigidezza tra trave e solido ($E = E_B/E_S$), il rapporto di volume tra trave e solido ($V$) e il rapporto di dimensione della mesh.
4. Limitazioni dei Metodi Esistenti: Il lavoro precedente degli autori utilizzava la regolarizzazione per penalità, che riduce la dimensione del sistema ma introduce un cattivo condizionamento e inconsistenza variazionale. Al contrario, gli approcci diretti con moltiplicatori di Lagrange richiedono strategie di precondizionamento specializzate che non sono ancora state pienamente esplorate per questo specifico contesto di elasticità vettoriale a dimensione mista.

Metodologia
Gli autori propongono una strategia di precondizionamento a blocchi Modified Augmented Lagrangian (MAL) per risolvere il sistema a punto di sella risultante. La metodologia procede come segue:

1. Formulazione Aumentata: Il sistema originale con moltiplicatore di Lagrange viene riformulato in un sistema di Augmented Lagrangian. Ciò comporta l'aggiunta di termini simili a penalità ($\epsilon M^T W^{-1} M$ e $\epsilon D^T W^{-1} D$) ai blocchi di rigidezza del solido e della trave. Fondamentalmente, questa aumentazione:

   • Preserva l'applicazione esatta dei vincoli di accoppiamento (a differenza dei metodi a pura penalità).
   • Regolarizza il sottoproblema singolare della trave, rendendo invertibile il blocco di rigidezza della trave anche per piccoli parametri di penalità.
   • Mantiene la struttura a punto di sella dell'intero sistema.

2. Costruzione del Precondizionatore a Blocchi: Partendo da un precondizionatore ideale di tipo Augmented Lagrangian, gli autori derivano varianti triangolari a blocchi pratiche. Questi precondizionatori permettono il trattamento indipendente dei blocchi solido, trave e dello Schur complement:

   • Blocco Solido: Precondizionato utilizzando metodi di Algebraic Multigrid (AMG) (specificamente Smoothed Aggregation e Plain Aggregation).
   • Blocco Trave: Poiché la matrice aumentata della trave è a blocchi diagonali (per fibre indipendenti) e invertibile, viene risolta utilizzando metodi diretti (fattorizzazione LU).
   • Schur Complement: Il cuore del metodo consiste nell'approssimare lo Schur complement aumentato. Gli autori introducono tre distinte varianti:
     • Variante I: Approssima lo Schur complement utilizzando inversi approssimati sparsi (SPAI) per il blocco della trave e un'approssimazione diagonale per il blocco del solido, assumendo un piccolo parametro di penalità ($\epsilon \ll E_B$).
     • Variante II: Utilizza un'approssimazione diagonale basata sulla matrice di scala $W$ (specificamente $\hat{S} \approx -\frac{1}{\epsilon}W$), richiedendo un parametro di penalità accuratamente calibrato per raggruppare gli autovalori.
     • Variante III: Introduce parametri di penalità separati per i contributi del solido ($\epsilon_S$) e della trave ($\epsilon_B$). Impostando $\epsilon_S \to 0$ e calibrando $\epsilon_B$, questa variante mira a disaccoppiare i problemi di scala, utilizzando un'approssimazione diagonale dello Schur complement.

3. Implementazione Numerica: I sistemi lineari vengono risolti utilizzando GMRES. L'implementazione sfrutta i framework software Trilinos e 4C, utilizzando il parallelismo in memoria distribuita (MPI).

Contributi Chiave

• Nuovo Framework di Precondizionamento: Il documento presenta la prima applicazione di precondizionatori a blocchi di tipo Modified Augmented Lagrangian specificamente adattati all'accoppiamento beam-solid a dimensione mista con moltiplicatori di Lagrange.
• Gestione dei Sottoproblemi Singolari: La strategia di aumentazione regolarizza efficacemente il sottoproblema singolare della trave senza ricorrere a tecniche complesse di proiezione dello spazio nullo o metodi di deflazione, consentendo l'uso di solver diretti standard sul blocco della trave.
• Varianti dello Schur Complement: La derivazione sistematica e il confronto di tre strategie di approssimazione dello Schur complement (Variante I, II e III), inclusa l'introduzione di parametri di penalità indipendenti per i domini solido e trave (Variante III).
• Analisi della Scalabilità: Valutazione completa della scalabilità forte e debole su cluster a memoria distribuita con milioni di gradi di libertà.

Risultati
Gli esperimenti numerici su elementi di volume rappresentativi (RVE) e un'applicazione a una piastra composita ibrida hanno prodotto i seguenti risultati:

• Robustezza: I precondizionatori proposti dimostrano robustezza rispetto al rapporto di rigidezza ($E$), al rapporto di volume ($V$) e al raffinamento della mesh. La Variante III, in particolare, ha mostrato una quasi-robustezza rispetto ai parametri per contrasti di rigidezza moderati o elevati.
• Indipendenza dalla Mesh: I metodi hanno raggiunto conteggi di iterazione quasi indipendenti dalla mesh (tipicamente tra 10 e 50 iterazioni per gli RVE, e circa 20 per la piastra composita) attraverso diverse densità di discretizzazione.
• Scalabilità: I metodi hanno mostrato una scalabilità forte e debole favorevole su cluster paralleli (fino a 512 processori). Le Varianti II e III hanno mostrato un'efficienza parallela superiore (fino al 97% di efficienza di scalabilità debole per il setup) rispetto alla Variante I, principalmente a causa del costo computazionale dell'approssimazione dello Schur complement nella Variante I.
• Applicazione Ingegneristica: Nell'applicazione a una piastra sandwich composita ibrida, il precondizionatore ILU ingenuo è fallito nel convergere. Al contrario, i precondizionatori MAL proposti (specialmente la Variante IIIc) hanno ottenuto la convergenza con bassi conteggi di iterazione e tempi di soluzione significativamente ridotti (fino a 4 volte più veloci della variante con le prestazioni peggiori) rispetto agli approcci standard.

Significatività
Il documento sostiene che i precondizionatori a blocchi di tipo Modified Augmented Lagrangian forniscano un'alternativa robusta ed efficiente sia alle formulazioni basate su pura penalità che agli approcci standard con moltiplicatori di Lagrange per simulazioni su larga scala di modelli a dimensione mista in meccanica dei solidi e delle strutture. Consentendo l'applicazione esatta dei vincoli e regolarizzando i sottoproblemi singolari, il metodo facilita la simulazione di strutture rinforzate con fibre complesse con alti contrasti di rigidezza e frazioni di volume variabili. Gli autori concludono che questo approccio è particolarmente adatto per applicazioni ingegneristiche che coinvolgono fibre incorporate, offrendo una via percorribile verso simulazioni ad alta fedeltà e scalabili. Il lavoro futuro prevede l'estensione di questo framework a formulazioni di trave più complesse (ad esempio, Simo–Reissner) che includono effetti di torsione e taglio.