Caveats on formulating finite elasto-plasticity in curvilinear coordinates

Questo articolo chiarisce le complessità nella formulazione dell'elasto-plasticità finita in coordinate curvilinee, fornendo una procedura pratica basata su cambiamenti di base espliciti per gestire correttamente i termini aggiuntivi e l'anelasticità necessari per una robusta implementazione agli elementi finiti in problemi assialsimmetrici.

L'essenziale

Immagina di dover costruire un ponte o di simulare come si comporta il terreno quando viene scavato un tunnel. Per farlo, gli ingegneri usano dei "supercomputer" che dividono il mondo in milioni di piccoli pezzi (come un puzzle) per calcolare le forze e i movimenti. Questo è il metodo degli Elementi Finiti.

Di solito, questi computer pensano in "linee rette" e "angoli retti", come se il mondo fosse fatto di scatole perfette (un sistema chiamato Cartesiano). È facile, come misurare una stanza con un righello.

Ma cosa succede se il mondo non è fatto di scatole, ma di tubi, cilindri o forme curve? Pensate a un serbatoio d'acqua, a un albero che cresce o a un tunnel rotondo. Se provate a usare il righello dritto su una superficie curva, le cose si complicano. È come cercare di stendere una mappa del mondo (che è sferica) su un foglio di carta piatto: i bordi si deformano e le distanze non tornano più.

Il Problema: La Mappa che Inganna

Questo articolo scientifico parla proprio di questo: come fare i calcoli corretti quando le forme sono curve, senza commettere errori che potrebbero far crollare il ponte virtuale.

Gli autori dicono: "Se usiamo le coordinate curve (come i cerchi e le linee radiali) senza fare attenzione, il computer pensa che le cose si stiano allungando o accorciando in modo sbagliato". È come se il computer dicesse: "Oh, questo tubo si è allargato!" mentre in realtà è solo che la nostra "riga di misura" ha cambiato forma.

Le Soluzioni: I "Trasformatori" Matematici

Per risolvere questo, gli autori introducono tre concetti chiave, che possiamo immaginare come strumenti magici:

1. Il "Gradiente di Deformazione" (Il Mago del Movimento):
   Immagina di avere un elastico. Quando lo tiri, ogni punto si sposta. In un mondo piatto, basta dire "si è spostato di 1 cm". In un mondo curvo, devi anche dire come è girato il tuo righello mentre lo spostavi. Questo "gradiente" è il mago che tiene traccia di tutto il movimento, anche se il righello cambia forma. Se sbagli qui, tutto il calcolo successivo è sbagliato.

2. Lo "Shifter" (Il Traduttore):
   Questo è il nostro "ponte". Quando passi da un sistema di coordinate all'altro (da uno curvo a uno dritto), serve un traduttore che dica al computer: "Ehi, questo punto qui corrisponde a quello lì, ma attenzione, le unità di misura sono diverse!". Senza questo traduttore, il computer confonderebbe il nord con l'est.

3. Il "Jacobian" (Il Contapassi del Volume):
   Quando deformi un oggetto (come schiacciare una pallina di argilla), il suo volume cambia. In un mondo curvo, il modo in cui calcoli questo volume è diverso rispetto a un mondo piatto. Il Jacobian è come un contapassi intelligente che ti dice: "Attenzione, perché siamo su una curva, questo piccolo pezzo di spazio è in realtà più grande o più piccolo di quanto sembri". Se sbagli questo, il computer potrebbe pensare che l'argilla sia sparita o apparsa dal nulla.

L'Esempio Pratico: Il Palloncino che Esplode

Per dimostrare che la loro teoria funziona, gli autori hanno simulato un problema classico: l'espansione di una cavità.
Immagina di avere un tubo di gomma e di gonfiarlo dall'interno come un palloncino.

• Hanno fatto il calcolo usando il loro nuovo metodo "curvo" (che è veloce perché vede solo la sezione del tubo, come un'ombra).
• Hanno fatto lo stesso calcolo usando il metodo "piatto" classico (che deve simulare l'intero tubo 3D, molto più lento).

Il risultato? I due metodi hanno dato esattamente lo stesso risultato.
Questo significa che il loro nuovo metodo è:

• Corretto: Non sbaglia i calcoli anche su forme curve.
• Veloce: Risparmia molto tempo di calcolo (come guardare un'ombra invece di girare intorno a tutto l'oggetto).

Perché è Importante?

Se un ingegnere usa i vecchi metodi "piani" su problemi "curvi" complessi (come la crescita di un vaso sanguigno o la perforazione di un terreno), potrebbe ottenere risultati sbagliati senza accorgersene. Potrebbe pensare che un materiale sia sicuro quando invece sta per rompersi.

Questo articolo è come un manuale di istruzioni aggiornato per i programmatori di questi computer. Dice: "Ehi, se lavorate su forme curve, non usate le vecchie formule! Usate queste nuove regole con lo 'Shifter' e il 'Jacobian' corretto, altrimenti rischiate di costruire castelli di sabbia che crollano appena tocca la prima onda".

In sintesi: È una guida per non farsi ingannare dalla geometria curva quando si fanno calcoli ingegneristici complessi.

Sommario tecnico

Ecco una sintesi tecnica dettagliata del paper "Caveats on formulating finite elasto-plasticity in curvilinear coordinates" di Prett et al., redatta in italiano.

Titolo

Avvertenze sulla formulazione dell'elasto-plasticità a grandi deformazioni in coordinate curvilinee

1. Problema e Contesto

Il lavoro affronta una sfida fondamentale nella meccanica computazionale dei solidi: l'implementazione numerica robusta di problemi di elasto-plasticità a grandi deformazioni (finite strain) in coordinate curvilinee, in particolare nel caso di simmetria assiale.

Sebbene l'analisi tensoriale fornisca una base invariante per la meccanica del continuo, le implementazioni numeriche (come il Metodo degli Elementi Finiti - FEM) si basano su rappresentazioni matriciali espresse in basi scelte dall'utente. Quando queste basi non sono cartesiane e ortonormali (ad esempio, coordinate cilindriche), sorgono termini aggiuntivi che sono spesso assenti o impliciti nelle formulazioni cartesiane.
Il paper evidenzia che estendere naive le formulazioni cartesiane ai sistemi curvilinei porta a errori nella definizione di grandezze cinematiche critiche, come il gradiente di deformazione, il Jacobian (fattore di trasformazione volumetrica) e il tensore spostatore (shifter). Questi errori possono compromettere la consistenza delle equazioni di equilibrio e la convergenza degli algoritmi numerici, specialmente in applicazioni ingegneristiche complesse come l'espansione di cavità, la perforazione di pozzi o la modellazione di strutture come silos e torri di raffreddamento.

2. Metodologia

Gli autori adottano un approccio pratico e computazionale, evitando il formalismo geometrico delle varietà riemanniane a favore di rappresentazioni cartesiane standard con cambi di base espliciti.

• Problema Modello: Viene utilizzato il problema classico dell'espansione di una cavità in un cilindro cavo per illustrare le differenze tra le formulazioni cartesiane e curvilinee.
• Cinematica e Grandezze Chiave:
  • Viene definita rigorosamente la relazione tra il gradiente di deformazione nelle basi cartesiane e curvilinee, dimostrando che il determinante del gradiente di deformazione in coordinate curvilinee non è direttamente uguale al rapporto dei volumi senza l'inclusione dei coefficienti metrici.
  • Viene introdotto e analizzato il ruolo del tensore spostatore (shifter tensor), che collega le basi vettoriali della configurazione di riferimento a quelle correnti. In coordinate cartesiane questo è spesso l'identità, ma in coordinate curvilinee ha componenti non banali.
  • Viene trattata la derivata covariante e i simboli di Christoffel per gestire correttamente le derivate di vettori e tensori in basi variabili.
• Elasto-plasticità a grandi deformazioni:
  • Viene adottata la decomposizione moltiplicativa del gradiente di deformazione ($F = F_e \cdot F_p$) per separare le parti elastica e plastica.
  • Viene affrontata la complessità legata alla definizione del Jacobian elastico ($J_e$) e plastico ($J_p$) in presenza di deformazioni plastiche isocore. Gli autori propongono una formula per la linearizzazione consistente di $J_e$ basata sulla traccia della deformazione logaritmica elastica, evitando la necessità di calcolare derivate complesse di tensori intermedi non noti.
  • Vengono definite le misure di sforzo appropriate (incluso il tensore $\zeta = J_e \sigma$) per garantire la validità della seconda legge della termodinamica e la conservazione di una configurazione intermedia priva di stress.
• Formulazione FEM:
  • Vengono derivati le forme forti e deboli dell'equilibrio della quantità di moto lineare sia in formulazione Lagrangiana Totale (TL) che Aggiornata (UL) in coordinate curvilinee.
  • Viene presentata una procedura passo-passo per la discretizzazione e la linearizzazione (algoritmo di Newton-Raphson) specifica per problemi assialsimmetrici.

3. Contributi Chiave

1. Chiarezza Cinematica: Dimostrazione che il gradiente di deformazione in coordinate curvilinee non è un semplice gradiente vettoriale e che la sua componente "fuori piano" (nel caso assialsimmetrico) non è necessariamente unitaria se non si considerano correttamente i coefficienti metrici.
2. Ruolo dello Shifter: Evidenziazione del fatto che il tensore spostatore non è trascurabile in coordinate curvilinee e deve essere incluso esplicitamente nel calcolo del gradiente di deformazione e delle sue inverse.
3. Linearizzazione Consistente: Sviluppo di una procedura per la linearizzazione consistente del Jacobian elastico in presenza di plasticità, essenziale per la convergenza quadratica degli schemi Newton-Raphson in FEM.
4. Guida Pratica per l'Implementazione: Fornitura di algoritmi dettagliati (Algoritmi 1 e 2 nel paper) che distinguono chiaramente tra le formulazioni TL e UL, mostrando come aggiornare i coefficienti metrici, i simboli di Christoffel e i gradienti delle funzioni di forma in ogni iterazione.
5. Validazione Numerica: Confronto diretto tra una simulazione 3D completa e la formulazione assialsimmetrica proposta, dimostrando l'equivalenza dei risultati con un costo computazionale drasticamente ridotto.

4. Risultati

Il paper presenta un esempio numerico sull'espansione di un cilindro cavo spesso con materiali elasto-plastici (modello Neo-Hookean decoupled con incrudimento e legge di snervamento Modified Cam-Clay).

• Confronto 3D vs 2D Assialsimmetrico: I risultati mostrano una sovrapposizione quasi perfetta tra le soluzioni ottenute con la formulazione 3D completa e quella assialsimmetrica proposta. Grandezze invarianti come la pressione di Cauchy, lo sforzo deviatorico e i Jacobiani elastico e plastico coincidono per punti di Gauss selezionati.
• Convergenza: L'algoritmo di Newton-Raphson dimostra una convergenza robusta e simile per entrambe le formulazioni (3D e 2D), confermando la correttezza della linearizzazione proposta.
• Efficienza: La formulazione assialsimmetrica riduce il problema a una discretizzazione 2D, offrendo un risparmio computazionale significativo senza perdita di accuratezza.

5. Significato e Implicazioni

Questo lavoro è significativo per la comunità della meccanica computazionale perché:

• Previene Errori Sottili: Avverte i ricercatori e gli sviluppatori di codici FEM che l'uso di coordinate curvilinee richiede un'attenzione maniacale alla definizione delle basi e dei coefficienti metrici; ignorare questi aspetti porta a errori fondamentali nella fisica del problema.
• Abilita Simulazioni Efficienti: Fornisce le basi teoriche e pratiche per implementare modelli di grandi deformazioni in geometrie assialsimmetriche complesse (comuni in geotecnica e ingegneria civile) in modo efficiente ed accurato.
• Ponte Teoria-Pratica: Colma il divario tra la teoria tensoriale astratta e l'implementazione pratica nel codice, offrendo una "ricetta" chiara per gestire la cinematica, la costituzione e la linearizzazione in contesti non cartesiani.

In sintesi, il paper stabilisce che per ottenere analisi robuste di problemi elasto-plastici a grandi deformazioni in simmetria assiale, è imperativo abbandonare le semplificazioni cartesiane e adottare una formulazione tensoriale completa che includa esplicitamente i coefficienti metrici, lo spostatore e le corrette relazioni tra le configurazioni di riferimento, corrente e intermedia.