混合变换法在无网格伽辽金方法中的应用

移动最小二乘近似的非插值特征给无网格伽辽金方法的应用带来了一定的的困难，本文将再生核质点法中的混合变换法与无网格伽辽金方法相结合，通过对移动最小二乘近似进行局部修正，实现了无网格伽辽金方法中本质边界条件的精确施加。对权函数、影响半径、积分阶等对计算精度的影响进行了有益的探讨。数值计算结果表明了方法的可行性。     

无网格方法求解稳定渗流问题

使用无单元伽辽金(EFG)方法求解圆形油藏中心井和矩形油藏裂缝井两种稳定渗流模型。在中心井模型计算过程中,观察到采用对数等分布置节点是最有效的;计算结果与理论解和有限元解相比较,表明无网格方法是一种比有限元更为精确的方法。裂缝井模型通过在初始节点基础之上加密节点,获得了比较好的结果,并且给出了等压力曲线图。     

无网格伽辽金法求解平面偶应力问题

提出采用无网格伽辽金法(EFGM)求解偶应力问题,以避免有限元求解中因C¹连续要求可能引起的不便。推导了基于二次基和移动最小二乘技术的EFGM计算公式,通过计算受轴向均匀拉伸的带中心小孔无限大板和细长杆的偶应力问题,对所提方法进行了数值验证,与解析解相比结果令人满意,此外还讨论了节点密度和权函数对计算结果的影响。数值结果表明,所提算法可有效地求解平面偶应力问题。     

完全变换法在无网格伽辽金方法中的应用

由于移动最小二乘形函数一般不具有常规有限元或边界元形函数所具有的插值特征,本质边界条件的处理成为无网格伽辽金法实施中的一个难点。本文通过建立节点位移和广义位移之间的关系对移动最小二乘形函数进行修正,给出了修正的移动最小二乘形函数;以二维问题为例,对完全交换法在无网格伽辽金方法中的应用进行了研究,实现了本质边界条件在节点处的精确施加。数值计算结果表明该方法不仅简单合理,而且具有较高的精度、收敛性和稳定性。     

一种改进的无单元方法

使用 1阶或 1阶以上最小滑动二乘法 ( MLS)形函数的无网格伽辽金法 ( EFGM) ,它们的主要缺点是形函数构造复杂、计算费用十分昂贵。本文提出了一种改进的无单元方法 ( IEFM) ,它通过采用 Shepard形函数 ( 0阶 MLS形函数 )对结点的覆盖位移函数加权求和来简化整体近似位移函数的构造 ,且能够避免 EFGM里求解结点形函数时矩阵的求逆及相乘计算。文中的数值算例表明 ,这种改进的 IEFM法收敛快、精度高 ,与标准的EFGM相比其计算时间得到了大幅度的减少     

无网格方法的研究进展与展望

目前正在发展的无网格方法采用基于点的近似,可以彻底或部分地消除网格,因此在处理不连续和大变形问题时可以完全抛开网格重构.无网格方法是目前科学和工程计算方法研究的热点,也是科学和工程计算发展的趋势.本文首先简单地阐述了无网格方法,然后详细叙述了目前提出的各种无网格方法的研究进展,最后对目前无网格方法存在的问题进行了探讨,提出了今后的研究方向.     

稳态斜轧延伸过程的刚塑性伽辽金法模拟

将无网格伽辽金法(Element-Free Galerkin Method,EFGM)与三维刚塑性流动理论相结合,对斜轧延伸过程进行了数值模拟。详细推导了斜轧延伸过程EFGM数值模型的刚度方程,给出了初始速度场和速度边界条件的建立方法和刚性区域的处理技术。得到的轧件的物理形态和金属流动的速度场均与实际情况相符,证明了采用EFGM计算斜轧延伸过程的可行性与正确性,将无网格算法引入到了斜轧领域。并将轧制力和壁厚的计算结果与实验结果进行了比较,结果发现：实验得到的轧制力结果偏小而壁厚值却偏大,其主要原因是由于在建模时做了简化和假定,同时也受轧机弹跳的影响。     

局部正交无网格伽辽金法的研究及其在含多裂纹多孔结构中的应用

通过对无网格法中正交基函数的研究,提出局部正交无网格伽辽金法,局部正交基函数保持原正交基函数的性质,但其导数具有了通式,简洁明了,易于编程实现,计算效率高,并将其应用到求解含多裂纹多孔均匀拉伸板的应力强度因子中,计算结果与用正交无网格伽辽金法和有限元法得到的结果进行比较,证明了局部正交无网格伽辽金法的可行性和正确性。     

影响无网格方法求解精度的因素分析

基于移动最小二乘法的无网格方法的计算精度除受到节点的分布密度和基底函数的阶次影响外,还受到其它因素的影响,其中权函数的选取、权函数影响域的大小及位移边界条件的引入对计算精度影响较大。本文分析了几种常用权函数在数值计算时的特点,包括计算精度、收敛情况、计算效率等,同时分析了影响域大小及边界条件的引入对计算精度的影响。通过分析给出了确定权函数及其影响域大小的方法。当受约束的自由度较多时,通过配点法引入位移边界条件会引起计算结果的振荡,通过施加稳定项可以消除振荡现象,通过对带孔方板的受力分析证明了其可行性。应用以上结论对J23—10曲柄压力机机身进行了受力分析,应力集中部位的计算结果得到了较高的精度。     

弹性力学的复变量无网格流形方法

为克服无网格流形方法配点过多、计算速度慢、容易形成病态方程组等缺点,将复变量移动最小二乘法与无网格流形方法相结合,提出了弹性力学的复变量无网格流形方法。分别采用线性基本与二次基进行计算,并与无网格流形方法相比。研究表明该方法计算量小、精度高。     

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

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

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.     

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

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

无网格法在弹塑性问题中的应用

在弹塑性问题中引入无网格法,得到了弹塑性无网格法的控制方法;计算了中心裂纹板的塑性区分布以及极限荷载;数值结果表明了该方法的可行性和有效性。     

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法分析非线性地基梁

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