Mohr-Coulomb准则在Mohr应力空间中具有最简形式,同时也因其非常可靠而在经典的极限分析或极限平衡法中得到了最广泛的应用.然而,应力空间中的Mohr-Coulomb屈服面是非光滑的,这给基于变形分析的弹塑性有限元法中的本构积分带来了巨大的麻烦.此外,在求解强度问题时,基于载荷控制法(LCM)的求解器很难将有限元模型带入极限平衡状态.针对这些问题,本研究给出了以下解决方案.首先,设计了适用于带非光滑屈服面的塑性本构积分算法GSPC.GSPC对任意大小的应变增量都收敛,数值特性远优于现有的返回-映射算法.还为弹塑性有限元分析定制了一个位移控制法(DCM)求解器,该DCM能将有限元模型带入极限平衡状态而不存在收敛性方面的问题,计算效率和鲁棒性都远优于现有的基于LCM的求解器.最后,结合强度折减法,建议了求解边坡安全系数的割线法,并给出了极限平衡状态下坡顶拉裂缝位置和深度的确定技术.