Niet-ordentelijk staat-gebaseerd peridynamica-model voor scheurvoortplanting in gesteente: een gecombineerde stress-energie breukmethode

Waarom het breken van gesteente ertoe doet

Van het aanleggen van tunneltjes onder steden tot het ondergronds opslaan van kooldioxide: veel moderne projecten zijn afhankelijk van hoe gesteente barst en breekt. Het in real time volgen van scheurgroei in een massief stuk gesteente is echter buitengewoon moeilijk en duur. Dit artikel introduceert een nieuwe computergebaseerde methode om te simuleren hoe scheuren beginnen, groeien, afbuigen en zich verbinden in gesteente, vooral onder de ingewikkelde combinatie van druk- en schuifkrachten die in echte technische omstandigheden optreden.

Een nieuwe manier om naar scheuren te kijken

De meeste traditionele computermodellen behandelen gesteente alsof het uit een continu rooster van punten bestaat die alleen hun directe buren 'voelen'. Dat werkt goed totdat er een scheur verschijnt, want een scheur is in wezen een plotselinge onderbreking van die continuïteit. De methode die in deze studie wordt gebruikt, peridynamica genoemd, vertrekt van een ander beeld: elk klein stukje materiaal wisselt rechtstreeks interacties uit met veel andere stukken binnen een bepaalde afstand, verbonden door onzichtbare verbindingen. Wanneer het gesteente wordt belast, rekken deze verbindingen uit, worden ze samengedrukt of schuiven ze; als ze te ver worden belast, breken ze, waardoor scheuren op natuurlijke wijze ontstaan zonder extra regels of hernetwerken.

Verbetering van verborgen tekortkomingen in eerdere modellen

Hoewel dit op verbindingen gebaseerde gezichtspunt krachtig is, leden eerdere geavanceerde versies aan een subtiele numerieke fout. Vanwege de manier waarop de vervorming werd gemiddeld, liet het model soms zogenaamde "nul-energietoestanden" toe — bewegingen die vrijwel geen energie kosten en fysiek onzin zijn. Deze uitten zich als schijnbare gerafelde verplaatsingen en konden de voorspellingen over waar scheuren zouden ontstaan aantasten. De auteurs lossen dit euvel op door elke verbinding een eigen zorgvuldig geconstrueerde maat voor lokale vervorming toe te kennen, opgebouwd uit een gemiddelde van de twee eindpunten en gecorrigeerd zodat de oorspronkelijke en vervormde posities van elke verbinding strikt compatibel zijn. Deze beschrijving op verbindingsniveau herstelt een consistente lokale voorstelling van rek en schuif en snijdt de niet-fysische oscillaties scherp weg.

Het model leren schuif- en trekbelastingen te onderscheiden

Scheuren in gesteente zijn niet allemaal hetzelfde. Bij zuivere trek neigen scheuren recht naar open te gaan, maar bij compressie gecombineerd met zijwaartse schuifbewegingen buigen ze af, vertakken ze en vormen ze complexe patronen van trek-, schuif- en gemengde scheuren. Gangbare engineeringcriteria zoals de Mohr–Coulomb-regel waren niet altijd volledig betrouwbaar wanneer ze direct in peridynamica werden toegepast, vooral omdat ze de intermediaire spanningscomponent grotendeels negeren. In dit werk verwerken de auteurs een verfijnder "drievoudig schuifenergie"-criterium op verbindingsniveau. In plaats van alleen de maximale en minimale spanningen te vergelijken, meet deze benadering de schuifenergie op drie mogelijke interne glijvlakken en houdt rekening met de invloed van de middelste spanning. Een verbinding faalt als óf de trekspanning ervan de treksterkte van het gesteente overschrijdt, óf de geaccumuleerde schuifenergie een drempel overschrijdt die gekoppeld is aan de cohesie van het gesteente. Hierdoor kan de simulatie in grotere overeenstemming met laboratoriumwaarnemingen onderscheid maken tussen openscheur- en schuiffalen.

De methode aan de praktijk toetsen

Om aan te tonen dat het verbeterde model niet alleen wiskundig elegant maar ook praktisch bruikbaar is, toetsen de auteurs het aan verschillende referentieproblemen. Ze trekken eerst een plaat met een cirkelvormig gat en vergelijken de verplaatsingen met standaard eindige-elementresultaten, waarbij ze vrijwel identieke gladde velden vinden zonder de eerdere nul-energie artefacten. Vervolgens simuleren ze een balk met twee randdefecten onder een vierpunts buigproef, een klassiek experiment in gemengde scheurvorming. De voorspelde gebogen scheurpaden, piekbelasting en de algehele belasting–verplaatsingscurve komen goed overeen met zowel experimenten als andere hoogwaardige numerieke studies. Daarna worden halfronde buigproeven met schuine inkepingen geladen zodat de breukmodus geleidelijk verschuift van zuiver openen naar gemengd openen–schuiven. Het model reproduceert de waargenomen initiatiehoeken en eindpaden van scheuren over een reeks inkepingshoeken. Tenslotte pakt het meer realistische gesteenteblokken aan met één of twee bestaande defecten onder compressie en schuifbelasting. De simulaties vangen bekende patronen goed, zoals vleugelscheuren die wegschieten vanaf defecttips, conjugate schuifbanden en complexe verbindingen over gesteentebruggen, opnieuw in overeenstemming met laboratoriumfoto's uit eerder werk.

Wat dit betekent voor geotechniek

Samengevat laat de studie zien dat een zorgvuldige beschrijving van vervorming en falen op de schaal van individuele verbindingen een groot verschil maakt voor het voorspellen van hoe echte gesteenten breken. Door niet-fysische numerieke modi te elimineren en een op schuifenergie gebaseerd faalcriterium te gebruiken dat de volledige driedimensionale spanningsstaat respecteert, kan het model volgen waar en hoe scheuren beginnen, hoe ze afbuigen bij veranderende spanningen en hoe afzonderlijke scheuren uiteindelijk verbinden. Hoewel het nog rekenkundig veeleisend is en uitgaat van geïdealiseerde gesteenseigenschappen, biedt dit verbeterde peridynamica-raamwerk een betrouwbaarder numeriek "laboratorium" om scheurevolutie in brosse materialen te onderzoeken, met duidelijke potentiële voordelen voor veiliger ontwerp in mijnbouw, ondergrondse bouw en geotechnische energiesystemen.

Bronvermelding: Gong, B., Song, Y., Zhang, L. et al. Non-ordinary state-based peridynamics model for rock crack propagation: a combined stress-energy fracture method. Sci Rep 16, 11386 (2026). https://doi.org/10.1038/s41598-026-40833-8

Trefwoorden: modellering van breuken in gesteente, peridynamica, voortplanting van schuifscheuren, numerale simulatie, bros materiaal