基于应变能的自适应EFG方法关键技术的研究.pdfVIP

基于应变能的自适应EFG方法关键技术研究1 张征，刘更，刘天祥 西北工业大学机电学院，西安 （710072 ） E-mail ：zhangzheng_1979@ 摘 要：本文研究了基于应变能的自适应无网格伽辽金方法。给出自适应无网格伽辽金方法 的原理、算法流程和程序实现，并对该方法中的关键技术进行了讨论，包括提出基于背景网 格的误差估计方法、由结点组成的四边形或三角形局部结点加密技术和EFG-FE 耦合方法及 自适应搜索半径算法，计算了带小圆孔平板拉伸和圆柱体与刚性平面的弹性接触的算例。计 算表明，该方法能在较少的计算规模下得到较正确的计算结果。 关键词：自适应；无网格伽辽金方法；无网格伽辽金-有限元耦合方法；应变能 1 自适应无网格方法 有限元方法在处理许多非线性和不连续问题，如冲击裂纹扩展、材料破坏及失效，材料 相变、腐蚀和渗透，以及大变形等遇到了困难。无网格方法(Meshless Methods 或Meshfree Methods)免除了定义在求解区域上的网格结构，不受网格约束，可以方便地在求解域内增加 和减少结点，从而可以极大地改善局部区域内的求解精度。1994 年，Belytschko等人[1]提出 了采用MLS近似函数的新方法，称为无网格伽辽金(Element-Free Galerkin Method, EFGM)方 法。 在无网格方法中，由于无网格方法排点灵活，且在计算中可根据需要插入结点，促成了 自适应算法的发展。其基本思想是初始时将结点排布得较为稀疏，得出初步结果后根据位移 或应变梯度、应变能密度 （Strain Energy Density ，SED）的变化，确定需要做进一步分析以 提高精确度的区域，在这些区域插入较多的结点重新计算，以此类推，直到计算结果满足要 求。此方法既保证了计算的精确度，而且减少了运算量，可看作有限元方法中h方法[2]在无 网格领域的推广。基于无网格伽辽金方法的 自适应方法可称为自适应无网格伽辽金方法 （Adaptive EFGM ，AEFGM ）。 基于 MLS 近似函数，EFG 方法位移近似的基本构造为 m u h (x) ∑p (x)a (x) ≡pT (x)a(x) ( 1) i i i 1 式中p (x)表示任意阶的基函数，m表示基函数的项数，而a(x)表示一个特定的系数矩阵 i a(x) A −1 (x)B(x)u (2) 式中 A PT W(x)P (3) B PT W(x) (4) 式中 P 是由基函数构成的矩阵，W 是由权函数构成的矩阵，分别为: 1 本课题得到 国家自然科学基金项目）；教育部高等学校博士学科点专项科研基金项 目 （20030699035 ）；陕西省自然科学基金项目 （2004E225 ）和西北工业大学 2003 年博士论文创新基金项目 （CX200312 ）的资助。 -1- ⎡p 1 (x1 ) p 2 (x1 ) K p m (x1 ) ⎤ ⎢ ⎥ ⎢p 1 (x2 ) p 2 (x2 ) K p m (x2 ) ⎥