A robust mixed finite element formulation for third medium contact

Sintesi Tecnica: Una Robusta Formulazione di Elementi Finiti Misti per il Contatto di Terzo Mezzo

Problematica
Le formulazioni classiche del contatto a elementi finiti si basano sulla discretizzazione esplicita delle interfacce di contatto, imponendo vincoli tramite metodi di penalità, moltiplicatori di Lagrange o barriere. Sebbene maturi, questi approcci affrontano carichi algoritmici significativi in scenari di grande scorrimento, auto-contatto, severa distorsione della mesh o ottimizzazione topologica dove i confini di contatto sono ignoti a priori. Il metodo del "contatto di terzo mezzo" (TMC) affronta questo problema sostituendo i vincoli espliciti con un mezzo fittizio altamente conforme che riempie lo spazio tra i corpi potenzialmente in contatto. Questo mezzo rimane morbido prima del contatto, ma si irrigidisce drasticamente sotto compressione, trasmettendo le forze di contatto attraverso il continuo.

Tuttavia, una sfida centrale nel TMC è la stabilizzazione. Il mezzo deve essere abbastanza morbido da evitare forze spurie pre-contatto, ma abbastanza stabile da resistere a estreme compressioni e tagli senza distorsione della mesh. Le strategie di stabilizzazione esistenti richiedono spesso elementi finiti di ordine superiore (per calcolare le derivate seconde dello spostamento), introducono campi ausiliari globali che aumentano la dimensione del sistema, o si affidano a discretizzazioni specializzate. Esiste la necessità di un meccanismo di stabilizzazione robusto compatibile con elementi finiti di primo ordine semplici e computazionalmente efficienti, che eviti la valutazione diretta delle derivate seconde dello spostamento.

Metodologia
Gli autori propongono una formulazione di elementi finiti misti che introduce una stabilizzazione tramite campo ausiliario per il terzo mezzo. Il nucleo del metodo prevede l'introduzione di:

1. Introduzione di un Campo Ausiliario: Un campo simile al gradiente di deformazione, denotato come $\Theta$, viene introdotto nel terzo mezzo e trattato come un incognita indipendente, disaccoppiata dal campo di spostamento fisico $u$.
2. Accoppiamento Debole: Il campo ausiliario $\Theta$ è debolmente accoppiato al gradiente di deformazione fisico $F$ tramite un termine di penalità ($W^p_{tm} \propto |\Theta - F|^2$). Ciò assicura che $\Theta$ segua la deformazione fisica.
3. Regolarizzazione del Gradiente: Inveve di penalizzare il gradiente del gradiente di deformazione fisico ($\nabla F$), che richiederebbe derivate seconde dello spostamento, la formulazione penalizza il gradiente del campo ausiliario ($\nabla \Theta$). L'energia di regolarizzazione è definita come $W^r_{tm} \propto |\nabla \Theta|^2$.
4. Discretizzazione Mista: La formulazione consente l'interpolazione indipendente di $u$ e $\Theta$. Gli autori indagano sia le approssimazioni continue che quelle discontinue per elemento per $\Theta$.
   • $\Theta$ Continuo: I valori nodali sono condivisi tra gli elementi, creando un accoppiamento inter-elemento.
   • $\Theta$ Discontinuo: I gradi di libertà sono locali a ciascun elemento. Questi possono essere eliminati tramite condensazione statica, garantendo che il sistema globale mantenga solo i gradi di libertà dello spostamento.

Contributi Chiave

• Compatibilità con Ordini Bassi: La formulazione raggiunge una regolarizzazione di tipo gradiente senza calcolare le derivate seconde del campo di spostamento, rendendola compatibile con elementi finiti standard di primo ordine (lineari) e secondo ordine.
• Meccanismo di Stabilizzazione: Il documento elucida che per elementi di spostamento di basso ordine (dove $F$ è costante per elemento), un'interpolazione continua di $\Theta$ è essenziale. Essa trasferisce i salti inter-elemento di $F$ in gradienti non nulli di $\Theta$, attivando la regolarizzazione. Al contrario, i campi $\Theta$ discontinui non possono generare questi gradienti inter-elemento, indebolendo l'effetto di stabilizzazione per gli elementi di basso ordine.
• Condensazione Statica: Per i campi ausiliari discontinui, il metodo dimostra che le incognite aggiuntive possono essere condensate localmente, evitando un aumento della dimensione del sistema globale pur mantenendo la stabilizzazione locale.

Risultati Numerici
La formulazione è stata valutata attraverso diversi problemi di benchmark:

• Compressione a Due Blocchi: Ha dimostrato che il parametro di scala del terzo mezzo $\gamma$ controlla il compromesso tra gap residuo (accuratezza del contatto) e robustezza della convergenza non lineare. Un $\gamma$ più piccolo produce un contatto più stretto ma richiede più iterazioni di Newton. Il metodo è rimasto stabile sotto compressione e taglio combinati.
• Auto-contatto all'interno di una Scatola:
  • Confronto di Interpolazione: I campi ausiliari continui ($T_u^1 T_\Theta^1$, $T_u^2 T_\Theta^2$) hanno fornito una chiusura del contatto robusta. La coppia mista $T_u^2 T_\Theta^1$ (spostamento quadratico, ausiliario lineare) ha offerto risultati comparabili a $T_u^2 T_\Theta^2$ con meno gradi di libertà ausiliari.
  • Discontinuo vs. Continuo: I campi ausiliari discontinui hanno mostrato una ridotta efficacia della regolarizzazione, in particolare per gli elementi lineari dove il termine del gradiente svaniva o non riusciva a accoppiare gli elementi.
  • Confronto con la Media del Gradiente di Deformazione: Rispetto al metodo di Faltus et al. [18], la proposta formulazione mista ha richiesto significativamente meno passi di carico e iterazioni di Newton, dimostrando una robustezza non lineare superiore. Sebbene la formulazione mista presenti un sistema algebrico più grande, il miglioramento della convergenza compensa il costo per mesh di dimensioni moderate, sebbene il costo dei gradi di libertà ausiliari diventi dominante per mesh estremamente fini ($>10^6$ incognite).
• Scatola a Forma di C: La formulazione ha catturato con successo grandi rotazioni, auto-contatto progressivo e severa compressione del terzo mezzo. La risposta forza-spostamento è coincisa strettamente con le formuliche classiche di contatto nodo-segmento, con transizioni fluide verso il contatto.
• Estensione 3D: Un esempio di auto-contatto in tre dimensioni ha confermato l'applicabilità del metodo oltre gli scenari di deformazione piana.

Significatività e Rivendicazioni
Il documento afferma che la stabilizzazione tramite campo ausiliario proposta fornisce un'alternativa robusta e compatibile con ordini bassi agli algoritmi di contatto classici e alle regolarizzazioni TMC esistenti. Evitando la valutazione diretta delle derivate seconde dello spostamento, essa consente l'uso di elementi finiti semplici e computazionalmente economici (lineari e quadratici) per problemi di contatto complessi che coinvolgono grandi deformazioni e auto-contatto.

Gli autori sottolineano che l'efficacia del metodo è intrinsecamente legata alla continuità del campo ausiliario; le interpolazioni continue sono necessarie per generare il meccanismo di stabilizzazione basato sul gradiente per gli elementi di spostamento di basso ordine. La formulazione è particolarmente adatta ad applicazioni dove i confini di contatto sono ignoti o evolvono, come nei meccanismi assistiti dal contatto, nella robotica soft e nell'ottimizzazione topologica, dove l'approccio basato sul continuo elimina la necessità di algoritmi espliciti di ricerca del contatto e strategie di set attivo. Il lavoro conclude che, sebbene il metodo scambi un sistema algebrico locale più grande con una migliore robustezza globale, offre una via percorribile per simulare scenari di contatto severo senza la complessità di elementi di ordine superiore o tecnologie di mesh specializzate.