刘 芸，李 林

(1.同济大学土木工程学院地下建筑与工程系，上海200092;2.岩土及地下工程教育部重点实验室(同济大学)，上海200092)

0 引 言

随着我国基础建设投资力度的增大，边坡工程的稳定问题也日渐突出.每年我国发生的边坡失稳灾害不计其数，严重威胁人民的生命财产安全［1］.因此，为加强对边坡及其作用形式的认识，正确理解边坡渐进性失稳的机理具有重大的现实指导意义.

就目前而言，分析边坡稳定性的主要方法有两种，分别为传统的基于极限平衡法和有限元数值模拟法［2～4］.但前者的分析方法中，需要做出许多近似假设［5～6］，因此结果中往往得不到滑移体内土体的变形特性以及应力迁移规律，也得不到支护结构对土坡变形及稳定性的影响.近年来，计算机技术不断发展，基于有限元的强度折减系数法也开始受到重视.它可以灵活地结合强度折减法和三维弹塑性有限元分析法，在已知评判指标的前提下，通过不断地调整稳定性折减系数，求出土坡发生破坏前提下的最小稳定安全系数，进而分析判定边坡稳定性.而边坡稳定性分析的变分法以变分原理和极限平衡法为基础研究边坡的稳定性问题［7］.上述两种方法在模拟边坡渐进破坏过程以及确定边坡滑裂面位置方面具有不可比拟的优势.

本文以上海地区某均质土坡工程为背景，构建二维有限元数值模型，模拟土坡的破裂区域的渐进性破坏发展过程，从而根据塑性区……