【2017年整理】弹性力学平面问题的有限元方法.ppt

第六章 用有限元法解平面问题;第六章 用有限元法解平面问题;第六章 用有限单元法解平面问题;简史; 1970年后，FEM被引入我国，并很快地得到应用 和发展。;导出方法;§6-1 基本量和基本方程的 矩阵表示 ;基本物理量：; 物理方程:; --结点虚位移; --对应的虚应变。; 3.整体分析。 ; 结构力学研究的对象是离散化结构。如桁架， 各单元（杆件）之间除结点铰结外，没有其他联 系（图（a））。; 将连续体变换为离散化结构（图（c））： 即将连续体划分为有限多个、有限大小的单元， 并使这些单元仅在一些结点处用绞连结起来，构 成所谓‘离散化结构’。; 图（c）与图( a）相比，两者都是离散 化结构；区别是，桁架的单元是杆件，而 图（c）的单元是三角形块体（注意：三角 形单元内部仍是连续体）。;2.单元分析;（1）应用插值公式, 由单元结???位移 ，求单元的位移函数;（4）应用虚功方程，由单元的应力 ， 求出单元的结点力，表示为;－－单元对结点的 作用力，与 数 值相同,方向相反， 作用于结点。;（5）将每一单元中的各种外荷载，按虚功 等效原则移置到结点上，化为结点荷 载，表示为 ; 为已知值, 是用结点位移表示的值。 通过求解联立方程，得出各结点位移值，从而求出各单元的应变和应力。;求解方法;思考题;应用插值公式，可由 求出位移 。; 插值公式（a）在结点 应等于结点位移值 。由此可求出 ;将式（a）按未知数 归纳为:;N －－ 称为形（态）函数矩阵。; A为三角形 的面积（图示坐标系中， 按逆时针编号），有：; 三结点三角形单元的位移模式，略去了2次以 上的项，因而其误差量级是 且其中只包含 了 的1次项，所以在单元中 的分布如图 （a）所示， 的分布如图（b）、（c）所示。 ; 所以当单元趋于很小时，即 时，为了使FEM之解逼近于真解。则为了保证FEM收敛性,位移模式应满足下列条件： ; 因为当单元 时，单元中的位移和应变都趋近于基本量－－刚体位移和常量位移。;收敛性条件;（3）位移模式应尽可能反映位移的连续性。 即应尽可能反映原连续体的位移连续 性。在三角形单元内部，位移为连续；在两单元边界ij 上， 之间均为线性变化，也为连续。; （1）和（2）是必要条件，而加上（3）就为充分条件。;思考题;§6-4 单元的应变列阵和应力列阵 ;应用几何方程，求出单元的应变列阵：;应变; 对于线性位移模式，求导后得到的应变和应力，均成为常量，因此，称为常应变（应力）单元。应变和应力的误差量级是 其精度比位移低一阶，且相邻单元的应力是跳跃式的。;思考题;§6-5 单元的结点力列阵与劲度矩阵 ;（2）单元与周围的单元在边界上已没有联 系，只在结点 互相联系。;假想将单元与结点i 切开，则： ; 按虚功方程，在虚位移上，外力的虚功等于应力的虚功。; 假设发生一组结点虚位移　　 则单元内 任一点（x,y）的虚位移为 单元内 任一点（x,y）的虚应变为 代入虚 功方程：在单元中，外力（结点力 ）在虚 位移（结点虚位移　 　）上的虚功，等于应 力 在虚应变 上的虚功，即： ;其中 与 无关，故式(a) 成为;式（c）是由结点位移求结点力的一般公式， －－称为单元的劲度矩阵;对于三角形单元，B 矩阵内均为常数， 有;（1） 是6×6的方阵， 中每一个元素都表示 单元各结点沿坐标方向发生单位位移时所 引起的结点力。;（4）由（3）可导出行列式 。; （书中P.117页），以直角三角形单元为例，计算了应力转换矩阵S和单元劲度矩阵 。 从例题中可以看出，将单元边界上的应力向结点移置，化为作用于结点上的力，正好就是结点力。在FEM中，单元边界之间的联系和相互作用力，都向结点简化，归结成为结点的铰结和结点力。 ;§6－6　荷载向结点移置　 单元的结点荷载列阵 ;（2）变形体静力等效原则－－在任意的虚位移上，使原荷载与移置荷载的虚功相等。 ; 假设发生一