第31卷 第12期 岩 土 工 程 学 报 Vol.31 No.122009年12月 Chinese Journal of Geotechnical Engineering Dec. 2009一种模拟节理岩体破坏的新方法：单元劈裂法张振南，陈永泉（上海大学土木工程系，上海 200072）摘 要：在有限单元法中，三结点三角单元在几何上具有独特的优点，即当一条裂纹贯穿单元（单元劈裂）时，总有一结点位于裂纹的一侧，而另外两个结点位于裂纹的另一侧。位于一侧的两个结点与位于另外一侧的一个结点可以潜在地构成两个接触点对，从而可以利用这两个接触点对来推导劈裂单元的刚度矩阵，以再现裂纹面之间的接触和摩擦效应。利用该三角单元的几何特点，提出了单元劈裂法，用于模拟节理或其它裂纹的扩展问题。由于劈裂单元刚度矩阵与原三角单元享有共同的结点，因而在模拟裂纹扩展过程中无需为设置节理单元而改变原网格划分方案，这为节理扩展的数值模拟提供了很大方便，提高了计算效率。在劈裂单元刚度矩阵推导过程中，没有考虑单元劈裂后所形成两个块体自身的弹塑性变形，所以该方法只是一种近似方法。对相应试验进行了数值模拟，结果表明该方法是有效的。关键词：单元劈裂法；节理扩展；裂纹；数值模拟；多维虚内键模型中图分类号：TU457 文献标识码：A 文章编号：1000–4548(2009)12–1858–08作者简介：张振南(1974–)，男，博士，主要从事岩石断裂破坏力学模型与数值模拟研究。E-mail:zhennanzhang@sina.com.cn。Novel numerical approach to jointed rock mass simulation:element partition methodZHANG Zhen-nan，CHEN Yong-quan(Department of Civil Engineering, Shanghai University, Shanghai 200072, China)Abstract: In the finite element method, the tri-node triangular element exhibits a geometrical characteristic, that is, when afracture runs through it, one node is always located at one side of the fracture and another two nodes located at the other side.The former node can potentially construct two contact pairs with the latter two nodes. Through the two contact pairs is thestiffness matrix of the partitioned element derived to represent the contact and friction effect between interfaces of joint orfracture. Based on the advantage of this geometrical characteristic of triangular element, an element partition method (EPM) isdeveloped to simulate the joint or fracture propagation of jointed rock mass. Due to the fact that the stiffness matrix of thepartitioned element shares the common nodes of the corresponding intact triangular element, it doesn’t need to modify theoriginal mesh configuration for setting up the joint element. The fracture propagates in the manner of successive elementpartition. To represent the newly generated and pre-existing fractures, it is just to displace the stiffness matrix of original intacttriangular element with that of the partitioned element. This makes the simulation of the fracture propagation highly convenientand efficient. However, the present method is only an approximate method for the elastic-plastic deformation of the two bodiesgenerated by the element partition. By simulating an experiment of fracture propagation and coalescence, good agreement isfound between the experimental and the simulated results, which suggests that the present method is validated and feasible.Key words: element partition method; joint propagation; fracture; numerical simulation; virtual multidimensional internal bondmodel0 引 言在应用有限单元法对节理岩体进行数值模拟过程中一个较为棘手的问题是节理扩展后网格处理问题及如何考虑新生节理面的接触和摩擦问题。对于已有节理裂隙，一般采用节理单元，如Goodman节理单元[1]以及后来发展起来的各种节理单元[2]或接触单元[3]等。这些方法可以很好地处理已有节理问题。根据已有节理的分布情况，可以不断地加密网格使计算达到任意精度。但是在荷载作用下，节理是不断变化的，───────基金项目：国家自然科学基金项目（50609013）收稿日期：2008–07–21
本文来源：https://www.2haoxitong.net/k/doc/abf2250c76c66137ee061901.html
