基于罚函数偶应力理论的土体应变局部化研究

与 Cosserat 理论相比,偶应力理论在一定程度上可以降低数值框架的复杂度,已逐渐应用于岩土体应变局部化分析中。然而,一般的偶应力有限元法需要满足 C¹连续性,即单元内部和单元交界面上的应变都需要具有连续性。为了避免开发较为复杂的C¹型偶应力单元,在 Cosserat 连续体理论框架下,通过借助罚函数方法对 C¹连续性进行松弛来获得偶应力理论的逼近解,建立了基于罚函数的偶应力有限元方法 PCS-FEM。通过平面应变条件下的弹性圆孔应力集中问题对 PCS-FEM 方法的有效性进行了验证,并应用于土体应变局部化分析中。通过对Ottawa砂的平面应变试验进行数值模拟,发现 PCS-FEM 方法获得的应力−应变曲线及剪切带破坏形态与试验结果基本一致,且能够克服经典连续体理论病态的网格敏感性问题,保证应变局部化问题的适定性;通过对承受偏心荷载作用下的土坡应变局部化经典算例进行分析,发现 PCS-FEM 方法同样可以克服土坡应变软化阶段的网格敏感性问题,展现土体的渐进破坏过程。     

基于二阶锥规划有限元增量加载法的条形浅基础极限承载力分析

地基极限承载力分析是土力学研究中的一个经典课题。基于Hellinger-Reissner混合变分原理和有限元方法,将岩土体弹塑性问题构造成基于有限元框架的二阶锥规划（second-order cone programming,SOCP）问题,进而提出一种基于二阶锥规划理论的增量有限元法,即FEM-SOCP法。将岩土体弹塑性问题构造成二阶锥规划的数学优化问题,可以避免采用传统弹塑性计算中复杂的应力点积分等算法和屈服面棱角的平滑处理。此外,对于二阶锥规划问题,可以采用具有原始?对偶内点求解法的标准数学规划求解器MOSEK进行求解。将增量加载FEM-SOCP法应用于经典的基底粗糙的条形浅基础地基极限承载力分析中,分别考虑了关联和非关联塑性条件下的Mohr-Coulomb屈服准则。数值结果表明：所提出的增量加载FEM-SOCP法获得的地基承载力系数及地基承载力与传统FEM计算结果基本一致,而与常规有限元计算结果相比,基于增量加载的FEM-SOCP法所获得的屈服区更加平滑。