Modellering van progressief falen in steile rotshellingen met de gecombineerde eindige‑elementen‑ en discrete‑elementenmethode

Waarom wankele bergwanden ertoe doen

Hoge, rotsachtige kloofwanden lijken solide, maar onder bepaalde omstandigheden kunnen ze plotseling bezwijken en miljoenen tonnen rots naar beneden laten donderen. Voor mensen stroomafwaarts van grote dammen of weggebruikers op bergwegen is het van levensbelang om te weten wanneer zulke hellingen kunnen falen; het gaat om veiligheid, planning en kosten. Deze studie onderzoekt een nieuwe manier om steile rotshellingen digitaal te "druktesten", zodat ingenieurs niet alleen kunnen vaststellen of een helling onveilig is, maar ook precies hoe en wanneer ze kan uiteenvallen.

Verborgen zwakheden in torenhoge kliffen

Berghellingen bestaan zelden uit uniform, ongeschonden gesteente. Lagen zandsteen en leisteen, oude breuken en verweerde zones vormen een lappendeken van sterke en zwakke banden. In westelijk China, vooral langs de diep uitgesneden kloven van de bovenste Lancangrivier, worden waterkrachtprojecten direct in zulke steile wanden gebouwd. Hoewel deze hellingen langzaam onder hun eigen gewicht vervormen, kunnen ze catastrofaal falen als scheuren zich verbinden langs kwetsbare lagen. Traditionele ingenieursinstrumenten kunnen een globale veiligheidsmarge inschatten, maar vatten niet volledig hoe kleine scheurtjes bij de voet van de helling ontstaan, omhoog groeien en uiteindelijk een bewegende massa rots loslaten.

Beperkingen van oudere voorspellingsmethoden

Ingenieurs gebruiken al lange tijd drie hoofdtypen numerieke methoden om de stabiliteit van hellingen te beoordelen. De limiet-evenwichtsmethode balanceert drijvende en remmende krachten langs een verondersteld glijvlak; ze is snel maar vereist dat het glijvlak van tevoren wordt verondersteld en beschrijft alleen het ogenblik van falen. De eindige‑elementenmethode volgt hoe massief gesteente vervormt, maar heeft moeite met het weergeven van het plotseling ontstaan en groeien van scheuren. De discrete‑elementenmethode stelt de helling voor als losse blokken die kunnen bewegen en botsen, maar kan de eerdere fase niet beschrijven waarin het gesteente nog als een continuum functioneert. Geen van deze methoden kan op zichzelf soepel een helling volgen van schijnbaar intact gesteente naar verspreide keien in het dal.

Een digitale helling die realistisch breekt

De auteurs combineren de sterke punten van deze benaderingen in een raamwerk dat de eindige–discrete elementenmethode (FDEM) wordt genoemd, en koppelen dit aan een "verhoging van de zwaartekracht"‑strategie. In hun virtuele helling wordt het gesteente voorgesteld als veel kleine massieve delen die langs onzichtbare voegen aan elkaar gebonden zijn. Naarmate de gesimuleerde spanningen toenemen, kunnen deze verbindingen verzwakken en breken, waardoor continu gesteente in afzonderlijke blokken verandert die schuiven en botsen. In plaats van een glijvlak te veronderstellen, verhoogt het model geleidelijk de effectieve zwaartekracht totdat de helling een plotselinge sprong in beweging en kinetische energie vertoont. Het zwaartekrachtsniveau op dat moment geeft een veiligheidsfactor, terwijl het evoluerende patroon van scheuren en beweging precies laat zien hoe het falen verloopt.

De methode testen in een echte kloof

Om te onderzoeken of deze digitale helling zich als in de echte wereld gedraagt, modelleerde het team een 796 meter hoge rotshelling naast een waterkrachtstation aan de bovenloop van de Lancangrivier. Ze bouwden een vereenvoudigde dwarsdoorsnede die verschillende verweerde zones en een belangrijke breuk omvatte, en stelden de modelparameters in op basis van gemeten gesteenteeigenschappen. Terwijl de zwaartekracht geleidelijk werd versterkt, reproduceerde de simulatie een realistische volgorde: eerst ontstonden kleine scheurtjes nabij de voet van de helling in het meest verweerde gesteente; deze scheurtjes voegden zich samen tot een continue zwakke band; vervolgens schoof een grote rotsmassa langs deze band, brak in blokken en kwam tot rust als een brede afzetting in het dal. De afgelegde afstanden, dieptes van scheurvorming en uiteindelijke puinbreedte kwamen nauw overeen met veldobservaties, met verschillen meestal onder de 20 procent.

Vergelijking met standaardmethoden

De onderzoekers vergeleken hun nieuwe raamwerk met de veelgebruikte limiet‑evenwichtsmethode en discrete‑elementenmethode door te berekenen hoe stabiel dezelfde helling in elk geval leek. Alle drie benaderingen gaven vergelijkbare veiligheidsfactoren en waren het eens over de algemene ligging van het belangrijkste glijvlak langs de zwakke, sterk verweerde laag. Het belangrijkste verschil was dat de nieuwe FDEM–zwaartekrachtbenadering meer deed dan bevestigen of de helling marginal stabiel was. Ze toonde ook het tijdstip van de eerste beweging, de rol van de breuk bij het uitlokken van secundaire instortingen hogerop, en hoe blokken interacteerden en fragmentiseerden tijdens hun afdaling—details die cruciaal zijn voor het ontwerpen van versterkingen en het plannen van monitoringsystemen.

Wat dit betekent voor veiligere bergprojecten

De studie concludeert dat dit gecombineerde modelleringsraamwerk betrouwbaar het hele levensverhaal van een aardverschuiving kan volgen, van de eerste haarfijne scheurtjes tot de uiteindelijke puinhoop. Voor de Lancang‑helling suggereert de berekende veiligheidsfactor dat de rotswand onder de huidige omstandigheden net stabiel is, wat impliceert dat extra ondersteuning en zorgvuldige ontgraving nodig zijn. Breder gezien biedt de methode ingenieurs een manier om te bepalen waar vroegtijdige schade waarschijnlijk begint, waar ankers en meetinstrumenten te plaatsen en hoe toekomstige triggers zoals aardbevingen of zware regenval het beeld kunnen veranderen. Hoewel de huidige simulaties tweedimensionaal en rekenkundig intensief zijn, zou het uitbreiden naar drie dimensies en het toevoegen van realistische triggers ze tot een krachtig onderdeel van het langdurig veiligheidsbeheer in steile berggebieden kunnen maken.

Bronvermelding: Xu, J., Deng, Z., Feng, Y. et al. Modeling progressive failure in steep rock slopes using the combined finite-discrete element method. Sci Rep 16, 11180 (2026). https://doi.org/10.1038/s41598-026-40966-w

Trefwoorden: aardverschuivingen, stabiliteit van rotshellingen, numerieke modellering, waterkrachtcentrales, gravitationeel falen