有限元(第六章)tg.pdfVIP

第六章 等参数单元 第1节 等参数单元的引入 第2节 四节点四边形等参数单元 第3节 等参数单元平面问题有限元格式 第4节 八节点曲边四边形等参数单元 第5节 数值积分及其应用 小 结 1 第1节 等参数单元的引入 线性单元的位移在单元内呈线性变化，应力、应变在单元内 是一个常量，因此在求解区域内应力和应变的变化都是不连续 的。所以采用线性位移插值函数的三角形单元的计算精度不 高，在许多情况下得不到正确结果，特别是在应力集中部位。 提高计算精度可以采用的方法有： (1) 单元分细； (2) 构造高精度新单元。 将单元分细无疑可以提高计算精度，因为有限元法的计算基 础是当单元无限分细时计算结果将收敛于精确解。但是单元分 细会增加单元数目和节点数目，从而大大增加了所要求解的方 程数目，占用和耗费大量计算机资源。所以，用细分单元的方 法来提高精度有时是不经济的。 2 如图（1）所示任意四节点四边形单元，采用位移插值函数 δ＝α＋αx ＋αy ＋αxy 1 2 3 4 单元不平行于坐标轴的任一条边的直 线方程可写成y ＝k x＋b （k不为0 ）， 将其代入位移插值函数公式，则该边 上的位移变化为δ＝Ax2 ＋B x ＋C，即 位移不再是线性变化的了。该边上的 位移值将不能由其上两个节点处的函 数值所惟一确定。从而在相邻两单元 的公共边将不能保证位移插值函数的 图1 任意四节点四边形单元 连续，也就是变形协调性得不到满 足，常用的解决办法是采用等参数单 元。 3 第2节 四节点四边形等参数单元 • 矩形单元是满足解的收敛性条件的(P86) 。如果 通过一个坐标变换将任意四边形单元变换成四 边形单元，只要 坐标变换中任意四边形单元与矩形单元之间的点是 一一对应的（称为是坐标变换的相容性）； 而变换后的位移插值函数又是满足解的收敛性条件 的； 这两条合在一起就能保证任意四边形单元在原 坐标系中满足解的收敛性条件。 4 首先看看局部坐标系下的位移插值函数、形状函数和收敛性条 件，再讨论坐标变换。 图（2 ）所示正方形单元的位移插值函数采 用双线性函数 + + + u α αξ αη αξη  1 2 3 4  （1） + + + v α αξ αη αξη  5 6 7 8 矩形单元采用双线性位移插值函数时的形 图2 局部坐标下正方形单元 状函数，将它用于正方形单元，a＝b ＝ 1，则形状函数为  1 ( ) ( )( ) N 1 ξ,η 4 1−ξ 1−η