关于无单元法中的插值基函数选取的探讨

无单元法不需要单元信息，它采用了一种基于移动最小二乘（MLS）的插值函数。插值基函数对插值函数以及无单元法的计算精度影响很大。本文就不同的基函数对插值函数及无单元法的计算精度的影响作了分析比较，得出了一些有益的结论，并用算例说明了这些结论的正确性。     

弹性力学的插值型边界无单元法

我们讨论了移动最小二乘插值法, 对Lancaster推导的公式进行了改进. 在边界无单元法的基础上, 将边界积分方程方法和本文改进的移动最小二乘插值法结合, 提出了弹性力学的插值型边界无单元法, 推导了相应的公式. 本文改进的移动最小二乘插值法中的形函数具有Kronecker函数的性质, 所以我们提出的插值型边界无单元法可以很容易地直接施加本质边界条件. 我们提出的插值型边界无单元法以节点的真实变量作为未知函数, 是无网格边界积分方程方法的直接解法, 具有较高的精度. 最后给出了数值算例说明了本文方法的有效性.     

Signorini问题的插值型边界无单元法

本文将改进的移动最小二乘插值法和边界积分方程结合,提出了求解Signorini问题的一种新的边界类型无网格方法——插值型边界无单元法.该方法用投影算子处理Signorini问题中的非线性边界不等式条件,然后将Signorini问题归化为边界积分方程,并用改进的移动最小二乘插值法近似未知的边界变量,然后本文分析了该方法的收敛性.数值算例表明该方法在求解Signorini问题时的可行性和有效性,相对于边界元方法也具有更好的精度和收敛速度.     

无单元伽辽金法新形函数技术

针对目前以移动最小二乘技术构造的无单元形函数需要大量的求逆运算,且在边界处无过点插值性质而给计算带来了困难的问题,以泰勒展开理论为基础,继承最小移动二乘法的高阶连续性,用Shepard插值实现"移动最小二乘法的由局部到整体区域的移动性"及"有限元法形函数过点插值性",旨在使无单元伽辽金法的形函数在满足高阶连续性的同时具有过点插值的性质,并避免了现有无单元伽辽金法形函数求解繁琐的缺点.     

地基板无单元法节点分布方案研究

针对基于滑动最小二乘法的板弯曲无单元法进行研究.通过弹性地基板的具体算例,探讨了不同节点的分布方案与指数型权函数影响域大小对板挠度和内力计算精度的影响,并给出了合理布置节点的方案.其对于板弯曲无单元法的进一步研究与应用,具有积极的意义.     

配点型点插值加权残值法

点插值法是一种新型的无单元法，该方法克服了EFGM法中形函数计算复杂、本质边界不容易处理等问题。利用点插值构造试函数，加权残值法求出试函数中的系数，进而得到定解问题的数值解。该方法简化了试函数的选择，适用于岩土工程中的各种数值计算。     

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

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

无单元法在对称叠层板稳定问题中的应用

基于Kircihhoff板理论和移动最小二乘近似函数,进一步研究无单元Galerkin(EFGM)方法在对称叠层板稳定问题中的应用.分析中,本质边界条件采用罚因子法施加,对称叠层板的无单元法几何刚度矩阵由最小二乘法和变分原理得到.通过各向同性板和对称叠层板的数值算例及与其他方法的比较,表明采用EFGM法求解对称叠层板稳定问题具有收敛性好、精度高等一系列优点.     

利用插值型无单元Galerkin方法求解KdV-B方程

首先讨论移动最小二乘插值法,并利用移动最小二乘插值法建立形函数,结合KdV-B方程的Galerkin积分弱形式,提出求KdV-B方程数值解的插值型无单元Galerkin方法(IEFG),并推导其相应的公式,跟无单元Galerkin方法相比,利用插值型无单元Galerkin方法计算时,本质边界条件可直接施加,从而可提高计算效率,并给出算例说明了该方法的有效性.     

二维各向异性位势问题的改进插值型边界无单元法

【目的】把边界积分方程方法和基于非奇异权函数的改进移动最小二乘插值法相结合,建立数值求解二维各向异性位势问题的改进插值型边界无单元法。【方法】在改进移动最小二乘插值法的基础上,讨论了非奇异权函数的改进移动最小二乘插值法,它的形函数满足Kroneckerδ函数的性质,因此可以直接施加边界条件。【结果】数值算例表明该方法求解二维各向异性位势问题是有效和可行的。【结论】与边界元方法相比,该方法精度和收敛性更好。     

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

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

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

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

无单元伽辽金法求解不可压Navier-Stokes方程

用无单元伽辽金法(EFGM)求解了不可压的Navier-Stokes方程,由加权残值法推导了系统无单元伽辽金法离散的Navier-Stokes方程,在时间域上采用分步格式计算,使速度和压力采用同阶线性插值并由相互独立的方程以解耦的形式求解,在每一时间步中,对压力解和速度解采用了Newton-Raphson迭代法进行修正、最后将所得到的方法应用到Couette流中,验证了本文方法的有效性。     

凹凸边界形状热弹性问题的自适应无网格法

为了解决具有凹凸边界形状的平面非定常拟静态热力耦合问题,采用无网格伽辽金算法(EFG)进行求解,使用移动最小二乘法构造形函数及拉格朗日乘子法处理本质边界条件(第一类边界条件),通过引入Voronoi邻接准则和后验误差式,对后续结果进行自适应优化;构建了一种新的适用于非定常拟静态热力耦合问题的EFG法自适应计算模型。为了验证计算模型的可行性,分别计算在二维混合边界条件下光滑与凹凸边界形状平面的温度场以及位移场的分布,并与有限单元法的计算结果进行对比,表征了有限单元法和无网格法计算结果的差异,验证了非定常拟静态热力耦合问题的EFG法计算模型的有效性和精确性。     

弹性大变形问题的复变量无单元Galerkin方法

基于复变量移动最小二乘法,建立了适合于大位移、大转动等弹性大变形问题的复变量无单元Galerkin方法.复变量移动最小二乘法的优点是采用一维基函数构造二维问题的试函数.将复变量移动最小二乘法应用于弹性大变形平面问题,结合大变形问题的Galerkin积分弱形式,采用罚函数法施加本质边界条件,建立了全Lagrange格式下的弹性大变形问题的复变量无单元Galerkin方法,推导了相应的计算公式,数值实现中采用了Newton-Raphson迭代法.最后通过数值算例证明了该方法的有效性.     

三维弹性静力问题的无网格局部Petrov-Galerkin法

将二维平面问题的无网格局部Petrov-Galerkin法拓展到三维的相应理论中,编制了该法相应的三维Fortran程序.分析了均匀受拉立方体和悬臂梁两个经典算例,将所得结果与有限元法和解析解对比.结果表明了无网格局部Petrov-Galerkin法在解决三维弹性静力问题时的可行性和有效性,相对于有限元方法在位移解和应力解上也具有更好的精度.     

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

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

伽辽金最小二乘无网格法在几何非线性问题中的应用

目的在不需要划分单元的情况下求解几何非线性问题。方法伽辽金最小二乘无网格法(MGLS)采用移动最小二乘近似函数作为试函数,并用罚函数法施加本质边界条件,内部区域用最小二乘域,边界区域用伽辽金域,是一种与单元划分无关的无网格方法。在求解几何非线性问题时,采用了增量和修正的Newton-Raphson迭代分析的方法,并在整个分析过程中所有变量的表达格式都采用更新的拉格朗日格式。结果通过对受均布载荷作用的悬臂梁用MGLS法进行内力分析,由于考虑大变形的影响,结构呈现出比线性分析结果刚硬的性质,结果与解析解符合的很好。结论算例表明:MGLS法在求解几何非线性问题时具有可行性,而且计算精度也较好。