有限单元法-上课.pptVIP

有限单元法(Finite Element Method) 2.2　局部坐标系中的杆单元分析 2.2.1　拉压杆单元 2.3　杆系结构的整体分析 2.3.1　平面问题坐标转换矩阵 2.4　等效结点荷载与边界条件的处理 2.4.1　非结点荷载的处理 2.5　杆系结构分析算例 3.1　概述 3.1.1　单元划分类型 第4章空间轴对称问题的有限元分析 4.1　概述 4.2　三角形单元 4.2.1　位移函数的选取 4.2.2　单元刚度矩阵 4.2.3　等效结点荷载 　　如图4.1所示的受均布内压作用的长圆筒，通过z轴的一个纵截面就是对称面。由于对称性，轴对问题共有4个应力分量： 其中　表示沿半径方向的正应力，称为径向应力；　表示沿　方向的正应力，称为环向应力或切向应力；　表示沿ｚ方向的正应力，称为轴向应力；　表示在圆柱面上沿ｚ方向作用的剪应力。 　　同样，轴对称问题共有4个应变分量： 其中　表示沿半径方向的正应变，称为径向正应变；　表示沿　方向的正应变，称为环向正应变或切向正应变；　表示沿ｚ方向的正应变，称为轴向正应变；　表示沿ｒ和ｚ方向的剪应变。 　　在轴对称问题中，弹性体内任意一点上，不存在切向位移，只存在径向位移 u 和轴向位移 w ，两个位移分量表示为， 　　轴对称问题分析中所使用的三结点单元，在对称面上是三角形，在整个弹性体中是三棱圆环，各单元中圆环形铰相联接。 参照平面问题的三角形单元位移函数，轴对称问题的三结点三角形单元位移函数取为， 　　按照平面问题三角形单元的分析过程，将结点坐标和结点位移代入（4-4）得到， 其中： 定义形函数为： 用矩阵表示的单元位移为： 轴对称问题的几何方程： 由（4-9）式得， 其中， 用几何矩阵表示单元的应变， 　　由于　在是坐标 r、z 的函数，　分量在单元中不为常量，其它三个应变分量在单元中仍为常量。 由轴对称问题的物理方程，得到弹性矩阵 D， 　　由弹性矩阵 D 和几何矩阵 B 可以得到应力矩阵 S，并计算出单元内的应力分量， 其中： 式中： 有了单元应力场和应变场，可以利用虚位移原理或最小势能原理建立单元刚度矩阵 单元刚度矩阵的分块矩阵为， 　　由于几何矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算，可以用三角形单元形心位置的坐标　　代替 B 矩阵中的变量　　。 应变矩阵（4-13）变成， 其中： 单元刚度矩阵的近似表达式为： 单元刚度矩阵的分块矩阵近似表达式为， 这里： 　　（1）集中力 　　集中力的处理很简单，一般直接把集中力作用点取为结点，不需要作特殊处理，就可以直接把集中力加入到结点荷载列阵中去。 　　（2）体积力 　　设单元内单位体积上作用的体积力为　，则移置到单元各结点的等效结点力为 3.2.4　单元刚度矩阵 　　（1）单元的应变能 　　 （3-23） 　　（2）单元上外力的势能 　　　　 （3-24） 式中　　、　　、　　分别表示单位体积的体积力、单元上的表面力、单元结点上的结点荷载。 将式（3-13）、（3-15）、（3-18）代入（3-23）、（3-24）后相加，得到单元的总势能为： 利用最小势能原理，取结点位移　　的变分，得到： 由　　的任意性，有： 考虑到　　　的对称性，对式（3-25）求偏导得到： 记： 则式（3-28）可写为： 这就是描述单元结点力和结点位移向量之间关系的平衡方程。其中　称为单元刚度矩阵。在3结点等厚三角形单元中　和 的分量均为常量，则单元刚度矩阵可以表示为： 其中t、Ａ分别为单元的厚度和面积。单元刚度矩阵　可以表示为分块矩阵的形式： 其中： 对于平面应力问题，其刚度矩阵的显式： 3.2.5　等效结点荷载 　　在式（3-28）中括号中的前两个量分别为体积力、表面力移置到结点上的等效结点力，依次定义为　　、　　，即： 、 那么（3-30）中的荷载列阵就等于体积力、表面力和单元结点荷载叠加。因此作用在弹性体上的外力，需要移置到相应的结点上成为结点荷载。荷载移置要满足静力等效原则。 3.2.6　整体刚度矩阵 　　得到了单元刚度矩阵后，需要将一系列的将单元组成一个整体结构，然后根据结点载荷平衡的原则进行分析，得到整体刚度矩阵。整体分析包括以下4个步骤： 　　（1）建立整体刚度矩阵， 　　（2）根据支承条件修改整体刚度矩阵， 　　（3）解方程组，求出结点的位移， 　　（4）根据结点位移，求出单元的应变和应力。 　　由单元刚度矩阵得到整体刚度矩阵的基本方法是刚度集成法，即整体刚度矩阵是单元刚度矩阵的集成。 单元编号 单元结点局部编号 单元结点整体编号 （1） 3 1 2 （2） 5 2 4 （3）