The fracture resistance of elastic networks increases with the density of defects like a random walk

Technische Samenvatting: De Breukweerstand van Elastische Netwerken Neemt Toe met de Dichtheid van Defecten zoals een Random Walk

Probleemstelling
Gestoorde veernetwerken dienen als een fundamenteel model voor het bestuderen van breuk in diverse materialen, variërend van keramiek en polymeernetwerken tot mechanische metamaterialen. Een centrale grootheid in de breukmechanica is de schijnbare breukenergie, $G_c$, die de weerstand tegen de voortplanting van een bestaande dominante scheur meet. Hoewel vaststaat dat wanorde $G_c$ kan verhogen door scheurgroei te stoppen via lokale inhomogeniteiten, blijft de kwantitatieve afhankelijkheid van deze taaiheid van de mate van wanorde (specifiek de defectdichtheid) slecht begrepen. Vorige studies hebben aangetoond dat afwezigheid van een vooraf bestaande scheur, defecten over het algemeen nucleatie bevorderen en de breukweerstand verlagen. Echter, in systemen met een vooraf aangebrachte scheur, is de relatie complexer: bij lage defectdichtheden kan toenemende wanorde de weerstand tegen scheurgroei versterken, terwijl bij hoge dichtheden het materiaal verzwakt. De specifieke mechanismen die microscopische defectstatistieken koppelen aan macroscopische breukenergie in het regime van vooraf bestaande scheuren vereisen verdere theoretische en kwantitatieve verduidelijking.

Methodologie
De auteurs maken gebruik van een tweedimensionaal driehoekig netwerk van lineaire veren onder quasi-statische belasting om het effect van willekeurig geplaatste ontbrekende bindingen op scheurgroei te onderzoeken.

• Modelstelsel: Het netwerk bestaat uit veren met rustlengte $\ell_r$ en stijfheid $k$. Een veen breekt abrupt wanneer zijn rek een drempelwaarde $\varepsilon_{max} = 1$ overschrijdt.
• Invoering van Defecten: Een fractie $\nu$ van de bindingen wordt willekeurig verwijderd uit een anders perfect rooster.
• Simulatieprotocol: Het stelsel wordt belast via uniforme verticale verplaatsingsstappen. Het netwerk wordt tot statisch evenwicht gebracht met het FIRE-minimalisatiealgoritme. Veren die de rek-drempel overschrijden, worden sequentieel verwijderd, beginnend bij de meest uitgerekte, volgens een specifiek protocol om scheurgroei te simuleren.
• Belangrijkste Metingen: De schijnbare breukenergie $G_c(a)$ wordt gedefinieerd als de kritische elastische energie-afgifte die nodig is om de scheur voort te planten van een initiële lengte $L_c$ naar $L_c + a$. De studie richt zich op de initiële energie $G_c(0)$ en de faalenergie $G_c(a_f)$ bij een specifieke scheurvordering $a_f$.
• Theoretisch Kader: De auteurs mappingen de elastische heterogeniteit veroorzaakt door ontbrekende bindingen af op een equivalent lokaal breukenergielandschap, $\Gamma_{loc}(a)$. Zij passen gevestigde theorieën voor planaire scheurstop toe, ervan uitgaande dat de schijnbare breukenergie op elk punt wordt bepaald door de maximale lokale breukenergie die langs het scheurpad tot dat punt is aangetroffen.

Belangrijkste Resultaten

1. Taaiheid door Defecten: Voor een vaste concentratie ontbrekende bindingen $\nu$ neemt de schijnbare breukenergie $G_c(a)$ toe met de scheurvordering $a$. In tegenstelling tot perfecte netwerken die direct falen bij het breken van de eerste veer, vertonen gestoorde netwerken intermitterende scheurgroei gekenmerkt door afwisselende instabiele groei en stopzetting.
2. Schaling met Defectdichtheid: De standaardafwijking van de fluctuaties in de equivalente lokale breukenergie $\Gamma_{loc}$ neemt toe met de wortel van de fractie ontbrekende bindingen, $\sqrt{\nu}$. Bijgevolg schalen ook de gemiddelde faalbreukenergie $\langle G_c(a_f) \rangle$ met $\sqrt{\nu}$.
3. Random Walk-mechanisme: Deze $\sqrt{\nu}$-schaling wordt verklaard door de verstoringen veroorzaakt door individuele ontbrekende bindingen te modelleren als onafhankelijke toevalsvariabelen. De totale verstoring van de lokale breukenergie is de superpositie van deze individuele effecten, analoog aan een random walk. De variantie van de totale verstoring groeit lineair met het aantal ontbrekende bindingen, wat leidt tot een standaardafwijking evenredig met $\sqrt{\nu}$.
4. Exponentiële Staart en Logaritmische Groei: De waarschijnlijkheidsdichtheidsfunctie (PDF) van de lokale breukenergie $\Gamma_{loc}$ vertoont een exponentiële staart in plaats van een puur Gaussische. Deze exponentiële staart leidt tot een Gumbel-extreemwaardeverdeling voor de schijnbare breukenergie. Als gevolg daarvan neemt de mediaan van $G_c(a)$ logaritmisch toe met de scheurvordering $a$ ($G_c \propto \ln a$).
5. Dominantie van Scheurstop: De studie toont aan dat scheurstop door lokale fluctuaties in $\Gamma_{loc}$ het dominante taaiheidsverhogende mechanisme is in het onderzochte regime. Hoewel scheurpaden licht kunnen afwijken van rechte lijnen (ruwheid), voorspelt het model voor planaire scheurstop het waargenomen gedrag nauwkeurig, wat suggereert dat padruwheid een ondergeschikte rol speelt ten opzichte van de door lokale heterogeniteit veroorzaakte stopzetting.

Betekenis en Beweringen
Het artikel beweert een kwantitatief verband te leggen tussen microstructurele wanorde en macroscopische breukenergie. Door aan te tonen dat de fluctuaties in lokale breukenergie schalen als $\sqrt{\nu}$ vanwege de random-walk-achtige superpositie van verstoringen van individuele defecten, bieden de auteurs een theoretisch kader voor het voorspellen van de taaiheid van gestoorde materialen. De identificatie van de exponentiële staart in het lokale energielandschap verklaart de logaritmische toename van breukenergie met scheurvordering, een gedrag dat verschilt van systemen zonder vooraf bestaande scheuren.

De auteurs merken op dat deze bevindingen naar verwachting van toepassing zijn op andere roosterstructuren en continuüm-systemen waarbij ontbrekende bindingen overeenkomen met microholtes of microscheurtjes, mits de defectdichtheid laag blijft en breuk plaatsvindt langs een bijna rechte weg. Zij geven expliciet beperkingen aan, waaronder de aanname van planaire scheurgroei en de onafhankelijkheid van bindingverstoringen, en erkennen dat deze aannames kunnen bezwijken bij zeer hoge defectdichtheden of in driedimensionale systemen waar middeling van de scheurfront plaatsvindt. Het werk beoogt de weg te effenen voor kwantitatieve voorspellingen van breukenergie in een breed scala aan gestoorde materialen door microscopische defectstatistieken te koppelen aan macroscopische taaiheid.