閆兆敏，張東俊, 溫激鴻，溫熙森

(1.聲學與振動研究組 裝備綜合保障技術重點實驗室 國防科技大學,長沙 410073；2.海軍92337部隊，大連 116023)

薄板及由其組成的復雜結構憑借良好的工程性質(zhì)，在機械、船舶與航空等領域得到了廣泛的應用，工程中的許多結構都是由薄板的典型結構組成，各國學者很早就開展了對其振動與聲輻射問題的研究。阻尼減振降噪技術(簡稱阻尼技術)是控制結構振動和聲輻射的一項重要技術，其原理是在結構表面貼敷阻尼材料，將結構的振動能轉(zhuǎn)化為熱能耗散從而達到減小振動和聲輻射的目的，在實際工程中經(jīng)常通過鋪設多層性質(zhì)不同的阻尼材料進一步增強阻尼減振降噪的效果[1]。

鋪設阻尼材料的薄板結構可作為復合材料層合板結構來進行研究。根據(jù)復合材料層合結構的特點，對層合板動態(tài)性能的求解常常需要利用一些物理假設來構建位移方程。根據(jù)位移方程的不同，層合板理論主要可分為經(jīng)典層合板理論(Classical Lamination Theory, CLD)[2-3]、一階剪切變形理論(First-Order Shear Deformation Theory, FSDT)[4]、高階剪切變形理論(Higher-Order Shear Deformation Theory, HSDT)[5-6]、分層理論(Layerwise Theory)[7-9]和三維彈性理論(3-D Elasticity Theory)[10]，Khandan等[11]對層合板理論的發(fā)展歷程進行了綜述。分層理論可看作是一種準三維彈性理論，其在板的厚度方向利用插值函數(shù)來描述每層內(nèi)位移沿厚度方向的變化規(guī)律，結合有限元法可憑借較小的計算開銷取得較好的計算精度，因而具有廣泛的適用性。

目前，許多學者已對貼敷單層阻尼的結構振動與聲輻射問題開展了廣泛的研究[12-15]，而對貼敷多層阻尼結構以及流動流體中貼敷阻尼結構的振動和聲輻射問題研究較少,與此同時，分層理論主要應用于單層板殼結構的振動特性分析[7-9]，將其應用于多層板殼結構聲輻射分析的研究也不常見。本文采用分層理論建立了四邊簡支約束條件下薄板聲輻射特性有限元模型，計算了其固有頻率、輻射聲功率以及振動均方速度，并與Comsol軟件計算結果進行對比，證明了該方法的有效性。其次，將分層理論拓展運用于貼敷阻尼薄板的聲輻射性能研究中，分析了單/多層阻尼對結構聲輻射的影響規(guī)律。最后通過分析貼敷阻尼薄板在輕/重流體、不同流速條件下的聲振特性，研究流體對貼敷阻尼薄板聲輻射性能的影響，為薄板結構貼敷阻尼以降低聲輻射水平奠定了理論基礎。

1 理論模型

1.1 基于分層理論的有限元模型

薄板貼敷多層黏彈性阻尼材料形成的層合板如圖1所示。

根據(jù)分層理論，可將各層阻尼沿厚度方向進行人為插值。以第k阻尼層為例，其沿x,y,z方向的位移可表示為[8]：

(1)

(2)

(3)

式中:n表示平面的節(jié)點數(shù)，Ni表示i節(jié)點的平面插值形函數(shù)。將式(3)代入式(1)可得：

(4)

由式(4)可以發(fā)現(xiàn)，分層理論將三維結構的位移利用二維平面的位移來表示的同時也考慮了厚度方向各層之間的連接關系。

1.2 結構聲輻射模型

圖2 聲輻射模型示意圖

聲輻射模型如圖2所示，有關流體聲學負載的計算，文獻[16-17]給出了詳細的推導過程，結構表面上一點M的聲壓可表示為

(5)

(-ω2M+iω(C+GAZGT)+K)w=f

(6)

值得注意的是，當考慮流體的流速時，式(5)只在較低的頻段內(nèi)有效[16]。

當求得結構表面的聲壓和法向振速后，結構的輻射聲功率可通過下式求得

(7)

式中:M是屬于面S上的一點，Re表示取實部，p(M)表示M點處的聲壓，un(M)表示M點處的法向振速，“*”表示復數(shù)共軛。

當結構發(fā)生振動時，結構會向周圍流體介質(zhì)進行壓縮產(chǎn)生聲波，聲波通過流體介質(zhì)向遠場輻射，形成聲輻射。聲輻射阻尼用以表征流體將振動能轉(zhuǎn)化為聲能的能力，其定義為[19]：

μf=ER/ES=WT/ES=2πW/ESω

(8)

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

2.1 數(shù)值算例

固有模態(tài)頻率是結構的重要參數(shù)，為驗證前述分層模型的準確性，構建四邊簡支約束矩形薄板模型，其物理參數(shù)見表1。

表1 薄板物理參數(shù)

通過前述分層模型理論求得的結構前5階固有頻率與Comsol軟件及文獻結果對比如表2所示：

表2 四邊簡支層合板結構前5階固有頻率對比(單位：Hz)

在四邊簡支邊界條件下，基于分層理論的有限元法與商用有限元軟件及文獻計算結果吻合度較高。此外，若沿Z方向施加簡諧激勵力F施加于(0.113 75,0.094 75,0)處，幅值為1 N，激勵頻率范圍為10-1 000 Hz，分別利用分層理論有限元法與Comsol計算在空氣介質(zhì)中(空氣的密度為1.225 kg/m3，聲速為344 m/s)的薄板輻射聲功率與振動均方速度結果對比如圖3所示，由結果可以發(fā)現(xiàn)，分層理論所得結果與商用有限元軟件結果較吻合，從而證明了利用該方法能夠較準確地計算出結構的聲振特性。

圖3 結構聲輻射特性對比

2.2 貼敷單/多層阻尼對薄板聲輻射的影響

以前述薄板結構為研究對象，在其表面貼敷阻尼材料以抑制其結構聲輻射，阻尼材料特性如表3所示。

表3 阻尼層材料物理參數(shù)

將厚度為6 mm的阻尼材料1均勻貼敷在薄板結構表面，結構周圍為空氣介質(zhì)，計算得到薄板結構的聲輻射特性如圖4所示。

從圖中可以發(fā)現(xiàn)，貼敷阻尼材料后，結構的等效剛度、等效密度以及等效阻尼等參數(shù)發(fā)生變化，使得結構的輻射聲功率級峰值頻率向高頻方向移動，在10-1 00 Hz范圍內(nèi)的輻射聲功率和振動均方速度均下降了約10 dB，效果顯著，由此說明貼敷阻尼材料對抑制薄板振動和聲輻射具有較明顯的作用。

進一步對表面貼敷不同種類阻尼組合而成的阻尼材料的結構聲輻射特性進行探討。以兩層阻尼結構為例，將表3中的兩種阻尼材料進行組合，構成混合阻尼貼附在前述薄板結構表面，每種阻尼層的厚度均為3 mm，此時共有2種貼敷方式，沿Z軸正方向排序，可將這2種貼敷方式按材料種類依次表示為1-2、2-1，外部激勵F與前述相同,計算結果如圖5所示。

由圖5可以看出，貼敷兩層阻尼的薄板結構聲輻射特性對阻尼材料的貼敷順序較為敏感。以第1階輻射峰為例，2-1貼敷方式下的輻射聲功率級和振動均方速度與1-2貼敷方式相比均降低了約20 dB，效果顯著，說明將彈性模量較大的阻尼材料貼敷在結構靠外側更有利于減小結構的振動和聲輻射。這是因為將彈性模量較大的阻尼材料貼敷在靠外側可以與內(nèi)側彈性模量較小的阻尼材料構成約束阻尼，當結構受到振動而發(fā)生彎曲變形時，內(nèi)側模量較小的阻尼層受上下兩個約束面的作用而不能有伸縮變形，各層之間發(fā)生剪切變形耗能，同時，外層阻尼還可作為自由阻尼層，利用交變的應力和應變來消耗振動能量，這樣可提高結構整體的損耗因子，進一步減小結構的振動，從而抑制聲輻射。

2.3 流體對貼敷阻尼層薄板聲輻射的影響

以2.2節(jié)所述貼敷單層阻尼薄板為研究對象，分別討論結構在輕、重流體中的聲振特性(本文以空氣作為輕流體，水作為重流體)，結果如圖6所示，由圖中可以發(fā)現(xiàn)流體對結構聲輻射性能有重要的影響。相比于輕流體，結構在重流體中的模態(tài)頻率向低頻方向移動，振動均方速度顯著降低，而輻射聲功率則明顯升高，這是由于流固耦合作用對結構引入了附加質(zhì)量，而重流體所引入的附加質(zhì)量遠遠高于輕流體，因此出現(xiàn)結構在重流體中的模態(tài)頻率向低頻方向移動，且振動受到明顯抑制。此外，重流體的密度和聲速較高，使得重流體的輻射聲阻抗較高，在較小振動的條件下能夠引起較高的聲壓，從而使得其輻射聲功率高于輕流體。

圖6 輕/重流體介質(zhì)中結構聲輻射特性

比較輕/重流體中的結構聲輻射情況還可以發(fā)現(xiàn)，在重流體中的輻射聲功率以及振動均方速度曲線更為平緩，結構阻尼效果更加明顯。通過比較兩種流體介質(zhì)的聲輻射阻尼(如圖7所示)可以發(fā)現(xiàn)，水中的聲輻射阻尼遠遠大于空氣，結構在水中所受的阻尼效果遠大于在空氣中的效果，因此會出現(xiàn)兩條曲線更為平緩的現(xiàn)象。

除流體的密度和聲速外，流體流速也是影響結構聲輻射的重要因素。前述貼敷單層阻尼薄板在不同流速重流體中的聲輻射情況如圖8所示，隨著流體流速的增大，結構的輻射峰向低頻方向移動，振動和聲輻射均有增大的趨勢，在峰值處的幅值變化更加平緩，且流速主要在非峰值頻率段的結構聲輻射以及峰值頻率段的結構振動方面產(chǎn)生較大影響，通過式(5)可知產(chǎn)生這種現(xiàn)象的主要原因是流體流動時會引入了負剛度和附加阻尼，改變結構的動態(tài)力學參數(shù)，從而使結構聲輻射性能發(fā)生變化。此外，由圖8還可以發(fā)現(xiàn)流體流速對結構聲輻射的影響是非線性的，以 10 Hz處的輻射聲功率及振動均方速度為例，其變化趨勢如圖9所示，當流速較低時，結構的聲振特性變化明顯，此后，隨著流體流速的增大，結構的聲振特性也趨于平穩(wěn)。

圖8 不同流速流體中結構聲輻射特性

3 結 論

本文基于分層理論建立了四邊簡支約束條件下貼敷阻尼薄板聲振特性有限元分析模型，計算了其固有頻率、輻射聲功率以及振動均方速度，并與Comsol軟件計算結果進行對比，證明了該方法用于分析結構聲振特性的有效性。其次，將分層理論用于貼敷阻尼薄板的聲輻射性能研究中，分析了單/多層阻尼對結構聲輻射的影響規(guī)律。最后，通過分析貼敷阻尼薄板在不同流體介質(zhì)、不同流速條件下的聲振特性，研究了流體對貼敷阻尼薄板聲輻射性能的影響。結果表明：

(1) 貼敷阻尼材料對結構等效密度、等效剛度以及等效阻尼影響顯著，通過貼敷阻尼后可使薄板在第1階峰值頻率處的輻射聲功率和振動均方速度降低約10 dB。貼敷多層阻尼的薄板結構聲輻射特性對阻尼材料的貼敷順序較為敏感，將彈性模量較大的阻尼材料貼敷在靠外側能夠增大約束阻尼層的剪切變形，有利于提高結構整體的損耗因子，降低結構的振動和聲輻射。

(2) 流固耦合會對結構產(chǎn)生附加阻尼作用，進一步加強阻尼層的貼敷效果，使得結構聲振特性變化更為平緩。結構在重流體介質(zhì)中的振動受到抑制，但由于重流體介質(zhì)中的聲輻射阻尼較大，導致結構在重流體中的輻射聲功率遠高于輕流體。

(3) 由于流體流速引入負剛度和附加阻尼作用，因而會使結構聲振特性發(fā)生較大變化且變化為非線性，當流速較低時，結構的聲振特性變化明顯，此后，隨著流體流速的增大，結構的聲振特性趨于平穩(wěn)。

[1]何作鏞. 結構振動與聲輻射[M]. 哈爾濱: 哈爾濱工程大學出版社, 2001.

[2]Petyt M. Introduction to finite element vibration analysis[M]. 2nd ed. New York: Cambridge University Press, 2010.

[3]Parthasarathy G, Reddy C V R. Partial coverage of rectangular plates by unconstrained layer damping treatments[J]. Journal of Sound and Vibration,1985, 102: 203-216.

[4]Bhar A, Phoenix S S, Satsangi S K. Finite element analysis of laminated composite stiffened plates using FSDT and HSDT: A comparative perspective[J]. Composite Structures,2010, 92: 312-321.

[5]Xiang S, Jiang S, Bi Z, et al. A nth-order meshless generalization of Reddy’s third-order shear deformation theory for the free vibration on laminated composite plates[J]. Composite Structures,2011, 93(2011): 299-307.

[6]Kumar J S, Raju T D, Reddy K V K. Vibration analysis of composite laminated plates using higher-order shear deformation theory with zig-zag[J]. Indian Journal of Science and Technology,2011, 4(8): 960-966.

[7]Desai Y M, Ramtekkar G S, Shah A H. Dynamic analysis of laminated composite plates using a layer-wise mixed finite element model[J]. Composite Structures,2003, 59(2): 237-249.

[8]Reddy J N. Mechanics of laminated composite plates and shells: theory and analysis[M]. 2nd ed. Boca Raton: CRC PRESS, 2003.

[9]胡明勇，王安穩(wěn). 纖維增強粘彈性復合材料層合板的自由振動和應力分析[J]. 工程力學,2010, 27(8): 10-20.

HU Ming-yong, WANG An-wen. Free vibration and stresses analysis of fiber-reinforced viscoelastic composite laminated plates[J]. Engineering Mechanics, 2010, 27(8):10-20.

[10]Reddy J N. Theories and computational models for composite laminate[J]. Applied mechanics reviews,1994, 47(6): 147-169.

[11]Khandan R, Noroozi S, Sewell P, et al. The development of laminated composite plate theories:a review[J]. Journal of Material Science, 2012, 47: 5901-5910.

[12]Laulagnet B, Guyader J L. Sound radiation from finite cylindrical coated shells, by means of asymptotic expansion of three-dimensional equations for coating[J]. J.Acoust.Soc.Am,1994, 96(1): 277-286.

[13]Laulagnet B, Guyader J L. Sound radiation from finite cylindrical shells partially covered with longitudinal strips of compliant layer[J]. Journal of Sound and Vibration,1995, 186(5): 723-741.

[14]王獻忠，孫龍泉，邱忠輝，等. 部分敷設阻尼材料的水下結構聲輻射分析[J]. 振動與沖擊, 2012, 31(18): 122-127.

WANG Xian-zhong, SUN Long-quan, QIU Zhong-hui, et al. Sound radiation analysis for a submerged structure with a viscoelastic damping layer partially covered[J]. Journal of Vibration and Shock,2012,31(18): 122-127.

[15]陳美霞，駱東平，周峰，等. 阻尼材料敷設方式對雙層殼體聲輻射性能的影響[J]. 聲學學報, 2005, 30(4): 296-301.

CHEN Mei-xia, LUO Dong-ping, ZHOU Feng, et al. Effect of viscoelastic layer cover type on characteristics of sound radiation from double cylindrical shell[J]. Acta Acustica, 2005, 30(4):296-301.

[16]Atalla N, Nicolas J. A Formulation for mean flow effects on sound radiation from rectangular baffled plates with arbitrary boundary conditions[J]. Journal of Vibration and Acoustics,1995, 117: 22-29.

[17]Sgard F, Atalla N, Nicolas J. Coupled FEM-BEM approach for mean flow effects on vibro-acoustic behavior of planar structures[J]. AIAA Journal,1994, 32(12): 2351-2358.

[18]姜弘道. 彈性力學問題的邊界元法[M]. 北京: 中國水利水電出版社, 2008.

[19]郭新毅，洪明. 結構振動聲輻射阻尼研究[J]. 中國造船,2004, 45: 189-196.

GUO Xin-yi, HONG Ming. Research on acoustic radiation damping by structual vibration[J]. Shipbuilding of China,2004,45:189-196.

[20]黎勝. 水下結構聲輻射和聲傳輸?shù)臄?shù)值分析及主動控制模擬研究[D]. 大連: 大連理工大學, 2001.