A multi-phase-field model for fiber-reinforced composite laminates based on puck failure theory

Questo articolo propone un modello multi-campo di fase basato sulla teoria del fallimento di Puck, che utilizza un metodo di sovrapposizione di mesh per prevedere accuratamente i danni da fibra e inter-fibra nei laminati compositi rinforzati con fibre, convalidando i risultati attraverso confronti quantitativi e qualitativi con dati sperimentali su diversi casi di carico.

L'essenziale

Immagina di dover costruire un ponte o l'ala di un aereo usando un materiale incredibilmente leggero e resistente, ma anche molto complicato: le fibre composite (come la fibra di carbonio). Questi materiali sono fatti di strati sottilissimi, come una lasagna, dove ogni strato ha le fibre orientate in una direzione diversa per resistere a forze diverse.

Il problema è: come facciamo a sapere esattamente quando e dove questo materiale si romperà?

Fino a poco tempo fa, gli ingegneri dovevano costruire centinaia di prototipi fisici, romperli in laboratorio e sperare di capire il perché. Questo è costoso, lento e non sempre sicuro.

Questo articolo presenta un nuovo software di simulazione (un "cervello digitale") che permette di prevedere la rottura di questi materiali con grande precisione, senza doverli rompere fisicamente. Ecco come funziona, spiegato con parole semplici e analogie.

1. Il Problema: La "Lasagna" che si spacca

Pensa a un pannello in fibra di carbonio come a una lasagna.

• Ogni strato (sfoglia) è fatto di fibre.
• Se tiri la lasagna, può rompersi in due modi principali:
  1. Le fibre si spezzano (come se rompessi i fili di spaghetto).
  2. La colla tra le fibre si rompe (come se la sfoglia si staccasse dal ripieno).

Spesso, quando un materiale si rompe, questi due eventi accadono insieme o uno dopo l'altro in modo molto complesso. I vecchi modelli di computer erano come bambini che vedono solo "tutto o niente": o il materiale è intero, o è rotto. Non capivano le sfumature.

2. La Soluzione: Il "Termostato della Rottura" (Teoria di Puck)

Gli autori di questo studio hanno creato un modello che usa una teoria chiamata Teoria di Puck.
Immagina che ogni strato della tua lasagna abbia un piccolo termostato interno che controlla due cose separate:

• Un termostato per le fibre.
• Un termostato per la colla (matrice).

Questi termostati non si limitano a dire "rottura". Calcolano esattamente quanto stress sta subendo il materiale in ogni direzione. Se la colla inizia a cedere, il termostato della colla si attiva. Se le fibre iniziano a cedere, quello delle fibre si attiva. In questo modo, il computer capisce esattamente quale tipo di danno sta avvenendo e in che ordine.

3. La Magia: Il "Metodo dell'Overlay" (Sovrapposizione)

Qui arriva la parte più creativa. Normalmente, per simulare una lasagna di 10 strati, dovresti costruire un modello 3D complesso con 10 livelli di mattoncini digitali. È lento e pesante.

Questi ricercatori hanno usato un trucco geniale chiamato "Metodo dell'Overlay" (Sovrapposizione):

• Immagina di prendere un unico foglio di carta digitale (la base).
• Poi, invece di costruire 10 strati sopra, sovrapponi 10 "fantasmi" digitali sullo stesso foglio.
• Ogni "fantasma" rappresenta uno strato della lasagna con la sua direzione di fibre specifica.
• Tutti i fantasmi condividono gli stessi punti di ancoraggio (i nodi), ma ognuno ha le sue regole interne.

È come se avessi 10 persone diverse che indossano lo stesso vestito, ma ognuna pensa e reagisce in modo diverso agli stress. Questo rende la simulazione velocissima e permette di vedere come si propagano le crepe all'interno di ogni singolo strato, anche se sono tutti nello stesso spazio.

4. Il Campo Diffuso: La "Neve che Scioglie"

Un altro concetto chiave è il Campo di Fase.
Nella realtà, una crepa è una linea netta e tagliente. Nel computer, trattare una linea netta è un incubo matematico.
Quindi, gli autori usano un'analogia con la neve che sta sciogliendo:

• Invece di una linea netta, la crepa è una "zona di transizione" sfumata, come un'area dove la neve sta diventando fango.
• Il computer calcola quanto è "fango" (danneggiato) ogni punto del materiale.
• Questo permette al software di gestire la rottura in modo fluido, senza che il programma si blocchi quando la crepa cambia direzione o si ramifica.

5. Cosa hanno scoperto? (I Test)

Hanno testato questo "cervello digitale" su quattro scenari reali, confrontando i risultati del computer con esperimenti fisici reali:

1. Trazione e Compressione: Hanno tirato e schiacciato strisce di materiale. Il modello ha previsto esattamente quando si sarebbero rotti.
2. Buco nel mezzo (Tensione con foro): Come quando hai un buco in una lastra di plastica. Il modello ha visto esattamente come le crepe partono dal buco e si allargano (spesso seguendo le fibre).
3. Trazione Compatta: Un test per vedere quanto il materiale resiste a essere strappato.
4. Doppio Taglio: Materiali con due tagli ai lati.

In tutti i casi, il modello ha imitato perfettamente il comportamento reale, sia nel modo in cui si rompe (qualitativamente) sia nella forza necessaria per romperlo (quantitativamente).

In Sintesi

Questo articolo ci dice che abbiamo finalmente un modo per "vedere il futuro" dei materiali compositi.
Grazie a questo modello:

• Possiamo progettare aerei e auto più leggeri e sicuri.
• Sappiamo esattamente dove si formerà la crepa prima che accada.
• Risparmiamo tempo e denaro perché non dobbiamo costruire e rompere così tanti prototipi fisici.

È come avere una palla di cristallo ingegneristica che ti dice: "Attenzione, se carichi questo pezzo in questo modo, la colla tra le fibre cederà qui, e poi le fibre si spezzeranno lì". Tutto questo grazie a un po' di matematica intelligente e a un trucco di sovrapposizione degli strati digitali.

Sommario tecnico

Ecco un riassunto tecnico dettagliato del documento in lingua italiana.

Titolo: Un modello multi-campo di fase per laminati in compositi rinforzati con fibre basato sulla teoria del fallimento di Puck

1. Problema e Contesto

I compositi rinforzati con fibre (FRC) sono ampiamente utilizzati in settori critici come l'aerospaziale e l'automotive grazie al loro elevato rapporto resistenza/peso. Tuttavia, la previsione della loro frattura rimane una sfida significativa a causa dei complessi meccanismi di fallimento, che includono:

• Criccatura della matrice intralaminare (trasversale e longitudinale).
• Delaminazione (fallimento interlaminare).
• Fallimento delle fibre (trazione e compressione).
• Interazione tra diversi modi di fallimento.

I metodi numerici esistenti, come gli elementi coesivi o le tecniche XFEM, spesso richiedono una discretizzazione complessa o non catturano adeguatamente la nucleazione delle cricche e l'interazione tra i diversi meccanismi di danno in configurazioni di laminati con sequenze di stratificazione complesse. Inoltre, i modelli di campo di fase (Phase-Field, PF) tradizionali applicati ai compositi sono spesso limitati a singoli strati o sequenze semplici (es. cross-ply) e faticano a gestire correttamente i carichi di compressione e la nucleazione delle cricche senza criteri aggiuntivi.

2. Metodologia Proposta

Gli autori propongono un modello multi-campo di fase (Multi-Phase-Field) bidimensionale che integra la teoria del fallimento di Puck per predire il danno nei laminati FRC.

• Variabili di Campo di Fase: Il modello introduce due variabili di campo di fase indipendenti:
  • $d_f$: per il fallimento dominato dalle fibre.
  • $d_m$: per il fallimento dominato dalla matrice (inter-fibra).
    Questo permette di attivare e modellare separatamente i meccanismi di danno, catturando la loro interrelazione.
• Teoria di Puck: La nucleazione del danno non è basata su una semplice energia di deformazione, ma è attivata dai criteri di fallimento di Puck. Questi criteri distinguono tra:
  • Fallimento delle fibre (trazione e compressione).
  • Fallimento della matrice in tre modalità (A: trazione trasversale, B: compressione trasversale + taglio, C: compressione trasversale dominante).
    Il modello calcola i "fattori di esposizione" ($F_i$) basati sugli stati di stress locali; quando $F_i \geq 1$, il danno inizia.
• Metodo di Sovrapposizione della Mesh (Mesh Overlay): Per gestire laminati con diverse orientazioni delle fibre senza dover modellare esplicitamente ogni strato con una mesh 3D complessa (e senza simulare la delaminazione), il modello utilizza una tecnica di sovrapposizione.
  • Ogni strato (ply) è rappresentato da una mesh separata ma condivisa (nodi comuni).
  • Le variabili di spostamento sono accoppiate tra gli strati, mentre le tensioni e i danni sono decoupled e dipendono dall'orientamento delle fibre specifico di ogni strato.
  • Questo approccio rende il modello computazionalmente efficiente per l'analisi 2D in piano.
• Formulazione Variazionale: Il sistema è formulato in base alla minimizzazione dell'energia potenziale totale, che include l'energia di deformazione degradata e l'energia di superficie (danno). La consistenza termodinamica è garantita attraverso le condizioni di Karush-Kuhn-Tucker (KKT) e l'uso di forze motrici basate sull'energia critica di fallimento di Puck.
• Implementazione: Il modello è stato implementato come un elemento utente (UEL) in Abaqus, utilizzando il metodo degli elementi finiti.

3. Contributi Chiave

• Integrazione Puck-Phase-Field: È uno dei primi modelli che applica rigorosamente la teoria di Puck all'interno di un framework di campo di fase per gestire la nucleazione e la propagazione del danno in laminati complessi.
• Gestione Separata dei Modi di Danno: L'uso di due campi di fase indipendenti permette di distinguere chiaramente tra la rottura delle fibre e quella della matrice, superando le limitazioni dei modelli a singolo campo.
• Efficienza Computazionale per Laminati: L'uso della tecnica "mesh overlay" permette di simulare laminati con sequenze di stratificazione complesse in 2D, evitando la complessità computazionale dei modelli 3D pur catturando l'effetto delle diverse orientazioni degli strati.
• Validazione Estensiva: Il modello è stato testato su quattro casi di riferimento (benchmark) confrontando i risultati numerici con dati sperimentali della letteratura.

4. Risultati e Validazione

Il modello è stato validato su quattro configurazioni geometriche e di carico:

1. Test su Campioni (Coupon Tests):

   • Simulazione di carichi di trazione e compressione su laminati angolari [(θ/-θ)8]s.
   • I risultati mostrano un buon accordo quantitativo con gli esperimenti (Flatscher et al.) per quanto riguarda i carichi massimi, sebbene il modello tenda a sovrastimare leggermente la resistenza a causa dell'assenza di plasticità.
   • Viene riprodotta correttamente la propagazione instabile delle cricche dopo l'iniziazione.

2. Trazione con Foro Aperto (Open Hole Tension - OHT):

   • Testati laminati [0°4/90°4]s e [45°4/-45°4]s.
   • Per il caso [0/90], il modello cattura la propagazione simmetrica del danno della matrice seguita dalla rottura delle fibre, con un picco di carico in ottimo accordo con l'esperimento (24.3 kN vs 25.2 kN sperimentali).
   • Per il caso [45/-45], il modello riproduce correttamente l'asimmetria della propagazione delle cricche dovuta all'effetto di "schermatura delle cricche" (crack shielding), un fenomeno difficile da catturare con modelli semplificati.

3. Tensione Compatta (Compact Tension - CT):

   • Confronto su laminati bloccati [0°2/90°2]2s e [0°4/90°4]2s.
   • Il modello riproduce bene il comportamento lineare iniziale e la risposta forza-spostamento.
   • Le discrepanze nei carichi finali sono attribuite alla mancata simulazione della delaminazione (che nel modello reale dissipa energia), ma la morfologia delle cricche è qualitativamente corretta.

4. Trazione con Doppia Intaccatura (Double-Edged Notched Tension):

   • Analisi di laminati cross-ply [0/90]s, [45/0/-45]s e quasi-isotropi [45/90/-45/0]s.
   • Il modello dimostra la capacità di simulare l'interazione tra fallimento delle fibre e della matrice, la nucleazione di cricche sub-critiche e la propagazione complessa in diverse direzioni.

5. Significatività e Conclusioni

Questo lavoro rappresenta un avanzamento significativo nella modellazione computazionale dei compositi strutturali.

• Capacità Predittiva: Il modello offre una valutazione sia qualitativa (morfologia delle cricche) che quantitativa (carichi di rottura) affidabile per una vasta gamma di sequenze di laminazione.
• Limitazioni: Essendo un modello 2D, non include la delaminazione (fallimento interlaminare) e non considera la plasticità o l'incrudimento del materiale. Di conseguenza, è più adatto per l'analisi di fallimento in piano e per la previsione della resistenza, mentre richiede calibrazione dei parametri di guida (dimensionless driving parameters) per la precisione quantitativa assoluta.
• Impatto: La disponibilità del codice e dei dati di ricerca favorisce la trasparenza e permette alla comunità scientifica di costruire su questa base per sviluppare modelli ancora più avanzati.

In sintesi, l'articolo presenta un framework robusto che combina la fisica del fallimento di Puck con l'efficienza del metodo del campo di fase, offrendo uno strumento potente per la progettazione e l'analisi di laminati compositi avanzati.