Peridynamic modeling of the crack velocity dependence via an incubation time fracture criterion

Questo studio utilizza un approccio peridinamico con un criterio di frattura basato su un tempo di incubazione per modellare gli esperimenti di Ravi-Chandar e Knauss su Homalite-100, rivelando che le variazioni del fattore di intensità di sforzo in modalità I a velocità di cricca costanti e l'inizio della micro-ramificazione a velocità più elevate offrono nuove prospettive sulla natura della dipendenza della velocità di cricca nella frattura dinamica.

L'essenziale

Immagina di osservare una crepa che si propaga attraverso un pezzo di plastica fragile, come un foglio di Homalite-100. Nei vecchi tempi della fisica, gli scienziati pensavano che, se si conosceva la velocità di propagazione della crepa, fosse possibile calcolare esattamente quanto "stress" (o pressione) la stesse spingendo in avanti. Era come pensare che, se un'auto viaggia a 100 km/h, il motore deve produrre esattamente 100 cavalli vapore. Semplice, vero?

Ma esperimenti negli anni '80 dimostrarono che non era così. A volte, la crepa si muoveva alla stessa identica velocità, ma la pressione che la spingeva era radicalmente diversa. Era come se due auto viaggiassero entrambe a 100 km/h, ma una avesse un minuscolo motore e l'altra un razzo booster. Gli scienziati erano perplessi: Perché la stessa velocità ha "spinte" diverse?

Questo articolo è una storia da detective in cui gli autori utilizzano un nuovo tipo di simulazione al computer per risolvere questo mistero.

Lo Strumento del Detective: La Peridinamica

La maggior parte dei modelli al computer delle crepe è come una catena di domino. Se un domino cade, spinge il successivo. Ma se manca un domino (una crepa), la catena si spezza e la matematica si blocca.

Gli autori hanno utilizzato un metodo chiamato Peridinamica. Immagina questo non come una catena, ma come uno sciame di api. Ogni ape può comunicare con ogni altra ape entro una certa distanza, anche se c'è un vuoto nel mezzo. Se un'ape vola via (si forma una crepa), le altre api smettono semplicemente di parlarle, ma il resto dello sciame continua a muoversi perfettamente. Questo permette al computer di gestire la rottura e la fessurazione senza confondersi.

L'Ingrediente Segreto: Il "Tempo di Incubazione"

La vera svolta in questo articolo è come hanno deciso quando una crepa dovrebbe effettivamente rompersi.

Nel vecchio modo, se la pressione diventava abbastanza alta, il materiale si rompeva istantaneamente. Ma gli autori hanno utilizzato una regola chiamata Criterio del Tempo di Incubazione.

Immagina di cercare di spezzare un rametto secco. Non lo tiri e si spezza istantaneamente. Lo tiri, lo mantieni lì per un istante mentre le fibre si allungano e si indeboliscono, e poi si spezza. Quel frazione di secondo è il "tempo di incubazione".

Gli autori hanno programmato il loro sciame al computer per ricordare gli ultimi microsecondi di pressione. Il materiale si rompe solo se la pressione media su quel breve periodo di "incubazione" è sufficientemente alta. Questo tiene conto del fatto che i materiali hanno bisogno di un piccolissimo lasso di tempo per "decidere" di rompersi.

Cosa Hanno Scoperto

Hanno eseguito simulazioni di lastre di plastica che venivano tirate, proprio come negli esperimenti reali. Ecco cosa hanno scoperto:

1. Il Puzzle Velocità vs Pressione: Proprio come negli esperimenti reali, il loro computer ha mostrato che per la stessa velocità di propagazione della crepa, la pressione (Fattore di Intensità delle Tensioni) non era un singolo numero. Era un intervallo. A volte era bassa, a volte alta.
2. L'Effetto "Micro-ramificazione": Quando la crepa si muoveva lentamente, procedeva dritta. Ma quando accelerava (oltre i 400 metri al secondo), iniziava a diventare irrequieta. Iniziava a germogliare minuscole crepe laterali microscopiche, come un ramo d'albero che si divide in rametti.
   • L'Analogia: Immagina un corridore che scatta. A un trotto costante, corre in linea retta. Ma quando scatta alla massima velocità, inizia a oscillare e zigzagare leggermente per mantenere l'equilibrio.
   • Il Risultato: Questi piccoli "oscillamenti" (micro-ramificazioni) facevano saltare la lettura della pressione su e giù in modo selvaggio. Questo spiegava perché la pressione non era unica per una data velocità; la crepa stava cambiando fisicamente forma leggermente mentre correva.

La Conclusione

L'articolo conclude che il motivo per cui vediamo diversi valori di pressione per la stessa velocità di crepa è che la crepa non è una linea liscia e perfetta. È una cosa caotica e vivente che fluttua.

• A velocità più basse: La crepa è stabile e la pressione è relativamente costante.
• A velocità più elevate: La crepa inizia a "micro-ramificarsi" (germogliare minuscole crepe laterali). Questo caos fa rimbalzare la pressione, creando la dispersione osservata negli esperimenti.

Utilizzando questo "sciame di api" (Peridinamica) combinato con il "periodo di attesa" (Tempo di Incubazione), gli autori hanno ricreato con successo la relazione disordinata e non unica tra velocità di propagazione della crepa e pressione che gli esperimenti reali avevano mostrato per decenni. Hanno dimostrato che il "rumore" nei dati non è un errore; è una caratteristica fisica reale di come si comportano le crepe che si muovono velocemente.

Sommario tecnico

Sintesi Tecnica: Modellazione Peridinamica della Dipendenza della Velocità di Cricca tramite un Criterio di Frattura a Tempo di Incubazione

Enunciato del Problema
Questo studio affronta una questione centrale irrisolta nella meccanica della frattura dinamica: la dipendenza del fattore di intensità degli sforzi istantaneo in Modo-I (SIF, $K_I$) dalla velocità di propagazione della cricca ($\dot{a}$). La Meccanica della Frattura Lineare Elastica (LEFM) classica e i criteri dinamici basati su Freund e Rosakis suggeriscono spesso una relazione univoca tra $K_I$ e $\dot{a}$. Tuttavia, i dati sperimentali di Ravi-Chandar e Knauss su lastre di polimero Homalite-100 fragile dimostrano una dipendenza non univoca, in cui una singola velocità di cricca corrisponde a un intervallo di valori SIF. Inoltre, a velocità più elevate, si verifica la micro-ramificazione, portando a una significativa dispersione nelle misurazioni SIF. La sfida risiede nel riprodurre numericamente questa non univocità e le relative fluttuazioni senza fare affidamento su proprietà materiali dipendenti dalla velocità aggiustate empiricamente.

Metodologia
Gli autori impiegano un quadro peridinamico non locale per simulare la propagazione dinamica della cricca nell'Homalite-100, replicando le condizioni sperimentali di Ravi-Chandar e Knauss. L'innovazione metodologica fondamentale è l'integrazione del Criterio di Frattura a Tempo di Incubazione (ITFC) nella formulazione peridinamica.

• Quadro Peridinamico: Le simulazioni utilizzano un modello peridinamico basato su stati ordinari implementato nel software open-source Peridigm. Il materiale è discretizzato in punti che interagiscono all'interno di un orizzonte sferico ($\delta$), consentendo la gestione naturale delle discontinuità come le cricche senza operatori differenziali spaziali.
• Criterio di Frattura a Tempo di Incubazione (ITFC): A differenza dei criteri classici che assumono un cedimento istantaneo, l'ITFC postula che la frattura richieda un tempo finito ($\tau$) affinché si verifichino processi preparatori (ad esempio, microcriccatura). Il criterio valuta lo sforzo mediato nel tempo durante il periodo di incubazione $\tau$ e una lunghezza caratteristica $d$.
  • La lunghezza caratteristica $d$ è derivata dal criterio quasi-statico di Irwin per garantire la coerenza con il SIF statico critico ($K_{Ic}$).
  • L'implementazione peridinamica adatta il criterio di frattura a Sforzo Remoto (RS). Il cedimento del legame è determinato integrando lo sforzo normale del legame durante il tempo di incubazione $\tau$ (Eq. 11).
  • Per gestire la complessità computazionale e focalizzarsi sulla micro-ramificazione piuttosto che sulla piena macro-ramificazione, il criterio è valutato solo su piani paralleli ($0^\circ$) e perpendicolari ($90^\circ$) all'asse principale della cricca.
• Configurazione della Simulazione: Il modello simula un provino di 300 mm $\times$ 500 mm con una pre-intaccatura di 50 mm. Il carico dinamico è applicato tramite un impulso volumetrico trapezoidale con un tempo di salita di 25 $\mu$s. Il SIF dinamico è calcolato utilizzando un contorno rettangolare mobile per l'integrale J dinamico, assicurando che il percorso di integrazione racchiuda la punta della cricca principale escludendo le micro-ramificazioni in sviluppo.

Risultati Chiave
Le simulazioni numeriche riproducono con successo i fenomeni osservati sperimentalmente riguardanti la relazione $K_I$-$\dot{a}$:

1. Dipendenza Non Univoca $K_I$-$\dot{a}$: I risultati confermano che per una velocità di propagazione della cricca quasi costante, il SIF istantaneo mostra fluttuazioni significative. Di conseguenza, il modello prevede un intervallo di valori SIF per una data velocità, rispecchiando la non univocità osservata negli esperimenti.
2. Stabilizzazione della Velocità e Transitori: Mentre gli esperimenti mostrano che la velocità della cricca si stabilizza quasi istantaneamente, le simulazioni esibiscono una fase iniziale di accelerazione della durata di circa 30–40 $\mu$s prima di raggiungere uno stato stazionario. Durante questa fase transitoria, la storia del SIF mostra variazioni pronunciate.
3. Micro-ramificazione e Dispersione: A tassi di carico più elevati (che risultano in velocità > 400 m/s), le simulazioni mostrano l'insorgenza della micro-ramificazione. Questo fenomeno correla con un marcato aumento dell'ampiezza delle oscillazioni SIF. La dispersione numerica nei valori SIF ad alte velocità (0,53–1,25 MPa$\sqrt{m}$) si allinea ragionevolmente bene con gli intervalli sperimentali (0,45–1,5 MPa$\sqrt{m}$).
4. Validazione del Meccanismo: Lo studio dimostra che la dispersione osservata nel SIF è fisicamente giustificata dalla dinamica della punta della cricca, specificamente la fase di accelerazione e l'insorgenza della micro-ramificazione, piuttosto che essere un artefatto di errore di misura.

Significato e Affermazioni
Il documento afferma che l'implementazione peridinamica dell'ITFC fornisce un approccio prezioso per valutare la tenacità alla frattura dinamica e modellare la propagazione della cricca su un'ampia gamma di tassi di carico. Il significato principale risiede nella capacità del modello di catturare la sensibilità al tasso della frattura fragile e la non univocità della relazione $K_I$-$\dot{a}$ senza introdurre parametri materiali arbitrari dipendenti dal tasso.

Gli autori notano modestamente che, sebbene il modello catturi qualitativamente le tendenze sperimentali, rimangono discrepanze quantitative, in particolare riguardo al tempismo della stabilizzazione della velocità e alla magnitudine precisa delle fluttuazioni SIF. Attribuiscono la previsione riuscita della dipendenza non univoca all'inclusione del tempo di incubazione, che tiene conto della storia temporale degli sforzi, e alla natura non locale della peridinamica, che gestisce naturalmente i campi singolari e le discontinuità associate alla frattura dinamica. Lo studio conclude che l'interazione tra la fase di accelerazione della cricca, le fluttuazioni SIF e l'insorgenza della micro-ramificazione sono i meccanismi governanti alla base del comportamento $K_I$-$\dot{a}$ osservato sperimentalmente.