水線以上機電設備對船舶水下輻射噪聲影響研究

王家文，陳新傳，劉云生，霍柏琦，葉林昌

(1. 海軍裝備研究院，北京 100161；2. 中國船舶重工集團公司 第七一一研究所，上海 200090)

水線以上機電設備對船舶水下輻射噪聲影響研究

王家文1，陳新傳1，劉云生1，霍柏琦1，葉林昌2

(1. 海軍裝備研究院，北京 100161；2. 中國船舶重工集團公司 第七一一研究所，上海 200090)

建立小水線面雙體船聲學有限元模型，以柴油發電機組振動噪聲為激勵，計算該船型水線以上振動噪聲源產生的水下輻射噪聲。探討船舶水線以上振動噪聲對其水下輻射噪聲的影響規律，為整船水下輻射噪聲評估提供參考。

水下輻射噪聲；有限元；小水線雙體船

0 引 言

某型科考船采用小水線面雙體船型，全電力推進，由于其任務使命的特點，要求最大限度地減小船體結構噪聲輻射，控制科考船自噪聲對測量區域被探測目標的影響。目前，對于單體船的振動和聲學預報已經相對比較成熟，船體結構受外載荷激勵而振動，在流體介質中輻射形成聲場，聲場反過來對結構施加反作用力，從而形成一個聲振耦合系統[1]；而對于小水線面雙體船，其大部分機電設備都布置在水上平臺，機械設備產生的振動噪聲除了通過振源-船體-水下噪聲這一條傳統路徑進行傳播外，由于其船型的特殊性，機械振動噪聲還通過振源-船體-空氣-水下噪聲路徑進行傳播，國內船舶水下輻射噪聲預報，大多都是基于常規船型水下輻射噪聲預報的思想，對振源-船體-水下噪聲這一傳統路徑進行研究，張維等[2]采用結構有限元法及聲學邊界元法，對小水線面雙體船不同工況下聲學特性進行數值分析；嚴斌，李廣等[3–4]對小水線面雙體船振動傳遞特性進行模型仿真及實船試驗研究，且對其水下輻射噪聲進行預報。

然而，目前對于小水線面雙體船機械振動噪聲通過振源-船體-空氣-水下噪聲這一傳播路徑的相關研究報道較少。雖然空氣-海水 2 種介質阻抗相差較大，水面以上產生的輻射噪聲大部分被海面反射，能夠透射進入海水中的能量很少，但是當空氣中的平面波入射到水中，雖然其聲強衰減，但是其聲壓并沒有衰減，反而增大到原來的 2 倍[5]。本文利用現代船舶水下輻射噪聲預報技術，對小水線面雙體船水線以上機電設備輻射噪聲形成的水下聲場進行初步研究，仿真得到的預報結果可為該類船舶水下輻射噪聲的預報評估及減振降噪提供一定的參考。

1 數學模型的建立

1.1 聲輻射模型

將船體看作是一個多自由度的振動系統，在其內部基座點 a 激勵下的外板點 b 振動響應為[6]：

式中：F 為設備在其基座處產生的激勵載荷；λi為復模態頻率；Ψi為復振型；αi為復模態質量與復模態剛度 Ki以及復模態阻尼 Ci有關；Hba（ω）為傳遞導納。

船體表面振動能量傳遞到其他周圍介質的能力是影響船體輻射噪聲的關鍵因素。估算船體表面輻射聲功率為[7]：

式中：WR為輻射聲功率；ρ 為流體密度；c 為流體介質中的聲速；v 為結構表面振動速度；S 為結構輻射的表面積；σrad為該輻射面上的聲輻射效率。

1.2 聲在空-水面透射模型

從空氣中聲源 O 出發的聲線，如圖 1 所示，在空-水界面 T 處折射入水，然后到達水下 S 處，入射角為θ，折射角 φ，另外 1 條聲線 OT'O 其入射角為 θ + dθ，折射角為 φ + dφ，假定空氣中距聲源單位距離處的聲強為 I0，聲線管橫截面積為 A0；S 處的聲強為 IS，聲線管橫截面積為 AS。根據射線聲學，在空氣中，聲線管內的聲功率保持為[8]：

在水中，聲線管內的聲功率保持為：

聲在空-水界面處發生反射和折射，一部分聲能反射回空氣中，另一部分聲能透射進入水中。進入水中的聲功率與入射波的聲功率之比近似等于平面波功率透射系數，即

式中 Tw 為聲功率透射系數。則由聲功率透射系數定義可以得到 S 處的聲強：

圖1 聲波透射入水示意圖Fig.1 Air-to-water sound transmission

將式（6）代入式（5），可得到：

其中 m = ρ1/ρ2為海水與空氣的密度比。P0為距聲源單位距離處的聲壓幅度。 Ps為在水下 S 處的聲壓幅度。則聲強 I0和 Is可以近似表示為:

將式（8）代入式（7）中，有：

2 模型計算

該型船在設計吃水時，使甲板距離海面 3.6 m，如圖 2 所示。根據圖紙進行必要的簡化[9]，建立結構有限元計算模型，船體的甲板、艙壁等板殼結構采用二維殼單元，加強筋采用一維梁單元，船體輪機設備用集中質量點模擬，油水用三維流體單元模擬，柴油發電機組下的隔振器采用一維彈簧單元模擬。借助 MSC. Patran & Nastran 軟件進行船體模型的前后處理和計算工作，建立船體的有限元模型包括船體各層甲板、艙壁、外板、各個肋位的加強筋、縱橫強橫梁及艙壁的扶墻材等，船體模型的質量由船舶有效質量和附加水質量組成，有限質量主要考慮船體動力設備、舾裝、外裝、電器設備以及油水重量，確保有限元船體的總重量、重心和實際相差在 0.5% 以內。附加水質量則可通過指定單元與水接觸，并由 Nastran 軟件自動識別。小水線面雙體船整船有限元模型如圖 3 所示。

圖2 水線以上輻射噪聲透射入水示意圖Fig.2 Air-to-water sound transmission

圖3 整船有限元模型Fig.3 Finite element model

該型船在工作工況下，航速較低，只需 1 組柴油發電機組運行工作，柴油發電機組位于船體中部的機艙，如圖 4 和圖 5 所示，單臺功率 1 380 kW，轉速1 500 r·min–1，重約 24 t，采用雙層隔振，柴油發電機組振動噪聲通過減振器-船體和空氣傳播等路徑，引發濕甲板振動，形成水上輻射噪聲源。在試驗室測得工作工況下的單臺柴油發電機組基座面板振動加速度級，如圖 6 所示。對獲取的載荷激勵進行處理和加載，將其作為輸入激勵，計算由單臺柴油發電機組激發的船體振動而產生的輻射噪聲。

通過在有限元軟件 MSC.Patran & Nastran 中建立該船型的結構有限元模型，按照該船的標準排水量進行整船模態和振動響應分析，由于模型較大，相應的自由度較多，隨著頻率的升高，模態越來越復雜，現給出幾階較為典型的振動模態，如圖 7 ～ 圖 10 所示。

圖4 船體模型半剖圖Fig.4 Finite element model

圖5 局部半剖放大圖-中部Fig.5 Finite element model

圖6 柴油發電機組基座面板振動加速度級Fig.6 Vibration acceleration level of diesel generating set

圖7 1 階扭轉Fig.7 Vibration mode

圖8 潛體 1 階反相位水平彎曲Fig.8 Vibration mode

圖9 整體 1 階水平彎曲Fig.9 Vibration mode

圖10 整體 1 階垂向彎曲Fig.10 Vibration mode

將柴油發電機組基座面板振動加速度級作為輸入激勵，求解出船體的振動響應，由于構建的有限元模型較大，為提高運算效率，本文采用邊界元軟件 Virtual.Lab 對輻射噪聲的中低頻部分進行運算分析，其基本思路如圖 11 所示。

統計能量法充分利用了在頻率較高的頻段上振動和聲模態密集的特性，采用統計量描述系統的振動規律，是解決復雜系統寬頻帶動力學問題的有效方法，對有規則結構的振動與聲學分析比較有效，計算量也小于有限元 + 邊界元方法，本文運用統計能量分析軟件 VA One 對高頻部分的輻射噪聲進行分析，在 VA One 中導入有限元模型，生成 SEA 分析模型，定義各

材料屬性及各種加筋板，并賦予給對應的子系統，定義設備的振動噪聲，約束到對應子系統。最后求解出濕甲板的輻射噪聲級，將該噪聲級看作是由濕甲板中心的一個點聲源產生，只考慮直達聲波，忽略側面聲波和海底、海面反射的影響，設海面處空氣密度約為1.225 kg/m3，聲速約為 340.3 m/s，海水密度為 1 000 kg/m3，聲速為 1 500 m/s，得到水線以上平臺在柴油發動機組單獨激勵下，通過空-水界面透射到水下，形成的水下輻射噪聲（1 m 處，1 μPa）如圖 12 所示。

圖11 中低頻水下噪聲預報思路Fig.11 Prediction of underwater noise in medium and low frequency

圖12 柴油發電機組激勵下透過空-水界面的水下輻射噪聲Fig.12 Underwater radiation noise of diesel generating set

通過圖 4 與圖 12 進行對比可以發現，圖 4 中振動噪聲低頻部分較高頻部分較弱，而圖 12 中計算的結果卻剛好相反，這說明噪聲透射到水下時高頻部分比低頻部分衰減較多，而且隨著頻率的增加，衰減量也在增加，可以看出機械噪聲在水-空界面透射入水時，噪聲高頻部分透射性較差，隨著頻率的增加，噪聲在水-空界面透射的衰減也在增大。

通過圖 13 可以看出，在濕甲板正下方時，聲壓最大，且聲壓隨深度的增大而減小，當偏離濕甲板正下方時，聲壓呈減小趨勢，其聲壓場有余弦方向性。特別是在較小水深時，隨著水平距離的增加聲壓迅速衰減，水平距離每增大1倍，聲壓近似下降約 11 dB；當不在濕甲板正下方時，聲壓隨深度的增加先增大，在某個深度達到最大后，隨深度的進一步增加而減小。所以在深海探測時，為減小該船振動輻射噪聲對探測效果的影響，并不是將探測器放置在較大深度上就好。

圖13 水下輻射噪聲分布Fig.13 Underwater radiation noise distribution

3 結 語

本文以某型小水線面雙體船為對象，通過 MSC. Patran & Nastran 軟件建立了該型船結構有限元模型，以水線以上柴油發電機組的振動噪聲為輸入激勵，運用聲學邊界元軟件 Virtual.Lab 和統計能量分析軟件 VA One 對小水線面雙體船水線以上平臺引起的水下輻射噪聲預報進行初步探究，仿真發現由小水線面雙體船水線以上機械設備產生的振動噪聲，通過水-空界面形成的水下輻射噪聲主要為低頻噪聲，高頻噪聲在水-空界面處衰減較大，且隨著水平距離的增加，輻射噪聲迅速衰減，為小水線面雙體船機電設備水下噪聲和減振降噪預報提供了參考。

[1]孫雪榮, 朱錫. 船舶水下結構噪聲的研究概況與趨勢[J]. 振動與沖擊, 2005, 24(1): 106–113.

[2]張維, 李明, 陳爐云, 等. 小水線面雙體船聲輻射特性研究[J].艦船科學技術, 2010, 32(11): 17–20.

[3]李廣, 嚴斌, 劉朋, 等. 小水線面雙體船振動傳遞特性研究[C]//中國造船工程學會船舶力學學術委員會水下噪聲學組成立三十周年第十五屆船舶水下噪聲學術討論會論文集. 鄭州:中國造船工程學會, 2015.

[4]嚴斌, 李廣. 小水線面雙體船水下輻射噪聲計算研究[C]//中國造船工程學會船舶力學學術委員會水下噪聲學組成立三十周年第十五屆船舶水下噪聲學術討論會論文集, 鄭州: 中國造船工程學會, 2015.

[5]馬遠良. 聲波從空氣進入水中的現象及研究[C]//水聲情報交流會論文集, 宜昌, 1970.

[6]傅志方, 華宏星. 模態分析理論與應用[M]. 上海: 上海交通大學出版社, 2000.

[7]孫進才, 王沖. 機械噪聲控制原理[M]. 西安: 西北工業大學出版社, 1993.

[8]HUDIMAC A A. Ray theory solution for the sound intensity in water due to a point source above it[J]. The Journal of the Acoustical Society of America, 1957, 29(8): 916–917.

[9]American Bureau of Shipping. Guidance notes on ship vibration[S]. New York: ABS, 2006.

Research on the influence of above the waterline equipment on underwaterradiated noise of SWAH

WANG Jia-wen1, CHEN Xin-chuan1, LIU Yun-sheng1, HUO Bai-qi1, YE Lin-chang2
(1. Naval Academy of Armament, Beijing 100161, China; 2. The 711 Research Institute of CSIC, Shanghai 200090, China)

This article established the acoustic finite element model of the SWAH, the vibration and noise of the diesel engine group as a source, calculated the underwater radiated noise of SWAH, made a preliminary research on the influence of above the waterline equipment on underwater radiated noise of SWAH, and provided a reference for underwater radiated noise of SWAH.

underwater radiated noise；finite element；SWAH

TB53

A

1672 – 7619(2016)11 – 0057 – 04

10.3404/j.issn.1672 – 7619.2016.11.011

2016 – 04 – 27；

2016 – 05 – 31

王家文(1986 – )，男，碩士研究生，研究方向為振動與噪聲控制。