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

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

功能梯度材料动态断裂力学的径向积分边界元法

采用径向积分边界元法分析功能梯度材料动态断裂力学问题. 该方法使用与弹性模量无关的弹性静力学开尔文基本解作为问题的基本解,在导出的边界-域积分方程中含有由材料的非均质性和惯性项引起的域积分,通过径向积分法将域积分转化为等效的边界积分,得到只含边界积分的纯边界积分方程;从而建立只需边界离散的无内部网格边界元算法. 采用候博特方法求解关于时间二阶导数的系统离散的常微分方程组. 最后通过数值算例验证本文方法的精度和有效性.      

二维弹性力学问题的光滑无网格伽辽金法

计算效率低的问题长期阻碍着无网格伽辽金法(element-free Galerkin method, EFGM)的深入发展.为了提高EFGM的计算速度,本文提出一种求解二维弹性力学问题的光滑无网格伽辽金法.该方法在问题域内采用滑动最小二乘法(moving least square, MLS)近似、在域边界上采用线性插值建立位移场函数;基于广义梯度光滑算子得到两层嵌套光滑三角形背景网格上的光滑应变,根据广义光滑伽辽金弱形式建立系统离散方程.两层嵌套光滑三角形网格是由三角形背景网格本身以及四个等面积三角形子网格组成.为了提高方法的精度,由Richardson外推法确定两层光滑网格上的最优光滑应变.几个数值算例验证了该方法的精度和计算效率.数值结果表明,随着光滑积分网格数目的增加,光滑无网格伽辽金法的计算精度逐步接近EFGM的,但计算效率要远远高于EFGM的.另外,光滑无网格伽辽金法的边界条件可以像有限元那样直接施加.从计算精度和效率综合考虑,光滑无网格伽辽金法比EFGM具有更好的数值表现,具有十分广阔的发展空间.     

Hermite梯度重构核近似配点法及其在功能梯度材料板的静力分析中的应用

由于直接配点法在求解边值问题时边界上的求解精度较低,本文提出了Hermite梯度重构核近似配点法（HGCM）来改进边界求解精度。重构核近似是无网格法中一种常用的近似函数,但是其在求解高阶导数时格式复杂且非常耗时。HGCM采用梯度重构核近似构建形函数的任意高阶导数,提高了计算效率;通过Hermite配点法构建离散方程,提高了边界求解精度。这种方法在求解对应变系数四阶偏微分方程的功能梯度材料板的静力问题时精度高,计算效率高,并可进一步推广应用于高阶偏微分方程描述的边值问题。     

两个新型的八结点块体单元—三维离散算子差分法

给出了弹性力学三维问题的离散算子差分法 ,讨论离散算子差分法在三维问题中的特点 ,意在为该方法的进一步发展提供依据 ,为应用弱形式进行数值求解的研究提供参考。本文从弹性力学平衡方程更为一般的弱形式出发 ,给出了含边界参数的弱形式方程。由该方程不仅可以得到有限元法 ,还可得到离散算子差分法。给出了两个八结点块体单元 ,虽然单元中位移函数是非协调的 ,不需特殊处理便可保证离散格式收敛 ,并对单元位移有十分好的反映能力。     

二维弹性力学问题的光滑无网格伽辽金法

计算效率低的问题长期阻碍着无网格伽辽金法(element-free Galerkin method, EFGM) 的深入发展. 为了提高EFGM 的计算速度, 本文提出一种求解二维弹性力学问题的光滑无网格伽辽金法. 该方法在问题域内采用滑动最小二乘法(moving least square, MLS)近似、在域边界上采用线性插值建立位移场函数; 基于广义梯度光滑算子得到两层嵌套光滑三角形背景网格上的光滑应变, 根据广义光滑伽辽金弱形式建立系统离散方程. 两层嵌套光滑三角形网格是由三角形背景网格本身以及四个等面积三角形子网格组成. 为了提高方法的精度, 由Richardson外推法确定两层光滑网格上的最优光滑应变. 几个数值算例验证了该方法的精度和计算效率. 数值结果表明, 随着光滑积分网格数目的增加, 光滑无网格伽辽金法的计算精度逐步接近EFGM 的, 但计算效率要远远高于EFGM的. 另外, 光滑无网格伽辽金法的边界条件可以像有限元那样直接施加. 从计算精度和效率综合考虑, 光滑无网格伽辽金法比EFGM具有更好的数值表现, 具有十分广阔的发展空间.      

高阶边界元奇异积分的一种通用高效计算方法

针对三维边界元法中曲面单元上的(弱、强、超)奇异积分提出了一种通用高效的计算方法。经极坐标变换,将奇异积分转化为常规积分;采用数值方法计算Cauchy主值积分和Hadamard有限项积分系数;引入保角变换和反曲变换消除因单元畸形或因积分点靠近单元边界而引起的周向积分奇异性。该方法可以统一处理(弱、强、超)奇异积分,并且只需要知道核函数的奇异阶数和少数几个点上的被积函数值,不依赖于积分和函数的具体选取;所需的积分点少,精度高,并且受单元畸形程度影响较小,稳定性好。采用该方法计算了声学和弹性力学中的典型奇异积分,并结合二阶Nystrm方法求解了弹性力学的边界积分方程,验证了方法的高精度和高效性。本文数值积分程序可向作者索取。     

无网格局部强弱法求解不规则域问题

无网格局部彼得洛夫-伽辽金(meshless local Petrov-Galerkin,MLPG)法是一种具有代表性的无网格方法,在计算力学领域得到广泛应用.然而,这种方法在边界上需执行积分运算,通常很难处理不规则求解域问题.为了克服MLPG法的这种局限性,提出了无网格局部强弱(meshless local strong-weak,MLSW)法.MLSW法采用MLPG法离散内部求解域,采用无网格介点(meshless intervention-point,MIP)法施加自然边界条件,并采用配点法施加本质边界条件,避免执行边界积分运算,可适用于求解各类复杂的不规则域问题.从理论上讲,这种结合式方法,既保持了MLPG法稳定而精确计算的优势,同时兼备配点型方法在处理复杂结构问题时简洁而灵活的优势,实现了弱式法和强式法的优势互补.此外,MLSW法采用移动最小二乘核(moving least squares core,MLSc)近似法来构造形函数,是对传统移动最小二乘(moving least squares,MLS)近似法的一种改进.MLSc使用核基函数代替通常的基函数,有利于数值求解的精确性和稳定性,而且其导数近似计算变得更为简单.数值算例结果初步表明:这种新方法实施简单,求解稳定、精确,表现出适合工程运用的潜力.     

Timoshenko梁谐响应求解新方法探索

在空间域上采用只与结点有关的无网格方法离散,在时间域上采用精细积分方法求 解. 无网格离散过程中,利用伽辽金积分等效弱形式代替微分形式的控制方程,并 用修正变分原理满足位移边界条件,采用移动最小二乘法求解离散的形函数,把形 函数代入等效积分弱形式得到离散的二阶方程;精细积分过程中非齐次项采 用Romberg积分. 同时给出了两种不同边界条件的谐响 应求解的两个数值算例,得到了精确的数值结果.     

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

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

加权最小二乘无网格法

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

Singular boundary method for 2D dynamic poroelastic problems

This paper extends a strong-form meshless boundary collocation method, named the singular boundary method (SBM), for the solution of dynamic poroelastic problems in the frequency domain, which is governed by Biot equations in the form of mixed displacement–pressure formulation. The solutions to problems are represented by using the fundamental solutions of the governing equations in the SBM formulations. To isolate the singularities of the fundamental solutions, the SBM uses the concept of the origin intensity factors to allow the source points to be placed on the physical boundary coinciding with collocation points, which avoids the auxiliary boundary issue of the method of fundamental solutions (MFS). Combining with the origin intensity factors of Laplace and plane strain elastostatic problems, this study derives the SBM formulations for poroelastic problems. Five examples for 2D poroelastic problems are examined to demonstrate the efficiency and accuracy of the present method. In particular, we test the SBM to the multiply connected domain problem, the multilayer problem and the poroelastic problem with corner stress singularities, which are all under varied ranges of frequencies.     

自由单元法及其在结构分析中的应用

通过吸收有限元与无网格法的优点,提出了一种新的数值方法------自由单元法.此方法在离散方面,采用有限元法中的等参单元,表征几何形状和进行物理量的插值;在算法方面,采用单元配点技术,逐点产生系统方程.主要特点是,在每个配置点只需要一个和周围自由选择的节点而形成的一个独立的等参单元,因而不需要考虑物理量在单元之间的相互连接关系与导数连续性问题. 本文介绍强形式与弱形式两种自由单元法,前者直接由控制方程和边界条件直接产生系统方程,后者通过在自由单元上建立控制方程的加权余量式产生弱形式积分式,并通过像传统有限元法中的积分过程建立系统方程组.本文提出的方法是一种单元配点法,对于域内点为了获得较高的导数精度,需要采用至少具有一个内部点的等参单元,为此除了可使用各阶次的拉格朗日四边形单元外, 还 给出了七节点三角形等参单元,用于模拟较为复杂的几何形状问题.      

改进的平均源边界节点法模拟平面位势问题

平均源边界节点法ASBNM是一种最近提出的边界型无网格法。该方法仅使用边界节点不涉及任何单元和积分的概念,具有方法简单和程序设计容易等特点。但是,对于依赖于边界积分方程的边界型无网格法,关键问题是如何准确高效地估计影响矩阵的对角元。本文提出直接计算影响矩阵对角元的方法,是已有ASBNM法的改进,将对角元的计算转化为一个纯几何问题,因此适用于任何二维边值问题。数值算例证明了本文方法的有效性和准确性。     

Kinetic meshless method for compressible flows

We present a grid‐free or meshless approximation called the kinetic meshless method (KMM), for the numerical solution of hyperbolic conservation laws that can be obtained by taking moments of a Boltzmann‐type transport equation. The meshless formulation requires the domain discretization to have very little topological information; a distribution of points in the domain together with local connectivity information is sufficient. For each node, the connectivity consists of a set of nearby nodes which are used to evaluate the spatial derivatives appearing in the conservation law. The derivatives are obtained using a modified form of the least‐squares approximation. The method is applied to the Euler equations for inviscid flow and results are presented for some 2‐D problems. The ability of the new scheme to accurately compute inviscid flows is clearly demonstrated, including good shock capturing ability. Comparisons with other grid‐free methods are made showing some advantages of the current approach. Copyright © 2007 John Wiley & Sons, Ltd.     

弹性力学的一种边界无单元法

首先对移动最小二乘副近法进行了研究，针对其容易形成病态方程的缺点，提出了以带权的正交函数作为基函数的方法－改进的移动最小二乘副近法，改进的移动最小二乘逼近法比原方法计算量小，精度高，且不会形成病态方程组，然后，将弹性力学的边界积分方程方法与改进的移动最小二乘逼近法结合，提出了弹性力学的一种边界无单元法，这种边界无单元法法是边界积分方程的无网格方法，与原有的边界积分方程的无网格方法相比，该方法直接采用节点变量的真实解为基本未知量，是边界积分方程无网格方法的直接解法，更容易引入界条件，且具有更高的精度，最后给出了弹性力学的边界无单元法的数值算例，并与原有的边界积分方程的无网格方法进行了较为详细的比较和讨论。     

基于单位分解法的无网格数值流形方法

在数值流形方法和单位分解法的基础上，提出了无网格数值流形方法. 无网格数值流形 方法在分析时采用了双重覆盖系统，即数学覆盖和物理覆盖. 数学覆盖提供的节点形成求解 域的有限覆盖和单位分解函数；而物理覆盖描述问题的几何区域及其域内不连续性. 与原有 的数值流形方法相比，无网格数值流形方法的数学覆盖形状更加灵活，可以用一系列节点的 影响域来建立数学覆盖和单位分解函数，具有无网格方法的特性，从而摆脱了传统的数值流 形方法中网格所带来的困难. 与无网格方法相比，由于采用了有限覆盖技术，试函数的构造 不受域内不连续的影响，克服了原有的无网格方法在处理不连续问题时所遇到的困难. 详细推导了无网格数值流形方法的试函数和求解方程，最后给出了算例，验证了该方法的正 确性.     

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

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

Enriched meshless manifold method for two-dimensional crack modeling

The meshless manifold method is based on the partition of unity method and the finite cover approximation theory which provides a unified framework for solving problems dealing with both continuum with and without discontinuities. The meshless manifold method employs two cover systems. The mathematical cover system provides the nodes for forming finite covers of the solution domain and the partition of unity functions. And the physical cover system describes geometry of the domain and the discontinuous surfaces in the domain. The shape functions are derived by the partition of unity and the finite covers approximation theory. In meshless manifold method, the mathematical finite cover approximation theory is used to model cracks that lead to interior discontinuities in the displacement. Therefore, the discontinuity is treated mathematically instead of empirically by the existing methods. However, one cover of a node is divided into two irregular sub-covers when the meshless manifold method is used to model the discontinuity. As a result, the method sometimes causes numerical errors at the tip of a crack. To improve the precision of the meshless manifold method, the enriched methods are introduced in this work for crack problems.