轴对称问题中的耦合多项式基的径向点插值法.pdfVIP

第 12期 贾 延，等：轴对称问题中的耦合多项式基的径向点插值法 1715 法很难得到其准确的应力分布，且解析法中常需 为单纯的RBFs，否则为添加了 个多项式基函 要重复构造应力函数，带有一定的偶然性和相当 数的 RBF。在 二维 问题 中，一般 采用线性 大的局限性。随着计算机技术的发展，数值方法 基P ()一71，X，]。 得到了人们的重视。在各种数值方法中，有限元 为避免采用多项式基 PIM所引起的力矩矩 方法是发展最为成熟以及应用最为广泛的一种， 阵奇异性问题，引入了径向基。本文采用复合 2 也为工程实践提供了大量宝贵的解决方案。但是 次(MQ)函数，其形式为： 随着现代工业的发展，有限元法遇到了一些挑战， R ()一R (，_y)一R：(r)一 如动态裂纹的扩展、材料的破坏与失效等不连续 [( + ( ))]，a≥0 (2) 问题以及一些大变形问题，采用有限元法解决是 其中，r7为高斯点影响域中相邻节点之间的距离； 有困难 的。 、 g为形状参数[。； 为与计算点 的局部支持 无网格法l1]的出现突破了传统有限元法基 域中节点间距有关的特征长度，该长度通常等于 于网格划分的限制，在分析动态裂纹扩展[3]、金属 局部支持域中的节点平均间距。 加工成型以及高速碰撞等涉及网格畸变、特大变 在径向基函数R ()中，仅有表示计算点 形和 自适应分析[4l_等方面具有突出优势，现已成 与节点z之间距离的变量为： 为国内外力学界的研究热点。诸多学者通过多种 r一~／(z—lz。)。+(3，～y。)。， 不同途径构造不同的近似场函数并进行积分，从 距离是无方 向的，所 以有 R ( )===R (rk)，即 而衍生出了多种具有各 自特点的无网格法 而耦 R (z，y)一R ( ， )成立 。 合多项式基的径向点插值法¨5 是其中之一。它 (1)式中的待定系数 以和b，可通过 个影响 在点插值方法[7]基础上，引入径 向基函数 (RBF) 域内节点的插值过程得到，如第k个点的插值方 构造近似函数，耦合的近似场函数兼并了多项式 程为： 基函数和径向基函数的优 良特性，确保了RBF的 一 (z ，Yk)一 air ( ，yk)+ 力矩矩阵可逆，同时其形函数满足 Kronecker— i= 1 Delta函数性质，形函数及其导数计算相对简单， m bjpj(xk，Yk)，k一 1，2，… ， 。 使得位移边界条件的施加变得容易，从而克服了 J—l 以往无网格法难 以实现位移边界条件的难点。 (1)式的矩阵形式为： 本文将耦合多项式基的径向点插值法应用于 U。一R0n+Pmb (3) 固体力学轴对称问题中，得到其无网格离散方程， 其中，函数值向量为 数值算例验证了该方法的有效性和适用性，并通 El。一Eul 2 U3 … ]T； 过与ANSYS有限元软件计算结果相比较，表明 RBFs的系数向量为 该方法是求解工程轴对称问题的一种较为有效的 n===[n1 a2 a3 …