无网格伽辽金法中本质边界问题处理.pdfVIP

______________________________________________________________________________ 无网格伽辽金法中本质边界问题的处理 1 2 3 王琥 李光耀 钟志华 1,2,3 湖南大学机械与汽车工程学院, 长沙 410082 Email: wanghuenying@ 摘 要： 无网格伽辽金法(EFG)是近几年应用较广泛的一种无网格数值算法。它采用移动 的最小二乘法构造形函数,从能量泛函的弱变分形式中得到控制方程,从而得到偏微分方程 的数值解。但是由于 EFG 方法是采用移动最小二乘法对形函数进行逼近，不满足 Kronecker δ条件，所以本质边界条件不能直接作为约束条件。为此，本文介绍对本质边界进行处理几 种的几种方法，并且编制了相应的程序对其效果进行了验证。 关键词： 无网格;移动最小二乘法;有限元 1. 引言 有限元、边界元等传统的计算方法对于大变形的冲压成型问题、裂纹扩展问题的求解有其局 限性，主要原因是网格的存在使处理与原始网格线不一致的非连续性以及大变形，这些基于 网格的计算方法，在处理随时间变化的非连续性和大变形问题时，需要进行的网格的重新构 形。导致计算的前处理困难，计算结果通常应力、应变不连续。为了减少或消除单元划分给 数值分析带来的不便，近二十年来出现了一些只需节点离散数据，而不需要单元信息的无网 格方法。无网格的种类十分繁多，历史最长的无网格方法是光滑粒子法 [1] （SPH, Smooth Particle Hydrodynamic Method ），最初被用于无边界的天体问题, 采用核函数对离散点进行 逼近或插值。SPH 对于边界点，以及不规则离散点的情况，不能精确的进行求解，为此 Liu [2 ] 提出了对核函数的修正方案，称为 RKPM （Reproducing Kernel Particle Method ），该方法是 SPH 方法的加强完备形式。目前应用最多的构造无网格方法的途径是采用移动最小二乘法 （Moving Least Square Approximation Method ，简称MLS ）进行逼近。Nayroles [3] 最早将 MLS 用于 Galerkin 方法，称之为扩散单元法（Diffuse Element Method,简称 DEM）。Belytschko [4 ] __________________________________ 本课题得到高等学校博士学科点专项科研基金资助课题（20020532021）资助 ______________________________________________________________________________中国科技论文在线 提出无单元的 Galerkin 方法（Element Free Galerkin 简称 EFG ）。Atluri [5] 提出局部 Petrov-Galerkin 方法（Meshless Local Petrov-Galerkin 简称MLPG）。此外， Duarte 和 Oden [6] 提出了单位分解法（Partition of Unity Finite Element Method 简称 PUFEM 法），并认识到MLS 实际上是 PU的特例。此外，近年来又有多种无网格方法问世，如 RFB,PIM 方法等。 由于采用移动最小二乘对形函数进行构造，所以不满足Kronecker δ 条件，所以本质边 界条件不能直接作为泛函的约束条件。这成为无网格方法在实际运用中最棘手的问题之一。 由此出现了有限单元耦合法，罚函数法以及拉格朗日乘子法用来解决无网格方法中的本质边 界条件问题。 2.