蘇常偉，朱海潮，毛榮富

(1. 海軍工程大學(xué) 船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室， 湖北 武漢 430033；2. 中國人民解放軍92578部隊(duì)， 北京 100161)

研究結(jié)構(gòu)振動(dòng)向封閉聲腔中輻射噪聲這一課題具有實(shí)際意義，例如船舶艙室及飛機(jī)座艙內(nèi)噪聲分析和控制。國內(nèi)外學(xué)者針對此類腔體結(jié)構(gòu)聲輻射問題開展了大量的研究[1-5]，其中Snyder等[5]通過對封閉聲腔聲勢能解耦，首次提出了耦合封閉聲腔的聲輻射模態(tài)(Acoustic Radiation Mode, ARM)概念并初步探討了利用其進(jìn)行有源結(jié)構(gòu)聲控制(Active Structural Acoustic Control, ASAC)的可能性。耦合封閉聲腔的聲輻射模態(tài)消除了利用結(jié)構(gòu)模態(tài)表示封閉聲腔聲勢能時(shí)復(fù)雜的耦合項(xiàng)，而且在中、低頻時(shí)通常只有一階或幾階模態(tài)對封閉聲腔聲勢能的貢獻(xiàn)起主要作用，這使得分析和控制此類聲輻射問題更為簡單，成為近年來結(jié)構(gòu)聲學(xué)領(lǐng)域的研究熱點(diǎn)[6-8]。

在實(shí)際設(shè)計(jì)基于聲輻射模態(tài)的耦合封閉聲腔ASAC系統(tǒng)時(shí)，如果能夠精準(zhǔn)確定對聲勢能有主要貢獻(xiàn)的聲輻射模態(tài)，那么在選取控制目標(biāo)時(shí)就具有非常明確的針對性，這直接影響控制效果，而且能夠有效減少控制通道數(shù)。然而，文獻(xiàn)[5-8]中僅依據(jù)輻射效率選擇前一階或前幾階聲輻射模態(tài)進(jìn)行控制，自由空間ASAC中的相關(guān)研究[9-10]也大都強(qiáng)調(diào)輻射效率的作用而忽視模態(tài)幅度的影響，這顯然是不全面的。李雙等[11]、和衛(wèi)平等[12]指出利用聲輻射模態(tài)和結(jié)構(gòu)模態(tài)之間的內(nèi)在聯(lián)系，可以有效確定主導(dǎo)聲輻射模態(tài)，但該方法需要用到結(jié)構(gòu)模態(tài)信息，而實(shí)際中要準(zhǔn)確獲得輻射體的結(jié)構(gòu)模態(tài)信息并不容易，因此這一方法在應(yīng)用上存在不便之處。截至目前，還沒有一種準(zhǔn)確的、便于工程應(yīng)用的確定方法。

本文在前期研究中，將封閉聲腔聲勢能直接表示為結(jié)構(gòu)振動(dòng)表面法向速度的二次型形式，重新定義了耦合封閉聲腔的聲輻射模態(tài)概念[13]。在此基礎(chǔ)上，本文進(jìn)一步給出了耦合封閉聲腔中主導(dǎo)聲輻射模態(tài)的定義；然后全面考慮輻射效率和模態(tài)幅度的作用，提出了一種基于“初選—預(yù)留—后驗(yàn)”的主導(dǎo)聲輻射模態(tài)確定方法；最后通過實(shí)驗(yàn)研究驗(yàn)證了上述方法的有效性和可行性。

1 耦合封閉聲腔的主導(dǎo)聲輻射模態(tài)

文獻(xiàn)[13]定義的耦合封閉聲腔的聲輻射模態(tài)是結(jié)構(gòu)振動(dòng)表面法向速度的一組可能分布形式。因此，結(jié)構(gòu)振動(dòng)表面法向速度向量v可由這組模態(tài)線性疊加表示，即

(1)

式中：D為聲輻射模態(tài)矩陣，滿足正交性(即DTD=E)，其列向量di(i=1,2,…,M)即為重新定義的耦合封閉聲腔的聲輻射模態(tài)；y為聲輻射模態(tài)(速度)幅度向量，其元素yi(i=1,2,…,M)對應(yīng)第i階聲輻射模態(tài)幅度；M為結(jié)構(gòu)振動(dòng)表面的離散單元數(shù)目。

封閉聲腔聲勢能Ep也可進(jìn)一步寫成

(2)

式中，A為對角矩陣，其對角線元素λi(i=1,2,…,M)為誤差加權(quán)矩陣Ω的特征值，均為正實(shí)數(shù)，與第i階模態(tài)的輻射效率成正比，稱之為第i階模態(tài)對應(yīng)的輻射效率系數(shù)。由式(2)可得，耦合封閉聲腔的聲輻射模態(tài)能夠?qū)崿F(xiàn)聲勢能解耦。因此，只要減小任意一階聲輻射模態(tài)幅度，就能保證總的封閉聲腔聲勢能降低。

將在某一頻率下或某一頻段內(nèi)對封閉聲腔聲勢能起主要作用的一階或幾階聲輻射模態(tài)定義為耦合封閉聲腔的主導(dǎo)聲輻射模態(tài)(Dominant Acoustic Radiation Mode, DARM)。在中、低頻時(shí)，一般只需抵消主導(dǎo)聲輻射模態(tài)貢獻(xiàn)的聲勢能，就能獲得較為理想的有源控制效果。

2 主導(dǎo)聲輻射模態(tài)確定方法

在實(shí)際設(shè)計(jì)ASAC系統(tǒng)時(shí)，首要環(huán)節(jié)就是要準(zhǔn)確認(rèn)定主導(dǎo)聲輻射模態(tài)，從而正確選取控制目標(biāo)，以保證控制效果。由式(2)可得，各階聲輻射模態(tài)單獨(dú)貢獻(xiàn)的聲勢能Epi為

Epi=λi|yi|2，i=1,2,…,M

(3)

可見，各階聲輻射模態(tài)單獨(dú)貢獻(xiàn)的聲勢能為該階模態(tài)幅度幅值的平方與輻射效率系數(shù)的乘積。已有的主導(dǎo)聲輻射模態(tài)確定方法大都強(qiáng)調(diào)輻射效率的作用而忽視模態(tài)幅度的影響，這顯然是不全面的。

為直觀說明，以一尺寸為Lx×Ly×Lz=0.65 m×0.50 m×0.45 m、頂面為四邊固支的彈性鋼板(尺寸為Lx×Ly=0.65 m×0.50 m，板厚h=0.002 6 m)、其余5個(gè)面均為剛性壁面的矩形聲腔為例進(jìn)行仿真分析。腔內(nèi)介質(zhì)為空氣，設(shè)彈性板左下角為坐標(biāo)原點(diǎn)，將彈性鋼板表面均勻離散為15×15個(gè)小面積元。假設(shè)彈性板受到集中點(diǎn)力激勵(lì)，激勵(lì)頻率為150 Hz，激勵(lì)位置為(0.152 m, 0.117 m)。圖1～3分別給出了該激勵(lì)下此耦合板-腔模型各階聲輻射模態(tài)的輻射效率系數(shù)、模態(tài)幅度的幅值及其單獨(dú)貢獻(xiàn)的聲勢能隨階數(shù)變化情況(為便于觀察，這里只顯示前20階)。分析圖1～3可知：如果某階聲輻射模態(tài)的模態(tài)幅度很小，即使該階模態(tài)的輻射效率很高，其輻射聲勢能還是很小，即該階模態(tài)不能作為主導(dǎo)聲輻射模態(tài)，比如第1階模態(tài)；反之亦然，比如第12階模態(tài)。

圖1 前20階聲輻射模態(tài)對應(yīng)的輻射效率系數(shù)Fig.1 Radiation efficiency coefficients of the first twenty ARMs

圖2 前20階聲輻射模態(tài)幅度幅值Fig.2 Magnitudes of modal amplitude of the first twenty ARMs

圖3 前20階聲輻射模態(tài)單獨(dú)貢獻(xiàn)的聲勢能Fig.3 Acoustic potential energy of the first twenty ARMs

基于此，本文全面考慮了輻射效率和模態(tài)幅度兩個(gè)因素，提出了一種基于“初選—預(yù)留—后驗(yàn)”的主導(dǎo)聲輻射模態(tài)確定方法，具體如下：

1)計(jì)算輻射效率，初選可能的主導(dǎo)聲輻射模態(tài)。單頻分析時(shí)，計(jì)算各階模態(tài)對應(yīng)的輻射效率系數(shù)，初步選擇輻射效率系數(shù)不小于最大值1%的聲輻射模態(tài)作為可能的主導(dǎo)聲輻射模態(tài)[6]；寬頻分析時(shí)，繪制各階輻射效率系數(shù)-頻率曲線，將在分析頻段內(nèi)出現(xiàn)峰值的最高模態(tài)階數(shù)作為主導(dǎo)聲輻射模態(tài)的初選階數(shù)，記為js。

2)考慮模態(tài)幅度作用，放大初選階數(shù)以預(yù)留余量。設(shè)預(yù)留余量后的主導(dǎo)聲輻射模態(tài)階數(shù)為JS，則JS=μ×js。其中，μ為放大因子，一般取整，根據(jù)經(jīng)驗(yàn)賦初值。

3)傳感前JS階聲輻射模態(tài)幅度，直接對比聲勢能后驗(yàn)初選結(jié)果。具體可分成以下兩個(gè)步驟：

①在結(jié)構(gòu)振動(dòng)表面布置傳感器陣列測量振速數(shù)據(jù)，傳感前JS階模態(tài)幅度。在結(jié)構(gòu)表面均勻離散的前提下，可類比自由空間聲輻射模態(tài)[14]，將不同自由度下的聲輻射模態(tài)向量均轉(zhuǎn)換為標(biāo)準(zhǔn)聲輻射模態(tài)(聲輻射模態(tài)函數(shù)的離散形式)，以保證振型的一致性。因此，可以利用低自由度下的標(biāo)準(zhǔn)聲輻射模態(tài)設(shè)計(jì)前JS階模態(tài)幅度傳感策略，即

(4)

(5)

②將利用式(4)～(5)獲得的前JS階模態(tài)幅度代入式(3)，計(jì)算前JS階模態(tài)單獨(dú)貢獻(xiàn)的聲勢能Epi(i=1,2,…,JS)，直接觀察Epi隨階數(shù)收斂情況，最終確定主導(dǎo)聲輻射模態(tài)。

上述確定方法的流程圖如圖4所示。不難發(fā)現(xiàn)，該方法綜合運(yùn)用計(jì)算和測試手段，本質(zhì)上是通過直接對比前JS階聲輻射模態(tài)單獨(dú)貢獻(xiàn)的聲勢能確定主導(dǎo)聲輻射模態(tài)，一方面保證了主導(dǎo)聲輻射模態(tài)的精準(zhǔn)認(rèn)定，另一方面通過“初選”環(huán)節(jié)進(jìn)行了模態(tài)階數(shù)截?cái)?，較大程度地減小了工作量，而且整個(gè)過程中不需要用到輻射體的結(jié)構(gòu)模態(tài)信息，便于工程實(shí)現(xiàn)。

圖4 主導(dǎo)聲輻射模態(tài)確定流程Fig.4 Determination process of DARM

值得一提的是，利用該方法確定主導(dǎo)聲輻射模態(tài)后，不僅能夠指導(dǎo)ASAC中控制目標(biāo)選取，還可以用于實(shí)現(xiàn)封閉聲腔聲勢能重構(gòu)，即

(6)

3 實(shí)驗(yàn)研究

為驗(yàn)證主導(dǎo)聲輻射模態(tài)確定方法的有效性及可行性，按照第2節(jié)中尺寸設(shè)計(jì)、制作固支彈性板-矩形腔耦合模型，開展實(shí)驗(yàn)研究。采用HEV-50型激振器向彈性板施加帶寬為30～300 Hz的白噪聲激勵(lì)，激勵(lì)位置為(0.155 m, 0.05 m)。將彈性板表面均勻離散為13×11個(gè)小面積元，圖5給出了該激勵(lì)頻段內(nèi)前6階聲輻射模態(tài)的輻射效率系數(shù)隨頻率變化曲線。

圖5 前6階模態(tài)的輻射效率系數(shù)隨頻率變化曲線Fig.5 Radiation efficiency coefficients of the first six ARMs as a function of frequency

首先，根據(jù)圖4中流程逐步確定該激勵(lì)下的主導(dǎo)聲輻射模態(tài)，具體為：觀察圖5中給出的輻射效率系數(shù)-頻率曲線可得，初選前2階模態(tài)作為可能的主導(dǎo)聲輻射模態(tài)，即js=2；為充分考慮模態(tài)幅度的作用，將初選階數(shù)放大3倍以預(yù)留余量，即令μ=3，則JS=6；根據(jù)式(4)～(5)給出的模態(tài)幅度傳感策略，假設(shè)將彈性板表面均勻離散為4×3個(gè)小面積元，利用均布在面元中心的12個(gè)加速度傳感器測量彈性板表面法向速度，傳感此耦合板-腔的前6階聲輻射模態(tài)幅度，然后代入式(3)計(jì)算得到前6階模態(tài)單獨(dú)貢獻(xiàn)的聲勢能。前6階聲輻射模態(tài)幅度的幅值及單獨(dú)貢獻(xiàn)的聲勢能分別如圖6和圖7所示。觀察圖7中各階聲勢能-頻率曲線可得：在不同頻率下，不同階次的聲輻射模態(tài)對總的封閉聲腔聲勢能起主導(dǎo)作用，比如第1階聲輻射模態(tài)在30～100 Hz時(shí)對聲勢能起主要貢獻(xiàn)；從整個(gè)分析頻段來看，對封閉聲腔聲勢能起主要貢獻(xiàn)的是第1～3階聲輻射模態(tài)，第4～6階模態(tài)的貢獻(xiàn)量則相對較小。因此，可最終確定該激勵(lì)下此耦合板-腔模型的主導(dǎo)聲輻射模態(tài)為第1～3階聲輻射模態(tài)。后續(xù)控制中，若能有效抵消第1～3階主導(dǎo)聲輻射模態(tài)貢獻(xiàn)的聲勢能，就可以獲得較為理想的寬頻(30～300 Hz)控制效果。

(a) 第1階(a) The first order (b) 第2階(b) The second order

(c) 第3階(c) The third order (d) 第4階(d) The fourth order

(e) 第5階(e) The fifth order (f) 第6階(f) The sixth order圖6 前6階聲輻射模態(tài)幅度的幅值隨頻率變化曲線Fig.6 Magnitudes of modal amplitude of the first six ARMs as a function of frequency

(a) 第1階(a) The first order (b) 第2階(b) The second order

(c) 第3階(c) The third order (d) 第4階(d) The fourth order

(e) 第5階(e) The fifth order (f) 第6階(f) The sixth order圖7 前6階聲輻射模態(tài)的聲勢能隨頻率變化曲線Fig.7 Acoustic potential energy of the first six ARMs as a function of frequency

利用4個(gè)布置在矩形腔底部四個(gè)角落的MPA416型參考傳聲器估算封閉聲腔聲勢能[15]，記為測量值。同時(shí)，將第1～3階主導(dǎo)聲輻射模態(tài)代入式(6)重構(gòu)封閉聲腔聲勢能，記為重構(gòu)值。將兩者對比，以驗(yàn)證上述確定結(jié)果的正確性，結(jié)果如圖8所示?？梢?，分別經(jīng)重構(gòu)和測量得到的聲勢能曲線除在波谷處由于信噪比低略有差別外，總體吻合良好，尤其在波峰處兩者基本重合。這就證明了上述確定結(jié)果的正確性，從而表明本文提出的主導(dǎo)聲輻射模態(tài)確定方法是有效且可行的。

最后，分析圖5～7中第3階聲輻射模態(tài)的輻射效率系數(shù)、模態(tài)幅度的幅值及其單獨(dú)貢獻(xiàn)的聲勢能隨頻率變化曲線，可以發(fā)現(xiàn)：在162～186 Hz頻段內(nèi)，雖然該階模態(tài)的輻射效率相比第1、2階模態(tài)低，但其模態(tài)幅度的幅值較大，因而此時(shí)該階模態(tài)單獨(dú)輻射的聲勢能在總聲勢能中的占比明顯大于第1、2階模態(tài)，成為該頻段內(nèi)的主導(dǎo)聲輻射模態(tài)。這就再次說明了在確定主導(dǎo)聲輻射模態(tài)時(shí)綜合考慮輻射效率和模態(tài)幅度兩個(gè)因素的必要性，這也正是本文所提方法的意義所在。

圖8 聲勢能對比Fig.8 Comparison of acoustic potential energy

4 結(jié)論

在前期研究的基礎(chǔ)上，定義了耦合封閉聲腔的主導(dǎo)聲輻射模態(tài)，提出了一種基于“初選—預(yù)留—后驗(yàn)”的主導(dǎo)聲輻射模態(tài)確定方法，并以固支板-腔耦合模型為例開展了實(shí)驗(yàn)研究，驗(yàn)證了上述方法的有效性和可行性，得到結(jié)論如下：①提出的主導(dǎo)聲輻射模態(tài)確定方法綜合運(yùn)用計(jì)算和測試手段、全面考慮輻射效率和模態(tài)幅度兩個(gè)因素，能夠精準(zhǔn)認(rèn)定主導(dǎo)聲輻射模態(tài)；②該方法通過“初選”環(huán)節(jié)較大程度地減小了工作量，且整個(gè)過程中不需要用到輻射體的結(jié)構(gòu)模態(tài)信息，便于工程實(shí)現(xiàn)；③該方法可用于指導(dǎo)ASAC中控制目標(biāo)選取或重構(gòu)封閉聲腔聲勢能，而且“初選—預(yù)留—后驗(yàn)”的主導(dǎo)聲輻射模態(tài)確定思路在自由空間中同樣適用。