配点型广义节点无网格法的基本原理及其应用.pdfVIP

第25卷第l期 计算力学学报 V01．25，NO．1 Chinese of Mechanics 2008 2008年2月 Journal February Computational 配点型广义节点无网格法的基本原理及其应用 栾茂田“一， 叶祥记1’2， 黎 勇1’2， 杨 庆1’2 (1．大连理工大学海岸和近海工程国家重点实验室，辽宁大连116024； 2．大连理工大学土木水利学院岩土工程研究所，辽宁大连116024) 摘 要：借鉴流形方法思想．引入广义节点的概念，对传统的无网格法进行了改进，建立了可具有任意高阶多项 式插值函数的广义节点无网格方法。同时采用径向插值函数构造具有插值特性的逼近函数；采用配点法建立系 统的离散方程。在阐述了这种方法基本原理的同时。针对线弹性力学问题给出了这种方法的数值计算列武。与 传统无网格方法相比，这种方法更具有一般性；同时由于采用了配点法而不需要背景积分网格，所以可以认为这 种方法是某种真正意义上的无网格法。当选取0阶广义节点位移插值函数时便可得到传统的无网格法；在不增 加支持域内节点数目的条件下，通过选取高阶广义节点位移插值函数可以提高计算精度。最后通过算例分析。 对0阶、1阶及2阶广义节点无网格法与现有的有关解答进行了对比，论证了其合理性。 关键词：无网格法；径向插值函数；广义节点无网格法 中图分类号：TU 文献标识码：A 方程的离散形式。在由Galerkin加权余量法所建 l 引 言 立的偏微分方程的弱形式中因积分而需要划分背 近30年来，无网格法得到了快速发展和广泛 景网格，而配置法可以得到控制方程的强形式而无 应用，其具体形式包括光滑粒子法(SPH)[1]，扩散 需背景网格。在石根华E7]所发展的数值流形方法 单元法(DEM)[21，无单元Galerkin法(EFGM)[副， Manifold (NumericalMethod)中，针对广义节点 再生核粒子法(RPKM)E4。，自然单元法(NEM)[5] 提出了广义自由度、数学覆盖和物理覆盖。通过提 和点插值法(PIM)[6]等。这些方法与基于网格的 高定义在覆盖的局部逼近函数即覆盖函数的阶数 有限单元法、有限体积法的最大差别在于逼近函数 可有效地提高流形方法总体解的逼近精度。梁国 的构造方式。在基于网格的传统数值方法中，试函 平[8]将广义节点思想对传统流形方法进行了改进， 数的构造依赖于网格，而无网格法只需采用离散点 并且从数学上证明了这种改进的可行性。栾茂田 直接构造试函数。 和田荣等[。3通过在传统有限单元法中引入广义节 一般地，无网格逼近形函数及其导数的构造和 点概念，建立了可具有任意高阶多项式插值函数的 偏微分控制方程的离散化数值列式是无网格法的 广义节点有限元方法。邵国建[10]等对广义四节点 两个关键步骤。移动最小二乘法(MLS)，再生