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

思考题： 请说明平面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 轴对称问题及空间问题有限元 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，可以选择二条直线，这二条直线通过除结