一种新的数值方法_无网格伽辽金法_EFGM_.pdfVIP

第 16 卷第 3 期 计 算 力 学 学 报 V o l. 16 N o. 3 1999 年 8 月 CH IN E SE JOU RN AL O F COM PU TA T ION AL M ECHAN IC S A u gu st 1999 一种新的数值方法——无网格 伽辽金法 EFGM 庞作会, 　葛修润, 　郑　宏, 　王水林 中国科学院武汉岩土力学研究所, 武汉　43007 1 摘　要　无网格伽辽金法 EF GM 是近几年发展起来的与有限元相似的一种数值算法。它采 用移动的最小二乘法构造形函数, 从能量泛函的弱变分形式中得到控制方程, 并用拉氏乘子 满足本征边界条件, 从而得到偏微分方程的数值解。该法只需节点信息, 不需将节点连成单 元; 此外, 还有精度高、后处理方便等优点。本文介绍其基本原理及实现过程, 并用算例表明, 该法具有一定的发展前景。 关键词　无网格; 移动的最小二乘法; 权函数; 拉氏乘子 分类号　 24 1 O 1　前　言 [ 1 ] 有限元、边界元等传统方法是求解边值问题强有力的数值方法 。这类方法以单元作为基 本概念, 一个缺点是, 每次计算前都要剖分网格, 数据准备工作量大, 尤其是对三维问题。当这 类方法用于自适应计算或模拟裂纹扩展时, 一般要不断更新网格 rem e sh in g 。虽然 目前已有 一些网格生成器, 但人们还是觉得准备数据 占用机时多, 不方便[2, 3 ] 。因此, 人们呼唤一种不划 网格的数值方法, 即无网格法 。 m e sh le ss m etho d o r g r id le ss m etho d 无网格法只需节点信息, 不需将节点连成单元。这类方法起源于约 20 年前[ 1, 2 ] , 按其出现 的先后顺序, 有代表性的大致有几下几种: Sm oo th ed P ar t ica l H ydro dyn am ics SPH 1977, 1988 [ 1, 2, 4, 5, 6 ]、 1992 [ 1, 2, 7 ]、 D iffu se E lem en t M etho d D EM E lem en t F ree Ca lerk in M etho d EF GM 1994 [ 1, 2, 3 ] 和 P ar t it ion o f U n ity F in ite E lem en t M etho d PU F EM 1995 [ 1, 2, 8 ] 。 EF GM 与其它几种方法相比具有数值稳定、精度高的优点[2, 9, 10 ] 。 移动的最小二乘法 通过几个互不相关的节点上的值, 拟合出 M ov in g L ea st Squ are M L S 一个函数, 该函数的光滑性好且导数连续[7, 11 ] 。 . 1992 在求解边值问题偏微 B N ay ro le s et a l 分方程数值解时, 用M L S 构造形函数, 这样形函数的形成及区域积分的实现都可以脱离单元 的概念。 等称这一方法为 。 等 1994 对此做了进一步改进, 这 B N ay ro le s D EM T B ely t schko [ 1, 2, 3, 10 ] 些改进包括 : 1 对形函数导数考虑得更全面; 2 采用高阶高斯积分完成区域积分; 3 引入拉氏乘子满足本征边界条件。这些改进使得 求解精度更高、更具发展前途[ 1, 3 ] 。 . D EM T B ely t schko 等称改进后的EDM 为 EF GM 。 EF GM 不仅不需要准备网格, 而且还有精度高、后处理方便、可消除体积闭锁现象 对不 [2, 3 ] 可压缩物体 、收敛快等优点 。 收稿 日期: 1997123 1. 作者简介: 庞作会 1969～ , 男, 博士 ? 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved. 3 期 庞作会等: 一种新的数值方法——无网格伽辽金法 EF GM 32 1 2　E F GM 基本原理 2 1 法 M L S 2 11 基本概念 M L S 在域 中, 数据 u 可拟出函数 uh x : m uh x p j x aj x p T x a x 1a ∑ j 1 式中p T x 是完备多项式基, 对二维域来说 线性基 p T x [ 1　x 　y ], 　m 3 1b 平方基 p T x [ 1　x 　y 　x 2 　x y 　y 2 ], 　m