湍流脈動(dòng)壓力下橢球聲吶腔水動(dòng)力自噪聲分析

彭華康 李天勻 朱 翔

華中科技大學(xué)船舶與海洋工程學(xué)院，湖北武漢430074

湍流脈動(dòng)壓力下橢球聲吶腔水動(dòng)力自噪聲分析

彭華康 李天勻 朱 翔

華中科技大學(xué)船舶與海洋工程學(xué)院，湖北武漢430074

運(yùn)用統(tǒng)計(jì)能量法建立水下橢球殼形狀的聲吶腔模型。利用計(jì)算流體動(dòng)力學(xué)軟件計(jì)算不同航速下聲吶罩表面湍流邊界層的分離點(diǎn)，考慮了湍流脈動(dòng)壓力空間分布的不均勻性。分析航速以及不同吸聲和阻尼處理方案對(duì)聲吶腔水動(dòng)力自噪聲的影響。理論評(píng)估吸聲處理的結(jié)果和數(shù)值計(jì)算的結(jié)果比較吻合。研究表明：航速越大，自噪聲越大；透聲窗面積所占總表面積的比例越小，吸聲和阻尼處理的降噪效果越明顯；吸聲和阻尼處理在中高頻段降噪效果更好。研究結(jié)果為水下聲吶腔水動(dòng)力自噪聲預(yù)報(bào)和控制提供了參考。

統(tǒng)計(jì)能量分析；聲吶腔；湍流脈動(dòng)壓力；水動(dòng)力自噪聲

1 引言

聲吶的探測(cè)性能與其自噪聲密切相關(guān)，降低自噪聲級(jí)，可以提高聲吶檢測(cè)信號(hào)的能力，增大聲吶的工作距離。聲吶罩是聲吶基陣外具有良好水動(dòng)力性能的外罩，其作用是避免水流對(duì)聲吶基陣的沖擊，減小湍流、抑制空化、降低偽聲對(duì)基陣的直接干擾。航行時(shí)，聲吶罩表面會(huì)產(chǎn)生湍流邊界層，一方面湍流脈動(dòng)壓力直接向聲吶罩內(nèi)輻射噪聲；另一方面湍流脈動(dòng)壓力激勵(lì)聲吶罩結(jié)構(gòu)振動(dòng)輻射噪聲。Lighthill［1］深入研究了自由湍流輻射噪聲，Curle［2］拓展了Lighthill的理論，研究了流體中有剛體結(jié)構(gòu)時(shí)湍流的輻射噪聲。研究表明，有結(jié)構(gòu)存在的流體直接輻射噪聲強(qiáng)度與馬赫數(shù)的6次方成正比。考慮到水下結(jié)構(gòu)物低馬赫數(shù)運(yùn)動(dòng)時(shí)，湍流邊界層脈動(dòng)壓力直接輻射噪聲可以忽略不計(jì)，因此，聲吶罩水動(dòng)力噪聲主要是湍流邊界層脈動(dòng)壓力激勵(lì)聲吶罩結(jié)構(gòu)振動(dòng)的輻射噪聲。

國內(nèi)外很多學(xué)者已經(jīng)對(duì)聲吶罩的水動(dòng)力噪聲開展了研究。Dowell［3］建立了彈性平板覆蓋在矩形腔上的聲吶罩模型，計(jì)算了彈性平板在湍流脈動(dòng)壓力激勵(lì)下矩形腔中的噪聲。Lauchle［4］采用邊界層過渡區(qū)湍流猝發(fā)的聲輻射理論，給出了回轉(zhuǎn)體艏部聲吶部位自噪聲的估算公式。俞孟薩［5］、彭臨慧等［6］運(yùn)用統(tǒng)計(jì)能量法分析了均勻湍流邊界層激勵(lì)下，聲吶自噪聲和水下結(jié)構(gòu)物的自噪聲，均沒有考慮湍流脈動(dòng)壓力的空間不均勻性。張娟等［7］利用統(tǒng)計(jì)能量法建立了圓柱殼形狀的聲吶罩模型，進(jìn)行了空間分布不均勻湍流邊界層對(duì)聲吶罩輸入功率的計(jì)算及統(tǒng)計(jì)能量參數(shù)的確定。本文在上述研究的基礎(chǔ)上，運(yùn)用振動(dòng)聲學(xué)軟件VAOne建立了形狀更接近工程實(shí)際的水下橢球殼聲吶罩模型。為了考慮湍流脈動(dòng)壓力空間分布的不均勻性，運(yùn)用CFD軟件計(jì)算不同航速下聲吶罩表面湍流邊界層的分離點(diǎn)。利用數(shù)值分析了航速對(duì)自噪聲的影響，同時(shí)分析了在相同和不同透聲窗面積所占總表面積的比例下，吸聲和阻尼處理對(duì)聲吶腔水動(dòng)力自噪聲的影響等。理論評(píng)估吸聲處理的結(jié)果和數(shù)值計(jì)算的結(jié)果比較吻合。研究結(jié)果為聲吶腔水動(dòng)力自噪聲的預(yù)報(bào)和控制提供了參考。

2 統(tǒng)計(jì)能量分析模型

2.1 統(tǒng)計(jì)能量分析基本理論

上世紀(jì)60年代，Lyon［8］提出了統(tǒng)計(jì)能量分析方法，將整個(gè)系統(tǒng)用多個(gè)物理子系統(tǒng)模擬，子系統(tǒng)之間可以相互耦合，其中的1個(gè)或多個(gè)子系統(tǒng)受到寬帶穩(wěn)態(tài)隨機(jī)力激勵(lì)。通過建立子系統(tǒng)間的功率流平衡方程，估算子系統(tǒng)的能量分布。對(duì)于保守弱耦合的系統(tǒng)，其穩(wěn)態(tài)功率流平衡方程為：

式中，Ei，Pi，ni，ηi分別為第 i個(gè)子系統(tǒng)的平均能量、平均輸入功率、模態(tài)密度以及能量損耗系數(shù)；ηij為第i個(gè)子系統(tǒng)和第j個(gè)子系統(tǒng)之間的能量傳遞損耗系數(shù)；ω為1/3倍頻程的中心頻率。

統(tǒng)計(jì)能量分析建立在以下基本假設(shè)之上：

1）子系統(tǒng)是弱耦合的；

2）激勵(lì)在統(tǒng)計(jì)上是獨(dú)立的；

3）在給定頻帶內(nèi)所有共振模態(tài)能量之間能量等分；

4）功率流與平衡耦合模態(tài)能量流之間的差成正比。

VAOne提供高頻區(qū)的振動(dòng)聲學(xué)分析，高頻區(qū)是指帶寬Δf內(nèi)振型數(shù)N大于或等于5。

2.2 聲吶罩幾何模型及流場(chǎng)計(jì)算

如圖1所示，橢球殼短半軸R1＝1 m，長半軸R2＝2 m，原點(diǎn)為O點(diǎn)，U0為自由來流速度，來流方向與長半軸方向平行。

運(yùn)用CFD軟件對(duì)橢球殼進(jìn)行流體計(jì)算。受到計(jì)算機(jī)硬件資源和計(jì)算時(shí)間的限制，本文采用RGN k－ε湍流模型，DES和LES湍流模型計(jì)算結(jié)果更好，但是對(duì)計(jì)算機(jī)要求更高。先進(jìn)行穩(wěn)態(tài)計(jì)算，待其充分收斂后，驗(yàn)證流場(chǎng)模擬的精度再進(jìn)行非穩(wěn)態(tài)計(jì)算，最后確定邊界層分離的位置。速度不同，分離點(diǎn)位置有所不同。如圖1所示，假設(shè)分離點(diǎn)與原點(diǎn)O的連線與長半軸的夾角為θ，則速度越小，θ越小；速度越大，θ越大，分離點(diǎn)越靠前，這是因?yàn)椋俣仍酱螅嗤恢玫睦字Z數(shù)越大，流動(dòng)越劇烈，橢球殼表面的逆壓梯度越靠前，邊界層分離越早。當(dāng)雷諾數(shù)達(dá)到（3.5～5）×105時(shí)，橢球殼表面會(huì)形成充分發(fā)展的湍流邊界層，橢球殼表面沿著來流方向勢(shì)流壓力分布從最大值很快下降到最小值，接著又迅速回升，如圖2所示圖，流體由順壓梯度向逆壓梯度變化，此時(shí)在切應(yīng)力消耗動(dòng)能和流體減速增壓的雙重作用下，邊界層迅速增大，由于流體的不可壓縮性和慣性，邊界層在橢球殼的某點(diǎn)脫離后，形成渦流區(qū)。如圖3所示，此時(shí)橢球殼表面邊界層分為2個(gè)區(qū)域：湍流邊界層和分離區(qū)，并形成不一致的湍流脈動(dòng)壓力分布。Han等［9］研究指出，邊界層的每一區(qū)域都對(duì)輸入功率有重大影響，必須考慮邊界層分布的空間不均勻性。

圖1 聲吶罩幾何模型Fig.1 Geometrical model of the sonar dome

圖2 橢球殼近壁區(qū)壓力分布Fig.2 The pressure distribution of near－wall zone

圖3 邊界層分離Fig.3 Turbulent boundary layer separation

2.3 聲吶罩統(tǒng)計(jì)能量分析計(jì)算模型

在給定的速度下，進(jìn)行流場(chǎng)計(jì)算后確定邊界層的分離點(diǎn)，然后在VAOne中建模。VAOne較之以前的舊版本AutoSEA2_2004有很大的改進(jìn)，即可以導(dǎo)入有限元幾何模型，建立更加精確的曲面模型。首先在ANSYS中建立橢球殼模型并劃分單元，再將有限元幾何文件導(dǎo)入VAOne中進(jìn)行一系列的操作。如圖4所示，在無限水域中，忽略聲吶基陣的影響，聲吶罩模型分為8個(gè)單曲殼子系統(tǒng)，模擬聲吶罩外殼，右邊的2個(gè)子系統(tǒng)在湍流邊界層分離區(qū)，左邊6個(gè)子系統(tǒng)在未分離區(qū)；內(nèi)部建立1個(gè)三維聲腔子系統(tǒng)，模擬罩內(nèi)聲吶腔；在每個(gè)殼子系統(tǒng)上定義湍流邊界層，作為輸入功率源；在聲吶罩外部上方和下方各建立1個(gè)半無限流場(chǎng)，模擬外部輻射聲場(chǎng)。然后建立點(diǎn)、線、面連接，半無限流場(chǎng)需要手動(dòng)連接到殼子系統(tǒng)。

在VAOne軟件中，邊界層分離區(qū)和未分離區(qū)中的子系統(tǒng)設(shè)置如圖5所示，X0為湍流邊界層起點(diǎn)至各子系統(tǒng)表面脈動(dòng)壓力中心的距離。將此模型作為保守弱耦合系統(tǒng)進(jìn)行分析，設(shè)定結(jié)構(gòu)損耗因子，其他各分析參數(shù)由軟件自動(dòng)計(jì)算。

圖4 橢球殼統(tǒng)計(jì)能量分析模型Fig.4 SEA model of the ellipsoid wall

圖5 邊界層分離區(qū)和未分離區(qū)流動(dòng)Fig.5 The attached and separated boundary layer flow

3 自噪聲計(jì)算分析

本文計(jì)算頻率范圍為500～10 000 Hz，采用1／3倍頻程中心頻率。子系統(tǒng)在每個(gè)中心頻率的頻率范圍內(nèi)振型數(shù)均大于5，滿足高頻區(qū)條件。聲壓參考值為10－6Pa。相關(guān)研究指出［10］，聲腔一定時(shí)，聲吶罩壁的材料對(duì)自噪聲的影響很大。首先，應(yīng)該選用阻尼因子較大的材料，玻璃鋼較不銹鋼的阻尼因子大；其次，應(yīng)該選用吻合頻率較高的材料，降低輻射效率。吻合頻率主要與材料的縱波波速和壁厚有關(guān)，壁厚一定時(shí)應(yīng)選擇縱波波數(shù)小的材料，玻璃鋼的縱波波數(shù)較不銹鋼小。壁厚對(duì)自噪聲也有影響，通常情況下，罩壁越厚，自噪聲越小，這是因?yàn)橥牧髅}動(dòng)壓力激勵(lì)罩壁振動(dòng)所產(chǎn)生的噪聲與面密度成反比，但是厚度增加會(huì)降低聲吶罩的透聲性，增大聲波束畸變。本文中聲吶罩壁透聲窗材料選用玻璃鋼，非透聲窗材料選用不銹鋼。參數(shù)設(shè)置如表1所示。

表1 罩壁材料參數(shù)Tab.1 Material coefficient of the sonar dome wall

聲吶在海水中工作，在其工作的頻段內(nèi)，假設(shè)透聲窗透聲系數(shù)為95%，透聲窗的位置在來流方向的下方。在非透聲窗界面上貼覆吸聲材料或者敷設(shè)阻尼材料等降噪措施。其他參數(shù)設(shè)置如表2。

表2 參數(shù)設(shè)置Tab.2 Coefficient settings

圖6所示為無聲學(xué)處理且透聲窗面積占總表面積比例為1／4時(shí)，不同航速下聲吶腔水動(dòng)力自噪聲隨頻率的變化曲線。從圖可見，隨著頻率的增大，聲吶腔內(nèi)的自噪聲整體上逐漸降低，只是在個(gè)別頻率點(diǎn)有所突變，這和模型的固有頻率點(diǎn)重合。當(dāng)航速增加時(shí)，聲吶腔內(nèi)的自噪聲也不斷變大，這是由于航速增大后，湍流脈動(dòng)的強(qiáng)度隨之增大，約與航速的5～6次方成正比，此時(shí)湍流邊界層激勵(lì)殼體振動(dòng)向腔內(nèi)輻射的噪聲也變大。當(dāng)航速增大1倍時(shí)，聲壓級(jí)大約增大20 dB。

圖6 不同航速下聲吶腔水動(dòng)力自噪聲Fig.6 Sonar cavity flow－induced self noise among different cruising speeds

圖7所示為透聲窗面積占總表面積的比例為1／4時(shí)，非透聲窗界面上貼覆不同的吸聲材料對(duì)聲吶腔水動(dòng)力自噪聲的影響。分析可知，貼覆聲學(xué)材料和未貼覆聲學(xué)材料的聲吶腔水動(dòng)力自噪聲隨頻率變化趨勢(shì)一致，均隨頻率增大而減小；隨著吸聲系數(shù)的增大，腔內(nèi)自噪聲聲壓級(jí)越來越小。此外還可知，吸聲處理在中高頻率段降噪效果更加明顯。

圖8所示為不同透聲窗面積占總表面積的比例下，不同吸聲處理對(duì)聲吶腔水動(dòng)力自噪聲的影響。橫軸為透聲窗所占的比例，縱軸為聲腔內(nèi)的總聲壓級(jí)。從圖可知，吸聲處理的降噪效果和透聲窗面積所占總表面積的比例有很大的關(guān)系。研究表明，只有當(dāng)透聲窗面積占總表面積的比例較小，非透聲窗壁面上的聲反射在罩內(nèi)形成較強(qiáng)的混響時(shí)，吸聲處理才會(huì)有較好的降噪效果。罩內(nèi)混響時(shí)間越長，吸聲處理的降噪效果就越好，而混響時(shí)間和透聲窗面積所占比例有很大的關(guān)系。根據(jù)水下封閉空間混響時(shí)間T60公式［11］：

圖7 相同透聲窗比例，不同吸聲處理（U0＝10 m／s）Fig.7 Same acoustic window proportion with different absorption coating

圖8 不同透聲窗比例，不同吸聲處理（U0＝10 m／s）Fig.8 Different acoustic window proportion with different absorption coating

式中，T60混響時(shí)間單位為s；V為封閉空間的總體積，m3；c為水中聲波的傳播速度，c＝1 500 m／s；Si為內(nèi)表面的總面積，m2；為內(nèi)表面的平均吸聲系數(shù)。當(dāng)聲腔一定時(shí)，V、S均為定值，那么55.26 V／（cS）也為定值，則混響時(shí)間T60與聲腔內(nèi)表面的平均吸聲系數(shù)成反比。由此可知透聲窗面積所占的比例越大，相同條件下內(nèi)表面上的平均吸聲也越大，即混響時(shí)間越短，所以吸聲處理的效果就越差。圖7、圖8中透聲窗面積占總表面積的比例為1／4，吸聲處理的降噪效果比較明顯。圖8中當(dāng)透聲窗面積占總表面積的比例為3／4時(shí)，吸聲處理的降噪效果較差。

如圖9所示，透聲窗面積占總表面積比例為1／4時(shí)，敷設(shè)聲學(xué)阻尼材料和未敷設(shè)聲學(xué)阻尼材料的聲腔內(nèi)水動(dòng)力自噪聲隨頻率變化的趨勢(shì)基本一致，均隨著頻率的增大而減小，只是在個(gè)別頻率點(diǎn)有所突變，這和模型的固有頻率點(diǎn)重合。敷設(shè)阻尼材料后，各頻率下聲吶腔聲壓級(jí)均有降低，這是由于在板子系統(tǒng)中設(shè)置阻尼損耗因子后，降低了板子系統(tǒng)的輸出能量。且隨著阻尼損耗因子的增大，聲腔內(nèi)自噪聲聲壓級(jí)越來越小。此外可見，阻尼處理在中高頻段降噪效果更加明顯。

圖10所示為不同透聲窗面積占總表面積的比例下，不同阻尼處理對(duì)聲吶腔水動(dòng)力自噪聲的影響。橫軸為透聲窗所占的比例，縱軸為聲腔內(nèi)的總聲壓級(jí)。分析圖中的1號(hào)線，當(dāng)無聲學(xué)處理時(shí)，透聲窗面積所占的比例越大，聲腔自噪聲越小，這時(shí)主要靠透聲窗起到降噪作用。對(duì)比分析圖10中任意2條線之間的差值，當(dāng)透聲窗面積占總表面積的比例增大時(shí)，總聲壓級(jí)的差值越來越小，阻尼的降噪效果越來越有限，這是因?yàn)榉笤O(shè)阻尼的界面面積所占總表面積的比例越來越小。所以只有當(dāng)非透聲窗面積所占總表面積的比例較大時(shí)，阻尼處理才會(huì)有明顯的效果。圖中當(dāng)透聲窗面積占總表面積的比例為1／4時(shí)，阻尼處理的降噪效果明顯。當(dāng)透聲窗面積占總表面積的比例為3／4時(shí)，阻尼處理有降噪效果，但明顯不及透聲窗所占比例較小的情況。

圖9 相同透聲窗比例，不同阻尼處理（U0＝10 m／s）Fig.9 Same acoustic window proportion with different damping coating

圖10 不同透聲窗比例，不同阻尼處理（U0＝10 m／s）Fig.10 Different acoustic window proportion with different damping coating

4 計(jì)算結(jié)果的理論評(píng)價(jià)

以吸聲處理為例，分析不同透聲窗比例下，不同吸聲處理的降噪效果。

已知橢球聲吶腔體積V和總內(nèi)表面積S。根據(jù)式（2），聲吶腔內(nèi)表面平均吸聲公式α表達(dá)式為：

以非透聲窗界面上無聲學(xué)處理時(shí)為方案1，當(dāng)吸聲系數(shù)α為0.1～0.9時(shí)分別為：方案2，α＝0.1；方案3，α＝0.5；方案4，α＝0.8；方案5，α＝0.9。當(dāng)透聲窗面積所占比例為1／4時(shí)，5個(gè)方案中的橢球聲吶腔內(nèi)的平均吸聲系數(shù)、混響時(shí)間如表3所示。當(dāng)透聲窗面積所占比例為3／4時(shí)，同理可得5個(gè)方案下橢球聲吶腔內(nèi)的平均吸聲系數(shù)、混響時(shí)間。

表3 不同聲學(xué)處理方案下聲吶腔內(nèi)相關(guān)參數(shù)（透聲窗面積所占比例為1／4）Tab.3 The coefficients in sonar cavity under different acoustical treatments

從計(jì)算結(jié)果可知，兩種透聲窗面積所占比例下，混響時(shí)間均較短，反射聲衰減較快，聲吶腔內(nèi)的聲場(chǎng)主要以直達(dá)聲為主，而且當(dāng)透聲窗所占比例增大時(shí)，相同條件下混響時(shí)間更短，直達(dá)聲的貢獻(xiàn)相對(duì)更大。

封閉空間內(nèi)聲壓平方為：

式中，W為聲源的聲功率；ρ為流體密度；Q為指向性指數(shù)；c為流體中聲波的傳播速度；r為距離。

式（4）中第1項(xiàng)表示直達(dá)聲的貢獻(xiàn)，第2項(xiàng)表示混響聲的貢獻(xiàn)。透聲窗結(jié)構(gòu)的聲輻射可簡化為若干個(gè)簡單聲源的疊加，簡要分析橢球聲吶腔表面附近的聲場(chǎng)。取Q＝1，r＝0.2 m，按照式（4），大致可得到各方案總聲場(chǎng)的聲壓平方值 （單位為Pa2），采取聲學(xué)處理措施方案后聲吶腔內(nèi)平均降低的聲壓級(jí)如表4所示。

對(duì)比表4和圖8可知，理論評(píng)估結(jié)果和VAOne計(jì)算結(jié)果比較吻合，表明吸聲處理措施后，聲吶腔內(nèi)聲壓級(jí)的變化結(jié)果基本可信，即當(dāng)透聲窗面積所占總表面積的比例增大時(shí)，吸聲處理的降噪效果會(huì)弱化。

表4 聲學(xué)處理方案后聲吶腔內(nèi)平均降低聲壓級(jí)Tab.4 The average reduction value of sound pressure

5 結(jié)論

本文基于振動(dòng)聲學(xué)軟件VAOne，建立了簡化的水下橢球殼形狀的聲吶腔模型，運(yùn)用CFD軟件計(jì)算不同流速下，聲吶罩表面湍流邊界層的分離位置，考慮了湍流脈動(dòng)壓力空間分布的不均勻性。通過數(shù)值計(jì)算，分析了在相同和不同透聲窗面積所占總表面積的比例下，吸聲和阻尼處理對(duì)聲吶腔水動(dòng)力自噪聲的影響等。理論評(píng)估吸聲處理措施的結(jié)果與數(shù)值計(jì)算結(jié)果比較吻合。研究表明：

1）航速越大，聲吶腔內(nèi)水動(dòng)力自噪聲越大。當(dāng)航速增大1倍時(shí)，聲壓級(jí)大約增大20 dB，符合水動(dòng)力噪聲的規(guī)律。

2）無聲學(xué)處理和有聲學(xué)處理時(shí)，聲吶腔水動(dòng)力自噪聲總體上隨頻率增高而減小。

3）吸聲材料的吸聲系數(shù)越大，降噪效果越好；阻尼材料的阻尼因子越大，降噪效果越好；吸聲和阻尼處理在中高頻段降噪效果較好。

4）透聲窗面積所占總表面積的比例大小，對(duì)聲學(xué)處理的降噪效果有較大的影響。如果透聲窗面積所占總表面積的比例較大，則非透聲窗界面的吸聲處理的降噪效果就會(huì)弱化。只有當(dāng)透聲窗面積占總表面積的比例較小時(shí)，非透聲窗壁面上的聲反射在罩內(nèi)形成較強(qiáng)的混響時(shí)，吸聲處理才會(huì)有較好的降噪效果。

5）當(dāng)非透聲窗界面面積所占總表面積的比例較大時(shí)，阻尼處理的降噪效果更加明顯。

［1］LIGHTHILL M J.On sound generated aerodynamically.I.general theory［J］.Proceedings of the Royal Society A Mathematical，Physical&Engineering Sciences，1952，211（1107）：564－587.

［2］CURLE N.The influence of solid boundaries upon aerodynamic sound ［J］.Proceedings of the Royal Society A Mathematical，Physical&Engineering Sciences，1955，231（1187）：505－514.

［3］DOWELL E H.Transmission of noise from a turbulent boundary layer through a flexible plate into a closed cavity［J］.J.Acoust.Soc.Am.，1969，46（1）：238－252.

［4］LAUCHLE G C.Noise generated by axisymmetric turbulent boundary layer flow［J］.J.Acoust.Soc.Am.，1977，61（3）：694－702.

［5］俞孟薩，李東升.統(tǒng)計(jì)能量法計(jì)算聲吶自噪聲的水動(dòng)力噪聲分量［J］.船舶力學(xué)，2004，8（1）：99－105.

［6］彭臨慧，陸建輝.水下結(jié)構(gòu)物自噪聲統(tǒng)計(jì)能量分析［J］.中國造船，2001，42（3）：58－62.

［7］張娟，李天勻，劉敬喜，等.空間不均勻湍流邊界層激勵(lì)下聲腔自噪聲統(tǒng)計(jì)能量分析［J］.中國艦船研究，2007，2（6）：24－29.

［8］LYON R H.Statistical analysis of power injection and response in structures and rooms［J］.J.Acoust.Soc.Am.，1969，45（3）：545－565.

［9］HAN F，BERNHARD R J，MONGEAU L G.Prediction of flow－induced structural vibration and sound radiation using energy flow analysis［J］.Journal of Sound and Vibration，1999，227（4）：685－709.

［10］俞孟薩，葉劍平，吳有生，等.船舶聲吶部位自噪聲的預(yù)報(bào)方法及其控制技術(shù)［J］.船舶力學(xué)，2002，6（5）：80－94.

［11］MORSE P M.振動(dòng)與聲［M］.北京：科學(xué)出版社，1974.

Self Noise of Hydrodynamic Component for Ellipsoid-Shaped Sonar Cavity Induced by Turbulent Fluctuating Pressure

Peng Hua－kang Li Tian－yun Zhu Xiang
College of Naval Architecture and Ocean Engineering，Huazhong University of Science and Technology，Wuhan 430074，China

An underwater ellipsoid－shaped sonar cavity model was created by SEA method.The separation position of turbulent boundary layer was calculated with different flow speeds by CFD software，and the spatial distribution of non－uniform for turbulent fluctuating pressure was considered.This paper analyzed the influence of flow speeds and different absorption and damping coatings on the sonar flow－induced self noise.The results of theoretical evaluation for absorption coating agreed with the numerical results.It shows that the greater flow speed is，the greater sonar self noise is.And the smaller percentage the acoustic window area accounts for，the better performance noise reduction for absorption and damping coating will.The results could support the prediction and control of flow－induced noise for underwater sonar cavity.

statistical energy analysis；sonar cavity；turbulent fluctuating pressure；flow－induced self noise

U661.44

：A

：1673－3185（2011）04－37－06

2010-07-13

國家自然科學(xué)基金資助項(xiàng)目（40976058）

彭華康（1984－），男，碩士研究生。研究方向：振動(dòng)與噪聲控制。E-mail：penghuakang＠sohu.com

李天勻（1969－），男，教授，博士生導(dǎo)師。研究方向：結(jié)構(gòu)力學(xué)、振動(dòng)與噪聲控制。E-mail：ltyz801＠m(xù)ail.hust.edu.cn

10.3969/j.issn.1673-3185.2011.04.008