Dynamic brittle fracture using Lip-field approach in an explicit dynamics context

Sintesi Tecnica: Frattura Fragile Dinamica mediante l'approccio Lip-field in un contesto di Dinamica Esplicita

Problematica
Questo articolo affronta la simulazione numerica della frattura fragile dinamica in materiali soggetti a carichi ad alto tasso di deformazione, come impatti o esplosioni. Mentre la meccanica della frattura classica (ad esempio, il modello di Griffith) fatica a gestire l'innesco e la ramificazione delle crepe, sono stati sviluppati approcci moderni come i Modelli a Zona Coesiva (CZM) e i modelli di danno diffuso. Tuttavia, i modelli di danno diffuso locale soffrono di una sensibilità patologica alla mesh, richiedendo una regolarizzazione per introdurre una scala di lunghezza e garantire risultati indipendenti dalla mesh. Le tecniche di regolarizzazione esistenti, come i metodi Phase-Field e Thick Level Set (TLS), comportano spesso funzionali di energia non locali che possono essere computazionalmente costosi. La sfida specifica affrontata qui è l'estensione dell'approccio Lip-field — un metodo variazionale che introduce scale di lunghezza tramite vincoli di Lipschitz sul campo di danno — dai suoi precedenti utilizzi in casi quasi-statici e dinamici monodimensionali alla frattura fragile bidimensionale all'interno di un framework di integrazione temporale esplicita.

Metodologia
Gli autori propongono una formulazione variazionale per la meccanica del danno dinamico basata sulla minimizzazione di un potenziale incrementale. Il nucleo della metodologia prevede:

1. Formulazione Variazionale: Il problema meccanico è definito dalla minimizzazione di un potenziale incrementale $F$ che include l'energia cinetica, l'energia libera e il potenziale di dissipazione. Le variabili di stato sono lo spostamento ($u$) e una variabile scalare di danno ($d$). Il modello costitutivo utilizza materiali standard generalizzati con specifiche funzioni di degradazione dell'energia ($g(d)$) e di addolcimento ($h(d)$).
2. Regolarizzazione Lip-field: A differenza dei modelli Phase-Field che modificano il funzionale di energia con un termine di gradiente, l'approccio Lip-field mantiene il funzionale di energia locale (non regolarizzato). Invece, impone vincoli di Lipschitz sul campo di danno, richiedendo che il gradiente del danno sia limitato da un valore critico ($1/l$). Ciò restringe il campo di danno a uno specifico spazio di Lipschitz $\mathcal{L}$, introducendo efficacemente una scala di lunghezza $l$ senza alterare la densità di energia locale.
3. Schema Staggered Esplicito: Il documento impiega uno schema di integrazione temporale esplicita di Newmark (metodo delle differenze centrali). Viene utilizzato un approccio staggered dove i campi di spostamento e di danno vengono risolti sequenzialmente. Per un dato passo temporale, lo spostamento viene aggiornato per primo, seguito dall'aggiornamento del danno.
4. Solutore di Danno Efficiente: Un importante contributo algoritmico è l'uso della stima dei limiti (derivata da lavori precedenti) per ottimizzare il calcolo del danno. Il processo prevede:
   • Risoluzione di un problema di minimizzazione puramente locale per trovare un campo di danno locale ($d_{loc}$).
   • Proiezione di $d_{loc}$ in limiti superiori e inferiori basati sul vincolo di Lipschitz.
   • Identificazione delle regioni in cui tali limiti differiscono.
   • Risoluzione del problema di minimizzazione non locale (vincolato da Lipschitz) solo nella piccola porzione del dominio in cui i limiti non sono uguali. Nelle regioni in cui i limiti coincidono, la soluzione locale viene accettata. Ciò riduce significativamente il costo computazionale rispetto alla risoluzione del problema non locale su l'intero dominio.
5. Discretizzazione: La discretizzazione spaziale utilizza elementi triangolari lineari per il campo di spostamento. Una "lip-mesh" separata (costruita dai centroidi degli elementi) viene utilizzata per memorizzare e vincolare il campo di danno, garantendo che il numero di vincoli di Lipschitz sia minimizzato.

Contributi Chiave

• Estensione alla Dinamica 2D: Il documento estende con successo l'approccio Lip-field, precedentemente validato in frammentazione dinamica 1D e casi quasi-statici, ai problemi di frattura fragile bidimensionale.
• Integrazione con la Dinamica Esplicita: Descrive un'implementazione specifica dell'approccio Lip-field utilizzando uno schema di integrazione temporale esplicita, che è computazionalmente efficiente per problemi di impatto ad alta velocità.
• Efficienza Computazionale tramite i Limiti: Gli autori dimostrano che l'utilizzo della stima dei limiti permette di restringere la minimizzazione non locale a piccoli sottodomini, rendendo l'approccio Lip-field computazionalmente competitivo per le simulazioni dinamiche.
• Distinzione Variazionale: Il lavoro chiarisce la distinzione tra approccio Lip-field e Phase-field: il Lip-field mantiene un funzionale di energia locale con vincoli non locali, mentre il Phase-field utilizza un funzionale di energia non locale.

Risultati della Simulazione
Il modello proposto è stato validato attraverso esempi numerici bidimensionali:

1. Single Edge Notched Tension (SENT): Il modello ha simulato con successo la propagazione e la ramificazione della crepa dinamica. I risultati hanno mostrato la convergenza della mesh nel potenziale e nell'energia dissipata. Il fenomeno della ramificazione della crepa è stato catturato senza criteri ad hoc, e le velocità della punta della crepa sono rimaste generalmente al di sotto della velocità dell'onda di Rayleigh, eccetto per picchi transitori durante gli eventi di ramificazione. I risultati sono comparabili alle simulazioni Phase-field e TLS trovate in letteratura.
2. Influenza dell'Intensità del Carico: Variando la trazione applicata, il modello ha dimostrato la capacità di catturare molteplici rami di crepa, con la ramificazione che avviene più vicino all'intaglio sotto carichi più elevati.
3. Test Kalthoff-Winkler: Il modello ha simulato l'innesco e la propagazione della crepa sotto carico di taglio dinamico.
   • A basse velocità di impatto ($33$ m/s), ha previsto un angolo di propagazione della crepa di Modo-I di $\approx 62.4^\circ$, coerente con le osservazioni sperimentali di $\approx 70^\circ$.
   • Ad alte velocità di impatto ($100$ m/s), ha catturato la successione di ramificazioni della crepa e molteplici crepe.
   • Limitazione: Gli autori osservano che, sebbene gli esperimenti mostrino una transizione verso la rottura duttile (frattura per taglio a $\approx -10^\circ$) ad alte velocità, l'attuale modello elastico-danno non cattura tale comportamento duttile, poiché manca della plasticità.

Significatività e Rivendicazioni
Il documento sostiene che l'approccio Lip-field offre un robusto framework variazionale per la frattura fragile che preserva la rilevanza matematica del problema introducendo al contempo le necessarie scale di lunghezza. La sua principale significatività risiede nella sua efficienza computazionale ottenuta attraverso il solutore basato sui limiti, che evita di risolvere equazioni non locali globalmente. Gli autori affermano che l'approccio gestisce efficacementamente fenomeni di frattura complessi come l'innesco, la ramificazione e la coalescenza senza criteri ad hoc.

Il documento conclude modestamente che, sebbene l'approccio Lip-field sia efficace per la frattura fragile e la frammentazione, esso è limitato al comportamento elastico-danno. Riconosce che la simulazione della rottura duttile ad alti tassi di deformazione richiederebbe un'estensione ai modelli di danno elastico-plastico, il che è fuori dall'ambito del presente lavoro. Il metodo è presentato come un'alternativa valida ai metodi Phase-field e TLS, particolarmente laddove l'efficienza computazionale nella dinamica esplicita è critica.