基于聲管測量的粘彈性介質(zhì)動態(tài)力學(xué)參數(shù)辨識方法

陶 猛 , 趙 陽

(1.貴州大學(xué) 機械工程學(xué)院，貴陽 550025；2.上海交通大學(xué) 機械系統(tǒng)與振動國家重點實驗室,上海 200240)

以粘彈性介質(zhì)為基底材料的聲學(xué)覆蓋層既可以主動聲吶的探測回波，還具有隔振隔聲的作用以降低噪聲的向外輻射，作為提高結(jié)構(gòu)隱身性能的一種重要技術(shù)手段已被廣泛采用。聲學(xué)覆蓋層的吸聲、隔聲等性能預(yù)報除了要采用合理的計算分析模型外，準(zhǔn)確的粘彈性動態(tài)力學(xué)參數(shù)也是必不可少的。

通常測定粘彈性材料參數(shù)的方法可分為力學(xué)方法和聲學(xué)方法兩大類。力學(xué)方法是利用測量材料樣品振動特性來計算其力學(xué)參數(shù)的方法，比較常用的有強迫非共振法[1-2]、振動梁法[3-4]、動態(tài)粘彈譜儀法[5-6]。盡管這些測試方法都相對簡單，但是這些方法都有各自的局限性，例如振動測試方法的頻段較低，或者基于時溫等效原理推導(dǎo)的楊氏模量還沒有直接的測試結(jié)果加以驗證等。與力學(xué)測試方法相比，聲學(xué)測試方法的優(yōu)點是直接測得材料的聲學(xué)性能參數(shù)，這樣可在頻段范圍上與機理研究同步。已有方法包括在水池中由待測材料的平板試樣測量斜向入射聲波的回聲降低或插入損失來反演材料參數(shù)[7-8]，但低頻測量時由于樣品邊緣的衍射干擾使得誤差較大。此外，也有利用做成球形的待測材料試樣，使實際測量的散射系數(shù)與理論計算的散射系數(shù)誤差達到最小的辦法來計算材料體積模量的做法[9]。

本文研究了一種基于水聲聲管測量粘彈性覆蓋層的反射系數(shù)，以此識別粘彈性材料動態(tài)力學(xué)參數(shù)的方法。該方法的原理是考慮平面波垂直入射到圓柱空腔覆蓋層時，覆蓋層單元-粘彈性圓柱管中傳播的軸對稱波與基底材料的縱波和剪切波滿足粘彈性圓柱管的特征方程。若其中兩個參數(shù)已知，求解特征方程即可獲得未知的第三個參數(shù)，并進而得到其它的粘彈性動態(tài)力學(xué)參數(shù)。然后，以某種橡膠制成的聲學(xué)覆蓋層樣品為例，進行了聲管測量和結(jié)果分析。

1 粘彈性材料的動態(tài)力學(xué)參數(shù)辨識方法

不論是均勻結(jié)構(gòu)的還是空腔結(jié)構(gòu)的、厚度為d的聲學(xué)覆蓋層，都可用四端網(wǎng)絡(luò)法建立覆蓋層前表面(下標(biāo)f表示)和后表面(下標(biāo)b表示)的聲壓p和質(zhì)點振速v的關(guān)系：

(1)

式中：對于均勻覆蓋層來說，z和k分別是特性阻抗zc=ρω/kl和縱波波數(shù)kl=ω/cl，ρ是材料密度；對于含有圓柱空腔結(jié)構(gòu)覆蓋層來說，z和k可以認(rèn)為是結(jié)構(gòu)等效阻抗ze=ρω/ke[10]和等效波數(shù)ke=ω/ce。

若波數(shù)kl和ke已知，可通過以下1.1和1.2節(jié)的參數(shù)辨識模型獲得粘彈性介質(zhì)的動態(tài)力學(xué)參數(shù)。因此，kl和ke是參數(shù)辨識過程中的關(guān)鍵中間變量，其獲取方法是通過測量聲學(xué)覆蓋層的反射系數(shù)并計算得到。

當(dāng)平面波正入射到聲學(xué)覆蓋層時，反射系數(shù)R可由下式計算：

(2)

式中：zf和z0分別是覆蓋層前表面阻抗和水介質(zhì)特性阻抗。

根據(jù)式(1)，在后表面阻抗(背襯條件)zb=pb/vb已知的前提下(如絕對硬背襯zb=∞)，前表面阻抗zf可由下式得到

(3)

顯然，對于絕對硬背襯，zf=t11/t21；對于絕對軟背襯，zf=t12/t22。由于t12和t21既是波數(shù)k的函數(shù)，也是阻抗z=ρω/k的函數(shù)，所以當(dāng)實驗測得覆蓋層的反射系數(shù)后，聯(lián)立求解式(2)和式(3)，即可得到波數(shù)kl(對于均勻結(jié)構(gòu))或ke(對于空腔結(jié)構(gòu))。需要指出的是，由于粘彈性介質(zhì)通常用復(fù)數(shù)形式表示其損耗特性，因而上述過程是通過在復(fù)平面上解算超越方程實現(xiàn)的。

1.1 復(fù)縱波聲速計算

(4)

式中：cl和ηcl是復(fù)縱波聲速的實部和損耗因子。需要指出的是，對于多數(shù)橡膠材料而言，盡管對縱波損耗因子來說很小，但是并不應(yīng)該忽略不計。這是因為在計算復(fù)剪切波聲速的過程中，復(fù)縱波聲速是作為已知參數(shù)代入計算的，如果忽略復(fù)縱波聲速的損耗因子，必然會在復(fù)剪切波聲速的計算中引入誤差。

1.2 復(fù)剪切波聲速計算

對于圓柱空腔按照正六邊形周期排列的吸聲覆蓋層，由于對稱性可以只取其中的一個單元分析，但是建立正六面棱柱體的理論模型相當(dāng)困難，用圓柱代替可以簡化模型的復(fù)雜性，這樣就得到粘彈性圓柱管的模型[11]，如圖1所示，其中粘彈性圓柱管的外半徑為a，內(nèi)半徑為b。

圖1 周期結(jié)構(gòu)圓柱空腔型吸聲覆蓋層的結(jié)構(gòu)示意圖

假設(shè)在聲波小振幅激勵的情況下，采用Kelvin-Voigt線性粘彈性模型來描述覆蓋層的基體材料，可以得到與彈性模型形式相同的方程和解，區(qū)別在于將彈性常數(shù)換成相應(yīng)的粘彈性常數(shù)[12]。圓柱管中軸對稱波的形式解一般可以寫成[11](省略時間因子ejωt)：

(5)

在圓柱坐標(biāo)系下，位移和應(yīng)力用兩個標(biāo)量勢函數(shù)可表示為:

(6)

將式(5)代入式(6)得到:

(7)

其中

(8)

對于按照正三角形周期排列的圓柱空腔覆蓋層，在平面波垂直入射條件下，根據(jù)相鄰單元平衡的條件取單元的外邊界(r=a)位移等于0，切應(yīng)力等于0；內(nèi)壁的邊界條件忽略空氣的影響，取正應(yīng)力和切應(yīng)力均為0。因此，內(nèi)、外壁邊界條件可寫為：(省略e-jkaz)

(9)

將位移和應(yīng)力的表達式代入邊界條件，得到方程組：

(10)

這就是粘彈性圓柱管的特征方程。要使該方程有解，則系數(shù)行列式必須為零，即:

(11)

(12)

(13)

2 動態(tài)力學(xué)參數(shù)識別的結(jié)果及討論

在水聲聲管中對以某種橡膠材料制作的實心覆蓋層和空心覆蓋層兩種樣品進行了反射系數(shù)的測量，測量的方法采用成熟的傳遞函數(shù)法[13]，根據(jù)測試結(jié)果對橡膠材料的粘彈性動態(tài)力學(xué)參數(shù)進行了計算，其中橡膠密度為1 120 kg/m3，圓柱空腔半徑為2 mm。

圖2 橡膠材料的復(fù)縱波聲速和復(fù)剪切波聲速

圖3 橡膠材料的復(fù)楊氏模量

此外，表1還比較了是否考慮縱波損耗因子對楊氏模量及其損耗因子計算結(jié)果的影響，其中定義了楊氏模量及其損耗因子的相對誤差

(14)

式中E和ηE分別是考慮縱波損耗因子的楊氏模量及其損耗因子，E*和ηE*分別是忽略縱波損耗因子的楊氏模量及其損耗因子。從表中可以看出，當(dāng)忽略較大的縱波損耗因子時，楊氏模量損耗因子計算的相對誤差也較大，如頻率1.4 kHz、ηcl=0.024時，相對誤差εηE=3.89%。需要指出的是，這里只是分析了忽略縱波損耗因子對復(fù)楊氏模量的影響關(guān)系，而對反射系數(shù)測試過程中引入的誤差環(huán)節(jié)暫不涉及。另一方面，通過比較得知，忽略縱波損耗因子對楊氏模量損耗因子的影響要大于對楊氏模量實部的影響。

表1 楊氏模量及其損耗因子計算的相對誤差

圖4 橡膠材料的復(fù)泊松比

3 結(jié) 論

本文研究了一種基于聲管測量粘彈性覆蓋層反射系數(shù)，并采用解析法計算和識別粘彈性介質(zhì)動態(tài)力學(xué)參數(shù)的方法，這種方法的特點是能夠與聲學(xué)性能測試同步進行，便于聲學(xué)覆蓋層的性能預(yù)報和優(yōu)化設(shè)計。通過對某種橡膠材料的測試分析可知，盡管橡膠類粘彈性材料的縱波損耗因子較小，但是仍對剪切波損耗因子的準(zhǔn)確識別有一定的影響，因此縱波損耗因子不應(yīng)簡單地忽略不計。最后需要指出的是，由于測試聲管的截止頻率限制，文中對于材料參數(shù)的分析局限于中頻范圍，而本文提出的方法原則上也適用于高頻情況，如可采用脈沖法測量覆蓋層的反射系數(shù)，進而拓寬粘彈性介質(zhì)動態(tài)力學(xué)參數(shù)的識別范圍。

參 考 文 獻

[1]Caracciolo R, Gasparetto A, et al. An experimental technique for complete dynamic characterization of a viscoelastic material[J]. Journal of Sound and Vibration, 2004, 272 (3): 1013-1032.

[2]Boiko A V, Kulik V M, et al. Measurement method of complex viscoelastic material properties[J]. International Journal of Solids and Structures, 2010, 47 (3): 374-382.

[3]Liao Y, Wells V. Estimation of complex modulus using wave coefficients[J]. Journal of Sound and Vibration, 2006, 295 (1-2):165-193.

[4]Liao Y, Wells V. Estimation of complex Young’s modulus of non-stiff materials using a modified Oberst beam technique[J]. Journal of Sound and Vibration, 2008, 316 (1-5):87-100.

[5]Li R Z. Time-temperature superposition method for glass transition temperature of plastic materials[J]. Materials Science and Engineering A, 1999, 278 (1-2): 36-45.

[6]Liu H W, Yao L, et al. A mixed method for measuring low-frequency acoustic properties of macromolecular materials[J]. Science in China Series G: Physics, Mechanics and Astronomy, 2006, 49(6): 729-739.

[7]Piquette J C. Shear property determination from underwater acoustic panel tests[J].Journal of the Acoustical Society of America, 2004, 115 (5): 2110-2121.

[8]Gartland G J, Radcliffeb C J, et al. Measurement of dilatational wave speed using an echo reduction test[J]. Journal of Sound and Vibration, 2009, 320 (3): 491-495.

[9]宋 揚. 中高頻段下粘彈性材料聲學(xué)參數(shù)測量[D].哈爾濱:哈爾濱工程大學(xué), 2007.

[10]Tao Meng. Simplified model for predicting acoustic performance of an underwater sound absorption coating[J]. Journal of Vibration and Control, doi: 10.1177/1077546312461027.

[11]湯渭霖,何世平,范 軍.含圓柱形空腔吸聲覆蓋層的二維理論[J].聲學(xué)學(xué)報,2005, 30 (4): 289-295.

TANG Wei-lin, HE Shi-ping, FAN Jun. Two-dimensional model for acoustic absorption of viscoelastic coating containing cylindrical holes[J]. Shengxue Xuebao, 2005, 30 (4): 289-295.

[12]Gaunaurd C. Methods for solving the viscoelasticity equations for cylinder and sphere problems[R]. NSWC/WOL Tech. Rep., ADA025302, 1976.

[13]朱蓓麗,肖今新. 雙水聽器傳遞函數(shù)法低頻測試及誤差分析[J]. 聲學(xué)學(xué)報, 1994, 19 (5): 351-360.

ZHU Bei-li, XIAO Jin-xin. A two-hydrophone transfer function method for measuring low-frequency acoustic properties and its error analysis[J]. Shengxue Xuebao, 1994, 19 (5): 351-360.

[14]Pritz T. The Poisson’s loss factor of solid viscoelastic materials[J]. Journal of Sound and Vibration, 2007, 306 (3): 790-802.