用于振动分析的有限元方法

《用于振动分析的有限元方法》由会员分享，可在线阅读，更多相关《用于振动分析的有限元方法（59页珍藏版）》请在文档大全上搜索。

1、用于振动分析的有限元方法指导老师：陈益 报告人：成志斌 韩宗彪 何瑜 宁鹏内容有限元介绍单个元素的运动方程整个系统的运动方程整个系统的边界条件的加载及质量矩阵 MATLAB实例及总结单个元素的质量矩阵、刚度矩阵、力矢量及其转化有限元法简介 有限元法是一种可用于精确地(但近似)解决许多复杂的振动问题的数值方法。 对于基本的一维元素进行有限元分析，能得到质量矩阵与刚度矩阵和所需的力矢量，对于二维三维，元素矩阵会转换成相关的更高维的空间。使用一致的和集中质量矩阵的有限元方程并结合边界条件能为复杂系统提供解释。 最后,使用MATLAB程序得到在轴向载荷下的指定节点位移,固有振动频率和特征值分析。本章目

2、的 *认识用于解决不同类型振动问题的刚度和质量矩阵。 *将矩阵元素从局部坐标系变换到全球坐标系。 *装配单元矩阵和应用边界条件。 *对杆、梁元素进行静态分析。 *对杆、梁元素进行动态分析来得到固有频率和振型。 *在有限元振动分析使用一致的集中质量矩阵。 *使用MATLAB解决振动问题。有限元思想 1，实际结构被一些元素所取代,这些元素都是被假定为一个连续的结构部件即有限元，这些元素在特定点即节点上互相关联。 2，如果解决方案的各方面都选择得当,那么它可以收敛到精确的解决方案，因为组成总体结构的元素很小，在节点上的力的平衡和元素之间的位移都令人感到满意,这样整个结构(组合的元素)表现为单一实体。

3、 3，因为得到准确解很难，所以得到一个方便且逼近的近似解很有价值。元素的运动方程龙门刨铣床有限元模型三角板元素梁元素元素的运动方程 位移函数 形状函数 各点对应位移 未知节点位移数 n 动能 应变能质量矩阵刚度矩阵 位移函数 形状函数 各点对应位移 未知节点位移数 n 动能 应变能主要内容： 一，单元的质量、刚度矩阵、等效节点力 二，单元矩阵的坐标变换 三，整个系统的运动方程一 单元的质量、刚度矩阵，等效节点力矢量一 杆单元一个杆单元是从杆上划分出的一个小段，如下图所示。由于单元很小，、A均视为常量。现在就以这最简单的杆单元，推导出它的质量、刚度矩阵，等效节点力。图12.1图12.1 (1)

4、求杆单元上任意点的位移u(x,t)本来，杆单元上任意点的位移u(x,t)与节点的位移u1(t)、u2(t)之间的关系是未知的，但是，只要单元划分的足够小，那么其间的关系就无关大局。所以可以假定它们之间有简单的线性关系，即根据节点位移对单元内任意点位移进行插值：)()(),(2211tututxu（1）式中，1、 2称为线性系数，与单元里点的位置有关，是x的函数。此函数与单元的形状有关，又叫形状函数。形状函数和插值函数一样是任意的，但必须边界条件：)(), 0(1tutu)(),(2tutlu只有满足此条件单元才能协调一致运动，而不致破坏系统的完整性，因此这两个条件实际上就是变形协调条件。将式（

5、1）带入（2）中，就可以得到形状函数1(x)、2(x)所满足的边界条件： （2） 1 , 00 0 , 102121ll （3） 以上边界条件确定了 xx21、 lxxlxx21 ,1由于这两个函数的任意性我们可以用简单的线性函数来近似，因此有：（4）代回（1）式中有： tulxtulxtxu21)1 (),(（5）为此，我们已经找到了用节点位移表示单元内任意一点位移的表达式。 （2）计算此单元的动能和势能杆单元的动能可表示成：（6） 上式中，是材料的密度，A是杆单元的横截面积。用矩阵形式表示（6）式为：（7）其中，（8）所以，质量矩阵可以认为是：（9） 杆单元的势能可以写成：（10）式中，E

6、是弹性模量，（10）表示成矩阵形式为：这里，所以刚度矩阵k可以表示成：（11）（12）(3)计算等效节点力设单元上x处作用有分布力f ( x , t)，现在要把它等效成节点力 tftf21,遵循等效原则，即原载荷和等效之后的节点载荷在虚位移上所做的虚功相等。其实，就是对应于广义坐标 tutu21,的广义力，为此，计算txf,所做的虚功：把上式写成矩阵形式：（13）所以等效节点力可以写成： tftuktum （14）二 梁单元 ,31tftf如下图所示，一个梁单元也是有两个节点，但是有四个自由度，每个节点处，有两种位移形式，一个是线位移，即挠度，一种是角位移。图中，tftf42, 是力，是力矩。

7、txf,是分布载荷 ,31tt tt42, 是对应的线位移，是对应的转角。tx,是梁单元上任意位移 x处的挠度。图12.2在静载弯曲条件下，梁单元上任意点出的挠度是x的三次方程，可写成：此方程必须满足下面的边界条件：由此可以求解处a (t)、b (t)、c (t)、d (t)，进而挠度方程为：（16）（17）（18）tan挠度的斜率上式可以写成形状函数的表示：其中，形函数分别为：梁单元的动能、势能、虚功表达式分别为：（19）式中I是横截面的惯性矩上式中：（20）（21）（22）通过上式，可以得到梁单元的质量、刚度矩阵，等效节点力：二 单元矩阵的坐标变换局部坐标系：以各个单元本身的轴线为基准所设