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

摘要：

自然单元法采用自然相邻点插值方法在全域构造近似函数和试函数,该方法基于整个求解域内离散节点的Voronoi结构.采用标准Galerkin法形成系统的平衡控制方程时,对弱形式的积分是在Voronoi图的对偶图Delaunay三角形内进行的;但自然单元法形函数其局部支撑域与背景积分域不一致,从而导致了相当的积分误差.单位分解积分方法利用无网格形函数满足单位分解的条件,从而将对弱形式的积分转化到形函数的紧支域内进行;但自然单元法中形函数紧支域的形状和大小由Voronoi结构所限定且形状较为复杂,实现单位分解积分算法需对其形函数紧支域进行分解映射,计算量较大.若在Delaunay三角形内采用基于三角形内平均应变的单点积分方案,利用散度定理可将三角形内平均应变的计算转化为在三角形边界环路上的积分,从而在形成系统矩阵的过程中无需形函数导数的计算,计算量小而精度高.采用单位分解积分方案,其计算精度和收敛性均好于基于平均应变的点积分方案,但综合计算精度和计算效率考虑,则基于平均应变的点积分方案较为理想.

关键词：岩土力学、自然单元法、自然相邻节点插值、弱形式、紧支域、单位分解积分、平均应变

所属期刊栏目：24

分类号：TU43;O241.8(土力学、地基基础工程)

在线出版日期：2005-08-04（万方平台首次上网日期，不代表论文的发表时间）

页数：共8页

页码：1917-1924