数值流形方法在有限元三维二十结点单元上实现.pdfVIP

第23 卷 第15 期 岩石力学与工程学报 23(15) ：2622～2627 2004 年8 月 Chinese Journal of Rock Mechanics and Engineering Aug .，2004 数值流形方法在有限元三维二十结点单元上的实现 彭自强 葛修润 ( 中国科学院武汉岩土力学研究所 武汉 430071) 摘要 在有限元三维20 结点单元构成的空间网格上构建流形覆盖和权函数，采用拉格朗日乘子法施加位移约束条 件，推导了分析静态问题的计算列式。无需细致网格划分即可更精确地分析具有曲线边界的区域，计算结果的误 差能量模较有限元法降低超过一个量级。对采用完全 1 阶覆盖函数时的线性相关性进行了分析，提出在单元角结 点上采用 1 阶覆盖函数基，其他结点上采用零阶覆盖函数基的方法。利用算例分析了用这种方法计算的两套覆盖 函数的收敛率。数值算例表明：效果明显，对解答有体积闭锁的问题均有很高精度。 关键词 数值分析，数值流形方法，拉格朗日乘子法，20 结点单元，体积闭锁 分类号 O 241 文献标识码 A 文章编号 1000-6915(2004)15-2622-06 IMPLEMENTATION OF NUMERICAL MANIFOLD METHOD BY USING 3D TWENTY-NODE ISO-PARAMETRIC ELEMENT MESHES OF FEM Peng Ziqiang ，Ge Xiurun (Institute of Rock and Soil Mechanics ，The Chinese Academy of Sciences， Wuhan 430071 China) Abstract The numerical manifold method (NMM) has been widely used for 2D problems. The construction of numerical manifold covers and weight functions on traditional 3D twenty-node iso-parametric element meshes of FEM is proposed in the paper. With Lagrange multiplier ，Dirichlet boundary conditions are imposed along the essential boundaries. Formulae for static analysis are given. Then it is easier to model 3D problems with curved boundary surfaces by using the NMM without fine meshing in common FEM. Linear dependency of the approximation space is studied，and a method of applying enhanced cover functions only at physical covers connected to corner vertices is proposed to resolve this problem. Thus the scale of discrete system decreases more significantly than that in conventional implementations of the NMM. Example analyses show that the error energy norms decrease dramatically to about ten percent of the co