⚡ electrical engineering

A discrete‑natured 3D DEM framework for spatial fracture analysis of unidirectional CFRP composites

Questo studio presenta un nuovo framework di metodo a elementi discreti 3D con uno schema di impacchettamento cubico che tratta le celle unitarie come elementi continui per modellare efficacemente la frattura trasversale nei compositi CFRP unidirezionali, rivelando che mentre la distribuzione delle fibre influenza minimamente la risposta meccanica globale, la frazione di volume e il raggio delle fibre governano criticamente la rigidezza, la resistenza e la localizzazione del danno.

Autori originali: Dechao Sun, Xiaoxuan Ding, Shuqi Zhang, Jianlin Liu, Zewen Gu, Xiaonan Hou

Pubblicato 2026-09-11

📖 5 min di lettura🧠 Approfondimento

CC BY 4.0

Autori originali: Dechao Sun, Xiaoxuan Ding, Shuqi Zhang, Jianlin Liu, Zewen Gu, Xiaonan Hou

Articolo originale sotto licenza CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Questa è una spiegazione generata dall'IA dell'articolo qui sotto. Non è stata scritta né approvata dagli autori. Per precisione tecnica, consulta l'articolo originale. Leggi il disclaimer completo

L'ingegneria moderna si affida pesantemente a materiali che siano resistenti ma leggeri, come i polimeri rinforzati con fibra di carbonio. Questi compositi sono realizzati inserendo sottili e rigidi filamenti di carbonio in una colla plastica più morbida, creando una struttura in grado di sopportare carichi pesanti nelle applicazioni aerospaziali e automobilistiche. Tuttavia, come ogni materiale, essi hanno un punto debole: quando vengono tirati lateralmente, la colla plastica tra le fibre può creparsi. Questo tipo di frattura, nota come frattura trasversale, spesso inizia come minuscole imperfezioni invisibili che crescono e si collegano tra loro, causando infine la perdita di resistenza dell'intera struttura. Comprendere esattamente come queste crepe si formino e si muovano attraverso il complesso labirinto di fibre e colla è difficile perché il danno avviene su una scala troppo piccola per essere vista facilmente con gli strumenti standard, e le crepe non percorrono linee rette. Esse si torcono, virano e interagiscono con le fibre in uno spazio tridimensionale, rendendo difficile prevedere quando un componente possa cedere.

Per risolvere questo enigma, i ricercatori Dechao Sun e il suo team della China University of Petroleum, lavorando con Xiaonan Hou della Lancaster University, hanno sviluppato un nuovo modo per simulare questi materiali su un computer. Invece di trattare il materiale come un blocco solido e liscio, lo hanno scomposto in una vasta collezione di minuscole particelle individuali che si toccano e si spingono a vicenda. Questo approccio, chiamato metodo degli elementi discreti, permette al computer di osservare il comportamento del materiale come se fosse composto da milioni di microscopici mattoni costruttivi. Il team ha creato un modello tridimensionale specifico di un cubo di questo materiale composito, riempiendolo con fibre di carbonio virtuali e resina epossidica. Hanno poi ideato una nuova regola matematica per collegare la rigidità delle minuscole connessioni tra queste particelle alla resistenza complessiva del materiale. Un'innovazione chiave del loro metodo è che il modello genera automaticamente la tendenza del materiale a gonfiarsi quando viene schiacciato — una proprietà nota come rapporto di Poisson — semplicemente osservando come si muovono le particelle, invece di costringere il computer a usare un numero preimpostato.

Quando i ricercatori hanno testato questo nuovo modello, hanno scoperto che poteva riprodurre accuratamente il comportamento reale del materiale. La simulazione ha previsto che il materiale si sarebbe allungato e rotto in modi che corrispondevano agli esistenti esperimenti di laboratorio e ai precedenti modelli informatici bidimensionali. Tuttavia, la visione tridimensionale ha rivelato qualcosa che i modelli più vecchi avevano tralasciato: le crepe non tagliano il materiale tutto in una volta come un coltello attraverso una pagnotta. Invece, il fronte della crepa avanza in modo irregolare. In alcuni punti, la crepa avanza rapidamente, mentre in altri rallenta, creando un percorso seghettato e tortuoso attraverso lo spessore del materiale. Questa progressione irregolare, che i ricercatori chiamano tortuosità del fronte della crepa, è un dettaglio cruciale che solo una visione tridimensionale completa può catturare.

Il team ha poi utilizzato questo strumento per investigare come diverse disposizioni delle fibre influenzino la resistenza del materiale. Hanno testato se il posizionamento casuale delle fibre all'interno della colla plastica cambiasse la resistenza del materiale. Hanno scoperto che, sebbene il percorso esatto seguito da una crepa dipendesse fortemente da dove le fibre si trovassero effettivamente, la resistenza e la rigidità complessive del materiale rimanevano quasi invariate indipendentemente dalla disposizione. Ciò suggerisce che, per prevedere le prestazioni generali del composito, il particolare schema casuale delle fibre conta meno di quanto precedentemente ipotizzato. Tuttavia, la quantità di fibra impacchettata nel materiale faceva una differenza significativa. Aumentare il contenuto di fibra dal 50 per cento al 68 per cento rendeva il materiale più rigido e resistente, ma lo rendeva anche più fragile, causandone il cedimento con un allungamento inferiore. I ricercatori hanno osservato che, quando le fibre erano impacchettate troppo strettamente, non c'era abbastanza colla tra di esse per assorbire lo stress, portando a una rapida fratturazione.

Infine, i ricercatori hanno esaminato le dimensioni delle fibre stesse, mantenendo costante la quantità di fibra. Hanno scoperto che il raggio delle fibre, ovvero quanto fossero spesse, influenzava il modo in cui il materiale gestiva lo stress. Un raggio della fibra di 3,5 micrometri si è rivelato la configurazione più efficiente, offrendo il miglior equilibrio tra il trasferimento del carico e la capacità della colla di allungarsi. Le simulazioni hanno mostrato che la dimensione delle fibre dettava la spaziatura tra di esse, che a sua volta controllava dove si concentrava lo stress e dove iniziavano le crepe. Mappando le forze tra le particelle virtuali, il team ha potuto vedere esattamente come il carico passasse dalla colla alle fibre e dove il materiale cedeva finalmente. Questo nuovo quadro tridimensionale fornisce uno strumento potente agli ingegneri per comprendere la meccanica nascosta dei materiali compositi, offrendo un quadro più chiaro di come queste strutture complesse falliscano e di come possano essere progettate per durare più a lungo.

Sintesi Tecnica: Un Framework DEM 3D a Natura Discreta per l'Analisi della Frattura Spaziale di Compositi CFRP Unidirezionali

Problematica
I polimeri rinforzati con fibra di carbonio (CFRP) sono critici nelle applicazioni aerospaziali e automobilistiche, ma soffrono di una limitata tolleranza al danno, in particolare per quanto riguarda la fessurazione trasversale. Questa modalità di cedimento avvia microfessure che si propagano attraverso lo spessore del ply, portando alla degradazione della rigidezza e al compromesso strutturale. Sebbene gli approcci numerici esistenti come la Meccanica del Danno Continuo (CDM), i Modelli a Zona Coesiva (CZM) e il Metodo degli Elementi Finiti Esteso (XFEM) abbiano fatto progressi nel campo, essi spesso dipendono da percorsi di fessurazione predefiniti o faticano con la convergenza e la complessità di molteplici fessure 3D intersecanti. Esistono inoltre studi sul Metodo degli Elementi Discreti (DEM) che sono prevalentemente limitati ad analisi di sezione trasversale 2D, fallendo nel catturare la continuità della fessura attraverso lo spessore, la connettività del danno 3D e la ridistribuzione spaziale dello stress inerente alle microstrutture eterogenee. Rimane un vuoto nei modelli DEM 3D espliciti capaci di simulare la fessurazione trasversale nei CFRP unidirezionali senza fare affidamento su proprietà anisotrope equivalenti o strutture omogeneizzate.

Metodologia
Gli autori propongono un nuovo framework a metodo degli elementi discreti tridimensionale (3D DEM) per investigare la frattura trasversale in compositi epossidici rinforzati con fibra di carbonio. La metodologia si basa sui seguenti pilastri:

Modellazione Discreta: Un Elemento Rappresentativo di Volume (RVE) cubico di 49 µm viene costruito utilizzando uno schema di impacchettamento cubico 3D. La microstruttura consiste in 38 fibre circolari (raggio 3,5 µm) immerse in una matrice epossidica, con interazioni modellate tramite molle a forza centrale lungo sei direzioni ortogonali.
Calibrazione a Natura Discreta: A differenza dei tradizionali modelli a particelle legate che trattano i legami come elementi del continuo richiedendo la calibrazione simultanea delle rigidezze normale ( $k_n$ $k_{n}$ ) e tangenziale ( $k_s$ $k_{s}$ ), questo studio introduce una innovativa relazione di equivalenza energetica. Il modello tratta ogni cella unitaria come un elemento continuo dove il coefficiente di Poisson macroscopico è generato automaticamente dalla deformazione della cella anziché essere prescritto.
- Il framework deriva relazioni analitiche esplicite (Eq. 41 e 42) che collegano la rigidezza normale microscopica ( $k_n$ ) al modulo di Young macroscopico ( $E$ ) e al coefficiente di Poisson ( $\nu$ ).
- Questo approccio elimina la necessità della rigidezza tangenziale ( $k_s$ ) nel processo di calibrazione, affidandosi esclusivamente alla rigidezza normale per soddisfare l'equivalenza energetica tra le forze di contatto microscopiche e la deformazione del continuo macroscopico.
Validazione: Il modello calibrato è validato rispetto a dati sperimentali e risultati DEM 2D. Le simulazioni sono condotte sotto carico di trazione quasi-statico a spostamento controllato per osservare l'evoluzione del danno, la propagazione delle fessure e le risposte sforzo-deformazione macroscopiche.

Risultati Chiave
Lo studio conduce indagini parametriche sugli effetti della distribuzione casuale delle fibre, della frazione di volume delle fibre ( $V_f$ ) e del raggio della fibra ( $r$ ):

Validazione del Modello: Il 3D DEM riproduce con successo il comportamento meccanico macroscopico dell'epossidica pura e degli RVE di CFRP, con il modulo di Young e il coefficiente di Poisson simulati che deviano di meno del 2% dai valori target. Il modello cattura la continuità della fessura attraverso lo spessore e la tortuosità del fronte della fessura, mostrando che la propagazione della fessura non è simultanea lungo la lunghezza della fibra, a differenza delle traiettorie planari dei modelli 2D.
Distribuzione delle Fibre: Mentre le distribuzioni casuali delle fibre alterano significativamente i siti locali di innesco delle fessure e i percorsi finali delle stesse (ad esempio, fluttuando tra fibre strettamente impacchettate), esse hanno un'influenza trascurabile sulle risposte meccaniche globali. Gli RVE statisticamente equivalenti esibiscono comportamenti sforzo-deformazione e l'evoluzione del modulo di Young quasi identici.
Frazione di Volume delle Fibre ( $V_f$ ): L'aumento di $V_f$ (dal 50% al 68%) migliora la rigidezza trasversale e la resistenza di picco, ma riduce la deformazione a rottura, portando a una risposta post-picco più fragile. Un $V_f$ più elevato induce carenza di resina, promuovendo concentrazioni di stress localizzate e un danno precoce dominato dalla matrice.
Raggio della Fibra: A una frazione di volume fissa ( $V_f = 60\%$ ), il raggio della fibra governa l'equilibrio tra il trasferimento di stress interfacciale e l'accomodamento della deformazione della matrice. La configurazione con un raggio della fibra di $r = 3,5$ µm ha dimostrato la capacità di carico più efficiente e la rigidezza più elevata rispetto a $r = 3,2$ µm e $r = 3,8$ µm.

Significatività e Rivendicazioni
L'articolo sostiene che il framework 3D DEM proposto offre uno strumento efficace per l'analisi del danno progressivo multiscala nei compositi CFRP unidirezionali. La sua principale significatività risiede in:

Calibrazione a Natura Discreta: Fornire una mappatura rigorosa, basata sull'energia, da parametri microscopici a proprietà macroscopiche che rispetti la natura discontinua del DEM pur generando naturalmente il coefficiente di Poisson.
Analisi della Frattura Spaziale: Catturare con successo complessi meccanismi di cedimento 3D, inclusa la penetrazione della fessura attraverso lo spessore, la tortuosità del fronte della fessura e la coalescenza del danno, che sono irraggiungibili con le semplificazioni 2D.
Intuizione Meccanicistica: Elucidare che, mentre la distribuzione delle fibre detta i pattern di frattura locale, la frazione di volume delle fibre è il fattore dominante che controlla le prestazioni meccaniche globali.

Gli autori concludono che, sebbene l'attuale modello affronti efficacemente la fessurazione trasversale nei laminati unidirezionali, sono necessari lavori futuri per incorporare layup di fibre multidirezionali, relazioni costitutive raffinate dell'interfaccia fibra-matrice e strategie di accoppiamento multiscala del danno.

Articoli simili