有限元法基础讲稿-第10讲新doc.ppt

《有限元法基础讲稿-第10讲新doc.ppt》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、结构静力学问题的有限元法梁结构问题比桁架复杂一些，也可用矩阵分析法（线性代数方程组）得到问题的精确解。在ANSYS软件中上述两类问题的建模和求解较为简单。定义单元用ET指令，考虑材料各种特性选不同的Link、Beam单元。有限单元法把杆系结构的矩阵分析方法推广应用于连续介质：把连续介质离散化，用有限个单元的组合体代替原来的连续介质，这样一组单元只在有限个节点上相互连接，因而包含有限个自由度，可用矩阵方法进行分析。主题:A.平面问题有限元法B.轴对称问题有限元法C.空间问题有限元法D.等参数有限元

2、法E.单元与整体分析结构静力学问题的有限元法平面问题有限元法对一些特殊情况可把空间问题近似地简化为平面问题，只须考虑平行于某个平面的位移分量、应变分量与应力分量，且这些量只是两个坐标的函数。平面问题分平面应力问题和平面应变问题两类。设有很薄的均匀薄板，只在板边上受有平行于板面并且不沿厚度变化的面力，同时，体力也平行于板面并且不沿厚度变化，记薄板的厚度为t，以薄板的中面为xy面，以垂直于中面的任一直线为z轴．由于板面上不受力，且板很薄，外力不沿厚度变化，可以认为恒有，，不为零的应力分量为、、，这种

3、问题就称为平面应力问题。设有无限长的柱形体，在柱面上受有平行于横截面而且不沿长度变化的面力，同时，体力也平行于横截面且不沿长度变化。以任一横截面为xy面，任一纵线为z轴，由于对称性（任一横截面都可以看做对称面），此时，结构静力学问题的有限元法…平面问题有限元法不为零的应变分量为、、，这种问题就称为平面应变问题。二维连续介质，用有限单元法分析的步骤如下：(1)用虚拟的直线把原介质分割成有限个三角形单元，这些直线是单元的边界，几条直线的交点即为节点。(2)假定各单元在节点上互相铰接，节点位移是基本的

4、未知量。(3)选择一个函数，用单元的三个节点的位移惟一地表示单元内部任一点的位移，此函数称为位移函数（位移模式）。(4)通过位移函数，用节点位移惟一地表示单元内任一点的应变；再利用广义虎克定律，用节点位移可惟一地表示单元内任一点的应力。(5)利用能量原理找到与单元内部应力状态等效的节点力；再利用单元应力与节点位移的关系，建立等效节点力与节点位移的关系。这是有限单元法求解应力问题的最重要的一步。结构静力学问题的有限元法…平面问题有限元法(6)将每一单元所承受的荷载，按静力等效原则移置到节点上。(7

5、)在每一节点建立用节点位移表示的静力平衡方程，得到一个线性方程组；解出这个方程组，求出节点位移；然后可求得每个单元的应力。单元的位移模式及插值函数由于三角形单元对复杂边界有较强的适应能力，因此很容易将一个二维域离散成有限个三角形单元。在边界上以若干段直线近似原来的曲线边界，随着单元增多，这种拟合将越精确。下面以3节点三角形单元为代表讨论平面问题的有限元格式。图2-13节点三角形单元结构静力学问题的有限元法…平面问题有限元法设三角形单元节点编码为i,j,m，以逆时针方向编码为正向，否则后面求出的面

6、积A为负值。每个节点有2个位移分量如图2-1所示，节点位移为δe=[ui,vi,uj,vj,um,vm]T节点的坐标分别为(xi,yi)、(xj,yj)、(xm,ym)。在有限单元法中单元的位移模式（也称位移函数和插值函数）一般采用多项式作为近似函数，因为多项式运算简便，并且随着项数的增多，可以逼近任何一段光滑的函数曲线。多项式的选取应由低次到高次。3节点三角形单元位移模式选取一次多项式单元内的位移是坐标x,y的线性函数。β1～β6是待定系数，称之为广义坐标。6个广义坐标可由单元的6个节点位移来

7、表示。在式（2-1-1）中代入三角形单元各节点的坐标然后解出(2-1-1)结构静力学问题的有限元法…平面问题有限元法式中，A为三角形的面积上式中(i，j，m)表示下标轮换，如i→j，j→m，m→i。(2-1-2)(2-1-3)结构静力学问题的有限元法…平面问题有限元法将求得的广义坐标代入式(2-1-1)，可将单元位移函数表示成节点位移的函数其中Ni，Nj，Nm称为单元的插值基函数或形函数，它是坐标x、y的一次函数。单元上任一点的形函数之和为1。式（2.1.4）写为矩阵的形式其中是二阶单位阵。(2

8、-1-4)(2-1-5)(2-1-6)结构静力学问题的有限元法…平面问题有限元法应变矩阵确定了单元位移后，可以很方便地利用几何方程和物理方程求得单元的应变和应力。作为平面问题，单元内具有3个应变分量εx、εy、γxy（各符号的意义见附录1），用矩阵表示为将（2-1-4）式代入上式中，得到结构静力学问题的有限元法…平面问题有限元法或式中B称为应变矩阵，写为分块形式，即B=[BiBjBm]而其子阵为3节点三角形单元的B是常量阵，所以称为常应变单元。在应变梯度较大(也即应力梯度较大)的部位，单元划分应