有限单元法的进一步讨论.ppt

4-1 矩形单元与插值函数 4-2 高次单元与插值函数(6结点三角形单元)4-3 等参单元和单元坐标变换4-4 数值积分方案 4-5 有限元解的收敛性4-6 轴对称问题有限元及空间问题有限元,4-1 矩形单元与插值函数,平面问题的有限元分析中，矩形单元也是常用的单元之一。,它采用了比常应变三角形单元更高次的位移模式，故可更好地反映弹性体中的位移状态与应力状态。,单元节点位移向量和单元位移向量分别为：,，,，,单元位移函数：,将4个节点的坐标值 代入上式，可解出上式的8个系数。,于是，单元位移可通过单元节点位移表示。,矩阵形式为：,其中,-,为型函数。,，,其中,它们的几何意义为,全局坐标,局部坐标,因此有：,，,单元坐标变换,此单元位移模式是具有完全一次式的非完全二次式。因此单元内位移沿x y方向都是线性变化的，但沿其他方向则不是线性的。,单元每边界上有两个节点，按线性变化的位移形式，单元变形后此两节点确定了惟一的单元直线边界，因而此单元边界间的位移是连续的，,所以矩形单元的这种位移模式又称为双线性位移模式。,正是由于在单元的边界上位移按线性变化，且相邻单元公共结点上有共同的结点位移，因此保证了两个相邻单元在其公共边界上位移的连续性。,这种单元的位移模式是完备的和协调的，满足解的收敛条件，因此4结点矩形元是协调元。,矩形单元中的应力分量都不是常量。正应力分量x的主项(即不与相乘的项)沿y方向按线性变化，而它的次要项(与相乘的项)沿x方向按线性变化。正应力y与此相反。剪应力分量则沿x y两个方向都呈线性变化。这种一个方向为常量，另一方向呈线性变化的情况通常并不能提高单元的精度。,矩形单元给出某种线性应力场，解题时可能给出比常应力三角形单元更好的结果，其数据准备也更方便。,但是，三角形单元的网格划分比较灵活，适应性更强，疏密网格的变化较为方便。,解决上述问题的方法之一是同时采用矩形单元和三角形单元混合使用，如图所示的计算模型，在直边区改用矩形单元,而在曲边附近采用三角形单元，以达到两者兼顾的效果。,矩形单元明显的缺点是不能很好地符合曲线边界，包括与坐标轴不平行的直线边界，因此直接应用受到限制。,更为一般的方法是通过等参变换将局部坐标内的矩形单元变换为总体坐标内的任意四边形(包括曲边四边形)的单元,4-2 6结点三角形单元,单元的6个结点是三角形的3个角点及3个边的中点。,六个节点，12个节点位移能确定12个多项式的系数。,4-2 6结点三角形单元,位移函数取完全的二次多项式。,它采用了比常应变三角形单元更高次的位移模式，故可更好地反映弹性体中的位移状态与应力状态：,单元边界上位移按二次抛物线分布。这种单元的应变在两个坐标方向上都呈线性变化。应力也呈线性变化。因此，单元精度较3结点三角形单元要高。,位移函数中的常数项和完整的一次式满足收敛条件中对完备性的要求。,六个节点，12个节点位移能确定12个多项式的系数。,如角点1，可利用2-3边的边方程以及过边中点4-6的直线方程，并使在结点1上N1=1,即,用面积坐标给出的6节点三角形单元的插值函数：,对于每个结点i，可以选择二条直线，这二条直线通过除结点i以外的所有结点，利用直线方程的左部作为线性函数来构造插值函数。,N1=L1(2L11),N2=L2(2L21),N3=L3(2L31),同理：,三个角点的插值函数可统一写作,Ni=Li(2Li1),三个边中点，可选择非此边中点所在的其他二条边作为直线方程，则可得,N4=4L1 L2,N5=4L2 L3,N6=4L1 L3,4-3 等参单元,对于一个给定问题的求解域，预期用较少的单元即可获得需要精度的解答。,但是用较少的形状规则的单元离散几何形状复杂的求解域常会遇到困难，因此需要寻找适当的方法将规则形状的单元转化为能够适应边界的单元。,4.3.1 等参变换的概念,4-3 等参单元,在有限单元法中最普遍采用的变换方法是等参变换，即单元几何形状的变换和单元内的场(位移)函数的变换采用相同数目的结点参数及相同的插值函数进行。,借助于等参单元可以对于一般的任意几何形状的工程问题方便地进行有限元离散，,等参单元的提出为有限单元法成为现代工程实际领域最有效的数值分析方法迈出了极为重要的一步。,采用等参变换的单元称之为等参单元。,4.3.1 等参变换的概念,为将局部坐标中几何形状规则的单元转换成总体(笛卡尔)坐标中几何形状扭曲的单元，以满足对一般形状求解域进行离散化的需要，要建立一个坐标变换：,为建立上述的变换最方便的方法是将上式也表示成插值函数的形式：,可以看到坐标变换关系式和场(位移)函数的插值表示式在形式上是相同的。,如果坐标变换和函数插值采用相同的结点m=n，并且采用相同的插值函数 Ni=Ni，则称这种变换为等参变换。,如果坐标变换结点数多于函数插值的结点数则称为超参变换。m大于n,反之，如果坐标变换结点数多于函数插值的结点数，则称为次(亚)参变换。m小于n,场函数是用局部坐标(,)表述,坐标变换：将局部坐标中几何形状规则的单元转换成总体(笛卡尔)坐标中几何形状扭曲的单元,局部坐标中几何形状规则的单元,总体(笛卡尔)坐标中几何形状扭曲的单元,n=4,4.3.2 两个坐标系内导数、面积微元间的变换关系,4.3.2 两个坐标系内导数、面积微元间的变换关系,为此需要建立两个坐标系内导数、体积微元、面积微元之间的变换关系，,在局部坐标内的积分限(-1,1)是规格化的，,因此我们希望能在局部坐标内并按标准化、规格化的数值积分方法对单元单刚进行积分。,局部坐标,4.3.2 两个坐标系内导数、面积微元间的变换关系,1.导数之间的变换,J 称为Jacobi矩阵，,这样一来，Ni对于x，y的偏导数可用节点坐标显式地表示为,J可以显式地表示为节点坐标的函数,Ni对于x，y的偏导数？,2.面积微元的变换,3.线段微元的变换,4.3.2 两个坐标系内导数、面积微元间的变换关系,4.单元刚度矩阵的积分变换,4.3.2 两个坐标系内导数、面积微元间的变换关系,4.3.3 等参变换的条件,从微积分学知识已知，两个坐标之间一对一变换的条件是Jacobi行列式|J|不得为0,等参变换作为一种坐标变换也必须服从此条件。,如|J|=0,|J|-1将不成立，所以两个坐标之间偏导数的变换就不可能实现。,如|J|=0，则表明笛卡儿坐标中体积微元(或面积微元)为0,即在自然坐标中的体积微元或面积微元对应笛卡儿坐标中的一个点，这种变换显然不是一一对应的。,4.4 数值积分方法,一维高斯积分：,二维高斯积分：,4.5 有限元的收敛性准则,在有限单元法中，场(位移)函数的总体泛函是由单元泛函集成的。如果采用完全多项式作为单元的插值函数(即试探函数)，则有限元解在一个有限尺寸的单元内可以精确地和真正解一致。,但是实际上有限元的试探函数只能取有限项多项式，因此有限元解只能是真正解的一个近似解答。,我们需要研究在什么条件下，当单元尺寸趋于零时，有限元解趋于真正解。,4.5 有限元的收敛性准则,如果出现在泛函中场函数的最高阶导数是m阶，则有限元解收敛的条件之一是单元内场函数的试探函数至少是m次完全多项式。或者说试探函数中必须包括本身和直至m阶导数为常数的项。,收敛准则1：完备性要求。,单元的插值函数满足上述要求时，我们称单元是完备的。,或者说试探函数中必须包括本身和直至m阶导数为常数的项。,4.5 有限元的收敛性准则,收敛准则2：协调性要求。,当单元的插值函数满足上述要求时，我们称单元是协调的。,简单地说，当选取的单元既完备又协调时，有限元解是收敛的，即当单元尺寸趋于零时，有限元解趋于真正解。,我们称这种单元为协调单元。,如果出现在泛函中的最高阶导数是m阶，则试探函数在单元交界面上必须具有Cm-1连续性，即在相邻单元的交界面上应有函数直至m-1阶的连续导数。,为了从物理意义上加深对收敛准则的理解，我们以平面问题为例加以说明。,收敛准则的物理意义,在平面问题中，泛函中出现的是位移u和v的一次导数，因此m=l。收敛准则I要求插值函数或位移函数至少是x、y的一次完全多项式。,所以完备性准则要求由插值函数所构成的有限元解必须能反映单元的刚体位移和常应变状态。,若不能满足上述要求，那么赋予结点以单元刚体位移(零应变)或常应变的位移值时，在单元内部将产生非零或非常值的应变，这样有限元解将不可能收敛于真正解。,我们知道位移及其一阶导数为常数的项是代表与单元的刚体位移和常应变状态相应的位移模式。,对于平面问题，协调性要求是C0连续性，即要求位移函数x、y的零阶导数，也就是位移函数自身在单元交界面上是连续的。如果在单元交界面上位移不连续，将在交界面上引起无限大的应变，这时必须有发生于交界面上的附加应变能补充到系统的应变能中去，而我们在建立泛函时，没有考虑这种情况，而只考虑了产生于各个单元内部的应变能。因此，倘若边界上位移不连续，有限元解就不可能收敛于真正解。,思考题：请说明平面4节点等参单元的收敛性？,4.6 轴对称问题及空间问题有限元,4.6.1 轴对称问题有限元：,4.6 轴对称问题及空间问题有限元,4.6.1 轴对称问题有限元：,请同学们自己证明,i节点的位移：,单元节点位移向量：,单元内部某点位移：,i节点的型函数：,4.6.1 轴对称问题有限元：,单元内部某点应变：,应变矩阵：,4.6.1 轴对称问题有限元：,节点外力：,1、如果 表示作用在圆周上的单位长度的线载荷的径向分量，,则相应的径向节点力(此时，一个点代表一个圆周)为：,如果 表示作用在圆周上的单位长度的线载荷的轴向分量，,则相应的轴向节点力(此时，一个点代表一个圆周)为：,2、如果 表示作用在圆周上的单位长度的面载荷，(此时，一条线代表一个圆柱面),则相应的节点力为：,3、如果 表示作用在圆周上的单位长度的体载荷，(此时，一个面代表一个圆环体),则相应的节点力为：,请同学们自己证明,4.6.2 空间问题有限元：,4.6 轴对称问题及空间问题有限元,4.6.2 空间问题有限元：,4.6 轴对称问题及空间问题有限元,4.6.2 空间问题有限元：,4.6 轴对称问题及空间问题有限元,