靜壓下考慮腔壓的吸聲覆蓋層吸聲性能分析

董文凱，陳美霞

華中科技大學(xué) 船舶與海洋工程學(xué)院，湖北 武漢 430074

0 引 言

隱身性是確保艦艇作戰(zhàn)能力的重要特性，敷設(shè)在潛艇外殼表面的水下吸聲覆蓋層是潛艇聲隱身應(yīng)用技術(shù)的重要組成部分。典型的水下吸聲覆蓋層由橡膠基體和內(nèi)部空腔構(gòu)成，這種結(jié)構(gòu)通過波形轉(zhuǎn)換、空腔諧振等效應(yīng)在某些頻段能夠很好地吸收聲波能量，使其聲反射量達(dá)到最小，從而降低回聲和目標(biāo)強(qiáng)度，提升潛艇隱身性能。因此，對(duì)靜壓下水下吸聲覆蓋層的吸聲性能開展理論研究有著重要意義。

目前，有關(guān)含空腔吸聲覆蓋層的研究取得了相當(dāng)多的成果，這些研究主要采用解析方法和有限元法。例如，Gaunaurd[1]提出一種含短圓柱形腔體的黏彈性介質(zhì)的一維分析方法，其根據(jù)腔體的不同物理尺寸及模態(tài)的共振情況，對(duì)覆蓋層的吸聲性能進(jìn)行了物理解釋。朱蓓麗和黃修長[2]采用彈性波傳播理論對(duì)含圓柱腔覆蓋層進(jìn)行了分析，通過對(duì)頻散方程的簡化和近似處理，將含圓柱腔的橡膠介質(zhì)模型轉(zhuǎn)化為具有等效波數(shù)的均勻材料層，再結(jié)合傳遞矩陣法[3]求解整體聲學(xué)覆蓋層的聲學(xué)系數(shù)。何世平和湯渭霖等[4-6]在Kelvin-Voigt的線性黏彈性理論模型框架內(nèi)，分析無限/有限長圓柱管單元在聲波垂直入射下激發(fā)的軸對(duì)稱波傳播及其衰減規(guī)律，發(fā)現(xiàn)低階空腔的穿孔率對(duì)傳播和衰減有著相當(dāng)大的影響，穿孔率越小，吸聲系數(shù)的峰值越大，峰值頻率也越高。在此基礎(chǔ)上，作者還分析了變截面空腔覆蓋層的吸聲性能，將變截面空腔分解為多段和均勻的圓柱腔，在各子段邊界處建立了應(yīng)力及位移連續(xù)方程，最終結(jié)合入射、透射端面的邊界條件求解反射和吸聲系數(shù)。

然而，在高靜壓下，潛艇外殼表面敷設(shè)的水下吸聲覆蓋層在受到擠壓后，其形狀、材料參數(shù)將發(fā)生改變，導(dǎo)致吸聲性能受到較大影響。目前，研究靜壓下的吸聲覆蓋層主要考慮覆蓋層形變對(duì)吸聲的影響。例如，姜聞文等[7]應(yīng)用有限元法計(jì)算靜壓下幾種橡膠吸聲結(jié)構(gòu)的形變及其吸聲性能，結(jié)果顯示，隨著靜壓的增大，幾種典型腔體結(jié)構(gòu)的吸聲頻帶都朝高頻移動(dòng)，與實(shí)測結(jié)果的變化趨勢(shì)一致。陶猛等[8]利用有限元軟件分析靜壓下的吸聲覆蓋層腔體單元形變，并將形變量引入空腔吸聲覆蓋層的二維解析理論和傳遞矩陣法中，分析了靜壓對(duì)水下吸聲覆蓋層聲學(xué)性能的影響。張沖等[9]和楊立軍等[10]先后使用有限元軟件研究靜壓下球形、橢球形空腔吸聲覆蓋層的吸聲性能，發(fā)現(xiàn)這兩種形狀的空腔在靜壓下的穿孔率越大并不代表該覆蓋層的低頻吸聲性能就越好。

此外，對(duì)于靜壓下的橡膠參數(shù)變化，仍主要集中在實(shí)驗(yàn)研究層面。例如，黃修長等[11-12]建立了利用水聲聲管進(jìn)行靜壓下橡膠材料動(dòng)態(tài)力學(xué)參數(shù)測量的方法。目前，多數(shù)研究仍然將靜壓下的腔形近似為規(guī)則的幾何體，較少考慮空腔內(nèi)部的反力，且未能對(duì)吸聲曲線上的谷值進(jìn)行機(jī)理分析。

因此，本文擬將受壓后的有限元模型導(dǎo)出后重新賦予材料參數(shù)進(jìn)行建模，考慮空腔壓力對(duì)潛艇外殼表面的水下吸聲覆蓋層性能的影響，以消除因幾何近似帶來的誤差。通過三維建模和二維軸對(duì)稱建模的對(duì)比，簡化有限元建模以提高計(jì)算效率?；诖?，對(duì)比在考慮空腔壓力的情況下對(duì)吸聲系數(shù)的影響，以揭示靜壓下含空腔的覆蓋層吸聲系數(shù)曲線上產(chǎn)生尖銳谷值的機(jī)理。

1 理論模型

1.1 物理分析

圖1所示為潛艇外殼表面敷設(shè)含圓柱形空腔的吸聲覆蓋層。敷設(shè)后的整體結(jié)構(gòu)模型由4個(gè)部分組成（圖1（a）），包括入射端的水介質(zhì)、吸聲覆蓋層、鋼結(jié)構(gòu)層和終端的空氣介質(zhì)。其中，吸聲覆蓋層結(jié)構(gòu)的平面上有交錯(cuò)和平行兩種周期性排布結(jié)構(gòu)（圖1（b）），分別可以拆解為正四棱柱單元和正六棱柱單元。在單元橫截面積與孔徑率相等的情況下，其吸聲性能的等效性已經(jīng)得到了包括文獻(xiàn)[3]在內(nèi)的較多文獻(xiàn)的驗(yàn)證。

圖1 水下吸聲覆蓋層單元結(jié)構(gòu)圖Fig.1 Structural diagram of anechoic coating layer unit

本文以兩種棱柱單元為例，將吸聲覆蓋層結(jié)構(gòu)分解后近似為圓柱形單元，如圖2所示。

圖2 吸聲覆蓋層單元等效圖Fig.2 Equivalent diagram of anechoic coating layer unit

基于此，下文將主要討論含該圓柱形空腔的吸聲覆蓋層在300 m海水深度以內(nèi)時(shí)的吸聲效果。因橡膠具有楊氏模量小、泊松比大的特點(diǎn)，在靜壓作用下其易發(fā)生擠壓形變，故已不可忽略。在靜壓下含空腔的橡膠形變行為相當(dāng)復(fù)雜，此時(shí)，根據(jù)材料特性表征空腔橡膠形變行為的模型包括有線彈性、超彈性等模型。目前，主要采用有限元法計(jì)算靜壓下的橡膠形變，即將靜態(tài)材料參數(shù)代入有限元軟件中求解得到形變。而對(duì)于聲傳播的計(jì)算，由于橡膠材料的特殊性，在不同頻率的聲波激勵(lì)下橡膠材料的楊氏模量、泊松比、損耗因子是不同的，因此，計(jì)算聲波激勵(lì)下橡膠材料的吸聲系數(shù)時(shí)必須采用動(dòng)態(tài)材料參數(shù)。

在上述物理分析的基礎(chǔ)上，本文給出了如圖3所示的水下吸聲覆蓋層吸聲系數(shù)計(jì)算流程。首先，將橡膠靜態(tài)參數(shù)導(dǎo)入有限元軟件COMSOL來計(jì)算靜壓下的形變，以求解得到覆蓋層形變量；然后，從軟件中導(dǎo)出每段圓柱形空腔的形變量；最后，基于含圓柱形空腔的吸聲覆蓋層聲傳播一維理論，求解理論吸聲系數(shù)，并將形變后的網(wǎng)格導(dǎo)入聲?固耦合模塊，再引入動(dòng)態(tài)材料參數(shù)及建立聲場，計(jì)算有限元吸聲系數(shù)。

圖3 靜壓下吸聲覆蓋層吸聲系數(shù)的計(jì)算流程圖Fig.3 Flowchart of sound absorption coefficient calculation for anechoic coating layer under hydrostatic pressure

1.2 數(shù)學(xué)分析

圖4所示為空腔內(nèi)徑為rb、外徑為ra、管段長度為d的均勻圓柱形管，其可視為含空腔吸聲覆蓋層的一段。其中，空腔外壁約束徑向位移，內(nèi)壁自由。

圖4 均勻圓柱形管段示意圖Fig.4 Illustration of uniform cylindrical cavity section

橡膠基體內(nèi)圓柱形空腔的等效復(fù)波數(shù)kz的表達(dá)式如下[2]。

式中，ε=rb/ra，為圓柱形空腔的孔徑率；kl，λ，μ分別為均勻橡膠的縱波波數(shù)和拉密常數(shù)。

首先，在得到圓柱形空腔的等效復(fù)波數(shù)kz和密度后，將各種軸對(duì)稱腔形的空腔圓柱體劃分為多段，每段近似為圓柱形管，求得其等效復(fù)波數(shù)和密度。等效為均勻?qū)又?，考慮均勻圓柱體徑向位移的約束，將空腔圓柱體聲傳播狀態(tài)等同于無限大多層均勻介質(zhì)。然后，利用多層介質(zhì)的傳遞矩陣法得出入射端面的反射系數(shù)和通過結(jié)構(gòu)的透射系數(shù)，如圖5所示。圖中：pin為入射壓力；pr為反射壓力；pt為透射壓力；ρi為第i層的密度；kiz為第i層的等效軸向波數(shù)；εi為第i層的孔徑率；di為第i層的厚度，i=1,2,···,n；為出射端面的聲壓；為出射端面的法向振速；Z0為入射端特性阻抗;Zn+1為出射端特性阻抗。

將吸聲覆蓋層單元分層后等效為多層介質(zhì)的過程如圖5（a）所示，等效為多層均質(zhì)結(jié)構(gòu)后用傳遞矩陣法求解聲傳播系數(shù)的過程如圖5（b）所示。相關(guān)參數(shù)的計(jì)算公式如下：

圖5 漸變腔等效為多層介質(zhì)及傳遞矩陣法示意圖Fig.5 Illustration of transition cavity equivalent to multilayer medium and transfer matrix method

式中： ρrubber為 均勻橡膠密度；ρa(bǔ)ir為空氣密度；f為入射波頻率；rbi為第i層空腔內(nèi)徑；ci為第i層等效聲速。

在聲波垂直入射情況下，水下吸聲覆蓋層多層介質(zhì)的相鄰兩層介質(zhì)的聲壓與法向速度傳遞關(guān)系如下：

式中：p(zi)為第i端面的聲壓；v(zi)為第i端面的法向振速；Zi為第i層等效特性阻抗，Zi=ρici。

將式（3）寫為如下傳遞矩陣形式：

根據(jù)各層介質(zhì)界面處的壓力和振速連續(xù)條件，將傳遞矩陣?yán)鄢?，得到總傳遞矩陣形式如下：

由此，可以分別求出等效后多層均質(zhì)結(jié)構(gòu)的輸入面阻抗Zin和復(fù)反射系數(shù)R，以及復(fù)聲壓透射系數(shù)tp：

在求得復(fù)聲壓透射系數(shù)后，由下式獲得聲強(qiáng)透射系數(shù)tI：

由式（8）可以計(jì)算吸聲覆蓋層的吸聲系數(shù)α：

水下吸聲覆蓋層通常敷設(shè)在艇體外殼板表面，外殼板鋼材料與橡膠、水介質(zhì)的阻抗相差極大，若將鋼板考慮為鋼性背襯，即Zn+1→∞，透射系數(shù)為0，則式（8）可以改寫為

對(duì)于靜壓和空腔內(nèi)壓的作用，僅考慮其對(duì)幾何的影響，利用有限元軟件計(jì)算形變和位移，再將形變量導(dǎo)入一維理論模型計(jì)算吸聲系數(shù)。在2.1節(jié)，本文運(yùn)用有限元計(jì)算得到靜壓下吸聲覆蓋層的形變后，將變截面含圓柱形空腔吸聲覆蓋層單元?jiǎng)澐譃閚層，每層前、后界面的內(nèi)徑分別記為rbi,f和rbi,b，則形變前每段圓柱形空腔的平均內(nèi)半徑為

因?yàn)殪o壓下的水下吸聲覆蓋層圓柱形空腔內(nèi)部被壓縮，其空腔單元結(jié)構(gòu)的內(nèi)半徑和厚度都會(huì)產(chǎn)生變化，所以受壓后每段圓柱形空腔單元結(jié)構(gòu)的平均內(nèi)徑為

式中， Δrbi,f和 Δrbi,b分別為第i層單元結(jié)構(gòu)前、后界面處內(nèi)徑節(jié)點(diǎn)的徑向位移。考慮到覆蓋層單元結(jié)構(gòu)的外邊界滿足徑向約束條件，即外壁節(jié)點(diǎn)的徑向位移為0，因此，本文以新的內(nèi)徑求得孔徑率，得到每層圓柱形空腔的等效波數(shù)和平均阻抗。此外，每層圓柱形空腔前、后界面的軸向坐標(biāo)分別為di,f和di,b，受壓后第i層圓柱形空腔單元結(jié)構(gòu)的厚度為

式中： Δdi,f和 Δdi,b分別為第i層圓柱形空腔單元結(jié)構(gòu)前、后界面處的軸向位移；di,f和di,b分別為每段圓柱形空腔單元結(jié)構(gòu)的前、后界面處的軸向坐標(biāo)。將靜壓下新求得的吸聲覆蓋層內(nèi)圓柱形空腔單元結(jié)構(gòu)平均內(nèi)徑rˉbi和每段第i層圓柱形空腔單元結(jié)構(gòu)的厚度dˉi代入式（2）～式（4），最終，可以求得靜壓下吸聲覆蓋層的吸聲系數(shù)。

2 有限元計(jì)算

2.1 靜壓形變

鑒于靜壓下圓柱形空腔內(nèi)的空氣與橡膠基體之間相互作用的問題較難解決，本文將空腔內(nèi)的空氣簡化為作用在空腔壁面上的壓力ΔP。此時(shí)，空腔受壓前、后腔內(nèi)的空氣與壓力之間的關(guān)系式表示如下：

式中：p0，p分別為壓縮前、后空腔內(nèi)的壓強(qiáng)，其中p0=1.01×105Pa；ρ0，ρ分別為壓縮前、后腔內(nèi)的空氣密度；V0，Vcavity分別為壓縮前、后空腔的體積；γ為氣體常數(shù)，γ=1.4。

在橡膠基體表面受到靜壓作用時(shí)，空腔內(nèi)部的空氣會(huì)受到擠壓，產(chǎn)生抵抗空腔壁面上的壓力ΔP，并作用于腔體表面。由式（14），可得到抵抗空腔壁面上的壓力ΔP。

在使用COMSOL有限元軟件進(jìn)行靜力學(xué)計(jì)算時(shí)，對(duì)于未定義物理場的空腔域，因不能計(jì)算其體積分，而采用常規(guī)方法又無法計(jì)算空腔體積，故本文使用了Stokes公式將體積轉(zhuǎn)化為邊界面上的面積分：

式中：?為梯度算子；ξ為某一特定函數(shù)；S為待積分域V的邊界面。

采用下式計(jì)算壓縮后的空腔體積：

式中：z為空腔內(nèi)點(diǎn)的z向坐標(biāo)；nz為z向坐標(biāo)的法向矢量；Scavity為空腔內(nèi)壁面。求得壓縮后的空腔體積后，在有限元軟件COMSOL中可定義ΔP變量，作為作用于空腔壁面上的外載荷。圖6所示為有限元軟件COMSOL中計(jì)算吸聲覆蓋層單元靜壓下的形變示意圖。

圖6 運(yùn)用COMSOL計(jì)算靜壓下吸聲覆蓋層單元形變示意圖Fig.6 COMSOL calculation of deformation of anechoic coating unit cell under hydrostatic pressure

因單胞在靜壓和腔壓的作用下發(fā)生形變后仍然為回轉(zhuǎn)體，故后續(xù)吸聲計(jì)算中繼續(xù)使用二維軸對(duì)稱模式進(jìn)行分析。

2.2 聲學(xué)有限元

將含周期性空腔結(jié)構(gòu)的吸聲覆蓋層分解為多邊形單胞，等效為圓柱形單胞，并在有限元軟件COMSOL中可選用三維單元進(jìn)行建模。然而，三維模型隨著入射波頻率的增高，網(wǎng)格數(shù)量也會(huì)迅速增大，使得節(jié)點(diǎn)數(shù)量多，計(jì)算效率不高。因此，考慮到將吸聲覆蓋層單元等效為圓柱形單元后，在垂直入射波的作用下具有軸對(duì)稱性，可將其簡化為二維軸對(duì)稱單元進(jìn)行建模。

圖7所示為在有限元軟件COMSOL中建立的典型吸聲覆蓋層軸對(duì)稱模型。該模型共分為4個(gè)域，其中：Ⅰ，Ⅲ和Ⅳ為壓力聲學(xué)域，Ⅱ?yàn)楣腆w力學(xué)域；邊界1，2，3，4為軸對(duì)稱邊界，6，7，8為硬聲場邊界。為了模擬無限大阻抗背襯，邊界9約束徑向和軸向位移，邊界5約束徑向位移。在壓力聲學(xué)域Ⅲ中定義入射聲場，Ⅳ定義為完美匹配層（PML），模擬無反射邊界。在固體力學(xué)域（Ⅱ）與壓力聲學(xué)域（Ⅰ，Ⅲ）交界處設(shè)置聲?固耦合界面，以便處理和考慮水?橡膠基體及空腔?橡膠基體的耦合作用。

圖7 吸聲覆蓋層單元的三維和二維軸對(duì)稱有限元模型Fig.7 Three and two-dimensional axi-symmetric finite element models of anechoic coating unit cell

當(dāng)入射聲場Ⅲ中的聲波沿對(duì)稱軸向空腔?橡膠基體結(jié)構(gòu)垂直入射時(shí)，將產(chǎn)生聲波反射和透射，而結(jié)構(gòu)因橡膠材料的高阻尼特性會(huì)產(chǎn)生吸聲效應(yīng)。考慮到阻抗無限大背襯下沒有透射聲波，可以通過在壓力聲學(xué)域Ⅲ、固定力學(xué)域Ⅱ交界面處定義平均值算子，求得反射系數(shù)后，再求出吸聲覆蓋層結(jié)構(gòu)的吸聲系數(shù)。

本文在COMSOL有限元軟件中建立模型時(shí)采用的尺寸和材料數(shù)據(jù)參考了文獻(xiàn)[13]，將圓柱形空腔內(nèi)徑取為rb=4 mm。文獻(xiàn)[13]中采用的是三維建模，在邊界使用了周期性邊界條件來模擬無限大周期性排列的吸聲覆蓋層結(jié)構(gòu)，而本文使用的是三維建模及二維軸對(duì)稱建立的吸聲覆蓋層單胞等效模型，圖8所示即為這兩種模型模擬的結(jié)果與文獻(xiàn)[13]的結(jié)果關(guān)于吸聲覆蓋層吸聲系數(shù)的對(duì)比。

圖8 兩種模型與文獻(xiàn)[13]關(guān)于吸聲覆蓋層單元吸聲系數(shù)對(duì)比Fig.8 Sound absorption coefficient comparison of anechoic coating unit calculated by the two models with ref.[13]

由圖8可見，本文采用兩種建模方法建立的模型得到的吸聲曲線與文獻(xiàn)[13]的結(jié)果完全重合。其中，采用三維建模方法建立的模型在計(jì)算200個(gè)頻率點(diǎn)時(shí)需要10 min，而采用二維軸對(duì)稱建模方法建立的模型在計(jì)算200個(gè)頻率點(diǎn)時(shí)僅需11 s，極大地提升了計(jì)算速率，故在后續(xù)有限元計(jì)算時(shí)本文采用了二維軸對(duì)稱單元建模。

3 結(jié)果討論

模擬采用的計(jì)算模型為含圓柱形空腔的水下吸聲覆蓋層單胞結(jié)構(gòu)，如圖9所示，其空腔中心與覆蓋層幾何中心重合。表1和表2分別給出了該吸聲覆蓋層及其空腔尺寸和材料參數(shù)[8]。

表1 吸聲覆蓋層單元幾何尺寸Table 1 Geometric dimensions of anechoic coating unit cell

表2 吸聲覆蓋層橡膠材料參數(shù)Table 2 Rubber material parameters of anechoic coating

3.1 無空腔內(nèi)壓

圖10所示為不考慮空腔內(nèi)壓時(shí)不同靜壓下含圓柱形空腔水下吸聲覆蓋層的形變情況。由圖可見，靜壓下的吸聲覆蓋層單胞產(chǎn)生了軸向與徑向收縮。

圖10 不考慮腔內(nèi)壓力時(shí)不同靜壓下含空腔吸聲覆蓋層的形變Fig.10 Deformation of anechoic coating layer under different hydrostatic pressures without considering cavity pressure

在計(jì)算了靜態(tài)形變后，假設(shè)計(jì)算模型中的空腔仍為圓柱形空腔。首先，將長度縮短、徑向收縮的幾何參數(shù)代入程序，計(jì)算吸聲覆蓋層的理論吸聲系數(shù)；然后，將形變后的幾何模型導(dǎo)入聲?固耦合模塊中，計(jì)算該覆蓋層的有限元吸聲系數(shù)，結(jié)果對(duì)比如圖11所示。

由圖11可見，在不考慮靜壓下橡膠材料的參數(shù)變化而僅考慮形變時(shí)，理論計(jì)算結(jié)果與有限元計(jì)算結(jié)果基本吻合。隨著靜壓的增大，吸聲曲線基本趨勢(shì)仍保持不變，在吸聲覆蓋層出現(xiàn)軸向、徑向收縮的情況下，曲線總體上向高頻移動(dòng)。與此同時(shí)，在0.2～3 kHz頻段內(nèi)，吸聲系數(shù)有所上升，且為諧振峰移動(dòng)所致。在中頻段，靜壓下吸聲覆蓋層的吸聲性能顯著降低。而在高頻段，由于靜壓下的穿孔率變小，吸聲峰頻率變高，吸聲系數(shù)峰值變大。

圖11 不同靜水壓下吸聲覆蓋層吸聲系數(shù)Fig.11 Sound absorption coefficients of anechoic coating layer under different hydrostatic pressures

此外，理論計(jì)算的吸聲曲線相對(duì)于有限元計(jì)算的吸聲曲線，在4.4和8.5 kHz附近少了兩個(gè)尖銳谷值，這是理論方法無法考慮空腔內(nèi)空氣與橡膠基體的耦合作用所致。

3.2 考慮空腔內(nèi)壓

當(dāng)考慮空腔內(nèi)壓（腔壓）時(shí)，在空腔壁面施加隨空腔體積變化的空腔壓力，重新計(jì)算吸聲覆蓋層的形變，以及重復(fù)3.1節(jié)吸聲系數(shù)計(jì)算步驟，可得到不同靜壓下考慮空腔壓縮反力時(shí)的吸聲系數(shù)，如圖12所示。

由圖12可見，靜壓下吸聲覆蓋層發(fā)生了形變，腔壓起到了抵抗靜壓下的收縮形變作用。當(dāng)靜壓為0時(shí)（圖12（a）），空腔無靜壓形變，不產(chǎn)生對(duì)應(yīng)的抵抗反力，此時(shí)，在考慮腔壓和無腔壓這兩種情況下的吸聲系數(shù)曲線出現(xiàn)重合。然而，隨著靜壓的增大，空腔受到擠壓，腔壓作用開始明顯，吸聲峰谷值向高頻移動(dòng)的幅度減小。

圖12 不同靜壓下考慮/無腔壓的吸聲覆蓋層吸聲系數(shù)對(duì)比Fig.12 Comparison of sound absorption coefficients of anechoic coating layer with and without considering cavity pressure under different hydrostatic pressures

3.3 聲腔模態(tài)

對(duì)于軸對(duì)稱型含空腔水下吸聲覆蓋層的大部分理論研究，都將空腔吸聲覆蓋層考慮為真空?qǐng)A柱形空腔，但僅能考慮真空腔對(duì)吸聲性能的影響，而無法考慮空腔內(nèi)空氣?橡膠基體的耦合作用。然而，運(yùn)用有限元法均能夠考慮真空腔、空氣腔這種兩種情況。圖13所示為考慮上述兩種情況下采用有限元法得到的吸聲覆蓋層吸聲系數(shù)的對(duì)比。

圖13 空腔中有/無空氣的吸聲覆蓋層吸聲系數(shù)對(duì)比Fig.13 Comparison of sound absorption coefficient of anechoic coating layer with and without air in cavity

由圖13可見，當(dāng)腔體內(nèi)部為真空時(shí)，吸聲系數(shù)曲線與空腔有空氣時(shí)的大體一致，但是在0～10 kHz內(nèi)少了兩個(gè)吸聲谷值，這兩個(gè)谷值是由空腔內(nèi)空氣?橡膠基體的耦合作用所致。由圖13中（a）和（b）兩個(gè)局部詳圖可見，在4 821 Hz聲波激勵(lì)下，腔內(nèi)聲壓分布呈現(xiàn)一階軸向簡正模式，而在8 561 Hz聲波激勵(lì)下，腔內(nèi)聲壓分布呈現(xiàn)二階軸向簡正模式。

在對(duì)應(yīng)頻率下，因激起了對(duì)應(yīng)頻率激勵(lì)下的腔體聲腔模態(tài)，產(chǎn)生了吸聲谷值。而通過理論公式能夠預(yù)報(bào)得出在空腔內(nèi)空氣?橡膠基體的耦合作用下形成的谷值頻率。如式（17）所示，對(duì)于lb長度的鋼性空腔，其z軸向聲腔模態(tài)固有頻率為

式中：c0為空氣聲速；pz為空腔軸向半波數(shù)。

利用式（17），可以計(jì)算出長度為lb的鋼性空腔z軸向聲腔模態(tài)頻率。靜壓下空腔軸向發(fā)生壓縮，聲腔模態(tài)頻率升高。由模型中提取靜壓下形變后的空腔平均長度，可得到軸向聲腔一階和二階模態(tài)頻率（分別記為f10和f20），并在對(duì)應(yīng)的吸聲曲線上找出對(duì)應(yīng)的吸聲第1谷值和第2谷值頻率，如表3所示。

表3 不同靜壓下腔體聲腔模態(tài)頻率和吸聲第1及第2谷值頻率Table 3 The cavity modal frequency and the 1st and 2nd valley frequencies of sound absorption under different hydrostatic pressures

由表3可見，吸聲系數(shù)曲線上的吸聲第1和第2谷值與相應(yīng)長度腔體的聲腔模態(tài)有明顯對(duì)應(yīng)關(guān)系。在靜壓為0 MPa時(shí)，吸聲第1谷值頻率為4 281 Hz，圓柱形空腔一階軸向聲腔模態(tài)頻率（f10）為4 287.50 Hz，而吸聲第2谷值頻率為8 561 Hz，圓柱形空腔二階軸向聲腔模態(tài)頻率（f20）為8 575.00 Hz，二者幾乎相等。然而，隨著靜壓的增大，圓柱形空腔均發(fā)生徑向、軸向收縮，產(chǎn)生了不規(guī)則形變，不再是規(guī)則的圓柱體，因此，基于軸向平均形變的鋼性圓柱體的聲腔模態(tài)頻率計(jì)算來預(yù)測谷值頻率不再準(zhǔn)確。在靜壓為2.5 MPa時(shí)，吸聲第1谷值頻率為4 478 Hz，而理論分析所得到的一階軸向聲腔模態(tài)頻率為4 606.64 Hz。

4 結(jié) 論

本文分別采用理論解析與有限元法計(jì)算和分析了靜壓下含空腔水下吸聲覆蓋層的吸聲性能，著重考慮了靜壓下覆蓋層形變對(duì)吸聲系數(shù)的影響，計(jì)算驗(yàn)證了一維理論對(duì)于計(jì)算靜壓下吸聲覆蓋層吸聲性能的適用性，為快速預(yù)報(bào)吸聲性能提供了可靠的路徑?；谘芯亢陀?jì)算，主要得到如下結(jié)論：

1） 若不考慮橡膠材料參數(shù)變化，靜壓下吸聲覆蓋層的空腔會(huì)發(fā)生軸向、徑向壓縮形變。而隨著靜壓的增加，吸聲系數(shù)曲線總體向高頻移動(dòng)，低頻吸聲系數(shù)略有上升，中頻吸聲性能明顯降低。

2） 靜壓下空腔發(fā)生壓縮形變，腔內(nèi)空氣在形變下會(huì)產(chǎn)生腔壓，以抵抗水壓作用下的收縮形變。靜壓較小時(shí)，對(duì)于考慮腔壓和無腔壓這兩種情況下得到的吸聲系數(shù)曲線，二者基本重合。但是，隨著靜壓的增大，腔壓作用開始明顯，吸聲系數(shù)曲線向高頻移動(dòng)的幅度減小。

3） 目前的理論方法解析計(jì)算吸聲覆蓋層的吸聲系數(shù)大多是將含空腔水下吸聲覆蓋層分解為真空?qǐng)A柱形空腔，無法考慮空腔內(nèi)空氣?橡膠基體的耦合作用，而有限元法則可以。結(jié)果表明，考慮上述耦合作用時(shí)，吸聲系數(shù)曲線上多出了若干尖銳谷值，此處即為空腔的聲腔模態(tài)頻率。

綜上所述，本文僅考慮了靜壓下水下吸聲覆蓋層的壓縮形變，未能考慮靜壓下材料參數(shù)變化對(duì)吸聲效果的影響。若后續(xù)研究有靜壓下橡膠的靜態(tài)參數(shù)與動(dòng)態(tài)實(shí)驗(yàn)參數(shù)，按照本文方法，可以分別將其引入到形變計(jì)算和吸聲計(jì)算中，從而獲得更接近于實(shí)際的吸聲覆蓋層吸聲性能。