二维弹性力学问题的光滑无网格伽辽金法

计算效率低的问题长期阻碍着无网格伽辽金法(element-free Galerkin method, EFGM)的深入发展.为了提高EFGM的计算速度,本文提出一种求解二维弹性力学问题的光滑无网格伽辽金法.该方法在问题域内采用滑动最小二乘法(moving least square, MLS)近似、在域边界上采用线性插值建立位移场函数;基于广义梯度光滑算子得到两层嵌套光滑三角形背景网格上的光滑应变,根据广义光滑伽辽金弱形式建立系统离散方程.两层嵌套光滑三角形网格是由三角形背景网格本身以及四个等面积三角形子网格组成.为了提高方法的精度,由Richardson外推法确定两层光滑网格上的最优光滑应变.几个数值算例验证了该方法的精度和计算效率.数值结果表明,随着光滑积分网格数目的增加,光滑无网格伽辽金法的计算精度逐步接近EFGM的,但计算效率要远远高于EFGM的.另外,光滑无网格伽辽金法的边界条件可以像有限元那样直接施加.从计算精度和效率综合考虑,光滑无网格伽辽金法比EFGM具有更好的数值表现,具有十分广阔的发展空间.     

质点积分无单元伽辽金法及其在金属挤压过程中的应用

无单元伽辽金法需要在背景网格上积分,计算量大.节点积分无单元伽辽金法把对求解域的积分转化为对节点的求和,效率高,但因零能模态不受控制而会产生不稳定现象,需要采取一定的稳定化方案.本文采用应力点思想,通过Newtor-Cotes法计算积分,建立了质点积分无单元伽辽金法,并通过小变形弹性静力学问题说明了该方法具有良好的稳定性,且计算效率远高于无单元伽辽金法.最后本文将质点积分无单元伽辽金法成功地应用于三维金属挤压成型过程的数值模拟,显示了该方法在分析此类问题时的优势和潜力.     

伽辽金型无网格法的数值积分方法

无网格法直接通过节点信息构造形函数,不依赖于节点之间的有序单元连接,能够建立任意高阶连续的整体协调形函数。与传统的有限元法相比,无网格法对大变形问题、移动边界问题和高阶问题的求解方面有比较明显的优势。伽辽金型无网格法是目前应用最为广泛的一类无网格法。虽然无网格形函数本身不依赖于单元,但伽辽金型无网格法需要采取合适的方法进行弱形式的数值积分。由于无网格形函数一般不是多项式,具有非插值性且影响域与背景积分网格通常不重合,伽辽金型无网格法通常需要采用高阶的高斯积分进行数值积分,导致了计算效率低下,难于求解大型实际问题。因此,如何通过建立高效积分方法提高无网格法的计算效率成为无网格法研究领域的一个核心问题。本文对无网格法中强制边界条件的施加方法进行了简要归纳,总结了伽辽金型无网格法中若干常用的数值积分方法,并对伽辽金型无网格法的数值积分方法领域存在的一些问题进行了探讨。     

二维弹性力学问题的光滑无网格伽辽金法

计算效率低的问题长期阻碍着无网格伽辽金法(element-free Galerkin method, EFGM) 的深入发展. 为了提高EFGM 的计算速度, 本文提出一种求解二维弹性力学问题的光滑无网格伽辽金法. 该方法在问题域内采用滑动最小二乘法(moving least square, MLS)近似、在域边界上采用线性插值建立位移场函数; 基于广义梯度光滑算子得到两层嵌套光滑三角形背景网格上的光滑应变, 根据广义光滑伽辽金弱形式建立系统离散方程. 两层嵌套光滑三角形网格是由三角形背景网格本身以及四个等面积三角形子网格组成. 为了提高方法的精度, 由Richardson外推法确定两层光滑网格上的最优光滑应变. 几个数值算例验证了该方法的精度和计算效率. 数值结果表明, 随着光滑积分网格数目的增加, 光滑无网格伽辽金法的计算精度逐步接近EFGM 的, 但计算效率要远远高于EFGM的. 另外, 光滑无网格伽辽金法的边界条件可以像有限元那样直接施加. 从计算精度和效率综合考虑, 光滑无网格伽辽金法比EFGM具有更好的数值表现, 具有十分广阔的发展空间.      

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

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

利用无网格伽辽金法分析材料稳态蠕变问题

稳态蠕变是高温环境下材料的一个重要的考虑问题，它也是材料破坏的一种重要形式。由于存在划分网格，利用传统有限元法模拟稳态蠕变有一定的不足之处。作为一种新兴的数值模拟方法，无网格法不需要划分单元，只需将求解问题离散成独立的结点，计算过程可以局部细化。利用连续介质的虚功原理，将无网格伽辽金法应用于稳态蠕变的数值模拟，推导了稳态蠕变的无网格伽辽金法控制方程。利用罚参数来实现稳态蠕变的不可压缩条件和本质边界条件，能够保证求解过程中刚度矩阵的对称正定性。通过实例的计算结果表明，无网格伽辽金法在求解稳态蠕变时具有较高的计算精度，结果与理论解结果吻合，而且前后处理较为简单。     

反平面断裂问题的无单元伽辽金比例边界法

将比例边界法与无单元伽辽金法相结合,建立了反平面断裂分析的无单元伽辽金比例边界法。这是一种边界型无网格法,在环向方向上采用无单元伽辽金法进行离散,因此计算时仅需要边界上的节点信息,不需要边界元所要求的基本解。为了便于施加本质边界条件,通过建立节点值和虚拟节点值之间的关系给出了修正的移动最小二乘形函数。在径向方向上,该方法利用解析的方法求解,因此是一种半解析的数值方法。最后,给出了数值算例,并验证了所提方法后处理简单和计算精度高的特点,适合于求解反平面断裂问题。     

基于无网格一图论方法的岩体边坡稳定性分析研究

将无网格伽辽金法应用于岩体边坡稳定性分析,发展了基于无网格模型和有向加权图Bellman—Ford最短路径搜索算法相结合的无网格-图论边坡滑移面搜索方法,以搜寻节理岩体边坡失稳时的临界滑移面并得出其相应的安全系数。区别于传统的滑移面搜索算法,本文方法无需假定滑移面形状,更适用于具有复杂滑移线形状的节理岩体边坡的稳定性分析与计算,具有稳定和高效的算法特点。文中详细论述了无网格-图论最短路径算法的理论、方法和程序实现,并通过算例说明该方法在岩体边坡稳定性分析中的适用性。     

基于无网格—图论方法的岩体边坡稳定性分析研究

将无网格伽辽金法应用于岩体边坡稳定性分析,发展了基于无网格模型和有向加权图Bellman-Ford最短路径搜索算法相结合的无网格-图论边坡滑移面搜索方法,以搜寻节理岩体边坡失稳时的临界滑移面并得出其相应的安全系数.区别于传统的滑移面搜索算法,本文方法无需假定滑移面形状,更适用于具有复杂滑移线形状的节理岩体边坡的稳定性分析与计算,具有稳定和高效的算法特点.文中详细论述了无网格-图论最短路径算法的理论、方法和程序实现,并通过算例说明该方法在岩体边坡稳定性分析中的适用性.     

基于刚（粘）塑性流动理论的自然单元法研究

将自然单元法与刚（粘）塑性流动理论相结合,对自然单元法在金属塑性成形过程数值模拟中的应用进行了研究。采用基于Voronoi图和Delaunay三角化结构的Non-Sibsonian插值方法构造近似速度场向量,实现无网格方法中速度边界条件的直接精确施加,提出了基于刚（粘）塑性流动理论的无网格自然单元法。运用不完全广义变分...     

无网格法研究进展及其应用

从加权残量法的角度出发，系统地总结了现有各种无网格法的基本格式，阐明了无网格法的特点，论述了无网格法的研究进展，给出了无网格法在碰撞、动态裂纹扩展、金属加工成型、流体力学以及其它领域中的应用。     

基于单位分解法的无网格数值流形方法

在数值流形方法和单位分解法的基础上，提出了无网格数值流形方法. 无网格数值流形 方法在分析时采用了双重覆盖系统，即数学覆盖和物理覆盖. 数学覆盖提供的节点形成求解 域的有限覆盖和单位分解函数；而物理覆盖描述问题的几何区域及其域内不连续性. 与原有 的数值流形方法相比，无网格数值流形方法的数学覆盖形状更加灵活，可以用一系列节点的 影响域来建立数学覆盖和单位分解函数，具有无网格方法的特性，从而摆脱了传统的数值流 形方法中网格所带来的困难. 与无网格方法相比，由于采用了有限覆盖技术，试函数的构造 不受域内不连续的影响，克服了原有的无网格方法在处理不连续问题时所遇到的困难. 详细推导了无网格数值流形方法的试函数和求解方程，最后给出了算例，验证了该方法的正 确性.     

无网格MLPG法求解轴对称弹性体扭转问题

应用无网格局部彼得洛夫-伽辽金法(MLPG)研究轴对称弹性体扭转问题,给出了矩阵形式的控制方程,发展了MLPG求解轴对称体弹性扭转问题的数值计算方法。算例分析表明:此方法对求解此类问题具有良好的适应性,数值解能达到理想的计算精度。     

无网格MPM法三维前处理系统设计

无网格法在解决爆炸冲击、裂纹扩展、高速碰撞以及大变形等问题时,由于与网格无关,能够避免有限元法因网格畸变计算精度降低或计算终止的矛盾,是目前求解固体力学问题的有效方法之一,也是近年来研究的热点。前处理对于无网格法数值计算是非常重要的环节,基于面向对象程序设计思想,应用C 和OpenGL开发了三维图形交互方式的MPM法前处理系统,提出了应用递归四叉树法、扫描表示法和CSG法划分任意形状非等间距的MPM法质点单元,为建立通用的无网格前处理平台提出了有效的解决方案。     

Ideal sheet forming with frictional constraints

The ideal forming theory, previously developed as a direct design method to guide iterative design practices based on analytic methods, does not properly account for frictional effects prevalent in real forming. A method for introducing the effects of the Coulomb friction in the ideal sheet forming theory was developed in this work by modifying the extremum work criterion. For numerical implementation, a rigid-plastic finite element method was used based on triangular membrane elements. Computational results were compared with experiments given in NUMISHEET'93 using a planar anisotropic strain rate potential proposed by Barlat et al. (Barlat, F., Chung, K., Richmond, O., 1993. Strain rate potential for metals and its application to minimum plastic work path calculations. Int. J. Plasticity 9, 51.) and the Coulomb friction model.     

HIGH PERFORMANCE SPARSE SOLVER FOR UNSYMMETRICAL LINEAR EQUATIONS WITH OUT-OF-CORE STRATEGIES AND ITS APPLICATION ON MESHLESS METHODS

A new direct method for solving unsymmetrical sparse linear systems(USLS) arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin (MLPG) method need to solve large USLS. The proposed solution method for unsymmetrical case performs factorization processes symmetrically on the upper and lower triangular portion of matrix, which differs from previous work based on general unsymmetrical process, and attains higher performance. It is shown that the solution algorithm for USLS can be simply derived from the existing approaches for the symmetrical case. The new matrix factorization algorithm in our method can be implemented easily by modifying a standard JKI symmetrical matrix factorization code. Multi-blocked out-of-core strategies were also developed to expand the solution scale. The approach convincingly increases the speed of the solution process, which is demonstrated with the numerical tests.     

基于核重构思想的最小二乘配点型无网格方法

介绍重构核点法的基本原理和近似函数的构造方法，并基于核重构思想，应用配点法和最小二乘原理，离散微分方程，建立求解的代数方程，提出了一种基于核重构思想的最小二乘配点型无网格方法．与一般配点法相比，该方法的系数矩阵是有对称正定的，计算精度高，稳定性好．该方法的实施不需要背景网格，不需要进行高斯积分，与Galerkin法相比，具有计算量小、边界条件处理简单的特点，是一种真正的无网格法．对该方法构造过程中的近似函数及其导数的计算、修正函数的计算及方法的实现等问题进行了探讨．文中结合若干典型算例，检验了该方法的有效性．     

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.     

Enriched meshless manifold method for two-dimensional crack modeling

The meshless manifold method is based on the partition of unity method and the finite cover approximation theory which provides a unified framework for solving problems dealing with both continuum with and without discontinuities. The meshless manifold method employs two cover systems. The mathematical cover system provides the nodes for forming finite covers of the solution domain and the partition of unity functions. And the physical cover system describes geometry of the domain and the discontinuous surfaces in the domain. The shape functions are derived by the partition of unity and the finite covers approximation theory. In meshless manifold method, the mathematical finite cover approximation theory is used to model cracks that lead to interior discontinuities in the displacement. Therefore, the discontinuity is treated mathematically instead of empirically by the existing methods. However, one cover of a node is divided into two irregular sub-covers when the meshless manifold method is used to model the discontinuity. As a result, the method sometimes causes numerical errors at the tip of a crack. To improve the precision of the meshless manifold method, the enriched methods are introduced in this work for crack problems.     

基于局部Petrov-Galerkin离散方案的无网格法

基于局部Petrov-Galerkin离散方案,选用自然邻近插值构造试函数,用Shepard函数作为权函数,提出了一种无网格方法(MNNPG),这种方法充分发挥了局部Petrov-Galerkin法的优势,并且结合了自然邻近插值的特点,方便引入边界条件,由于以Shepard函数的圆形支集作为积分子域,用分片中点插值来完成区域积分,无需额外背景网格,是一种真正的无网格法。本文将该无网格方法用于求解二维弹性力学边值问题,算例结果很好地吻合了精确解,表明该方法具有良好的数值精度和稳定性。