華大成,彭臨慧,于小濤

(中國(guó)海洋大學(xué)信息科學(xué)與工程學(xué)院,山東青島266100)

水中彈性圓柱體聲散射實(shí)驗(yàn)研究

華大成,彭臨慧,于小濤

(中國(guó)海洋大學(xué)信息科學(xué)與工程學(xué)院,山東青島266100)

對(duì)水中彈性圓柱進(jìn)行聲散射實(shí)驗(yàn)研究,結(jié)合彈性體的共振散射理論驗(yàn)證實(shí)驗(yàn)結(jié)果。在實(shí)驗(yàn)中得到較好的時(shí)域散射波形,與理論得出的時(shí)間間隔吻合,并利用準(zhǔn)諧波法得到目標(biāo)的形態(tài)函數(shù),該結(jié)果與理論一致。文章使用維格納威利(Wigner-Ville)分布來(lái)分析散射信號(hào)的時(shí)頻關(guān)系,得到時(shí)頻聯(lián)合分布圖,從中看出散射信號(hào)中存在頻散現(xiàn)象和時(shí)間延遲現(xiàn)象。

共振散射;散射時(shí)間序列;形態(tài)函數(shù);時(shí)頻分布

目標(biāo)散射一直是水聲學(xué)領(lǐng)域關(guān)注的熱點(diǎn)問(wèn)題,目標(biāo)散射問(wèn)題的主要研究方向有目標(biāo)探測(cè)、目標(biāo)特性分析以及目標(biāo)識(shí)別等。本文的研究對(duì)象是水中懸浮目標(biāo),其主要目的是通過(guò)實(shí)驗(yàn)對(duì)目標(biāo)散射的時(shí)域信號(hào)進(jìn)行分析。首先通過(guò)理論得出相關(guān)的解析解,并通過(guò)計(jì)算機(jī)進(jìn)行數(shù)值模擬,分析其所的到的結(jié)果,再以表面硬度較大的彈性目標(biāo)為實(shí)驗(yàn)對(duì)象,并與理論和數(shù)值結(jié)果進(jìn)行比較,最后得出相關(guān)結(jié)論。

關(guān)于聲波入射下的目標(biāo)散射問(wèn)題,國(guó)內(nèi)外都進(jìn)行了很多研究。Flax[1]首次將量子力學(xué)中的共振散射理論(RST)應(yīng)用到彈性體的散射,在理論上指出了共振現(xiàn)象的存在,他推導(dǎo)了圓柱和球體的形態(tài)函數(shù)公式,并將形態(tài)函數(shù)看成各階分波形態(tài)函數(shù)的疊加,其中每一階形態(tài)函數(shù)可以看成共振項(xiàng)和背景項(xiàng)的疊加。Veksler N D[2]和W hite R M[3]對(duì)不同參數(shù)的圓柱散射目標(biāo)進(jìn)行了理論分析和數(shù)值模擬。Goodman[4]分析了彈性球殼的共振散射情況,對(duì)球殼內(nèi)外流體一致時(shí)進(jìn)行了討論。Diarmuid J[5]在Goodman工作上進(jìn)一步討論了球殼內(nèi)外流體密度不同的時(shí),共振散射的情況,并給出了形態(tài)函數(shù)中共振項(xiàng)和背景項(xiàng)的表達(dá)式。Tesei[6]使用T矩陣法計(jì)算得出了散射信號(hào),并結(jié)合RST對(duì)散射信號(hào)進(jìn)行了分析,得到了很好的結(jié)果。Gaunaurd G C[7]與B reitenbach E D[8]分析了彈性圓柱殼的共振散射理論,并解出了形態(tài)函數(shù)的表達(dá)式以及各階形態(tài)函數(shù)的分波級(jí)數(shù),并對(duì)其進(jìn)行了分析。

Williams KL與Marston PL[9]使用Sommerfeld-Waston變換(SW T)分析了彈性球殼,并且進(jìn)行了實(shí)驗(yàn)驗(yàn)證。Sommerfeld-Waston變換是將分波形態(tài)轉(zhuǎn)化到復(fù)空間中的積分,通過(guò)求解積分得到分波級(jí)數(shù)。Marston[10]還提出幾何繞射理論,在ka不大的情況下,即相對(duì)目標(biāo)在低頻情況下可使用這種方法。Kargl[11]用該方法計(jì)算得到散射信號(hào)強(qiáng)度和散射信號(hào)各種成分的到達(dá)時(shí)間,并在實(shí)驗(yàn)室中比較了高頻聲脈沖入射小型號(hào)目標(biāo)的時(shí)間序列,從而驗(yàn)證了該理論對(duì)于目標(biāo)散射的正確性。時(shí)頻分析工具很多,如常用的如短時(shí)傅里葉變換,小波變換等,LatifR[12]使用了了維格納威利分布(Wigner-Ville)分布,以圓柱為目標(biāo),從而對(duì)散射信號(hào)進(jìn)行了分析,并且得到了彈性波的群速度求解方法,求解了群速度。劉孝斌[13]碩士學(xué)位論文通過(guò)數(shù)值模擬方法,從目標(biāo)散射的時(shí)域信號(hào)中得到了形態(tài)函數(shù)。本文對(duì)水中彈性實(shí)心圓柱進(jìn)行了聲散射實(shí)驗(yàn)研究(實(shí)驗(yàn)裝置詳見圖1),通過(guò)對(duì)彈性圓柱發(fā)射聲脈沖,得到其聲散射信號(hào)的時(shí)域波形,并從該信號(hào)提取出形態(tài)函數(shù),所得結(jié)果與共振散射理論所得的數(shù)值結(jié)果相符合。文章又對(duì)所得的時(shí)域信號(hào)進(jìn)行了時(shí)頻分析,得出了接收信號(hào)中存在頻散現(xiàn)象的結(jié)論。

1 彈性圓柱的共振散射理論

設(shè)平面波入射到圓柱上,入射聲壓為[1]

根據(jù)邊界條件[2]得到圓柱散射聲場(chǎng)的表達(dá)式為

ρ為水中水的密度,ρ1為圓柱的密度,為了簡(jiǎn)單起見,這里用引入函數(shù)

函數(shù)Sn稱做分波散射函數(shù),把δn稱做散射相位躍變

考慮到公式(4),圓柱的散射聲壓可以寫成下面的形式

在遠(yuǎn)場(chǎng)情況下krμ1

這樣在遠(yuǎn)場(chǎng)式(7)寫成下面形式

引入由分波形態(tài)函數(shù)fn(<)組成的形態(tài)函數(shù)f(<)

這樣式(8)可以表示成

分波形態(tài)函數(shù)fn(<)可以寫成

其中分波形態(tài)函數(shù)由硬共振項(xiàng)g(ns)(<)(或軟共振項(xiàng)和硬背景項(xiàng)(或軟背景項(xiàng))構(gòu)成[1]。

在彈性圓柱內(nèi)激發(fā)的彈性波相速度與共振頻率[2]關(guān)系式為:

n對(duì)應(yīng)共振模態(tài),l對(duì)應(yīng)共振族。群速度表達(dá)式為:

如圖1所示,入射波耦合到彈性圓柱體的角度[10]為

利用聲傳播的幾何關(guān)系可以求得鏡反射和圓柱體再輻射彈性波之間的時(shí)間間隔

其中,m=0,1,2……對(duì)應(yīng)于表面環(huán)繞波繞行圓柱的次數(shù),見圖1。

根據(jù)共振散射理論,散射信號(hào)可以表示為共振項(xiàng)和剛性背景項(xiàng)的疊加,在時(shí)間序列中背景項(xiàng)和共振項(xiàng)會(huì)有時(shí)間上的間隔,這樣就可以將背景項(xiàng)提取出來(lái),從而計(jì)算得到形態(tài)函數(shù)。

圖1 圓柱截面散射示意圖Fig.1 Cylindrical section scattering Fig.re

2 數(shù)值計(jì)算和分析

數(shù)值計(jì)算所用參數(shù)為:圓柱的密度8 900 kg/m3,球殼壓縮波速度為4 700 m/s,剪切波速度為2 322m/s。球殼內(nèi)外流體密度均為1 000 kg/m3,流體中聲速為1 450 m/s,半徑為0.035 m。

根據(jù)文章后續(xù)中實(shí)驗(yàn)的需要,這里選取20,36,49,60,71,81k作為入射頻率,由式(14)可知,這些頻率對(duì)應(yīng)著前6階的瑞利波,下圖列出了n=2～7階的共振項(xiàng)形態(tài)函數(shù),第1個(gè)峰值為瑞利波。在這些頻率的激發(fā)下,彈性波較強(qiáng)。

圖2 圓柱的階的共振項(xiàng)形態(tài)函數(shù)Fig.2 Cylindrical resonance form function(FF)n=2～7

根據(jù)式(17)可知,瑞利波的情況下,隨著n的增大,第1次到達(dá)時(shí)間是在不斷延長(zhǎng)的,這就是群速度的頻散現(xiàn)象和時(shí)間延遲。瑞利波的傳播速度在低頻的情況下較高頻的情況快,這就是所謂的頻散現(xiàn)象。下圖繪出了前6階瑞利波的到達(dá)時(shí)間圖,可以明顯的看出頻散現(xiàn)象的存在,低頻時(shí)到達(dá)時(shí)間較短,高頻到達(dá)時(shí)間較長(zhǎng)。

圖3 圓柱n=2～7階彈性波到達(dá)時(shí)間頻散圖Fig.3 The arrival time of elastic wave n=2～7

3 實(shí)驗(yàn)研究與計(jì)算

實(shí)驗(yàn)連接見圖4,發(fā)射換能器發(fā)射脈沖信號(hào),接收水聽器首先接收到直達(dá)聲信號(hào),再接收目標(biāo)散射信號(hào),接收水聽器將聲信號(hào)經(jīng)過(guò)測(cè)量放大器放大,再由信號(hào)采集器接收,輸入計(jì)算機(jī)。

由于實(shí)驗(yàn)所用的銅圓柱不是無(wú)限長(zhǎng)的,而且圓柱側(cè)面與底面之間不是平滑的,因此在棱上包裹了厚橡膠層,從而抑制棱散射波的影響。圓柱的中軸與聲波入射方向垂直,且與水面平行。在數(shù)值模擬中,假定圓柱為無(wú)限長(zhǎng),這里選用圓柱長(zhǎng)度為15 cm,大于所發(fā)射聲波波長(zhǎng)的3～5倍,因此該實(shí)驗(yàn)數(shù)據(jù)具有可信性,可與數(shù)值結(jié)果比較驗(yàn)證。

圖4 實(shí)驗(yàn)連接示意圖Fig.4 Experiment connection diagram

下圖為聲源發(fā)射信號(hào),以f=60k,,即ka=x51為例

圖5 f=60k發(fā)射信號(hào)時(shí)域與頻域波形圖Fig.5 f=60k Hz,the time domain and the frequency domain of the source signal

圖6 不同頻率下散射信號(hào)與理論時(shí)間對(duì)比Fig.6 The comparison between the scatter signal and the theoretic time,w hen the frequency is different

圖6對(duì)比了實(shí)驗(yàn)所測(cè)得的散射信號(hào),頻率為49k,60k,71k,81k(n=4～7階),與圖3所對(duì)應(yīng)的到達(dá)時(shí)間的對(duì)比結(jié)果。可以在時(shí)域上較為清晰地分辨出鏡反射波(第1條豎線)與再輻射波(第2、第3條豎線),尤其在式(17)m=0的情況下,即第一次到達(dá)再輻射波。而m=1時(shí),由于振幅的衰減較為嚴(yán)重,不容易在時(shí)域上進(jìn)行分辨,但在之后的時(shí)頻域分析中(文章第4部分),可以在偽彩圖上分辨。從圖上看,49k與60k的散射信號(hào)較為清晰,而較高頻率71k和81k的散射信號(hào),特別是彈性波部分干擾較為嚴(yán)重,其原因是由于較高的頻率能激發(fā)更多的共振頻率,如其它模態(tài)上的瑞利波以及回音廊波,從而在時(shí)域上受到干擾?？偟膩?lái)說(shuō),實(shí)驗(yàn)所得到的彈性波的時(shí)間位置與理論得到的彈性波位置基本吻合。

這里發(fā)射的信號(hào)是1個(gè)填充正弦的脈沖信號(hào),頻率帶寬不大,。因此,這里使用準(zhǔn)諧波法來(lái)繪制形態(tài)函數(shù),即發(fā)射多個(gè)不同頻率的脈沖信號(hào),每次得到1個(gè)區(qū)間的形態(tài)函數(shù)值,再進(jìn)行疊加,從而得到1個(gè)有一定寬度的形態(tài)函數(shù)。這里的發(fā)射頻率與前文所使用的6個(gè)頻率一致?？傻眯螒B(tài)函數(shù),見圖7。

圖7 實(shí)驗(yàn)所得形態(tài)函數(shù)(虛線)與理論所得形態(tài)函數(shù)比較(實(shí)線)Fig.7 The comparision between the experiment FF and theoretic FF

圖7的虛線為實(shí)驗(yàn)所得到的形態(tài)函數(shù),實(shí)線為理論所得到得形態(tài)函數(shù)。實(shí)驗(yàn)的形態(tài)函數(shù)通過(guò)疊加計(jì)算以后,在中間頻率部分與理論所得的形態(tài)函數(shù)無(wú)論在幅值和頻率上都有較好的吻合度,即在415的情況下,也就是在低頻和高頻段,以及發(fā)射頻帶的旁瓣覆蓋的區(qū)域形態(tài)函數(shù)都不理想,符合度較差?？傮w來(lái)說(shuō),在中頻部分所得到的形態(tài)函數(shù)較為理想,可以明顯的看到共振的峰值。

4 散射信號(hào)時(shí)頻分析

這里使用維格納威利分布(Wigner-Ville distribution)對(duì)散射信號(hào)進(jìn)行時(shí)頻分析。維格納威利分布定義式如下[14]它表示信號(hào)的瞬時(shí)自相關(guān)函數(shù)的傅里葉變換。在短時(shí)傅里葉變換中,信號(hào)只出現(xiàn)一次,為信號(hào)的線性表示;而在維格納威利分布中,信號(hào)在其計(jì)算中出現(xiàn)2次,所以維格納威利分布為信號(hào)的雙線性表示。

由于維格納威利分布進(jìn)行了自相關(guān),因此對(duì)于2個(gè)信號(hào)和的維格納威利分布等于每個(gè)信號(hào)的自維格納威利分布之和,再加上他們的交叉維格納威利分布的實(shí)部,且交叉項(xiàng)總是出現(xiàn)在2個(gè)維格納威利分布位置的中間。

圖8為f=49k時(shí),散射信號(hào)根據(jù)維格納威利分布(W igner-Ville distribution)公式(18)所得到的時(shí)頻聯(lián)合分布圖。

圖8 f=49k,散射信號(hào)維格納威利分布圖Fig.8 f=49k,WVD of the scattering signal

圖8的黑色橢圓標(biāo)記從左到右分別代表了鏡反射波,第1次環(huán)繞所輻射能量及第2次環(huán)繞所輻射能量,所標(biāo)記的具體時(shí)間與圖6一致?？梢钥吹?圖中的鏡反射波最為明顯,第1次環(huán)繞后所輻射的彈性波也可以分辨,而第2次環(huán)繞的彈性波則無(wú)法再此圖中辨別,這主要是因?yàn)樵诰S格納威利分布中存在干擾項(xiàng)(圖中箭頭標(biāo)出),鏡反射波與第一次輻射的彈性波中間的成分就是兩者的干擾項(xiàng)。這也就直接導(dǎo)致無(wú)法從維格納威利分布的時(shí)域坐標(biāo)中分辨散射信號(hào)成分,但卻得到較好的頻域分辨率。從圖8的中間橢圓信號(hào)可以明顯的看出信號(hào)在頻域上有峰值,這些峰值正好對(duì)應(yīng)著左圖總頻譜的峰值,且可以看出第1次輻射信號(hào)的底部明顯向右傾斜,這就是前面所提到的頻散現(xiàn)象,低頻部分群速度比高頻部分的群速度大。

為了降低干擾項(xiàng)的影響,這里引入偽維格納威利分布(Pseudo W igner-Ville distribution)

由(20)可知,偽維格納威利分布像短時(shí)傅里葉變換一樣,通過(guò)對(duì)信號(hào)加窗,有截?cái)嘈盘?hào)的維格納威利分布來(lái)分析信號(hào)的局部特性。這樣使得維格納威利分布的頻率分辨率降低,但對(duì)于多種分量的信號(hào),這種平滑可以降低交叉項(xiàng)的影響。

使用偽維格納威利分布對(duì)f=49k的散射信號(hào)進(jìn)行時(shí)頻分析。

圖9 f=49k,散射信號(hào)偽維格納威利分布圖Fig.9 f=49k,PWVD of the scattering signal

圖9可以看到,在偽維格納威利分布圖中,可以明顯的看到鏡反射波和2個(gè)彈性波的信號(hào),且在時(shí)域上有著良好的區(qū)分度,雖然3個(gè)信號(hào)中間中依然有干擾項(xiàng)的存在(圖中箭頭標(biāo)出),但已不影響在時(shí)域上分辨信號(hào)。在時(shí)域上得到了好的分辨率后,頻域上必然會(huì)變的模糊,可以看到2個(gè)彈性波信號(hào)在頻域上無(wú)法分辨各個(gè)共振峰值,但是依然可以看到這2個(gè)信號(hào)的底部是向右偏移的,頻散現(xiàn)象依然觀察到。

5 結(jié)語(yǔ)

文章通過(guò)實(shí)驗(yàn)得到了目標(biāo)的散射信號(hào)的時(shí)間序列,將時(shí)域波形分為了鏡反射波部分和彈性波部分,鏡反射部分對(duì)應(yīng)著形態(tài)函數(shù)中的背景項(xiàng),彈性波部分則為目標(biāo)的再輻射波形,對(duì)應(yīng)著形態(tài)函數(shù)的共振項(xiàng)。文章在時(shí)域上進(jìn)行了詳細(xì)的分析,所測(cè)得的信號(hào)的到達(dá)時(shí)間與理論推導(dǎo)基本吻合。在頻域上,通過(guò)準(zhǔn)諧波法,實(shí)驗(yàn)得到了目標(biāo)的形態(tài)函數(shù),在中頻部分吻合較好。最后又對(duì)實(shí)驗(yàn)的散射信號(hào)進(jìn)行了時(shí)頻分析,通過(guò)對(duì)幾種維格納威利分布的討論,在時(shí)域上得到了較好的分辨率,明確的觀察到了散射波的時(shí)域成分,頻域上也論證了共振峰值以及頻散現(xiàn)象的存在。

[1] Flax L,Dragonette L R,Uberall H.Theory of elastic resonance excitation by sound scattering[J].J Acoust Soc Am,1978,63(3):723-731.

[2] Veksler N D.Information Analysis in Hydroelasticity[M].Russian:Valgus Tallinn,1982.

[3] White R M.Elastic wave scattering at a cylindrical discontinuity in a solid[J].J Acoust Soc Am,1958,30:771-785.

[4] Ralph R.Goodman and raya stern,reflection and transmiss-ion of sound by elastic spherical shells[J].J Acoust Soc Am,1962,34(3):338-344.

[5] JDiarmuid Murphy,Jacob George,Anton Nagl,et al.Isolation of the resonant component in acoustic scattering from fluid-loaded elastic spherical shell[J].J Acoust Soc Am,1979,65(2):368 -373.

[6] Tesei A,Maguer W L,Fox J.Measurem-ents and modeling of acoustic scattering from partially and completely buried spherical shells[J].J Acoust Soc Am,2002,112(5):1817-1830.

[7] Guanaued G C,Brill D.Acoustic spectrogram and complex-frequency poles of a resonantly excited elastic tube[J].J Acoust Soc Am,1984,75:1680-1693.

[8] Breitenbach ED,Uberall H,Yoo KB.Resonant scattering from elastic cylindrical shells[J].J Acoust Soc Am,1983,74:1267-1273.

[9] Williams K L,Marston P L.Backscattering from an elastic sphere:Sommerfeld-Watson transformation and experimental confirmation[J].Acoustical Society of America Journal,1985,78:1093 1102.

[10] Marston P L.GTD for backscattering from elastic spheres and cylinders in water and the coupling of surface elastic waves with the acoustic field[J].Acoustical Society of America Journal,1988,83:25-37.

[11] Steven G,Kargl,Philip L.Observations and modeling of the backscattering of short tone bursts from a spherical shell:Lamb wave echoes,glory,and axial reverberations[J].The Journal of the Acoustical Society of America,1989,85(3):1014-1028.

[12] Latif R,Aassif E,Maze G,et al.Analysis of the circumferential acoustic waves backscattered by a tube using the time-frequency rep resentati on of Wigner-Ville[J].Measurement Science and Technology,2000,11:83-88.

[13] 劉孝斌.流體中彈性球殼聲散射研究[D].青島:中國(guó)海洋大學(xué),2009.

[14] 唐向宏,李齊良.時(shí)頻分析與小波變換[M].北京:科學(xué)出版社,2008.

Underwater Acoustic Scattering Study of Elastic Cylinder

HUA Da-Cheng,PENG Lin-Hui,YU Xiao-Tao
(College of Information Science and Engineering,Ocean University of China,Qingdao 266100,China)

In the paper,we performed acoustic scattering experiment of the elastic cylinder in water,combined with RST(Resonance Scattering Theory)to verify the experiment results.In the experiment,the time-domain scattering signals have been obtained,which coincide with the theoretical results.The aut her used the quasi-harmonic method to obtain the fo rm function,the results consistent with the theory.The aut her make use of Wigner-Ville distribution to analyze the time-frequency scattering relations,and obtain the Fig.res of the time-frequency distribution.which suggested that there is frequency dispersion phenomenon and the phenomenon of time delay.

resonance scattering;scattering time-series;fo rm function;time-frequency distribution

O427.2

A

1672-5174(2010)12-141-06

2010-03-03;

2010-06-02

華大成(1984-),男,碩士生。E-mail:huadacheng@sohu.com

責(zé)任編輯 陳呈超