第17卷第4期 2013年4月 文章编号:1007—7294(2013)04—0430—09 船舶力学 Journal of Ship Mechanics V01.17 No.4 Apr.2013 基于水弹性力学的SWATH船结构 振动与噪声分析 叶永林,吴有生,邹明松,倪其军 (中国船舶科学研究中心,江苏无锡214082) 摘要:文章以三维水弹性力学理论为基础,提出了综合考虑波浪与设备机械激励下的船体结构振动分析方法,并 对小水线面双体船的船体振动响应特性进行了分析,预报了该船体的固有频率特性,及在试验工况与航行时受发 电机激励下的局部结构振动响应,并与实船测量结果进行了对比。结果表明,该方法对于船体总振动、及船体受机 械激励作用下的强迫振动分析预报是有效的,可应用于工程设计。同时,利用水弹性方法对SwATH船在波浪中 的辐射噪声与近场噪声进行了预报,数值预报结果与试验进行了对比,表明该方法在噪声分析中也是可行的。 关键词:水弹性力学;小水线面双体船;机械激励;振动响应;辐射噪声 中图分类号:U661.44 Tv131.2 文献标识码:A doi:10.39690.issn.1007—7294.2013.04.013 Analysis of the structural vibration and noise radiation of a SWATH ship based on hydroelastic method YE Yong-lin. You-sheng,ZOU Ming-song,NI Qi-jun (China Ship Scientiifc Research Center,Wuxi 214082,China) Abstract:The structural vibration and noise radiation of the ship in wave is achieved synthetically to con— sider the wave force and mechanical exciting force based on the 3D hydroelasticity method.The structural vibration response characteristics of SWATH ship are analyzed.The naturl farequency characteristics of the hull and the local structural vibration response excited by diesel generator are predicted when in 0 kn and in navigation.The predicted results are compared with the test resuhs of the ship.and illustrate that the analysis method is effective for the prediction of overal1 vibration and forced vibration of the hull excited by mechanical force.The near and far field underwater noise radiations of the SWATH ship are predicted by the hydroelasticity method also.The predictions are compared with the measured results.showing that the present method is also rational for the noise radiation analysis. Key words:hydroelasticity;SWATH ship;mechanical excitation;vibration response;radiation noise 1引 言 小水线面双体船(SWATH)具有独特的水动力船型,其耐波性能非常优秀,适合作为远洋科学考 察的平台。在海洋环境中航行时,波浪和船上机械设备均会引起船体振动,包括总体振动与局部振 动,需要引起足够的重视。过大的振动会导致局部结构的疲劳失效、设备故障,以及影响乘员的舒适 收稿日期:2012—12—08 作者简介:叶永林(1974-),男,中国船舶科学研究中心高级工程师;吴有生(1942一),男,中国船舶科学研究 中心研究员,博士生导师,中国工程院院士。 第4期 叶永林等:基于水弹性力学的SWATH船结构…431 性[】1。且如果作为海洋科学考察船,过大的振动会影响船上的试验,同时产生的结构噪声,也会影响 到其声学测量系统的性能。 水弹性力学分析方法l2- ,是把流场和结构的运动与变形作为一个完整的系统来进行分析。可应 用解决的问题涉及各类船舶结构在波浪激励、砰击与甲板上浪下的稳态和瞬态响应;海洋工程结构 物的运动、外载荷、承载能力、波激振动、疲劳、可靠性与安全性的评估l4_5 ;水中结构在流场中的振动 与噪声辐射 ]、水翼的颤振、水下爆炸及快速运动结构出入水的瞬时冲击响应;各类薄壁储液结构和 容器的晃荡、管道或管群的涡激振动等。 本文采用线性三维水弹性理论,对波浪中船体在机械激励下的强迫振动特性与水下噪声进行了 分析。对于结构振动,文献[7]尝试了采用水弹性力学方法进行计算零航速工况时的局部振动。分析 可以采用文献[8]所介绍的软件。本文预报了航行工况的振动响应,包括下潜体与工作甲板的结构振 动响应。结果与实船测试结果进行了比较,结果表明该预报方法合理、可行。同时,利用声介质中的 三维船舶水弹性分析方法,对该船的水下噪声进行了分析,其中6 kns航速的分析结果与测试值进 行了对比,结果比较吻合。说明该方法可适用于水面船的噪声预报,具有较强的工程性。同时,对零 航速时的近场噪声与辐射噪声也进行了预报,以考察其在进行相关试验时的声学噪声环境。 2三维水弹性力学理论 假定船舶在波浪中作微幅运动和变形,根据模态叠加的原理,船体结构相对其平衡位置的运动 和变形可以表达为: =∑H-q),p,:∑( , , ,)p r=l r=l (1) 其中:p,(r=l,2,…,171,)表示相对于第r阶船体干模态位移U,的主坐标分量。 (1)式中的主坐标p 满足广义线性水弹性力学运动方程: [口 ]{b}+Ib+8] )+[c+c] }={F}+{ )+{△)+{G) (2) 式中: )={P (t),P (t),…,P (t))为广义主坐标列阵,[a]、[b]和[C]分别为结构广义质量矩阵、广 义阻尼矩阵和广义刚度矩阵;[A]、[ ]和[C]分别为广义流体附加质量矩阵、广义流体附加阻尼矩 阵和广义流体恢复力矩阵;fF}、{ }、{△}和fG)为广义波浪激励力、广义静态力、广义集中力和广义 体积力列阵。 其中,广义集中激励力{△)可以表示为: f△)={A。,A:,…,h } 当外部激励力为简谐力时,有 △,(3) ( ):Zf/(t)T~, e 一 (4) 3 SWATH船及主要激励源 3.1船体主要参数 该SWATH船的图片见图1,其主要技术参数见表1。 3.2船上主要激励频率 该船采用综合电力推进系统,螺旋桨及主机柴油机、推进电机的激励频率为: a.螺旋桨轴频率:250 r/min(全速)、163 r/m(巡航速度) 激励频率:250/60=4.167 Hz;1 630/30=2.717 Hz 第4期 4.2水动力特性 叶永林等:基于水弹性力学的SWATH船结构…433 图4、图5分别是该试验船在水中的附加质量和阻尼。由图4可见,在波浪中的船体有一定的附加 质量,对船体结构的响应会产生影响,尤其是第l2阶振型模态,其附加质量比较大,达到船体质量的 l0之 10一-量级,因此总体振动计算时,考虑附连水的影响是必要的。同时,可以看到,仅在低频段(1 Hz 以内)波浪频率对船体的附加质量影响较大,频率较大(大于1 Hz)时,附加质量趋于稳定。由图5可 知,阻尼仅在1 Hz以内有一些波动,频率较大(大于1 Hz)时,将逐步趋于零。 lO’’ 1 0・t 10 10・_ 01 0‘ 菌,0 1O 7 10’8 1 O 10.,口 ,(Ik) ,( 图4附加质量 Fig.4 Additional mass 图5阻尼 Fig.5 Additional damp 4.3湿谐振频率 船体响应的主坐标频率,即湿谐振频率见图6。 10 1 E 06 1E.O6 1 O。7 BE.O7 亡 6E.0 7 1 0 龟 tO一9 4E.O7 2E.0 7 l0’1- 0 … ,(H 图6主坐标响应 Fig.6 Principal coordinate response 根据机械设备激励频率与上述湿船体固有频率,可以发现,设备激励与船体固有频率之问有一定 的间距,不会激起船体结构总体振动。 5船体结构强迫振动分析 船舶航行时,主要受到波浪激励、机械激励及螺旋桨的激励。SWATH船耐波性好,首阶湿谐振频 率较大,且船体尾部不受桨的表面力作用。本文仅考查二甲板作业区的振动情况,此处主要考虑受柴 油发电机的激励下的垂向结构振动。 5.1局部振动分析区域 结构的局部振动分析区域为艉二甲板的工作区与试验室,布置图见图7。 根据前述位置提取局部模态,在低频率部分,多数与总体模态相同,在60 Hz以后,节点垂向位移 较大,提取出来的模态数量相对较多一些,可能与结构模型的精细度有关。 第4期 输入机械阻抗。 叶永林等:基于水弹性力学的SWATH船结构… 435 主发、辅发的等效激励力为: Ⅳ2(厂). 测/V 同样,在主机、辅机基座上分别施加垂向单位正弦激励力,并求解(2)式中的广义水弹性力学运动 方程,可分别得到基座在单位垂向力作用下的垂向位移加 、 。 5-3航行工况下的局部振动 当船舶航行时,求解非稳态流场速度势的移动脉动源Green函数表示形式为[61: G( ;£ ) 赢+ 三坤书2 畚J… ( ,z,k)k-sgn(c0s , 式中: 0( ,Y,z,k)=exp【kZ+jk(XcosO+YsinO)] sgn(cosO): 一]、/1+d4O (5) 1, 一号≤Re( )≤号 1,一1T≤Re( )≤一号 ) X=K( ),Y=K(y一叩),Z=K( K=g/U,r=Uo@  ̄p=aretan(Y/X) a=arecos(1/4 ̄")  ̄=sinh I/ +y ) {外 ( 佩)] 航行时船体二甲板所受激励力,可根据该节点的实测加速度与设备基座加速度的对比来分析,如 图9(b)。由该图可见,二甲板的局部振动频率响应特征,主要受到发电机组的激励影响。因此,对二甲 板上的局部振动问题进行分析时,采用发电机组作为激励源是合理的。 为计算船体任意节点在机械激励下的振动响应,可先计算单位激励下的节点位移、加速度,然后 倍乘上述等效激励力。采用在关心节点上施加如(4)式所述的单位正弦激励力,求解(2)式,可得到节 '0O 峰值频率(Hz) 23.9 37.5 50 误差(dB) —12.6 12.9 一 。0 81) l3.6 9.2 14.2 皂7o 々 62.5 75 87.4 60 50 7.8 -99.9 9.5 ●O ,I (a)各峰值频率下的峰值加速度 a)Peak value of the acceleration (b)节点振动加速度 (b)Acceleration of the focus node 图l0航行(12 kns)时二甲板艉部节点振动计算与实测对比 Fig.10 Comparison of calculation results and test results of the node vibration response at main desk produced by the main generator at 12 kns 436 船舶力学 第17卷第4期 峰值频率(Hz) 23.9 37.5 49.8 误差(dB) —15I2 -6.6 21.2 0.8 18.1 62-2 74.4 81.8 87.1 -5-3 6.0 (a)各峰值频率下的峰值加速度 (b)节点振动加速度 (a)Peak value of the acceleration (b)Acceleration of the focus node 图11航行(12 kns)时潜体外侧节点振动计算与实测对比 Fig.1 1 Comparison of calculation results and test results of the node vibration response on the sub—hull produced by the main generator at 12 kns 点在机械激励力下的强迫振动的主坐标响应等,再根据实船在试验工况下所测试出的发电机基座加 速度(见图9),按上节所述方法所计算得出的激励力作为发电机的激励力,可求得局部振动结果,见图 10。 同时,本文对下潜体结构振动作了分析。图1l为该船航行时下潜体振动预报结果。分析该节点振 动时,考虑到推进电机工作时会对结构产生激励,采用了柴油发电机组与推进电机作为激励源,以其 航行时的实测加速度换算为激励力。 总的来说,在航行中二甲板与潜体结构的局部振动加速度计算结果与试验结果吻合,响应频率均 能一一对应,振动响应加速度幅值误差也较小。振动加速度满足中国船级社的相关要求…】。 6结构噪声分析 6.1声介质中水弹性力学频域分析方法 假定浮体周围为可压理想声介质。辐射速度势满足带航速的自由液面声学边界条件为: 2 ( 一 吾)咖=0 为: (6) 其中: 为作用在浮体上外载荷的激励角频率。引入与该自由液面条件相适应的频域脉动Green函数 G(P,Q; =击e h1一 1 e 其中:k=to/c。为流体中声波波数。 (7) 辐射势可用平均湿表面 及水线C上的源汇分布边界积分方程表示: ( 』so-(Q)G( 一 )dsQ+吉 f 。(Q) (Q)G(弛一 )dz 声介质中广义水弹性力学运动方程可表示为如下形式: (8) (9) a ]{q}+Ib+8 }+【c+c]{q}={ } 流场中的辐射声压可表示为: p( )= 。[暑+ ・ ] 如果浮体是扁的或狭的或细长的结构,略去浮体对均匀稳态流场的扰动,此时有: 0o) 第4期 叶永林等:基于水弹性力学的SWATH船结构… 437 =一Ui (11) 可得: ^ = (ttrn1+ n2+ ,n3)+ [n36Ir2一n2 3 J (12) , 本文分析零航速及6 kns航速航行时机械激励引起的结构噪声,包括近场噪声与辐射噪声。航速 大于10 kns时,流体噪声的影响较大,本文不作分析。 6.2航行时的近场噪声与辐射噪声分析 利用前述三维船体有限元模型,并对设备基座有限元模型进行了细化,计算了基座的机械阻抗。 SWATH船航行时,产生结构噪声的机械设备主要为发电机和推进电机。因此,本文考虑了主发电 机、辅助发电机及推进电机的影响,激励力计算如5.2节所示。利用对应的基座振动加速度测试结果和 上述基座的机械阻抗计算结果,可以计算出相应的设备激励力。 利用6.1节所述的水弹性分析方法,可以计算单位激励力作用下的流场中的辐射声压,考虑上述 激励力,从而可以得到该船行时的辐射噪声和近场噪声。 在6 kns航速的工况下,对水深7.5 m,横向水平距离65 in的水下位置进行了辐射噪声的预报。 图l2为三种激励辐射声总能量三分之一倍频程谱级与测试值比较曲线。 图13为近场噪声计算结果。计算场点为水深5.29 m,重心后11 nl处的纵舯面处。 ▲ l’ I一一一一一Te 【】土‘ I l 0 a1el at4nn r lJ1t‘ l 肇 . ~击 、i ,, , 厂 \ i聿 V1, VlI 、- ‘、\ ¨} j、 , { , , l , 蹲!r I I} , , \ \ /- } , , \ ,_/ , f, , \ / l \ ~《 ,t, , 7 7 T I、 、、 A 、 j1 i 、 / — 一1 1 ≠一 X ,… ,J l ,、 _,l、 ~一一rf t_{ ¨ \ '£ 一— 凸I l lf j ,(Hz】 /(Hz) 图14辐射噪声1/3OCT谱级 图15近场噪声I/3OCT谱级 Fig.1 4 I/3OCT of the noise radiation Fig.15 1/3OCT of noise of the Real"field 438 船舶力学 第17卷第4期 由图15可知,静止时该船的近场噪声水平较小。由于零航速时,机械激励力较小,因此近场噪声 与辐射噪声较小是合理的。 7结 论 通过采用三维水弹性方法,利用三维有限元模型和模态叠加原理,将波浪与机械设备的激励力有 效地结合,形成了航行中船体结构的强迫振动分析方法。该方法能有效地预报船舶在航行时的船体总 振动谐振频率,为控制船舶设备引起结构谐振提供参考。同时,可以预报在发电机等机械设备引起的 强迫振动下的船体局部振动情况,为船舶舱室及工作区的振动水平控制提供有力的参考。同样,利用 三维水弹性方法,可有效地预报船体结构噪声。但该方法对有限元模型建模以及模态的选取等提出了 比较严格的要求。 参考文献: [1】朱石坚,何琳.船舶机械振动控制[M].北京:国防工业出版社,2006. [2]Wu Y S.Hydroelasticity of lfoating bodies[D].Ph.D.Thesis,London,Brunel University,1984. [3]杜双兴.完善的三维航行船体线性水弹性力学频域分析方法【D】.无锡:中国船舶科学研究中心,1996. [4]崔维成,杨建民,吴有生,刘应中.水弹性理论及其在超大型浮式结构物上的应用【M】.上海:上海交通大学出版社, 2007. 【5]Ye Yonglin,Wu Yousheng,et a1.Hydroelastic analysis of SWATH ship wave loads ̄].Journal of Ship Mechanics.2010。 14(12):1340-1348. 【6】Zou Mingsong,Wu Yousheng,Ye Yonglin.Three-dimensional hydroelascity analysis of acoustic responses of ship stmc— tures[C]g 9th International Conference on Hydrodynamics,October l1—15,2010.Shanghai,China,2010:844—851. [7】叶永林,吴有生,等.水弹性力学在SWATH船体结构振动响应分析中的应用研究[J】.船舶力学,2011,15(5):506— 512. Ye Yonglin,Wu Yousheng,et a1.Research on application of hydroelasticity method to vibration response analysis of the SWATH ship hull[J].Journal f oShip Mechanics,2011,15(5):506—512. 【8】金建海,叶永林,田 超,等.三维水弹性集成软件的设计与实现【J].船舶力学,2011,15(5):521—529. Jin Jianhai,Ye Yongling,Tian Chao,et 1.Desaign and implementation for 3D hydroelastieity integrated software[J].Jour— nal of Ship Mechanics,201 1,15(5):521—529. 【9】席亦农,叶永林,等.SWATH船结构振动检测报告【R].无锡:中国船舶科学研究中心检测报告,2011. [10】熊草根,等.SWATH船航行结构噪声检测报告[R].大连:中船重工第760研究所检测报告,2008. [11】中国船级社.船上振动控制指南2000[M].北京:人民交通出版社,2000.