基于六面体覆盖的三维数值流形方法的理论探讨与应用-PolyUIRA.PDF

第23 卷 第 10 期 岩石力学与工程学报 23( 10) ：1745～1754 2004 年5 月 Chinese Journal of Rock Mechanics and Engineering May ，2004 基于六面体覆盖的三维数值流形方法的 理论探讨与应用 1 1，2 郑榕明 张勇慧 (1 香港理工大学土木与结构工程系 香港) (2 中国科学院武汉岩土力学研究所 武汉 430071) 摘要 介绍了基于六面体单元的三维数值流形方法的基本理论，并在六面体单元中采用了有限元中C8 型等参单元 的形函数作为流形元中的覆盖函数。推导了相应平衡方程的刚度、初始应力、点荷载、体荷载、惯性力、固定点、 法向接触、切向接触和摩擦力子矩阵。在二维接触中，切向应力和摩擦力的方向很容易沿着接触边得到，但在三 维问题中要在接触面上确定该方向则非常困难，这也是所有计算方法在遇到三维接触时的难题。引入了向量理论 和迭代方法，避开了相应困难的干扰。最后，给出了一个简单的地下失稳算例，证明了该三维方法的正确性。 关键词 岩石力学，三维数值流形方法，六面体，覆盖函数，三维接触 分类号 TU 458.4 文献标识码 A 文章编号 1000-6915(2004)10-1745-10 APPLICATION AND DISCUSSION OF THREE DIMENSIONAL NUMERICAL MANIFOLD METHOD BASED ON HEXAHEDRON ELEMENT 1 1，2 Cheng Rongming ，Zhang Yonghui 1 ( Department of Civil and Structural Engineering ，Hong Kong Polytechnic University ，Hong Kong ，China) 2 ( Institute of Rock and Soil Mechanics ，CAS， Wuhan 430071 China) Abstract In this paper ，three-dimensional numerical manifold method (NMM) based on hexahedron element is proposed and the shape functions of C8 isoparametric element are used as the cover weight functions. All sub-matrices of equilibrium equations are derived according to this cover weight function. In 2D contact ，the direction of shear stiffness and friction force is easily defined along the contact edge. For the corresponding 3D contact problem，it is difficult to determine similar relations in simple way.