A 2D discrete model with a phase-field regularisation applied to dynamic crack propagation

Quando un materiale fragile come il vetro o l'acciaio temprato si rompe, non si limita a spezzarsi semplicemente in due. Al contrario, una crepa corre attraverso il solido, spesso ramificandosi in una complessa rete di fratture che frantuma l'oggetto in molti pezzi. Prevedere esattamente dove queste crepe viaggeranno e quanto velocemente si muoveranno è una delle sfide più difficili dell'ingegneria. Se un modello utilizzato per simulare questa distruzione è difettoso, potrebbe suggerire che un ponte crollerà in una direzione sicura quando invece collasserà catastroficamente. Per decenni, gli scienziati hanno cercato di costruire modelli informatici che imitino la struttura atomica dei materiali, trattandoli come una griglia di minuscoli punti collegati da molle. Sebbene questi modelli siano potenti, soffrono di un difetto peculiare: il percorso che segue la crepa dipende spesso da come la griglia del computer viene ruotata, piuttosto che dalla fisica del materiale stesso. È come se una crepa in un vero pezzo di vetro cambiasse improvvisamente direzione semplicemente perché l'ingegnere ha ruotato lo schermo del computer.

Per risolvere questo problema, Matías Braun, un ricercatore della Universidad Nebrija in Spagna, ha sviluppato un nuovo modo per simulare la frattura dinamica. Il suo approccio combina una griglia tradizionale di punti con un metodo a "campo di fase" (phase-field), una tecnica che tratta la crepa non come una linea netta che salta da un punto all'altro, ma come una zona di danno fluida e sfumata che si diffonde attraverso il materiale. Fondendo queste due idee, Braun ha creato un modello che permette alle crepe di viaggiare in qualsiasi direzione, indipendentemente dall'orientamento della griglia informatica sottostante. Ha testato questo nuovo metodo su tre scenari distinti: una piastra d'acciaio colpita da un impatto ad alta velocità, una piastra con un intaglio tirata fino a spezzarsi e un cilindro metallico che esplode dall'interno verso l'esterno. In ogni caso, il nuovo modello ha prodotto risultati che corrispondevano agli esperimenti del mondo reale ed eliminavano il pregiudizio artificiale causato dalla griglia del computer, offrendo una visione più chiara e affidabile di come i materiali falliscono realmente.

Il cuore del problema affrontato da Braun risiede nel modo in cui i computer gestiscono il mondo invisibile degli atomi. Nel suo lavoro precedente, utilizzava un modello discreto in cui il materiale era rappresentato da un reticolo triangolare, un insieme di punti collegati da linee. Quando il materiale veniva sollecitato, i collegamenti tra i punti si rompevano, simulando una crepa. Tuttavia, poiché i punti erano disposti in un modello fisso, la crepa era costretta a seguire le linee della griglia. Se la griglia veniva ruotata, la crepa cambiava il suo percorso, anche se il materiale stesso non era cambiato. Questo è un problema significativo perché, nella realtà, una crepa dovrebbe seguire il percorso di minor resistenza determinato dalle proprietà del materiale, non l'angolo arbitrario di una mesh informatica.

La soluzione di Braun è stata l'introduzione di una regolarizzazione a campo di fase. Invece di forzare la rottura di specifici collegamenti tra i punti, ha assegnato a ogni punto un valore che rappresenta quanto sia danneggiato. Questo valore varia da zero, che significa che il materiale è perfettamente intatto, a uno, che significa che è completamente rotto. La crepa diventa una regione in cui questo valore aumenta gradualmente, creando una transizione fluida piuttosto che una linea frastagliata. Ciò consente alla crepa di curvarsi e ramificarsi in qualsiasi direzione, ignorando efficacementamente la struttura rigida della griglia. Fondamentalmente, Braun ha calibrato questo sistema utilizzando due proprietà fondamentali del materiale: la sua resistenza alla trazione, ovvero la forza necessaria per iniziare una crepa, e l'energia di frattura, ovvero l'energia necessaria per far crescere la crepa. Ancorando la simulazione a queste costanti fisiche reali, il modello assicura che i risultati dipendano dal materiale, non dalla matematica.

Per dimostrare che questo nuovo metodo funzionasse, Braun ha eseguito tre test rigorosi. Il primo è stato il test Kalthoff–Winkler, un esperimento classico in cui una piastra d'acciaio con un intaglio pre-tagliato viene colpita da un proiettile ad alta velocità. In questo scenario, una crepa scatta dall'intaglio con un angolo specifico. I modelli precedenti spesso prevedevano l'angolo errato a seconda di come la griglia veniva ruotata. Nelle simulazioni di Braun, la crepa ha viaggiato costantemente a un angolo di circa 70 gradi, corrispondendo perfettamente alle osservazioni sperimentali. Questo risultato è rimasto valido sia che la griglia fosse allineata con la piastra, sia che fosse ruotata di 15 o 30 gradi. Anche la velocità della crepa è rimasta stabile, viaggiando a circa il 70 percento della velocità del suono nel materiale, indipendentemente dall'orientamento della griglia.

Il secondo test ha coinvolto una piastra con un intaglio a bordo singolo sotto tensione, una configurazione utilizzata per studiare come le crepe crescano lentamente sotto una trazione costante. Qui, l'obiettivo era vedere se la crepa avrebbe viaggiato in linea retta dall'intaglio al bordo opposto. I modelli a griglia tradizionali spesso costringevano la crepa a procedere a zigzag o a deviare in base all'allineamento della griglia. Il modello di Braun, tuttavia, ha prodotto una crepa dritta per ogni orientamento della griglia. Anche la forza necessaria per rompere la piastra è stata coerente, variando di meno del due percento indipendentemente da come la griglia veniva ruotata. Ciò ha dimostrato che il modello può prevedere accuratamente sia il percorso che la resistenza del materiale senza l'influenza artificiale della mesh informatica.

Il terzo test, il più complesso, è stato la frammentazione di un cilindro spesso sottoposto a un'esplosione interna. In questo scenario, nessuna singola crepa è guidata da un intaglio; invece, molte crepe devono nucleare spontaneamente e correre verso l'esterno, frantumando il cilindro in dozzine di pezzi. Questo è un evento caotico in cui la simmetria della griglia potrebbe facilmente dettare il pattern dei frammenti. Quando Braun ha eseguito la simulazione con la griglia centrata sull'asse del cilindro, le crepe hanno formato un pattern simmetrico. Tuttavia, quando ha spostato leggermente la griglia in modo che non fosse più perfettamente centrata, le crepe si sono distribuite in modo irregolare e il numero di frammenti è cambiato. Ciò ha confermato che il modello era sensibile alla configurazione fisica piuttosto che alla semplice simmetria della griglia. La simulazione ha previsto tra i 12 e i 19 frammenti, un numero che rientra nel range di altri metodi avanzati, ma che è inferiore ad alcuni, probabilmente perché la zona di danno fluida nel modello impedisce la formazione di crepe troppo vicine tra loro.

Durante tutti questi test, i ricercatori hanno monitorato attentamente l'energia del sistema. Hanno confrontato l'energia persa mentre il materiale si rompeva rispetto all'energia necessaria per creare le nuove superfici di frattura. In alcuni casi, c'è stata una piccola discrepanza tra l'energia calcolata dal modello di danno e l'energia effettivamente rilasciata dai legami rotti. Questa è una limitazione nota dell'approccio ibrido, in cui la forza motrice della crepa è calcolata diversamente dall'energia dei legami rotti. Tuttavia, il bilancio dell'energia meccanica è rimasto accurato entro lo 0,5%, confermando che la simulazione è fisicamente solida. I risultati suggeriscono che, sebbene il modello non sia perfetto, rappresenta un passo avanti significativo nella comprensione della frattura dinamica.

Le implicazioni di questo lavoro vanno oltre la semplice migliore grafica computerizzata. Rimuovendo la dipendenza dall'orientamento della griglia, gli ingegneri possono fidarsi del fatto che le loro simulazioni di crash test, cedimenti strutturali o test sui materiali riflettano la vera fisica del materiale. Il modello riesce a colmare il divario tra la natura discreta degli atomi e il comportamento continuo degli oggetti solidi. Dimostra che, con gli strumenti matematici corretti, è possibile simulare il processo caotico e violento della rottura senza lasciare che la griglia del computer detti l'esito. Lo studio conclude che questo approccio offre un modo robusto per prevedere come i materiali falliscono in condizioni estreme, fornendo una via più chiara per progettare strutture più sicure e resilienti.