Undrained stability of a trapdoor

Abstract

This paper examines the stability of a purely cohesive layer resting on a trapdoor. Although this problem has been studied for a number of years and is of considerable interest for analysing the stability of tunnels and mine workings, its exact collapse load is unknown. Rigorous bounds on the load needed to support the trapdoor against active failure are derived using two numerical techniques which employ finite elements in conjunction with the limit theorems of classical plasticity. Both of the methods assume a perfectly plastic soil model with a Tresca yield criterion and lead to large linear programming problems. The solution to the lower bound linear programming problem defines a statically admissible stress field whereas the solution to the upper bound linear programming problem defines a kinematically admissible velocity field. The upper and lower bounds derived from the numerical methods typically bracket the exact collapse load to within ten per cent over the range of trapdoor geometries considered, and are substantially tighter than existing analytical bounds. They are shown to agree well with the collapse loads that are predicted from a traditional dis-placment type of finite element analysis. L'article examine la stabilité d'une couche purement cohérente reposant sur une trappe. Bien qu'on étudie depuis quelques années ce problème de grande importance pour analyser la stabilité des travaux miniers, la charge de rupture reste encore inconnue. Des limites strictes sur le chargement nécessaire pour la rupture active au droit de la trappe contre la rupture active sont derivées à l'aide de deux techniques numériques qui emploient des éléments finis conjointement avec les théorèmes limites de la plasticité classique. Les deux méthodes admettent un modèle de sol parfaitement plastique avec un critère d'ecoulement de Tresca et conduisent à des problemes considérables en ce qui concerne la préparation de programmes linéaires. La solution du problème de la préparation des programmes linéaires à la limite inférieure définit un champ de contrainte statiquement admissible, tandis que la solution de ce problème à la limite supérieure définit un champ de vitesse cinématiquement admissible. De façon typique, les limites supérieure et inférieure dérivées à partir des méthodes numériques se situent au voisinage du chargement de rupture précis à moins de dix pour cent près pour la gamme considérée de géometries de trappe et elles sont essentiellement plus précises que les limites analytiques à present acceptées. L'article démontre qu'elles s'accordent bien avec les chargements de rupture prédits à partir d'une analyse à elements finis traditionnelle basée sur le déplacement.