叶轮机械计算中动静叶排间的非定常效应

为了在数值分析中正确反映叶轮机械内部动静叶排流场之间的相互干涉，本文解决了两个关键问题。首先，推导了二维、无粘、可压缩、多级流场的通道平均方程并对其中的确定应力项进行了模化。其次，在Ｇｉｌｅｓ理论的基础上，提出一种与时间推时法相匹配的，适用于多级流场计算的无反射边界条件，并从物理本质上解释了它的含义；发展了一种用时间推进法求解多级流场的连续界面方法。采用无反射边界条件比较了连续界面法和混合平面法两种级间参数传递模型的实际效果，因为考虑了动静叶间的干扰，连续界面法的计算结果更接近实际值。     

无网格伽辽金法在板弯曲问题中的应用

无网格伽辽金法采用移动最小二乘近似试函数,形函数一般不具有插值特性,本质边界条件需要特殊处理.本文采用替换式拉格朗日乘子法施加本质边界条件,为提高精度,对修正泛函使用罚函数法再次施加本质边界条件.此方法没有增加未知量的数目,而且刚度矩阵仍具有对称正定带状特点.数值算例表明了该方法的合理性及数值稳定性.     

轴对称几何非线性问题中的无网格算法

应用无网格伽辽金法对轴对称几何非线性问题进行了分析。在小变形假设的条件下，利用几何非线性的应变-位移关系，基于线性弹性本构关系，推导了无网格法的计算控制方程，并采用Newton—Raphson迭代法来求解非线性方程，初步计算了压力管道的几何非线性问题。由于无网格方法中的形函数不具备Kroneckerdelta性质，采用罚方法来实现本质边界条件。数值实例表明．无网格伽辽金法在处理轴对称几何非线性问题时，具有较高的计算精度，是一种有效的数值计算方法。     

基于时段混合能量的Riccati微分方程的求解

基于计算结构力学与最优控制的模拟关系，引入时段混合能的概念， 研究了两类边界条件（对状态约束的边界条件和弹性约束的边界条件）下 Ｒｉｃｃａｔｉ微分方程的求解．通过变分原理，使问题成为两点边值问题；采用惩 罚函数法，对状态约束边界条件用弹性支承边界条件来代替．这样就把两 类边界条件化为了一类进行处理，给出了相应的例题．     

夹层环形板轴对称非线性弯曲和屈曲

由于等厚度夹层环形板大挠度计算边界条件复杂，迄今仅有一种特殊情形的数值解答．本文以三次B样条函数为试函数，用配点法计算夹层环形板的非线性弯曲和屈曲．夹层板采用Reissner模型．讨论了16种边界条件，荷载可为多项式型的分布荷载、内缘均布线荷载、均布边缘力矩或均布径向压力及它们的组合．用非线性理论计算了夹层环形板的压曲临界荷载．样条配点法具有精度高、编写程序通用、输入的数据量少以及节省计算时间的优点．     

第一、二、三类非齐次线性边界条件的齐次化

本征函数法和格林函数法是求解数学物理定解问题的常用方法，这些方法要求问题具有齐次线性边界条件，若定解问题具有非齐次线性边界条件，则必须选择辅助函数，使之化为齐次边界条件，本文首先对稳定的非齐次边界条件进行讨论，得到选择辅助函数的普遍方法和结果；然后对其结果应用参数变易法进行改造，使之适应于不稳定的非齐次边界条件，从而得到选择辅助函数的普遍方法和结论。     

加权残值法求解矩形薄板大挠度问题

采用加权残值法探讨大挠度力学问题时，用5次B样条作为未知函数的试函数，使其直接满足边界条件，这种方法与传统的有限元法相比，具有未知量少，自由度少，连续性强，边界条件容易满足等优点．     

四点法曲面造型的自然边界条件

研究四点插值细分法在曲线曲面造型应用中的边界条件确定问题．针对四点法的特点，对传统几何造型技术中常用的非线性自然边界条件决定方法进行改进，基于等腰梯形四点共圆的事实，提出通过求轴对称点而得自然边界点的线性方法，很好地解决了四点法在曲面造型应用中的关键问题。     

非均质材料热传导问题的扩展无单元伽辽金法

采用滑动克里金(Kriging)插值法构造单位分解函数,并对扩展无单元伽辽金(Galerkin)方法进行了改进.与移动最小二乘法对比,其形函数具备克罗内克(Kronecker)δ函数插值特性,克服了移动最小二乘逼近难以直接准确施加本质边界条件的不足.进一步将该方法应用于非均质材料稳态热传导问题的求解,单夹杂和多夹杂数值结果可以看出:改进的扩展无单元伽辽金法易于施加本质边界条件,只需考虑夹杂几何界面进行节点增强,求解更为方便.     

一种改进的无网格局部Petrov-Galerkin方法

本文提出了一种改进的无网格局部Petrov-Galerkin方法来分析平面弹性力学问题．这种无网格方法采用移动最小二乘近似函数（MLS）来近似试函数,采用Heaviside函数作为加权残值法中的权函数,采用直接插值法来施加本质边界条件．最后通过数值实例表明平面弹性力学问题中改进的无网格Petrov-Galerkin方法具有收敛快、稳定性好、精度高和简单有效的特点．     

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

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

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

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

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

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

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

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

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

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

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

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

无网格法在沥青路面瞬态温度场分析中的应用

为了提高沥青路面温度场的计算精度,采用一种新方法——无网格法对其进行计算分析。根据传热学理论,基于变分原理推出了沥青路面瞬态温度场的无网格法计算公式,采用罚函数法处理本质边界条件。针对工程实例,编制无网格法程序,建立了数值计算模型并进行模拟分析。结果表明:无网格法计算结果与现场实测数据最大误差不超过2.6%,且小于有限元方法最大误差的3.5%。     

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.     

两种点插值法在瞬态涡流分析中的比较

运用多项式点插值法（PPIM）和径向基点插值法（RPIM）构造形函数,推导了适合于工程电磁场瞬态涡流问题的多项式基点插值边界无单元方法（BPPIM）和径向基点插值边界无单元方法（BRPIM）,这两种方法的空间插值形函数满足Kronecker delta条件,从而强加边界条件可以直接施加在边界点上．以金属长方柱的瞬态涡流分析作为数值算例,证实了两种方法的正确性和有效性,并对两种基类的点插值法进行了精度分析和比较．     

无网格法计算误差分析

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