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

用于振动分析的有限元方法 指导老师：陈益 报告人：成志斌 韩宗彪 何瑜 宁鹏 内容 有限元法简介 有限元法是一种可用于精确地(但近似)解决许多复杂的振动问题的数值方法。 对于基本的一维元素进行有限元分析，能得到质量矩阵与刚度矩阵和所需的力矢量，对于二维三维，元素矩阵会转换成相关的更高维的空间。使用一致的和集中质量矩阵的有限元方程并结合边界条件能为复杂系统提供解释。 最后,使用MATLAB程序得到在轴向载荷下的指定节点位移,固有振动频率和特征值分析。 本章目的 *认识用于解决不同类型振动问题的刚度和质量矩阵。 *将矩阵元素从局部坐标系变换到全球坐标系。 *装配单元矩阵和应用边界条件。 *对杆、梁元素进行静态分析。 *对杆、梁元素进行动态分析来得到固有频率和振型。 *在有限元振动分析使用一致的集中质量矩阵。 *使用MATLAB解决振动问题。 有限元思想 1，实际结构被一些元素所取代,这些元素都是被假定为一个连续的结构部件即有限元，这些元素在特定点即节点上互相关联。 2，如果解决方案的各方面都选择得当,那么它可以收敛到精确的解决方案，因为组成总体结构的元素很小，在节点上的力的平衡和元素之间的位移都令人感到满意,这样整个结构(组合的元素)表现为单一实体。 3，因为得到准确解很难，所以得到一个方便且逼近的近似解很有价值。 元素的运动方程 元素的运动方程 位移函数 形状函数 各点对应位移 未知节点位移数 n 动能 应变能 主要内容： 一，单元的质量、刚度矩阵、等效节点力 二，单元矩阵的坐标变换 三，整个系统的运动方程 一 单元的质量、刚度矩阵，等效节点力矢量 (1) 求杆单元上任意点的位移u(x,t) 形状函数和插值函数一样是任意的，但必须边界条件： 以上边界条件确定了 （2）计算此单元的动能和势能 上式中，ρ是材料的密度，A是杆单元的横截面积。 杆单元的势能可以写成： 其实，就是对应于广义坐标 所以等效节点力可以写成： 二 梁单元 在静载弯曲条件下，梁单元上任意点出的挠度是x的三次方程，可写成： 上式可以写成形状函数的表示： 通过上式，可以得到梁单元的质量、刚度矩阵，等效节点力： 二 单元矩阵的坐标变换 局部坐标系：以各个单元本身的轴线为基准所设立的坐标系。便于计算节点位移。 缺点：如果整个系统里各个单元取向各异，各个节点位移方向不一致，如下图。如 何使汇交于一个节点的各个杆件的节点位移真正相等? 解决方法：进行坐标变换 如下图，节点位移在局部、全局坐标系中的关系： 得到在全局坐标下的单元质量、刚度矩阵为： 三 全系统运动方程 经过坐标变换，各个单元的节点位移方向被统一起来，但是不同的节点有不同的节点位移，为了便于综合出全系统的运动方程，首先要建立全系统的节点位移向量。 例如，图的12.5中的1单元，方程，变为： 这样就得到了整个系统的质量矩阵: 12.8 MATLAB应用举例 例12.5 阶梯轴的有限元分析 图12.11中的阶梯轴满足一下条件：A1=16×10-4 m2 ， A2=9×10-4 m2 ， A3=4×10-4 m2 ，Ei=20×1010Pa，i=1,2,3，pi=7.8×103Kg/m3， i=1,2,3，l1=1m，l2=0.5m，l3=0.25m。编写一个MATLAB程序解决以下问题。 a，在在载荷p3=1000N下u1，u2，u3的位移 b，阶梯轴的固有频率和模态 12.8 MATLAB应用举例 解决方案：阶梯轴的刚度矩阵和质量矩阵如下所示： 在载荷p3的作用下系统的平衡方程如下所示： b.由下面方程可求得特征值 例12.6 阶梯梁的特征值分析的程序 如图12.12所示的两端固定的阶梯梁，编写一个MATLAB程序，命名为Program17.m，对它进行特征值分析。 Program17.m程序的输入数据如下： xl(i)=阶梯梁的长度 i xi(i)=阶梯梁的转动惯量 i a(i)=阶梯梁的横截面积 i b(i，j)=对应本地的自由度 j 和阶梯梁的自由度 i e=杨氏模量 rho=质量密度 程序给出的阶梯梁的自然频率和模态的输出结果如下： 本章小结： 有限元法是一种流行的数值计算方法，准确的找到解决复杂的实际系统的方案。我们介绍的方法适用于振动问题。我们介绍简单的结构元素，如条，杆，梁，全局坐标系变换的矩阵，装配的单元矩阵，和解决方案的有限元方程推导的刚度和质量矩阵的方法。通过几个静态和