船舶动力装置轴系扭转振动计算课程设计

1、船舶动力装置轴系扭转振动计算课程设计班级：轮机0801班学号：u200812201姓名： 一 设计任务及意义：在推进装置中，从主机到推进器之间，用传动轴及保证推进装置正常工作所需的全部设备连接在一起的中间机构成为轴系。船舶轴系是船舶动力装置的重要组成部分之一。轴系的工作好坏，将直接影响船舶的推进特性和正常航行，并对船舶主机的正常工作也有直接的影响。如果轴系设计质量欠佳，将会引起机体振动、传动系统零部件损坏、轴承过度磨损、甚至轴件折断等事故，不仅会中止机械系统的正常运行，也会危急工作人员的生命安全。因此对轴系必须进行深入的研究，以利于其正确的设计、制造、安装和检验。船舶轴系振动控制就是设计及安装

2、中采取措施，以保证动力装置的振动限制在容许的范围内。这次设计主要是针对简化实际系统后的理想的轴系当量系统图进行分析，采用其参数，通过各种方法（矩阵特征值特征向量、holzer法、专门解微分方程的软件等）求出系统的各阶频率及其主阵型，通过对着2个参数进行分析，得出所需的数据，并总结归纳出轴运转过程中要注意的问题，以保证轴能够安全有效的运转。二 柴油机推进轴系布置图：图1所选主机的型号为6350zc-1，其额定功率为661kw，额定转速为350r/m。三 轴系当量系统图：为了方便对船舶的推进轴系进行分析和振动计算，将实际的船舶推进轴系简化成当量系统，如下图：图2其中：1.空气压缩机 2.水泵 3.

3、变速齿轮 4-8.柴油机气缸 9.飞轮 10.减速器 11.联轴节 12.螺旋浆 各当量参数如下表：序号1234789101112转动惯量(kgm2)5.981.081.042.9132.91351.4630.61.1153.944扭转刚度10-5(nm/rad)8.2392.2150112.78169.660.50.550.29表1转动惯量与扭转刚度的等效计算原理：a，转动惯量：轴系作扭转振动时，其运动部件可分为旋转运动件和往复式运动件，其中，旋转运动件的转动惯量一般都是对圆盘这类有规则几何形状的物体进行积分：j=.比如真空心圆轴的转动惯量为j=（）l（kgm）。往复运动件的等效转动惯量按照

4、动能守恒的原则进行等效，一般，往复式运动件的等效转动惯量近似等于一半的往复质量集中在曲柄销中心时所具有的转动惯量，即：=.b，扭转刚度：轴段扭转刚度k=g/l， 其中，g为轴段材料剪切弹性模数，l为轴段长度。 所以，决定扭转刚度的三个参数都是轴段的结构与材料参数，对于一根具体的轴来说，刚度是不随运转条件变化的定植。一次，扭转刚度的等效过程比较容易计算。四 轴系转动惯量矩阵、刚度矩阵及振动微分方程：转动惯量矩阵为：刚度矩阵为: 根据表1的当量参数，将数据带入两矩阵，得：转动惯量矩阵：j=刚度矩阵为：k=x系统微分方程为：j+k=f(t)其中，f(t)整个系统受到的激励力。五 用matlab求解微

5、分方程： 我们知道系统的固有特性与系统受到的激励力无关，即忽略f(t)项，则有一下形式：j+k=0两边同乘以，则有 j+ k=0令w= k，有i+w=0, i为单位矩阵，w称为刚度动力矩阵，假设系统自由振动的主阵型为简谐振动，则，设=cos(wt-),假设=，则最后有w=.从上式可以看出系统固有频率就是系统特征值的开二次方，系统的主阵型阵列就是系统的特征向量。按照这个思路，很容易用matlab求解。建立的m文件代码为：function x,a=tttj=diag(5.98,1.08,1.04,2.913 ,2.913,2.913,2.913 ,2.913,51.463,0.6,1.115 ,3

6、.944);m=8.2,400.4,542.2,262.78,225.56,225.56,225.56,282.44,170.16,1,50.79,50.29;n=-8.2,-392.2,-150,-112.78,-112.78,-112.78,-112.78,-169.66,-0.5,-0.5,-50.29;k=diag(m)*105+diag(n,1)*105+diag(n,-1)*105;b=k*inv(j);x,a=eig(b);a=a.0.5;end 通过上面的程序可以求出系统的12阶固有频率与相对应的主阵型，但我们只要求前4阶，为了更好地对应和观察，将各将前4阶固有频率及其主阵型列

7、成表2：阶 次1234固有频率11.4051.4765.94107.69主振型11.000001.000001.00001.0000020.962570.23743-0.25168-2.3389730.961650.22081-0.27666-2.3792940.958910.17574-0.33868-2.4091850.953980.11104-0.40616-2.1640460.947790.04334-0.45563-1.6629770.94036-0.02553-0.48490-0.9625080.93167-0.09370-0.49268-0.1533890.92507-0.137

8、34-0.483330.3983710-6.20064-0.1630988.07290-0.1145111-12.94440.01581-4.767940.00174712-12.99680.017722-5.509550.0272527表2六 运用holzer法求各主阵型：在此，不对进行试根，即：不通过二分法求解个阶固有频率。仅仅利用上面求出的4阶固有频率（即：71.644，323.37,414.29,676.65），通过迭代求出各主阵型: 从轴段弹性力矩的表达式可以得出第k质量振幅的表达式而由第k质量所受力矩的平衡关系可以得到设a1=1，并由边界条件，根据以上两个递推公式对试算频率有：如果

9、所选的是系统的一个固有频率，则应有即 按照上面原理编写matlab的m文件：function a=tdba=1:12;u=0:11;j=5.98,1.08,1.04,2.913 ,2.913,2.913,2.913 ,2.913,51.463,0.6,1.115,3.944;k=8.2 392.2 150 112.78 112.78 112.78 112.78 169.66 0.5 0.5 50.29*105;w=input(w=);for i=1:11u(i+1)=u(i)+j(i)*w2*a(i);a(i+1)=a(i)-u(i+1)/k(i);endend输入（input）各频率：71.

10、644，323.37,414.29,676.65；计算结果如表3：阶 次1234固有频率71.644323.37414.29676.65主振型11.00001.00001.00001.000020.96260.2374-0.2517-2.339030.96160.2208-0.2767-2.379340.95890.1757-0.3387-2.409250.95400.1110-0.4062-2.164160.94780.0433-0.4556-1.663070.9404-0.0255-0.4849-0.965280.9317-0.0937-0.4927-0.153490.9251-0.137

11、4-0.48330.398410-6.2007-0.160588.0747-0.122311-12.94450.0178-4.76900.028812-12.99680.0192-5.51060.0274表3从表2和表3的结果来看，数据相当吻合的，证明了结果的正确！七 画出主阵型、计算“转速禁区”： 用matlab画出4阶谐振频率下轴系扭转振动得到主阵型：=71.644=323.37=414.29=676.65“转速禁区”计算如表4（其中）：(rad/s)(r/min)r=/“转速禁区”：16/（18-r）-（18-r）71.644684.201.90679.95688.48r/min323.

12、373088.188.58 1818.175245.32 r/min414.293956.4710.99 1733.439030.46 r/min676.656462.0117.95 20.192067840 r/min表4可以看出，系统的频率越大，“转速禁区”越大，对系统的限制越大。八 小结：这次计算通过对轴系当量系统图进行分析，由其转动惯量和扭转刚度2参数列出微分方程，在用课本的转化方法，将所求频率和主阵型与矩阵的特征值、特征向量关联起来，最后用matlab轻易的求出频率和主阵型。在课设过程中我很好的熟悉了matlab软件，并由逐渐摸索转为较为熟练的运用，这是我最大的收获。希望能把这次经验总结积累起来，努力运用到今后的学习工作之中去。参考文献：1