強(qiáng)化水下航行器聲隱身性能一直是提高海軍綜合突防能力、生存能力和作戰(zhàn)效能，并取得戰(zhàn)略、戰(zhàn)術(shù)優(yōu)勢的重要途徑和核心舉措。由于重流體作用，水下結(jié)構(gòu)的聲源分布和聲輻射機(jī)理特別復(fù)雜，這給水下聲輻射預(yù)報和低噪聲設(shè)計帶來不少的難題。本文首先對水下聲輻射機(jī)理進(jìn)行了梳理；然后簡要介紹了Simcenter Acoustics聲仿真工具；最后，分別針對不同水下聲源給出了聲輻射仿真方法和流程，同時也分享了一些仿真案例，為相關(guān)的水下聲研究提供仿真經(jīng)驗和數(shù)據(jù)參考。
水下航行器噪聲的主要聲源有：機(jī)械噪聲、螺旋槳噪聲和水動力噪聲。這三種聲源根據(jù)產(chǎn)生的部位和機(jī)理，相互之間相互獨立也相互有所關(guān)聯(lián)。在低速隱蔽航行工況，機(jī)械噪聲是其最主要的噪聲源，其譜結(jié)構(gòu)特征也最容易被敵方探測到。機(jī)械噪聲的主要來源是：動力設(shè)備與管路系統(tǒng)和艉部推進(jìn)傳動系統(tǒng)；在巡航和高速工況，螺旋槳噪聲和水動力噪聲的貢獻(xiàn)量逐漸增大以至于占主導(dǎo)地位。有時為了簡化，水聲研究人員也會將這兩種噪聲統(tǒng)稱為流制噪聲，將螺旋槳作為最重要的流噪聲聲源。本文研究內(nèi)容是對水下聲輻射機(jī)理進(jìn)行詳細(xì)論述，并針對這5種機(jī)理采用不同的方法和流程來進(jìn)行聲學(xué)仿真分析。
二、水下聲輻射機(jī)理
2.1 結(jié)構(gòu)振動輻射聲
結(jié)構(gòu)振動輻射聲的聲源特征可以視為結(jié)構(gòu)輻射面上一個個具有一定相位關(guān)系的活塞輻射，結(jié)構(gòu)表面振動引起附近流體的壓縮和擴(kuò)張，密度變化而形成聲波傳播出去。因此，在考慮煤質(zhì)振動速度時需要考慮煤質(zhì)對結(jié)構(gòu)振動的影響，在水聲學(xué)中稱為附連水的影響。在進(jìn)行水下結(jié)構(gòu)振動聲輻射時，通常將結(jié)構(gòu)振動與聲輻射分開處理（這與流固耦合中聲振區(qū)分是兩個不同的概念）。在結(jié)構(gòu)動力學(xué)模型中，基于結(jié)構(gòu)的固有屬性（干模態(tài)）建立附連水的模型從而獲得濕模態(tài)，并在此基礎(chǔ)上加載載荷、流體力或傳遞的振動來計算獲得結(jié)構(gòu)表面的振動位移、速度或加速度。這些參量就是聲輻射仿真的聲源量，它已經(jīng)考慮了聲煤質(zhì)對結(jié)構(gòu)的影響，于是聲仿真就可以采用聲學(xué)波動仿真直接求解聲場。2.2流體繞流噪聲
流體動力噪聲屬于流體力學(xué)和聲學(xué)的交叉學(xué)科領(lǐng)域。當(dāng)飛機(jī)、艦艇、魚雷等運(yùn)載武器在流體介質(zhì)中運(yùn)動時，物體邊界層由層流發(fā)展為湍流，即在時間上和空間上都會有隨機(jī)變化的不穩(wěn)定流動狀態(tài)。這種非定常的流動一方面直接產(chǎn)生輻射噪聲，另一方面，湍流邊界層內(nèi)隨機(jī)的速度擾動產(chǎn)生隨機(jī)的脈動壓力，激勵物面彈性結(jié)構(gòu)振動并產(chǎn)生二次輻射噪聲。上述兩部分不同來源的噪聲我們統(tǒng)稱為流致噪聲。這里的艦船流體繞流噪聲分析指的是非定常湍流引起的直發(fā)聲。在進(jìn)行流體力引起的聲學(xué)仿真時，首先需要對聲源的產(chǎn)生機(jī)理和聲源特性要有一定的理解。通常繞流不會帶來流體質(zhì)量流量的變化（反例：水泵喘振），而且流動處于低馬赫數(shù)下，此時線速度較大的推進(jìn)器葉片表面的非定常力為主要聲源，其次是具有結(jié)構(gòu)引起流場邊界層變化較大的區(qū)域：如圍殼、尾舵、定子等位置，而且這些位置通常受到非均勻來流的影響。也是因為馬赫數(shù)較低的緣故，流體繞流聲仿真中的流場非定常計算與聲輻射仿真可以分開，其中不考慮流場對聲場的散射效應(yīng)。在艦船流體繞流噪聲仿真中，最關(guān)鍵的是要獲得精確的聲源信息，了解聲源的成分、聲輻射效率和聲源分布。這部分工作是需要聲仿真工程師與CFD工程師聯(lián)合進(jìn)行完成的：進(jìn)行網(wǎng)格適應(yīng)性研究，兼顧宏觀量如阻力、力矩平衡和微觀非定常力加密區(qū)域的聲源提?。贿€需要權(quán)衡計算量與計算精度，提取主要噪聲源。2.3 聲散射噪聲
根據(jù)FW-H方程中的繞流噪聲特性分析中，湍流脈動力在界面上形成偶極子為主要輻射聲源，其中在許多脈動力噪聲源中，與螺旋槳有關(guān)的噪聲源通常占主導(dǎo)地位。于是，該方案在此給出一種一種快速預(yù)報推進(jìn)器噪聲的方法：首先考慮主要聲源的輻射場，通過CFD計算獲得螺旋槳中的表面壓力脈動和軸向脈動力（通常CFD獲得的軸向一階葉頻脈動力精度能達(dá)到30%）；其次，將軸向推力等效為偶極子源；最后，將螺旋槳誘發(fā)的艇體輻射場可以視為是由螺旋槳直接輻射，入射聲經(jīng)過艇體激發(fā)艇體的反射聲場共同疊加形成的聲場。其實在聲仿真中，都不可避免的遇到聲散射問題，如結(jié)構(gòu)表面、自由液面、海底淤泥等等，在這里單獨拿出來分析主要是為了更好地詮釋艦船局部結(jié)構(gòu)的聲場特性和散射體的聲特性（又稱為固壁因子）。當(dāng)散射體邊界條件與聲源位置一定，其固壁因子不會發(fā)生變化，該參數(shù)在聲源結(jié)構(gòu)布置安裝和匹配中具有一定的參考價值。2.4 流激勵結(jié)構(gòu)振動輻射聲
在研究水下結(jié)構(gòu)二次輻射聲之前，首先需要對邊界層紊流壓力脈動激勵源做一個簡單的區(qū)分。脈動邊界層紊流壓力脈動通常區(qū)分為聲壓和偽聲兩部分：聲壓起因于湍流中的起伏Reynolds應(yīng)力產(chǎn)生的密度起伏，它服從波動方程；偽聲起因于湍流速度起伏的動量起伏，它直接平衡與動量起伏，滿足泊松方程。對比邊界層中的聲壓成分可以發(fā)現(xiàn)：偽聲壓力脈動能量占絕大部分，可見它是二次輻射噪聲的激勵源。繞流引起的二次聲仿真方法主要針對激勵力和結(jié)構(gòu)動力系統(tǒng)兩個方面。在激勵力方面，由于結(jié)構(gòu)受到繞流激勵的強(qiáng)度不僅取決于繞流表面指定點的自功率譜，還取決于壓力脈動沿空間的隨機(jī)關(guān)聯(lián)性（相關(guān)性）的程度，在給定頻率處這一脈動相關(guān)性的程度恰好表征互譜。隨機(jī)壓力脈動在空間的關(guān)聯(lián)性越強(qiáng)，結(jié)構(gòu)的振動激勵也就越有效。采用互譜矩陣的方式來加載激勵力是隨機(jī)振動的最直接有效的方法。2.5其他全頻段噪聲機(jī)理
古金噪聲：螺旋槳的勻速旋轉(zhuǎn)運(yùn)動導(dǎo)致的聲輻射，以此來強(qiáng)調(diào)與不導(dǎo)致聲輻射的勻速平移運(yùn)動的區(qū)別。代表了旋轉(zhuǎn)機(jī)構(gòu)的最低聲輻射。在水下低馬赫數(shù)時，這類噪聲聲能量很低，且被不均勻和非定常來流引起的聲輻射所掩蓋。
螺旋槳推、扭力噪聲（旋轉(zhuǎn)聲）：在不均勻來流中，螺旋槳周期掃掠流場，使得周圍流體介質(zhì)承受一種周期起伏力而輻射噪聲，這些噪聲的頻譜都是線狀譜。通常呈現(xiàn)軸頻、葉頻及倍頻特征，這是低頻線譜噪聲產(chǎn)生機(jī)理之一。
厚度噪聲：在非定常來流中，由于紊亂脈動的隨機(jī)特性，使得槳盤面速度場的周期性和重復(fù)性遭到破壞——使得槳面的非定常力的嚴(yán)格周期性喪失，在周期性分量外出現(xiàn)了非周期性的分量，不斷改變頻率和幅值的諧調(diào)分量，而形成的隨機(jī)分量的連續(xù)能量譜，這是低頻寬帶噪聲分量產(chǎn)生機(jī)理之一。有時，在葉頻附近的線譜來源于圍殼尾流卡門渦街導(dǎo)致的頻率和幅值的諧調(diào)。
渦旋噪聲：渦旋噪聲是流體流經(jīng)螺旋槳葉片、舵、各種突起物的后緣所產(chǎn)生的渦旋噪聲。特別是當(dāng)渦脫落頻率與物體的某階固有頻率相同且兩者的空間波數(shù)吻合時激發(fā)物體的共振，即唱音。通常情況下，槳葉是三維翼型，不同半徑處厚度不一樣，渦脫落頻率也不一樣，而且由于葉片間的葉柵效應(yīng)會破壞脫出渦的周期性，因此，它更接近于寬帶噪聲。
邊緣噪聲：由于槳葉后緣的半平面屏蔽效應(yīng)造成輻射源從偶極子源降級至單級子源，而輻射增大。該噪聲有明顯的指向性特征，垂直于流向。該噪聲頻率也由渦脫落頻率確定，屬于高頻范圍。
唱音：槳后緣的渦脫落形成周期力，當(dāng)頻率與彈性振動某個固有頻率一致時，會發(fā)生諧頻放大，并使得渦脫落增強(qiáng)。當(dāng)速度變化不大時，還會發(fā)生鎖頻自激振動。甚至有時唱音還會不同頻段發(fā)生，例如某螺旋槳300-400Hz、550-600Hz和650-700Hz。弱唱音產(chǎn)生的機(jī)理是來源于后緣產(chǎn)生的渦旋，這些渦旋引起的結(jié)構(gòu)振動僅僅增強(qiáng)了渦旋強(qiáng)調(diào)和調(diào)節(jié)渦脫落頻率，也就是前面提到的渦旋噪聲和邊緣噪聲，其表征的是聲輻射能量不僅集中在葉片固有振動頻率上，而且還存在于附件區(qū)域，它隨著流速變化比較平緩。
運(yùn)動激勵噪聲：由潛艇的軸系振動和艉部振動傳播至葉片而引起的聲輻射。頻率特征主要體現(xiàn)低頻窄帶（如艇體模態(tài)），也有中高頻窄帶（如軸系縱振）。
空化噪聲：物體與水的相對運(yùn)動引起局部壓力下降導(dǎo)致空化。螺旋槳空炮又分成稍渦空化、葉表面空炮和轂渦空泡。潛艇到達(dá)臨界轉(zhuǎn)速，空化噪聲強(qiáng)度急劇上升，其強(qiáng)度正比于空泡體積，空泡尺度逐漸從微米級增加到毫米，強(qiáng)度增大的同時，能量峰值逐步向低頻移動。在實際工程中，空泡通常會收到非均勻來流的影響，因此出現(xiàn)軸頻調(diào)制。
三、Simcenter Acoustics 仿真工具介紹
·Simcenter 3D繼承了Virtual lab高級聲學(xué)模塊，幾乎涵蓋了所有的成熟聲學(xué)仿真方法，針對大尺度模型、全頻段、多種聲源類型問題均具有詳細(xì)解決方案。
·Simcenter 3D沿用NX高級幾何建模功能，基于NX Nastran強(qiáng)大的結(jié)構(gòu)強(qiáng)度、剛度、模態(tài)、振動分析功能，通過整合流體、聲學(xué)、電磁等工程仿真軟件，實現(xiàn)從仿真驗證設(shè)計到仿真驅(qū)動設(shè)計轉(zhuǎn)變。
·Simcenter 3D軟件提供工作平臺導(dǎo)航功能，通過相應(yīng)菜單可以在設(shè)計、建模、和仿真功能模塊之間進(jìn)行任意切換，還能調(diào)用諸多其它仿真軟件內(nèi)核進(jìn)行求解。
·Simcenter 3D平臺集成了結(jié)構(gòu)動力學(xué)模塊和聲學(xué)模塊，能直接進(jìn)行結(jié)構(gòu)振動聲輻射聯(lián)合仿真。Simcenter 3D平臺提供多種仿真接口，通過調(diào)用結(jié)構(gòu)、流體仿真結(jié)果，可實現(xiàn)繞流噪聲和流固耦合聲輻射的計算。
·Simcenter 3D平臺的開放體系架構(gòu)(支持多種編程語言C/ C++、VB、C#、Java和Python) 為客戶化定制提供便利。
·Simcenter 3D平臺繼承了Virtual lab中的實驗和仿真混合建模及仿真分析功能，使其同時具有實驗的可靠性和仿真的便捷性。
四、水下聲輻射仿真方案
4.1結(jié)構(gòu)振動輻射聲仿真
由于聲源是結(jié)構(gòu)輻射面的振動信息，因此艦船結(jié)構(gòu)振動聲輻射仿真比較合適的方法是聲學(xué)邊界元法。針對內(nèi)外聲場封閉或聯(lián)通的問題，分別提供直接聲學(xué)邊界元法和間接聲學(xué)邊界元法，針對實際尺度模型提供多級子快速邊界元方法。本文考慮的是一般的外場聲輻射問題，所采用的方法為直接聲學(xué)邊界元法。4.2流體繞流噪聲分析
流體噪聲是通過CFD與FW-H聲類比、聲學(xué)軟件相結(jié)合計算聲場，即采用CFD模擬聲源周圍的區(qū)域，聲類比等效聲源，應(yīng)用聲學(xué)軟件求解聲場，該方法求解精度良好。聲學(xué)軟件中的聲學(xué)分析方法主要有兩種：有限元方法和邊界元方法。聲學(xué)有限元要求把所計算的聲場離散成實體網(wǎng)格，因此聲學(xué)有限元的計算空間通常是有限的，例如封閉空間的內(nèi)聲場問題。邊界元與聲學(xué)有限元相比，它有許多的靈活性，邊界元既可以計算封閉空間中的聲場，也可以計算非封閉空間中的聲場，邊界元需要的是面網(wǎng)格（二維網(wǎng)格），而不是實體網(wǎng)格，通過在面網(wǎng)格上積分，得到場點的聲場分布。邊界元相比有限元而言，雖然降了一個維度，但由于邊界元矩陣是滿秩矩陣，其計算所用的時間相差不大。4.3 聲散射噪聲計算
艦船的聲散射噪聲仿真重點在于聲源獲取和散射體建模。流體水動力獲得聲源的方法在前面導(dǎo)管槳算例中已經(jīng)有介紹，在這里要特別提的是，如何將其等效成簡單的聲源，當(dāng)然如果不考慮計算能耗仍然可以采用艦船流體繞流噪聲仿真的方法。散射體模型除了結(jié)構(gòu)外表面輪廓外，還有就是需要設(shè)定聲邊界，即穩(wěn)態(tài)聲場的邊界條件三類邊界條件：Dirichlet邊界條件(給定聲壓)，Neumann邊界條件(給定，為法向單位矢量)或Robin邊界條件(給定聲學(xué)阻抗，其中，和為給定的參數(shù))，例如剛性邊界為Neumann邊界條件=0。4.4 流激勵結(jié)構(gòu)振動輻射聲仿真
從多物理場仿真的角度來說，艦船的流激勵結(jié)構(gòu)振動輻射聲仿真只是將振動輻射聲中的激勵力換成由流場CFD獲得的脈動力，而且該脈動力具有遷移性特征。本文中的脈動力通過時域激勵力互功率譜來表征該激勵力特性。結(jié)構(gòu)振動仿真在前面章節(jié)中已經(jīng)講過了，就是利用結(jié)構(gòu)有限元軟件進(jìn)行干模態(tài)計算，并導(dǎo)入聲學(xué)軟件Simcenter中，采用邊界元將結(jié)構(gòu)干模態(tài)與聲場進(jìn)行耦合計算獲得結(jié)構(gòu)的濕模態(tài)。聲輻射計算在輻射表面振動信息已知的情況下，就是通過聲學(xué)邊界元或有限元來進(jìn)行求解。4.5 其他全頻段噪聲仿真
經(jīng)過對推進(jìn)器噪聲的產(chǎn)生機(jī)理進(jìn)行梳理之后，結(jié)合工程實際我們不難發(fā)現(xiàn)，唱音和空化噪聲的聲仿真幾乎是無法精確實現(xiàn)的。而其它的噪聲機(jī)理都可以用前面章節(jié)中介紹的方法來進(jìn)行仿真。唱音的仿真難點在于很難定義入流邊界，而且與結(jié)構(gòu)的制造工藝有關(guān)（同一型號的槳，工況一致，其中就有一兩條槳發(fā)生唱音）。然而，通過對流場仿真和槳葉結(jié)構(gòu)仿真以及唱音的機(jī)理分析可以有效地預(yù)防唱音的發(fā)生。空化噪聲仿真難點在于：1、聲源為非穩(wěn)態(tài)聲源，且只具有統(tǒng)計規(guī)律；2、聲源頻率高達(dá)10kHz，聲源尺度為幾毫米，將給聲仿真計算量巨大；最重要一點，CFD計算無法較精確的定量計算出聲源。在這些認(rèn)知基礎(chǔ)上，本方案尋求一些定性的仿真方法，如CFD+虛擬面FW-H方法和CFD+經(jīng)驗公式法。·CFD+虛擬面FW-H方法
·CFD+經(jīng)驗公式法
根據(jù)單個空泡的噪聲特性研究，空化的輻射聲功率是每個氣泡輻射的平均能量于每秒氣泡崩潰數(shù)目的乘積。由于每次崩潰輻射的能量正比于崩潰壓力與最大氣泡體積的乘積。因此，輻射功率正比于單位時間所產(chǎn)生的全部空化體積，即測量獲得的空化頻譜圖是在峰值前約9dB倍頻增加，峰值后約6dB倍頻降低。而峰值頻率可利用最大氣泡的崩潰關(guān)系式獲得：利用CFD計算獲得空化體積變化率和最大空泡半徑，利用聲功率公式計算并定性分析。另外，針對螺旋槳空化噪聲表達(dá)式，楊勇和熊鷹在文獻(xiàn)顯示噪聲頻率在高頻時空化噪聲比非空化噪聲聲功率級增加約20dB，獲得螺旋槳聲壓級經(jīng)驗表達(dá)式?？梢詫⑦@些經(jīng)驗公式封裝到聲學(xué)仿真軟件中，作為一種快速空化噪聲的預(yù)估仿真方法。總結(jié)
本文對水下聲輻射機(jī)理進(jìn)行了詳細(xì)的論述，幾乎涵蓋所有的聲源和頻段，為水下聲學(xué)研究提供理論指導(dǎo)。分別針對多種不同聲源機(jī)理提供了聲仿真的方法和流程，并通過相關(guān)的案例分享，能有效提高讀者的仿真實踐能力。能科股份介紹能科科技股份有限公司，成立于2006年12月，2016年10月在上交所上市（股票代碼：603859.SH），致力成為業(yè)內(nèi)領(lǐng)先的智能制造與智能電氣先進(jìn)技術(shù)提供商。堅持智能制造、智能電氣雙輪驅(qū)動戰(zhàn)略，依托先進(jìn)的工業(yè)軟件和電力電子技術(shù)，為客戶定制專屬的，以工業(yè)互聯(lián)網(wǎng)為核心的，數(shù)字化、網(wǎng)絡(luò)化、智能化系統(tǒng)解決方案。
