首页 理论教育无网格计算方法:伽辽金法与有限元法的比较

无网格计算方法:伽辽金法与有限元法的比较

【摘要】:鉴于这种缺陷,近几年来国际上许多著名的计算力学学者,如T.Belytschko、O.C.Zienkiewicz、S.N.Atluri、J.T.Oden、W.K.Liu等都对无网格方法表现出了极大的兴趣,并进行了大量的研究工作。无网格伽辽金法与有限元法的相似之处在于,两者都是将边值问题等价为一个条件变分问题,然后由条件变分问题通过数值积分离散为代数方程组。

有限元法已成为计算力学中解决工程问题的主要数值手段,然而随着其应用范围的扩展,其固有的一些缺陷也日益突出。在高速碰撞、流固耦合等涉及特大变形的领域中,基于拉格朗日法的有限元网格可能产生严重的扭曲,甚至使得单元的雅可比行列式为负值,不仅在计算中需要网格重构,而且严重地影响解的精度。由于有限元近似基于网格,因此必然难于处理与原始网格线不一致的不连续性和变形。网格重构不仅计算费用昂贵,而且会损害计算精度。鉴于这种缺陷,近几年来国际上许多著名的计算力学学者,如T.Belytschko、O.C.Zienkiewicz、S.N.Atluri、J.T.Oden、W.K.Liu等都对无网格方法表现出了极大的兴趣,并进行了大量的研究工作。无网格方法采用基于点的近似,可以彻底或部分地消除网格,不需要网格的初始划分和重构,不仅可以保证计算的精度,而且可以大大减小计算的难度。然而,由于目前的无网格近似一般没有解析表达式,且大都基于伽辽金原理,因此计算量很大,要超出传统的有限元法;另外,无网格近似大都是拟合,因此对于位移边界的处理比较困难,多采用拉格朗日乘子法处理。

目前已提出了十余种无网格法,其主要区别在于离散微分方程的方法(如伽辽金法、配点法、最小二乘法、彼得洛夫—伽辽金法等)和建立近似函数的方法(移动最小二乘近似、核近似、重构核质点近似、单位分解法、hp云团法、径向基函数法、点插值法等)。

无网格伽辽金法(Element Free Galerkin Method,EFGM)是近20年发展起来的,其理论基础是移动最小二乘法,其基本思想是将计算场域离散成若干个点,由移动最小二乘法来构造形函数,再获得问题的弱变分形式,通过施加边界条件,从而得到数值解。在EFGM中,只需要知道整个分析区域的节点分布以及边界条件,而并不要求节点之间有任何连接。因此,这种算法非常适合在计算过程中场域几何形状和参数发生变化这类问题的求解。无网格伽辽金法与有限元法的相似之处在于,两者都是将边值问题等价为一个条件变分问题,然后由条件变分问题通过数值积分离散为代数方程组。不同之处在于,有限元法是对逐个有限单元进行数值积分,形成单元矩阵,然后将其叠加到单元节点所对应的方程中;而无网格伽辽金法是在积分单元上进行数值积分,然后将每个高斯点上的积分值叠加到该高斯点所支撑的若干节点所对应的方程中。