基于二阶锥规划理论的有限元强度折减法及应用

针对岩土体稳定性问题,常用的方法有极限平衡法和有限元强度折减法等。传统的有限元强度折减法通常需设置很大的最大允许非线性迭代次数(如200或500),计算耗时严重,此外,采用的平衡迭代和应力积分算法可能导致岩土体塑性区计算不够准确,进而影响稳定性分析结果。提出一种二阶锥规划有限元强度折减法,该方法基于Hellinger–Reissner混合变分原理和有限元法,将岩土体弹塑性问题构造成基于有限元框架的二阶锥规划问题,结合强度折减技术来分析岩土体稳定性。将该新方法应用于平面应变岩土体稳定性分析,结果表明:与传统的有限元强度折减法相比,新方法结果可靠,但其计算效率更高,所获得的塑性区更加平滑。