一种改进的无单元流形方法-中国科技论文在线.PDF

计算力学学报 2002 年 3 月 19 卷 1 期 一种改进的无单元流形方法 张 湘 伟， 蔡 永 昌 (汕头大学工学院，汕头大学，515063） 摘 要 使用 1 阶或 1 阶以上最小滑动二乘法(MLS)形函数的无网格伽辽金法(EFGM)和无单元流形 方法(EFMM)，它们的主要缺点是形函数构造复杂、计算费用十分昂贵。本文提出了一种改进的EFMM 法，它通过采用Shepard形函数(0阶MLS 形函数)对结点的覆盖位移函数加权求和来简化整体近似 位移函数的构造，且能够避免EFGM和EFMM里求解结点形函数时矩阵的求逆及相乘运算。文中的数 值算例表明，这种改进的 EFMM 法收敛快、精度高，与标准的 EFGM 相比其计算时间得到了大幅度 的减少。 关键词 无网格；流形方法；覆盖函数；最小滑动二乘法 分类号 O241 1 前 言 为了减少或消除单元划分给数值分析带来的不便，近年来出现了一些只需结点离散数据， 不需单元信息的无单元方法。按照不同的理论和应用范围，无单元法的种类十分繁多，其中尤 以 Element-free Galerkin Method(EFGM)[1,2]和在此基础上发展起来的 Node-By-Node Meshless Method(NBNM) [3] 、Meshless Boundary Integral Equation(MBIE) [4]等方法因具有求解精度高、 收敛快、能消除体积闭锁和数值解稳定等优点而应用最为广泛。然而即使是这类已被大量研究 的无单元法也仍然存在着许多尚待解决的问题。例如，基于最小滑动二乘法(MLS)的无网格伽辽 金法(EFGM)其整体位移函数的构成十分复杂，形函数的形成需要求矩阵的逆和相乘运算，更为 重要的是现阶段的 EFGM 仍然没有有限单元法(FEM)稳定和快速，其计算花费超过对应的低次有 限单元法 4~ 10 倍以上，这在那些需要迭代计算的弹塑性、大变形、冲击动力学等本来就十分 耗时的问题里，矛盾就更为突出了。因此正如文献[2]所指出的那样，目前无单元法最大的挑战 应该是发展能够如低次有限单元法一样快捷的无单元方法。 本文作者在文献[13]里提出了一种新的无单元流形方法 (EFMM)，它能用较小的结点影响半 径和很少的离散结点得到所求问题的收敛解答，从一定程度上缓解了以往无单元法计算时间长 的矛盾，但是由于 EFMM 仍采用 1 阶MLS 形函数来构造位移函数，因而并没有从根本上解决无单 元法面临的问题。实际上，无单元法的计算效率和结构形式仍能得到进一步的改善和提高。本 文基于流形方法[11,12] [5] 和Partition of Unity(PU) 的理论对文[13]的 EFMM进行了改进，它采用 Shepard形函数(即0阶MLS形函数)对结点的覆盖位移函数加权求和来构造新的EFMM位移函数， 在位移边界上则与基于有限单元网格的数值流形方法(NMM)相耦合来施加位移边界条件，从而较 好地克服了原来无单元法的缺点。以下将详细讨论该方法的原理。 2 整体近似位移函数 2.1 覆盖位移函数 采用 EFMM 寻求边值问题的近似解，首先在多边形域Ω上建立互相迭合的覆盖 C (图1)，其 I 覆盖图形有扇形、圆盘或半圆盘，用数个覆盖图形去覆盖多边形域Ω就形成了一个求解的流形 作者简介：张湘伟(1950~)，男，博士，教授，博士生导师 26 中国科技论文在线 计算力学学报 2002 年 3 月 19 卷 1 期 结构。对各个基本覆盖独立定义覆盖上的位移函数，