基于Mohr-Coulomb准则和二阶锥规划技术的轴对称自适应下限有限元法

摘要

轴对称Mohr-Coulomb准则屈服面的角点问题导致其在数值计算中存在困难,如何高效处理该屈服准则一直是极限分析下限有限元法的重要内容。首先,引入完全塑性假定将轴对称Mohr-Coulomb准则转化为1组不等式约束和3个线性等式约束;然后,将不等式约束直接转化为二阶锥约束,避免了对角点进行光滑近似处理;最后,将基于Mohr-Coulomb准则的轴对称极限分析下限有限元计算模型转化为具有较高计算效率的二阶锥规划数学优化模型。极限分析下限有限元法采用的线性应力单元难以精确模拟破坏区域的应力变化,单元的分布形式对计算精度存在较大影响。因此提出一种基于单元应力的网格自适应加密策略,通过判断单元内节点应力接近屈服的程度,自动识别破坏区域待加密的单元,实现对破坏区域应力分布的精确模拟,进而能够以较少单元获得高精度下限解。通过分析圆形基础承载力及竖向锚板极限抗拔力等典型轴对称岩土工程稳定性问题,表明了所提方法具有较高计算效率及计算精度,具有一定的理论价值和应用前景。