移动最小二乘形函数插值精度

移动最小二乘近似作为无网格法中广泛采用的形函数构造方法,其插值精度直接决定数值分析的质量.移动最小二乘形函数的性质通过编写的程序进行计算验证和讨论,重点分析了形函数插值精度对各影响因素的敏感性,并对已知函数的表面拟合进行检验,给出了合理的参数取值与选择范围.研究结果表明,权函数形状、支持域尺度、基函数形式和插值点密度等,对移动最小二乘形函数的插值稳定性和插值精度均有重要影响.     

插值样条δ-序列求解非线性对流扩散方程

提出了一种用广义函数δ序列求解偏微分方程的数值方法.首先对一阶B样条函数N1(x)进行卷积得到四阶B样条函数N4(x),用N4(x)的线性组合构造出三次样条插值基函数;然后用样条插值基序列逼近δ函数,利用δ函数的性质构造插值样条δ序列,该δ序列具有对称、Riesz基和插值性质.以非线性对流扩散方程(伯格方程)为例,用插值样条δ序列离散该方程的空间形式,用四阶龙格库塔方法描述发展过程,取得了较好的精度.为减少计算量,加快插值函数的收敛速度,进一步提高求解精度,对δ序列进行了改进,对同一算例进行数值实验,结果表明,改进后的算法求解过程稳定发展,能够有效描述局部快速变化的情况.     

二次基函数的构造及其应用

本文构造了具有二次代数精度的基函数，这些基函数无论在计算方面还是在应用方面都比Ｅ．Ｌ．Ｗａｃｈｓｐｒｅｓｓ所建立的方便．利用这些基函数，得到了三角剖分下插值多元样条ｕ（ｘ，ｙ）∈Ｓ（△，Ｄ）存在的充要条件，并以数值便阐明了求解插值多元作条空间的维数问题．     

利用点插值函数作为试函数的最小二乘配点法

采用径向基函数与多项式基函数作为耦合的基函数,并利用点插值法构造加权残值法中的近似试函数,试函数中的形函数具有狄拉克-δ函数性质,因此可以直接施加本质边界条件．利用这种试函数和采用最小二乘配点法求解了一维二阶微分方程和薄板的弯曲问题,并与理论结果进行对比;同时还检验了配点数以及节点支持域半径对计算精度的影响．数值结果表明：这是一种与单元划分无关的无网格方法,具有模拟简单,计算精度高,收敛快的优点．     

基于Taylor展开的无单元插值形函数及应用

在无单元伽辽金法的基础上,构造了基于Taylor展开的具有过点插值的无单元形函数,它可以和有限元法一样处理边界条件,克服了传统的无单元伽辽金法遇到的瓶颈问题;对非凸边界的处理,提出了新的准则--弧弦准则(arc-string criterion).这样,可大大减少了无单元法的计算工作量,提高了边界处理的精度,并且继承了无单元法及有限元法的优点.     

用高阶有限元素法模拟双侧向测井仪器在层状介质中的响应

利用有限元法模拟双侧向测井响应,基函数的特性对有限元计算效率和精度影响较大.当基函数为线性函数时,要达到理想的仿真精度需要采用十分稠密的网格剖分,这将导致计算量增加从而影响仿真效率.采用高阶基函数插值,模拟双侧向测井响应,结果表明,稀疏的网格剖分在采用高阶基函数时仍然可以保证较好的仿真精度,同比线性插值基函数,高阶基函数具有高效率和高精度特点.针对深侧向截断误差较大的问题,通过优化网格划分和边界处理方法可提高深侧向结果精度.     

三角域上C1连续的插值样条函数

通过引进三角域上的插值基函数，给出了一种新的三角域上的二元三次插值样条函数，这种插值样条函数整体达到C^1连续，且在各网格点处的参数可由递推公式得到。文中给出的插值样条函数较之Farin提出的分裂三角形方法，具有计算方便、待定系数少且参数易于确定等优点，更易于在CAD中应用。     

基函数的阶数对无单元伽辽金方法精度的影响

无单元伽辽金(Element-Free Galerkin)方法是无网格方法中很重要的一种数值计算方法,利用无单元伽辽金方法求解二维稳态热传导方程,当选取基函数为线性基、二次基时,分别将数值解和解析解对比,分析了基函数的阶数对无单元伽辽金方法精度的影响,并说明无单元伽辽金方法是一种高精度的数值计算方法 .     

基函数应用于插值多元样条

本文借助于代数几何技巧构造了具有三次代数精度的基函数。这些基函数计算起来要比E.L.Wachspress定义的楔函数简便且可用于确定多边形域D在三角剖分Δ下空间S′3(Δ,D)中插值多元样条U(x,y)存在的充分必要条件。     

Multi-Quadric函数与Gauss函数的插值比较

径向基函数在散乱数据插值中有着广泛应用.本文对Multi-Quadric函数与常用的Gauss函数列举多个实例,使用Mathematica编程计算,对其参数进行分析与比较,获得效果较好的误差图.     

无网格伽辽金法(EFGM)求解接触问题

本文考虑到无网格伽辽金法只需节点信息,不需将节点连成单元,并且有精度高,后处理方便等优点,从而用它求解接触问题。这里采用的方法是将Ｋａｔｏｎａ界面单元引入ＥＦＧＭ,迭代求得两物体间的接触状态。算例表明,本文方法基本可行。     

无网格法计算误差分析

探讨了无网格法中形函数的性态及对计算结果的影响 ,讨论了无网格法产生误差的原因 .主要分析了无网格伽辽金法 (EFGM )节点不良分布以及采用一般高次多项式基构造形函数时 ,致使形函数中矩阵A(X)病态 ,从而导致全局数值解振荡的原因 .就不同的基函数对插值函数及无网格法的计算精度的影响作了分析比较 ,得出了基函数的选取标准 ,算例说明使用三次基函数计算精度最高 .     

Meshless Least-Squares Method for Solving the Steady-State Heat Conduction Equation

The meshless weighted least-squares (MWLS) method is a pure meshless method that combines the moving least-squares approximation scheme and least-square discretization. Previous studies of the MWLS method for elastostatics and wave propagation problems have shown that the MWLS method possesses several advantages, such as high accuracy, high convergence rate, good stability, and high computational efficiency. In this paper, the MWLS method is extended to heat conduction problems. The MWLS computational parameters are chosen based on a thorough numerical study of 1-dimensional problems.Several 2-dimensional examples show that the MWLS method is much faster than the element free Galerkin method (EFGM), while the accuracy of the MWLS method is close to, or even better than the EFGM These numerical results demonstrate that the MWLS method has good potential for numerical analyses of heat transfer problems.     

轴对称问题中的无网格Galerkin法

采用无网格Galerkin法分析轴对称问题，得到弹性力学中的对称问题的无网格离散方程．将这一方法与有限元耦合，即在边界处布置有限单元，这样就可以用传统有限元方法方便地处理力学边界条件．算例考察表明：本文方法通过了分片检验，计算结果达到了较高的精度，最大误差不超过5％．     

采用无网格伽辽金法求解电容层析成像正问题

由于有限元法求解电容层析成像正问题的计算准备及后处理非常费时，对正问题的三维求解造成了瓶颈，为此，提出采用无网格伽辽金法求解电容层析成像正问题，获得正问题的弱变分形式，并用拉格朗日乘子法施加边界条件，从而得到数值解．在同样的仿真条件下，2种方法的计算时间分别为14．046S和5．078S．对5种典型流型进行仿真，结果表明，2种方法计算结果的最大相对误差为2．25％．因此，无网格伽辽金法与有限元法具有相当的精度，且计算速度有较大提高．     

模拟平板轧制的新方法——无网格伽辽金法

将一种新的数值方法无网格伽辽金法(EFGM)用于刚塑性可压缩材料稳态轧制过程的模拟,由于形函数不满足插值条件,采用罚函数法满足本质边界条件;为提高精度,选用矩形影响域的张量积核函数;利用有限元背景网格作为积分单元,对求解域内和边界上采用不同的高斯积分方案·数值计算结果与刚塑性有限元的计算结果和文献中的实验数据吻合较好,说明无网格伽辽金法用于刚塑性可压缩材料轧制过程的可行性和正确性·     

无网格伽辽金法在裂纹扩展计算中的应用

采用了一种基于t-分布的新型权函数，提高了无网格伽辽金法的计算精度；采用完全变换法处理本质边界条件，实现了本质边界条件在节点处的精确施加；针对裂纹扩展中的实际情况，对动态影响半径法作了进一步的补充和改进．算例验证了方法的正确性和有效性．     

金属三维塑性成形过程无网格伽辽金法数值模拟技术

将无网格伽辽金法(EFGM)与三维刚(粘)塑性流动理论相结合,对EFGM在金属三维塑性成形过程数值模拟中的应用技术进行了研究.分别采用边界奇异权法和修正的罚函数法处理速度边界条件和体积不可压缩条件,采用反正切摩擦模型处理摩擦边界条件,推导了金属三维塑性成形过程EFGM法数值模拟的刚度方程,给出了关键算法.对长方体金属镦粗过程进行了数值模拟,并将数值结果与三维刚塑性有限元体积成形商品软件Deform3D计算结果作了比较.发现两者吻合良好,表明了本文方法的正确性和有效性.     

无网格伽辽金法在二维结构问题中的应用研究

无网格伽辽金法(EFGM)是一种新型的求解偏微分方程的数值计算方法,不需对结构进行有限元网格的离散化,只需节点信息而不需将节点连成单元.本文论述和研究了EF-GM的基本原理与实现过程,主要包括用移动最小二乘法(MLS)构造形函数、用变分原理推导控制方程、用拉格朗日乘子法增强本征边界条件和域的高斯积分4个主要过程.基于MAT-LAB平台,实现了二维弹性结构的EFGM算法,并将典型算例的EFGM求解结果与有限元近似解、解析解结果进行了比较,结果表明了EFGM算法的正确性和有效性.     

基于基本解和径向基函数插值的热弹性无网格分析

采用基于基本解方法和径向基函数插值的无网格算法(MFS-RBF)分析了广义的热弹性问题.位移分为齐次解和特解两部分,径向基函数被用来插值任意温度变化和体力分布,对应于径向基函数的特解集被用来构造位移特解部分,而基本解方法用来计算相应的位移齐次解.最后,常见的重力荷载、惯性荷载和温度荷载算例验证了算法的有效性和简单性.