有限单元法.ppt

有限单元法(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　等效结点荷载 　　（2）框图设计 　　计算方案选定以后，需先画出总框图。平面问题的有限元法的总框图如图3.11所示。 　　（3）设计程序 　　根据总框图的规划和有限元法的基本原理，用自己所使用的程序语言把这个解决方案严格地描述出来，也就是编写出程序代码。它就是用于求解问题的源程序。 　　（4）程序的编译和调试 　　将编写好的源程序输入计算机，进行编译，改正其中的语法等错误。并用实际的输入数据对编好的程序进行调试，分析所得到的运行结果，进行程序的测试和调整，直至获得预期的结果。 　　教学提示：空间轴对称结构是工程中一类特殊的结构。对于这类结构，由于其几何形状、约束情况、所受荷载都对称于某一根轴（如图4.1所示），因此，通常是在柱坐标系中进行分析。本章以三角形单元为例介绍了这一类结构的有限元分析过程。 　　教学要求：本章要求学生了解柱坐标系中三角形单元的建立方法，熟悉轴对称结构有限元分析的基本方法。 　　在工程中有许多结构，如活塞、厚壁容器等，他们的几何形状、约束情况及所受的荷载都对称于空间的某一根轴，因此在物体中通过该轴的任何平面都是对称面，所有应力、应变和位移也对称于该轴，这类问题称为空间轴对称问题。研究轴对称问题时通常采用圆柱坐标系　　　　，以　轴为对称轴。 　　如图4.1所示的受均布内压作用的长圆筒，通过z轴的一个纵截面就是对称面。由于对称性，轴对问题共有4个应力分量： 其中　表示沿半径方向的正应力，称为径向应力；　表示沿　方向的正应力，称为环向应力或切向应力；　表示沿ｚ方向的正应力，称为轴向应力；　表示在圆柱面上沿ｚ方向作用的剪应力。 　　同样，轴对称问题共有4个应变分量： 其中　表示沿半径方向的正应变，称为径向正应变；　表示沿　方向的正应变，称为环向正应变或切向正应变；　表示沿ｚ方向的正应变，称为轴向正应变；　表示沿ｒ和ｚ方向的剪应变。 　　在轴对称问题中，弹性体内任意一点上，不存在切向位移，只存在径向位移 u 和轴向位移 w ，两个位移分量表示为， 　　轴对称问题分析中所使用的三结点单元，在对称面上是三角形，在整个弹性体中是三棱圆环，各单元中圆环形铰相联接。 参照平面问题的三角形单元位移函数，轴对称问题的三结点三角形单元位移函数取为， 　　按照平面问题三角形单元的分析过程，将结点坐标和结点位移代入（4-4）得到， 其中： 定义形函数为： 用矩阵表示的单元位移为： 轴对称问题的几何方程： 由（4-9）式得， 其中， 用几何矩阵表示单元的应变， 　　由于　在是坐标 r、z 的函数，　分量在单元中不为常量，其它三个应变分量在单元中仍为常量。 由轴对称问题的物理方程，得到弹性矩阵 D， 　　由弹性矩阵 D 和几何矩阵 B 可以得到应力矩阵 S，并计算出单元内的应力分量， 其中： 式中： 有了单元应力场和应变场，可以利用虚位移原理或最小势能原理建立单元刚度矩阵 单元刚度矩阵的分块矩阵为， 　　由于几何矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算，可以用三角形单元形心位置的坐标　　代替 B 矩阵中的变量　　。 应变矩阵（4-13）变成， 其中： 单元刚度矩阵的近似表达式为： 单元刚度矩阵的分块矩阵近似表达式为， 这里： 将竖向位移和结点坐标代入（3-4）中的第二式，可以得到： 　　　 　　 （3-8） 将（3-7）、（3-8）代回（3-4）整理后可得： 　　　　　　　　　　　　　　　　　　　　　　　 （3-9） 　　　　 　　　　　　　　　　　　　　　　　　　（3-10） 其中系数， 　　　　　　　　　　　　　　　（下　　　　标轮换） 令　　　　 （3-11） 可得： 　　　 　 　　　　　（3-12） 单元内的位移场函数可以简写成： 　　　　 　 （3-13） 把　　称为形函数矩阵，　称为形函数。根据形函数的定义，　　具有以下性质： 　　（1）在单元结点上形函数的值为1或为0。