基于核重构思想的配点型无网格方法的研究--一维问题

无网格方法按其离散原理可分为Galerkin型、配点型等。其中Galerkin型无网格方法的实施需要背景网格，不属于真正的无网格法；配点型无网格方法的实施不需要背景网格，是真正的无网格法。本文首先介绍了重构核点法的基本原理，然后基于核重构思想，与配点法相结合，以一维问题为例，研究了配点型无网格方法，对该方法构造过程中的近似函数及其导数的计算、修正函数的计算及方法的实现等问题进行了探讨。并结合若干典型算例，检验了其计算精度与收敛姓。     

无网格全局介点法

提出了一种新型无网格法,即无网格全局介点(MGIP)法。该方法采用移动最小二乘核近似来构造形函数,有利于提高数值方法的计算稳定性,而且算法更为简单。该方法需要引入全局介点进行数值离散,基于有限点的广义变分法导出求解系统方程,并采用罚系数法来保证边界条件,数值实现较为简洁。数值算例结果表明:MGIP法的计算耗时不到无网格局部彼得洛夫-伽辽金法的1%,具有较高的计算效率;相比于一般配点法,本文方法的计算稳定性更好,计算精度更高。     

双重点最小二乘配点无网格法

配点类无网格法需要计算近似函数的二阶导数,因而在移动最小二乘(MLS)近似中至少要采用二次基函数。本文利用Voronoi图对双重点移动最小二乘近似法进行了改进,建立了基于Voronoi图的双重点移动最小二乘近似(VDG),并利用加权最小二乘法离散微分方程,导出了双重点最小二乘配点无网格法(MD GLS)。该方法将求解域用节点离散,并以节点为生成点建立Voronoi图,取Voronoi多边形的顶点为辅助点。近似函数及其二阶导数的计算过程可分解为两个步骤:首先用场函数节点值拟合辅助点处近似函数的一阶导数,再以辅助点处近似函数的一阶导数值拟合节点处近似函数的二阶导数。由于在每一步中只需计算MLS形函数及其一阶导数,这种近似方法需要较少的影响点和较小的影响域。同时借助于Voronoi结构的优良几何性质,可以快速地搜索影响点。研究表明,与基于MLS的加权最小二乘无网格法(MWLS)相比,这种方法可以显著提高计算效率,并且在精度和收敛性方面也有所改善。     

加权最小二乘无网格法

在最小二乘法和移动最小二乘近似的基础上提出了加权最小二乘无网格法．该方法除节点外又引入了一些辅助点，控制方程在所有节点和辅助点处的残差用最小二乘法予以消除，边界条件用罚函数法引入．另外对移动最小二乘近似进行了改进，并给出了最小二乘法中泛函的简化格式，因而提高了计算效率．与配点法相比，新方法精度高，稳定性好，并且系数矩阵是对称正定矩阵．与Galerkin法相比，该方法不需要进行高斯积分，因而计算量小．算例表明该方法具有效率高、精度高和稳定性好等优点，并且易于实现．     

求解Helmholtz方程基于核重构思想的最小二乘配点法

基于核重构思想构造近似函数，将配点法和最小二乘原理相结合对微分方程进行离散， 建立了Helmholtz方程的最小二乘配点格式，并分别研究了Helmholtz方程的波传播问题和 边界层问题. 通过数值算例可以发现，给出的数值计算结果非常接近于精确解，计算精度明显高于SPH 法的数值结果，且随着节点数目的增加，其精确度越来越高，具有良好的收敛性.     

基于径向基函数配点法的梁板弯曲问题分析

采用径向基函数配点法分析考虑剪切效应的梁板弯曲问题,该方法利用径向基函数作为近似函数,基于配点法离散方程,通过最小二乘法求解。径向基函数配点法在离散和计算过程中不需要任何形式的网格划分,是一种真正的无网格法;径向基函数可以用一元函数来描述多元函数,存在明显的储存和运算简单的特点;而基于配点法求解不需要积分,提高了计算效率。分析考虑剪切效应的薄梁板问题时,传统的有限元法或无网格法求解均会存在剪切锁闭问题,而径向基函数在全域内存在无限连续性,能够准确地满足Kirchhoff约束条件,因此径向基函数配点法能够消除剪切锁闭现象,而且不会出现应力波动。该方法的优势在于,其不仅易于离散、精度高,而且具有指数收敛率,计算效率高。数值算例验证了上述结论和该方法的稳定性。     

带源参数的二维热传导反问题的无网格方法

利用无网格有限点法求解带源参数的二维热传导反问题，推导了相应的离散方程. 与 其它基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不 需要划分网格，用配点法离散控制方程，可以直接施加边界条件，不需要在区域内部求积分. 用有限点法求解二维热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点. 最后通过算例验证了该方法的有效性.     

无网格法的理论及应用

详细论述了近年来迅速发展的无网格法的理论基础及其在各个领域内的应 用. 无网格法网格依赖性弱, 避免了传统的有限元、边界元等基于网格的数值方法 中可能出现的网格畸变和扭曲, 在一些有限元、边界元等方法难以较好处理的领域体现 出独特的优势. 以加权余量法为主线归纳了已有的30多种无网格法, 各类 无网格法的主要区别在于使用了不同的加权余量法和近似函数. 详尽介绍 了各种无网格近似方案(包括移动最小二乘近似、核近似和重构核近似、单位分 解近似、径向基函数近似、点插值近似、自然邻接点插值近似等)和无网格法 中常用的各类加权余量法(伽辽金格式、配点格式、局部弱形式、加权最小二乘 格式和边界积分格式等), 并讨论了数值积分方法和边界条件的处理等问题. 在 此基础上较系统地总结了无网格法在冲击爆炸、裂纹传播、超大变形、结 构优化、流固耦合、生物力学和微纳米力学等领域的应用, 展示了无网格法相 对于传统数值方法的优势.      

瞬态热传导问题的加权最小二乘无网格法

加权最小二乘无网格法是一种基于节点信息的纯无网格法,该方法使用最小二乘法建立系统的变分原理,通过移动最小二乘法构造近似函数,控制方程在节点处的残量使用最小二乘法予以消除,边界条件通过罚函数法引入。本文推导了瞬态热传导问题的加权最小二乘无网格计算格式,编制了相应的计算程序,算例结果表明,该方法具有精度高、前后处理简单的优点,是一种高效的的新型无网格法。     

紧支试函数加权残值法

将紧支函数引入加权残值法中，提出了紧支试函数加权残值法，其数值格式具有和有限元相似的窄带系数矩阵，提高了加权残值法的计算效率．在紧支试函数加权残值的基础上，导出了紧支试函数直接配点法、紧支试函数Hermite配点法和紧支试函数最小二乘配点法的具体格式，并且对几个典型算例进行了分析．与配点法相比，这些方法精度高，稳定性好，而与Galerkin法相比，这些方法效率高．     

A meshless collocation method for the plane problems of functionally graded material beams and plates using the DRK interpolation

A meshless collocation method is developed for the static analysis of plane problems of functionally graded (FG) elastic beams and plates under transverse mechanical loads using the differential reproducing kernel (DRK) interpolation, in which the DRK interpolant is constructed by the randomly distributed nodes. A point collocation method based on this DRK interpolation is developed for the plane stress and strain problems of homogeneous and FG elastic beams and plates. It is shown that the present DRK interpolation-based collocation method is indeed a truly meshless approach with excellent accuracy and has a fast convergence rate.     

一种新的求解对流扩散边值问题的无网格方法

基于配点法和楔形基函数,提出了一种新的求解对流扩散边值问题的无网格方法。通过一维和二维的问题验证了该数值方法的可行性;并根据数值算例和分析,可以看到该数值方法能达到满意的收敛效果。该数值方法的隐格式形式能够有效地消除对流占优问题的数值振荡现象,是一种真正的无网格方法。     

无网格MPM法三维前处理系统设计

无网格法在解决爆炸冲击、裂纹扩展、高速碰撞以及大变形等问题时,由于与网格无关,能够避免有限元法因网格畸变计算精度降低或计算终止的矛盾,是目前求解固体力学问题的有效方法之一,也是近年来研究的热点。前处理对于无网格法数值计算是非常重要的环节,基于面向对象程序设计思想,应用C 和OpenGL开发了三维图形交互方式的MPM法前处理系统,提出了应用递归四叉树法、扫描表示法和CSG法划分任意形状非等间距的MPM法质点单元,为建立通用的无网格前处理平台提出了有效的解决方案。     

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

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

再生核质点法研究进展

概述了无网格再生核质点法的近似方案和本质边界条件的处理;综述了再生核质点法在结构动力学、大变形分析、计算流体力学、断裂及损伤力学、微机电结构分析中的应用;从工程应用和改进再生核质点法的角度提出了未来的若干研究方向.      

Local Petrov-Galerkin method for a thin plate

The meshless local Petrov-Galerkin (MLPG) method for solving the bendingproblem of the thin plate were presented and discussed. The method used the moving least-squares approximation to interpolate the solution variables, and employed a local symmetricweak form. The present method was a truly meshless one as it did not need a finite elementor boundary element mesh, either for purpose of interpolation of the solution, or for theintegration of the energy. All integrals could be easily evaluated over regularly shapeddomains ( in general, spheres in three-dimensional problems ) and their boundaries. Theessential boundary conditions were enforced by the penalty method. Several numericalexamples were presented to illustrate the implementation and performance of the presentmethod. The numerical examples presented show that high accuracy can be achieved forarbitrary grid geometries for clamped and simply-supported edge conditions. No postprocessing procedure is required to computer the strain and stress, since the originalsolution from the present method, using the moving least squares approximation, is already smooth enough.     

势问题的重构核粒子边界无单元法

将重构核粒子法和势问题的边界积分方程方法结合,提出了势问题的重构核粒子边界无单元法. 推导了势问题的重构核粒子边界无单元法的公式,研究其数值积分方案,建立了重构核粒子边界无单元法的离散化边界积分方程,并推导了重构核粒子边界无单元法的内点位势的积分公式. 重构核粒子法形成的形函数具有重构核函数的光滑性,且能再现多项式在插值点的精确值,所以该方法具有更高的精度. 最后给出了数值算例,验证了所提方法的有效性和正确性. }