張?jiān)娍?朱海潮 毛榮富 蘇俊博 蘇常偉
(海軍工程大學(xué)振動(dòng)與噪聲研究所1) 武漢 430033) (船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室2) 武漢 430033)
強(qiáng)干擾環(huán)境下雙面Patch近場(chǎng)聲全息技術(shù)*
張?jiān)娍?,2)朱海潮1,2)毛榮富1,2)蘇俊博1,2)蘇常偉1,2)
(海軍工程大學(xué)振動(dòng)與噪聲研究所1)武漢 430033) (船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室2)武漢 430033)
針對(duì)目標(biāo)聲源在強(qiáng)干擾環(huán)境下進(jìn)行聲場(chǎng)重建易受到干擾源、散射效應(yīng)影響及測(cè)量孔徑效應(yīng)的雙重問(wèn)題,提出強(qiáng)干擾聲場(chǎng)環(huán)境下基于聲輻射模態(tài)聲場(chǎng)分離算法的雙面Patch近場(chǎng)聲全息技術(shù).通過(guò)理論分析和公式推導(dǎo),建立了基于聲輻射模態(tài)聲場(chǎng)分離算法的內(nèi)插外推數(shù)學(xué)模型,不僅清除干擾源聲場(chǎng)及散射聲場(chǎng)的影響,而且減少測(cè)量工作量,加快計(jì)算速度;數(shù)值仿真及實(shí)驗(yàn)表明,強(qiáng)干擾環(huán)境下忽略散射聲場(chǎng)的影響很難重建目標(biāo)聲源的聲場(chǎng),且隨著干擾源強(qiáng)度增大分離誤差逐漸增大,而利用文中究的理論同時(shí)剔除了散射聲場(chǎng)及干擾源聲場(chǎng)的影響,提高了聲場(chǎng)重建的精度,降低了測(cè)量成本.
聲全息;強(qiáng)干擾環(huán)境;聲輻射模態(tài);聲場(chǎng)分離
近場(chǎng)聲全息技術(shù)在聲場(chǎng)重建過(guò)程中保留了高波數(shù)空間的倐逝波成分,提高了重建精度,從全息領(lǐng)域脫穎而出.然而,常規(guī)的近場(chǎng)聲全息技術(shù)在測(cè)量過(guò)程中不僅對(duì)聲場(chǎng)環(huán)境要求嚴(yán)格,即要求所有聲源位于測(cè)量面一側(cè),以滿足測(cè)量面另一側(cè)為自由場(chǎng)條件,而且為了保證重建精度,通常要求測(cè)量面的測(cè)量孔徑必須大于實(shí)際聲源的尺寸[1].但是,實(shí)際的測(cè)量環(huán)境比較復(fù)雜,測(cè)量面的兩側(cè)或都存在聲源,或存在其他物體的散射、反射,若直接進(jìn)行測(cè)量會(huì)產(chǎn)生較大的重建誤差.此時(shí),聲源識(shí)別及聲場(chǎng)預(yù)測(cè)的關(guān)鍵在于能否將目標(biāo)聲源從復(fù)雜聲場(chǎng)環(huán)境中分離出來(lái).另外,對(duì)于大尺寸結(jié)構(gòu)聲源的測(cè)量,不僅測(cè)量工作量較大、計(jì)算速度慢,而且測(cè)量成本較高.于是,在針對(duì)以上2個(gè)問(wèn)題的研究中形成了基于近場(chǎng)聲全息的聲場(chǎng)分離技術(shù)及基于內(nèi)插外推的Patch技術(shù).
目前在非自由聲場(chǎng)環(huán)境下聲場(chǎng)分離技術(shù)的應(yīng)用研究中形成了基于空間聲場(chǎng)變換理論、邊界元法、波疊加法、統(tǒng)計(jì)最優(yōu)法等多種算法聲場(chǎng)分離技術(shù).畢傳興等[2]利用空間Fourier變換算法提出基于質(zhì)點(diǎn)振速的雙面聲場(chǎng)分離技術(shù),在近場(chǎng)聲全息聲場(chǎng)重建中取得很好的效果.但是,空間Fourier變換算法在計(jì)算中不僅要求測(cè)量面的形狀與聲源面共形,而且測(cè)量面尺寸容易帶來(lái)窗效應(yīng)和卷繞誤差.Valdivia等[3]利用基于邊界元法的雙面聲場(chǎng)分離技術(shù)實(shí)現(xiàn)了目標(biāo)聲源的聲場(chǎng)分離及重建.基于邊界元法的聲場(chǎng)分離技術(shù)雖然能夠?qū)崿F(xiàn)任意形狀聲源的聲場(chǎng)分離及重建,但是邊界源法不僅存在特征波數(shù)處解的非唯一性問(wèn)題,而且隨著聲源尺寸增大,表面網(wǎng)格數(shù)目劃分較多,計(jì)算量增大.畢傳興等[4]基于等效源法利用雙測(cè)量面上的質(zhì)點(diǎn)振速數(shù)據(jù)實(shí)現(xiàn)目標(biāo)聲源的速度場(chǎng)的分離及重建.等效源法雖然克服了奇異積分和解的非唯一性問(wèn)題,但是要求測(cè)點(diǎn)數(shù)目必須大于等效源數(shù)目并存在逆問(wèn)題處理問(wèn)題.Fernandez等[5]采用基于統(tǒng)計(jì)最優(yōu)的雙面質(zhì)點(diǎn)振速聲場(chǎng)分離技術(shù)實(shí)現(xiàn)目標(biāo)聲源的速度場(chǎng)分離及重建,并與單面聲壓-質(zhì)點(diǎn)振速聲場(chǎng)分離技術(shù)作對(duì)比,證明其優(yōu)越性.統(tǒng)計(jì)最優(yōu)算法雖避免了窗效應(yīng)及卷繞誤差,十分適合局部聲場(chǎng)的重建,但是無(wú)法實(shí)現(xiàn)復(fù)雜結(jié)構(gòu)聲源的聲場(chǎng)重建.BraiKia等[6]利用球面波疊加方法不僅去除了球體聲源外干擾聲源的影響,而且去除了干擾聲源在目標(biāo)聲源表面散射聲的影響,大大提高了聲場(chǎng)重建的精度.球面波疊加方法適合于類球形聲源,對(duì)于其他形狀聲源聲場(chǎng)重建效果不好.另外,綜合考慮聲場(chǎng)分離技術(shù)與基于數(shù)據(jù)內(nèi)插與外推Patch技術(shù)的研究并不多,畢傳興等[7]研究了非自由場(chǎng)環(huán)境下基于空間聲場(chǎng)變換理論的雙面Patch聲場(chǎng)分離技術(shù),對(duì)非自由聲場(chǎng)環(huán)境下局部聲源的識(shí)別進(jìn)行了初步研究,減少了測(cè)量工作量,為大尺寸結(jié)構(gòu)聲源內(nèi)部聲場(chǎng)的重建提供了理論基礎(chǔ),但是仍需要進(jìn)一步研究各參數(shù)與誤差間的關(guān)系.文中綜合考慮以上所提2個(gè)問(wèn)題及各種算法的優(yōu)缺點(diǎn),將聲輻射模態(tài)應(yīng)用于非自由場(chǎng)環(huán)境下目標(biāo)聲源的聲場(chǎng)分離與重建,既具有等效源法和邊界元法適用于復(fù)雜結(jié)構(gòu)聲源聲場(chǎng)重建的優(yōu)點(diǎn),又克服了它們的缺陷,而且利用聲輻射模態(tài)的空間濾波特性,實(shí)現(xiàn)利用聲場(chǎng)中少量測(cè)點(diǎn)的數(shù)據(jù)準(zhǔn)確剔除干擾源影響的同時(shí)為進(jìn)一步的聲場(chǎng)重建提供較多測(cè)點(diǎn)全息數(shù)據(jù)的目的.
在強(qiáng)干擾聲場(chǎng)環(huán)境下,干擾聲源在目標(biāo)聲源表面產(chǎn)生的散射效應(yīng)明顯,實(shí)際的測(cè)量環(huán)境比較復(fù)雜,此時(shí)的聲場(chǎng)環(huán)境應(yīng)包含目標(biāo)聲源的聲場(chǎng)、干擾源產(chǎn)生的聲場(chǎng)及散射聲場(chǎng).目標(biāo)聲源的聲場(chǎng)分離時(shí)只考慮干擾源聲場(chǎng)而忽視散射聲場(chǎng)的影響,會(huì)給目標(biāo)聲源聲場(chǎng)的分離帶來(lái)很大的誤差,尤其在散射聲場(chǎng)影響較為明顯時(shí),會(huì)導(dǎo)致目標(biāo)源聲場(chǎng)分離的失敗.所以,在強(qiáng)干擾聲場(chǎng)環(huán)境下實(shí)現(xiàn)目標(biāo)源聲場(chǎng)的分離,不得不同時(shí)剔除干擾源聲場(chǎng)及散射聲場(chǎng)的影響.文中充分考慮強(qiáng)干擾環(huán)境下聲場(chǎng)的特點(diǎn),從確保聲場(chǎng)重建精度,減少測(cè)量成本入手,運(yùn)用源強(qiáng)聲輻射模態(tài)理論及Patch近場(chǎng)聲全息技術(shù),推導(dǎo)出了強(qiáng)干擾環(huán)境下目標(biāo)源聲場(chǎng)分離公式,同時(shí)剔除了干擾源聲場(chǎng)及散射聲場(chǎng)的影響,為聲場(chǎng)重建提供準(zhǔn)確數(shù)據(jù).數(shù)值仿真及實(shí)驗(yàn)研究對(duì)研究理論進(jìn)行了正確性驗(yàn)證.
1.1 目標(biāo)聲源分離技術(shù)
置于密度為ρ0、聲速為c0均質(zhì)流體中的復(fù)雜結(jié)構(gòu)以角頻率ω振動(dòng),向結(jié)構(gòu)外空間輻射聲能量.在邊界表面S0上滿足Neumann邊界條件,在無(wú)窮遠(yuǎn)處滿足Sommerfeld輻射條件,向外輻射的聲場(chǎng)滿足Helmhotlz方程,在省略時(shí)間相關(guān)項(xiàng)e-jωt后,由聲輻射模態(tài)理論[8],可得振源在全息面處的聲壓矩陣與振速矩陣
P=GΦC=ΨC
(1)
(2)
式中:G,D為聲場(chǎng)中N個(gè)點(diǎn)與聲源表面上M個(gè)點(diǎn)之間的傳遞矩陣;Φ為聲輻射模態(tài);Ψ,Ψv分別為聲壓聲場(chǎng)分布模態(tài)及振速聲場(chǎng)分布模態(tài);C為聲源的聲場(chǎng)分布模態(tài)對(duì)應(yīng)的展開系數(shù)矩陣.當(dāng)全息面位置及聲源結(jié)構(gòu)形狀確定的條件下,G、Φ、Ψ、Ψv都是可直接求取的量,聲壓計(jì)算時(shí)可認(rèn)為是已知的量,見圖1.
圖1 雙全息面與目標(biāo)聲源、干擾聲源的空間位置關(guān)系圖
由圖1可知,目標(biāo)聲源和干擾聲源在測(cè)量面兩側(cè),采用雙測(cè)量面測(cè)量聲場(chǎng)聲壓,測(cè)量面SH1上測(cè)得的聲壓數(shù)據(jù)記為p1;測(cè)量面SH2上測(cè)得的聲壓數(shù)據(jù)記為p2.在強(qiáng)干擾聲場(chǎng)環(huán)境下,干擾聲源在目標(biāo)聲源表面產(chǎn)生的散射效應(yīng)比較顯著,全息面上測(cè)量的聲壓由目標(biāo)聲源的聲場(chǎng)、干擾源產(chǎn)生的聲場(chǎng)及干擾聲源在目標(biāo)聲源表面產(chǎn)生的散射聲場(chǎng)組成.其中,目標(biāo)聲源聲場(chǎng)與散射聲場(chǎng)向外輻射聲壓記為pout;干擾聲源向內(nèi)入射的聲壓記為pin.可以得到
(3)
(4)
聯(lián)立式(3),(4)可得聲輻射模態(tài)對(duì)應(yīng)的展開系數(shù)為
(5)
(6)
根據(jù)求得的系數(shù)C1可以求出測(cè)量面SH1上向外聲場(chǎng)產(chǎn)生的聲壓為
(7)
(8)
在目標(biāo)聲源表面,聲壓和質(zhì)點(diǎn)振速的關(guān)系為
(9)
(10)
(11)
(12)
(13)
將式(10)~(13)與式(9)聯(lián)合可以得到C3的求解公式.
(15)
當(dāng)干擾源的強(qiáng)度較小時(shí),可忽略散射聲場(chǎng)的影響,則全息面上的聲壓由目標(biāo)源聲場(chǎng)及干擾源聲場(chǎng)聲壓組成,從而計(jì)算得到目標(biāo)聲源在全息面SH1上的聲壓為
(16)
1.2 基于內(nèi)插與外推的Patch技術(shù)
全息面上的聲壓與聲場(chǎng)分布模態(tài)有關(guān),而聲場(chǎng)分布模態(tài)又與全息面上的測(cè)點(diǎn)分布有關(guān).倏逝波與高階聲輻射模態(tài)有關(guān),在不影響獲取充足倏逝波的條件下,利用最優(yōu)模態(tài)截止方法對(duì)模態(tài)階數(shù)進(jìn)行截取,從而可對(duì)基于數(shù)據(jù)內(nèi)插與外推的Patch近場(chǎng)聲全息技術(shù)進(jìn)行研究.
假設(shè)目標(biāo)聲源在全息面SH1上的模態(tài)截止階數(shù)為M0,全息面上的測(cè)點(diǎn)數(shù)目為N1,則目標(biāo)聲源在全息面上的聲壓可表示為
p1(N1×1)=Ψ1(N1×M0)C1(M0×1)
(17)
式中:Ψ1(N1×M0)分別為目標(biāo)聲源在全息面上的聲壓聲場(chǎng)分布模態(tài),模態(tài)截止階數(shù)為M0;C1(M0×1)為對(duì)應(yīng)的聲場(chǎng)分布聲輻射模態(tài)展開系數(shù).
擴(kuò)大全息面的測(cè)量點(diǎn)數(shù),使之具有N個(gè)測(cè)點(diǎn),且N大于N1,則可得目標(biāo)聲源在N個(gè)測(cè)點(diǎn)測(cè)量面上聲壓為
ppatch=Ψ1(N×M0)C1(M0×1)
(18)
式中:Ψ1(N×M0)為目標(biāo)聲源在N個(gè)測(cè)點(diǎn)全息面上的聲壓聲場(chǎng)分布模態(tài).
聯(lián)立式(17)與式(18)消去C1(M0×1)可得
ppatch=Ψ1(N×M0)(Ψ1(N1×M0))+p1(N1×1)
(19)
式中:ppatch為全息面上N個(gè)測(cè)點(diǎn)的聲壓數(shù)據(jù);p1(N1×1)為全息面上N1個(gè)測(cè)點(diǎn)的聲壓數(shù)據(jù);Ψ1(N×M0),Ψ1(N1×M0)為可以計(jì)算的量,若已知全息面上N1個(gè)測(cè)點(diǎn)的p1(N1×1)的數(shù)據(jù)就可以預(yù)估全息面上N個(gè)測(cè)點(diǎn)ppatch的數(shù)據(jù),即利用全息面上少量測(cè)點(diǎn)的聲壓數(shù)據(jù)預(yù)估更多測(cè)點(diǎn)聲壓的數(shù)據(jù).
為驗(yàn)證本文研究理論的正確性,在測(cè)量面兩側(cè)采用不同聲源進(jìn)行仿真分析,假定目標(biāo)聲源為受簡(jiǎn)諧激勵(lì)的簡(jiǎn)支板,干擾源為剛性脈動(dòng)球,半徑為5 cm.簡(jiǎn)支平板參數(shù)如下:平板尺寸為0.5 m×0.5 m,厚度為8 mm,泊松比υ取0.28,楊氏模量E為2×1011Pa,密度ρ為7 800 kg/m3.數(shù)值仿真計(jì)算中,測(cè)量面SH1和SH2大小皆為1×1 m,其上均布13×13個(gè)測(cè)量點(diǎn),目標(biāo)聲源與全息測(cè)量面SH1距離為d1=0.02 m,干擾源與測(cè)量面SH2距離為d2=0.05 m,2全息面之間的間距為Δd=0.03 m.以簡(jiǎn)支板中心為坐標(biāo)原點(diǎn)建立空間直角坐標(biāo)系,干擾源小球的球心坐標(biāo)位于點(diǎn)(0,0,d1+d2+Δd)處,仿真中假設(shè)鋼板表面是剛性邊界,則脈動(dòng)球在鋼板表面產(chǎn)生的散射聲場(chǎng)相當(dāng)于在(0,0,-d1-d2-Δd)處相同的脈動(dòng)小球產(chǎn)生的聲場(chǎng).
現(xiàn)定義干擾聲源與目標(biāo)聲源在全息面SH1上的輻射聲功率比為
(20)
目標(biāo)源聲場(chǎng)聲壓的分離誤差定義為
(21)
2.1 理論正確性驗(yàn)證
取激振頻率f=600 Hz,作用在平板(x0,y0)=(0.125, 0.125)的位置處.為了使仿真結(jié)果接近真實(shí)測(cè)量情況,對(duì)測(cè)量面的聲壓數(shù)據(jù)加入了30 dB的隨機(jī)白噪聲.此次仿真為模擬強(qiáng)干擾的聲場(chǎng)環(huán)境,采用的輻射聲功率比為0 dB,即干擾聲源與目標(biāo)聲源在全面SH1上的聲源強(qiáng)度相同.
圖2 測(cè)量面SH1上的聲壓幅值分布圖
從圖2可知,直接測(cè)量所得聲壓的高聲壓區(qū)域位置發(fā)生偏移,所有區(qū)域聲壓數(shù)據(jù)皆被放大,測(cè)量誤差高達(dá)126.45%;忽略散射聲場(chǎng)進(jìn)行聲場(chǎng)分離所得聲壓幅值分布圖像的中部區(qū)域出現(xiàn)偽聲壓,分離誤差達(dá)42.37%;剔除散射聲場(chǎng)進(jìn)行聲場(chǎng)分離所得聲壓幅值分布圖像與理論聲壓幅值分布圖像非常相似,較好的反映目標(biāo)源聲場(chǎng)的實(shí)際情況,分離誤差僅為12.96%;對(duì)剔除散射聲場(chǎng)分離所得聲壓數(shù)據(jù)進(jìn)行2倍內(nèi)插與外推處理,patch誤差為13.47%.
數(shù)值仿真表明,剔除散射聲場(chǎng)進(jìn)行聲場(chǎng)分離能夠進(jìn)一步提高目標(biāo)源聲場(chǎng)的重建精度;對(duì)分離數(shù)據(jù)進(jìn)行內(nèi)插外推后,依然能夠反映目標(biāo)源聲場(chǎng)的實(shí)際情況,圖像漸變區(qū)域更加圓滑,誤差并增大幅度可以忽略不計(jì).
強(qiáng)干擾環(huán)境下的干擾源強(qiáng)度、激勵(lì)頻率等關(guān)鍵參數(shù)對(duì)分離誤差會(huì)產(chǎn)生一定的影響,下面分析參數(shù)變化對(duì)分離結(jié)果產(chǎn)生的具體影響.
2.2 干擾源強(qiáng)度的影響
調(diào)整剛性脈動(dòng)小球表面振速,改變干擾聲源的輻射聲功率,實(shí)現(xiàn)全息面SH1上干擾聲源與目標(biāo)聲源的輻射聲功率比調(diào)節(jié),從而模擬不同干擾源強(qiáng)度下的聲場(chǎng)環(huán)境,見圖3.
圖3 誤差曲線圖
如圖3所示,剔除散射和忽略散射2種情況下,不同干擾源強(qiáng)度環(huán)境下的誤差曲線圖.由圖可以看出,當(dāng)輻射聲功率比較小,即目標(biāo)源的輻射聲功率遠(yuǎn)大于干擾源的輻射聲功率時(shí),此時(shí)散射聲場(chǎng)的影響并不顯著,剔除散射及忽略散射進(jìn)行聲場(chǎng)分離都能取得很好的效果.當(dāng)輻射聲功率比較大時(shí),干擾源的輻射強(qiáng)度增大,散射聲場(chǎng)的影響開始凸顯,忽略散射進(jìn)行聲場(chǎng)分離誤差變大,且隨著聲功率比的增大而增大,若聲功率比較大則很難將目標(biāo)聲源的聲場(chǎng)分離出來(lái);而剔除散射進(jìn)行聲場(chǎng)分離時(shí),隨著聲功率比的增大,誤差略有增大,但增加幅度非常?。徊逯岛髉atch誤差比插值前誤差略小,但兩者相差不大,變化趨勢(shì)相似,patch技術(shù)能夠有效實(shí)現(xiàn)強(qiáng)干擾聲場(chǎng)環(huán)境下目標(biāo)聲源的聲場(chǎng)分離.由此可以看出,強(qiáng)干擾聲場(chǎng)環(huán)境下,散射聲場(chǎng)影響顯著,聲場(chǎng)分離時(shí)必須同時(shí)剔除干擾源聲場(chǎng)及散射聲場(chǎng)的影響才能有效重建目標(biāo)源的聲場(chǎng).
2.3 激勵(lì)頻率的影響
輻射聲功率比為-1 dB,其他參數(shù)與上文仿真中保持一致.頻率在較大范圍內(nèi)進(jìn)行遍歷選取,采用雙面Patch技術(shù)對(duì)陣列信息進(jìn)行處理,見圖4.
圖4 誤差隨激勵(lì)頻率的變化曲線
不考慮干擾源影響直接測(cè)量誤差大部分在75%以上,忽略散射聲場(chǎng)分離誤差大部分在40%以上,剔除散射聲場(chǎng)分離誤差始終低于20%,可以看出該方法對(duì)干擾聲源的影響具有很強(qiáng)的抑制作用;從整個(gè)曲線變化趨勢(shì)上可以看出,隨著激勵(lì)頻率的增高誤差有增大趨勢(shì),中低頻段的分離效果較好而高頻段的分離效果較差,這是因?yàn)楦哳l段的波長(zhǎng)較短,測(cè)量所需測(cè)點(diǎn)間距更小,而實(shí)際仿真中要保持測(cè)點(diǎn)間距不變,所以采集到的聲場(chǎng)信息并不充分,測(cè)量誤差在逆運(yùn)算過(guò)程中被放大,致使分離誤差增大,但分離效果仍然處于可接受范圍之內(nèi).
以固支平板為目標(biāo)聲源,音箱為干擾聲源進(jìn)行實(shí)驗(yàn)研究,驗(yàn)證理論的正確性.固支平板是大小為0.5 m×0.5 m、厚度為0.001 m的鋼板,由信號(hào)發(fā)生器產(chǎn)生一路信號(hào)送入JZK-2型激振器激勵(lì)鋼板,產(chǎn)生穩(wěn)態(tài)聲場(chǎng).音箱聲源尺寸為0.13 m×0.14 m×0.23 m,信號(hào)發(fā)生器經(jīng)同一通道同時(shí)發(fā)出激振信號(hào)和音響信號(hào).全息面的測(cè)量網(wǎng)格點(diǎn)數(shù)為13(X向)×13(Y向),2個(gè)方向上的傳聲器的間距均為5 cm.測(cè)量時(shí)設(shè)有15個(gè)通道,前13通道為陣列傳聲器,第14通道為參考傳聲器,第15通道為參考加速度傳感器.以平板中心為坐標(biāo)原點(diǎn)建立空間直角坐標(biāo)系,激勵(lì)點(diǎn)坐標(biāo)為(0.125,-0.125,0),板對(duì)角布置參考加速度傳感器.平板、音箱、測(cè)量面布置時(shí)保持三者的中心點(diǎn)在一條直線上.第一次測(cè)量時(shí),線陣從右至左進(jìn)行掃描(以測(cè)量者為基準(zhǔn),測(cè)量者面對(duì)目標(biāo)聲源),測(cè)量面與目標(biāo)聲源面的距離為0.04 m,目標(biāo)聲源面與干擾聲源面的距離為0.5 m,采樣頻率2 048 Hz,采樣時(shí)間間隔為1 s,利用單參考源傳遞函數(shù)法計(jì)算各測(cè)點(diǎn)復(fù)聲壓并作為全息面SH1上的聲壓數(shù)據(jù).第二次測(cè)量時(shí),測(cè)量面與目標(biāo)聲源面的距離為0.08 m,其他保持不變,測(cè)量計(jì)算的復(fù)聲壓作為全息面SH2上的聲壓數(shù)據(jù).同理,設(shè)置全息面的測(cè)量網(wǎng)格點(diǎn)數(shù)7(X向)×7 (Y向),其他設(shè)置方法及參數(shù)保持不變,研究patch技術(shù).實(shí)驗(yàn)現(xiàn)場(chǎng)見圖5.
圖5 實(shí)驗(yàn)現(xiàn)場(chǎng)圖
圖6為137 Hz時(shí)測(cè)量面SH1上的聲壓幅值分布圖.實(shí)驗(yàn)時(shí),首先用激振器對(duì)平板進(jìn)行單獨(dú)激勵(lì),激勵(lì)頻率為137 Hz,利用線陣進(jìn)行掃描,設(shè)置測(cè)量面與目標(biāo)聲源面距離為0.04 m,得到13×13測(cè)點(diǎn)下平板聲場(chǎng)的理論聲壓數(shù)據(jù)幅值分布圖,見圖6b);然后打開音箱,頻率與目標(biāo)聲源的激勵(lì)頻率保持一致,再用線陣進(jìn)行掃描,得到13×13測(cè)點(diǎn)下干擾環(huán)境下混合聲場(chǎng)的聲壓數(shù)據(jù)p1幅值分布圖,見圖6a),直接測(cè)量偽聲壓區(qū)域增多,很難重建目標(biāo)聲源聲場(chǎng),測(cè)量誤差達(dá)328.56%;調(diào)整測(cè)量面與目標(biāo)聲源面距離為0.08 m,得到13×13測(cè)點(diǎn)下干擾環(huán)境下混合聲場(chǎng)的聲壓數(shù)據(jù)p2.利用p1,p2聲壓數(shù)據(jù)進(jìn)行聲場(chǎng)分離,得到忽略散射進(jìn)行聲場(chǎng)分離全息面SH1上的聲壓幅值分布圖,見圖6c),偽聲壓區(qū)域增加,并不能確定目標(biāo)聲源的峰值準(zhǔn)確信息,此時(shí)分離誤差達(dá)76.34%;剔除散射進(jìn)行聲場(chǎng)分離全息面SH1上的聲壓幅值分布圖,見圖6d),重建目標(biāo)聲源的聲場(chǎng)與目標(biāo)聲源的理論聲場(chǎng)十分吻合,分離誤差僅為23.38%.同理,利用7×7個(gè)測(cè)點(diǎn)重復(fù)上述步驟,剔除散射進(jìn)行聲場(chǎng)分離并利用patch技術(shù)預(yù)估13×13測(cè)點(diǎn)下的全息面SH1上的聲壓數(shù)據(jù),幅值分布圖見圖6e),與理論聲場(chǎng)基本吻合,較好反映了目標(biāo)源聲場(chǎng)的位置信息,patch誤差為26.47%.
圖6 137 Hz時(shí)測(cè)量面SH1上的聲壓幅值分布圖
對(duì)實(shí)驗(yàn)結(jié)果的有效性進(jìn)行進(jìn)一步的研究,分別對(duì)不同的頻率進(jìn)行了研究,各類誤差見表1.
表1 不同頻率下的相關(guān)誤差
1) 強(qiáng)干擾聲場(chǎng)環(huán)境下,利用直接測(cè)量的聲壓數(shù)據(jù)很難重建目標(biāo)聲源的聲場(chǎng);忽略散射進(jìn)行聲場(chǎng)分離依然很大,不能反映目標(biāo)聲源的聲場(chǎng)信息;只有同時(shí)剔除散射聲場(chǎng)及干擾源聲場(chǎng)的影響才能比較準(zhǔn)確的重建目標(biāo)聲源的聲場(chǎng).
2) 隨著干擾源強(qiáng)度增大,即干擾聲源與目標(biāo)聲源的聲功率比增加,剔除散射及忽略散射下的分離誤差都會(huì)增大,但剔除散射下的分離誤差增加幅度非常小,而忽略散射下的分離誤差則快速增加;隨著激勵(lì)頻率的增高誤差有增大趨勢(shì),中低頻段的分離效果較好而高頻段的分離效果較差.
3) 數(shù)值仿真及實(shí)驗(yàn)表明,強(qiáng)干擾聲場(chǎng)環(huán)境下基于聲輻射模態(tài)聲場(chǎng)分離理論的patch技術(shù)在減少測(cè)量工作量,提高計(jì)算速度的同時(shí),實(shí)現(xiàn)了目標(biāo)聲源的聲場(chǎng)重建.
[1]陳心昭,畢傳興.近場(chǎng)聲全息技術(shù)及其應(yīng)用[M].北京:科學(xué)出版社,2013.
[2]畢傳興,張永斌,徐亮,等.基于雙面質(zhì)點(diǎn)振速測(cè)量的聲場(chǎng)分離技術(shù)[J].聲學(xué)學(xué)報(bào),2010,35(6):653-658.
[3]VALDIVIA N P, WILLIAMS E G, HERDIC P C. Approximations of inverse boundary element methods with partial measurements of the pressure field[J]. Journal of the Acoustics Society of America,2008(1):109-120.
[4]畢傳興,胡定玉,張永斌,等.基于等效源法和雙面質(zhì)點(diǎn)振速測(cè)量的聲場(chǎng)分離方法[J].物理學(xué)報(bào),2013,62(8):284-301.
[5]FERNANDEZ G E F,JACOBSEN F. Sound field separation with a double layer velocity transducer array[J].Journal of the Acoustics Society of America,2013(1):5-8.
[6]BRAIKIA Y, MELON M, LANGRENNE C, et al. Evaluation of a separation method for source identification in small spaces[J]. Journal of the Acoustics Society of America,2013(1):323-331.
[7]畢傳興,唐琦,張永斌,等.雙面Patch近場(chǎng)聲全息實(shí)驗(yàn)研究[J].合肥工業(yè)大學(xué)學(xué)報(bào),2011,34(3):321-325.
[8]聶永發(fā),朱海潮.利用源強(qiáng)密度聲輻射模態(tài)重建聲場(chǎng)[J].物理學(xué)報(bào),2014,63(10):256-267.
[9]王冉,陳進(jìn),董廣明.基于改進(jìn)HELS方法的局部近場(chǎng)聲全息技術(shù)研究[J].振動(dòng)與沖擊,2014,33(10):157-161.
Double-layer Patch Near-field Acoustic Holography under Strong Interference Environment
ZHANG Shike1,2)ZHU Haichao1,2)MAO Rongfu1,2)SU Junbo1,2)SU Changwei1,2)
(InstituteofNoise&vibration,NavalUniversityofEngineering,Wuhan430033,China)1)(NationalKeyLaboratoryonShipVibration&Noise,Wuhan430033,China)2)
Aiming at the twin problems that sound field reconstruction of the target sound source is easy to be disturbed by the interference sound source, scattering sound field and the measuring aperture effect, the Double-layer Patch near-field acoustic holography based on acoustic field separation algorithm of acoustic radiation modes under strong interference environment is presented. Firstly, the internal difference extrapolation mathematical model based on acoustic field separation algorithm of acoustic radiation modes is established by theoretical analysis and formula deduction. It not only eliminates the influence of the interference sound source and the scattered sound field, but also reduces the workload of the measurement and speeds up the calculation speed. Then the numerical simulation and experimental results verify the theoretical correctness. Numerical simulation and experimental results show that neglecting the influence of the scattered field would make it difficult to reconstruct the sound field of the target source under strong interference environment and the separation error increases with the increase of the intensity of the interference source. By using the theory of this paper, the influence of scattering sound field and the sound field is eliminated, the reconstruction accuracy is improved and the cost is reduced.
acoustic holography; strong interference environment; acoustic radiation modes; acoustic field separation
2016-10-13
*國(guó)家自然科學(xué)基金項(xiàng)目資助(51305452)
TB532
10.3963/j.issn.2095-3844.2016.06.017
張?jiān)娍?1990—):男,碩士生,主要研究領(lǐng)域?yàn)檎駝?dòng)與噪聲控制