自然单元法数值积分方案研究.pdfVIP

维普资讯 第24卷 第 11期 岩石力学与工程学报 V_01．24 NO．11 2005年 6月 ChineseJournalofRockMechanicsandEngineering June，2005 自然单元法数值积分方案研究 卢 波 ，葛修润 ，王水林 (1．中国科学院 岩土力学重点实验室，湖北 武汉 430071；2．中国科学院 武汉岩土力学研究所 ，湖北 武汉 430071) 摘要 ：自然单元法采用 自然相邻点插值方法在全域构造近似函数和试函数，该方法基于整个求解域内离散节点的 Voronoi结构。采用标准 Galerkin法形成系统的平衡控制方程时，对弱形式的积分是在 Voronoi图的对偶图Delaunay 三角形 内进行的；但 自然单元法形函数其局部支撑域与背景积分域不一致 ，从而导致了相当的积分误差。单位分 解积分方法利用无 网格形函数满足单位分解 的条件 ，从而将对弱形式的积分转化到形函数的紧支域内进行；但 自 然单元法中形函数紧支域的形状和大小由 Voronoi结构所限定且形状较为复杂，实现单位分解积分算法需对其形 函数紧支域进行分解映射，计算量较大。若在 Delaunay三角形内采用基于三角形 内平均应变的单点积分方案，利 用散度定理可将三角形内平均应变的计算转化为在三角形边界环路上的积分，从而在形成系统矩阵的过程中无需 形函数导数的计算，计算量小而精度高。采用单位分解积分方案，其计算精度和收敛性均好于基于平均应变的点 积分方案，但综合计算精度和计算效率考虑，则基于平均应变的点积分方案较为理想。 关键词：岩土力学；自然单元法；自然相邻节点插值；弱形式；紧支域；单位分解积分；平均应变 中图分类号：TU43；O241．8 文献标识码：A 文章编号：1000—6915(2005)11—1917—08 STUDY oN THENUMERICALINTEGRATIoN SCHEMEW ITH NATURAL ELEM ENTM ETHoD LU Bo ’，GEXiu—mn 2 ， WANG Shui—lin ’ (1．KeyLaboratoryofRockandSoilMechanics，InstituteofRockandSoilMechanics，ChineseAcademyofSciences， Wuha n 43007l， ina： 2．InstituteofRockandSoilMechanics，ChineseAcademyofSciences，Wuhan430071，China) Abstract：Innaturalelementmethod(NEM)，thetrialandtestfunctionareconstructedwiththenaturalneighbor interpolation(NNT)method．Theinterpolation isconstructedwiht respecttotheVoronoitessellationofhte scatterednod esinhteproblem domain．Th eintegrationofhteweakform isperfomr edintheDelauna