聲激勵簡單導(dǎo)管的聲輻射特性研究

鄭 晗，王 瑋，周其斗

（1.海軍裝備研究院，北京100161；2.海軍工程大學(xué) 艦船工程系，武漢430033）

聲激勵簡單導(dǎo)管的聲輻射特性研究

鄭 晗1，王 瑋1，周其斗2

（1.海軍裝備研究院，北京100161；2.海軍工程大學(xué) 艦船工程系，武漢430033）

采用NASTRAN無限元技術(shù)，對簡單導(dǎo)管受聲激勵的聲輻射特性進行了研究，建立了導(dǎo)管聲波模態(tài)和輻射聲能量之間的直接聯(lián)系，并考慮了導(dǎo)管半徑、長度、長徑比三種因素對管內(nèi)聲波模態(tài)、管外輻射聲功率和遠場聲指向性的影響。研究結(jié)果表明：影響導(dǎo)管聲波模態(tài)和輻射聲功率的參數(shù)主要有兩個，一是無因次波數(shù)，二是導(dǎo)管長徑比；導(dǎo)管的長徑比越小，總輻射聲能量變小，但側(cè)方屏蔽效應(yīng)也變差，因此需綜合考慮長徑比變化對總輻射聲能量和聲指向性的影響；和導(dǎo)管的側(cè)方屏蔽效應(yīng)主要作用于中高頻，低頻時導(dǎo)管外的聲場與入射聲場基本相同。

導(dǎo)管聲學(xué)；聲源激勵；單極子

0 引 言

導(dǎo)管受聲激勵下的聲輻射特性分析，對于優(yōu)化航空發(fā)動機、泵噴推進器和導(dǎo)管螺旋槳等導(dǎo)管結(jié)構(gòu)[1]的聲學(xué)設(shè)計，具有十分重要的意義。以往的研究中，對于導(dǎo)管聲學(xué)計算方法的關(guān)注較多[2-6]，對管內(nèi)外聲場和聲輻射特性的關(guān)注卻比較少。

Doak[7-8]推導(dǎo)了聲源分布于有限長導(dǎo)管內(nèi)的聲場聲壓和聲功率公式，指出對于剛性有凸緣導(dǎo)管，決定聲場的主要是三個因素：導(dǎo)管橫截面幾何、管內(nèi)聲源分布的時域和空間域特性和導(dǎo)管端口的邊界條件。Wang和Tszeng[9]針對有限長導(dǎo)管，推導(dǎo)了一種新的計算開口廣義輻射阻抗的公式，并指出，導(dǎo)管端口對整個聲場的影響在高頻、旋轉(zhuǎn)聲源及長導(dǎo)管情況下更加突出。這些研究主要從數(shù)學(xué)和公式推導(dǎo)的角度分析了導(dǎo)管聲場，對于管內(nèi)外聲場如何發(fā)生聯(lián)系，以及導(dǎo)管尺寸因素如何影響聲輻射特性，并沒有給出具體的分析。

本文以簡單剛性導(dǎo)管為研究對象，采用基于NASTRAN無限元技術(shù)的導(dǎo)管聲學(xué)計算方法，詳細闡述了管內(nèi)聲波模態(tài)與管外輻射聲功率的關(guān)系，并研究了導(dǎo)管尺寸對其降噪效果和屏蔽特性的影響，給出了具體的結(jié)論和建議。

1 模型和方法說明

如圖1（a）所示，一簡單剛性導(dǎo)管長為L，半徑為a，長徑比σ=L/D=L/2a，流體介質(zhì)為水，單極子聲源置于導(dǎo)管正中心，源強度A=1 kg/s2，計算頻段為0-4 000 Hz，有限元模型如圖1（b）。文獻[10]中已經(jīng)詳細說明了有限元建模的具體方法和計算參數(shù)設(shè)置，并驗證了采用NASTRAN無限元技術(shù)計算導(dǎo)管聲場的正確性，指出其優(yōu)勢在于適用于任意多連通域問題的聲場預(yù)報，適應(yīng)任意復(fù)雜的導(dǎo)管截面和形狀，并能考慮導(dǎo)管壁的彈性振動。

本文將基于該方法，分析簡單導(dǎo)管的聲輻射特性，并討論導(dǎo)管尺寸變化對管內(nèi)外聲場的影響規(guī)律。除特別說明，在研究半徑的影響時，固定導(dǎo)管長度L=1m，導(dǎo)管半徑分別取a=0.125m，0.25m，0.5 m，1m；研究長度的影響時，固定導(dǎo)管半徑a=0.5m，導(dǎo)管長度分別取L=0.5m、1m、2m和4m。

圖1 導(dǎo)管的幾何與有限元模型Fig.1 Geometry and FEMmodel of duct

2 導(dǎo)管內(nèi)聲場

表1 Jmn′（xmn）=0的解Tab.1 Solutions of Jmn′（xmn）=0

首先對傳統(tǒng)的無限長導(dǎo)管聲波模態(tài)理論進行簡要闡述，對一半徑為a的無限長導(dǎo)管，管內(nèi)聲壓的一個特解為[11]

其中：m和n分別為導(dǎo)管的周向和徑向模態(tài)數(shù)；Amn為聲壓幅值；Jm（kmnr ）表示宗量為（kmnr ）的m階柱貝塞爾函數(shù)，kmn=xmn/a，xmn是Jm′（xmn）=0的第n個零點（表1）；ω為聲源角頻率，kx為周向聲波數(shù)，并且滿足

當聲波圓頻率低于（m, n ）次簡正波對應(yīng)的簡正頻率，即k＜kmn時，kx為虛數(shù)，管內(nèi)只能傳播（m, n）階以下的聲波，此時也稱（m, n ）階聲波模態(tài)被截止，k＞kmn時稱為被打開。將臨界條件k=kmn轉(zhuǎn)化為ka=kmna=xmn，右邊的xmn為常數(shù)，可見決定無限長導(dǎo)管聲波模態(tài)打開或截止（即管內(nèi)聲場分布）的參數(shù)為ka。

公式（1）中，貝塞爾函數(shù)z=Jm（kmnr）=Jm（xmnr/a ）表示無限長導(dǎo)管內(nèi)任一截面處的聲壓分布，取周向模態(tài)數(shù)m=0，借助表1中的xmn值，就可以得到無限長導(dǎo)管內(nèi)第1-6階徑向聲波模態(tài)的云圖，如圖2所示，該云圖的形狀與半徑a無關(guān)。

圖2 無限長導(dǎo)管內(nèi)的聲場云圖Fig.2 Sound field contour inside the infinite-length duct

2.1 半徑的影響

下面研究有限長導(dǎo)管內(nèi)的聲場，對a=0.125 m，0.25m，0.5 m，1 m四種不同半徑導(dǎo)管，計算單極子聲源激勵下的管內(nèi)聲場。假定導(dǎo)管長為單位1，中心截面為x=0，左右兩端分別為x=-0.5和0.5，取x= 0.25截面處的導(dǎo)管內(nèi)聲場側(cè)視圖進行觀察。圖3給出了四種導(dǎo)管在取不同ka值時的管內(nèi)聲場截面云圖，它們的周向模態(tài)數(shù)均為0。

比較圖2和圖3可以發(fā)現(xiàn)，具有相同ka值的不同半徑導(dǎo)管的管內(nèi)聲場云圖形狀基本一致，并且與無限長導(dǎo)管內(nèi)聲場具有對應(yīng)關(guān)系，ka=0.63對應(yīng)J0（k00r）的平面波形式，ka=1.88對應(yīng)J0（k01r）的徑向（0,1）聲波模態(tài)……ka=7.54對應(yīng)J0（k04r）的徑向（0,4）聲波模態(tài)。這說明，有限長導(dǎo)管內(nèi)的聲場主要由參數(shù)ka決定，與導(dǎo)管半徑a和聲源頻率f均有關(guān)系，有限長導(dǎo)管內(nèi)的聲波模態(tài)仍然遵循傳統(tǒng)的無限長導(dǎo)管聲波模態(tài)理論。

需要說明的是，由于導(dǎo)管端口邊界的影響，在有限長導(dǎo)管中很難找到一個明顯的分界頻率點，以區(qū)分一個聲波模態(tài)和另一個聲波模態(tài)，這也正是有限長導(dǎo)管與無限長導(dǎo)管不同的地方，盡管管內(nèi)聲波模態(tài)的基本形式相同，但出現(xiàn)的頻率位置卻并不相同。

圖3 不同半徑的有限長導(dǎo)管內(nèi)聲場云圖Fig.3 Sound field contour inside the finite-length ductwith different radius

2.2 長度的影響

如圖4所示，給出了導(dǎo)管取四種不同長度L=0.5m、1m、2m和4m時的管內(nèi)聲場，四種導(dǎo)管的聲場云圖形狀沒有差別，這與Tyler和Sofrin在文獻[2]中的論述一致。他們指出，導(dǎo)管長度對于管內(nèi)聲場并不十分重要，它只在一個很窄的頻段內(nèi)發(fā)揮作用，當激勵頻率遠小于截止頻率，即f＜＜fq時，衰減波數(shù)為

壓力的衰減因子為p（x）=e-2πx/λq。這意味著，在長度等于截止頻率對應(yīng)波長λq=c/fq的導(dǎo)管里，壓力將下降20log（e2π），大約50 dB。

圖4 不同長度的導(dǎo)管內(nèi)聲場云圖ka=4.62Fig.4 Sound field contour inside the finite-length ductwith different length ka=4.62

3 聲源在導(dǎo)管中的輻射聲功率

剛性導(dǎo)管受聲源激勵時不會產(chǎn)生振動，因此不會直接向外輻射聲，本文中的“導(dǎo)管聲輻射”指的是聲源置于導(dǎo)管中時產(chǎn)生的輻射聲。通常，在計算結(jié)構(gòu)的輻射聲功率時，要在一個封閉面上進行數(shù)值積分。然而，導(dǎo)管是一個具有開口的多連通域結(jié)構(gòu)，不具有封閉結(jié)構(gòu)面，因此，在導(dǎo)管外人為地建立一個封閉的聲場球面，通過NASTRAN計算求出聲場節(jié)點的速度和聲壓，再在聲場球面上積分，就可以得到導(dǎo)管的輻射聲功率和聲功率級

式中：N為聲場球面的單元數(shù)目，pi和vi為第i個聲場單元上的平均聲壓和平均速度，“*”表示取共軛值，聲功率的參考值為Wref=10-12（W）。將計算得到的聲場結(jié)果文件（含有節(jié)點位移和速度）進行處理，采用FORTRAN語言編程，就可以得到導(dǎo)管的輻射聲功率，程序的正確性作者已通過與自由聲場中點聲源輻射聲功率的理論值進行了比較驗證，這里不再贅述。

3.1 半徑的影響

對a=0.125m、0.25m、0.5m和1m四種不同半徑的導(dǎo)管，計算單極子聲源激勵下的導(dǎo)管輻射聲功率值，聲功率圖以無因次波數(shù)ka為橫坐標，如圖5所示。同時給出的還有單極子聲源入射聲場的聲功率曲線，其理論值[11]為，聲功率級66.22 dB。觀察聲功率曲線，可以發(fā)現(xiàn)：

圖5 不同半徑導(dǎo)管的輻射聲功率Fig.5 Radiated power of ductswith different radius

（1）從整體上看，以ka為橫坐標的四條聲功率曲線形狀和趨勢十分相似，均經(jīng)歷了三個主波峰和兩個主波谷，且波峰和波谷的位置基本重合，好像可以用一條曲線（如a=1 m曲線）擬合出四種不同半徑導(dǎo)管聲功率隨ka變化的趨勢。這說明，ka是決定導(dǎo)管輻射聲功率的一個非常重要的參數(shù)，因為ka決定了導(dǎo)管內(nèi)聲波模態(tài)的基本形式和第（m,n）階模態(tài)出現(xiàn)的具體位置，所以也就決定了導(dǎo)管輻射聲功率主波峰和波谷出現(xiàn)的位置。

（2）再看曲線的局部特征，以a=0.125 m曲線為例說明，對三個主峰和前兩個主峰之間的4個次峰作垂線相交于橫軸，找到對應(yīng)的頻率值，并提取相應(yīng)的導(dǎo)管軸向聲波模態(tài)，如圖6所示。可以看到，導(dǎo)管輻射聲功率的峰值與軸向聲波模態(tài)數(shù)（記為p），即管內(nèi)聲場的軸向波峰數(shù)目有著直接關(guān)系。模態(tài)數(shù)為1時，輻射聲功率出現(xiàn)主峰；模態(tài)數(shù)為3、5、7等奇數(shù)時，輻射聲功率出現(xiàn)次峰，并且峰值隨著模態(tài)數(shù)的增大而減小。這是因為，有限長導(dǎo)管內(nèi)的聲波要受到管口邊界反射的影響，因此管內(nèi)軸向聲波模態(tài)會產(chǎn)生不同數(shù)目的波峰，當模態(tài)數(shù)為1時，管內(nèi)聲波受邊界反射的影響最小，聲波被反射回來的比例最小，輻射效率最高，因而輻射聲功率的峰值最大；當模態(tài)數(shù)增大為3、5、7時，聲波被反射回來的比例增大，輻射效率降低，因而輻射聲功率峰值也依次減小。波谷位置對應(yīng)的管內(nèi)軸向聲波模態(tài)數(shù)為2、4等偶數(shù)，如圖7所示，此時管內(nèi)正中位置的聲壓值最小，與聲源發(fā)出的聲波相互抵消，輻射效率最低，輻射聲功率也最小，因此出現(xiàn)聲功率曲線上的波谷。

圖6 輻射聲功率波峰對應(yīng)的導(dǎo)管軸向聲波模態(tài)Fig.6 Axial acoustic wavemodes of duct corresponding to the peak of the radiated power curve

圖7 輻射聲功率波谷對應(yīng)的導(dǎo)管軸向模態(tài)Fig.7 Axial acoustic wavemodes of duct corresponding to the peak of the radiated power curve

繪圖表示聲功率曲線波峰、波谷與軸向模態(tài)數(shù)p的對應(yīng)關(guān)系，如圖8所示，隨著聲源頻率的增大，導(dǎo)管軸向聲波模態(tài)數(shù)p先從1開始增大到最大值，然后又逐漸減小到1，對應(yīng)的聲功率曲線從一個主峰值轉(zhuǎn)到另一個主峰值，中間還經(jīng)歷了多個峰值和谷值。進一步的觀察還可以發(fā)現(xiàn)，在相鄰的兩個主波峰之間，盡管管內(nèi)聲波的軸向模態(tài)數(shù)p不盡相同，但徑向和周向模態(tài)數(shù)（m,n）卻是相同的，這也再次印證了上一條結(jié)論。

（3）比較四條輻射聲功率曲線的主峰值大小，可以看到，導(dǎo)管的半徑越小，聲功率的峰值就越大。這是因為，點聲源的輻射聲功率要受到環(huán)境的影響，根據(jù)文獻[12]中的公式推導(dǎo)，當一個點聲源十分靠近無限大剛性面時，點聲源的輻射聲功率是自由空間中輻射聲功率的兩倍，聲功率級高約3 dB，即聲功率的大小受環(huán)境的影響。由此可以解釋這個現(xiàn)象：聲源在半徑小的細導(dǎo)管中受到端口邊界條件的影響更大，聲功率級相比入射聲場的增加值更大。

（4）比較四條輻射聲功率曲線的峰值數(shù)目，可以看到，導(dǎo)管的半徑越小，聲功率曲線的次波峰和波谷數(shù)目越多，曲線的振蕩越厲害。這是因為半徑小的細導(dǎo)管相比粗導(dǎo)管具有更為豐富的軸向聲波模態(tài)，軸向模態(tài)數(shù)p的取值更多，由于每個p值都對應(yīng)一個波峰或者波谷，因而峰值的數(shù)目更多。

3.2 長度的影響

圖9給出了四種不同長度導(dǎo)管的輻射聲功率曲線，它與圖5中的曲線形狀十分相似，L=0.5 m對應(yīng)前面的a=1 m，L=1 m對應(yīng)前面的a=0.5m……這些對應(yīng)的曲線都具有同樣的長徑比，由此猜測，直接影響輻射聲功率的導(dǎo)管參數(shù)不是半徑，也不是長度，而是導(dǎo)管的長徑比。下面將驗證這個猜測。

3.3 長徑比的影響

對前面兩個系列計算中具有相同長徑比的導(dǎo)管輻射聲功率進行比較，如圖10所示，可以看到：具有相同長徑比的兩導(dǎo)管輻射聲功率曲線形狀完全一致，由此驗證長徑比是真正影響輻射聲功率的導(dǎo)管參數(shù)，3.1節(jié)中關(guān)于半徑對導(dǎo)管聲輻射影響的結(jié)論可全部轉(zhuǎn)換為長徑比的結(jié)論，對主要結(jié)論歸納總結(jié)如下：

（1）聲源在剛性導(dǎo)管中的輻射聲功率是一個隨頻率變化的量，其輻射聲功率曲線與管內(nèi)聲場的內(nèi)在聯(lián)系表現(xiàn)為：無因次波數(shù)ka的大小決定了管內(nèi)聲場的形態(tài)，即管內(nèi)聲波模態(tài)數(shù)（m,n,p），其中徑向和周向模態(tài)數(shù)（m,n）決定了輻射聲功率出現(xiàn)主峰值的位置，即曲線的整體特征；而軸向模態(tài)數(shù)p決定了輻射聲功率出現(xiàn)次波峰和次波谷的位置，即曲線的局部特征。（2）導(dǎo)管的長徑比越大，其輻射聲功率的主峰值越大，峰值的數(shù)目也越多。

圖8 輻射聲功率曲線的波峰波谷與模態(tài)數(shù)p的對應(yīng)關(guān)系Fig.8 The relationship between the peak of radiated power curve and mode number p

圖9 不同長度導(dǎo)管的輻射聲功率Fig.9 Radiated power of ductswith different length

圖10 不同長徑比導(dǎo)管的輻射聲功率Fig.10 Radiated power of ductswith different ratio of length to diameter

4 導(dǎo)管外聲場指向性

導(dǎo)管的輻射聲功率反映的是聲源在導(dǎo)管中輻射聲能量的總大小，而聲指向性則反映了聲源在導(dǎo)管中發(fā)聲形成聲場的方向性特征。設(shè)計一個聲學(xué)性能優(yōu)化的導(dǎo)管，不僅要抑制總能量的大小以降低噪聲，還要關(guān)注一些特定方向的聲壓級，發(fā)揮屏蔽作用。本節(jié)將討論管外聲場的指向性，固定導(dǎo)管長度l= 1m，導(dǎo)管半徑分別取為a=0.125m、0.25m、0.5m和1m，研究長徑比σ變化對管外聲場指向性的影響。

選取圖5中第1～3個主波峰和第1～2個主波谷對應(yīng)的無因次波數(shù)ka=0.2、4、7.3和ka=3.5、7，繪出不同長徑比的導(dǎo)管在10m半徑處的遠場聲指向性圖，如圖11所示。其中，無導(dǎo)管存在時的入射聲場聲壓），聲壓級96.99 dB。

觀察可以發(fā)現(xiàn)：（1）輻射聲功率主波峰位置對應(yīng)的聲場指向性圖（a）、（c）、（e）中，聲壓級大于入射聲場聲壓級；而主波谷位置對應(yīng)的聲場指向性圖（b）、（d）則恰好相反。由此可見，采用輻射聲功率和遠場聲壓級來判斷導(dǎo)管的降噪效果在一定程度上是等效的。（2）聲場指向性圖（b）、（d）中的90°方向聲壓級小于自由聲場聲壓級，這一現(xiàn)象被稱之為“側(cè)方屏蔽效應(yīng)”[3]。顯然，σ=4導(dǎo)管的側(cè)方屏蔽效應(yīng)最好，并且依次好于σ=2，1，0.5導(dǎo)管。這說明，導(dǎo)管的長徑比越大，其側(cè)方屏蔽效應(yīng)也越好，因此可以通過增大長徑比的方法抑制導(dǎo)管90°方向的噪聲。（3）比較相同ka值下的不同長徑比導(dǎo)管外聲場指向性，長徑比越大，管外聲場指向性越復(fù)雜，波瓣數(shù)目也越多，這是因為細導(dǎo)管相比粗導(dǎo)管擁有更豐富的軸向聲波模態(tài)。（4）當聲源頻率較低時，四種導(dǎo)管的聲場指向性圖十分接近，并且與點聲源的入射聲場一樣沒有指向性。這是因為聲源的頻率很小，其波長比導(dǎo)管半徑大很多，導(dǎo)管的存在對聲傳播的影響很小，大部分入射聲波可以繞過導(dǎo)管傳播，發(fā)生聲波的繞射（或衍射）現(xiàn)象，因此表現(xiàn)出很弱的指向性。隨著頻率的增加，導(dǎo)管端口邊界的影響逐漸增強，輻射聲場表現(xiàn)出較強的指向性，此時對噪聲的預(yù)報和控制都將變得更加困難。

圖11 不同波數(shù)激勵下的聲指向性圖Fig.11 Directivity of sound with differentwave number at sound field radius 10m

5 結(jié)論與展望

（1）有限長導(dǎo)管內(nèi)的聲波模態(tài)仍然遵循傳統(tǒng)的無限長導(dǎo)管聲波模態(tài)理論，但由于端口輻射阻抗的影響，同一階聲波模態(tài)出現(xiàn)的頻率位置不同。導(dǎo)管半徑對管內(nèi)聲場的影響表現(xiàn)為，具有相同ka值的不同半徑導(dǎo)管的管內(nèi)聲場云圖形狀基本一致，且與無限長導(dǎo)管內(nèi)聲場云圖一一對應(yīng)，導(dǎo)管長度對管內(nèi)聲場的影響較小，主要表現(xiàn)在遠低于截止頻率的頻段。

（2）聲源在剛性導(dǎo)管中的輻射聲功率是一個隨頻率變化的量，其輻射聲功率曲線與管內(nèi)聲場的內(nèi)在聯(lián)系表現(xiàn)為：無因次波數(shù)ka的大小決定了管內(nèi)聲場的形態(tài)，即管內(nèi)聲波模態(tài)數(shù)（m,n,p），其中徑向和周向模態(tài)數(shù)（m,n）決定了輻射聲功率曲線出現(xiàn)主峰值的位置，即曲線的整體特征；而軸向模態(tài)數(shù)p決定了輻射聲功率曲線出現(xiàn)次波峰和次波谷的位置，即曲線的局部特征。

（3）導(dǎo)管的長徑比越小，輻射聲功率的主峰值越小，峰值數(shù)目也越少，對于抑制總輻射聲能量的大小越有利；但是，長徑比小的導(dǎo)管側(cè)方屏蔽效應(yīng)差，對于抑制導(dǎo)管90°方向的聲壓級卻不利；因此，在導(dǎo)管的聲學(xué)設(shè)計中，必須兼顧長徑比變化對總輻射聲能量和某方向聲壓級大小的影響，使導(dǎo)管總的降噪效果和特定方向上的屏蔽特性達到相對最優(yōu)。

（4）導(dǎo)管內(nèi)外聲場形狀的復(fù)雜程度與兩個因素有關(guān)，一是聲源頻率，二是導(dǎo)管的長徑比，頻率和長徑比的增大都將使管內(nèi)聲波模態(tài)數(shù)和管外聲場的波瓣數(shù)增多。

（5）導(dǎo)管的側(cè)方屏蔽效應(yīng)主要作用于中高頻，低頻時導(dǎo)管外的聲場與入射聲場基本相同。

下一步，可利用簡單導(dǎo)管聲學(xué)特性的結(jié)論，研究復(fù)雜導(dǎo)管（如導(dǎo)管螺旋槳、泵噴推進器、分支型導(dǎo)管）的聲學(xué)特性，并提出具體的降噪措施。

[1]劉忠族,孫玉東,吳有生.管道流固耦合振動及聲傳播的研究現(xiàn)狀及展望[J].船舶力學(xué),2001,5(2):82-90. Liu Zhongzu,Sun Yudong,Wu Yousheng.Current situation and trends on the study of coupled fluid-structure vibration and sound propagation of pipeline systems[J].Journal of Ship Mechanics,2001,5(2):82-90.

[2]Tyler JM,Sofrin TG.Axial flow compressor noise studies[J].SAE Transaction,1962,70:309-332.

[3]MyersM K.Radiation of sound from a point source in a short duct[C].NASA CP-3352,1997:19-26.

[4]楊迎春,周其斗,潘雨村.有限長導(dǎo)管聲場預(yù)報的一種邊界元方法[J].聲學(xué)學(xué)報,2011,36(1):60-65. Yang Yingchun,Zhou Qidou,Pan Yucun.Boundary elementmethod for sound field prediction of finite duct[J].Acta A-custica,2011,36(1):60-65.

[5]楊迎春,周其斗,盧東慶.一種計算內(nèi)外域聲場的簡單源方法[J].海軍工程大學(xué)學(xué)報,2011,23(3):67-71. Yang Yingchun,Zhou Qidou,Lu Dongqing.Simple sourcemethod for sound field prediction in coupled interior/exterior domain[J].Journal of Naval University of Engineering,2011,23(3):67-71.

[6]柳貢民,李艷華.考慮流固耦合的多分支管道系統(tǒng)振動研究[J].船舶力學(xué),2012,16(5):532-541. Liu Gongmin,Li Yanhua.Research on vibration ofmulti-branched pipe system considering fluid-structure interaction[J]. Journal of Ship Mechanics,2012,16(5):532-541.

[7]Doak P E.Excitation,transmission and radiation of sound from source distributions in hard-walled ducts of finite length (I):The effects of duct cross-section geometry and source distribution space-time pattern[J].Journal of Sound and Vibration,1973,31(1):1-72.

[8]Doak PE.Excitation,transmission and radiation of sound from source distributions in hard-walled ducts of finite length (II):The effects of duct length[J].Journal of Sound and Vibration,1973,31(2):137-174.

[9]Wang K S,Tszeng T C.Propagation and radiation of sound in a finite length duct[J].Journal of Sound and Vibration, 1984,93(1):57-59.

[10]鄭 晗,周其斗,謝志勇.一種基于NASTRAN無限元技術(shù)的導(dǎo)管聲學(xué)計算方法[J].海軍工程大學(xué)學(xué)報,2015,27(1): 68-73. Zheng Han,Zhou Qidou,Xie Zhiyong.A calculation method for duct acoustics based on NASTRAN infinite element technology[J].Journal of Naval University of Engineering,2015,27(1):68-73.

[11]杜功煥,朱哲民,龔秀芬.聲學(xué)基礎(chǔ)[M].南京:南京大學(xué)出版社,2004.

[12]Ver IL,Beranek L L.Noise and vibration control engineering[M].USA:Wiley,2005.

Study of sound radiation of a simple duct excited by sound source

ZHENG Han1,WANGWei1,ZHOU Qi-dou2
(1.Naval Academy of Armament,Beijing 100161,China;2.Dept.of Naval Architecture Engineering, Naval Univ.of Engineering,Wuhan 430033,China)

Based on NASTRAN infinite element technology,the sound radiation of a simple duct excited by sound source was studied detailedly.The relationship between wavemode and radiated energy was founded, and the influence of the duct’s radius,length,length-to-diameter ratio on wavemode,radiated power and directivity of far field was considered.The results show that:there are two parameters that affect duct’s wavemode and radiated power,the non-dimensional wave number and the length-to-diameter ratio;and with the ratio becomes smaller,the radiated power becomes smaller,and the lateral shielding effect becomes weaker.It should be considered comprehensively to choose proper size of duct for the best effects of sound power reduction and lateral shielding；the lateral shielding effects act on high frequency band rather than low frequency band.

ductacoustics;acoustic excitation;monopole

U661.44

：Adoi:10.3969/j.issn.1007-7294.2017.02.013

2016-07-08

總裝備部“十二五”預(yù)研項目（9140A14080512JB11165）

鄭 晗（1988-），男，博士，工程師，E-mail:zhenghanhust@qq.com；

王 瑋（1981-），男，碩士；周其斗（1962-），男，教授，博士生導(dǎo)師。

1007-7294（2017）02-0228-09