机械振动大作业——简支梁的各情况分析2.docVIP

机械振动大作业——简支梁的各情况分析2 机械振动大作业——简支梁的各情况分析2 PAGE 机械振动大作业——简支梁的各情况分析2 机 械 振 动 大 作 业 姓 名：徐 强 学 号：SX1302106 专 业：航空宇航推进理论与工程 能源与动力学院 2013年12月  简支梁的振动特性分析 题目:针对简支梁、分别用单、双、三、十个自由度以及连续体模型，计算其固有频率、固有振型。单、双、三自由度模型要求理论解；十自由度模型要求使用李兹法、霍尔茨法、矩阵迭代法、雅可比法、子空间迭代法求解基频；连续体要求推导理论解，并通过有限元软件进行数值计算。 解答： 单自由度简支梁的振动特性 如图1，正方形截面（取5mm×5mm）的简支梁，跨长为=1m，质量m沿杆长均匀分布，将其简化为单自由度模型，忽略阻尼，则运动微分方程为,固有频率ωn=，其中k为等效刚度，为等效质量。因此，求出上述两项即可知单自由度简支梁的固有频率。 根据材料力学的结果，由于横向载荷F作用在简支梁中间位置而引起的变形为（）, 为最大挠度，则： == 梁本身的最大动能为： = Tmax=2×= 如果用表示简支梁的质量等效到中间位置时的大小，它的最大动能可表示为： Tmax= 所以质量为m的简支梁，等效到中间位置的全部质量为： 故单自由度简支梁横向振动的固有频率为： ωn== 图1 简支梁的单自由度模型 双自由度简支梁的振动特性 如图2，将简支梁简化为双自由度模型，仍假设在简支梁中间位置作用载荷，根据对称性，等效质量相等,因此只要求出在处的等效质量即可。在至之间积分，利用最大动能进行质量等效，略去小量得： 所以，质量矩阵为： 双自由度简支梁的柔度矩阵： 在b=处作用单位力，挠曲线方程为：则处的变形为：，同理可求：，，其中。 所以，柔度矩阵为： 动力矩阵： 令特征行列式为零，得到频率方程为： 其中，，将上式整理得： 其中，。 解上述方程的根为： , , 由式， 其中，分别将、代入上式，得 第一、二阶主振型分别为： ， 图2 简支梁的双自由度模型 三自由度简支梁的振动特性 如图3，将简支梁简化为三自由度模型，按照双自由度类似的等效思想，可得等效质量： 因此，质量矩阵为： 由机械振动中文教材例可知，系统的柔度矩阵为： 其中，。 动力矩阵： 令特征行列式为零，得到频率方程为： 其中，，将上式整理得： 其中，。 利用Matlab软件,求解上述方程的根为： , , , 由式， 其中，分别将、、代入上式，得 第一、二、三阶主振型分别为： ， , 图3 简支梁的三自由度模型 十自由度简支梁的数值方法 将简支梁简化为十自由度模型（如图4）。 图4 简支梁的十自由度模型 通过在一点施加单位力，计算其余点的挠度，可得柔度矩阵： 表1 十自由度挠度变形矩阵 十自由度简支梁为十个集中质量的振动模型，每个质量都近似等于，因此，质量矩阵为： 动力矩阵为： 下面，用如下几种方法计算十自由度简支梁的固有频率与振型。 1、邓克莱法 利用邓克莱法求基频（比准确值小）： 因此，将柔度矩阵主对角线上各元素相加并乘以，可求得： 2、瑞利法 （1）瑞利第一商 柔度矩阵求逆得刚度矩阵： ，其中，矩阵见表2。 表2 矩阵各元素 假设力作用在简支梁中间位置而得到各点的静变形，可以表示为： 其中，。 因此，可以假设振型： 则由瑞利第一商公式：，可得： （2）瑞利第二商 同样假设力作用在简支梁中间位置，由瑞利第二商公式： 可得: 瑞利法中，代表质量矩阵，代表刚度矩阵，代表柔度矩阵，为模态向量。 3、李兹法 将十自由度简支梁缩减为三自由度，假设振型为： 则可求出： 0 0 18 0 0 由式 ，得： ,, 其中，，因此可得： 以及： 所以系统的前三阶主振型的近似为： × × × 4、矩阵迭代法 单位力作用在简支梁中间位置得到各点的挠度变形，将首项化一，得： 其中，。 因此可假设振型： 利用矩阵迭代求第一阶固有频率和主振型： 由上，仅3次矩阵迭代后，与基本相等，因此可以认为系统的第一阶主振型为： 第一阶固有频率为： 5、雅可比法 根据雅可比