Modélisation de la rupture progressive dans les versants rocheux raides par la méthode combinée éléments finis–éléments discrets

Pourquoi les parois montagneuses instables comptent

Des parois abruptes et rocheuses peuvent sembler solides, mais, dans certaines conditions, elles peuvent céder soudainement et envoyer des millions de tonnes de roche dévaler la pente. Pour les populations en aval de grands barrages ou les usagers des routes de montagne, comprendre quand ces versants peuvent lâcher est une question de sécurité, de planification et de coûts. Cette étude présente une nouvelle façon de « soumettre à l’essai » numériquement des pentes rocheuses raides afin que les ingénieurs voient non seulement si une pente est dangereuse, mais aussi comment et quand elle peut se fragmenter.

Faiblesses cachées dans des falaises imposantes

Les versants de montagne sont rarement constitués d’un bloc de roche homogène et parfait. Des alternances de grès et d’ardoise, des fractures anciennes et des zones altérées forment un patchwork de bandes plus ou moins résistantes. Dans l’ouest de la Chine, notamment le long des gorges profondément encaissées du haut cours du Lancang, des ouvrages hydroélectriques sont implantés directement dans de telles parois raides. Bien que ces versants puissent se déformer lentement sous leur propre poids, ils peuvent céder de manière catastrophique si des fissures se relient le long de couches vulnérables. Les outils d’ingénierie traditionnels estiment une marge de sécurité globale, mais ne rendent pas complètement compte de la manière dont de petites fissures commencent à la base de la pente, progressent vers le haut et finissent par libérer une masse rocheuse en mouvement.

Limites des outils de prédiction anciens

Les ingénieurs utilisent depuis longtemps trois principales familles de méthodes informatiques pour évaluer la stabilité des pentes. La première, la méthode des équilibres limites, équilibre forces motrices et résistantes le long d’une surface de glissement supposée ; elle est rapide mais oblige à deviner à l’avance la forme de cette surface et ne décrit que l’instant de la rupture. La seconde, la méthode des éléments finis, suit la déformation de la roche solide, mais peine à représenter l’apparition soudaine et la croissance des fissures. La troisième, la méthode des éléments discrets, modélise la pente comme un assemblage de blocs séparés qui peuvent se déplacer et entrer en collision, mais elle ne peut pas décrire la phase antérieure où la roche se comporte encore comme une masse continue. Aucune de ces méthodes, prise isolément, ne peut suivre de façon fluide une pente depuis un état apparemment intact jusqu’à des blocs éparpillés au fond de la vallée.

Une pente numérique qui se rompt de façon réaliste

Les auteurs combinent les points forts de ces approches dans un cadre appelé méthode éléments finis–éléments discrets (FDEM) et l’associent à une stratégie d’« augmentation de la gravité ». Dans leur modèle de pente virtuel, la roche est représentée par de nombreux petits éléments solides reliés entre eux par des joints invisibles. À mesure que les contraintes simulées augmentent, ces liaisons peuvent s’affaiblir et se rompre, transformant la roche continue en blocs séparés qui glissent et entrent en collision. Plutôt que de supposer une surface de rupture, le modèle augmente lentement l’attraction gravitationnelle effective jusqu’à ce que la pente manifeste un saut soudain de déplacement et d’énergie cinétique. Le niveau de gravité à ce moment-là fournit un facteur de sécurité, tandis que le schéma évolutif des fissures et des mouvements montre précisément comment la rupture se déroule.

Tester la méthode dans une gorge réelle

Pour vérifier si cette pente numérique se comporte comme dans la réalité, l’équipe a modélisé un versant rocheux de 796 mètres de haut près d’une centrale hydroélectrique sur le haut Lancang. Ils ont construit une coupe simplifiée incluant différentes zones d’altération et une faille clé, puis utilisé des propriétés mécaniques mesurées pour paramétrer le modèle. À mesure que la gravité était progressivement amplifiée, la simulation a reproduit une séquence réaliste : d’abord, de très petites fissures se sont formées près du pied de la pente dans la roche la plus altérée ; ces fissures se sont reliées en une bande continue faible ; ensuite une grande masse rocheuse a glissé le long de cette bande, s’est fragmentée en blocs et s’est déposée en une large accumulation dans la vallée. Les distances parcourues, les profondeurs des fissurations et la largeur finale des débris ont tous été en bon accord avec les observations de terrain, avec des écarts généralement inférieurs à 20 %.

Comparaison avec les méthodes standard

Les chercheurs ont comparé leur nouveau cadre avec les méthodes des équilibres limites et des éléments discrets largement utilisées en calculant la stabilité apparente de la même pente dans chaque cas. Les trois approches ont fourni des facteurs de sécurité similaires et se sont accordées sur la position générale de la principale surface de glissement le long de la couche fortement altérée et faible. La différence clé est que la nouvelle approche FDEM–gravité a fait plus que confirmer si la pente était marginalement stable. Elle a également révélé le moment du premier mouvement, le rôle de la faille dans le déclenchement d’effondrements secondaires provenant d’altitudes supérieures, et la manière dont les blocs interagissent et se fragmentent lors de leur descente — des détails cruciaux pour concevoir des renforcements et planifier des systèmes de surveillance.

Ce que cela signifie pour des projets montagnards plus sûrs

L’étude conclut que ce cadre de modélisation combiné peut suivre de manière fiable l’histoire complète d’un glissement de terrain, de la première microfissure au tas final de décombres. Pour le versant du Lancang, le facteur de sécurité calculé suggère que la paroi rocheuse est seulement tout juste stable dans les conditions actuelles, ce qui implique que des soutiens supplémentaires et une excavation prudente sont nécessaires. Plus généralement, la méthode offre aux ingénieurs un moyen de localiser où des dommages précoces sont susceptibles d’apparaître, où placer des ancrages et des instruments de surveillance, et comment des déclencheurs futurs comme les séismes ou de fortes pluies pourraient modifier la situation. Bien que les simulations actuelles soient bidimensionnelles et coûteuses en calcul, leur extension en trois dimensions et l’ajout de déclencheurs réalistes pourraient en faire un outil puissant pour la gestion de la sécurité à long terme dans les régions montagneuses raides.

Citation: 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

Mots-clés: glissements de terrain, stabilité des pentes rocheuses, modélisation numérique, barrages hydrauliques, rupture gravitationnelle