Meshless_1无网格方法介绍.docVIP

的过程，在最近十年倍受注意。无网格法的主要目标是摆脱，至少是减轻在整个 结构划分和重分网格的困难，代之以只要在整个结构中增加或消除节点。无网格 法还能减轻某些与有限元法有关的其它问题，例如闭锁、单元畸变等等。 无网格法的原始思想可追溯到用于模拟天体物理现象的平稳粒子流体力学 （SPH）法(Gingold, Monaghan, 1977). 仅当Nayroles, Touzot, Villon (1992)发表了 扩散单元法之后，无网格法的研究才变得非常活跃。其后又发表了许多所谓无网 格法：Belytschko, Lu, Gu (1994)的无单元伽辽金(EFG)法，Liu, Chen, Uras, Chang (1996)的再生核粒子法(RKPM)，Babuska, Melenk (1997)的单位分解有限元法 (PUFEM)，Duarte, Oden(1996)的hp-云法，Suku-mar, Moran, Belytschko (1998)的 自然单元法(NEM)，Wendland (1999)的利用径向基函数的无网格伽辽金法(RBF)。 在本专著的书末提供了对于无网格法详细的文献目录。这些无网格法的主要差别 只是来源于试探函数插值采用的技术。虽然在这些方法中对于求解变量的试探函 数和检验函数的插值都不需要网格，但在这些方法中对于弱形式或能量的积分采 用影子单元是不可避免的。因此，这些方法不是真正的无网格。 最近，对于求解线性和非线性边值问题，Zhu, Zhang, Atluri (1998a, b), Atluri, Zhu(1998a, b) and Atluri, Kim, Cho (1999)发展了两种真正的无网格法，即无网格 局部边界积分方程(LBIE)法和无网格局部彼得洛夫－伽辽金(MLPG)法。在这两 种方法中无论为了对求解变量的试探函数和检验函数的插值，或者为了对弱形式 的积分，在区域或边界都不需要网格，因此这两种方法是真正的无网格法。所有 相关的积分很容易在互相覆盖的、形状规则的区域（在三维问题中通常是球）上 求积。实际上，LBIE法可以简单地看作MLPG法(Atluri, Kim, Cho, 1999)的一种特 例。Zhang, Yao (2001)综合LBIE和边界点法的优点，提出了一种新的规则杂交边 界点法。MLPG法令人注目的成功在一些问题中已经得到了报道，其中包括：求 解对流－扩散问题 [Lin Atluri (2000)]，断裂力学问题 [Kim, Atluri (2000), Ching, Batra (2001)]，Navier-Stokes流动问题 [Lin, Atluri (2001)]，剪切变形梁[Cho, Atluri (2001)] 以及板弯曲问题 [Gu Liu (2001), Long Atluri (2002)] 等。 总之，MLPG是一种真正的无网格法，它不仅包含对于试探函数的一种无网 格插值（诸如移动最小二乘MLS、单位分解PU、Shepard函数或RBF），而且包 含一种弱形式的无网格积分（即所有积分总是在规则形状的子域，例如三维问题中的球、平行六面体或椭球上求积）。 本专著的目的是展示MLPG法，包括它的一般原理和各种实施方案。可以看 到，MLPG法的某些方案不仅完全消除了区域划分网格的人工成本，而且其计算 成本也比有限元法和边界元法低得多。 1.2 整体弱形式和加权余量法(WRM) 为了搭建展示MLPG法的舞台，我们通过一个例题讨论利用整体弱形式的各 种建模策略来开始。可以看到，支配工程中许多线性边值问题的数学描述可以写 成如下形式： 其中 Lm 是一个m阶的线性微分算子。 Lm 可以是一个一维问题的常微分算子，或 者一个二维、三维物理域的偏微分算子。在一般的三维域中，设想一条穿过该域、 连结边界两点的直线。沿着这条直线，可以设想一个m阶的常微分方程。然后， 在边界每点有m/2个边界条件，每个边界条件包含一个不高于m?1阶的微分算子。 图1.1 一个实际的边值问题 例如，我们考虑一个用笛卡尔坐标 xi 定义的二维空间内的线性Poisson方程。 考虑域 ? 如图1.1所示。域 ? 的边界记作 Γ (或者有时记作 ?? )，由曲线A-B-C-D-E-F-A, 以及曲线GH定义。沿着AB段，给定边界条件φ = φ 。类似地， 沿着其它各段的边界条件为：沿着BC段?φ ?n = q ；沿着CD段一部分为φ = φ ， 另一部分为?φ ?n = q ；沿着切口（或裂纹）D-E-F ?φ ?n = 0 ；沿着FA段一部分 为φ = φ ，另一部分为 ?φ ?n = q ；沿着GH孔边 ?φ ?n = 0 。在图1.1中定义的边 值问题是相当实际的工程问题。因此，