压应力状态下多裂纹扩展过程数值模拟.pdfVIP

32卷 2期 岩 土 工 程 学 报 Vol.32 No.2 2010 年 .2 月 Chinese Journal of Geotechnical Engineering Feb. 2010 应力状态下多裂纹扩展过程数值模拟 周小平，杨海清，董 捷 (重庆大学土木工程学院，重庆400045) 摘 要：基于扩展有限元法，建立了压应力状态下多裂纹扩展过程数值模拟方法。该方法采用罚函数法考虑裂纹面的 摩擦，然后对刚度矩阵和荷载列阵进行修正，通过迭代计算可以得到压应力状态下裂纹的扩展路径及应力场。单元形 函数中包含有间断函数和裂尖富集化函数，因 裂纹的扩展路径与网格无关，从而避免了裂纹扩展计算过程中不断对 网格重新划分。由于考虑了裂纹闭合的影响， 数值模拟方法适用于岩土工程中常见的压剪应力状态多裂纹扩展。最 后，将 数值模拟结果与相似模型试验结果进行对比，分析表明与相似模型试验结果接近，验证了 数值模拟方法的 可行性。 关键词：压应力；多裂纹；裂纹扩展；数值模拟；罚函数法 中图分类号：TU452 文献标识码：A 文章编号：1000 –4548(2010)02 –0192 –06 作者简介：周小平(1970 – )，男，江西瑞金人，教授，博士生导师，从事岩石力学与工程研究方面的工作。E-mail: zhouxiaopinga@ 。 Numerical simulation of multiple-crack growth under compressive loads ZHOU Xiao-ping，YANG Hai-qing，DONG Jie (School of Civil Engineering, Chongqing University, Chongqing 400045, China) Abstract: A numerical simulation method for modeling the growth of multiple frictional cracks in linear elastic media within the framework of the extended finite element method is presented. The contact constraint on crack surface is embedded within the total stiffness matrix by the penalty method, in which the load vectors are modified too. The path of crack growth and stress field are determined through iterative computations. Because additional discontinuous functions and enriched tip elements are added in the displacement field, the geometry of cracks is independent of the finite element mesh. As a result, remeshing is not necessary to model crack propagation. Because the frictional contact of crack surface is taken into account, the present method is fit for modeling the growth of multiple cracks in geomaterials under compressive loads. F