一种新数值方法_无网格伽辽金法_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