An improved local radial point interpolation method for transient heat conduction analysis

The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method （FEM） as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline （TPS） radial basis functions.

Factors affecting accuracy of radial point interpolation meshfree method for 3-D solid mechanics

Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.

A Local Meshless Method for Two Classes of Parabolic Inverse Problems

A local meshless method is applied to find the numerical solutions of two classes of inverse problems in parabolic equations. The problem is reconstructing the source term using a solution specified at some internal points;one class is that the source term is time dependent, and the other class is that the source term is time and space dependent. Some numerical experiments are presented and discussed.

Three-dimensional magnetotelluric forward modeling using structural-point EFGM and FEM coupling method

In this paper,the authors propose a method of three-dimensional(3D)magnetotelluric(MT)forward modeling algorithm based on the meshfree and finite element coupling method.The model is discretized by regular nodes in the central area,and the radial point interpolation method(RPIM)based on the global weakness is utilized to construct the meshfree shape function.The Governing equations in each background gird are solved by Gaussian integration.In the extended area where the points are sparsely distributed,to avoid the instability of the meshfree method,finite element method(FEM)with regular grids is used to solve the governing equation.Finally,the meshfree and finite element governing equations are coupled by the continuity of the field at the interfaces,and the direct solution technique is used to realize the 3D MT forward modeling.Numerical experiments of several typical electrical models are used to verify the effectiveness of the method.

极限下限分析的区域光滑径向点插值法

极限分析的高效数值计算方法在结构设计和安全评定中具有非常重要的作用.为了更加有效地求解极限分析问题,将区域光滑径向点插值法与二阶锥规划相结合,提出了理想弹塑性结构极限下限分析的一种新方法.将问题域离散为简单的三角形背景单元,每个单元进一步划分成若干个光滑域.为了将复杂的域积分转化为简单的边界积分,并且避免计算形函数的导数,采用广义梯度光滑技术对每个光滑域进行应变光滑处理.由于径向点插值法构造的形函数满足Kronecker delta性质,本质边界条件可以直接施加.依据下限定理,在满足以等效积分弱形式表达的自平衡应力场平衡条件的基础上,用二阶锥规划成功构建了极限分析下限法的计算模型,从而可方便地通过基于原始-对偶内点法的数学规划求解器MOSEK直接求解该问题.数值算例结果表明,本文所提方法有效地克服了维数障碍问题,具有较高的计算精度,并且计算结果对网格畸变十分不敏感.

基于局部区域对应点匹配的MRI图像层间插值方法

目的医学影像层间插值对提高图像分辨率十分重要。基于灰度插值的传统算法在处理结构复杂,层间灰度变化大的区域时常会造成组织细节丢失或边缘模糊。本研究旨在建立一种局部区域对应点匹配插值算法,从而有效改善插值图像质量。方法构建一种局部区域对应点匹配算法。通过将线性插值图像与同位置原始扫描图像之间的对比构建误差图像并确定出高误差区域。对非高误差区域使用线性插值。对于各高误差区分别采用对应点匹配方法处理:首先根据每个区域的局部信息确定其对应于点对匹配度误差函数最小误差的最优参数组,然后为区域内每个插值点找到两侧图像上的最佳匹配点对用于该插值点灰度计算。将此方法应用于轴、矢、冠状位方向系列扫描图像的层间插值,将插值图像与同位置原始扫描图像的均方差作为评价插值质量的标准。结果与其他算法插值结果相比,局部区域对应点匹配算法可显著降低插值误差,矢状位图像插值结果的改善尤为突出。与整体对应点匹配算法插值结果相比,其误差平均值降低达19%,有效改善了原高误差区域的插值质量。结论本研究建立的局部区域对应点匹配插值算法可以有效地解决以往算法在处理复杂组织边界区域插值时的不足,显著提高图像层间插值的精度和质量。

二维水下声辐射问题的改进径向基点插值法研究

径向基点插值法是一种典型的无网格数值计算方法,在分析声学问题时,相比于传统有限元法能更好地抑制频散误差,且在相同的节点分布下通常可以得到更精确的数值解。本文提出一种改进的节点选取方案用于构造插值形函数,即改进径向基点插值法。该方案采取一个简单而直接的格式,可确保在进行数值积分时同一背景积分单元中的被积函数是连续可微的,从而减小数值积分误差,得到比原始径向基点插值法更精确的数值解。同时,为了处理外声场问题,本文采用DtN映射技术将无限域截断为有界计算域,满足索默菲尔德辐射条件。数值试验表明,相比于传统有限元法和原始径向基点插值法,本文改进方法具有更高的计算精度和计算效率,在研究水下声辐射问题时具有良好的应用前景。

ELASTIC DYNAMIC ANALYSIS OF MODERATELY THICK PLATE USING MESHLESS LRPIM

A meshless local radial point interpolation method （LRPIM） for solving elastic dy-namic problems of moderately thick plates is presented in this paper. The discretized system equation of the plate is obtained using a locally weighted residual method. It uses a radial basis function （RBF） coupled with a polynomial basis function as a trial function,and uses the quartic spline function as a test function of the weighted residual method. The shape function has the properties of the Kronecker delta function,and no additional treatment is done to impose essen-tial boundary conditions. The Newmark method for solving the dynamic problem is adopted in computation. Effects of sizes of the quadrature sub-domain and influence domain on the dynamic properties are investigated. The numerical results show that the presented method can give quite accurate results for the elastic dynamic problem of the moderately thick plate.

Coupling Magneto-Electro-Elastic node-based smoothed radial point interpolation method for free vibration and transient analysis of Functionally Graded Magneto-Electro-Elastic structures

A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(FGMEE)structures.By introducing the modified Newmark method,the displacement,electrical potential and magnetic potential of the structures under transient mechanical loading were obtained.Based on G space theory and the weakened weak(W2)formulation,the equations of the multi-physics coupling problems were derived.Using triangular background elements,the free vibration and transient responses of three numerical examples were studied.Results proved that CM-NS-RPIM performed better than the standard FEM by reducing the overly-stiff of structures.Moreover,CM-NS-RPIM could reduce the number of nodes while guaranteeing the accuracy.Besides,triangular elements could be generated automatically even for complex geometries.Therefore,the effectiveness and validity of CM-NS-RPIM were demonstrated,which were valuable for the design of intelligence devices,such as energy harvesters and sensors.

The analysis of composite laminated beams using a 2D interpolating meshless technique

Laminated composite materials are widely implemented in several engineering constructions. For its relative light weight, these materials are suitable for aerospace, military, marine, and automotive structural applications. To obtain safe and economical structures, the modelling analysis accuracy is highly relevant. Since meshless methods in the recent years achieved a remarkable progress in computational mechanics, the present work uses one of the most flexible and stable interpolation meshless technique available in the literature—the Radial Point Interpolation Method(RPIM).Here, a 2 D approach is considered to numerically analyse composite laminated beams. Both the meshless formulation and the equilibrium equations ruling the studied physical phenomenon are presented with detail. Several benchmark beam examples are studied and the results are compared with exact solutions available in the literature and the results obtained from a commercial finite element software. The results show the efficiency and accuracy of the proposed numeric technique.

On the Factors Affecting the Accuracy and Robustness of Smoothed-Radial Point Interpolation Method

In order to overcome the possible singularity associated with the Point Interpolation Method(PIM),the Radial Point Interpolation Method(RPIM)was proposed by G.R.Liu.Radial basis functions(RBF)was used in RPIM as basis functions for interpolation.All these radial basis functions include shape parameters.The choice of these shape parameters has been and stays a problematic theme in RBF approximation and interpolation theory.The object of this study is to contribute to the analysis of how these shape parameters affect the accuracy of the radial PIM.The RPIM is studied based on the global Galerkin weak form performed using two integration technics:classical Gaussian integration and the strain smoothing integration scheme.The numerical performance of this method is tested on their behavior on curve fitting,and on three elastic mechanical problems with regular or irregular nodes distributions.A range of recommended shape parameters is obtained from the analysis of different error indexes and also the condition number of the matrix system.All resulting RPIM methods perform very well in term of numerical computation.The Smoothed Radial Point Interpolation Method(SRPIM)shows a higher accuracy,especially in a situation of distorted node scheme.

MESHLESS METHOD OF DUAL RECIPROCITY HYBRID RADIAL BOUNDARY NODE METHOD FOR ELASTICITY

Combining the radial point interpolation method （RPIM）, the dual reciprocity method （DRM） and the hybrid boundary node method （HBNM）, a dual reciprocity hybrid radial boundary node method （DHRBNM） is proposed for linear elasticity. Compared to DHBNM, RPIM is exploited to replace the moving least square （MLS） in DHRBNM, and it gets rid of the deficiency of MLS approximation, in which shape functions lack the delta function property, the boundary condition can not be applied easily and directly and it＇s computational expense is high. Besides, different approximate functions are discussed in DRM to get the interpolation property, in which the accuracy and efficiency for different basis functions are compared. Then RPIM is also applied in DRM to replace the conical function interpolation, which can greatly improve the accuracy of the present method. To demonstrate the effectiveness of the present method, DHBNM is applied for comparison, and some numerical examples of 2-D elasticity problems show that the present method is much more effective than DHBNM.

一种预报水下声辐射的无网格弱式径向点插值方法

提出了一种预报结构水下声辐射的径向点插值法(radial point interpolation method, RPIM)耦合MDtN(modified Dirichlet-to-Neumann)边界条件的无网格弱式方法。在RPIM-MDtN方法中,采用人工边界将结构外部无限问题域截断为有限计算域,利用RPIM构造声形函数,并在截断边界处施加MDtN边界条件,以提高计算精度和确保声波自由衰减及解的唯一性。在该方法中,声场量的插值无需使用网格或场点的连接属性。研究了该方法性能的影响因素,采用算例对该方法的有效性进行了验证。结果表明:与有限元法相比,该方法具有计算精度高和收敛速度快的优势,对声波数的敏感度显著降低;在要求高精度时,该方法具有计算效率的优势。

有限覆盖径向点插值方法理论及其应用

数值流形方法能够统一地处理连续与非连续变形问题,有限覆盖技术是这种方法的核心。无网格方法前处理过程比较简单,径向点插值法是其中的一种计算格式。本文将有限覆盖技术与径向点插值方法相结合发展了有限覆盖径向点插值无网格方法,综合了数值流形方法与点插值方法的各自优点,能够有效地处理连续与非连续性问题,由此所构造的形函数具有Kronecker δ-函数属性,能够有效地处理位移边界条件。本文在阐述了这种方法基本原理的基础上,通过算例分析与数值计算论证了本文所建议方法的可靠性及其有效性。

起伏地形下各向异性的2D大地电磁无网格法数值模拟

作为一种基于节点的计算方法,无网格法具有构造高阶导数方便,自适应分析便利的优点,特别适合复杂地质构造的数值模拟.本文针对实际地球物理勘探中存在的起伏地形和各向异性的地电结构,提出用无网格法来模拟大地电磁响应,采用复合二次径向基函数构造了形函数,推导了大地电磁无网格法等价线性方程组,研究了系数矩阵的压缩存储方法以及大型稀疏复线性方程组快速求解算法,实现了起伏地形下各向异性的2D大地电磁高精度数值模拟.基于层状模型验证了算法的正确性,计算结果表明:无论是TM模式还是TE模式,计算相对误差均小于1%;通过对地垒和地堑模型的模拟,得出了起伏地形对视电阻率和相位的影响规律;对起伏地形下含有不同各向异性系数异常体的模型进行了数值模拟,为开展复杂地质结构的电磁场特征研究以及地形校正奠定了理论基础.随着计算科学的快速发展,无网格法必将成为新的高精度电磁场数值模拟方法.

基于点插值的配点型无网格法解Helmholtz问题

基于点插值法的思想,用三角函数作为基函数在局部支持域内构造具有Kroneckerδ函数性、单位分解性、高阶连续性、再生性和紧支性的形函数。用配点法离散微分方程,得到了具有稀疏带状性的系数矩阵,用GMERS方法求解代数方程组,分别研究了Helmholtz问题的边界层问题和波传播问题。通过数值算例可以发现,给出的数值结果非常接近于精确解,且随着节点的增加,其精确度越来越高,具有良好的收敛性。

功能梯度材料中的无网格局部径向点插值法

提出一种新的无网格局部径向点插值法来分析功能梯度材料．这种无网格方法采用径向基函数耦合多项式基函数来近似试函数,采用Heaviside函数作为加权残值法中的权函数．构造成的形函数具有Kronecker Delta性质,不再需要额外的处理来施加本质边界条件．若不考虑体力,则所形成的整体刚度矩阵只包含局部边界积分,而不包含局部域积分和奇异积分．在计算过程中,取局部边界积分中的高斯点的材料参数来模拟问题域材料特性的变化．结果表明,这是一种真正的无网格方法,具有模拟简单,计算精度高等优点．

有限覆盖径向点插值法及其在土工问题中的应用

将有限覆盖技术与径向点插值方法相结合发展了有限覆盖径向点插值无网格方法,从而综合了数值流形方法与点插值方法的各自优点,能够有效地处理连续与非连续性问题。用该方法构造的形函数具有Kroneckerδ-函数属性,方便了位移边界条件的处理。在简要阐述了这种方法基本原理的基础上,将该方法用用于地基附加应力计算分析,结果证明本文方法精度高,方法简单有效。

用无网格径向点插值法分析中厚板的弯曲问题

利用无网格径向点插值方法对四边固支和四边简支中厚方板以及悬臂中厚梯形板的挠度和应力进行了分析和计算.编制了该方法的计算机程序,研究了计算结果的精度和收敛性.由于该方法是采用径向基函数耦合多项式基函数来构造形函数,其插值函数具有Kronecker Delta函数性质,可以和有限元法一样很方便地施加本质边界条件.而且该方法是基于节点信息而不是基于单元或网格信息,所以用该方法求解薄板问题时也可以避免剪切自锁现象.算例结果表明,用无网格径向点插值法分析中厚板的挠度和应力问题所得计算结果与已有文献解以及有限元解都十分地吻合,并且具有效率高、精度高、收敛性好和易于实现等优点.

粘性流模拟中基于点云重构的最小二乘无网格法

为改善最小二乘无网格法的粘性流模拟能力,该文提出了一种利用点云重构拟合空间导数的方法。对边界层内的高度各向异性的点云结构进行重构,形成适合最小二乘无网格法计算的近似各向同性的点云结构。在此基础上,采用拉格朗日方法对新点进行插值,AUSM+-UP格式计算点云中心点的数值通量,最小二乘法拟合空间导数,三阶SSP型Runge-Kutta法进行时间迭代。为验证基于点云重构的最小二乘无网格法的准确性和鲁棒性,分别对NACA0012翼型在不同来流条件下的绕流进行了数值模拟,并将计算结果与传统最小二乘无网格法结果进行了对比。结果表明:与传统最小二乘无网格法相比,基于点云重构的最小二乘无网格法获得的分离点位置、涡心位置和表面压力系数更接近于参考结果。