《有限元分析第五章(第一部分).doc》由会员分享,可在线阅读,更多相关《有限元分析第五章(第一部分).doc(20页珍藏版)》请在金锄头文库上搜索。
1、第五章 等参数单元(Isoparametric Elements)在前面的章节中我们已经认识了三角形单元和矩形单元。这两种单元的边均为直边,用直边单元离散曲边的求解域势必要用更多的单元数才能较准确地描述实际边界。本章将要介绍的等参数单元是目前应用最广的一类单元,可用这类单元更精确的描述不规则的边界。这类单元的出现不仅系统的解决了构造协调位移单元的问题,而且自然坐标系的描述方法也广泛为其他类型的单元所采用。等参数单元在构造形函数时首先定义一个规则的母体单元(参考单元),在母体单元上构造形函数,再通过等参数变换将实际单元与母体单元联系起来。变换涉及两个方面:几何图形的变换(坐标变换)和位移场函数的
2、变换,由于两种变换采用了相同的函数关系(形函数)和同一组结点参数,故称其为等参数变换。5-1四结点四边形等参数单元、 母体单元自然坐标和形函数(1,1)(-1,-1)342图5-1母体单元:边长为的正方形,自然坐标系,示于图5-1。取四个角点为结点,在单元内的排序为、。仿照矩形单元,可定义出四个形函数(5-1-1)显然有如下特点:(i)是,的双线性函数(ii)(iii)(5-1-2)、 实际单元与母体单元之间的坐标变换() 坐标变换设xy平面上的实际单元e由母体单元经过变换F得到,即 且规定结点(i,i)与结点(xi, yi)对应(i=14)。这样的变换不只一个,利用(5-1-1)定义的形函数
3、即可写出这种变换中的一个(5-1-3)x,uy,v3241=-1=1=-1/2=0=1/2=1=1/2=0=-1/2=-1图5-2(5-1-3)所定义的变换有如下特点:x, y是,的双线性函数。沿母体单元中常数的直线(坐标线),x, y是的线性函数,对应于单元e中的一组直线,特别,单元e的一组对边1-2、3-4为直线。类似,中常数的另一组坐标线对应于单元e中的另一组直线。特别,e的另一组对边2-3、4-1也是直线,单元e为直边四边形。单元的其他直线(例如对角线1-3),变换到单元e中将是一条曲线(图5-2)()Jacobi矩阵Jacobi行列式矩阵(5-1-4)称为变换的Jacobi矩阵。de
4、tJ称为变换的Jacobi行列式。一般情况下,J的元素和detJ都是,的函数。若detJ恒不为零(一般使它恒正),则J-1存在,变换F存在逆变换F-1。使单元e内的任一点(x, y)对应于单元内的一确定点(,)。此时称变换F为非奇异的。detJ称为变换特征量。detJ还具有明显的几何意义,如图5-3所示。设在(,)处detJ0在(,)附近取一边长为d,d的长方形。设此长方形与单元e内的一个小子区域d对应,可以证明,此小子域的面积d在略去高阶微量后有ed(,)dd图5-3(x,y)例图5-4所示的实际单元e为边长分别为2a、2b的矩形。结点坐标为:2a 2b 1 (c,d)4 2 3 xy 图5
5、-4则由(5-1-3),可得出坐标变换为同样得到:表明:当实际单元e为矩形时,经坐标变换得到的x, y是,的线性函数。Jacobi矩阵Jacobi行列式在单元内是常数。当结点序号按图5-4的转向排列时,detJ恒正。、 单元内假设的位移场对于平面问题,设沿总体坐标系的位移为u、v,结点(xi, yi)的位移为ui,vi实际单元e内的假设位移场(Trial function)取为(5-1-5)注意,这里u、v虽然是用点的自然坐标,表述的,但位移u、v (以及后面的单元刚度矩阵)却是对总体坐标系的。这与第二章中在单元局部坐标系下定义位移场的作法有区别。在坐标变换(5-1-3)和假定的位移场(5-1
6、-5)中使用的是同一套变换关系(形函数),同一套变换参数(与(xi, yi)对应的结点位移(ui,vi))满足这一特征的单元称为等参数单元。这样定义单元有不少优点,但也对我们提出了一些新问题。假定的位移场是,的双线性函数,当实际单元为矩形时,可表示成x, y的线性函数,假定的位移场u、v是x,y的多项式。但对一般单元而言,不能表示成x,y的多项式,因而位移场u、v不再是x,y的多项式,不能直接利用第四章的结果进行收敛性分析。、 收敛性分析() 单元内位移场连续x、y、u、v都是,的双线性函数(连续函数)。只要Jacobi行列式detJ0,u、v就是x,y的连续函数。即在实际单元内u、v连续。(
7、) 刚体位移和常应变条件对于二阶问题,这个条件归结为假定的位移场中包括总体坐标的完全一次多项式。或者换一个提法:假定的位移场可以精确地表述任何一种线性变化的真实位移场。当试探函数直接用总体坐标的多项式描述时(像第四章所做的那样)采用前面一种提法是方便的。现在试探函数是用自然坐标表述的,则用后一种提法更合适一些。我们定义的形函数满足:(5-1-6)设真实位移场为x,y的线性函数将x,y按(5-1-3)代入,并利用 有注意到 (结点位移的真实值)y,vx,ueeMs图5-54则有类似有上述论证表明:只要所定义的形函数满足(5-1-6)(不管形函数的具体表达式如何),且坐标变换和假定的位移场使用同一
8、组形函数(等参数单元总是如此),那么这样假设的位移场一定能够精确地表述任何一种线性位移场,即刚体位移和常应变条件总可以得到满足。() 协调性对于二阶问题要求穿过单元边界时位移连续。如图5-5所示,考虑一个实际单元e,它的母体单元为。以1-2边为例。沿1-2边常数,x、y、u、v都是的线性函数。设e边界上的M点与边界上的点对应,则M到结点的距离将是的线性函数。反过来也是的线性函数,因而u,v也是的线性函数,完全由这个边界上两个结点1、2的位移值u1、u2、v1、v2所决定。从另一相邻单元 e 看来,沿边1-2, u、v也是的线性函数。完全被结点、处的位移值所决定。从单元e和e 看来沿共同边界1-
9、2上的位移处处相同,即在边界上位移是连续的。对其他边界可用类似的方法加以证明。四结点四边形等参元的形状有较大灵活性,巧妙地解决了单元形状的灵活性和收敛条件(主要是协调条件)之间的矛盾。但是一般的四边形单元只能精确地再现线性变化的位移场,有限元空间Sh的次数k1=1。虽然能保证有限元解的收敛性,但精度不够满意。当实际单元是矩形时,是x、y的线性函数,假定的位移场将是x、y的二次多项式,但只完全到一次多项式,二次项不完全。这不完全的二次项有时可能改善精度,有时则不能。例如,在分析图5-6的“纯弯曲”应力场时,图(a)中的单元将比图(b)中的单元效果好,尽管还不能说满意。提高单元精度的一个途径是增加
10、结点个数,提高插值函数阶次。图5-6(b)(a)5、四结点单元的应用实例及相关限制条件 某求解域如图5-7(a)所示,若将该区域用个四结点等参元进行离散,母体单元如图5-7(b)所示。xy0.01124342332.03.05.05.03.02.00.0142y1432x(c)(a)12341(d)(b)图5-7从图中可以看出:、2号单元与母体单元的结点编号顺序一致,均为逆钟向,而号单元的编号顺序为顺钟向;、号单元为凸形单元,即连接任意两点结点的线段均在单元内部,而单元是非凸形单元,如连接结点、的线段不在单元内。下面讨论这些差别在母体单元与实际单元进行映射时的影响。在母体单元中形函数如式(5-
11、1-1),坐标变换关系如式(5-1-3)。首先,计算出Jacobi矩阵中的各元素如下下面计算出各单元具体的变换关系及Jacobi行列式的值单元:各结点的坐标为,则Jacobi行列式是的线性函数,对所有的值()Jacobi行列式的值恒为正,因此,母体单元与单元的变换是可逆的。单元:各结点的坐标为,则 Jacobi行列式的值沿着直线为零,母体单元中的阴影部分将映射到实际单元的阴影部分,这部分显然在实际单元之外。例如,母体单元中的点落在阴影部分,该点映射到了实际单元的。因此,母体单元与单元的变换不是可逆的。所以内角大于网格在任何单元中都是不允许的。一般来说,有限元网格中内角过大或过小都是不合适的。单
12、元:各结点的坐标为。则Jacobi行列式的值小于零表示:右手坐标系映射到左手坐标系,这种变换关系在有限元方法中也是不允许的。若将单元3的结点编号顺序改为逆钟向,即:等参变换成为可逆变换。5-2八结点四边形(直边或曲边)等参数单元8342图5-765、 母体单元形函数母体单元仍为边长为的正方形(图5-8)。自然坐标系,如图所示。在单元中配置八个结点,其中仍位于角点上,则位于各边中点。构造出八个形函数N1N8。如下:(5-2-1)验证可知,具备以下性质(5-2-2)构造出(5-2-1)的八个形函数是不太容易的,下一节将介绍构造这类形函数的一般方法。y76543218x0图5-、 实际单元和坐标变换
13、实际单元e的八个结点坐标为(xi, yi)(图5-9),则母体单元到实际单元e的坐标变换取为(5-2-3)Jacobi矩阵当Jacobi行列式 detJ0时,(5-2-3)规定的变换是非奇异的。(5-2-1)所定义的形函数对于变量或变量来说次数都不超过。沿母体单元中常数的直线(坐标线),根据(5-2-3),x、y将是的二次函数,因而对应于实际单元e中的一族曲线(图5-3)。母体单元中常数的直线将对应于实际单元e中的另一族曲线。在一般情况下单元e将是曲边四边形。当实际单元e为矩形,且结点位于各边中点时,变换(5-2-3)的右端退化为、的一次多项式,反过来、也可表示为x、y的线性函数。、 单元内假设的位移场单元内的位移场(即试探函数)采用与坐标变换相同的一套形函数(5-2-1),对平面问题则认为单元e内有(5-2-4