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

Deze studie maakt gebruik van een peridynamische aanpak met een incubatietijd-breukcriterium om de experimenten van Ravi-Chandar en Knauss op Homalite-100 te modelleren, waarbij wordt aangetoond dat variaties in de Mode-I spanningsintensiteitsfactor bij constante krasnelheden en het optreden van microvertakking bij hogere snelheden nieuwe inzichten bieden in de aard van de snelheidsafhankelijkheid van krasnelheid bij dynamische breuk.

De kern

Stel je voor dat je een scheur ziet racen door een stuk broos plastic, zoals een vel Homalite-100. In de oude dagen van de natuurkunde dachten wetenschappers dat als je wist hoe snel de scheur bewoog, je precies kon berekenen hoeveel "spanning" (of druk) erop werd uitgeoefend om haar vooruit te duwen. Het was als denken dat als een auto 96 km/u rijdt, de motor precies 100 paardenkracht moet leveren. Eenvoudig, toch?

Maar experimenten in de jaren 1980 toonden aan dat dit niet waar was. Soms bewoog de scheur met exact dezelfde snelheid, maar was de druk die haar voortduwde enorm verschillend. Het was alsof twee auto's allebei 96 km/u reden, maar de ene een piepklein motor had en de andere een raketmotor. Wetenschappers waren in de war: Waarom heeft dezelfde snelheid verschillende "duwkrachten"?

Dit artikel is een detectiveverhaal waarin de auteurs een nieuw soort computersimulatie gebruiken om dit mysterie op te lossen.

Het detectivegereedschap: Peridynamica

De meeste computermodellen van scheuren zijn als een rij dominostenen. Als één steen valt, duwt hij de volgende. Maar als een steen ontbreekt (een scheur), breekt de keten en blijft de wiskunde steken.

De auteurs gebruikten een methode genaamd Peridynamica. Denk hierbij niet aan een keten, maar aan een zwerm bijen. Elke bij kan met elke andere bij binnen een bepaalde afstand praten, zelfs als er een gat in het midden zit. Als een bij wegvliegt (een scheur vormt zich), stoppen de andere bijen gewoon met praten met hem, maar de rest van de zwerm beweegt perfect verder. Hierdoor kan de computer breken en scheuren afhandelen zonder in de war te raken.

Het geheime ingrediënt: De "Incubatietijd"

De echte doorbraak in dit artikel is hoe ze beslisten wanneer een scheur daadwerkelijk zou breken.

Op de oude manier brak het materiaal direct als de druk hoog genoeg werd. Maar de auteurs gebruikten een regel genaamd het Incubatietijdcriterium.

Stel je voor dat je een droge tak probeert te breken. Je trekt niet zomaar en hij breekt direct. Je trekt, houdt het een fractie van een seconde vast terwijl de vezels rekken en verzwakken, en dan breekt hij. Die fractie van een seconde is de "incubatietijd".

De auteurs programmeerden hun computervermenging om de druk van de laatste paar microseconden te onthouden. Het materiaal breekt alleen als de gemiddelde druk over die korte "incubatie"-periode hoog genoeg is. Dit houdt rekening met het feit dat materialen een klein beetje tijd nodig hebben om te "beslissen" om te breken.

Wat ze vonden

Ze voerden simulaties uit van de plastic platen die uit elkaar werden getrokken, net als de echte experimenten. Dit is wat ze ontdekten:

1. Het Snelheid versus Druk Raadsel: Net als bij de echte experimenten toonde hun computer aan dat voor dezelfde scheursnelheid de druk (Spanningsintensiteitsfactor) geen enkel getal was. Het was een bereik. Soms was het laag, soms hoog.
2. Het "Micro-vertakking"-effect: Wanneer de scheur langzaam bewoog, ging hij recht. Maar toen hij versnelde (boven de 400 meter per seconde), begon hij te trillen. Hij begon kleine, microscopische zij-scheurtjes te vertonen, alsof een tak van een boom zich splitst in twijgjes.
   • De Analogie: Stel je een hardloper voor die sprint. Bij een steady jog rennen ze in een rechte lijn. Maar wanneer ze sprinten op topsnelheid, beginnen ze te wiebelen en lichtjes te zigzaggen om het evenwicht te bewaren.
   • Het Resultaat: Deze kleine "wiebels" (micro-vertakkingen) zorgden ervoor dat de drukmeting wild omhoog en omlaag sprong. Dit verklaarde waarom de druk niet uniek was voor een gegeven snelheid; de scheur veranderde fysiek zijn vorm lichtjes terwijl hij racete.

De Conclusie

Het artikel concludeert dat de reden waarom we verschillende drukwaarden zien voor dezelfde scheursnelheid is dat de scheur geen gladde, perfecte lijn is. Het is een chaotisch, levend ding dat fluctueert.

• Bij lagere snelheden: De scheur is stabiel en de druk is relatief stabiel.
• Bij hogere snelheden: De scheur begint te "micro-vertakken" (kleine zij-scheurtjes laten groeien). Dit chaos zorgt ervoor dat de druk heen en weer springt, wat de spreiding veroorzaakt die in de experimenten werd gezien.

Door deze "zwerm bijen" (Peridynamica) te combineren met de "wachtperiode" (Incubatietijd), slaagden de auteurs erin de rommelige, niet-unieke relatie tussen scheursnelheid en druk die decennialang door werkelijke experimenten was aangetoond, succesvol te reconstrueren. Ze bewezen dat de "ruis" in de data geen fout is; het is een echt fysiek kenmerk van hoe snel bewegende scheuren zich gedragen.

Technische samenvatting

Technische Samenvatting: Peridynamische Modellering van de Snelheidsafhankelijkheid van Scheurgroei via een Incubatietijd-Bruchcriterium

Probleemstelling
Deze studie adresseert een centraal onopgelost probleem in de dynamische breukmechanica: de afhankelijkheid van de instantane Mode-I spanningsintensiteitsfactor (SIF, $K_I$) van de scheurgroeisnelheid ($\dot{a}$). Klassieke Lineair Elastische Breukmechanica (LEFM) en dynamische criteria gebaseerd op Freund en Rosakis suggereren vaak een unieke relatie tussen $K_I$ en $\dot{a}$. Echter, experimentele data van Ravi-Chandar en Knauss aan breekbare Homalite-100 polymeerplaten tonen een niet-unieke afhankelijkheid, waarbij één enkele scheurgroeisnelheid correspondeert met een bereik van SIF-waarden. Bovendien treedt bij hogere snelheden micro-vertakking op, wat leidt tot aanzienlijke spreiding in SIF-metingen. De uitdaging ligt in het numeriek reproduceren van deze niet-uniekheid en de bijbehorende fluctuaties zonder te vertrouwen op empirisch aangepaste, snelheidsafhankelijke materiaaleigenschappen.

Methodologie
De auteurs maken gebruik van een niet-lokaal peridynamisch kader om dynamische scheurgroei in Homalite-100 te simuleren, waarbij de experimentele omstandigheden van Ravi-Chandar en Knauss worden gerepliceerd. De kernmethodologische innovatie is de integratie van het Incubatietijd-Bruchcriterium (ITFC) in de peridynamische formulering.

• Peridynamisch Kader: De simulaties maken gebruik van een gewone staatsgebaseerd peridynamisch model, geïmplementeerd in de open-source software Peridigm. Het materiaal is gediskretiseerd in punten die interageren binnen een sferische horizon ($\delta$), wat een natuurlijke behandeling van discontinuïteiten zoals scheuren mogelijk maakt zonder ruimtelijke differentiaaloperatoren.
• Incubatietijd-Bruchcriterium (ITFC): In tegenstelling tot klassieke criteria die een onmiddellijk falen veronderstellen, stelt ITFC dat breuk een eindige tijd ($\tau$) vereist voor voorbereidende processen (bijv. microscheurvorming) om plaats te vinden. Het criterium evalueert de tijd-gemiddelde spanning over de incubatietijd $\tau$ en een karakteristieke lengte $d$.
  • De karakteristieke lengte $d$ is afgeleid uit het quasi-statische Irwin-criterium om consistentie te waarborgen met de kritische statische SIF ($K_{Ic}$).
  • De peridynamische implementatie past het Remote Stress (RS) breukcriterium aan. Bondfalen wordt bepaald door de normale bondspanning te integreren over de incubatietijd $\tau$ (Vergelijking 11).
  • Om de computationele complexiteit te beheersen en te focussen op micro-vertakking in plaats van volledige macro-vertakking, wordt het criterium alleen geëvalueerd op vlakken evenwijdig ($0^\circ$) en loodrecht ($90^\circ$) op de hoofdscheuras.
• Simulatieopstelling: Het model simuleert een specimen van 300 mm $\times$ 500 mm met een voorafgaande inkeping van 50 mm. Dynamische belasting wordt aangebracht via een trapeziumvormige volumetrische krachtpuls met een opstijgtijd van 25 $\mu$s. De dynamische SIF wordt berekend met behulp van een bewegend rechthoekig contour voor de dynamische J-integraal, waarbij wordt gewaarborgd dat het integratiepad de hoofdscheurtip omsluit en ontwikkelende micro-vertakkingen uitsluit.

Belangrijkste Resultaten
De numerieke simulaties reproduceren succesvol de experimenteel waargenomen fenomenen met betrekking tot de $K_I$-$\dot{a}$-relatie:

1. Niet-Unieke $K_I$-$\dot{a}$-Afhankelijkheid: De resultaten bevestigen dat voor een bijna constante scheurgroeisnelheid de instantane SIF aanzienlijke fluctuaties vertoont. Bijgevolg voorspelt het model een bereik van SIF-waarden voor een gegeven snelheid, wat de in experimenten waargenomen niet-uniekheid weerspiegelt.
2. Snelheidsstabilisatie en Transiënten: Waar experimenten aantonen dat de scheurgroeisnelheid bijna onmiddellijk stabiliseert, vertonen de simulaties een initiële acceleratiefase van ongeveer 30–40 $\mu$s voordat een steady-state wordt bereikt. Tijdens deze transiëntfase toont de SIF-historie uitgesproken variaties.
3. Micro-vertakking en Spreiding: Bij hogere belastingssnelheden (wat resulteert in snelheden > 400 m/s) vertonen de simulaties het begin van micro-vertakking. Dit fenomeen correleert met een duidelijke toename in de amplitude van SIF-oscillaties. De numerieke spreiding in SIF-waarden bij hoge snelheden (0,53–1,25 MPa$\sqrt{m}$) sluit redelijk goed aan bij experimentele bereiken (0,45–1,5 MPa$\sqrt{m}$).
4. Validatie van het Mechanisme: De studie toont aan dat de waargenomen spreiding in SIF fysisch gerechtvaardigd is door de dynamiek van de scheurtip, specifiek de acceleratiefase en het begin van micro-vertakking, en niet het gevolg is van een artefact van meetfouten.

Betekenis en Aanspraken
Het artikel stelt dat de peridynamische implementatie van het ITFC een waardevolle aanpak biedt voor het beoordelen van dynamische breuktaaiheid en het modelleren van scheurgroei over een breed scala aan belastingssnelheden. De primaire betekenis ligt in het vermogen van het model om de snelheidsgevoeligheid van breekbare breuk en de niet-uniekheid van de $K_I$-$\dot{a}$-relatie vast te leggen zonder willekeurige snelheidsafhankelijke materiaalparameters in te voeren.

De auteurs merken bescheiden op dat, hoewel het model de experimentele trends kwalitatief vastlegt, kwantitatieve discrepanties blijven bestaan, met name wat betreft het tijdstip van snelheidsstabilisatie en de precieze grootte van SIF-fluctuaties. Zij schrijven het succesvolle voorspellen van de niet-unieke afhankelijkheid toe aan de opname van de incubatietijd, die rekening houdt met de temporele geschiedenis van de spanning, en het niet-lokale karakter van peridynamica, dat op natuurlijke wijze omgaat met de singuliere velden en discontinuïteiten die geassocieerd zijn met dynamische breuk. De studie concludeert dat het samenspel tussen de scheuracceleratiefase, SIF-fluctuaties en het begin van micro-vertakking de regulerende mechanismen zijn achter het experimenteel waargenomen $K_I$-$\dot{a}$-gedrag.