边坡失稳是一种严重的地质灾害,灾害一旦发生就会造成大量人员伤亡与巨大财产损失,因此边坡稳定分析对于土工结构设计与灾害防治至关重要。同时,一些学者也认为,土压力与地基承载力问题亦可采用边坡稳定分析思想进行分析。目前工程应用中的边坡稳定分析方法大都存在一些不足,如,刚体极限平衡法需要引入一定假定条件才能求解;极限分析法计算过程较为复杂;有限元强度折减法在某些情况下可能会存在违背材料自身性质等问题。然而,有限元极限平衡法可以避免上述问题,该方法理论在任何工况下土体结构的稳定性分析中都可取得较好的效果。基于此,本文在前人研究的基础上,详细地介绍了有限元极限平衡法的理论,进而开发出了相应的计算程序,并将该程序应用于土体学三大经典问题的分析,同时与传统理论得到的结果进行对比。本文主要完成了以下工作:(1)讨论了土体内一点的极限平衡状态与土体沿整个滑动面的极限平衡状态之间的关系,在此基础上定义了土体沿整体与局部滑动面的安全系数,并以此作为土体整体与局部稳定性评价标准。(2)采用粒子群算法与模式搜索法结合进行临界滑动面的搜素,提高了局部和全局搜索能力,通过生成局部微小滑动面的方式实现了局部破坏区的搜索,并基于Qt框架与Python语言开发了相应程序。(3)本文中基于有限元极限平衡法开发的程序与传统理论方法在边坡稳定分析中得到的结果进行对比,有限元极限平衡法计算得到的临界滑动面上土体剪应力与抗剪强度分布符合实际特征,其位置和安全系数与瑞典圆弧法、简化Bishop法、有限元强度折减法得到的结果一致性较好。在土压力问题中,有限元极限平衡法得到的局部破坏区与临界滑动面符合Rankine土压力理论。在地基承载力分析中,采用逐渐增加地基载荷的方式,得到了地基的渐进破坏过程,有限元极限平衡法得到的极限载荷与Prandtl理论计算的结果同样具有较好的一致性。