一种改进的无单元方法

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

无单元法（EFGM）——在岩土工程上有限元法的有力补充

作为一种新的岩土工程数值计算方法，无网格伽辽金法具有局部化技术的无单元特性。该法与传统FEM相比，在精度和后处理等方面有着明显的计算优势。本文全面介绍了EFGM，并结合岩土工程领域中的土体固结变形问题，推导了EFGM的一维、二维固结问题的离散方程；同时对算例进行计算，将其结果与精确解和传统FEM解进行比较后，得出结论：精度是较高的，处理边界是准确的。尤其当与FEM联合使用时，既能发挥FEM的经验优质，又能体现EFGM的无网格特点，计算效果十分显著，因此，EFGM是FEM的有力补充。     

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

计算效率低的问题长期阻碍着无网格伽辽金法(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具有更好的数值表现, 具有十分广阔的发展空间.      

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

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

无网格Galerkin法与非协调元方法的耦合

无网格Galerkin法(EFGM)处理不可压缩问题时不存在自锁现象，有限元方法(FEM)也常被用来与其耦合以方便地施加边界条件和提高计算效率。在有限元方法中使用等参元，EFGM与FEM的耦合方法在处理不可压缩问题时仍然存在自锁现象。本文在有限元方法中，采用非协调元，将无网格kGalerkin法与非协调元耦合，保留了耦合方法的优点，且避免了求解不可压缩问题时的自锁现象。算例显示本文方法在分析平面应变不可压缩问题时能得到合理的结果。     

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

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

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

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

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

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

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

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

一种新的数值方法——无网格伽辽金法（EFGM）

无网格伽辽金法（ＥＦＧＭ）是近几年发展起来的与有限元相似的一种数值算法,它采用移动的最小二乘法构造形函数,从能量泛函的弱变分形式中得到控制方程,并用拉氏乘子满足本征边界条件,从而得到偏微分方程的数值解中得到该法只需节点信息,不需将节点连成单元,此外,还有精度高,后处理方便等优点,本文介绍其基本原理及实现过程,并用算例表明,该法具有一定的发展前景。     

TOPOLOGY SYNTHESIS OF GEOMETRICALLY NONLINEAR COMPLIANT MECHANISMS USING MESHLESS METHODS

This paper presents a new method for topology optimization of geometrical nonlinear compliant mechanisms using the element-free Galerkin method （EFGM）. The EFGM is employed as an alternative scheme to numerically solve the state equations by fully taking advantage of its capability in dealing with large displacement problems. In the meshless method, the imposition of essential boundary conditions is also addressed. The popularly studied solid isotropic material with the penalization （SIMP） scheme is used to represent the nonlinear dependence between material properties and regularized discrete densities. The output displacement is regarded as the objective function and the adjoint method is applied to finding the sensitivity of the design functions. As a result, the optimization of compliant mechanisms is mathematically established as a nonlinear programming problem, to which the method of moving asymptotes （MMA） belonging to the sequential convex programming can be applied. The availability of the present method is finally demonstrated with several widely investigated numerical examples.     

基于无网格Galerkin方法的整体式柔性机构的多准则拓扑优化设计

基于一种新型的数值计算技术—无网格伽辽金法,提出了一种整体式柔性机构拓扑优化设计的新方法.利用移动最小二乘形函数和伽辽金的弱变分形式建立弹性问题的控制方程,用Lagrange乘子法增强本质边界条件.在优化问题中,同时综合考虑机构的柔性和结构的刚度要求,用折衷规划法建立了柔性机构拓扑优化的多准则优化模型,这样对于非凸的优化问题也能保证搜索到Pareto解集所有的有效解.基于SIMP密度函数惩罚模型和优化准则法,建立了一种设计变量的显式迭代格式.运用经典算例证明了文中方法的正确性和有效性.     

无网格方法在平面粘弹性力学问题中的应用

本文综合应用无网格方法(EFGM)、线性粘弹性与弹性力学之间的对应原理,Laplace变换和逆变换等方法求解了拟静态平面弹性和粘弹性力学问题。首先,利用Laplace变换和逆变换推导了平面问题的粘弹性本构关系,建立了拟静态粘弹性平面问题的边值问题;其次,利用粘弹性与弹性力学之间的对应原理得到了Laplace变换域中平面问题的基本方程,在Laplace变换域中建立了相应的泛函,并得到了用无网格方法离散的控制方程;同时,求解了几个拟静态弹性和粘弹性平面问题,给出了它们的表达式和数值结果;最后,采用Laplace逆变换和数值逆变换,得到了粘弹性力学平面问题在物理空间中的解,并比较了由解析解和无网格数值方法所得到的数值结果,可以看到它们是非常吻合的。说明本文方法的正确性和有效性。     

RESEARCH ON THE COMPANION SOLUTION FOR A THINPLATE IN THE MESHLESS LOCAL BOUNDARYINTEGRAL EQUATION METHOD

The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method ( GFEM ), boundary element method (BEM) and element free Galerkin method (EFGM), and is a truly meshless method possessing wide prospects in engineeringapplications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.     

无单元法求解任意边界条件下的中厚板弯曲问题

本文用无单元法进行不同边界条件下的中厚板弯曲问题的求解,提出了构造其近似位移函数的三种形式的权函数,从变分原理出发导出了Mindlin-Reissner中厚板弯曲问题的控制方程,并编制了相应的计算程序。数值算例表明,无单元法用于中厚板弯曲问题是合理可行的,其结果具有相当高的精度。