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

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

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

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

边坡弹塑性稳定性分析的无网格法

将无单元伽辽金法(EFGM) 推广到求解弹—— 黏塑性问题, 自行编制了相应的计算程序, 并应用于土质边坡的弹塑性稳定分析. 通过无单元伽辽金法计算结果与有限元法的计算结果的对比分析, 可以发现用黏弹塑性问题的无单元伽辽金法计算程序去求解弹塑性问题是方便可行的, 本文编制的计算程序稳定性好, 收敛速度快.     

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

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

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

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

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

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

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

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

加权最小二乘无网格法

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

基于无网格法的弹性地基加肋斜板频率数值解

为了给实际地基工程中的经济效益提供技术指标参考,应用无网格法计算加肋斜板在地基上的自由振动行为,应用移动最小二乘近似MLS和一阶剪切变形原理描述加肋斜板的位移场,分别建立斜板与肋条的势能泛函,使用Winkler弹性地基模型处理加肋板与地基之间的接触势能。通过斜板与肋条的位移协调关系寻找斜板和肋条节点参数转换公式,确立加肋斜板的自由振动控制方程。本文运用了无网格法的优势,即使肋条位置改变也不需重置网格。与有限元解的对比验证了本文方法的有效性和精确性。     

无网格全局介点法

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

A MESHLESS LOCAL PETROV-GALERKIN METHOD FOR GEOMETRICALLY NONLINEAR PROBLEMS

Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.     

动力弹塑性分析的无网格自然单元法

基于无网格自然单元法,提出了结构动力弹塑性响应分析的一条新途径.自然单元法是一种新兴的无网格数值计算方法,其实质是基于自然邻近插值的伽辽金法.自然单元法在本质边界条件的施加上较采用移动最小二乘法的无网格法具有明显的优势.在空间域上采用自然单元法离散,并运用加权余量法推导了动力弹塑性分析的离散控制方程.然后,采用预校正形式的Newmark法在时间域上进行求解.最后给出了数值算例,并验证了所提方法的有效性和正确性.     

非均质问题中的无网格边界单元法

介绍了一种不需要内部网格计算非均匀介质问题的边界元算法.该算法是建立在一种能将任何区域积分转换成边界积分的径向积分转换法基础上,首先用对应各向同性问题的基本解来建立以正规化位移表示的非均质问题的积分方程,然后用径向积分转换法将出现在积分方程中的区域积分转换成边界积分,从而形成不需要使用内部网格来计算区域积分的纯边界元算法.与其它无网格法相比,此方法需要很少的内部点,有些问题甚至不需要内部点都能得到满意的结果,因此,可以计算大型的三维非均匀介质工程问题.由于此方法继承了边界元和无网格算法的优点,因而具有广阔的发展前景.     

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

无网格局部彼得洛夫-伽辽金(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使用核基函数代替通常的基函数,有利于数值求解的精确性和稳定性,而且其导数近似计算变得更为简单.数值算例结果初步表明:这种新方法实施简单,求解稳定、精确,表现出适合工程运用的潜力.     

ELASTO－PLASTICITY ANALYSIS BASED ON COLLOCATION WITH THE MOVING LEAST SQUARE METHOD

A meshless approach based on the moving least square method is developed for elasto-plasticity analysis, in which the incremental formulation is used. In this approach, the displacement shape functions are constructed by using the moving least square approximation, and the discrete governing equations for elasto-plastic material are constructed with the direct collocation method. The boundary conditions are also imposed by collocation. The method established is a truly meshless one, as it does not need any mesh, either for the purpose of interpolation of the solution variables, or for the purpose of construction of the discrete equations. It is simply formulated and very efficient, and no post-processing procedure is required to compute the derivatives of the unknown variables, since the solution from this method based on the moving least square approximation is already smooth enough. Numerical examples are given to verify the accuracy of the meshless method proposed for elasto-rdasticity analysis.     

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

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

最小二乘法在薄板塑性动力响应中的应用

本文通过拉普拉斯变换,将最小二乘法应用于薄板的塑性动力响应问题。文中分别以承受均布冲击荷载(中载情形)的简支圆板和方板的塑性动力响应问题为例,说明本文方法的计算方法及步骤,并显示出本文方法的优越性。     

弹塑性结构安定下限分析的无网格局部Petrov-Galerkin法

将基于Voronoi结构的无网格局部Petrov-Galerkin法与减缩基技术相结合,建立了一种安定下限分析的新方法.为了克服移动最小二乘近似难以准确施加本质边界条件的缺点,采用了自然邻近插值构造试函数.通过引入基准载荷域上载荷角点的概念,消除了安定下限分析中由时间参数所引起的求解困难.利用减缩基技术,将安定分析问题化为一系列未知变量较少的非线性规划子问题.在每个非线性规划子问题中,自平衡应力场由一组带有待定系数的自平衡应力场基矢量的线性组合进行模拟,而这些自平衡应力场基矢量可应用弹塑性增鼍分析中的平衡迭代结果得到.算例结果汪明了提出的分析方法的有效性.     

A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years .It uses the moving least square (MLS) approximation as a shape function . The smoothness of the MLS approximation is determined by that of the basic function and of the weight function, and is mainly determined by that of the weight function. Therefore, the weight function greatly affects the accuracy of results obtained. Different kinds of weight functions, such as the spline function, the Gauss function and so on, are proposed recently by many researchers. In the present work, the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method. The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed. Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and a in Gauss and exponential weight functions are in the range of reasonable values, respectively, and the higher the smoothness of the weight function, the better the features of the solutions.     

计算含大位移动边界非定常流动的无网格算法

研究了可以计算处理包含大位移动边界非定常流动的无网格算法. 创建了无网格方法中由于大位移边界运动造成的不合要求点的判断标准, 采用高效的填充方法实现点云重构, 运用线性插值方法得到新点参数, 实现了局部重构处理大位移动点. 流场计算方面, 在计算域自动布点基础上, 采用曲线逼近计算导数及HLLC格式计算数值通量, 发展了求解基于无网格的ALE方程组的算法. 最后, 对活塞问题进行了模拟, 结果与解析解相吻合, 验证了算法的准确性; 另外, 计算比较了气流流过静止圆柱以及圆柱在静止流场运动流场, 结果表明方法是成功的.