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

Wanneer een bros materiaal zoals glas of gehard staal breekt, knapt het niet simpelweg doormidden. In plaats daarvan raast een scheur door de vaste stof, waarbij deze vaak vertakt in een complex web van breuken die het object in vele stukken versplintert. Het precies voorspellen waar deze scheuren zullen reizen en hoe snel ze zullen bewegen, is een van de moeilijkste uitdagingen in de techniek. Als een model dat deze vernietiging simuleert gebrekkig is, kan het suggereren dat een brug op een veilige manier zal falen, terwijl deze in werkelijkheid catastrofaal instort. Decennialang hebben wetenschappers geprobeerd computermodellen te bouwen die de atomaire structuur van materialen nabootsen, waarbij ze deze behandelen als een rooster van minuscule punten verbonden door veren. Hoewel deze modellen krachtig zijn, lijden ze aan een vreemde tekortkoming: het pad dat de scheur aflegt, hangt vaak af van hoe het computerraster is geroteerd, in plaats van op de fysica van het materiaal zelf. Het is alsof een scheur in een echt stuk glas plotseling van richting zou veranderen simpelweg omdat de ingenieur het computerscherm heeft gedraaid.

Om dit probleem op te lossen, heeft Matías Braun, een onderzoeker aan de Universidad Nebrija in Spanje, een nieuwe manier ontwikkeld om dynamische breuk te simuleren. Zijn aanpak combineert een traditioneel rooster van punten met een "faseveld"-methode, een techniek die de scheur niet behandelt als een scherpe lijn die van punt naar punt springt, maar als een gladde, vage zone van schade die zich door het materiaal verspreidt. Door deze twee ideeën te mengen, creëerde Braun een model waarmee scheuren in elke gewenste richting kunnen reizen, ongeacht hoe het onderliggende computerraster is georiënteerd. Hij testte deze nieuwe methode op drie verschillende scenario's: een stalen plaat die geraakt wordt door een hogesnelheidsimpact, een ingesneden plaat die uit elkaar wordt getrokken tot hij knapt, en een dikke metalen cilinder die van binnenuit explodeert. In elk geval produceerde het nieuwe model resultaten die overeenkwamen met experimenten in de echte wereld en elimineerde het de kunstmatige bias veroorzaakt door het computerraster, wat een duidelijker en betrouwbaarder beeld geeft van hoe materialen werkelijk falen.

De kern van het probleem dat Braun aanpakte, ligt in de manier waarop computers de onzichtbare wereld van atomen verwerken. In zijn eerdere werk gebruikte hij een discreet model waarbij het materiaal werd gerepresenteerd door een driehoekig rooster, een patroon van stippen verbonden door lijnen. Wanneer het materiaal onder spanning kwam te staan, zouden de verbindingen tussen de stippen breken, wat een scheur simuleert. Echter, omdat de stippen in een vast patroon waren gerangschikt, werd de scheur gedwongen om de lijnen van het rooster te volgen. Als het rooster werd geroteerd, veranderde de scheur van pad, ook al was het materiaal zelf niet veranderd. Dit is een aanzienlijk probleem omdat een scheur in werkelijkheid het pad van de minste weerstand moet volgen, bepaald door de eigenschappen van het materiaal, en niet door de willekeurige hoek van een computermesh.

Brauns oplossing was het introduceren van een faseveld-regularisatie. In plaats van te eisen dat specifieke verbindingen tussen de stippen breken, kende hij een waarde toe aan elke stip die representatief is voor hoe beschadigd deze is. Deze waarde varieert van nul, wat betekent dat het materiaal perfect intact is, tot één, wat betekent dat het volledig gebroken is. De scheur wordt een regio waar deze waarde geleidelijk toeneemt, wat een gladde overgang creëert in plaats van een grillige lijn. Dit stelt de scheur in staat om in elke richting te buigen en te vertakken, waardoor het de rigide structuur van het rooster effectief negeert. Cruciaal is dat Braun dit systeem heeft gekalibreerd met behulp van twee fundamentele eigenschappen van het materiaal: de treksterkte, wat de kracht is die nodig is om een scheur te laten ontstaan, en de breukenergie, wat de energie is die nodig is om de scheur te laten groeien. Door deze simulatie te verankeren aan deze echte fysische constanten, zorgt het model ervoor dat de resultaten afhangen van het materiaal, en niet van de wiskunde.

Om te bewijzen dat deze nieuwe methode werkte, voerde Braun drie rigoureuze tests uit. De eerste was de Kalthoff–Winkler-test, een klassiek experiment waarbij een stalen plaat met een vooraf aangebrachte inkeping door een projectiel op hoge snelheid wordt geraakt. In dit scenario schiet een scheur vanuit de inkeping onder een specifieke hoek naar buiten. Eerdere modellen voorspelden vaak de verkeerde hoek afhankelijk van hoe het rooster was geroteerd. In Brauns simulaties reisde de scheur consequent onder een hoek van ongeveer 70 graden, wat de experimentele waarnemingen perfect mat. Dit resultaat bleef standhouden, of het rooster nu uitgelijnd was met de plaat of met 15 of 30 graden was geroteerd. De snelheid van de scheur bleef ook stabiel en bewoog zich op ongeveer 70 procent van de geluidssnelheid in het materiaal, onafhankelijk van de oriëntatie van het rooster.

De tweede test betrof een plaat met een enkele inkeping aan de rand onder trekspanning, een opstelling die wordt gebruikt om te bestuderen hoe scheuren langzaam groeien onder een constante trekkracht. Hier was het doel om te zien of de scheur in een rechte lijn van de inkeping naar de tegenoverliggende rand zou reizen. Traditionele roostermodellen dwongen de scheur vaak om te zigzaggen of af te wijken op basis van de uitlijning van het rooster. Brauns model produceerde echter een rechte scheur voor elke oriëntatie van het rooster. De kracht die nodig was om de plaat te breken, was ook consistent en varieerde met minder dan twee procent, ongeacht hoe het rooster werd geroteerd. Dit toonde aan dat het model zowel het pad als de sterkte van het materiaal nauwkeurig kon voorspellen zonder de kunstmatige invloed van de computermesh.

De laatste en meest complexe test was de fragmentatie van een dikke cilinder die onderhevig is aan een interne explosie. In dit scenario wordt geen enkele scheur door een inkeping gestuurd; in plaats daarvan moeten veel scheuren spontaan ontstaan en naar buiten racen, waardoor de cilinder in tientallen stukken uiteenspat. Dit is een chaotisch evenement waarbij de symmetrie van het rooster gemakkelijk het patroon van de fragmenten zou kunnen dicteren. Wanneer Braun de simulatie draaide met het rooster gecentreerd op de as van de cilinder, vormden de scheuren een symmetrisch patroon. Echter, wanneer hij het rooster iets verschoof zodat het niet langer perfect gecentreerd was, verdeelden de scheuren zich onregelmatig en veranderde het aantal fragmenten. Dit bevestigde dat het model gevoelig was voor de fysieke opstelling in plaats van alleen de symmetrie van het rooster. De simulatie voorspelde tussen de 12 en 19 fragmenten, een aantal dat binnen de range lag van andere geavanceerde methoden, maar lager was dan sommige andere, waarschijnlijk omdat de gladde schadezone in het model voorkwam dat scheuren die te dicht bij elkaar lagen, zich vormden.

Gedurende al deze tests hielden de onderzoekers de energie van het systeem zorgvuldig bij. Ze vergeleken de energie die verloren ging terwijl het materiaal brak met de energie die nodig was om de nieuwe scheuroppervlakken te creëren. In sommige gevallen was er een kleine discrepantie tussen de energie berekend door het schademodel en de energie die daadwerkelijk vrijkwam bij het breken van de verbindingen. Dit is een bekende beperking van de hybride aanpak, waarbij de drijvende kracht voor de scheur anders wordt berekend dan de energie van de gebroken verbindingen. Echter, de mechanische energiebalans bleef nauwkeurig tot binnen een halve procent, wat bevestigt dat de simulatie fysisch solide was. De resultaten suggereren dat, hoewel het model niet perfect is, het een significante stap voorwaarts is in het begrijpen van dynamische breuk.

De implicaties van dit werk reiken verder dan alleen betere computergraphics. Door de afhankelijkheid van de roosteroriëntatie te verwijderen, kunnen ingenieurs erop vertrouwen dat hun simulaties van crashtests, structurele defecten of materiaaltesten de ware fysica van het materiaal weerspiegelen. Het model overbrugt succesvol de kloof tussen de discrete aard van atomen en het continue gedrag van vaste objecten. Het toont aan dat het met de juiste wiskundige instrumenten mogelijk is om het chaotische, gewelddadige proces van breken te simuleren zonder dat het rooster van de computer de uitkomst dicteert. De studie concludeert dat deze aanpak een robuuste manier biedt om te voorspellen hoe materialen falen onder extreme omstandigheden, wat een duidelijk pad biedt voor het ontwerpen van veiligere en veerkrachtigere structuren.