动力机械振动与噪声测量课程设计.docxVIP

《动力机械振动与噪声测量》 课程设计报告书 船舶动力装置轴系扭转振动计算 作者：陈瑞爽 班级：轮机1302班 华中科技大学船舶与海洋工程学院 2016年7月 一：设计目的 船舶轴系振动控制就是要在设计及安装中采取措施，以保证动力装置的振动限制在容许的范围内。在船舶轴系设计阶段应使轴系的振动临界转速尽量不在工作转速范围内，或者不发生激烈的振动，以保证轴系工作正常。因此轴系扭转振动计算是船舶轴系设计时必不可少的项目。 二：已知条件 1.船舶轴系布置草图： 图1 船舶推进轴系布置草图 主机额定功率为1000r/min 2.船舶轴系当量系统图： 图2.轴系当量系统图 1.空气压缩机 2.水泵 3.变速齿轮 4-8.柴油机气缸 9.飞轮 10.减速器 11.联轴节 12.螺旋浆 3. 船舶推进轴系当量系统参数 表1 船舶推进轴系当量系统参数 序号 1 2 3 4~7 8 9 10 11 12 转动惯量 (kg·m2) 5.98 1.08 1.04 2.913 2.913 51.463 0.6 1.115 3.944 扭转刚度×105 (N·m/rad) 9.5 400.2 180 120.78 190.66 0.9 1.5 60.29 三：设计内容： 1.运动方程： 图2所示的扭振系统，第质量的运动方程为 式中——分别为第k质量的扭振角位移、角速度和角加速度； ——第k质量的转动惯量（kg·m2）; ——第k质量的外阻尼系数（ N·m·/rad）； ——第k质量与第k+1质量间轴段的内阻尼系数（N·m·/rad）； ——第k质量与第k+1质量间轴段的刚度（N·m·/rad）； ——作用在第k质量上的激振力矩（N·m）。 由上式可推得整个系统的振动微分方程为 由于要研究的系统为无阻尼自由振动的系统，所以该系统的振动微分方程为： 其中转动惯量矩阵为对角阵，即 刚度矩阵为三对角阵，即 2. 计算轴系的固有频率及相应的主振型 (1):算法 ①：霍尔茨法是求解轴系扭振系统自由振动特性的一种有效的方法。该法通过不断调整试算频率，得到满足所给定边界条件的真值，同时也求得了相应的振型。 对如图2所示的轴系当量扭振系统，从轴段弹性力矩的表达式 可以得出第k质量振幅的表达式 而由第k质量所受力矩的平衡关系可以得到 设A1=1，并由边界条件，根据以上两个递推公式对试算频率有： 如果所选的是系统的一个固有频率，则应有 =0 即 ②：求解的过程可通过剩余力矩曲线图来说明，见图3。图中曲线为随的变化规律，曲线与横轴的交点即为系统的固有频率。 图3 剩余力矩曲线图 具体的试根过程可采用二分法弦截法或切线法等，这里介绍二分法。 设第一次试算频率为，得到的余量为，第二次试算频率取为为试算频率增量。只要每次试算所得扭矩余量与的符号相同，那么下次试算频率就取，如果余量与的符号相反，则下一次试算频率取这样试算频率将收敛于所求的固有频率。实际计算中，当两次得到符号相反的余量的试算频率之差小于允许误差时就以最后一次的试算频率作为固有频率的计算值： 求出一阶固有频率后，以为起算频率可求得二阶固有频率。以此类推，直至所需阶数的固有频率。 (2):计算过程和结果 根据以上方法，可以编程计算出所需要的固有频率和阵型，具体的程序，请老师批阅附带的MATLAB文件。 根据程序，其中计算的结果包括:角频率 各阶频率 剩余力矩误差 阵型矩阵 转速禁区 其计算结果为：（共计算到四阶） 1阶----------- 角频率为 106.4596 频率为 16.9436 剩余力矩的误差值为 0.063473164613541839915424575091151 阵型矩阵为 [ 1.0, 0.92866, 0.92668, 0.92168, 0.9117, 0.89923, 0.88431, 0.86696, 0.85447, -7.3287, -11.906, -11.995] 避开的转速范围为 145.1367894---180.4010514 2阶----------- 角频率为 344.1434 频率为 54.7721 剩余力矩的误差值为 -1.4397945265409839843467823870334 阵型矩阵为 [ 1.0, 0.25449, 0.23597, 0.1932, 0.12395, 0.051146, -0.023114, -0.096715, -0.14159, -0.05927, 0.0182, 0.019729] 避开的转速范围为 478.8702297---571.3575645 3阶----------- 角频率为 644.1149 频率为 102.5141 剩余力矩的误差值为 -1.5012001867784