拋物面鏡球面波反射聲場特性的理論研究

張　軍，陳　鵬，陳正武，趙　云，曾新吾

(1.中國空氣動力研究與發(fā)展中心 氣動噪聲控制重點實驗室，四川 綿陽　621000；2.國防科學(xué)技術(shù)大學(xué) 光電科學(xué)與工程學(xué)院，長沙　410073)

?

拋物面鏡球面波反射聲場特性的理論研究

張軍1，陳鵬1，陳正武1，趙云2，曾新吾2

(1.中國空氣動力研究與發(fā)展中心 氣動噪聲控制重點實驗室，四川 綿陽621000；2.國防科學(xué)技術(shù)大學(xué) 光電科學(xué)與工程學(xué)院，長沙410073)

根據(jù)Kirchhoff-Helmholtz聲衍射理論，獨立推導(dǎo)了沿拋物面鏡軸線的球面波反射聲場的時域理論解?；诶碚摻忸A(yù)測近場反射波中存在三種子波，即“中心波”、“邊緣波”和“尾波”，中心波的相位與邊緣波和尾波相反。遠(yuǎn)場反射波與聲源波形的導(dǎo)數(shù)形式相同，聲壓幅值和傳播距離成反比而與聲波頻率成正比。以典型的正弦波為例給出了反射聲場的演化形成過程，并通過COMSOL軟件進(jìn)行數(shù)值模擬驗證了三種子波的存在及理論解的正確性。最后研究了球面波反射聲場的特性，如果拋物面鏡的口徑固定，則存在一個最優(yōu)的深焦比參數(shù)d/zF=3.92使得遠(yuǎn)場的反射波聲功率密度最大。

拋物面鏡；球面聲波；反射聲場

拋物面鏡具有特殊的幾何性質(zhì)，平行射線經(jīng)拋物面反射后將會聚焦或者從焦點發(fā)出的射線經(jīng)拋物面反射后將成為平行射線。因此，拋物面鏡多被用作信號的接收器或發(fā)射器，例如在工程上常使用拋物面天線來接收電磁波信號以提高信噪比，或者將拋物面鏡作為高能激光的中繼鏡以提高機(jī)載激光器的傳輸距離[1]。在聲學(xué)領(lǐng)域，拋物面鏡本質(zhì)上屬于一種曲面聲障板，常規(guī)換能器使用聲障板可以提高聲波的發(fā)射效率。與相控陣技術(shù)相比，聲障板的方法易于實現(xiàn)且成本更低。張軍等[2]研究了電弧放電式等離子體聲源(Plasma Sound Source)的聚焦原理，指出橢球面鏡雖可以將球面波由近焦點聚焦到遠(yuǎn)焦區(qū)域，但是在過遠(yuǎn)焦點之后聲能量的衰減非常迅速。所以橢球面鏡聚焦系統(tǒng)適合用于治療人體結(jié)石等對傳輸距離要求不高的場合，而要在較遠(yuǎn)的距離上提高聲學(xué)換能器的發(fā)射效率則需要使用拋物面鏡。

王鴻樟[3]使用拋物面鏡作為聲反射障板提高了聚焦換能器的發(fā)射效率，他利用頻域Kirchhoff公式推導(dǎo)了一個發(fā)射聲壓的表達(dá)式。WAHLSTROM[4]使用一個傳聲器和拋物面鏡，從遠(yuǎn)處成功接收到了生物類聲源的聲信號。他在平面波入射條件下建立了拋物面鏡的線性聲波反射模型，推導(dǎo)了反射聲場的頻域解析解，并據(jù)此研究了反射聲場的特性。DAI等[5]使用拋物面鏡作為接收器提高了S1型-零群速度(S1ZGV)蘭姆波的信噪比，S1ZGV蘭姆波在薄板中很容易被激發(fā)，因此其常被用于超聲無損檢測。后來TSAI等[6-7]推導(dǎo)了平面波及軸對稱非平面波入射條件下拋物面鏡反射聲場的時域理論解，討論了拋物面鏡幾何參數(shù)及入射波頻率參數(shù)變化對反射聲場的影響。王鴻樟和Wahlstrom得到的都是頻域結(jié)果，而Tsai等人得到的是時域結(jié)果。與頻域結(jié)果相比，時域結(jié)果具有以下優(yōu)點：一是可以更清楚的揭示聲場形成過程；二是受聲波的頻帶范圍限制小。本文根據(jù)Kirchhoff-Helmholtz聲衍射理論，獨立推導(dǎo)了拋物面鏡球面波反射聲場的時域理論解，建立了將拋物面鏡作為聲發(fā)射器的時域理論方法，從而完善了王鴻樟等的理論研究。

本文的內(nèi)容分為四個部分：第一節(jié)建立了拋物面鏡球面波反射的理論模型，并給出了反射聲場沿拋物面鏡軸線的解析解和遠(yuǎn)場解；第二節(jié)以正弦波為例給出了沿拋物面鏡軸線的反射聲場的形成過程；第三節(jié)利用多物理場有限元軟件COMSOL進(jìn)行了數(shù)值模擬，驗證了理論結(jié)果的正確性；第四節(jié)研究了拋物面鏡反射聲場的特性，討論了開角和深焦比等參數(shù)對反射聲場的影響。

1　理論

從拋物面鏡焦點發(fā)出的球面波在傳播過程中將分成兩個部分，一部分仍按球面波形式擴(kuò)散，另一部分入射到壁面后發(fā)生反射。根據(jù)惠更斯原理，反射波可以視為新的次級波源，次級波相干疊加便形成了反射聲場。因此空間中的聲場包含直達(dá)波和反射波。本文主要研究球面波的反射聲場。如圖1所示，從焦點發(fā)出的射線經(jīng)拋物面鏡反射后將成為平行于其軸線的射線。

圖1　拋物面鏡球面波反射的幾何示意圖Fig.1　The reflection of a spherical wave from a parabolic mirror

在直角坐標(biāo)系下，旋轉(zhuǎn)拋物面的方程可以表示為，

(1)

式中：(x,y,z)中的下標(biāo)s代表“Surface”，表示位于拋物面鏡上的點。R’代表x坐標(biāo)或y坐標(biāo)，zF表示拋物面鏡頂點到焦點之間的距離，即焦距；d表示拋物面鏡頂點到開口的距離，即深度；S為拋物面鏡上的一點，r1s=zs+zF為焦點到S之間的距離；r2s=d-zs為S到拋物面鏡開口之間的固定距離；r2為拋物面鏡開口到 r2s上任意一點之間的變化的距離。R為點S與拋物面鏡軸線上一點P之間的距離，n表示點S處垂直于拋物面鏡表面的內(nèi)法線方向單位向量。

假設(shè)傳播介質(zhì)無耗散、各向均勻，Gold’berg數(shù)(即非線性系數(shù)與耗散系數(shù)之比)小于1，則聲波的傳播由線性波動方程控制。

(2)

式中，c0表示介質(zhì)中的聲速。

根據(jù)封閉空間中的能流守恒原理，若已知封閉面上的聲場的分布則可得到封閉面內(nèi)任意一點處的聲場[8]。忽略封閉面上二次以上的聲反射，由式(2)可以得到非單色波的Kirchhoff-Helmholtz時域一般解。

(3)

Kirchhoff-Helmholtz公式是基于衍射理論的一般解，對具有波動特性的波(如光波、聲波)的傳播都適用。為了便于書寫，式(3)省略了物理量的時間變量和空間變量。其中，p表示目標(biāo)點處的反射聲壓，pg表示通過射線聲學(xué)得到的封閉面上的聲壓。?/?n表示沿拋物面內(nèi)法線方向的方向?qū)?shù)，S表示由拋物面鏡和從拋物面鏡出口到無窮遠(yuǎn)處的假想曲面組成的封閉面?！癧·]”是一種記號，表示其中的物理量經(jīng)歷了R/c0的時間延遲。根據(jù)Somerfield輻射條件[9]，可以忽略從無窮遠(yuǎn)處來的聲場的貢獻(xiàn)，式(3)的積分只在拋物面上進(jìn)行。

2.1邊界條件

假設(shè)聲波的波長與拋物面鏡的曲率半徑相比很小(ka?1，k表示波數(shù)，a表示曲率半徑)，則射線聲學(xué)的理論可用于聲場計算。球面波聲壓與傳播距離r之間有如下關(guān)系

(4)

(5)

(6)

在鏡面上一點S處，可以得到如下的邊界條件，

(7)

(8)

式中：ρ=0表示拋物面鏡開口圓的中心點，ρ=1表示邊緣點?？梢园l(fā)現(xiàn)，固定焦距，隨著拋物面鏡深度d增加反射波的壓力幅值減小。因為深度d增加將使得從焦點zF處發(fā)出的球面波到達(dá)相同鏡面點所經(jīng)歷的傳播距離延長，而球面波的聲壓幅值與傳播距離成反比。

2.2沿拋物面鏡軸線球面波反射聲場的理論解

邊界條件式(7)包含了任意形式的源函數(shù)f(t)，f(t)在積分過程中難以處理。根據(jù)線性系統(tǒng)理論，系統(tǒng)的輸出可以通過輸入函數(shù)與系統(tǒng)的沖擊響應(yīng)函數(shù)相卷積得到?？梢岳脹_擊函數(shù)δ(t)代替f(t)先求解系統(tǒng)的沖擊響應(yīng)函數(shù)h(t)，然后通過h(t)與f(t)相卷積便可以得到反射聲場的解。邊界條件(7)相應(yīng)改寫為，

(9)

結(jié)合拋物面鏡方程式(1)，經(jīng)過繁雜的數(shù)學(xué)推導(dǎo)得到了沖擊響應(yīng)函數(shù)。

(10)

式(10)右端各項的表達(dá)式如下

(11)

進(jìn)一步通過卷積便可以得到沿拋物面鏡軸線任意位置處的反射聲場的解。

(12)

式(12)右端包含三項，分別對應(yīng)從拋物面鏡頂點、出口邊緣和其余表面發(fā)出的反射波。根據(jù)Blackstocks[10]及Hamilton[7]，本文將其稱為“中心波”、“邊緣波”和“尾波”。從聲波的相位來看，邊緣波和尾波的相位與中心波是相反的。

2.3遠(yuǎn)場解

隨著觀察點遠(yuǎn)離拋物面鏡，中心反射波和邊緣反射波的到達(dá)時間差縮短，于是尾波項的積分時間趨近于0，尾波趨近于δ脈沖。在遠(yuǎn)場(z/zF?1)，式(12)退化為更加簡單的形式，

(13)

根據(jù)式(13)，遠(yuǎn)場反射波將退化為聲源波形的導(dǎo)數(shù)形式，例如從拋物面鏡焦點處發(fā)出的正弦波在遠(yuǎn)場將退化為余弦波。當(dāng)拋物面鏡深度遠(yuǎn)小于焦距時，上式與圓形平面活塞聲源的遠(yuǎn)場解形式上一致[11]。不難發(fā)現(xiàn)，拋物面鏡遠(yuǎn)場反射聲壓主要受幾何參數(shù)(d/zF)及聲源波形參數(shù)(df/dt)控制。定義遠(yuǎn)場反射波的聲壓幅值與p0之比為拋物面鏡的遠(yuǎn)場增益，當(dāng)焦距固定而深度變化時有如下關(guān)系

(14)

當(dāng)拋物面鏡的口徑固定而深焦比變化時，

(15)

對于簡諧變化的聲波，df(t)/dt～jω，則有

(16)

式中:λ=2πc0/ω為聲波的波長，ω為聲波的角頻率。可以發(fā)現(xiàn)，遠(yuǎn)場增益G3與聲波的波長λ成反比，但由于實際介質(zhì)中存在黏性和熱傳導(dǎo)耗散，高頻聲波更容易被吸收。

2　反射聲場的形成過程

下面以正弦波為例研究反射聲場的形成過程。取拋物面鏡的焦距zF=0.05 m，深度d=0.2 m；拋物面鏡的口徑ra=0.2 m。將球面波聲源置于拋物面鏡的焦點處。距離焦點r=zF處的球面波壓力幅值為P0，單個周期正弦波脈沖的持續(xù)時間T0=0.1 ms；傳播介質(zhì)為空氣，密度ρ0=1.25 kg/m3，聲速c0=340 m/s。計算結(jié)果如圖2所示。

圖2(a)為輸入的正弦波，圖2(b)～2(f)為沿拋物面鏡軸線的反射波計算結(jié)果。根據(jù)反射波形特點，可以將反射聲場分為近場和遠(yuǎn)場兩個區(qū)域：當(dāng)反射波形與初始波形的導(dǎo)數(shù)形式相同時表明觀察點到達(dá)遠(yuǎn)場區(qū)。當(dāng)拋物面鏡的深焦比d/zF較小時，近場和遠(yuǎn)場的臨界距離為平面活塞聲源的瑞利距離(lR=πa2/λ)[11]。

圖2　拋物面鏡軸線上球面波反射聲場的形成過程Fig.2　The evolution ofan on-axis sinusoidal waveform

在近場區(qū)，反射波形包含中心波“C”、邊緣波“E”和尾波“W”等復(fù)雜的波形結(jié)構(gòu)(見圖2(b)～2(c))。在到達(dá)觀察點的時間順序上，邊緣波滯后于中心波，“尾波”夾在中心波和邊緣波之間。根據(jù)幾何關(guān)系，當(dāng)深焦比d/zF較大時，邊緣波的傳播距離大于中心波的傳播距離，因此其達(dá)到時間也越晚。隨著距離增加，中心波和邊緣波到達(dá)觀察點的時間差將縮短。在幅值關(guān)系上，中心波比邊緣波大，兩者之間的關(guān)系受拋物面鏡的深焦比參數(shù)控制。深焦比越大，邊緣波的幅值越小。

在遠(yuǎn)場區(qū)，中心波、邊緣波和尾波將疊加在一起，反射波演化成聲源波形的導(dǎo)數(shù)形式。根據(jù)式(13)，遠(yuǎn)場反射波的聲壓幅值與傳播距離成反比。反射波的聲壓幅值隨距離增加下降較快。這在反射聲束上體現(xiàn)了一種“近場集中，遠(yuǎn)場發(fā)散”的特性。顯然，這是由聲波的衍射效應(yīng)引起的，頻率越低聲衍射效應(yīng)越明顯。這種性質(zhì)在具有波動性的光波和電磁波中也存在，例如手電筒發(fā)出的準(zhǔn)直光束到了較遠(yuǎn)距離處就逐漸發(fā)散了。盡管如此，在100倍焦距處，本算例的反射波聚焦增益超過了20 dB。

3　數(shù)值模擬驗證

王鴻樟等人的研究結(jié)果無法與本文的研究結(jié)果直接進(jìn)行對比，本節(jié)使用多物理場有限元軟件COMSOL對拋物面鏡的聲波反射過程進(jìn)行了模擬，希望通過模擬結(jié)果對推導(dǎo)結(jié)果進(jìn)行數(shù)值驗證。

COMSOL軟件包含了聲學(xué)、流體力學(xué)、結(jié)構(gòu)力學(xué)等多個求解模塊，該軟件的優(yōu)點是可以有效求解多個物理場的耦合問題。COMSOL的求解包括物理場選擇、幾何建模、網(wǎng)格劃分、求解器配置和結(jié)果后處理等幾個步驟。COMSOL軟件已由最初的3.X版本發(fā)展到了目前的5.X版本，本文使用的是5.0版本。首先選擇瞬態(tài)壓力聲學(xué)模型，并建立起簡化的二維軸對稱幾何模型，如圖3所示。

圖3　COMSOL軟件中建立的幾何模型和聲源輸入Fig.3　The geometric model of a parabolic mirror in COMSOL

幾何模型由參數(shù)化曲線(y2=4zFx)和貝氏曲線轉(zhuǎn)化成的實體區(qū)域組成。區(qū)域內(nèi)的材料為空氣，材料參數(shù)使用軟件內(nèi)置參數(shù)，區(qū)域的邊界采用剛性邊界條件，區(qū)域的總體尺寸如圖3所示。網(wǎng)格劃分采用自由三角形網(wǎng)格，最大網(wǎng)格尺寸hmax為聲波波長的1/6。最大網(wǎng)格尺寸和時間步長之間滿足CFL穩(wěn)定性條件，即tstep=CFL*hmax/Cair，CFL=0.05。聲源設(shè)置在拋物面鏡的焦點上，聲源波形為單個周期的正弦波f=Asin(f0*t)*(t

圖4給出了單個周期的正弦波在四個不同時刻的二維傳播云圖。

圖4　不同傳播時刻正弦波的二維云圖Fig.4　Snapshots of a sinusoidal wave at different times

圖4(a)，在t1=0.15 ms時刻從焦點發(fā)出的球面波的波陣面正好達(dá)到拋物面鏡的頂點。圖4(b)，t2=0.45 ms時刻聲場分成了三個部分：以球面波形式向前傳播的直達(dá)波“D”，反射中心波“C”，以及鏡面附近的反射尾波“W”。中心波的聲壓幅值高于直達(dá)波和尾波?？梢园l(fā)現(xiàn)，球面波經(jīng)拋物面鏡反射后變成了準(zhǔn)直平面波。圖4(c)，在t3=0.9 ms時刻反射中心波之后出現(xiàn)了邊緣波E。邊緣波起源于拋物面鏡的開口邊緣，在直達(dá)波到達(dá)鏡面邊緣的時刻產(chǎn)生，并以邊緣點為中心按球面波的形式擴(kuò)散。在向拋物面鏡軸線方向的傳播匯聚過程中，邊緣波的相位和中心波相反。尾波的幅值較小，云圖中難以看出尾波的存在。圖4(d)，在t4=1.05 ms時刻可以清楚的看出中心波和邊緣波的傳播情況。根據(jù)模擬結(jié)果，理論預(yù)測的近場三種子波確實是存在的。

假設(shè)單個周期正弦波的持續(xù)時間T=0.1 ms，即聲波的中心頻率f0=10 kHz。按照文獻(xiàn)[6]中的計算方法得到ka=19，表明聲波在鏡面附近的反射可以使用射線理論求解。沿拋物面鏡軸線z/zF=5，10，50，100將理論解和軟件的模擬結(jié)果進(jìn)行對比，如圖5所示。

圖5　COMSOL軟件模擬結(jié)果與理論解的對比Fig.5　Results of the present solution and COMSOL simulation

圖5中藍(lán)色實線代表理論解，紅色虛線代表使用COMSOL軟件得到的模擬結(jié)果?？梢园l(fā)現(xiàn)，理論解的中心波和COMSOL的模擬結(jié)果之間符合得很好，而理論解的邊緣波幅值比模擬結(jié)果略小。兩者之間的差異可能是由網(wǎng)格剖分在拋物面鏡邊緣處不細(xì)致引起的。對有限元軟件而言，網(wǎng)格剖分的好壞將直接影響計算結(jié)果。對于遠(yuǎn)場反射波，理論解和COMSOL模擬結(jié)果之間符合得非常好。從而通過數(shù)值模擬方法使第二節(jié)的理論解的正確性得到了驗證。

4　拋物面鏡球面波反射聲場的特性

反射鏡的形狀主要由深焦比d/zF或開角θ等參數(shù)控制。θ是由焦點與拋物面鏡出口邊緣的連線與軸線之間組成的夾角，

(16)

當(dāng)θ=90°時，拋物面鏡的焦點正好位于出口面上。如果焦距固定而深焦比增加，則拋物面鏡的口徑和開角也增加；如果口徑固定，則隨深度變小開角也變小。

分別將拋物面鏡的焦距和口徑固定，遠(yuǎn)場增益隨深焦比的變化如圖6所示。固定拋物面鏡焦距，遠(yuǎn)場增益將隨深度的增加而緩慢增加，這是因為深度增加將使得更多的聲線受到拋物面鏡的反射。固定拋物面鏡口徑，則存在一個最優(yōu)的深焦比參數(shù)d/zF=3.92使得遠(yuǎn)場增益達(dá)到最大值。Wahlstrom在平面波入射條件下也得到了一個相同的深焦比參數(shù)，使得拋物面鏡焦點處的反射聲壓達(dá)到最大值[3]。不難理解，這正是線性條件下聲波的傳播滿足互易原理的體現(xiàn)。

圖6　遠(yuǎn)場增益與深焦比的關(guān)系Fig.6　The relationship between depth-to-focal-length ratio and far field gain

取焦距zF=0.05 m，單個周期正弦波的持續(xù)時間分別為T1=0.1 ms和T2=0.02 ms，拋物面鏡的開角θ分別為60°，90°，120°和150°。根據(jù)式(12)和(17)，得到反射波的聲功率密度的變化，如圖7所示。

圖7　不同開角，拋物面鏡軸線上聲功率密度的變化Fig.7　The relationship between on-axis sound power density and axial propagation distance with different rim angles

由圖7可見，反射波的聲功率密度隨傳播距離呈現(xiàn)先增加而后減小的趨勢，這是因為反射波形的演化使得拋物面鏡的軸線上存在這樣一個特殊的觀察點：在該點之前，中心波和邊緣波是分離的(見圖2(b)、2(c))；在該點處，中心波的負(fù)壓部分和邊緣波的負(fù)壓部分恰好完全重疊而使得聲信號的負(fù)壓幅值達(dá)到極大；在該點之后，中心波的正壓部分與邊緣波的負(fù)壓部分或中心波的負(fù)壓部分與邊緣波的正壓部分重合使得有效聲壓的幅值減小。滿足遠(yuǎn)場條件后，反射波演化為聲源波形的導(dǎo)數(shù)形式，聲壓幅值與傳播距離成反比。

從圖7中可以看出，隨著拋物面鏡開角增大，遠(yuǎn)場反射波的聲功率密度也越大。因為開角增加，拋物面鏡的深焦比d/zF也增加，這就相當(dāng)于使拋物面鏡的焦點向底部移動，使得更多的聲線受到鏡面的反射并使得遠(yuǎn)場的反射波聲功率密度增加。拋物面鏡開角增大，近場反射波的聲功率密度反而減小。其原因在于：當(dāng)開角較小(θ<90°)時，拋物面鏡較“淺”，焦點位于開口面之外。此時從焦點發(fā)出的球面波到達(dá)拋物面鏡頂點和邊緣點的聲程差很小，反射中心波與邊緣波的幅值接近，并且兩者在達(dá)到觀察點的時間上較為一致，于是中心波的負(fù)壓部分與邊緣波的負(fù)壓部分相互疊加形成了與聲源波形導(dǎo)數(shù)形式相同的波形。如前文所述，這種反射波形對應(yīng)的聲功率是最大的。如果開角增大，焦點向拋物面鏡的底部移動，從焦點發(fā)出的球面波到達(dá)邊緣點的聲程將比到達(dá)頂點的聲程大，從而使得邊緣波的幅值比中心波的小，并且在到達(dá)觀察點的時間上邊緣波滯后于中心波，于是反射波對應(yīng)的聲功率也逐漸減小。當(dāng)觀察點位于遠(yuǎn)場時，反射波的傳播距離很遠(yuǎn)，由拋物面鏡開角變化引起的聲程差可以忽略不計，反射波的波形及聲功率的大小遵循式(13)的規(guī)律。

在固定的觀察點處，本文進(jìn)一步研究了反射波聲功率隨拋物面鏡深焦比的變化，如圖8所示。

圖8　拋物面鏡反射聲場的最優(yōu)深焦比參數(shù)Fig.8　Optimum depth-to-focal-length ratios for the reflection sound of a parabolic mirror

分析結(jié)果發(fā)現(xiàn)，在拋物面鏡軸線上存在一個最優(yōu)的深焦比參數(shù)值使得該處反射波的聲功率最大(圖8(a))。深焦比參數(shù)超過該值，反射波的聲功率將出現(xiàn)一定程度的下降。例如，當(dāng)聲波頻率f0=100 kHz，拋物面鏡焦距zF=0.05 m時，在z/zF=100處最優(yōu)深焦比參數(shù)d/zF=1.4；在z/zF=200處時，最優(yōu)深焦比參數(shù)d/zF=2.6。上述最優(yōu)深焦比參數(shù)對于拋物面鏡的設(shè)計具有實際意義，它可以避免在設(shè)定傳播距離上將拋物面鏡設(shè)計得過大。隨之我們計算了不同觀察點處的最優(yōu)深焦比參數(shù)(圖8(b))，發(fā)現(xiàn)最優(yōu)深焦比參數(shù)與傳播距離之間呈線性關(guān)系。并且隨著聲信號頻率的升高，最優(yōu)深焦比參數(shù)曲線的斜率減小。這說明對于高頻聲波不需要將拋物面鏡設(shè)計得過大；而對于低頻聲波，想要提高換能器的發(fā)射效率則需要將拋物面鏡設(shè)計得很大。

5　結(jié)　論

工程上使用拋物面鏡可提高換能器的發(fā)射效率。假設(shè)聲波波長與拋物面鏡的曲率半徑相比很小，傳播介質(zhì)無耗散、各向均勻，首先根據(jù)Kirchhoff-Helmholtz聲衍射理論獨立推導(dǎo)了沿拋物面鏡軸線球面波反射聲場的時域理論解，預(yù)測近場反射波包含 “中心波”、“邊緣波”和“尾波”，中心波的相位與邊緣波和尾波相反。遠(yuǎn)場反射波的波形與聲源波形的導(dǎo)數(shù)形式相同，聲壓幅值與傳播距離成反比而與聲波的頻率成正比。理論解的正確性得到了COMSOL軟件數(shù)值模擬的驗證。其次，以單個周期的正弦波為例分析了反射聲場的形成過程。最后，基于推導(dǎo)的理論結(jié)果研究了球面波反射聲場的特性，討論了拋物面鏡開角和深焦比參數(shù)對反射聲場的影響。開角越大，遠(yuǎn)場反射波的聲功率越大而近場反射波的聲功率越小。固定拋物面鏡的口徑，則存在一個最優(yōu)的深焦比參數(shù)d/zF=3.92使得遠(yuǎn)場反射波聲功率最大。固定拋物面鏡的焦距，在特定傳播距離上也存在一個最優(yōu)的深焦比參數(shù)值。該值與傳播距離之間成線性關(guān)系，并且隨聲波頻率的變化而改變。

[1]吳慧云,吳武明,陳金寶,等.100 kW功率固體激光中繼鏡系統(tǒng)對1 km高度目標(biāo)作用效果模擬[J].光學(xué)學(xué)報,2008,28(10):1967-1970.

WU Huiyun,WU Wuming，CHEN Jinbao,et al.Simulation of effect of 100 kW solid-state laser relay mirror system on 1 km altitude target[J].ACTA Optical Sinca,2008,28(10):1967-1970.

[2]張軍，曾新吾，陳聃，等.水下強(qiáng)聲波脈沖負(fù)壓的產(chǎn)生和空化氣泡運動[J].物理學(xué)報，2012,18:184302.

ZHANG Jun,ZENG Xinwu,CHEN Dan,et al.Generation of negative pressure in water an cavitation bubble dynamics [J].ACTA Physica Sinca,2012,18:184302.

[3]王鴻樟.換能器與聚焦系統(tǒng)[M].上海交大出版社,1995,243.

[4]WAHLSTROM S.The parabolic reflector as an acoustical amplifier[J].J.Audio Eng.Soc.,1985,33(6):41-429.

[5]DAI X,ZHU J,TSAI Y T,et al.Use of parabolic reflector to amplify in-air signals generated during impact-echo testing[J].J.Acoust.Soc.Am.,2011,130:167-172.

[6]HAMILTON M F.Transient axial solution for the reflection of a spherical wave from a concave ellipsoidal mirror[J].J.Acoust.Soc.Am.1993,93(3):1256-1266.

[7]TSAI Y T,ZHU J,HABERMAN M R.Transient axial solution for plane and axisymmetric waves focused by a paraboloidal reflector[J].J.Acoust.Soc.Am.,2013,133(4):2025-2035.

[8]FRANCO E E,ANDRADE M A,ADAMOWSKI J C.Acoustic beam modeling of ultrasonic transducers and arrays using the impulse response and the discrete representation methods[J]J.of the Braz.Soc.of Mech.Sci.& Eng.,2011,33(4):408-416.

[9]PIWAKOWSKI B,SBAI K.A new approach to calculate the field radiated from arbitrarily structured transducer arrays[J].IEEE Transactions on Ultrasonics,Ferroelectrics,and Frequency Control,1999,46(2):422-440.

[10]BLACKSTOCK D T,MORFEY C L.Reflection and transmission of spherical waves incident on a concentric spherical interface[J].J.Acoust.Soc.Am.,1991,89:1971.

[11]ROBINSON D E,LEES S,BESS L.Near field transient radiation patterns for circular pistons[J].IEEE Transactions on Acoustics,Speech,and Signal Processing,1974,22(6):395-403.

Reflection sound field’s characteristics for spherical sound wave from a parabolic mirror

ZHANG Jun1,CHEN Peng1,CHEN Zhengwu1,ZHAO Yun2,ZENG Xinwu2

(1.Key Laboratory of Aerodynamic Noise Control,China Aerodynamics Research and Development Center,Mianyang 621000,China；2.College of Opto-electronic Science and Engineering,National University of Defense Technology,Changsha 410073,China)

Based on Kirchhoff-Helmholtz sound diffraction theory,a transient axial solution to reflection sound field of spherical acoustic wave from a parabolic mirror was derived.It was assumed that the spherical ware’s wavelength be small compared to the dimensions of the mirror and the geometric acoustic theory could be applied.The theoretical solution indicated that the near field reflection wave consists of three parts,i.e.,center wave,edge wave and wake wave; the phase of the center wave is opposite to those of the edge wave and wake one; the wave form of the far field reflection wave is consistent with that of the derivative of the sound source wave,and its sound pressure amplitude is inversely proportion to propagation distance and direct proportion to frequency.Taking a sinusoidal wave as an example,evolution and formation processes of its reflection sound field were presented.The existence of three wavelets and the correctness of the theoretical solution were verified with simulations using the software COMSOL.Lastly,the characteristics of the reflection sound field were studied,and the effects of geometric parameters of the mirror on the characteristics were discussed.It was shown that if the aperture of the mirror is fixed,an optimal depth-to-focal-length ratio d/zF=3.92 exists to maximize the sound power density of the far field reflection wave.

parabolic mirror; spherical sound wave; reflection sound field

國家自然科學(xué)基金(11504417；11172007)

2015-02-11修改稿收到日期：2015-08-06

張軍 男，博士，助理研究員，1983年11月生

曾新吾 男，教授，博士生導(dǎo)師，1963年4月生

TB54

A DOI:10.13465/j.cnki.jvs.2016.17.022