《基于解析法的船舶推进轴系纵向强迫振动计算方法.pdf》由会员分享，可在线阅读，更多相关《基于解析法的船舶推进轴系纵向强迫振动计算方法.pdf（8页完整版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 103177171 A (43)申请公布日 2013.06.26 CN 103177171 A *CN103177171A* (21)申请号 201110433478.X (22)申请日 2011.12.22 G06F 19/00(2006.01) (71)申请人 中国舰船研究设计中心 地址 430064 湖北省武汉市武昌区紫阳路 268 号 05 分箱 (72)发明人 周瑞 管文生 杨义顺 王丹 彭戈 (74)专利代理机构 上海航天局专利中心 31107 代理人 冯和纯 (54) 发明名称 基于解析法的船舶推进轴系纵向强迫振动计 算方法 (57) 摘要 本发明公开了。
2、基于解析法的船舶推进轴系纵 向强迫振动计算方法。 根据原始数据输入， 构建系 统的质量及刚度矩阵 ； 计算系统的激振力和质量 阻尼系数， 构建激振力列向量和质量阻尼系数矩 阵 ； 计算轴段的滞后阻尼系数， 修正系统刚度矩 阵 ； 求解微分方程， 计算推进轴系纵向强迫振动 响应。本发明方法解决了采用解析法进行推进轴 系纵振响应计算时轴段滞后阻尼系数的处理和计 算问题， 运算简便、 易于实现， 提高了计算精度。 (51)Int.Cl. 权利要求书 1 页 说明书 4 页 附图 2 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书1页 说明书4页 附图2页 (10)申请公布。
3、号 CN 103177171 A CN 103177171 A *CN103177171A* 1/1 页 2 1. 基于解析法的船舶推进轴系纵向强迫振动计算方法， 其特征在于包括以下步骤 ： （a） 根据原始数据输入， 构建系统的质量及刚度矩阵 ； （b） 计算系统的激振力和质量阻尼系数， 构建激振力列向量和质量阻尼系数矩阵 ； （c） 计算轴段的滞后阻尼系数， 修正系统刚度矩阵 ； （d） 求解微分方程， 计算推进轴系纵向强迫振动响应。 2. 根据权利要求 1 所述的基于解析法的船舶推进轴系纵向强迫振动计算方法， 其特征 在于 ： 所述步骤 1-c） 中， 计算轴段的滞后阻尼系数， 修正系统。
4、刚度矩阵包括以下步骤 ： （a） 计算轴段的滞后阻尼功 轴段纵振时， 单位体积材料在一个振动循环中的阻尼功为 ： 式中 ：为纵振时轴段的交变拉压应力 整个轴段在一个振动循环中的滞后阻尼功为 ： 式中 ：A和L分别为轴段的横截面积和长度 根据和， 代入上式得到 式中 ：F为纵振交变推力，K为轴段刚度，E为材料弹性模量 ； （b） 计算轴段的滞后阻尼系数 轴段的滞后阻尼系数按照复刚度处理， 即在考虑阻尼力时， 在刚度系数项前乘以复常 数iV，V定义为复阻尼系数， 轴段滞后阻尼系数CZ=iKV， 则轴段在一个振动循环中的滞后阻 尼功为 ： 令WA=WZ， 可以得到， 则轴段的滞后阻尼系数CZ为 ： 。
5、（c） 修正系统刚度矩阵 根据计算得到的轴段滞后阻尼系数CZ， 对系统刚度矩阵K进行如下修正 。 权 利 要 求 书 CN 103177171 A 2 1/4 页 3 基于解析法的船舶推进轴系纵向强迫振动计算方法 技术领域 0001 本发明涉及船舶推进系统技术领域， 具体涉及一种基于解析法的船舶推进轴系纵 向强迫振动计算方法。 背景技术 0002 船舶推进轴系是一个弹性系统， 因此在推进主机产生的径向力以及在不规则伴流 中螺旋桨产生的周期性变化的纵向推力作用下， 推进轴系将产生纵向振动现象。当轴向激 振力频率与纵向振动固有频率相同时， 轴系会发生纵向共振。此时， 推力轴承、 传动齿轮及 推进主。
6、机都将承受附加的轴向交变负荷和弯曲应力， 极易引起疲劳破坏。 因此， 在船舶推进 轴系设计时， 需要进行纵振计算。 0003 船舶推进轴系纵振计算分为自由和强迫振动两个方面内容， 前者计算轴系的固有 频率和固有振型， 后者计算轴系在激振力作用下的振幅和交变推力。轴系纵向自由振动计 算通常采用离散化的集总参数模型， 采用的计算方法包括解析法、 有限元法和传递矩阵法 等， 利用这些方法可以精确地对推进轴系的固有频率和振型进行计算分析。纵向强迫振动 计算采用的模型与纵向自由振动相同， 但是需要考虑系统中激励和阻尼的影响。 0004 现有的船舶轴系纵振计算指导性技术文件 （CB*/Z 337-1984。
7、） 对系统中的激励和 阻尼均以功的形式给出， 即主机激励功、 螺旋桨激励功、 轴承阻尼功、 螺旋桨阻尼功和轴段 滞后阻尼功等， 没有明确地给出激振力和阻尼系数的计算公式， 因此只能采用近似算法， 如 能量法或放大系数法， 对强迫振动进行分析， 而不能采用解析法、 有限元法和传递矩阵法等 精度更高的计算方法。 0005 近年来， 国内相关学者给出了主机和螺旋桨激振力以及轴承和螺旋桨阻尼系数的 计算公式， 比如参考文献 船舶推进轴系振动 （陈之炎， 上海交通大学出版社， 1987） 和 动 力装置振动数值计算 （张志华， 哈尔滨工程大学出版社， 2007） 。但是， 关于轴段滞后阻尼 （内阻尼） 。
8、系数的计算和处理方式尚未见报道。 发明内容 0006 本发明的目的在于针对现有技术存在的不足， 提供基于解析法的船舶推进轴系纵 向强迫振动计算方法。通过推导纵向振动轴段的滞后阻尼系数， 使得解析法可应用于推进 轴系纵向强迫振动响应的计算， 从而提高了计算结果的精度。 0007 为达到上述目的， 本发明是通过以下的技术方案实现的。基于解析法的船舶推进 轴系纵向强迫振动计算方法， 包括如下步骤 ： 步骤一、 根据原始数据输入， 构建系统的质量及刚度矩阵 ； 步骤二、 计算系统的激振力和质量阻尼系数， 构建激振力列向量和质量阻尼系数矩 阵 ； 步骤三、 计算轴段的滞后阻尼系数， 修正系统刚度矩阵 ；。
9、 步骤四、 求解微分方程， 计算推进轴系纵向强迫振动响应。 说 明 书 CN 103177171 A 3 2/4 页 4 0008 所述的步骤三包括如下步骤 ： 1、 计算轴段的滞后阻尼功 ； 2、 计算轴段的滞后阻尼系数 ； 3、 修正系统刚度矩阵。 0009 本发明提供的基于解析法的船舶推进轴系纵向强迫振动计算方法运算简便， 易于 实现， 解决了采用解析法进行推进轴系纵振响应计算时轴段滞后阻尼系数的处理问题， 提 高了计算精度。 附图说明 0010 图 1 为本发明方法的程序流程图。 0011 图 2 为本发明一个实施例的推进轴系纵向振动计算模型。 0012 图 3 为轴系中推力轴承处纵向。
10、强迫振动的位移幅值。 0013 图 4 为轴系中推力轴承处纵向强迫振动的交变推力幅值。 具体实施方式 0014 以下将结合附图和一个具体实施例对本发明作详细说明。应该强调的是， 下述说 明仅仅是示例性的， 而不是为了限制本发明的范围及其应用。 0015 如图 1、 2、 3 和 4 所示， 本发明的具体方法步骤如下 ： 第一步 ： 根据原始数据输入， 构建系统的质量及刚度矩阵 本实施例中， 轴系纵振离散化集总参数模型如图 2 所示， 包括从螺旋桨到推力轴承之 间的轴段。其中，M1为螺旋桨的集中质量，M2M9为轴段简化后的集中质量，K1K8为轴段的 纵向刚度，K9为推力轴承及其基座的合成刚度。 。
11、0016 考虑到计算误差， 推力轴承的基座刚度分 6 档进行计算， 即K9分别取 6 个不同的 刚度值。 式中 ：n为集中质量个数。 0017 第二步 ： 计算系统的激振力和质量阻尼系数， 构建激振力列向量和质量阻尼系数 矩阵。 0018 本实施例中， 只有螺旋桨为激励源， 其纵向激振力计算公式为 ： 说 明 书 CN 103177171 A 4 3/4 页 5 式中 ：为叶频及其倍叶频， 为旋转角速度，k为正整数，为桨叶数 ； 为额定转速下的推力 ； 为相位角， 共振时取； 为推力变化系数。 0019 系统的激振力列向量构建为 ： 本实施例中， 螺旋桨质量阻尼系数的计算公式为 ： 式中 ： 。
12、为水的密度 ； 为旋转角速度 ；为螺旋桨直径 ；为螺旋桨盘面比。 0020 本实施例中， 轴段质量阻尼系数的计算公式为 ： 式中 ：Mi为第i个轴段集中质量， 本实施例中，i=29。 0021 系统的质量阻尼系数矩阵构建为 ： 第三步 ： 计算轴段的滞后阻尼系数， 修正系统刚度矩阵 3a. 计算轴段的滞后阻尼功 轴段纵振时， 单位体积材料在一个振动循环中的阻尼功为 ： 式中 ：为纵振时轴段的交变拉压应力 整个轴段在一个振动循环中的滞后阻尼功为 ： 式中 ：A和L分别为轴段的横截面积和长度 根据和， 代入上式得到 式中 ：F为纵振交变推力，K为轴段刚度，E为材料弹性模量 ； 3b. 计算轴段的滞。
13、后阻尼系数 轴段的滞后阻尼系数按照复刚度处理， 即在考虑阻尼力时， 在刚度系数项前乘以复常 数iV，V定义为复阻尼系数， 轴段滞后阻尼系数CZ=iKV， 则轴段在一个振动循环中的滞后阻 尼功为 ： 说 明 书 CN 103177171 A 5 4/4 页 6 令WA=WZ， 可以得到， 则轴段的滞后阻尼系数CZ为 ： 3c. 修正系统刚度矩阵 根据计算得到的轴段滞后阻尼系数CZ， 对系统刚度矩阵K进行如下修正 第四步 ： 求解微分方程， 计算推进轴系纵向强迫振动响应 构造矩阵H， 使得 ： 式中 ：F * 为激振力的幅值，A * 为纵向振动复幅值，H=K- 2M+iC。 0022 设定了转速ne和谐次v后， 则直接可以求解出纵向强迫振动的复振幅A *， 即 式中 ：=2vne 纵向强迫振动的振幅A为A * 的模值， 第i轴段纵向强迫振动交变推力的幅值为 Ki(Ai+1-Ai)。 0023 在本实施例中， 采用本发明方法对图 2 所示的推进轴系进行叶片次纵向强迫振动 响应分析， 计算结果如图 3、 4 所示。 说 明 书 CN 103177171 A 6 1/2 页 7 图 1 图 2 图 3 说 明 书 附 图 CN 103177171 A 7 2/2 页 8 图 4 说 明 书 附 图 CN 103177171 A 8 。