船舶推进轴系纵向振动共振转换器的优化设计胡泽超;何琳;徐伟;李正民;赵兴乾【摘要】[目的]在推力轴承上集成共振转换器(RC)可以改变轴系纵向振动的传递路径,衰减传递到基座的响应使轴系的固有频率避开螺旋桨叶频及其倍叶频激励力,从而实现减振、调频的目的.[方法]为此,建立推进轴系纵向振动的力学模型,基于传递矩阵法计算桨轴系统的振动响应,以力传递率为指标,分析RC的主要参数对推进轴系隔振效果的影响,分别采用最大值最小化方法和曲线面积最小的参数修正方法,对RC的主要参数进行优化设计.[结果]研究结果表明:加装RC后,轴系的隔振效果得到了明显的改善,采用曲线面积最小修正的优化设计方法可使RC的减振调频效果更佳.[结论]通过对RC结构参数的合理设计能使减振系统获得优良的隔振效果.【期刊名称】《中国舰船研究》【年(卷),期】2019(014)001【总页数】7页(P107-113)【关键词】共振转换器;推进轴系;传递矩阵法;纵向振动【作者】胡泽超;何琳;徐伟;李正民;赵兴乾【作者单位】海军工程大学振动噪声研究所,湖北武汉430033;船舶振动噪声重点实验室,湖北武汉430033;海军工程大学振动噪声研究所,湖北武汉430033;船舶振动噪声重点实验室,湖北武汉430033;海军工程大学振动噪声研究所,湖北武汉430033;船舶振动噪声重点实验室,湖北武汉430033;海军工程大学振动噪声研究所,湖北武汉430033;船舶振动噪声重点实验室,湖北武汉430033;海军工程大学振动噪声研究所,湖北武汉430033;船舶振动噪声重点实验室,湖北武汉430033【正文语种】中文【中图分类】U664.20 引言在不均匀伴流场下，螺旋桨的周期性运转产生的脉动激励力是船舶在中高速航行时的主要噪声源。
该纵向激励力通过推力轴承传递至船体，会引起轴系及船体的振动，影响船舶的运行安全，降低船体的声学性能。
为了减小纵向激励力向船体的传递，可在轴系上安装减振器。
考虑到推力轴承传递大推力、小位移的特性，需要设计一种刚度低、阻力大的隔振装置。
液压减振装置利用流体的可压缩性来调整系统的刚度和阻尼，通过合理的设计，可使推进轴系的动态特性满足指标要求。
在推力轴承处安装共振转换器（Resonance Changer，RC）不仅能调节推进轴系的纵向固有频率以偏离螺旋桨的脉动激励频率，还能降低船体艉部的纵向振动响应，达到实现隔振的目的。
Goodwin［1］认为，RC可等效为一种质量―弹簧―阻尼单元，据此设计了一种能在特定频段内降低轴系纵向振动的液压减振装置。
Dylejko等［2］和Li等［3］利用传递矩阵法建立了桨―轴―艇体系统的数学模型，分析了RC的主要参数对桨轴系统力传递率的影响。
李良伟等［4］和王珺等［5］运用动力谐调消振理论对RC进行了优化设计，得到了RC的最优固有频率比和阻尼比，但分析模型较为简单，应用范围有限。
本文拟建立桨轴系统纵向振动的力学模型，采用传递矩阵法，计算螺旋桨激励力传递到壳体的振动响应；以力传递率为指标，分析RC的活塞缸直径d0，连接管长度l1和直径d1及油箱体积V1的变化对桨轴系统隔振效果的影响，分别采用力传递率最大值最小化方法和力传递率与坐标轴围成的面积最小修正法，对RC的主要参数进行优化设计。
1 RC的动力学模型RC由充满油液的油箱、外接管系和活塞缸组成，装置内的工作流体可以改变轴系的纵向刚度和阻尼［6］。
图1为RC的原理模型，其中，P为活塞两侧的压力差，x0和x1为活塞缸两端的位移。
图1 RC结构简图Fig.1 Structure diagram of RC为了便于推导RC的动力学方程，需作出如下假设［7］：1）油箱壁是刚性的，流体的压缩全部发生在油箱内；2）连接管内的流体处于层流状态；3）管内的流体可视为集中质量；4）流体的水力有效长度等于连接管的实际长度；5）不考虑管道中的压缩效应。
根据假设条件，由D'Alembert原理可知，活塞缸内作用于连接管上的压力等于连接管中油受到的惯性力、连接管内的阻尼力以及压缩油箱内的油所需力之和，则RC的动力学方程可描述为式中：分别为活塞缸与连接管的截面积；B1为油的体积模量；μ1和ρ1分别为油的粘度和密度；0,0和1,1分别为活塞缸两端的速度、加速度。
将式（1）两端同乘A0/A1，得其中，Mh，Kh和Ch分别为RC的质量、刚度和阻尼，则可将式（2）转换成质量―弹簧―阻尼的数学模型：式中，F0为活塞受到的外力。
2 桨轴系统纵振数学模型桨轴系统纵向振动的力学模型如图2所示，模型可分解为5个子系统，每个子系统均可以用传递矩阵来表示元件左右两端纵向振动的传递关系。
图2中，下标p，t，c，b，h分别为螺旋桨、推力盘、联轴器、基座和RC，Mp，Mt，Mb，Mc分别为螺旋桨、推力盘、基座、联轴器的质量，K0和Kb分别为油膜和基座的刚度，C0为油膜的阻尼，Ls和Lse分别为艉轴的实际长度和有效长度，L为中间轴的长度，从螺旋桨至联轴器分别为1～5号单元，Ti（i=1，2，3，4，5）为其对应的单元传递矩阵，分别代表j单元左、右端面的位移响应和力响应，下标 j可用b，h，c，t，p替代。
图2 桨轴系统纵向振动力学模型Fig.2 Mechanical model of longitudinal vibration of propeller shafting2.1 轴系纵振的点传递矩阵考虑螺旋桨脉动激励力Fp的作用，桨轴系统的纵向传递矩阵应改写为T3×3的形式，式（4）～式（7）为桨轴系统各子系统的传递矩阵：式中：k=ω c，为轴（尾轴和中间轴的材料属性及截面积均相同）的纵向波数，其中ω为角频率,，为轴的纵向波速，ρ为轴的密度；E和A分别为艉轴和中间轴的弹性模量和截面积。
螺旋桨和联轴器可视为集中质量块。
由于艉轴较长，计算时一般需考虑其有效长度，相应的传递矩阵为T2。
集成RC的推力轴承可进一步分解为3个单元：推力盘、油膜和RC。
推力盘右端面至船体的传递矩阵方程可表示为式（8）：与基座连接的壳体刚度较大，可视为桨轴系统的刚性边界条件将其代入式（8），推导出的关系：式中，Ke为推力盘与壳体之间的等效刚度。
根据式（9），可将推力轴承—基座—壳体模型简化为如图3所示的等效力学模型，则与集成RC的推力轴承对应的传递矩阵T3可简化为式（10）：图3 推力轴承和基座的等效力学模型Fig.3 Equivalent mechanical model of thrust bearing and base2.2 轴系纵振响应计算在得到T1～T5各单元的传递矩阵后，依据边界条件，可对桨轴系统的振动响应进行求解。
式（11）为螺旋桨输入端至联轴器输出端的传递矩阵，其中 Tr，q表示T的第r行q列元素（1≤r，q≤3）。
当脉动激励力 Fp作用在螺旋桨上时，螺旋桨的左端面与联轴器的右端面可视为自由边界条件，将约束条件代入式（11），得式（13）和式（14）为螺旋桨输入端由基座至壳体的传递矩阵T′，其中T′r，q表示T′的第 r行 q列元素（1≤r，q≤3）。
将式（14）代入式（9），得到激励力通过桨轴系统传递到壳体的响应力Fb，由此可计算出桨轴系统的力传递率Tf=Fb/Fp。
3 RC结构的优化设计方法为了使桨轴系统具有良好的隔振效果，需要对RC进行参数设计。
合理的结构参数设计能使RC装置吸收桨轴系统的大部分振动能量，可将力传递率限制在一定范围内，使共振峰值较小；还能调节桨轴系统的固有频率从而避开螺旋桨叶频激励，实现减振和调频的目的。
本文以RC的l1，d0，d1和V1参数为设计目标，分别采用力传递率最大值最小化方法［9］和力传递率曲线面积最小修正法［10］，对RC的结构进行优化设计。
3.1 最大值最小化优化方法RC结构参数优化可表述为力传递率曲线峰值最大值最小化的设计问题，即通过某种算法，搜寻一组设计变量 (l1，d0，d1，V1)，使桨轴系统力传递率的多个峰值在分析频率范围内取最小值。
该优化方法可描述为式（16）：式中：a为纵向模态阶数；m为分析频带内力传递率峰值的个数。
最大值最小化的优化方法更关注力传递率峰值的变化情况，适用于对传递到壳体的力的幅值或对轴系强度有严格限制的约束条件的优化问题，因此，实质上是一种局部优化方法。
3.2 面积最小修正法对于每一组设计变量 (l1，d0，d1，V1)，力传递率曲线 Tf(l1，d0，d1，V1)与坐标轴 f 围成的面积记为，该方法可描述为即搜索在整个取值区间内的最小值所对应的 (l1，d0，d1，V1)。
该方法优化的是整个分析频带内的力传递率曲线，尽管曲线的局部峰值可能会有所增大，但能够在分析频段内改善桨轴系统的整体隔振效果，故是全局优化方法。
4 算例分析以某型船舶桨轴系统为例，计算的参数为：Mp=7 000 kg，Mc=1 000 kg，Mb=4 000 kg，Mt=500 kg，ρ1=860 kg/m3，ρ=7 850 kg/m3，E=200 GPa，A=0.02 m2，Ls=14.6 m，Lse=14 m，L=2 m，Kb=5×109N/m，B1=1.38 GPa，μ1=0.23 Pa·s，螺旋桨桨叶数m=7。
根据实验数据，油膜刚度K0和阻尼C0与转速n和螺旋桨所受载荷F有关。
当F=200 kN时，K0和C0随转速的变化曲线如图4所示。
取n=220 r/min（对应叶频为25.7 Hz），对应的K0=1.4×1010N/m，C0=6.5×108N·s/m。
初步设计的RC参数为l1=1 m，d0=0.06 m，d1=0.01 m，V1=1.6 L。
图4 油膜刚度和阻尼随转速的变化曲线Fig.4 Variation curve of oil film stiffness and damping with rotating speedRC的l1，d0，d1，V1参数对桨轴系统隔振效果的影响如图5所示。
由图可知，未加RC时，桨轴系统在30和132 Hz附近出现系统前两阶共振峰，且第1阶共振峰值较大。
与未加装RC系统的力传递率曲线对比，加装RC后桨轴系统的第1阶模态峰值转化成了两阶峰值更小的低阶模态，分别在未加装RC的第1阶模态频率两侧。
显然，在0～200 Hz分析频带内，RC对桨轴系统的前两阶模态影响较大。
图5 RC主要参数对桨轴系统隔振效果的影响Fig.5 Influence of RC's main parameters on vibration isolation effect of propeller shafting system根据RC的 Mh，Kh，Ch和ωh频率随参数l1，d0，d1，V1的变化趋势，可得出以下结论：1）随着l1增大，RC等效阻尼Ch增大，桨轴系统力传递率的峰值相应减小。