改进广义移动最小二乘近似的无网格法

无网格法是求解微分方程定解问题的一种新数值方法.移动最小二乘近似只要求近似函数在各节点处的误差的平方和最小,对近似函数导数的误差没有任何约束.而广义移动最小二乘近似要求近似函数及其导数在所有节点处的误差的平方和最小.为了降低计算工作量,本文构造了要求近似函数在全部节点处和任意阶导数在部分节点处误差的平方和最小的改进广义移动最小二乘近似.数值计算显示本文提供的方法关于函数值和各阶导数值都具有很高的精度.     

加权最小二乘无网格法

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

基于径向基函数的加权最小二乘无网格法

径向基函数插值是一种新型的无网格插值方法,具有形式简单、空间维数无关等优点.这种插值方法具有δ函数的性质,易于满足本质边界条件,且插值函数的导数求解过程比通常的移动最小二乘插值(MLS)简单,精度也较高.另一方面,通过加权最小二乘法离散控制方程不需要积分,具有效率高,精度高等优点.本文试图将两者的优点结合起来,发展一种...     

矩形加肋板线性弯曲分析的移动最小二乘无网格法

提出了一种求解矩形加肋板线性弯曲问题的移动最小二乘无网格法。将矩形加肋板模拟成平板与肋条组成的复合结构。先基于一阶剪切变形理论,由移动最小二乘近似建立板和肋条的位移场,再利用板与肋条叠合处的位移协调条件,推导肋条与板的节点参数转换方程,最后利用此方程将板与肋条的应变能叠加,推导出整个加肋板的刚度方程。由于本文提出的加肋板无网格模型中不涉及到网格,肋条不必像有限元那样必须沿网格线布置,肋条位置的改变也不会导致板网格的重新划分。文末算例表明由本文方法得到的解与采用实体单元得到的ANSYS有限元解吻合良好,证明了本文方法的准确性。     

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

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

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

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

基于核重构思想的最小二乘配点型无网格方法

介绍重构核点法的基本原理和近似函数的构造方法，并基于核重构思想，应用配点法和最小二乘原理，离散微分方程，建立求解的代数方程，提出了一种基于核重构思想的最小二乘配点型无网格方法．与一般配点法相比，该方法的系数矩阵是有对称正定的，计算精度高，稳定性好．该方法的实施不需要背景网格，不需要进行高斯积分，与Galerkin法相比，具有计算量小、边界条件处理简单的特点，是一种真正的无网格法．对该方法构造过程中的近似函数及其导数的计算、修正函数的计算及方法的实现等问题进行了探讨．文中结合若干典型算例，检验了该方法的有效性．     

一种求解三维声腔声压分布问题的无网格法

利用传统有限元法求解声压分布问题常常受到污染误差和色散误差的困扰.加权最小二乘无网格法(MWLS)是一种基于移动最小二乘(MLS)近似的无网格方法,求解声腔声压分布问题具有低色散、高精度的特点.然而传统的MLS近似有时容易产生病态矩阵,利用加权正交基函数构建改进的移动最小二乘(IMLS)近似,得到的系统方程为非病态的.论文基于改进的加权最小二乘无网格法(IMWLS)求解三维声腔内部声压分布.计算得到的声压分布和声压频响曲线都与参考值十分吻合,峰值误差和污染误差都比FEM的小,计算成本相比无单元伽辽金法显著降低.计算结果表明IMWLS相比传统的FEM,能在更高的频段内达到高精度,并且相比EFGM能大幅提高计算效率.     

无网格全局介点法

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

虚边界元最小二乘配点法

本方法是在虚边界上进行数值积分，在实边界上有限个点处满足给定问题边界条件的一种数值算法。在虚边界上积分是为能寻求到使原问题得到正确解的较好的分布虚体力；在实边界上配点，是依据加权残数法中超额配点，即最小二乘配点法的思想。文中给出了本文方法的基本思想，并给出了壳体的算例。由数值结果表明本文方法的计算精度是令人满意的。     

用无网格局部Petrov-Galerkin法分析非线性地基梁

利用无网格局部Petrov-Galerkin法求解了非线性地基梁。在Petrov-Galerkin方法中，采用移动最小二乘（MLS）近似函数作为场主量挠度的试函数并取移动最小二乘近似函数中的体验函数作为近似场函数的加权函数，采用罚因子法施加本质边界条件。文末给出了两个计算实例，算例的结果表明，Petrov-galerkin法不仅能成功地分析线性地基梁，而且也适用于解非线性地基梁，在分析非线性地基梁时具有收敛快，稳定性好的优点。     

Meshless method based on collocation with consistent compactly supported radial basis functions

Based on our previously study, the accuracy of derivatives of interpolating functions are usually very poor near the boundary of domain when Compactly Supported Radial Basis Functions (CSRBFs) are used, so that it could result in significant error in solving partial differential equations with Neumann boundary conditions. To overcome this drawback, the Consistent Compactly Supported Radial Basis Functions (CCSRBFs) are developed, which satisfy the predetermined consistency conditions. Meshless method based on point collocation with CCSRBFs is developed for solving partial differential equations. Numerical studies show that the proposed method improves the accuracy of approximation significantly. The project supported by the National Natural Science Foundation of China (10172052)     

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

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

应力高梯度问题的无网格分析

基于移动最小二乘法的无网格计算，采用线性基函数即可得到C^1连续位移场，使得应力，应变场在整个求解域内保持连续；节点之间脱离了单元的约束，对求解域进行离散和加密节点时变得十分灵活，因此适合分析应力高梯度问题。本文简要介绍了无网格方法的基本原理，给出了确定节点影响域大小的方法，应用无网格方法对带有V型缺口的受拉方板J23-10曲柄压力机机身进行了受力分析，得到的应力集中部位的计算结果与实际值更为接近。     

MESHLESS ANALYSIS FOR THREE-DIMENSIONAL ELASTICITY WITH SINGULAR HYBRID BOUNDARY NODE METHOD

The singular hybrid boundary node method (SHBNM) is proposed for solving three-dimensional problems in linear elasticity. The SHBNM represents a coupling between the hybrid displacement variational formulations and moving least squares (MLS) approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the later. The rigid movement method was employed to solve the hyper-singular integrations. The 'boundary layer effect', which is the main drawback of the original Hybrid BNM, was overcome by an adaptive integration scheme. The source points of the fundamental solution were arranged directly on the boundary. Thus the uncertain scale factor taken in the regular hybrid boundary node method (RHBNM) can be avoided. Numerical examples for some 3D elastic problems were given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method were studied through the numerical examples.     

一种新的数值方法——无网格伽辽金法（EFGM）

无网格伽辽金法（ＥＦＧＭ）是近几年发展起来的与有限元相似的一种数值算法,它采用移动的最小二乘法构造形函数,从能量泛函的弱变分形式中得到控制方程,并用拉氏乘子满足本征边界条件,从而得到偏微分方程的数值解中得到该法只需节点信息,不需将节点连成单元,此外,还有精度高,后处理方便等优点,本文介绍其基本原理及实现过程,并用算例表明,该法具有一定的发展前景。     

NUMERICAL ANALYSIS OF MINDLIN SHELL BY MESHLESS LOCAL PETROV-GALERKIN METHOD

The objectives of this study are to employ the meshless local Petrov-Galerkin method （MLPGM） to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many factors, including the dimension of compact support domain, the dimension of quadrture domain, the number of integral cells and the number of Gauss points. These factors＇ sensitivity analysis is to adopt the Taguchi experimental design technology and point out the dimension of the quadrature domain with the largest influence on the computational accuracy of the present MLPGM for shells and give out the optimum combination of these factors. A few examples are given to verify the reliability and good convergence of MLPGM for shell problems compared to the theoretical or the finite element results.