Introduction :
Compréhension générale
Les méthodes d’analyse de la stabilité des talus se divisent en deux grandes familles : les
méthodes analytiques classiques, aussi appelées méthodes d’équilibre limite (MEL), et les
méthodes numériques, plus récentes.
Les méthodes analytiques sont simples et rapides, mais elles ne prennent pas toujours bien en
compte les effets complexes comme la présence d’eau, les séismes ou les discontinuités dans
le sol.
Les méthodes numériques, elles, reposent sur les équations de la mécanique des milieux
continus. Elles permettent de calculer la répartition des contraintes et des déformations
dans un talus et d’obtenir une vision plus réaliste de son comportement.
Aujourd’hui, la tendance est de combiner ces deux approches.
La partie numérique sert à calculer précisément les contraintes internes, tandis que la
méthode d’équilibre limite est utilisée pour déterminer la surface de rupture critique et le
facteur de sécurité global.
Cette combinaison permet d’obtenir une analyse plus complète et plus fiable de la stabilité
des ouvrages en terre.
Principales contributions de la recherche
De nombreux chercheurs ont travaillé à l’amélioration des méthodes d’analyse de la stabilité
des talus en combinant approches numériques et équilibre limite.
En 1980, Baker a proposé une méthode de programmation dynamique associée à la
méthode de Spencer, permettant de trouver simultanément la surface de rupture et le facteur
de sécurité minimal.
En 1981, Martins et ses collègues ont utilisé la programmation linéaire pour optimiser la
recherche de la surface de rupture, en l’appliquant à des milieux fracturés intégrés à un code
d’éléments finis.
En 1988, Takuo et al. ont développé un algorithme de programmation dynamique basé sur
les contraintes calculées par la méthode des éléments finis, sans imposer de forme de
surface de glissement.
En 1993, Thompson a montré que le couplage entre les méthodes numériques (éléments
finis ou différences finies) et la programmation dynamique permet d’obtenir un facteur de
sécurité plus réaliste et une meilleure localisation de la surface de rupture.
En 1996, Ennour et al. ont proposé d’utiliser la valeur maximale du cisaillement comme
indicateur d’instabilité et ont observé que, dans les sols cohérents à faible angle de frottement,
les ruptures se produisent en profondeur.
Les années 1996 à 1998 ont vu un grand progrès avec les logiciels ITASCA (UDEC, FLAC,
etc.), qui introduisent une méthode simple pour calculer le facteur de sécurité en réduisant
progressivement l’angle de frottement jusqu’à la rupture. Cette approche permet aussi de
repérer les zones plastifiées et les surfaces potentielles de glissement.
Toujours en 1996, Stanley a proposé une méthode combinant les éléments finis et l’équilibre
limite, prenant en compte l’interaction sol-structure et la non-linéarité des propriétés du sol. Il
a défini un facteur de sécurité global calculé à partir de la moyenne pondérée des facteurs
locaux de chaque élément du maillage.
Takuo et Jiang (1996) ont ensuite développé une méthode 3D sans surface prédéfinie,
utilisant la programmation dynamique et la méthode de Janbu simplifiée.
En 1998, Law et al. ont prolongé la méthode de Morgenstern et Price en exprimant les
équilibres sous forme d’équations intégrales résolues à l’aide d’un tableur Excel, rendant la
recherche du facteur de sécurité plus rapide et souple.
En 1999, Korini a combiné la méthode des éléments finis et celle de l’équilibre limite en
supposant une surface circulaire, tandis que Fry a utilisé la méthode des différences finies
pour calculer des facteurs de sécurité locaux et en déduire un facteur global moyen.
La même année, Goh (1999) a introduit les algorithmes génétiques pour déterminer
automatiquement la surface critique de rupture, ouvrant la voie à des approches
d’optimisation plus avancées.
📊 3. Synthèse et conclusion
L’ensemble de ces travaux montre que le couplage entre méthodes numériques et
méthodes d’équilibre limite constitue une voie efficace pour améliorer l’analyse de stabilité
des talus.
Les méthodes numériques permettent de calculer les contraintes, déformations et zones
plastifiées, tandis que les méthodes d’équilibre limite servent à déterminer le facteur de
sécurité et la surface critique de rupture.
Ce couplage offre ainsi une analyse plus réaliste du comportement des sols, notamment pour
les talus, pentes naturelles et ouvrages en terre.
La méthode proposée dans le mémoire s’appuie sur ce principe : elle calcule l’état de
contrainte numérique, détermine un facteur de sécurité local le long de la surface de
glissement, puis utilise la méthode de Bishop pour obtenir le facteur global de sécurité.
Cette approche cherche donc à améliorer les limites des méthodes classiques d’équilibre
limite et à mieux adapter l’analyse aux conditions réelles de stabilité des ouvrages
géotechniques.