弹性力学第六章_用有限单元法解平面问题教材教学课件.ppt

第六章 用有限元法解平面问题;第六章 用有限元法解平面问题;第六章 用有限单元法解 平面问题;简史; 1956年，特纳等人提出了FEM。 20世纪50年代，平面问题的FEM建立， 并应用于工程问题。 1960年提出了FEM的名称。 20世纪60年代后，FEM应用于各种力学 问题和非线性问题，并得到迅速发展。 1970年后，FEM被引入我国，并很快地得 到应用和发展。;导出方法;§6-1 基本量和基本方程的 矩阵表示 ;面力 位移函数 应变 应力 结点位移列阵 结点力列阵 ; 结力研究的对象是离散化结构。如桁架，各单元（杆件）之间除结点铰结外，没有其他联系（图（a））。; 将连续体变换为离散化结构（图（c））：即将连续体划分为有限多个、有限大小的单元，并使这些单元仅在一些结点处用绞连结起来，构成所谓‘离散化结构’。; 与 相比，两者都是离散化结构；区别是，桁架的单元是杆件，而图（c）的单元是三角形块体（注意：三角形单元内部仍是连续体）。;2.应用结构力学方法(位移法)进行求解: ; （2) 应用插值公式, 由单元结点位 移 ，求单元的位移函数;（5）应用虚功方程，由单元的应力 ，求出 单元的结点力，表示为; ——结点对单元的作用力，作用 于单元，称为结点力，以正标向为正。 ——单元对结点 的作用力，与 数 值相同,方向相反， 作用于结点。;（6）将每一单元中的各种外荷载，按虚功 等效原则移置到结点上，化为结点荷 载，表示为 ;各单位移置到i 结点上的结点荷载 其中 表示对围绕i 结点的单元求和；;结力法求解; 思考题 1. 桁架的单元为杆件，而平面体的单元为三角形块体，在三角形内仍是作为连续体来分析的。试考虑后者在用结构力学方法求解时，将会遇到什么困难？ 2. 在平面问题中，是否也可以考虑其它的单元形状，如四边形单元？ ;FEM是取结点位移 为基本未知数的。但其中每一个单元仍是连续体，所以按弹力公式求应变、应力时，必须首先解决：如何由单元的结点位移 来求出单元的位移函数 应用插值公式，可由 求出位移d。这个插值公式表示了单元中位移的分布形式，因此称为位移模式。; 泰勒级数展开式中，低次幂项是最重要的。∴三角形单元的位移模式，可取为 插值公式 在结点 应等于结点位移值 由此可求出 ; 其中 包含 将式 按未知数 归纳，可表示为 或用矩阵表示为 ;N — 称为形（态）函数矩阵。; A为三角形 的面积（图示坐标系中， 按逆时针编号），; 三结点三角形单元的位移模式，略去了二次以上的项，因而其误差量级是 且其中只包含了 的一次项，所以在单元中 的分布如图（a）所示， 的分布如图 所示。 ; FEM中以后的一系列工作，都是以位移模式为基础的。 所以当单元趋于很小时，即 时，为了使FEM之解逼近于真解，即为了 保证FEM收敛性,位移模式应满足下列 条件： ;（1）位移模式必须能反映单元的刚体位移。 （2）位移模式必须能反映单元的常量应变。 因为当单元 时，单元中的位移和应 变都趋近于基本量——刚体位移和常量 位移。 ;收敛性条件;（3）位移模式应尽可能反映位移的连续性。 即应尽可能反映原连续体的位移连续 性。 在三角形单元内部，位移为连续； 在两单元边界ij 上， 之间均为线 性变化，也为连续。; 为了保证FEM的收敛性