自然单元法数值积分方案研究-岩石力学与工程学报.pdf

第24 卷 第11 期 岩石力学与工程学报 Vol.24 No.11 2005 年6 月 Chinese Journal of Rock Mechanics and Engineering June ，2005 自然单元法数值积分方案研究 1，2 1，2 1，2 卢 波 ，葛修润 ，王水林 (1. 中国科学院 岩土力学重点实验室，湖北 武汉 430071 ；2. 中国科学院 武汉岩土力学研究所，湖北 武汉 430071) 摘要：自然单元法采用自然相邻点插值方法在全域构造近似函数和试函数，该方法基于整个求解域内离散节点的 Voronoi 结构。采用标准Galerkin 法形成系统的平衡控制方程时，对弱形式的积分是在Voronoi 图的对偶图Delaunay 三角形内进行的；但自然单元法形函数其局部支撑域与背景积分域不一致，从而导致了相当的积分误差。单位分 解积分方法利用无网格形函数满足单位分解的条件，从而将对弱形式的积分转化到形函数的紧支域内进行；但自 然单元法中形函数紧支域的形状和大小由 Voronoi 结构所限定且形状较为复杂，实现单位分解积分算法需对其形 函数紧支域进行分解映射，计算量较大。若在Delaunay 三角形内采用基于三角形内平均应变的单点积分方案，利 用散度定理可将三角形内平均应变的计算转化为在三角形边界环路上的积分，从而在形成系统矩阵的过程中无需 形函数导数的计算，计算量小而精度高。采用单位分解积分方案，其计算精度和收敛性均好于基于平均应变的点 积分方案，但综合计算精度和计算效率考虑，则基于平均应变的点积分方案较为理想。 关键词：岩土力学；自然单元法；自然相邻节点插值；弱形式；紧支域；单位分解积分；平均应变 中图分类号：TU 43 ；O 241.8 文献标识码：A 文章编号：1000 –6915(2005)11 –1917 –08 STUDY ON THE NUMERICAL INTEGRATION SCHEME WITH NATURAL ELEMENT METHOD 1，2 1，2 1，2 LU Bo ，GE Xiu-run ，WANG Shui-lin (1. Key Laboratory of Rock and Soil Mechanics ，Institute of Rock and Soil Mechanics ，Chinese Academy of Sciences， Wuhan 430071 ，China； 2. Institute of Rock and Soil Mechanics ，Chinese Academy of Sciences，Wuhan 430071 ，China) Abstract ：In natural element method(NEM) ，the trial and test function are constructed with the natural neighbor interpolation(NNT) method. The interpolation is constructed with respect to the Voronoi tessellation of the scattered nodes in the problem domain. The