力电耦合扩展无网格伽辽金法在压电柔性臂断裂分析中的研究

为提高含裂纹压电柔性臂在断裂分析中的求解精度,基于压电材料断裂力学理论,建立了压电柔性臂的力学分析模型。将描述裂纹尖端奇异性的位移场函数和电场函数引入到传统无网格伽辽金法中,提出了力电耦合扩展无网格伽辽金法。与传统无网格伽辽金法相比,本方法只需要较小的影响域来描述裂尖奇异场,并且节点影响域不会被裂纹线影响,不要可视准则和衍射准则,提高了计算精度。由虚拟裂纹闭合法推导了压电材料能量释放率,讨论了不同高斯点密度对强计算结果的影响。与解析解、有限元法的计算结果进行了比较,在高斯点个数为6×6时,扩展无网格伽辽金法的计算误差为1.88%,明显好于有限元的计算误差2.5%。数值算例结果表明本方法具有较高的计算精度。     

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

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

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

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

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

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

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

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

无网格法的理论及应用

详细论述了近年来迅速发展的无网格法的理论基础及其在各个领域内的应 用. 无网格法网格依赖性弱, 避免了传统的有限元、边界元等基于网格的数值方法 中可能出现的网格畸变和扭曲, 在一些有限元、边界元等方法难以较好处理的领域体现 出独特的优势. 以加权余量法为主线归纳了已有的30多种无网格法, 各类 无网格法的主要区别在于使用了不同的加权余量法和近似函数. 详尽介绍 了各种无网格近似方案(包括移动最小二乘近似、核近似和重构核近似、单位分 解近似、径向基函数近似、点插值近似、自然邻接点插值近似等)和无网格法 中常用的各类加权余量法(伽辽金格式、配点格式、局部弱形式、加权最小二乘 格式和边界积分格式等), 并讨论了数值积分方法和边界条件的处理等问题. 在 此基础上较系统地总结了无网格法在冲击爆炸、裂纹传播、超大变形、结 构优化、流固耦合、生物力学和微纳米力学等领域的应用, 展示了无网格法相 对于传统数值方法的优势.      

波前法在无网格伽辽金法中的应用

有别于有限元法,无网格法采用基于点的近似,可彻底或部分地去除网格（只保留积分所需的背景网格）,在保证计算精度同时降低计算难度。无网格伽辽金法（Element Free Galerkin method, EFG）是一种基于移 动最小二乘近似（Moving Least-Squares, MLS）的全局弱式无网格法,广泛应用于计算力学等领域,该方法的一个缺点是：计算过程中产生的系数矩阵含有的非零元数量比有限元法多,即使处理中等规模模型时,也要求计算机有很大的存储空间,并且计算时间长。波前法在有限元法中已有很成熟的应用,但至今没有应用于无网格方法。本文介绍了波前法在无网格伽辽金法中的应用方法,编写了相应的计算程序,并以弹性力学为例做了验算。     

基于整体位移模式的平面裂纹结构数值模拟方法

对于平面裂纹问题,针对扩展有限元法和无网格伽辽金法的不足,从结构的整体位移模式出发,提出了一种新的数值模拟方法。在整个求解域内构造其试探函数,并引入裂纹修正项描述裂尖处的奇异性和裂纹面的强间断特性;同时,提出了一种新的强制边界条件施加方法,通过引入位移边界水平集函数,将位移边界条件包含在近似位移场的表达式中,有效地解决了位移边界条件问题,减小了刚度矩阵的阶数,非常方便地消除了刚度矩阵的奇异性,降低了线性方程组的求解难度。含裂纹矩形平板结构的数值算例验证了该方法的有效性。     

几何非线性分析的无网格伽辽金算法

利用几何非线性的应变-位移关系,在小应变假设的条件下,推导出二维几何非线性问题中的无网格伽辽金法的计算格式。由于无网格方法中的形函数不具备Kronecker delta性质,文中采用罚方法来实现本质边界条件。数值实例表明,无网格伽辽金法在处理几何非线性问题时,具有较高的计算精度,是一种有效的数值计算方法。     

复合材料层合板的谱随机无网格伽辽金法

针对基于摄动理论的计算方法不适用于分析随机结构参数大变异问题,提出谱随机无网格伽辽金法,该方法基于随机场正交分解理论,将随机场采用Karhunen-Loève级数展开为一系列不相关随机变量,再引入结构位移随机响应的混沌多项式分解,结合无网格伽辽金法,从而导出含随机变量的复合材料层合板的谱随机无网格伽辽金法,给出结构响应的统计特征值的计算公式,该方法既适用于随机参数大变异情况,又具有无网格法的优势;数值算例结果表明该方法是正确有效的。     

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

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

A 3D shell-like approach using element-free Galerkin method for analysis of thin and thick plate structures

A new efficient meshless method based on the element-free Galerkin method is proposed to analyze the static deformation of thin and thick plate structures in this paper. Using the new 3D shell-like kinematics in analogy to the solid-shell concept of the finite element method, discretization is carried out by the nodes located on the upper and lower surfaces of the structures. The approximation of all unknown field variables is carried out by using the moving least squares (MLS) approximation scheme in the in-plane directions, while the linear interpolation is applied through the thickness direction. Thus, different boundary conditions are defined only using displacements and penalty method is used to enforce the essential boundary conditions. The constrained Galerkin weak form, which incorporates only displacement degrees of freedom (d.o.f.s), is derived. A modified 3D constitutive relationship is adopted in order to avoid or eliminate some self-locking effects. The numeric efficiency of the proposed meshless formulation is illustrated by the numeric examples.     

A high-order Runge-Kutta discontinuous Galerkin method with a subcell limiter on adaptive unstructured grids for two-dimensional compressible inviscid flows

A robust, adaptive unstructured mesh refinement strategy for high-order Runge-Kutta discontinuous Galerkin method is proposed. The present work mainly focuses on accurate capturing of sharp gradient flow features like strong shocks in the simulations of two-dimensional inviscid compressible flows. A posteriori finite volume subcell limiter is employed in the shock-affected cells to control numerical spurious oscillations. An efficient cell-by-cell adaptive mesh refinement is implemented to increase the resolution of our simulations. This strategy enables to capture strong shocks without much numerical dissipation. A wide range of challenging test cases is considered to demonstrate the efficiency of the present adaptive numerical strategy for solving inviscid compressible flow problems having strong shocks.     

基于Delaunay三角化的无网格法计算结果后处理

最近新发展起来的无网格方法由于不需要显式网格，节省了网格生成所需的大量时间，并且避免了网格畸变问题，所以在处理一些特殊问题如移动边界、大变形、高梯度等方面显示出特殊的优越性。但另一方面也使得计算结果的全域后处理遇到困难。提出了一种基于Delaunay三角形背景网格的实用无网格计算结果后处理方法，以无网格离散节点为顶点生成Delaunay三角形，将无网格法计算得到的节点应力值映射插值得到三角形内的应力场云图颜色。给出的二维线弹性应力分析算例表明方法可靠实用。     

大变形问题分析的局部Petrov-Galerkin法

在微机电系统(MEMS)的建模和模拟研究中,大变形或大移动要充分予以考虑.用有限元法分析这类问题,由于难以避免的网格畸变,使模拟效率精度降低甚至失效,无网格方法(Meshless Method)则能在分析这类问题时显示出明显的优势,无网格局部Petrov-Galerkin(MLPG)法被誉为是一种有发展前景的真正无网格法.本文进一步发展了MLPG法,通过对任意的离散分布节点采用局部径向基函数构造插值形函数和Heaviside权函数,分析方程采用局部加权弱形式离散,建立了变量仅依赖于初始构型的完全Lagrange分析格式,最后用Newton-Raphson法迭代求解.文中分析了悬臂梁典型算例和微机电开关非线性大变形问题,通过与有限元结果的比较,表明本文提出的大变形问题无网格局部Petrov-Galerkin法具有稳定性好及收敛性快等优点.     

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

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

平面裂纹问题的h, p, hp型自适应无网格方法的研究

无网格方法以其独特的优点：不需“网格”（即节点间的连接信息）划分,特别适合自适应的分析,在分析中只需要高梯度域简单地插入离散点（ｈ型）或保持模型节点数、分布、覆盖大小均不变,中增加高误差覆盖上的函数的多项式阶次（ｐ型）,便可以得到更高精度的数值模型。针对平面弹性问题发展和推导一种显式后验误差指示公式,对平面裂纹实例进行了ｈ型,ｐ型,ｈｐ型三种不同类型的无网格自适应分析,数值分析结果表明了这种自适应     

Time discontinuous Galerkin finite element method for dynamic analyses in saturated poro-elasto-plastic medium

A time-discontinuous Galerkin finite element method for dynamic analyses in saturated poro-elasto-plastic medium is proposed. As compared with the existing discontinuous Galerkin finite element methods, the distinct feature of the proposed method is that the continuity of the displacement vector at each discrete time instant is automatically ensured, whereas the discontinuity of the velocity vector at the discrete time levels still remains. The computational cost is then obviously reduced, particularly, for material non-linear problems. Both the implicit and explicit algorithms to solve the derived formulations for material non-linear problems are developed. Numerical results show a good performance of the present method in eliminating spurious numerical oscillations and providing with much more accurate solutions over the traditional Galerkin finite element method using the Newmark algorithm in the time domain. The project supported by the National Natural Science Foundation of China (19832010, 50278012, 10272027) and the National Key Basic Research and Development Program (973 Program, 2002CB412709)     

基于不连续场修正权的无网格断裂分析

本文采用了一种基于不连续场修正权函数的无网格方法来处理二维平面问题中的有限长裂纹。相较于目前常用的无网格裂纹不连续性处理方案，采用修正权函数处理裂纹附近不连续场时只需要对原权函数进行修正，算法简便易实现。本文采用基于不连续场修正权函数的无单元Galerkin方法（EFGM），对在边界上施加I-II混合型裂纹位移场的斜裂纹板进行了数值分析。并与可视性准则、衍射法和透射法等不连续准则对比了裂尖位移场、应力场和应力强度因子解的数值精度。另外，本文还对这四种不连续准则形函数的计算效率进行了分析和比较。