彭茜蕤，周彥玲，范軍

?

水中規(guī)則波紋表面球回聲特性研究

彭茜蕤，周彥玲，范軍

(上海交通大學(xué)海洋工程國家重點實驗室，上海 200240)

運用聲吶目標(biāo)回聲特性預(yù)報的板塊元方法，研究了平面波入射水中規(guī)則波紋(正弦平方曲線波紋)表面球體的目標(biāo)回聲特性。分析了表面波紋高度、周期、入射平面波頻率和入射方位角對聲目標(biāo)強(qiáng)度空間分布和頻率響應(yīng)特性的影響規(guī)律。討論了不同頻段目標(biāo)強(qiáng)度空間分布特性的形成機(jī)理，得到了表面波紋參數(shù)與目標(biāo)強(qiáng)度頻率響應(yīng)和空間分布特性的近似定量關(guān)系。通過水槽實驗測量得到了3D打印的規(guī)則波紋表面尼龍球目標(biāo)強(qiáng)度的頻率響應(yīng)和空間分布特性，測量結(jié)果與理論計算結(jié)果吻合較好。掌握規(guī)則波紋表面球體回聲特性，有助于開展利用聲學(xué)超構(gòu)材料對散射聲場調(diào)控的研究。

波紋表面；回聲特性；目標(biāo)強(qiáng)度

0 引言

粗糙表面聲散射研究是一個既有理論意義又有實際意義的問題[1]。對粗糙表面回聲(反向散射)特性的研究可應(yīng)用于醫(yī)療超聲、超聲無損檢測、海底和海面混響建模與抑制、海洋浮游生物分類和探測等方面。從Lord Rayleigh開始，對粗糙表面聲散射已經(jīng)進(jìn)行了長期的研究。大多數(shù)水中粗糙表面聲散射的研究都以起伏海底、海面聲散射和混響建模為背景[2-4]，針對粗糙平面聲散射進(jìn)行研究，其中，粗糙表面假設(shè)為滿足高斯統(tǒng)計分布的隨機(jī)起伏平面或正弦起伏平面。研究的方法包括Kirchhoff近似方法、基于瑞利假設(shè)的微擾方法、積分方程方法以及小斜率近似方法等。部分針對水下粗糙曲面目標(biāo)聲散射開展研究，T. K. Stanton針對水下浮游生物聲散射，基于Sommerfeld-Watson變換和變形柱方法，研究了滿足高斯統(tǒng)計分布的隨機(jī)粗糙表面的拉長彈性柱體目標(biāo)平均散射聲場[5-7]，P. Jansson基于從Helmholtz 積分方程推導(dǎo)得到的零場方法，研究了滿足高斯統(tǒng)計分布的隨機(jī)粗糙表面球形目標(biāo)的聲散射特性[8]。

確定性規(guī)則波紋(正弦波紋、正弦平方波紋只是其中的子類)表面目標(biāo)是隨機(jī)粗糙表面目標(biāo)的一個特例，無論表面起伏幅度和斜率的大小程度如何改變，其散射聲場總是相干的，這就使人們可以更為深入地研究粗糙表面目標(biāo)的聲散射特征，如J. A. Fawcett[9]采用Kirchhoff/衍射混合方法研究了具有正弦平方波紋起伏的剛性球和剛性圓柱的寬帶脈沖信號聲散射時域特征，但針對具有規(guī)則波紋表面的水下曲面目標(biāo)聲散射特性的研究較少。近年來，以目標(biāo)表面微結(jié)構(gòu)設(shè)計為基礎(chǔ)的人工聲學(xué)超材料的研究發(fā)現(xiàn)，表面規(guī)則波紋(狹縫)可以產(chǎn)生許多新的物理現(xiàn)象[10-11]，如人工牛眼結(jié)構(gòu)就是平面表面上由多個周期性同心凹槽組成的，這也可以看成是一類具有規(guī)則波紋結(jié)構(gòu)的目標(biāo)，這類人工牛眼聲學(xué)結(jié)構(gòu)可以引起聲波異常透射和聲波能量的單向傳輸。掌握規(guī)則波紋表面目標(biāo)回聲特性有助于開展利用聲學(xué)超構(gòu)材料對散射聲場調(diào)控的研究。

本文應(yīng)用聲吶目標(biāo)回聲特性預(yù)報板塊元方法，研究了平面波入射到水中規(guī)則波紋(正弦平方曲線波紋)表面(剛性和非剛性)實心球體目標(biāo)的回聲特性；分析了表面波紋高度、周期、入射平面波頻率和入射方位角對聲目標(biāo)強(qiáng)度空間分布和頻率響應(yīng)特性的影響規(guī)律，討論了不同頻段目標(biāo)強(qiáng)度空間分布特性的形成機(jī)理；并開展了3D打印的規(guī)則波紋表面尼龍球回聲目標(biāo)強(qiáng)度的頻率響應(yīng)和空間分布特性的水聲水槽測量實驗。

1 規(guī)則波紋表面球體目標(biāo)

本文采用平面內(nèi)滿足式(1)的正弦平方波紋曲線繞軸旋轉(zhuǎn)得到的三維規(guī)則波紋表面球作為典型目標(biāo)開展研究，其表面方程見式(2)。

在計算仿真中假設(shè)表面具有以下兩種聲學(xué)特性[12-13]，假設(shè)1是剛性表面，表面滿足式(4)的剛性邊界條件，假設(shè)2是非剛性表面，表面滿足式(5)的邊界條件：

2 回聲特性預(yù)報的板塊元方法

基于Kirchhoff 近似所建立的板塊元方法是一種相對較為成熟的對聲吶目標(biāo)中高頻回聲特性進(jìn)行預(yù)報的數(shù)值方法[12]。對于收-發(fā)合置的情況，基于物理聲學(xué)或Kirchhoff 近似的非剛性聲吶目標(biāo)的遠(yuǎn)場目標(biāo)強(qiáng)度(TS)為[13]

板塊元方法首先將目標(biāo)的三維幾何表面進(jìn)行空間板塊離散劃分，用一組三維平面板塊來近似復(fù)雜形狀的目標(biāo)曲面，一般采用三角形板塊進(jìn)行劃分。所有板塊的散射聲場之和就是總散射聲場的近似值。算法中將單個三維板塊變化到局部二維平面，式(10)中的面積分則轉(zhuǎn)化為與板塊頂點坐標(biāo)相關(guān)的代數(shù)運算，如式(11)所示，從而避免了面積分運算，提高了計算速度。

假如劃分為個平面板塊，則式(10)可以寫成：

板塊元方法比通常的面元積分法計算速度快，不僅適用于遠(yuǎn)場、剛性目標(biāo)的回聲特性預(yù)報，還可推廣應(yīng)用于近場和非剛性表面的回聲特性預(yù)報，目前板塊元方法已經(jīng)較為廣泛應(yīng)用于聲吶目標(biāo)中高頻回聲特性預(yù)報[15-18]。

本文將采用板塊元方法對規(guī)則波紋表面球回聲特性進(jìn)行仿真計算，需要說明的是目前考慮的波紋球波紋起伏相對較小，計算中忽略了聲波在波紋凹槽之間的多次散射[19-20]。

3 數(shù)值計算與分析

以正弦平方波紋球為例，采用板塊元方法對其回聲特性進(jìn)行數(shù)值仿真。波紋球的半徑=0.3 m，在仿真中假設(shè)球面為剛性邊界條件。

3.1 目標(biāo)強(qiáng)度的頻率響應(yīng)特性

圖2 不同ξ條件下的波紋球目標(biāo)強(qiáng)度頻率響應(yīng)(a) φ=0°；(b) φ=90°

3.2 目標(biāo)強(qiáng)度的空間分布特性

圖3 波紋球目標(biāo)強(qiáng)度的空間分布特性

由圖3可以看出，在低頻(2 kHz)處波紋球目標(biāo)強(qiáng)度的空間分布基本是均勻的，如圖中虛線所示，接近光滑球的目標(biāo)強(qiáng)度空間分布，即波紋效應(yīng)沒有體現(xiàn)。隨著頻率的升高，波紋球目標(biāo)強(qiáng)度隨入射角變化呈現(xiàn)“花瓣狀”峰谷空間分布，并且頻率越高“花瓣狀”峰谷結(jié)構(gòu)分支越密集。相對來說，=8的波紋球目標(biāo)強(qiáng)度呈現(xiàn)的“花瓣狀”峰谷分支結(jié)構(gòu)比=5的更為密集。

值得注意的是，當(dāng)頻率升高到波紋效應(yīng)呈現(xiàn)的頻段，目標(biāo)強(qiáng)度“花瓣狀”峰谷空間分布的峰谷與波紋球表面波紋的峰/凹槽結(jié)構(gòu)基本對應(yīng)。并且隨著頻率升高到不同的頻段，目標(biāo)強(qiáng)度“花瓣狀”分支峰谷結(jié)構(gòu)對應(yīng)表面波紋結(jié)構(gòu)的不同部分。圖4為=8的波紋球局部放大圖。

圖4 表面波紋局部結(jié)構(gòu)

按如下原則得到頻段劃分近似公式：

圖5 頻率-角度譜

因正弦平方波紋球具有局部多鏡面反射(波峰、谷凹槽、凹槽斜線)，使得其散射聲場在空間形成多峰結(jié)構(gòu)，這種多峰結(jié)構(gòu)與入射波頻率、入射角度緊密相關(guān)，可以通過調(diào)整表面波紋周期、起伏幅度對散射聲場的空間分布特性進(jìn)行調(diào)控。

4 聲散射實驗

4.1 實驗布局

圖6 3D打印的波紋實驗?zāi)Ｐ秃蛯嶒灢贾脠D

4.2 實驗數(shù)據(jù)分析

圖7 目標(biāo)強(qiáng)度的實驗與數(shù)值計算結(jié)果對比

5 結(jié)論

通過板塊元數(shù)值計算和實驗，研究了平面波入射水中規(guī)則波紋(正弦平方曲線波紋)表面球體的目標(biāo)回聲特性。結(jié)果表明：低頻時波紋高度、周期對聲目標(biāo)強(qiáng)度基本沒有影響，隨著頻率的升高，多鏡面反射(波峰、波谷凹槽和凹槽斜面)的影響清晰顯示出來，且波紋高度越高和周期越大，聲目標(biāo)強(qiáng)度起伏越大、峰值越密集。通過水槽實驗測量得到了3D打印的規(guī)則波紋表面尼龍球目標(biāo)強(qiáng)度的頻率響應(yīng)和空間分布特性，測量結(jié)果與數(shù)值計算結(jié)果吻合較好。規(guī)則波紋表面球體回聲特性的掌握有助于開展利用聲學(xué)超構(gòu)材料對散射聲場調(diào)控的研究。

[1] OGILVY J A. Theory of wave scattering from random rough surfaces[M]. Bristol: Hilger, 1991.

[2] JACKSON D R, RICHARDSON M D. High-frequency seafloor acoustics[M]. New York: Springer, 2007.

[3] WELTON P J. Three dimensional analyses of scattering by pressure release sinusoidal surfaces[J]. J. Acoust. Soc. Am., 2012, 131 (1): 54-68.

[4] CHOO Y M, SONG H C, SEONG W J. Numerical study of three-dimensional sound reflection from corrugated surface waves [J]. J. Acoust. Soc. Am., 2016, 140(4): 2290-2296.

[5] STANTON T K. Sound scattering by rough elongated elastic objects. I: Means of scattered field[J]. J. Acoust. Soc. Am., 1992, 92(3): 1641-1664.

[6] STANTON T K, CHU D. Sound scattering by rough elongated elastic objects. II: Fluctuations of scattered field[J]. J. Acoust. Soc. Am., 1992, 92(3): 1665-1678.

[7] GURLEY J V, STANTON T K. Sound scattering by rough elongated elastic objects. III. Experiment[J]. J. Acoust. Soc. Am., 1993, 94(5): 2746-2755.

[8] JANSSON P. Acoustic scattering from a rough sphere[J]. J. Acoust. Soc. Am., 1993, 93 (6): 3032–3042.

[9] FAWCETT J A. Modeling of high-frequency scattering from objects using a hybrid Kirchhoff/diffraction[J]. J. Acoust. Soc. Am., 2001, 109(4): 1312-1319.

[10] MEI J, HOU B, KE M Z, et al. Acoustic wave transmission through a bull’s eye structure[J]. Appl Phys Lett, 2008, 92: 1-3.

[11] WANG J W, YUAN B G, CHENG Y, et al. Unidirectional acoustic transmission in asymmetric bull’s eye structure[J]. 2015, 58(2): 1-5.

[12] 范軍, 湯渭霖, 卓琳凱. 聲吶目標(biāo)回波特性預(yù)報的板塊元方法[J].船舶力學(xué), 2012, 16(1-2): 171-180.

FAN Jun, TANG Weilin, ZHUO Linkai. Planar elements method for forecasting the echo characteristics from sonar targets[J]. Journal of Ship Mechanics, 2012, 16(1-2): 171-180.

[13] 湯渭霖. 用物理聲學(xué)方法計算非硬表面的聲散射[J]. 聲學(xué)學(xué)報, 1993, 18(1): 45-53.

TANG Weilin. Calculation of acoustic scattering of a nonrigid sur- face using physical acoustic method[J]. Acta Acustica, 1993, 18(1): 45-53.

[14] 布列霍夫斯基赫著. 分層介質(zhì)中的波[M]. 楊訓(xùn)仁譯. 北京: 科學(xué)出版社: 1960.

BREKHOVSKIKH L. Wave in layered medium[M]. YANG Xunren, translate. Beijing: Science Press: 1960.

[15] 鄧偉, 范軍. 基于魚雷導(dǎo)引彈道的目標(biāo)回聲特性研究[J]. 聲學(xué)技術(shù), 2014, 33(5): 402-407.

DENG Wei, FAN Jun. Research on target echo characteristics in torpedo guidance trajectory[J]. Technical Acoustics, 2014, 33(5): 402-407.

[16] 卜考英, 范軍. 基于多波束系統(tǒng)的水下目標(biāo)近程回波亮點聲圖像[J].聲學(xué)技術(shù), 2007, 26(4): 570-573.

BU Kaoyin, FAN Jun. Acoustic image of bright spots of near-field underwater targets based on multi-beam system[J]. Technical Acoustics, 2007, 26(4): 570-573.

[17] 王二慶, 卓琳凱, 王華奎. 單雙層殼體潛艇的高頻回波特性研究[J].聲學(xué)技術(shù), 2012, 31(5): 486-489.

WANG Erqing, ZHUO Linkai, WANG Huakui. High-frequency echo characteristics of single and double shell submarine[J]. Technical Acoustics, 2012, 31(5): 486-489.

[18] 王汝夯, 李冰, 魏強(qiáng). 基于板塊元的水面艦船遠(yuǎn)場聲目標(biāo)強(qiáng)度預(yù)報方法[J].聲學(xué)技術(shù), 2015, 34(1): 6-10.

WANG Ruhang, LI Bing, WEI Qiang. A method for predicting far-field target strength of surface ship based on planar elements theory[J]. Technical Acoustics, 2015, 34(1): 6-10.

[19] 鄭國垠, 范軍, 湯渭霖. 考慮遮擋和二次散射的修正板塊元算法[J]. 聲學(xué)學(xué)報, 2011, 36(4): 377-383.

ZHENG Guoyin, FAN Jun, TANG Weilin. A modified planar elements method considering occlusion and secondary scattering[J]. Acta Acustica, 2011, 36(4): 377-383.

[20] 張陽, 李桂娟, 王振山, 等. 改進(jìn)的板塊元遮擋算法[J]. 聲學(xué)技術(shù), 2016, 35(4): 319-324.

ZHANG Yang, LI Guijuan, WANG Zhenshan, et al. An improved plate element sheltering algorithm[J]. Technical Acoustics, 2016, 35(4): 319-324.

[21] 劉伯勝, 雷家煜. 水聲學(xué)原理[M]. 哈爾濱: 哈爾濱工程大學(xué)出版社, 1993.

LIU Bosheng, LEI Jiayu. Principles of underwater acoustics[M]. Harbin:Harbin Engineering University Press, 1993.

Echo characteristics of a regular corrugated surface ball in water

PENG Xi-rui, ZHOU Yan-ling, FAN Jun

(State Key Laboratory of Ocean Engineering, Shanghai Jiao Tong University, Shanghai 200240, China)

Echo characteristics of a regular corrugated (sine squared curve) surface sphere in water with an incident plane wave are investigated by using the planar elements method which is used for forecasting echo characteristics of sonar targets. The effects of surface corrugation’s height, period, incident plane wave frequency and incident azimuth angle on the spatial distribution and frequency response characteristics of acoustic target strength are analyzed. The formation mechanism of the spatial distribution characteristics of target strength in different frequency bands is discussed. The approximate quantitative relationship between the parameters of corrugated surface and the spatial distribution as well as the frequency response characteristics of target strength is obtained. The target strength of the 3D printed regular corrugated surface nylon sphere is measured in acoustic tank. The measured spatial distribution and frequency response characteristics of target strength are in good agreement with the theoretical results. A thorough understanding of the echo characteristics of regular corrugated surface sphere is helpful for conducting research on controlling the scattering sound field by using the acoustic metamaterials.

corrugated surface;echo characteristics; target strength

TB556

A

1000-3630(2018)-06-0528-07

10.16300/j.cnki.1000-3630.2018.06.003

2017-12-13;

2018-03-15

國家自然科學(xué)基金(11774229)資助項目。

彭茜蕤(1993－), 女, 湖北武漢人, 碩士, 研究方向為水下目標(biāo)聲散射。

范軍,E-mail: fanjun@sjtu.edu.cn