基于水平集坐标及二维节理动态扩展过程无网格法模拟.docVIP

基于水平集坐标的二维节理动态扩展 庄晓莹1,2，黄润秋3，朱合华1,2 (1. 同济大学 岩土及地下工程教育部重点实验室, 上海 200092；2. 同济大学 地下建筑与工程系, 上海 200092； 3. 成都理工大学 地质灾害防治与地质环境保护国家重点实验室, 成都） 摘要：本文发展了模拟岩体节理张开闭合力学行为的无网格法近似函数构造方法，采用水平集坐标描述二维节理几何信息，沿任意几何形状节理两侧布置跳跃结点。在节理张开时引入断裂线，采用衍射准则构造无网格近似差值函数当节理闭合时，发展了类似于Goodman单元的无网格函数构造方法，考虑节理面上下两侧的相对滑动和法向接触。与不同之处在于本文构造方法基于材料点，对于几何形状复杂的节理更灵活。最后，应用本文方法对二维节理受拉、节理受压的静态和动态破裂扩展过程做了分析，并通过对照解验证了本文方法的正确性。 关键词：岩石力学；扩展无网格法；节理闭合； 中图分类号：O 319.56 文献标识码：A 文章编号：1000–6915(2005)01–0001–02 2D JOINT PROPAGTION IN ROCK USING THE MESHLESS METHODS AND LEVEL SETS ZHUANG Xiao-ying1,2 HUANG Run-qiu3 ZHU He-hua1,2 (1. Key Laboratory of Geotechnical Engineering of China Ministry of Education, Tongji University, Shanghai 200092, China; 2. Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China; 3. State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu, China) Abstract：In this paper, the meshless method, element-free Galerkin method is used to model 2D crack propagation in jointed rock mass where the level set method is used to describe the crack geometry. To capture the displacement jump across the joint when the join is open, the diffraction rule is used to modify the weight function, and hence construction of the discontinuous meshless approximation. The joint closure is modeled by a new meshless approximation, which is similar to the well known Goodman element in the family of finite elements to model the sliding and contact between joint surfaces. Different from Goodman element, the construction of stiffness matrix is based on material points instead of nodes and thereby is more flexible to deal with joint of complex geometry. Test examples are tested to model the crack propagation in rock samples subjected to compression and tension showing the correctness of the method. Key words：rock mechanics; propagation; meshless methods; joint closure; level sets; contact and friction; 1 引 言 岩体稳定性分析一直是岩土工程