自然单元法中的误差分析.pdfVIP

科技信息 博士·专家论坛 自然单元法巾的误差分析 安徽理工大学土木建筑学院 李武 [摘要]自然单元法基于整个求解域内离散结点的Voronoi结构，在全域内构造近似函数和试函数。当采用标准伽辽金法建立系统 存在明显误差。分析自然单元法积分产生误差的各种的原因，并找出新的积分方法解决这一问题。该通过分片试验和悬臂梁等算例 验证新积分方法解决这些误差的可行性和有效性。 [关键词]自然单元法 自然相邻结点插值 蒙特卡罗方法 误差分析 O．引言 方法直接在积分区域上投点。采用随机方法在任意有界区域上投点，该 无网格方法是一种新兴的数值方法，位移插值函数只需要结点信 点函数值的数学期望就是函数在区域上的平均值。虽然随机投点的函 息．不需要单元信息，具有前后处理简单、计算精度高度的优点。自从 数值的数学期望是函数的平均值，但是由于投一点的随机性较大，可以 20 采取投多点的方式，然后计算其平均值，作为函数在区域的平均值。设 I鼓纪90年代初Belysehko进一步对DEM做了改进，并用拉格朗日乘 Galerkin该平均值为m，则积分可以这样计算 子法引入边界条件，称为无单元Galerkin法(Element—free Method，简称EFG)t“．引起研究无网格热潮，各种无网格方法纷纷涌现 出来，它们的区别主要是在位移插值函构造形式上。但是这些研究重点 Jnf(x，y)dfl=mB (1) 都集中在二维空间上，比如自然单元法、伽辽金法、配点法及径向基法 式中，B为区域面积。 等无网格方法[2-31。在二维空间的无网格数值方法研究的比较成熟，并已 还可以把上式右边变化为每一点函数值乘以该区域面积除以投的 应用到实际工程问题中，如大变形问题和裂纹扩展等工程问是酽。 点数n的和，即 自然单元法(NEM)中采用基于求解域内离散结点Voronoi结构的 自然相邻结点插值方法(NNI)构造形函数。采用标准Galerkin方法建立 (2) 』。“x，y)dn=萋f(xi，yJ詈 整体求解的系统平衡控制方程时，在Voronoi图的直线对偶图Delaunay 这种积分方法是真正的完全无网格积分方法。它能很好地解决由 三角形内完成对弱形式的全域积j归。一般采用有限元方法的积分方案 形函数是有理式引起的误差、形函数支持域与积分区域不一致引起的 进行积分计算，这种计算方法对自然单元法形函数的计算误差有很大 误差。 影响。 4．数值算例 本文分析了自然单元法形函数的特性，采用蒙特卡罗积分解决这 4．1不同积分方法的分片试验 些误差。 通过分片试验说明各种积分方法的误差大小。由图2可知，H—mer 1．自然单元法(NEM) 积分在积分点增加时，误差在同一量级上波动，未趋向真实解。然而蒙 自然相邻结点的概念建立在Voronoi图的基础上。图1中左图所示 特卡罗方法在积分点数达到50个前随着积分点数增加，误差出现波 为平面7个结点的Voronol划分及待插值点x的二阶Voronoi结构。自 动，但在50个后随着积分点数增加，数值解趋近真实值。 然相邻结点根据空圆规则确定——若Delaunay三角形的外接圆包含待 插值点，则该三角形的3个顶点即为待插值