基于两级光滑边域混合单元的弹塑性二阶锥规划方法

弹塑性增量分析的数学规划方法是分析岩土工程变形和强度问题的有效途径之一,在处理非光滑屈服面、接触和多屈服面等复杂问题时具有独特的优势。为进一步简化计算和克服体积锁定问题,在广义Hellinger-Reissner(GHR)变分原理的基础上,提出一种新型混合常应力-两级光滑边域的三节点三角形单元,在关联流动法则和Mohr-Coulomb屈服准则条件下,将弹塑性增量问题转化为标准的二阶锥规划问题,并将黏-摩擦接触条件转化为锥约束条件引入到弹塑性增量分析的二阶锥规划问题中,随后采用高效的原对偶内点算法对其进行求解。最后将新方法用于岩土工程中两类经典问题的弹塑性数值分析。结果表明:新方法在计算效率、收敛性和精度方面均优于传统六节点三角形混合单元。