王海芝， 刘海洋， 韩振华*

（1. 北京市地质研究所， 北京 100120； 2. 中国科学院 地质与地球物理研究所， 北京 100029）

0 引言

我国多山， 滑坡灾害频繁发生， 尤其在每年的雨季， 滑坡灾害严重威胁危险区内的生命财产安全。 随着我国山区城镇化发展、 交通建设、 能源开发、 旅游开发等国民经济建设的全面展开， 潜在的滑坡灾害越来越多， 所造成的经济损失和人员伤亡也逐渐加大， 相应的稳定性分析对防灾减灾来说变得尤为重要。

边坡稳定性分析中最常规的是极限平衡分析法[1]， 其特点是方便快捷， 很多工程单位采用该方法来计算、 设计。 在工程实践中， 可根据边坡滑动面的形态来选择相应的极限平衡法。目前常用的极限平衡法有瑞典条分法、 Bishop法、 Janbu 法、 Sarma 法、 Morgenstern-Price 法和剩余推力法等[2-6]。 然而， 极限平衡分析法不能解决分析边坡应力和应变的问题。 后来出现了有限元分析边坡的稳定性， 有限元在解决小变形方面有其优越性， 但在解决大变形方面较为不便。 强度折减法由Zienkiewicz[7]于1975 年提出， Ugai[8]、 Duncan[9]、 郑颖人[10-11]、 赵尚毅[12]、徐平[13]等学者在有限元强度折减法应用方面做了大量的工作， 丰富了边坡稳定性分析的内容。

数值计算方法和计算机技术的飞速发展为边坡稳定性分析提供了有利的工具， 尤其是有限元、 有限差分、 离散元等数值模拟方法的广泛应用， 使深入研究滑坡岩土体内部的应力分布状态， 以及对坡体内外各点的变形模拟成为可能， 为复杂的滑坡稳定性分析开创了一条新的途径。……