壳结构的无网格局部Petrov-Galerkin方法

无网格近似函数具有高度光滑性,能够很好的逼近曲壳表面及其位移场。无网格局部Petrov-Galerkin方法不论插值还是离散都不需要单元,是一种真正的无网格方法。本文基于无网格局部Petrov-Galerkin方法的基本原理,采用移动最小二乘插值,利用控制微分方程弱形式,建立了Mindlin壳结构的无网格局部Petrov-Galerkin分析方法,用屋顶壳、受夹圆柱壳、几何非线性圆柱壳作为计算实例分析了求解精度、收敛性和稳定性,并与精确解和有限元计算结果进行了对比,表明该方法计算精度高及收敛性好。     

用局部Petrov-Galerkin法分析薄板自由振动

利用薄板振型方程的等效积分弱形式和对振型函数采用移动最小二乘近似函数进行插值，本文进一步研究了无网格局部Petrov-Galerkin方法在薄板自由振动问题中的应用。它不需要任何形式的网格划分，所有的积分都在规则形状的子域及其边界上进行。在插值近似时，采用虚拟-实际节点值变换方法直接引入本质边界条件。通过数值算例和与其他方法的结果进行比较，表明无网格局部Petrov-Galerkin法求解弹性薄板自由振动问题具有收敛性好、精度高等一系列优点。     

基于Voronoi结构的无网格局部Petrov-Galerkin方法

基于自然邻结点近似位移函数提出了一种用于求解弹性力学平面问题的无网格局部局部Petrov-Galerkin方法。这种方法在结构求解域Ω内任意布置离散的结点，并且利用需求结点的自然邻结点和Voronoi结构来构造整腐朽 求解的近似位移函数，对于构造好的近似位移函数，在局部Petrov-Galerkin方法建立整体求解的平控制方程，这样平衡方程的积分可在背景三角积分网格的形心上解析计算得到，而采用标准Galerkin方法的自然单元法需要三个数值积分点。该方法能够准确地施加边界条件，得到的系统矩阵是带状稀疏矩阵，对软件用户来说，这它学是一种安全的，真正的无网格方法，所得计算结果表明，该方法的计算精度与有限元四边界单元相当，但计算和形成系统平衡方程的时间比有限元法四边界单元提高了将近一倍，是一种理想的数值求解方法。     

轴对称弹性体扭转问题的无网格自然邻接点Petrov-Galerkin法

本文将无网格自然邻接点Petrov-Galerkin 法应用于轴对称弹性体扭转问题的求解．无网格自然邻接点Petrov-Galerkin 法采用自然邻接点插值构造试函数，并且采用三角形线性单元的形函数作为加权残值法的加权函数．自然邻接点插值构造的试函数满足Kronecker delta 函数性质，因此本质边界条件的施加十分方便．由于几何形状和边界条件的轴对称特点，原来的空间问题简化为二维问题求解，因此计算时只需要横截面上离散节点的信息．数值算例结果表明，所提出的方法对求解轴对称弹性体扭转问题是行之有效的．     

带源参数热传导问题的基于滑动Kriging插值的MLPG法

利用基于滑动Kriging插值的无网格局部Petrov-Galerkin(MLPG)法来求解带源参数的二维热传导问题,推导了相应的离散方程。由于滑动Kriging插值法构造的形函数满足Kronecker Delta特性,因此可以直接施加本质边界条件。在离散过程中采用Heaviside分段函数作为局部弱形式的权函数,时间域则通过向后差分法进行离散,这一处理过程中刚度矩阵只涉及到边界积分,而没有涉及到区域积分。最后通过算例验证了本方法的有效性。     

一种高效的局部径向基点插值无网格方法

提出了一种弹性动力分析的高效局部径向基点插值无网格方法(MLRPI).该方法采用径向基点插值形函数近似解变量,运用局部Petrov-Galerkin法推导出了相应的离散方程,并根据波动模拟的精度要求,得到某一结点的动力方程.然后采用Newmark常平均加速度法和中心差分法相结合的显式积分格式进行时域积分,得到每个自由度的一种解耦递推格式.最后,对一平面应变问题进行了求解,比较了该文提出的解耦MI.RPI方法、常规MLRPI方法和ANSYS有限元方法的精度和计算时间,结果表明解耦MLRPI方法与常规MLRPI方法的精度相当,但计算效率大大提高.     

用无网格LRPIM分析中厚板弯曲问题

用径向基函数构造无网格点插值法的形函数,插值函数具有Kronecker delta函数性质,因此可以很方便地施加本质边界条件.利用无网格局部径向点插值方法分别对一个对边固支另对边简支中厚板和一个悬臂中厚板的弯曲进行了分析计算.该方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行,是一种真正的无网格方法.算例表明:将无网格局部径向点插值法应用于计算中厚板的弯曲问题,所求得的位移场和应力场都是光滑的;在径向基函数的基础上,附加多项式大大提高了插值精度;所得结果与弹性力学理论解以及有限元解都十分吻合.     

大变形问题分析的局部Petrov-Galerkin法

在微机电系统(MEMS)的建模和模拟研究中,大变形或大移动要充分予以考虑.用有限元法分析这类问题,由于难以避免的网格畸变,使模拟效率精度降低甚至失效,无网格方法(Meshless Method)则能在分析这类问题时显示出明显的优势,无网格局部Petrov-Galerkin(MLPG)法被誉为是一种有发展前景的真正无网格法.本文进一步发展了MLPG法,通过对任意的离散分布节点采用局部径向基函数构造插值形函数和Heaviside权函数,分析方程采用局部加权弱形式离散,建立了变量仅依赖于初始构型的完全Lagrange分析格式,最后用Newton-Raphson法迭代求解.文中分析了悬臂梁典型算例和微机电开关非线性大变形问题,通过与有限元结果的比较,表明本文提出的大变形问题无网格局部Petrov-Galerkin法具有稳定性好及收敛性快等优点.     

无网格局部Petrov-Galerkin法求解板壳弹塑性大变形

无网格局部Petrov-Galerkin法构造的高阶光滑的形函数非常适合建立板壳结构场函数的逼近函数,是一种比较理想的研究板壳问题的方法。基于Mindlin板壳理论,采用更新拉格朗日原理和大变形条件下场量的无网格表达形式,实现了率型无网格局部Petrov-Galerkin方法对板壳弹塑性大变形的求解,算例分析表明了方法的有效性和较高的分析精度。     

复合材料层合板自由振动分析的无网格自然邻接点Petrov-Galerkin法

基于一阶剪切变形理论,提出了复合材料层合板自由振动分析的无网格自然邻接点Petrov-Galerkin法。计算时在复合材料层合板中面上仅需要布置一系列的离散节点,并利用这些节点构建插值函数。在板中面上的局部多边形子域上,采用加权余量法建立复合材料层合板自由振动分析的离散化控制方程,并且这些子域可由Delaunay三角形方便创建。自然邻接点插值形函数具有Kronecker delta函数性质,因而无需经过特别处理就能准确地施加本质边界条件。对不同边界条件、不同跨厚比、不同材料参数和不同铺设角度的复合材料层合板,由本文提出的无网格自然邻接点Petrov-Galerkin法进行自由振动分析时均可得到满意的结果。数值算例结果表明,本文方法求解复合材料层合板的自由振动问题是行之有效的。     

Local Heaviside‐weighted LRPIM meshless method and its application to two‐dimensional potential flows

In this paper, the local radial point interpolation meshless method (LRPIM) is used for the analysis of two‐dimensional potential flows, based on a local‐weighted residual method with the Heaviside step function as the weighting function over a local subdomain. Trial functions are constructed using radial basis functions. The present method is a truly meshless method based only on a number of randomly located nodes. Integration over the subdomains requires only a simple integration cell to obtain the solution. No element matrix assembly is required and no special treatment is needed to impose the essential boundary conditions. The novelty of the paper is the use of a local Heaviside weight function in the LRPIM, which does not need local domain integration and integrations only on the boundary of the local domains are needed. Effects of the sizes of local subdomain and interpolation domain on the performance of the present method are investigated. The behavior of shape parameters of multiquadrics has been systematically studied. Two numerical tests in groundwater and fluid flows are presented and compared with closed‐form solutions and finite element method. The results show that the use of a local Heaviside weight function in the LRPIM is highly accurate and possesses no numerical difficulties. Copyright © 2008 John Wiley & Sons, Ltd.     

基于局部搜索算法的自然邻接点方法

自然邻接点方法(NNM)采用自然邻接点形函数进行插值，其插值形函数具有严格定义，且与 有限元形函数一样形式简洁、性能优良，因而避免了EFG法里难以准确施加位移边界条件和 材料不连续条件等诸多主要困难. 但是从形式上看自然邻接点方法仍然属于有网格的方法， 其研究和应用受到了较大的限制. 为了克服这个缺点，对于任意给定的数值积分点，提出了 一种基于局部搜索自然邻接点的寻找算法对NNM进行改进. 改进后的NNM与无单元伽辽金法 (EFG)的插值和求解过程类似，兼具有EFG的真正无网格特性及NNM的便于处理边界和材料 不连续条件等优点. 所得计算结果表明，改进后的NNM的计算精度和计算时间与NNM相当， 是一种比较理想的数值求解方法.     

弹塑性结构安定下限分析的无网格局部Petrov-Galerkin法

将基于Voronoi结构的无网格局部Petrov-Galerkin法与减缩基技术相结合,建立了一种安定下限分析的新方法.为了克服移动最小二乘近似难以准确施加本质边界条件的缺点,采用了自然邻近插值构造试函数.通过引入基准载荷域上载荷角点的概念,消除了安定下限分析中由时间参数所引起的求解困难.利用减缩基技术,将安定分析问题化为一系列未知变量较少的非线性规划子问题.在每个非线性规划子问题中,自平衡应力场由一组带有待定系数的自平衡应力场基矢量的线性组合进行模拟,而这些自平衡应力场基矢量可应用弹塑性增鼍分析中的平衡迭代结果得到.算例结果汪明了提出的分析方法的有效性.     

Numerical simulation of ground water flow via a new approach to the local radial point interpolation meshless method

A novel approach to local radial point interpolation meshless (LRPIM) method is introduced to investigate the influence of leakage on tidal response in a coastal leaky confined aquifer system, based on a local weighted residual method with the Heaviside step function as the weighting function over a local sub-domain. The present approach is a truly meshless method based only on a number of randomly located nodes. In this approach, neither global background integration mesh nor domain integration is needed. Radial basis functions (RBFs) interpolation is employed in shape function and its derivatives construction for evaluating the local weak form integrals. Due to satisfaction of kronecker delta property in RBF interpolation, no special treatment is needed to impose the essential boundary conditions. In order to obtain the optimum parameters, shape parameters of multiquadrics (MQ)-RBF are tuned and studied. The leakage has a significant impact on the tidal behaviour of the confined aquifer. The numerical results of this research indicate that both tidal amplitude of groundwater head in the aquifer and the distance over which the aquifer can be disturbed by the tide are considerably reduced by leakage. The novelty of the approach is the use of a local Heaviside weight function in the LRPIM which does not need local domain integration and only integrations on the boundary of the local domains are needed. Therefore, in this research a new local Heaviside weight function has been proposed. Numerical results are presented and compared with the results of analytical solution. It is observed that the obtained results agreed very well with the results of analytical solution. The numerical results show that the use of a local Heaviside weight function in the LRPIM is highly accurate, fast and robust. It is also noticed that this novel meshless approach using MQ radial basis is very stable.     

meshless local Petrov-Galerkin method for static shakedown analysis of elasto-plastic structures

将基于Voronoi结构的无网格局部Petrov-Galerkin法与减缩基技术相结合，建立了一 种安定下限分析的新方法．为了克服移动最小二乘近似难以准确施加本质边界条件的缺点，采 用了自然邻近插值构造试函数.通过引入基准载荷域上载荷角点的概念，消除了安定下限分析 中由时间参数所引起的求解困难．利用减缩基技术，将安定分析问题化为一系列未知变量较 少的非线性规划子问题．在每个非线性规划子问题中，自平衡应力场由一组带有待定系数的 自平衡应力场基矢量的线性组合进行模拟，而这些自平衡应力场基矢量可应用弹塑性增量分析中 的平衡迭代结果得到．算例结果证明了提出的分析方法的有效性．     

h-ADAPTIVE ANALYSIS BASED ON MESHLESS LOCAL PETROV-G ALERKIN METHOD WITH B SPLINE WAVELET FOR PLATES AND SHELLS

Using the two-scale decomposition technique, the h-adaptive meshless local Petrov- Galerkin method for solving Mindlin plate and shell problems is presented. The scaling functions of B spline wavelet are employed as the basis of the moving least square method to construct the meshless interpolation function. Multi-resolution analysis is used to decompose the field variables into high and low scales and the high scale component can commonly represent the gradient of the solution according to inherent characteristics of wavelets. The high scale component in the present method can directly detect high gradient regions of the field variables. The developed adaptive refinement scheme has been applied to simulate actual examples, and the effectiveness of the present adaptive refinement scheme has been verified.     

基于自然单元法的轴对称结构极限上限分析

自然单元法是一种以自然邻近插值为试函数的新兴无网格数值方法,其形函数的计算不涉及矩阵求逆,也不需要任何人为参数。为了充分发挥自然单元法的优势,本文基于极限分析上限定理建立了轴对称结构极限上限分析的整套求解算法。轴对称结构的位移场由自然邻近插值构造,并且采用罚函数法处理材料的不可压条件。为了消除目标函数非光滑所引起的数值困难,采用逐步识别刚性区和塑性区,并对两者用不同方法进行处理。数值算例结果表明,本文提出的轴对称结构极限上限分析方法是行之有效的。     

无网格局部Petrov-Galerkin方法在弹塑性断裂力学问题中的应用

采用无网格局部Petroy-Galerkin方法来分析弹塑性断裂力学问题.这种无网格方法采用移动最小二乘法(MLS)来构造近似试函数和采用Heaviside函数作为加权残值法中的权函数,由于近似函数不满足KroneckerDelta条件,因此采用直接插值法来施加本质边界条件.如果不考虑体力,所形成的整体刚度矩阵只包含局部边界积分,而不包含局部域积分和奇异积分.采用增量Newton-Raphson迭代法来求解弹塑性增量形式的局部Petrov-Galerkin方程.数值算例结果表明,该文方法对于弹塑性断裂力学问题的求解是可行的和有效的,并且所得到的结果具有较好的精度.