基于径向基函数的局部边界积分方程方法.pdf

N0．11 第42卷第11期 机 械 工 程 学 报 V01．42 2006年11月 Nov． 2O06 CHINESEJOURNALOFMECHANICALENGINEERING 基于径向基函数的局部边界积分方程方法水 戴保东．1’2程玉民1 (1．上海大学上海市应用数学和力学研究所上海200072： 2．太原科技大学工程力学系太原030024) 摘要：将基于径向基函数构造的具有插值特性的逼近函数应用于弹性力学问题的局部边界积分方程，推导出相应 离散方程的计算公式，建立基于径向基函数的局部边界积分方程方法。与原有的局部边界积分方程方法相比，该 方法不需要虚拟节点变量，而是采用节点变量的真实解作为基本未知量，是局部边界积分方程无网格法的直接解 法。由于形函数及其导数的构造相对简单，并且满足Delta函数性质，故该方法具有计算量小、精度高，可以像有 限元法一样直接施加边界条件等优点。算例证明了该方法的有效性。 关键词：径向基函数多项式基函数紧支域无网格法局部边界积分方程 中图分类号：0343 光滑性好、计算精度高等优点，同时也带来计算量 O前言 大、边界条件处理不便等问题。尤其是这两种方法 均采用节点变量的近似解为未知量，是边界积分方 无网格法是一种新的求解偏微分方程边值问题 程无网格法的间接列式法。 的数值方法，也是目前计算力学领域的研究热点之 basis 将径向基函数(Ihdial 一【l】。由于该方法具有无网格特性，在处理大变形、 项式基函数耦合，可构造出具有Delta函数性质的 裂纹扩展及动态不连续边界的问题时比有限元法具 形函数，与移动最小二乘法相比，其优点是形函数 有明显优势。目前发展的无网格方法有无单元 构造简单、计算量小，并且如果采用一次多项式基 Galerkin Galerldn法(Element．仔ee memod，EFG)【21、 函数，则形函数具有一阶一致性，同时计算插值函 ken坞1 再生核粒子法(R印roducingpanicle，鼬@)【3J、 数时不容易出现矩阵奇异。这里介绍这种形函数的 有限点法(Finitepointme也od，FPM)04I、无网格局部 构造过程，并应用于局部边界积分方程，对弹性力 localPe仃ov-Galerkin， Petrov—Galerki方法(Meshless 学问题，推导基于径向基函数的局部边界积分方程， MLPG)例、单位分解法(Partitionof吼ity)帕J、点插值 建立相应的离散方程，提出基于径向基函数的局部 法(PominterpolationmeⅡlod，PIM)o川、最小二乘配 边界积分方程的无网格法。最后通过数值算例证明 点无网格法[8】、复变量无网格法[9】以及边界积分方程