完全隐式返回映射算法对岩土地基问题的求解.pdfVIP

第 30 卷第 8 期 Vol.30 No.8 工 程 力 学 2013 年 8 月 Aug. 2013 ENGINEERING MECHANICS 83 文章编号：1000-4750(2013)08-0083-07 完全隐式返回映射算法对岩土地基问题的求解 王军祥，姜谙男 (大连海事大学道路与桥梁工程研究所，辽宁，大连 116026) 摘 要：针对岩土工程中的复杂力学问题，在弹塑性力学理论框架和非线性有限元理论基础上，采用非关联等向 硬化 Drucker-Prager 模型的完全隐式积分算法—返回映射算法(Return Mapping Algorithm)编制了有限元求解程序。 该算法可以避免预测应力漂移屈服面的现象，对准静态变形条件下的本构方程可以获得准确解，在迭代中使用 Newton-Raphson 法获得近似平方的收敛速率，具有较高的精确性和稳定性。对岩土工程中的地基问题进行求解， 计算得出位移、应力等结果，模拟了塑性区随载荷步增加的演化过程，对地基极限承载力进行了解析解和数值解 的对比。结果表明了算法的优越性、程序的正确性和实用性。 关键词：岩土力学；有限元；非关联等向硬化；Drucker-Prager 模型；完全隐式积分算法；返回映射算法 中图分类号：TU452 文献标志码：A doi: 10.6052/j.issn.1000-4750.2012.04.0253 FULLY IMPLICIT RETURN MAPPING ALGORITHM FOR SOLVING THE PROBLEMS OF GEOTECHNICAL FOUNDATION WANG Jun-xiang , JIANG An-nan (Institute of Road and Bridge Engineering, Dalian Maritime University, Dalian, Liaoning 116026, China) Abstract: To solve complicated mechanics problems in geotechnical engineering, the finite element program is compiled with fully implicit integration algorithm which is return mapping algorithm of non-associative isotropic hardening Drucker-Prager criterion based on the elastic-plastic mechanics theory frame and the nonlinear finite element theory. Return mapping algorithm can avoid the drift phenomenon of the trial stress, and can achieve the accurate solution of the constitutive equation on the condition of the quasi-static deformation, a quadratic convergence rate when using the Newton-Raphson iteration scheme, higher accuracy and stability. To solve the problems of foundations in geotechnical engineering, displacements and stresses are calculated, and the evolution process of plastic zones is simulated. The ultimate bearing capacity of the analytical s