本文以ProNas能量有限元理論為基礎，建立了船舶的ProNas能量有限元計算模型，採用大型商用正版軟體ProNas對複雜激勵在船舶各艙室產生的中高頻結構噪聲及空氣噪聲進行模擬計算，得到船舶各艙室聲壓級，並利用ProNas軟體後處理功能確定激勵源及傳遞路徑處的能量分佈雲圖，對不滿足噪聲目標的艙室進行了聲學最佳化，最終解決了大型實際船舶工程的中高頻噪聲預測與控制問題。
1、引言複雜結構的中高頻噪聲的控制一直以來都是各工業領域研究的重點與難點問題，尤其對於大型船舶其內環境相比其它工業產品更加獨特[1]：結構形式縱橫交錯，艙室眾多，噪聲誘因複雜，聲源品種繁多密集，噪聲強度較大，噪聲頻域頻寬且持續穩定，結構噪聲與空氣噪聲相互轉化。以上這些特點，就使得船舶噪聲控制起來更加困難。2014年7月國際海事組織（IMO）簽訂生效的新的《船上噪聲等級規則》要求居住區部分艙室聲壓級在舊規範的基礎上降低5dB(A),這就要求船舶工程設計人員需要採取更加有效的控制手段來降低船舶噪聲。
目前預測結構噪聲及空氣噪聲的理論體系相對完善，並已將理論應用於大量實際工程中。按激勵源頻率及具體工業產品，可將振動和噪聲的問題劃分為低頻、中頻及高頻問題。低頻結構的響應具有隨機和確定性，且為全域性模態，工程中常用的數值方法為有限元（FEM）、邊界元（BEM）等。
理論上，有限元方法可計算任意結構的振動場，但是為了能準確反映結構的振動特性，劃分單元的網格長度要求要遠低於結構中彎曲波的波長。這樣隨著頻率的增加，結構彎曲波波長變小，結構網格需要進一步細化，對計算資源要求較高，且結構模態密集，重疊交錯，響應表現出不確定性，具有了統計的概念，因此，有限元和邊界元就不再適合解決中高頻問題。
統計能量分析（SEA）用於解決中高頻問題，且模態越密集，其計算精度就越高。但統計能量分析不能保證子系統的空間變數資訊的完整性[2]，難以精確預示子系統內能量密度分佈[3]且子系統的劃分需要一定的經驗，不易進行實際結構形態的設計與最佳化[4]，模態密度及耦合因子的準確與否直接影響結果的準確度[5-6]。所以為了更好的控制中高頻噪聲，就需要有更好的理論方法。
能量有限元法[7](EFEA)是一種預測中高頻動響應的新方法，它是以波動理論為基礎，將結構離散化，在單元之間建立能量密度關係式，從而求解得到所有節點的能量密度。能量有限元法以有限單元為物件，保證了模型的完整性，但是在結構突變處引入了大量重複節點，求解效率得不到平衡，能量有限元理論目前可以解決一些簡單的有限元噪聲問題，對於複雜問題的應用依然比較困難。
ProNas能量有限元是在統計能量分析及能量有限元理論的基礎上，以有限單元為研究物件，混合了SEA及EFEA理論，利用有限體積法及差分法推匯出得類似於SEA的理論方程，聯立求得每個有限單元的能量密度。ProNas能量有限元法避免了在結構突變處引入重複節點，保證了求解精度和計算效率，是一種具有較高研究價值及應用價值的中高頻噪聲控制方法。
本文著重介紹了ProNas能量有限元理論原理且應用ProNas商業軟體對大型實際船舶的中高頻噪聲進行了模擬計算，求解得出船舶各艙室聲壓級。
2、ProNas能量有限元基本理論ProNas能量有限元法可有效解決結構振動系統中高頻噪聲問題，它包含了能量有限元法（EFEA）及統計能量分析法（SEA）的優點，具有更大優勢。
ProNas能量有限元方法是在波動理論上建立的一種功率流方法，透過波動形式求解結構微元體運動方程，並將波動形式下的動能密度、勢能密度及功率流進行週期內時間平均和區域性空間平均，得到能量密度與功率流的關係，代入穩態下能量平衡方程中可以得到能量密度控制方程；採用一定的數值離散方法對控制方程進行離散，得到能量有限元方程；在結構存在耦合的地方考慮波的反射和折射，求出耦合處能量密度與功率流的關係；最後進行總係數矩陣的組集，得到總體的方程，從而求出能量密度。
有限體積法的基礎是將模型離散成若干個控制體，而這裡的控制體可理解為有限元的“單元”，這樣，ProNas能量有限元就是以單元為研究物件。
3、船舶模型建立① 船舶前處理模型
本文以某客箱船為例。利用Patran前處理軟體對某客箱船進行有限元建模，如圖2所示，該模型共有單元114358個，節點94357個。該客箱船全長約200米，船寬約25米，採用的是雙機雙漿配置，主機為兩臺低速機。該船的前部為載貨區域，後部為居住區域和機械處所。因對居住區域及機械處所噪聲控制要求較高，同時為了節省計算時間，本案例擷取區域性模型，即僅對居住區域和機械處所進行研究。
圖2.某客箱船有限元模型
② ProNas軟體模型處理
本案例透過ProNas聲學模擬軟體分別對船舶的結構噪聲與空氣噪聲進行了模擬計算。將Nastran格式模型匯入ProNas軟體中，如圖3所示。透過ProNas軟體聲腔自動識別功能自動劃分聲腔子系統，如圖4所示。軟體自動耦合聲腔子系統和其臨近貼合的結構單元，這樣能量即可在結構單元與聲腔子系統之間進行交換。然後對各艙室結構單元進行聲學材料設定，如甲板敷料、絕緣佈置材料，設定門窗材料等，如圖5所示。
圖3.某客箱船ProNas有限元模型
圖4.某客箱船各艙室聲腔子系統
圖5.某客箱船聲學材料施加
4、船舶模型激勵輸入船舶噪聲[9-10]主要包括主機噪聲、螺旋槳噪聲及水動力噪聲。其中主機噪聲是船舶噪聲中最強的噪聲源，主機噪聲主要分為空氣輻射噪聲、排氣噪聲和結構噪聲，空氣輻射噪聲有燃燒噪聲和機械往復運動碰撞噪聲，成因複雜，一般實測得到。螺旋槳噪聲也是船舶的主要噪聲源，螺旋槳直接產生的噪聲有空泡噪聲，空泡噪聲會引起船體結構噪聲。水動力噪聲主要是由於高速海流的不規則起伏作用於船體溼表面，激起船體的區域性振動並向周圍介質(空氣與水)輻射的噪聲。
該客箱船模型根據實際情況所受激勵載荷主要有24個，分別為主機、螺旋槳、發電機、空壓機、空調、風機等。其中主機、發電機、螺旋槳與船體作用，船體會振動產生結構噪聲。其它主要激勵源透過空氣傳遞到船舶艙室內，產生空氣噪聲。結構噪聲源透過ProNas軟體分別在各裝置的有限單元表面以速度級及加速度級方式載入，如圖6所示。空氣噪聲源透過ProNas軟體分別在各激勵源所在聲腔子系統內以聲功率級方式載入，如圖7所示。主機速度級激勵如圖8所示，螺旋槳加速度級激勵如圖9所示。
圖6.某客箱船結構噪聲載荷分佈圖
圖7.某客箱船空氣噪聲載荷分佈圖
圖8. 主機激勵頻譜圖
圖9.螺旋槳激勵頻譜
5、船舶結構噪聲及空氣噪聲分析結果與最佳化① 船舶噪聲分析結果
在倍頻程中心頻率63-8000Hz下分別將該客箱船結構噪聲模型與空氣噪聲模型提交ProNas軟體求解器，經過ProNas軟體計算，可得到各艙室聲壓級水平，圖10、圖11分別為該客箱船在中心頻率為500Hz時結構噪聲模擬得到的結構速度雲圖與艙室聲壓雲圖，圖12與圖13分別為該客箱船在中心頻率為500Hz時的空氣噪聲模擬得到的結構速度雲圖和艙室聲壓雲圖。
圖10.結構速度雲圖（500Hz-結構噪聲）
圖11.艙室聲壓雲圖（500Hz-結構噪聲）
圖12.結構速度雲圖（500Hz-空氣噪聲）
圖13. 艙室聲壓雲圖（500Hz-空氣噪聲）
該客箱船艙室聲壓分析結果見表1，部分艙室聲壓級不滿足目標值。
② 最佳化方案
對於結構噪聲超標的艙室，常見的最佳化方法為敷設阻尼。阻尼材料是將結構振動板的振動能量快速轉化為熱能，從而減弱金屬板的彎曲振動，阻尼材料透過這種方式有效的控制金屬板的結構輻射噪聲。
透過雲圖分析可知，本案例結構噪聲超標，多為主機與螺旋槳結構噪聲引起。由於該客箱船採用低速機，主機與船體鋼板直接螺接，船體底部振動區域較大，且底部鋼板厚度約為25mm-30mm，在實際應用中，阻尼層厚度一般為金屬板厚度2-4倍，因此如果採用常規的在激勵源處阻尼敷設方法，阻尼用量相對較大，成本較高。
又由於船舶結構複雜，結構噪聲和空氣噪聲相互轉化，因此常規阻尼敷設方法對超標艙室噪聲控制效果一般。針對以上特點，本案例透過ProNas軟體後處理介面可直觀精確顯示傳遞路徑處能量分佈，即在不達標艙室處直接敷設阻尼，採用這種敷設方法，阻尼用量及成本在可控範圍內，且噪聲控制效果顯著。圖14為常規敷設阻尼方法，圖15為本案例阻尼敷設方法。
圖14.常規敷設阻尼區域
圖15.本案例敷設阻尼區域
對於空氣噪聲不達標艙室，結合ProNas軟體後處理結果雲圖，本案例採用以下最佳化方案：
l 部分娛樂室及放映室等地面浮動地板均改為50mm岩棉+鋼板+阻尼+鋼板型，再在上側加10mm流平甲板敷料。
l 部分甲板廚房抽風機室的內部，艉、中、艏面增加4mm鋼板做成雙牆艙壁，形成密封夾層空間，兩層鋼壁都設阻尼塗料。
l 甲板房間合計26間處所的天花板需穿孔（穿孔面向室內），孔徑按1.8mm。壁板不穿孔。
l 在部分甲板機艙棚兩側直接面向公共處所和房間的艙壁，在外側增加4mm壓筋板做成的艙壁，形成密封的夾層空間，夾層內設50mm隔音棉。
l 第五、部分甲板均設9mm厚阻尼塗料處理。
經過以上方案最佳化並多輪模擬計算，最佳化以後各艙室噪聲控制效果較好，將以上最佳化方案應用到實船設計中，並對該客箱船實船艙室進行了噪聲測試，如表2所示，大多艙室噪聲數值滿足設計要求。
6、結論本文簡要介紹了ProNas能量有限元方法，在此基礎上應用ProNas軟體建立了某客箱船聲學模型，進行了整船的結構噪聲與空氣噪聲模擬計算，對結構噪聲及空氣噪聲不達標艙室進行模擬最佳化，並與實船測試結果進行對比，得出了以下結論：
l ProNas能量有限元方法是一種用於預測和解決船舶及其它工業產品中高頻噪聲的可行性及有效性方法，具有很大的理論意義與現實價值。
l 採用ProNas能量有限元方法，可避免常規中高頻方法子系統建模操作的複雜性與專業性，提高模擬計算精度。
l ProNas能量有限元方法不需要判斷子系統模態密度，不需要確定系統之間的耦合因子，既保證模型的完整性，又提高模擬的時效性。
l ProNas軟體為解決中高頻噪聲提供了新的方法，為改進船舶設計引數提供了科學依據，提高了船舶的噪聲設計質量，縮短研發週期，節約開發成本。
l ProNas軟體模擬得到的船舶結果與實測結果的最大誤差值的絕對值在5dB(A)內，滿足模型的預測精度。
參考文獻
[1] 陳小劍. 船舶噪聲控制技術[M]. 上海：上海交通大學出版社，2012:12-107.
[2] 林志立. 能量有限元法的模型及其在車身中的應用[D]. 中國科學院大學博士學位論文，2018，寧波.
[3] 原凱，王建民，韓麗，等. 能量有限元在振動與噪聲預示中的研究進展[J]. 強度與環境 2015. 42(3)：12-19.
[4] 祝丹暉，解妙霞，孔祥傑，等.複雜機械結構中高頻動響應能量有限元方法研究[J]. 中國工程科學，2013，15(1): 107-112.
[5] 龐劍，諶剛，何華. 汽車噪聲與振動----理論與應用[M]. 北京：北京理工大學出版社，2005: 54-62.
[6] 陳書明，王登峰，馬一功，等. 模態密度計算精度對車內噪聲預測精度的影響[J]. 吉林大學學報，2010，40(2):335-340.
[7] 孫麗萍，聶武. 能量有限元法在船舶結構中的應用[J]. 哈爾濱工業大學學報，2008，40(9):1491-1494.
[8] BERNHARD R J. The family of EFA equations and their relationship to SEA[C]//Proceedings International Conference Noise and Vibration Energy Methods．Lyon：[s.n.] 2000：248-257．
[9] 翁長儉. 船體振動學[M]. 大連：大連海運學院出版社，1992
[10] 邱斌. 高速船全頻段艙室噪聲預報與控制方法的研究[D]. 武漢理工大學碩士學位論文，2010，武漢.