加速对称刚度矩阵技术及其在岩土有限元分析中的应用

陈曦; 程勇刚

加速对称刚度矩阵技术及其在岩土有限元分析中的应用 English Version

陈曦¹,
程勇刚²

1. 北京交通大学土木建筑工程学院，北京 100044 ； 2. 武汉大学水资源与水电工程科学国家重点实验室，湖北武汉 430072

基金项目: 教育部留学回国人员科研启动基金资助项目（教外司留 [2010]1561 号）；中央高校基本科研业务费专项资金资助项目（ 2011JB-M070 ）

详细信息

作者简介:
陈曦 (1977 – ) ，男，博士，副教授，主要从事岩土工程数值计算方法等方面的研究。

中图分类号: TD853.34
计量
- 文章访问数: 1239
- HTML全文浏览量: 1
- PDF下载量: 441
出版历程
- 发布日期: 2011-08-14

Accelerated symmetric stiffness matrix technique and its applications in geotechnical finite element analysis

CHEN Xi¹,
CHENG Yong-gang²

1. School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China; 2. State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China

摘要

摘要: 非关联塑性岩土工程问题的有限元离散会产生非对称刚度矩阵，对于大规模岩土工程问题，大量非对称线性方程组的求解会显著增加内存需求和计算耗时。基于加速初始刚度法的思想，提出了加速对称刚度矩阵技术来求解具有非关联塑性模型的岩土工程问题。进一步评价和比较了两个基于最小二乘法的加速技术以及两个近似对称刚度矩阵，即弹性刚度矩阵 K e 和投影到塑性势面的具有关联塑性流的等效材料的刚度矩阵 K G 。通过一个二维孔洞扩张算例和一个二维边坡算例，评价和比较了这些加速刚度矩阵技术和 Newton-Raphson 迭代法。数值结果表明加速 K G 具有较好的计算性能，例如，加速 K G 技术只需要对称线性求解器，收敛速度较快，计算耗时有时会少于 Newton-Raphson 迭代法。
- 非关联塑性流 /
- 非对称刚度矩阵 /
- 加速对称刚度 /
- 初始刚度 /
- 最小二乘法 /
- 剪切强度折减法
Abstract: Finite element discretization of non-associated plastic geotechnical problems may result in non-symmetric stiffness matrices, and solving the nonsymmetric linear systems of equations arising from large-scale geotechnical problems may significantly increase the computer memory storage requirement and the computer runtime. Based on the accelerated initial stiffness matrix technique, the accelerated symmetric stiffness matrix techniques are proposed to solve the non-associated plastic geotechnical problems. Furthermore, two acceleration techniques based on the least-square minimization and two approximate symmetric stiffness matrices (i.e. the elastic stiffness matrix K _e and the stiffness matrix K _G obtained from an equivalent material with associated plastic flow by mapping the yield surface to the plastic potential surface) are assessed and compared. By using a 2-D cavity expansion example and a 2-D slope example, the accelerated symmetric stiffness matrix techniques and the Newton-Raphson iteration are evaluated and compared, and numerical results show that the accelerated K _G techniques possess better computational performances. For example, they only resort to symmetric linear solvers, and have faster convergence rates, and the consumed computer runtime may be less than that of the Newton-Raphson iteration.
- non-associated plastic flow /
- nonsymmetric stiffness matrix /
- accelerated symmetric stiffness /
- initial stiffness /
- least square method /
- shear strength reduction method

HTML全文

参考文献(16)

[1]	VERMEER P A, DE BORST R. Non-associated plasticity for soils, concrete and rock[J]. Heron, 1984, 29 (3): 1 – 64.
[2]	杨强, 陈新, 周维垣. 基于 D—P 准则的三维弹塑性有限元增量计算的有效算法 [J]. 岩土工程学报, 2002, 24 (1): 16 – 22. (YANG Qiang, CHEN Xin ， ZHOU Wei-yuan. A practical 3D elasto-plastic incremental method in FEM based on D-P yield criteria[J]. Chinese Journal of Geotechnical Engineering, 2002, 24 (1): 16 – 22. (in Chinese))
[3]	COLLIN I F. Associated and non-associated aspects of the constitutive laws for coupled elastic/plastic materials[J], International Journal of Geomechanics, 2002, 2 (2): 259 – 267.
[4]	熊文林. 非关联塑性切线刚度矩阵的对称表示 [J]. 应用数学和力学, 1986, 7 (11): 983 – 992. (XIONG Wen-lin. Symmetric formulation of tangential stiffnesses for non-associated plasticity[J]. Applied mathematics and mechanic, 1986, 7 (11): 983 – 992. (in Chinese))
[5]	ZIENKIEWICZ O C, VALLIAPPAN S, KING I P. Elasto-plastic solutions of engineering problems ‘initial stress’, finite element approach[J]. International Journal for Numerical Methods in Engineering, 1969, 1 (1): 75 – 100.
[6]	郑宏, 刘德富. 弹塑性矩阵 D _ep 的特性和有限元边坡稳定性分析中的极限状态标准 [J]. 岩石力学与工程学报, 2005, 24 (7): 1099 – 1105. (ZHENG Hong, LIU De-fu. Properties of elasto-plastic matrix D_ep and a criterion on limiting state of slope stability by FEM[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24 (7): 1099 – 1105. (in Chinese))
[7]	CHEN C. Efficient and reliable accelerated constant stiffness algorithms for the solution of non-linear problems[J]. International Journal for Numerical Methods in Engineering, 1992, 35 (3): 481 – 490.
[8]	THOMAS J N. An improved accelerated initial stress procedure for elasto-plastic finite element analysis[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1984, 8 (4): 359 – 379.
[9]	SLOAN S W, SHENG D, ABBO A J. Accelerated initial stiffness schemes for elastoplasticity[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2000, 24 (6): 579 – 599.
[10]	FAMIYESIN O O R. Robust symmetric formulations for nonassociated plasticity problems[J]. Journal of Engineering Mechanics, 1999, 125 (9): 1071 – 1080.
[11]	YU H S. Cavity expansion methods in geomechanics[M]. Dordrecht: Kluwer Academic, 2000.
[12]	赵尚毅, 郑颖人, 时卫民, 等. 用有限元强度折减法求边坡稳定安全系数 [J]. 岩土工程学报, 2002, 24 (3): 343 – 346. (ZHAO Shang-yi, ZHENG Ying-ren, SHI Wei-min, et al. Analysis on safety factor of slope by strength reduction FEM[J]. Chinese Journal of Geotechnical Engineering, 2002, 24 (3): 343 – 346. (in Chinese))
[13]	陈祖煜. 土质边坡稳定分析 — 原理·方法·程序 [M]. 北京：中国水利水电出版社， 2003. (CHEN Zu-yu. Soil slope stability analysis—theory, methods and programs[M]. Beijing: China Water Power Press, 2003. (in Chinese))
[14]	陈曦, 刘春杰. 有限元强度折减法中安全系数的搜索算法 [J]. 岩土工程学报, 2010, 32 (9): 1443 – 1447. (CHEN Xi, LIU Chun-jie. Search algorithms for safety factor in finite element shear strength reduction method[J]. Chinese Journal of Geotechnical Engineering, 2010, 32 (9): 1443 – 1447. (in Chinese))
[15]	MANZARI M T, NOUR M A. Significance of soil dilatancy in slope stability analysis[J]. Journal of Geotechnical and Geoenvironmental Engineering, ASCE, 2000, 126 (1): 75 – 80.
[16]	WANG Y J, YIN J H, LEE C F. The influence of a non-associated flow rule on the calculation of the factor of safety of soil slopes[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2001, 25 (13): 1351 – 1359.