广义扩展有限元法的研究.docxVIP

广义扩展有限元法的研究 1 广义扩展有限元法 自出版以来，xferm法被广泛用于制定传统的有限元法框架（gfm）和xferm法的两个新方法，自成立以来，它引起了学术界的广泛兴趣。广义有限元法将传统有限元结点自由度广义化,通过增加插值函数的阶次,有效提高了计算结果的精度,且在非线性分析方面较传统有限元法更为可靠。广义有限元法的概念由Babushka首次提出,随后Babuska、Strouboulis、梁国平和栾茂田等对广义有限元法进行系统研究,取得了许多重要成果。 扩展有限元法的优势则在于分析不连续问题,1999年,Belytschko教授在研究常规有限元框架下富集裂尖渐进位移场形函数的主要项时,首次提出了扩展有限单元法的思想。该方法以单位分解法为基础,将非连续结构直接嵌入单元内部,通过增加描述非连续性的附加函数,模拟裂纹的不连续性及裂纹尖端的奇异性,避免了传统有限元法分析不连续问题时需要重新剖分网格的缺点。扩展有限元法问世后在国际上引起了极大关注,得到了快速发展和广泛的应用。 在裂纹扩展分析中,一般是依据断裂力学的相关准则判断裂纹是否稳定,解的精度往往成为数值模拟的关键;与传统有限元法相比,扩展有限元法本身对解的精度并没有提高。本文结合上述两种方法各自的特点,提出了一种新的数值方法———广义扩展有限元法(GXFEM)。文中对扩展有限元法中的非连续附加位移函数进行了改进,阐述了广义扩展有限元法的基本原理,推导了相关公式,编写了广义扩展有限元法程序,应用所提出的广义扩展有限元法对裂纹扩展过程进行了计算模拟,有效提高了数值解的精度。 2 广义差分和扩展差分 2.1 广义有限元分析 传统有限元中位移插值基函数取为Φ={φ1,φ2,…,φn},单元结点位移向量δe=[u1,u2,…,un]T,则单元任一点的位移可表示为 而广义有限元将传统有限元的结点位移广义化,认为各结点位移具有任意多个广义位移的展开式,即将结点位移表示为广义自由度di的函数,即 式中de=[d1,d2,…,dn]T,F是由结点位移插值函数组成的矩阵。 将式(2)代入式(1)得到新的插值逼近空间,即 式中di=[di,1,di,2,…,di,kmi]T为结点的广义位移,mi为结点位移展开式的项数,k为维数,Ni称为广义形函数。其与传统形函数之间的关系为 式中F(i)为结点位移的插值函数,可任意选取,一般取为坐标的多项式函数。 以平面四结点等参单元为例,若结点位移插值函数取为零阶,有 则每个结点具有两个自由度,即退化为传统有限元法。 若结点位移插值函数取为一阶,有 则每个结点具有6个自由度。 若结点位移插值函数取为二阶,有 则每个结点具有12个自由度。 有了单元的位移模式,便可进行单元分析,类似传统有限元一样,形成广义单元劲度矩阵和广义单元荷载列阵,再进行整体分析,形成整体结构的广义支配方程: 式中K为广义整体劲度矩阵,R为广义整体荷载列阵,D为所求各结点的广义位移。 因为广义形函数中包含传统有限元中的形函数,故该方法的完备性与协调性自然满足。并且由于广义有限元法不同结点位移的多项式之间的协调性没有要求,因此可以按照实际问题的需要去选取多项式的阶次,具有很强的灵活性和适应性。 2.2 附加非连续函数 式中ui为节点i的自由度,Ni为传统有限元单元中节点i的形函数,ai是节点i的附加自由度,φi(x)为附加非连续函数,是考虑非连续性结构的影响而给节点i的丰富函数。 在处理裂纹扩展问题时,附加非连续函数的构造是其中的一个重要环节。裂纹在有限元网格中有两种情况:一种是把单元截成两块,另一种则是裂纹尖端终止于单元内部如图1所示。 对于横跨单元的裂纹,构造的附加非连续函数必须体现位移在裂纹处的不连续性。此外,单元中的点距离裂纹越远,裂纹对该点的影响就越弱。Belytschko教授在提出扩展有限元法时,采用的是单位阶跃函数,难以反映裂纹影响的衰减效应。本文将指数间断函数作为附加非连续函数,先定义最短距离d(x)为 式中f(x)=0为裂纹的曲线方程,然后构造的附加非连续函数φi(x)如下。 该指数间断函数的曲线如图2所示。 对于尖端终止于单元内部的裂纹,构造的附加非连续函数除了要求能反映裂纹处的位移不连续性外,还必须能体现裂纹尖端渐进位移场的奇异性。无论是哪一种裂纹形式(Ⅰ,Ⅱ,Ⅲ),裂纹尖端位移场均可由下面几个基函数组成的函数形式来表达,即 根据上述位移模式,便可类似传统有限元法,进行单元分析和整体分析,建立扩展有限元的支配方程,求出主变量u和附加变量a。 3 广泛扩展限制法 3.1 不连续变量 构造广义扩展有限元的位移模式为 与扩展有限元相比,此处仅对形函数作了变动,其中的Niu和Nia分别为主变量的广义形函数及附加不连续变量的广义形函数。 式中Fu(i)和Fa(i)分别为主变量