螺旋桨船舶振动辐射噪声分析

《螺旋桨船舶振动辐射噪声分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、推进器激励船舶振动辐射声计算方法谢基榕徐利刚沈顺根吴有生摘要；文章在铁木辛柯梁的基础上，建立螺旋桨、轴系及船体耦合振动的数学物理模型，采用声极子模型计算各船体截面振动引起的辐射声压。计算分析了激励力直接作用到船体和经螺旋桨、轴系激励船体振动产生辐射的差别，指出轴系振动对螺旋桨激励力激起船体振动并产生声辐射有着关键的影响作用。（中国船舶科学研究中心，江苏无锡21408）一,引言二,艇体低频振动计算方法三,辐射声计算方法四,计算结果分析五,讨论六,结论引言潜艇螺旋桨工作于非均匀艉流场中，除了提供动力外还产生不定常激励力作用到船艇结构上，是引起艇体振动的一个重要激励源。螺旋桨通过轴系

2、激励艇体振动引起辐射噪声是潜艇水下航行辐射噪声的重要组成部分，国外很早就在这方面研究了并建立了计算方法和控制措施。本文针对螺旋桨低频宽带力谱的频率特性，在艇体梁模型基础上建立了简化的适合工程设计应用的螺旋桨非定常力经轴系激励艇体振动辐射噪声计算方法，在满足工程计算精度的条件下，极大的简化了辐射噪声计算方法。艇体低频振动辐射声计算方法对于潜艇螺旋桨、轴系及艇体耦合振动系统建立上图所示的数学物理模型，该模型包括两大部分；船体梁模型和轴系子系统，支撑轴承和推力轴承作为轴系结构的约束条件划归为轴系子系统。由于艉流作用在螺旋桨上的合力具有低频特性，而在低频0~f(f为舱段周向的首阶共振频

3、率)范围内，采用艇体梁模型描述艇体振动便可得满足工程要求的解。在艇体梁模型中，船体结构被当做截面不发生变形的变截面梁结构。由于实际船体在各截面处的弯曲刚度中心与质心不重合，因此，垂直面内的弯曲振动与纵向振动耦合，而水平面内的弯曲振动与扭转振动耦合，而从，需对这两组无法解耦的振动进行同时解耦。在推导艇体体梁的刚度矩阵和质量矩阵时，把参考坐标系的原点建立在刚度中心，因此梁单元做垂向弯曲运动时，梁的单位截面的转角位移将产生明显的纵向位移，如上图。由于截面转角位移引起的纵向位移是微幅运动，所以引起的横向位移可以忽略，只考虑纵向位移。于是梁单位的实际纵向位移可以表示为梁单位的实际横向位移

4、可以表示为应用结构有限元理论，分别将船体梁垂向弯曲和纵向耦合振动、水平弯曲和扭转耦合的振动的单元位移向量记为并将梁单位元内部的横向位移和转角位移表达成梁单元节点的二次多项式关系，将纵向位移表示成节点位移的线性关系，从而可以得到任意节点位移下的梁单元应变势能和动能，进而得到梁单元的刚度和质量矩阵。以垂向弯曲与纵向振动耦合为例，梁的单元应变势能与动能可分别表示为；由上两式并根据矩阵知识转变可得梁单位元作用于垂向弯曲与纵向耦合振动时的刚度矩阵和质量矩阵。在低频段，连接轴系结构与艇体之间的轴承装置可以使用等效的阻尼弹簧模拟。对于径向支撑轴承，考虑横向刚度，在轴系与艇体间起弯曲振动的耦合

5、作用；主推力轴承同时具有横向刚度和纵向刚度，其在轴系与船体间起到弯曲振动及纵向振动耦合作用。将船艏到艉将船体梁单元和轴系梁单元的节点逐一编号，并将船体梁单元和轴系梁单元的刚度、质量矩阵填入对应自由度，从而组装形成总刚度阵和总质量阵，针对分析频率计算得广义刚度矩阵得式；再将用于描述各轴承动力性能的阻抗参数换算到等效刚度并分别叠加到广义刚度矩阵中，即可得到描述螺旋桨、轴系系统和船体耦合振动的结构有限元刚度矩阵。再考虑流场作用在螺旋桨上的激励力，从而得到了包括螺旋桨、轴系系统和船体结构在内的全系统振动方程；完成方程组求解后即可得到螺旋桨、轴系与船舶结构耦合振动的位移解。艇体梁模型的辐

6、射声计算方法针对推进器、轴系及艇体耦合振动具有低频特性，此频率段内的声波波长远大于船体截面的尺度，因此，无论是弯曲振动还是纵向振动，艇体梁单元的声辐射都满足致密条件，因而可以使用声极子模型计算各艇体截面振动引起的辐射声压。船体水下振动声辐射在辐射机理上可以分解为两部分，并分别使用不同声辐射模型计算。船体弯曲振动引起船体段元在水介质中横向平动且不产生船体截面变化，因此可用偶极模型计算各梁单元的引起辐射声压。声压公式为；当船体做纵向振动时，纵向运动经泊松效应引起截面产生体积变化，截面变化部位(含艏、艉端部分)的纵向运动也将产生局部的体积变化，见下图所示其中a（x）为各截面处的艇体平

7、均半径，阴影面积部分表示艇体纵向振动引起的局部体积变化。因此，可用单极子模型计算因船体纵向振动引起的水下辐射声压级，三部分体积变化产生的辐射噪声总和为；式中，S,u分别为梁单元的横截面积和纵向弯曲横向位移，0为声场点与艇轴线的夹角，纵向振动产生声压与声场点在艇体横截面上的方位角无关。当不同激励力间不存在固定相位差时，其振动响应及后续的声压也不具有固定的相位差，因此，多激励引起的辐射声场也不存在稳定的干涉作用，从而适合使用能量叠加方法计算由两方向振动引起的总辐射声压级。在本文中，由于考虑了因质