第三章 平面问题的有限元法ppt课件.ppt

《第三章 平面问题的有限元法ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第3章平面问题的有限元法平面问题的有限单元法不仅可以用于计算分析具有平面特征的实际机械结构，还可以通过它掌握有限单元法的基本思想和基本步骤。本章详细介绍平面三角形单元的基本原理、用平面三角形单元进行结构有限元分析的主要步骤，并给出了示例。本章概述3.1平面三角形单元矩阵推导平面弹性问题可以分为两类：平面应力问题：在平面应力问题中，连续体的二维坐标尺寸远大于第三维尺寸（如板），平面的法向应力可以忽略。平面应变问题：在平面应变问题中，连续体的二维坐标尺寸远小于第三维尺寸，加载平面的法向应变可以假设为零,因此分析该连续

2、体的应力和位移时，可以通过分析其法向应变为零的一个横断面来完成。3.1平面三角形单元矩阵推导平面问题可以用最简单的平面三角形常应变单元加以分析。在这节中讨论平面三角形单元的构造方法，给出用平面三角形单元求解平面应力问题的详细过程。平面三角形单元刚度矩阵的推导包括如下6个步骤：1.选择合适的单元，建立坐标系统，进行结构离散2.选择合适的位移函数3.用节点位移表示单元内部各点位移4.用节点位移表达单元内任一点的应变5.用应变和节点位移表达单元内任一点的应力6.单元刚度矩阵的形成3.1平面三角形单元矩阵推导1.选择合适

3、的单元，建立坐标系统，进行结构离散用平面三角形分析时，可以只建立一个整体坐标系OXY。将所有作用在单元上的载荷，都按虚功等效的原则移置到节点上，成为等效节点载荷。由此得到平面问题的有限元计算模型，如图3-1所示。图3-1弹性体和离散化后的有限元计算模型如图3-2所示，节点编号1、2、3按逆时针顺序编排，三个节点的位置坐标分别是,和。图3-2直角坐标系下平面三角形单元的节点位移和节点力对于平面问题，每个节点有x和y两个方向的自由度，即对应的位移是u和v。可以认为三角形单元共有6个自由度，即{}T，相应的单元节点力分

4、量分别为{}T。3.1平面三角形单元矩阵推导1.选择合适的单元，建立坐标系统，进行结构离散三角形单元的6个节点位移分量用列阵表示为三角形单元的节点载荷列阵表示为单元节点载荷列阵和节点位移列阵之间的关系可用下式表示其中，为单元刚度矩阵。对于平面三角形单元，节点位移列阵和节点载荷列阵都是6阶的，单元刚度矩阵是一个6×6阶的矩阵。在有限元中，将利用函数插值、弹性力学几何方程和物理方程、最小势能原理，建立上式具体关系。(3.1)(3.2)(3.3)3.1平面三角形单元矩阵推导2.选择合适的位移函数考虑建立以单元节点位移表

5、示的单元内各点位移的表达式，选择一个简单的单元位移模式，单元内各点的位移可按此位移模式由单元节点位移通过插值得到。设平面三角形单元的位移模式为由于在x和y方向的位移都是线性的，从而保证了沿接触面方向相邻单元间任意节点位移的连续性。(3.4)3.1平面三角形单元矩阵推导3.用节点位移表示单元内部各点位移三角形单元的三个节点必定满足位移模式的要求。将单元三个节点的坐标和三个节点位移都代入位移模式方程可以求解。已知单元三个节点的坐标分别为，，。对于节点1有类似的，节点2、3也按上述方法处理，三个节点的位移模型表达式可以

6、组成方程组利用上式就可求出未知的多项式系数，即，可以求得，(3.5)3.1平面三角形单元矩阵推导得到单元内任意一点（x，y）的位移为式中,N为形函数矩阵(Shapefunctionmatrix)，。平面三角形单元的形函数矩阵具体表达式如下其中,I为2阶单位矩阵,；为三角形单元的面积;(3.6)(3.7)(3.8)3.1平面三角形单元矩阵推导式中的其它各个系数为（1、2、3轮换）(3.9)式（3.4）经过整理可以写成如下展开形式(3.10)3.1平面三角形单元矩阵推导三角形单元用于解决弹性力学平面问题，单元内任一

7、点的应变列阵满足几何方程式中，和是线应变，是剪应变。上式中的分别用位移模式方程式(3.6)或(3.10)代入可求解应变分量。由于N是x,y的函数，对其进行偏微分处理后，可得4.用节点位移表达单元内任一点的应变(3.11)(3.12)3.1平面三角形单元矩阵推导其中,B为单元应变矩阵，其表达式为由于和b1、b2、b3、c1、c2、c3等都是常量，所以平面三角形单元的应变矩阵B中的诸元素都是常量，因而平面三角形单元中各点的应变分量也都是常量。(3.13)上式可简记为(3.14)3.1平面三角形单元矩阵推导对于平面应

8、力问题，一点的应力状态可以用、、这三个应力分量来表示，应力应变关系为式中，D为弹性矩阵，其表达式为其中，E为杨氏模量，为泊松比。把步骤(4)中导出的应变表达式(3.12)代入式(3.15)，可得用节点位移表示的应力。令S为应力矩阵，它是5.用应变和节点位移表达单元内任一点的应力(3.16)(3.17)(3.18)(3.15)3.1平面三角形单元矩阵推导用虚位移原理对图3-