有限元-4-等参

第四章等参单元本章包括以下内容：4.1等参单元的基本概念4.2四边形八节点等参单元4.3等参单元的单元分析4.4六面体等参单元4.1等参单元的基本概念在进行有限元分析时，单元离散化会带来计算误差，主要采用两种方法来降低单元离散化产生的误差：1）提高单元划分的密度，被称为h方法（h-method）；2）提高单元位移函数多项式的阶次，被称为p方法（p-method）。在平面问题的有限单元中，我们可以选择四结点的矩形单元，如图5-1所示，该矩形单元在x及y方向的边长分别为2a和2b。图4-1四结点矩形单元同第三章的方法类似，将单元的位移模式选为，（4-1）可得到，（4-2）形态函数为，xyayaxaau4321xyayaxaav8765ppmmjjiiuNuNuNuNuppmmjjiivNvNvNvNv)1)(1(41byaxNi)1)(1(41byaxNj)1)(1(41byaxNm)1)(1(41byaxNp上述单元位移模式满足位移模式选择的基本要求：1）反映了单元的刚体位移和常应变，2）单元在公共边界上位移连续。在矩形单元的边界上，坐标x和y的其中一个取常量，因此在边界上位移是线性分布的，由两个结点上的位移确定。与三结点三角形单元相比，四结点矩形单元的位移模式是坐标的二次函数，能够提高计算精度，但也有显著的缺点，两种单元的比较如下。表5-1三结点三角形单元与四结点矩形单元比较单元类型优点缺点三结点三角形单元适应复杂形状，单元大小过渡方便计算精度低四结点矩形单元单元内的应力、应变是线性变化的，计算精度较高不能适应曲线边界和非正交的直线边界如果任意形状的四边形四结点单元采用矩形单元的位移模式，则在公共边界上不满足位移连续性条件。为了既能得到较高的计算精度，又能适应复杂的边界形状，可以采用坐标变换。图4-2四结点任意四边形单元图4-3四结点正方形单元在图4-2所示的任意四边形单元上，用等分四条边的两族直线分割四边形，以两族直线的中心为原点，建立局部坐标系，沿及增大的方向作为轴和轴，并令四条边上的及值分别为。为了求出位移模式，以及局部坐标与整体坐标之间的变换式，在局部坐标系中定义一个四结点正方形单元，如图4-3所示。),(1参照矩形单元，四结点正方形单元的位移模式为，（4-4）其中，（4-5）44332211uNuNuNuNu44332211vNvNvNvNv)1)(1(411N)1)(1(412N)1)(1(413N)1)(1(414N四个结点的坐标为，定义新的变量，，(i=1，2，3，4)(4-6)形态函数表示为，（i=1，2，3，4）（4-7）把及作为任意四边形单元的局部坐标，把（4-4）的位移模式和（4-7）的形态函数用于任意形状的四边单元，可得：1）在四个结点处可以得到结点的位移；2）在单元的四条边上，位移线性变化，保证了单元公共边界上位移的连续性。因此给出任意四边形单元的结点位移就能得到整个单元上的位移，（4-4）的位移模式就是所要找的正确的位移模式。),(iii0i0)1)(1(4100iN把局部坐标与整体坐标的变换式也取为，将坐标变换式用于任意四边形单元，可得：1）在四个结点处给出结点的整体坐标，2）在四条边上的整体坐标是线性变化的。只要给出任意四边形单元四个结点的整体坐标，用（4-8）式就可以建立局部坐标系中的正方形单元和整体坐标系中的任意四边形单元之间的坐标变换关系。把图4-3中的局部坐标系中的正方形单元称为基本单元。把图4-2中的在整体坐标系中的任意四边形单元看作由基本单元通过坐标变换得来的，称为实际单元。单元几何形状和单元内的未知量采用相同数目的结点参数以及相同的插值函数进行变换，称为等参变换。采用等参变换的单元，称为等参单元。由于形态函数，正好反映了单元形状的变化，也称为形函数（Shapefunction）。（4-8）44332211xNxNxNxNx44332211yNyNyNyNy采用等参单元，使我们可以在局部坐标系中的规则单元上进行单元分析，然后再映射到实际单元上。等参单元同时具有计算精度高和适用性好的特点，是有限元程序中主要采用的单元形式。4.2四边形八节点等参单元为了更好地反映物体内的应力变化，适应曲线边界，在弹性力学平面问题的分析中经常使用四边形八节点等参单元。如图5-4所示，由于每条边上增加了一个结点，单元的边是一条二次曲线，可以更好地适应曲线边界，图4-4四边形八结点单元图4-5八结点基本单元对于等参单元，先在图5-5所示的八结点基本单元上进行分析。八结点单元一共有16个已知的结点位移分量，基本单元中取如下的位移模式：该位移模式实际上是一个双二次函数，待定系数由结点位移分量确定。在单元的每条边上，位移是局部坐标的二次函数，完全由边上的三个结点的位移值确定，所以这个位移模式满足位移连续性条件。282726524321aaaaaaaau282726524321bbbbbbbbv