楊德森，孫玉，胡博，韓闖，靳仕源

（哈爾濱工程大學(xué)水聲技術(shù)重點實驗室，黑龍江哈爾濱150001）

阻抗梯度變化介質(zhì)的聲學(xué)特性

楊德森，孫玉，胡博，韓闖，靳仕源

（哈爾濱工程大學(xué)水聲技術(shù)重點實驗室，黑龍江哈爾濱150001）

阻抗梯度變化介質(zhì)是從多層復(fù)合介質(zhì)領(lǐng)域脫穎而出的一種聲學(xué)介質(zhì)，它是通過選取具有不同特征阻抗的材料，從阻抗匹配的角度組成特性阻抗成梯度變化的介質(zhì)，從而獲得更好的聲學(xué)性能。針對水中聲學(xué)介質(zhì)的聲傳播問題，建立了阻抗梯度變化介質(zhì)中的聲傳播模型，對其聲學(xué)特性進(jìn)行了研究，推導(dǎo)了阻抗梯度變化介質(zhì)的反射、吸聲系數(shù)的計算公式。通過仿真計算，將阻抗梯度變化介質(zhì)與傳統(tǒng)均勻介質(zhì)的聲學(xué)性能相對比，討論離散計算中介質(zhì)層數(shù)的選取原則。制備了阻抗梯度變化介質(zhì)，進(jìn)行了水下聲學(xué)性能測量實驗，實驗結(jié)果表明：實驗介質(zhì)的阻抗具有良好的連續(xù)性分布能夠滿足實驗要求，測量結(jié)果驗證了阻抗梯度變化介質(zhì)更小的反射系數(shù)和更大的吸聲系數(shù)，為阻抗梯度變化介質(zhì)在水聲工程中的應(yīng)用提供了理論依據(jù)。

聲阻抗；梯度變化；聲學(xué)特性；反射系數(shù)；阻抗匹配；聲學(xué)介質(zhì)

眾所周知，提高分層介質(zhì)聲學(xué)性能最直接的方法就是使介質(zhì)表面的輸入阻抗與自由空間波阻抗相匹配，同時要求介質(zhì)具有很高的聲衰減性能［1-2］，但這2個條件本身就是矛盾的，對于均勻介質(zhì)無法同時滿足。深入研究發(fā)現(xiàn)，當(dāng)介質(zhì)的阻抗或結(jié)構(gòu)呈梯度變化時能夠滿足以上2個條件。據(jù)此，通過改變介質(zhì)的結(jié)構(gòu)、物理和化學(xué)特性的方法，人們設(shè)計了多種聲學(xué)介質(zhì)，并對介質(zhì)的聲學(xué)性能進(jìn)行了研究。

Beretti［3］研究了具有4層結(jié)構(gòu)的消聲層的聲學(xué)特性，發(fā)現(xiàn)在不增加消聲材料厚度的情況下，由多種材料組成的多層結(jié)構(gòu)具有消聲和減振功能。Emer-y［4］對阻抗?jié)u變的多層結(jié)構(gòu)進(jìn)行了研究，指出在高頻時阻抗匹配可以減小反射系數(shù)，低頻時通過采用不同厚度、阻抗的多層結(jié)構(gòu)可以解決反射系數(shù)的峰值問題。Andrew［5］針對不同高度會引起空氣中聲速和阻抗的變化這一現(xiàn)象，對多個高度下大氣層的聲傳播問題進(jìn)行了研究。何祚鏞［6-7］對聲垂直入射時水下均勻和非均勻材料復(fù)合層的聲學(xué)特性進(jìn)行了研究，提出了計算包含任意非均勻?qū)拥亩鄬訌?fù)合結(jié)構(gòu)反射、吸聲系數(shù)的近似方法。朱金華［8］選取3種不同材料組合成梯度聚氨酯材料，研究了梯度聚氨酯介質(zhì)中的聲吸收，但結(jié)果發(fā)現(xiàn)這種制作方法并非十分理想，不同阻抗的材料分界面會影響聲波的發(fā)射和透射。楊雪［9-10］采用澆注的方法制備了梯度聚氨酯，研究在全入射角下梯度聚氨酯的吸聲性能，但這種逐層澆注的方法仍然存在每層材料的阻抗不連續(xù)的問題。王源升［11-12］設(shè)計合成了多種不同配方的梯度高分子溶液，通過聲波在梯度層中被多次反射吸收實現(xiàn)聲波有效衰減，但這種溶液在實用性和可控性方面還存在很大問題。目前，人們的研究只限于對不同材料組成的梯度聲學(xué)介質(zhì)的可行性研究，而對于阻抗梯度變化介質(zhì)中聲傳播機(jī)理問題，人們沒有作更深入的研究，已有的方法也只是將介質(zhì)分層，再將每層介質(zhì)中的阻抗看成均勻的進(jìn)行求解，而這種計算方法會存在很大誤差。在制備方面，理想的阻抗梯度變化介質(zhì)的阻抗應(yīng)連續(xù)梯度變化，內(nèi)部不存在明顯的界面，從而使介質(zhì)的性質(zhì)和功能沿厚度方向也成梯度變化，但目前制備的梯度介質(zhì)仍存在阻抗不連續(xù)和實用性差等問題。

本文以水中消聲材料為研究背景，設(shè)計了阻抗連續(xù)梯度分布的不均勻介質(zhì)，達(dá)到介質(zhì)表面阻抗匹配和內(nèi)部高損耗的目的。建立了阻抗梯度變化介質(zhì)中的聲傳播模型，對介質(zhì)中聲傳播規(guī)律進(jìn)行了研究，推導(dǎo)了任意阻抗分布條件下介質(zhì)的吸聲系數(shù)、反射系數(shù)公式。進(jìn)行仿真研究，將阻抗梯度變化介質(zhì)結(jié)構(gòu)與傳統(tǒng)均勻介質(zhì)的聲學(xué)性能比較。最后設(shè)計制作了阻抗連續(xù)梯度變化介質(zhì)，采用脈沖聲管法對其水下聲學(xué)性能進(jìn)行了測量。

1 阻抗梯度變化介質(zhì)中聲傳播理論

1.1 聲傳播規(guī)律

假設(shè)有2個半空間無限均勻介質(zhì)，分別為介質(zhì)I和介質(zhì)III，阻抗分別為Z1和Z3，密度分別為ρ1和ρ3，聲速分別為c1和c3。在它們之間有介質(zhì)II，介質(zhì)II的特性阻抗Z（x）設(shè)計為在迎聲面處與Z1匹配，在介質(zhì)另一側(cè)與Z3匹配，介質(zhì)II的厚度為d。為滿足介質(zhì)阻抗匹配和高衰減的要求，可以將介質(zhì)II的阻抗值設(shè)計成在厚度方向上按一定規(guī)律呈連續(xù)變化。圖1為介質(zhì)II的阻抗分布圖。設(shè)計時主要考慮以下幾點：1）介質(zhì)表面阻抗匹配。只有先將聲波“引”進(jìn)來，才有可能在介質(zhì)內(nèi)部對其進(jìn)行吸收和衰減，因此應(yīng)保證介質(zhì)表面的輸入阻抗與自由空間的阻抗匹配；2）保證阻抗連續(xù)變化。應(yīng)保證介質(zhì)內(nèi)部阻抗變化的連續(xù)性，減小介質(zhì)內(nèi)部聲波的反射；3）實現(xiàn)高衰減高吸收。當(dāng)聲波被引入介質(zhì)內(nèi)部后，必須具有足夠大的聲能量損耗才能對聲波進(jìn)行衰減吸收，因此將介質(zhì)的阻抗設(shè)計成逐漸增大變化。最終使阻抗梯度變化介質(zhì)獲得更小的反射系數(shù)更小和更大的吸聲系數(shù)。

圖1 介質(zhì)阻抗梯度分布圖Fig.1 The gradient change of acoustic impedance

設(shè)平面波從介質(zhì)I垂直入射，經(jīng)介質(zhì)II，從介質(zhì)III透射出，除介質(zhì)III外所有介質(zhì)中都有2個波，其中一個波向正x方向傳播，另一個波向負(fù)x方向傳播，在介質(zhì)III中只有一個向正x方向傳播的波。在阻抗梯度變化介質(zhì)中阻抗是連續(xù)變化的，因此，為了便于計算需要進(jìn)行離散化，將介質(zhì)II分為N層，每層介質(zhì)厚度為Δxi。需要注意的是，模型中介質(zhì)II的阻抗是連續(xù)變化的，在離散后考慮到實際情況，每層介質(zhì)中的阻抗、密度和聲速也應(yīng)當(dāng)是變化的。為保證離散后介質(zhì)模型的聲學(xué)性能不發(fā)生改變，給出幾點假設(shè)：1）聲壓和法向質(zhì)點振速在每層介質(zhì)的分界面上是連續(xù)的，即分界面兩側(cè)介質(zhì)的聲壓和法向質(zhì)點振速分別相等；2）介質(zhì)II的左側(cè)和右側(cè)表面分別處于半無限空間中；3）為方便計算不考慮二次反射及多次反射的情況。

由理想流體媒質(zhì)中小振幅聲波運動方程可知：

式中：p為介質(zhì)質(zhì)點的聲壓，v為介質(zhì)質(zhì)點的速度，ρ為密度。當(dāng)考慮介質(zhì)具有粘滯性時［13］，可以得到聲壓p、質(zhì)點振速v與位移ξ的關(guān)系為

式中：K為復(fù)體積模量，表示為

式中：KS為體積彈性系數(shù)，η為粘滯系數(shù)。已知理想流體介質(zhì)中一維波動方程的解為

式中：c～為復(fù)聲速，p＋（x，t）代表了沿正x方向行進(jìn)的波，p－（x，t）代表了沿負(fù)x方向行進(jìn)的波。定義聲場中某位置的波阻抗Z（x）為

根據(jù)歐拉方程，式（5）可以得到質(zhì)點振速為

一維波動方程的解為關(guān)于時間的正弦函數(shù)，利用分離變量法，方程解寫為

其中，

將式（8）中的v（x，t）代入式（1）中，由于分離變量后公式左右兩側(cè)都有時間因子，因此可以消去時間因子得到

將式（5）和式（7）代入式（10）得到

分別對式（2）和式（8）中的p（x，t）對t求偏導(dǎo)，聯(lián)立式（3）可得

將式（5）和式（7）代入到式（12）中得到

利用式（6），式（12）可簡化為

由式（11）與式（14），可分別得到

將式（5）寫成向量的形式：

對式（17）求x偏導(dǎo)再代入式（15）和式（16），簡化可以得到

式中：α（x）為介質(zhì)阻抗的相對變化。

如果要求得微分方程式（18）的解，必須保證矩陣C（x）為常數(shù)，但在本模型中聲速c和α（x）都是隨x變化的，無法得到微分方程的解，因此必須對式（18）進(jìn)行進(jìn)一步變化。設(shè)聲波透過阻抗梯度變化介質(zhì)的總時間為T，考慮到聲速不僅和傳播距離x有關(guān)，而且還與傳播時間t有關(guān)，因此可以進(jìn)行變量替換，把C（x）轉(zhuǎn)化為隨時間t變化的函數(shù)。將式（18）變?yōu)?/p>

其中，當(dāng)t＝0時，x＝0。

式中：

為求得微分方程式（21）的解，則β（t）一定為常數(shù)。當(dāng)β（t）為常數(shù)時對式（22）積分，得到

式（23）說明，在本介質(zhì)模型中為了保證式（21）有解，介質(zhì)的阻抗應(yīng)滿足關(guān)于t的指數(shù)分布。而對于式（23）可以這樣理解，即介質(zhì)阻抗是密度和聲速的乘積，而介質(zhì)中的聲速是關(guān)于聲波傳播時間t的函數(shù)，因此，也可以把介質(zhì)阻抗看成是關(guān)于時間t的函數(shù)。

1.2 介質(zhì)聲學(xué)性能

設(shè)離散后介質(zhì)II中任意一層介質(zhì)分界面位置為x＝xM和x＝xN，xN－xM＝Δx，每層介質(zhì)的入射聲壓和透射聲壓可以通過傳遞函數(shù)G（xM，xN）聯(lián)系起來，即

對于式（24）將xM用x代替，可以得到最為普遍的表達(dá)式為

將式（25）對t求導(dǎo)可得

將式（25）代入式（21）中得到

對比式（26）和（27）可以得到

由1.1節(jié)內(nèi)容可知，當(dāng)介質(zhì)阻抗為關(guān)于t的指數(shù)分布時微分方程有解，同理當(dāng)介質(zhì)阻抗為關(guān)于t的指數(shù)分布時上式的解可以表示為

式中：ω*＝ω2－β2。介質(zhì)II總傳遞函數(shù)Gtotal為

其中，Gn是第n層介質(zhì)的傳遞函數(shù)。設(shè)介質(zhì)中聲速和密度也為關(guān)于t指數(shù)分布，即

由式（31）可以得到

整理得到

對于式（28），只有介質(zhì)阻抗?jié)M足隨指數(shù)分布時微分方程才有解，而對于阻抗為一般分布時，卻不能得到解析解。為解決這一問題，使得任意阻抗分布也滿足微分方程的條件，可以采用分段近似擬合的方法。該方法是將阻抗梯度變化介質(zhì)分為N層，當(dāng)N的取值足夠大時，每層介質(zhì)內(nèi)的阻抗分布可由指數(shù)函數(shù)進(jìn)行近似代替，并保證真實值與近似值的誤差非常小，整個介質(zhì)阻抗分布再由每層介質(zhì)中的阻抗進(jìn)行擬合。因此對于第n層介質(zhì)，已知介質(zhì)中任意兩點的阻抗Z（xn）和Z（xn＋1），再根據(jù)式（23）就可以得到近似的指數(shù)函數(shù)為

式中：bn為第n層介質(zhì)中近似指數(shù)函數(shù)的系數(shù)，Δt為第n層及之內(nèi)聲波傳播時間。整理得到bn為

這樣通過分段近似擬合，任意阻抗分布都滿足式（27）微分方程解的存在條件，因此第n層介質(zhì)的傳遞函數(shù)Gn表示為

以上給出了每層介質(zhì)中的傳遞函數(shù)，考慮到因為在介質(zhì)III中只有一個透射整個層系而向正x方向傳播的波，并且考慮梯度介質(zhì)左右界面處的邊界條件，可以得到

式中：Mtotal＝Q1GtotalQ2，Q1和Q2為

定義

則介質(zhì)II的透射系數(shù)和反射系數(shù)為

同樣，可以得到整個介質(zhì)II的反射系數(shù)為

以上推導(dǎo)了阻抗梯度變化的非均勻介質(zhì)中聲傳播特性公式。計算時首先確定聲波通過介質(zhì)的總時間，再求得每個時間點所對應(yīng)的位置xn，然后計算每層介質(zhì)中近似指數(shù)函數(shù)的系數(shù)bn和傳遞矩陣Gn，最終得到介質(zhì)的反射系數(shù)和透聲系數(shù)。

2 介質(zhì)聲學(xué)性能仿真

2.1 方法驗證

利用Matlab軟件編制程序進(jìn)行數(shù)值仿真。設(shè)圖1中的介質(zhì)I為水，介質(zhì)II為阻抗梯度變化的介質(zhì)，介質(zhì)III為鋼板，鋼板足夠厚可以看成半無限空間，水的特性阻抗為Z1＝1.5×106kg／（s·m2），鋼板的特性阻抗為Z3＝46.25×106kg／（s·m2）。將阻抗梯度變化介質(zhì)結(jié)構(gòu)與3種不同阻抗的均勻介質(zhì)比較，對比其聲學(xué)性能。設(shè)平面波垂直入射，介質(zhì)的厚度d均為0.1 m，離散層數(shù)為200層，均勻介質(zhì)分別選取為橡膠、環(huán)氧樹脂和鉛，阻抗分別為1.67×106kg／（s·m2）、3.2×106kg／（s·m2）和16.9×106kg／（s·m2）。已知介質(zhì)II的阻抗分布滿足指數(shù)分布：

圖2給出了4種介質(zhì)的反射系數(shù)和吸聲系數(shù)對比結(jié)果?？梢钥吹?，3種均勻介質(zhì)的聲學(xué)性能都具有頻率選擇性，聲學(xué)性能曲線作振幅相等的周期振蕩，而阻抗梯度變化介質(zhì)的聲學(xué)性能曲線并沒有出現(xiàn)振蕩現(xiàn)象。在低頻段，阻抗梯度變化介質(zhì)的聲學(xué)性能曲線與4種介質(zhì)的基本相同；在1～10 kHz頻率范圍內(nèi)，隨頻率增加，阻抗梯度變化介質(zhì)的反射系數(shù)單調(diào)遞減至0.8，而吸聲系數(shù)單調(diào)增大至0.4；在10 kHz以上的高頻段，阻抗梯度變化介質(zhì)的反射系數(shù)由0.8減小至0.2，反射系數(shù)要遠(yuǎn)小于均勻介質(zhì)，而其吸聲系數(shù)由0.4增大至0.9，吸聲系數(shù)要遠(yuǎn)大于均勻介質(zhì)，均好于3種均勻介質(zhì)。分析原因：1）在低頻段，聲波波長遠(yuǎn)大于介質(zhì)厚度，無論是均勻介質(zhì)還是阻抗梯度變化介質(zhì)，介質(zhì)都如同“不存在”，因此低頻段二者的聲學(xué)性能基本相同；2）在中高頻段，對于均勻介質(zhì)，其阻抗接近于水的特性阻抗使得聲波容易進(jìn)入介質(zhì)內(nèi)，但在分界面處，鋼板與均勻的特性阻抗相差很大，最大值達(dá)到30倍，此時聲波將發(fā)生強(qiáng)烈反射。而對于阻抗梯度變化介質(zhì)，它的表面阻抗接近于水的特性阻抗，有利于聲波進(jìn)入介質(zhì)內(nèi)部，在聲波前進(jìn)方向上，每一層的阻抗都比前一層略大一個微小量，起到阻抗過渡作用，并利用逐漸增大的阻抗對聲波進(jìn)行衰減，同時在分界面處，鋼板與阻抗梯度變化介質(zhì)的阻抗比值較小，僅為Z3／Z（0.1）≈1.1，能夠?qū)Ψ瓷渎暡ㄟM(jìn)行抑制，實現(xiàn)減少反射增加吸聲的目的；3）阻抗梯度變化介質(zhì)是由不同阻抗介質(zhì)進(jìn)行合理疊加而成，不同阻抗的介質(zhì)能夠在聲學(xué)性能上形成互補(bǔ)，因此整體具有較寬的頻帶寬度。

圖2 Z3＝46×106kg／（s·m2）時4種介質(zhì)聲學(xué)性能對比Fig.2 Comparison of acoustic characteristics of four media when Z3＝46×106kg／（s·m2）

為驗證模型的一般性，減小介質(zhì)III的阻抗，使得Z3＝20×106kg／（s·m2），圖3為4種介質(zhì)對比結(jié)果。通過對比可以發(fā)現(xiàn)，均勻介質(zhì)的透聲面處的聲波反射減小，吸聲性能有所提高，但在中高頻段阻抗梯度變化介質(zhì)仍能夠獲得較小的反射系數(shù)和較大的吸聲系數(shù)?？偨Y(jié)以上結(jié)果表明：與傳統(tǒng)均勻介質(zhì)相比，在不增加吸聲介質(zhì)厚度的情況下，在中高頻段阻抗梯度變化介質(zhì)具有更低的反射系數(shù)和更大的吸聲系數(shù)。

為了進(jìn)一步驗證算法的正確性，設(shè)阻抗梯度變化介質(zhì)阻抗分布的阻抗變化率為0，此時分別采用分層介質(zhì)傳統(tǒng)算法和本文方法對同一均勻介質(zhì)進(jìn)行計算，結(jié)果如圖4所示。從圖中可知，當(dāng)阻抗梯度變化介質(zhì)的阻抗變化率為0時，阻抗梯度變化介質(zhì)變?yōu)榫鶆蚪橘|(zhì)，此時2種方法的計算結(jié)果是吻合的，這說明本文方法對于均勻介質(zhì)也同樣適用，能夠驗證本文算法的正確性。

圖3 Z3＝20×106kg／（s·m2）時4種介質(zhì)聲學(xué)性能對比Fig.3 Comparison of acoustic characteristics of four media when Z3＝20×106kg／（s·m2）

圖4 傳統(tǒng)方法和本文算法對比Fig.4 Comparison of traditional method and this method

2.2 介質(zhì)層數(shù)影響

在離散計算中，為保證任意阻抗分布都能滿足微分方程求解要求，需要保證N的取值足夠大或每層介質(zhì)厚度Δx足夠小，使得每層介質(zhì)阻抗分布可由指數(shù)分布近似，進(jìn)而介質(zhì)的總阻抗分布可由每層介質(zhì)的阻抗擬合而成。這就引出人們最為關(guān)心的一個問題，即N選取為多大或每層介質(zhì)厚度為多小時才能保證近似值與阻抗真實分布的誤差足夠小，不會影響介質(zhì)的聲傳播規(guī)律和聲學(xué)性能。以下對這一問題進(jìn)行深入研究。

計算時首先將傳播時間均勻分為N份，再計算每個時間間隔內(nèi)介質(zhì)的厚度Δxn，由于每層介質(zhì)中聲速是逐漸增大的，在相同時間間隔內(nèi)所對應(yīng)的Δxn值也是逐層增大的，如果采用分析每層介質(zhì)厚度與波長比的方法將無法對其變化影響進(jìn)行有效說明，因此這里考慮介質(zhì)層數(shù)N的選取問題。圖5給出了N變化對介質(zhì)聲學(xué)性能的影響?？梢钥吹?，當(dāng)N＜100時反射系數(shù)曲線和吸聲系數(shù)曲線有明顯的變化；當(dāng)N＞100時反射系數(shù)和吸聲系數(shù)曲線非常平坦，介質(zhì)層數(shù)N由100增大到200時，介質(zhì)的反射系數(shù)僅發(fā)生了0.005的變化。

圖5 介質(zhì)層數(shù)對介質(zhì)聲學(xué)性能影響Fig.5 The relationship between N and acoustic characteristics of media

這說明：當(dāng)N＞100時指數(shù)近似的結(jié)果與實際阻抗分布之間的誤差已經(jīng)很小，對介質(zhì)的聲學(xué)性能影響也很小，此時N的選擇能夠滿足模型假設(shè)條件中足夠大的要求，介質(zhì)的阻抗分布可由每層介質(zhì)的阻抗近似擬合而成。

3 實驗研究與分析

為驗證理論和仿真研究的正確性，設(shè)計制作了阻抗梯度變化介質(zhì)實驗樣品，采用脈沖聲管法對介質(zhì)的水下聲學(xué)性能進(jìn)行測量。

3.1 實驗介質(zhì)的制備

首先，基于沉降法理論制備了阻抗梯度變化實驗介質(zhì)。制備裝置主要包括：沉降玻璃管、沉降模具、干燥箱、萬能拉伸試驗機(jī)、保溫箱、混合器皿，選用鋁粉和鐵粉作為制備材料。具體制備步驟如下：1）將沉降玻璃管和沉降模具內(nèi)表面擦洗干凈，并在沉降模具上均勻涂抹適量的脫模劑；2）稱量待包覆粉末和高分子粘結(jié)劑，高分子粘結(jié)劑為酚醛樹脂與烏洛托品的混合物；3）使粉末與粘結(jié)劑混合均勻，再將混料容器加熱繼續(xù)熱混；4）將包裹粉末迅速置入有機(jī)玻璃管中，管中盛滿水，粉末進(jìn)行自然沉降，待粉末完全沉降后，將沉降管中的水排出；5）將沉降模具和沉積粉末置于干燥箱中進(jìn)行干燥處理；6）將干燥后模具置于萬能拉伸試驗機(jī)下進(jìn)行壓制，壓制出相應(yīng)厚度；7）將模具固定后，放入保溫箱中在200℃下進(jìn)行固化處理。制備后的阻抗梯度變化介質(zhì)實物如圖6所示。

需要說明的是，本文目前掌握的制備技術(shù)和知識水平，很難解決由于材料互溶性不好導(dǎo)致固化后介質(zhì)非常脆、容易斷裂的問題。因此，在制作時綜合考慮材料之間的互溶性、材料的阻抗差異、固化效果、硬度等因素，在保證介質(zhì)阻抗連續(xù)變化的前提條件下提高了介質(zhì)迎聲面的阻抗值，制作完畢后對介質(zhì)的阻抗進(jìn)行測量，介質(zhì)的阻抗分布情況如圖7所示?？梢钥吹剑苽涞膶嶒灲橘|(zhì)的阻抗值在6.3×106～16.05×106kg／（s·m2）范圍之間，介質(zhì)具有明顯的梯度分布，隨著距離的增加介質(zhì)阻抗分布并沒有出現(xiàn)突變，基本上呈連續(xù)變化，能夠滿足實驗要求。

圖6 制備的阻抗梯度變化介質(zhì)Fig.6 The media with gradient change of impedance

圖7 介質(zhì)的阻抗分布Fig.7 Acoustic impedance value of media

3.2 實驗數(shù)據(jù)處理

依據(jù)國家標(biāo)準(zhǔn)GB／T14369-93，采用脈沖聲管法對介質(zhì)的水下聲學(xué)性能進(jìn)行測量。整個測試系統(tǒng)由脈沖聲管、收發(fā)換能器、信號源、B＆K2713功率放大器、收發(fā)轉(zhuǎn)換裝置、B＆K2636測量放大器、B＆K1617帶通濾波器和示波器組成，實驗測量系統(tǒng)如圖8所示。脈沖聲管的聲硬末端采用大質(zhì)量的圓形鋼塊，鋼塊厚度為5 cm，直徑為5.7 cm；選用聚氨酯作為對比樣品，直徑為5.6 cm，厚度為1 cm，梯度介質(zhì)的厚度為1 cm。介質(zhì)的特點，將制備的1號阻抗梯度變化介質(zhì)與相同厚度的聚氨酯均勻介質(zhì)相比較，對比結(jié)果如圖9所示。

圖8 實驗測量系統(tǒng)Fig.8 Experimental measurement system

圖9 1號梯度介質(zhì)與聚氨酯對比結(jié)果Fig.9 Comparison of No.1 media and polyurethane

在圖9中，阻抗梯度變化介質(zhì)與阻抗均勻分布的聚氨酯的聲學(xué)性能有很大不同，聚氨酯的聲學(xué)性能曲線呈周期振蕩，反射系數(shù)在0.65～0.95內(nèi)波動變化，吸聲系數(shù)0.1～0.58內(nèi)波動變化；與聚氨酯相比，在低頻段阻抗梯度變化介質(zhì)的反射系數(shù)更小，吸聲系數(shù)更大，并且反射系數(shù)隨頻率的增大而減小，吸聲系數(shù)隨頻率的增大而增大，這與仿真結(jié)果得到的規(guī)律是相同的。值得注意的是，圖中阻抗梯度變化介質(zhì)的聲學(xué)性能曲線的變化規(guī)律并沒有仿真結(jié)果給出的那么理想，曲線中會出現(xiàn)不規(guī)則的凹陷和凸起，分析原因，這可能是因為仿真模型是建立在阻抗連續(xù)變化的基礎(chǔ)上的，而在樣品介質(zhì)制備過程中，由于制備工藝的原因會存在阻抗局部連續(xù)性不好的問題，此外實驗測量中某些頻率下反射波波形很不穩(wěn)定，信號幅值的讀取也會帶來誤差。

圖10為2號阻抗梯度變化介質(zhì)與聚氨酯的測量結(jié)果?？梢钥吹?，圖10的實驗結(jié)果與圖9中具有很好的一致性，雖然阻抗梯度變化介質(zhì)的阻抗分布發(fā)生了變化，但與均勻聚氨酯介質(zhì)相比，阻抗梯度變化介質(zhì)仍然獲得更小的反射系數(shù)，更大的吸聲系數(shù)，得到的實驗規(guī)律與仿真結(jié)果相同。以上結(jié)果表明，在不增加介質(zhì)厚度的情況下，與傳統(tǒng)均勻介質(zhì)相比，阻抗梯度變化介質(zhì)具有更大的中高頻吸聲特性。

圖10 2號介質(zhì)與聚氨酯對比結(jié)果Fig.10 Comparison of No.2 gradient media and polyurethane

4 結(jié)束語

本文對阻抗梯度變化介質(zhì)的聲學(xué)性能進(jìn)行了研究，建立了阻抗連續(xù)梯度變化介質(zhì)中聲波傳播模型，推導(dǎo)了阻抗任意分布條件下介質(zhì)中聲學(xué)性能公式。通過數(shù)值仿真在多種條件下驗證了本文方法，制備了阻抗梯度變化介質(zhì)，對介質(zhì)的阻抗分布情況進(jìn)行了測量，進(jìn)行了水中阻抗梯度變化介質(zhì)的聲學(xué)性能測量實驗。研究發(fā)現(xiàn)：

1）在不增加吸聲介質(zhì)厚度的情況下，在中高頻段阻抗連續(xù)梯度變化介質(zhì)具有更小反射系數(shù)和更大的吸聲系數(shù)；

2）當(dāng)N＞100時，介質(zhì)層數(shù)能夠滿足模型假設(shè)條件中足夠大的要求，以保證阻抗變化的連續(xù)性；

3）制備的阻抗梯度變化介質(zhì)能夠滿足實驗要求，證明了本文方法的可行性，實驗結(jié)果驗證了理論研究的正確性。

可以說本文的研究為阻抗梯度變化介質(zhì)在實際工程中的應(yīng)用提供了理論依據(jù)和實踐基礎(chǔ)，對水下減振降噪具有比較積極的意義。

參考文獻(xiàn)：

［1］劉志宏，盛美萍.水下多層均勻材料垂直入射的聲特性分析［J］.噪聲與振動控制，2005，4（2）：11-13.LIU Zhihong，SHENG Meiping.Analysis of sound characteristics of a normal incidence plane wave on isotropic symmetric plate in water［J］.Noise and Vibration Control，2005，4（2）：11-13.

［2］李海濤，朱錫，王林，等.水下均勻復(fù)合層結(jié)構(gòu)的聲特性研究［J］.武漢理工大學(xué)學(xué)報，2008，30（6）：105-109.LI Haitao，ZHU Xi，WANG Lin，et al.Investigation of sound field in homogeneous multiple-layer structures in water［J］.Journal of Wuhan University of Technology，2008，30（6）：105-109.

［3］BERETTI S，VASSAS M.Multifunction coating and enhanced anechoic coating by using multi-layer optimization［C］／／Conference Proceedings of UDT.London，1996：300-303.

［4］EMERY P A.New cladding material［C］／／Conference Proceedings of UDT.London，1995：527-531.

［5］PEPLOW A.Acoustic propagation in variable sound speed profiles［C］／／8th European Conference on Noise Control.Edinburgh，Scotland，2009：302-307.

［6］何祚鏞，王曼.水下均勻材料復(fù)合層吸聲特性研究［J］.應(yīng)用聲學(xué)，1996，15（5）：6-11.HE Zuoyong，WANG Man.Investigation of the sound absorption of homogeneous composite multiple-layer structures in water［J］.Applied Acoustics，1996，15（5）：6-11.

［7］何祚鏞，王曼.水下非均勻復(fù)合層結(jié)構(gòu)吸聲的理論研究［J］.應(yīng)用聲學(xué)，1996，15（5）：12-19.HE Zuoyong，WANG Man.Investigation of the sound absorption of non-homogeneous composite multiple-layer structures in water［J］.Applied Acoustics，1996，15（5）：12-19.

［8］朱金華，劉巨斌，姚樹人，等.分層高分子介質(zhì)中的聲吸收［J］.高分子材料科學(xué)與工程，2001，17（3）：34-38.ZHU Jinhua，LIU Jubin，YAO Shuren，et al.The acoustic absorption in layered polymer materials［J］.Polymer Materials Science and Engineering，2001，17（3）：34-38.

［9］楊雪，王源升，朱金華，等.梯度聚氨酯水聲吸聲性能的研究［J］.高分子材料科學(xué)與工程，2007，23（1）：169-172.YANG Xue，WANG Yuansheng，ZHU Jinhua，et al.Study on the underwater sound absorption properties of gradient polyurethanes［J］.Polymer Materials Science and Engineering，2007，23（1）：169-172.

［10］楊雪，王源升，余紅偉.梯度聚氨酯吸聲性能研究［J］.武漢理工大學(xué)學(xué)報，2008，30（3）：21-23.YANG Xue，WANG Yuansheng，YU Hongwei.Research on sound absorption properties of graded polyurethanes［J］.Journal of Wuhan University of Technology，2008，30（3）：21-23.

［11］王源升，楊雪，朱金華，等.水溶性高分子梯度溶液吸聲機(jī)理的研究［J］.聲學(xué)學(xué)報，2006，31（1）：14-18.WANG Yuansheng，YANG Xue，ZHU Jinhua，et al.Study on the sound absorption mechanism in gradient water-soluble polymer solution［J］.Acta Acustica，2006，31（1）：14-18.

［12］王源升，楊雪，朱金華，等.梯度高分子溶液的聲衰減［J］.高分子材料科學(xué)與工程，2005，21（5）：129-132.WANG Yuansheng，YANG Xue，ZHU Jinhua，et al.The underwater sound attenuation of gradient polymer solutions［J］.Polymer Materials Science and Engineering，2005，21（5）：129-132.

［13］杜功煥，朱哲民，龔秀芬.聲學(xué)基礎(chǔ)［M］.南京：南京大學(xué)出版社，2010：169-172.

Acoustic characteristics of the media with gradient change of impedance

YANG Desen，SUN Yu，HU Bo，HAN Chuang，JIN Shiyuan
（Science and Technology on Underwater Acoustic Laboratory，Harbin Engineering University，Harbin 150001，China）

The media with gradient change of impedance is an acoustic media which is developed from multiple layer media.Using the principle of impedance matching method，better acoustic properties based on the media with gradient change of impedance with different characteristic impedance of materials can be obtained.Applying this media to underwater acoustic engineering，the model of acoustic propagation in the media was established and the sound properties were investigated.The calculation formulas for the reflection and sound absorption coefficients of the media with gradient change of impedance were derived.Numerical simulations were also conducted.The acoustic properties of the media with gradient change of impedance were compared with traditional homogeneous media.The rules of the number of layers were discussed.The media with gradient change of impedance were prepared and the experiments in a pulse acoustic tube were carried out.The experimental results showed that the preparation media with continuous distribution of impedance can meet the need of experiment.The measure results verified smaller reflection coefficient and larger absorption coefficient of the media with gradient change of impedance，which provide a theoretical basis for the application of the media with gradient change of impedance in acoustic engineering.

acoustic impedance；gradient change；acoustic characteristic；reflection coefficient；impedance matching；acoustic media

10.3969／j.issn.1006-7043.201401002

http：／／www.cnki.net／kcms／doi／10.3969／j.issn.1006－7043.201401002.html

TB535

A

1006-7043（2014）12-1458-09

2014-01-02.網(wǎng)絡(luò)出版時間：2014-12-02.

長江學(xué)者和創(chuàng)新團(tuán)隊發(fā)展計劃資助項目（IRT1228）；國家自然科學(xué)基金資助項目（11204049，11204050）；高等學(xué)校博士學(xué)科點專項科研基金資助項目（20122304120023）.

楊德森（1957-），男，教授，博士生導(dǎo)師；胡博（1980-），男，副研究員.

胡博，E-mail：kidd1105＠sina.com.