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

04030102对于平面问题，协调性要求是C0连续性，即要求位移函数x、y的零阶导数，也就是位移函数自身在单元交界面上是连续的。如果在单元交界面上位移不连续，将在交界面上引起无限大的应变，这时必须有发生于交界面上的附加应变能补充到系统的应变能中去，而我们在建立泛函时，没有考虑这种情况，而只考虑了产生于各个单元内部的应变能。因此，倘若边界上位移不连续，有限元解就不可能收敛于真正解。4-1矩形单元与插值函数4-2高次单元与插值函数(6结点三角形单元)4-3等参单元和单元坐标变换4-4数值积分方案4-5有限元解的收敛性4-6轴对称问题有限元及空间问题有限元平面问题的有限元分析中，矩形单元也是常用的单元之一。它采用了比常应变三角形单元更高次的位移模式，故可更好地反映弹性体中的位移状态与应力状态。02014-1矩形单元与插值函数单元节点位移向量和单元位移向量分别为：单元位移函数：其中为型函数。矩阵形式为：将4个节点的坐标值代入上式，可解出上式的8个系数。于是，单元位移可通过单元节点位移表示。单元坐标变换06局部坐标04其中01全局坐标03因此有：05它们的几何意义为02此单元位移模式是具有完全一次式的非完全二次式。因此单元内位移沿xy方向都是线性变化的，但沿其他方向则不是线性的。单元每边界上有两个节点，按线性变化的位移形式，单元变形后此两节点确定了惟一的单元直线边界，因而此单元边界间的位移是连续的，所以矩形单元的这种位移模式又称为双线性位移模式。正是由于在单元的边界上位移按线性变化，且相邻单元公共结点上有共同的结点位移，因此保证了两个相邻单元在其公共边界上位移的连续性。这种单元的位移模式是完备的和协调的，满足解的收敛条件，因此4结点矩形元是协调元。矩形单元中的应力分量都不是常量。01正应力分量σx的主项(即不与μ相乘的项)沿y方向按线性变化，而它的次要项(与μ相乘的项)沿x方向按线性变化。正应力σy与此相反。02剪应力分量则沿xy两个方向都呈线性变化。03这种一个方向为常量，另一方向呈线性变化的情况通常并不能提高单元的精度。04矩形单元给出某种线性应力场，解题时可能给出比常应力三角形单元更好的结果，其数据准备也更方便。但是，三角形单元的网格划分比较灵活，适应性更强，疏密网格的变化较为方便。解决上述问题的方法之一是同时采用矩形单元和三角形单元混合使用，如图所示的计算模型，在直边区改用矩形单元,而在曲边附近采用三角形单元，以达到两者兼顾的效果。矩形单元明显的缺点是不能很好地符合曲线边界，包括与坐标轴不平行的直线边界，因此直接应用受到限制。更为一般的方法是通过等参变换将局部坐标内的矩形单元变换为总体坐标内的任意四边形(包括曲边四边形)的单元4-26结点三角形单元单元的6个结点是三角形的3个角点及3个边的中点。六个节点，12个节点位移能确定12个多项式的系数。位移函数取完全的二次多项式。单元边界上位移按二次抛物线分布。这种单元的应变在两个坐标方向上都呈线性变化。应力也呈线性变化。因此，单元精度较3结点三角形单元要高。六个节点，12个节点位移能确定12个多项式的系数。它采用了比常应变三角形单元更高次的位移模式，故可更好地反映弹性体中的位移状态与应力状态：位移函数中的常数项和完整的一次式满足收敛条件中对完备性的要求。123454-26结点三角形单元如角点1，可利用2-3边的边方程以及过边中点4-6的直线方程，并使在结点1上N1=1,即用面积坐标给出的6节点三角形单元的插值函数：对于每个结点i，可以选择二条直线，这二条直线通过除结点i以外的所有结点，利用直线方程的左部作为线性函数来构造插值函数。N1=L1(2L1－1)N2=L2(2L2－1)N3=L3(2L3－1)同理：三个角点的插值函数可统一写作Ni=Li(2Li－1)三个边中点，可选择非此边中点所在的其他二条边作为直线方程，则可得N4=4L1L2N5=4L2L3N6=4L1L34-3等参单元对于一个给定问题的求解域，预期用较少的单元即可获得需要精度的解答。但是用较少的形状规则的单元离散几何形状复杂的求解域常会遇到困难，因此需要寻找适当的方法将规则形状的单元转化为能够适应边界的单元。4.3.1等参变换的概念在有限单元法中最普遍采用的变换方法是等参变换，即单元几何形状的变换和单元内的场(位移)函数的变换采用相同数目的结点参数及相同的插值函数进行。01借助于等参单元可以对于一般的任意几何形状的工程问题方便地进行有限元离