对进一步发展点插值插值函数地研究.pdfVIP

对进一步发展点插值插值函数的研究 马文涛, 刘学文 王燕昌 (宁夏大学物理与电气信息工程系, 宁夏银川, 750021) −1 摘 要: 目前,点插值方法面临的最大困难在于如何避免p Q 奇异 修改基函数是解决这个难题的 −1 有效手段 本文提出了固定节点法,使每个子域内的结点数目保持一致 该方法不但可以避免p Q 奇异的问题,还简化了形函数的计算过程 同时指出当节点分布相同时 点插值与有限元二者的 形函数具有相同的表达形式;通过分片试验检验了固定节点法的正确性 关键字 点插值 插值函数 径向函数,固定节点法 1 . 概述 点插值是近期由G.R.LIU 和Y.T.GU[1]提出的一种无网格方法 在该方法中问题域由一系列 节点离散 某点的位移值由该点影响子域内的节点值插值得到 由于该方法具有Delta 函数性质 从而使本征边界条件的施加变得很方便 点插值方法与无单元法都是利用插值函数构造近似解来 逼近真实解 二者有一定的类似 但又有本质的区别 最本质的区别在于二者使用不同的插值函 数 无单元法使用了一种基于MLS 的插值函数 形函数不具有Delta 函数性质 点插值方法则 采用过点插值 使插值过程变得简单 2 . 点插值方法(A point interpolation method) 2.1 点插值公式推导 在域内分布一系列节点 一点 的插值函数u(x) 由该点周围一系列节点插得到: x Q n T u(x , x ) ∑p (x)a (x ) p a(x ) (1) Q i i Q Q i 1 其中 p (x) 是空间坐标x T (x , y ) 的单项式 n 是x 影响域内节点数 a (x ) i Q i Q 是p i (x) 的系数 对已给定的点xQ 有: T a (x ) [a , a , a La ] (2) Q 1 2 3 n 方程 1 中p i (x) 采用PASCAL 三角来建立 因此是完备的 T [ 2 3 n ] ( ) 1, , , , L (3) 对一维情况 p x x x x x 固体力学会议论文集 对二维情况 T [ 2 2 2 2 3 3 ] (4) p (x) 1,x ,y ,xy ,x ,y ,x y ,xy ,x ,y LL 方程(1)中的系数 由 影响域内满足该方程的n 个节点值求得 在节点 处