基于波函數(shù)法的自由阻尼薄板與聲學(xué)耦合響應(yīng)

夏小均，徐中明，賴詩洋，張志飛, 賀巖松

(1.重慶大學(xué) 機械傳動國家重點實驗室,重慶 400030；2.重慶大學(xué) 汽車工程學(xué)院,重慶 400030;3.重慶工程職業(yè)技術(shù)學(xué)院 機械工程學(xué)院,重慶 402260)

基于波函數(shù)法的自由阻尼薄板與聲學(xué)耦合響應(yīng)

夏小均1,2，徐中明1,2，賴詩洋3，張志飛1,2, 賀巖松1,2

(1.重慶大學(xué) 機械傳動國家重點實驗室,重慶 400030；2.重慶大學(xué) 汽車工程學(xué)院,重慶 400030;3.重慶工程職業(yè)技術(shù)學(xué)院 機械工程學(xué)院,重慶 402260)

對敷設(shè)自由阻尼薄板振動響應(yīng)進行了分析，以復(fù)剛度為基礎(chǔ)得到了在Kirchhoff理論下復(fù)合薄板振動控制方程。基于波函數(shù)法理論，推導(dǎo)了自由阻尼薄板振動分析模型以及包含自由阻尼結(jié)構(gòu)與聲腔的三維耦合模型的建模方法;以四邊固支自由阻尼矩形板及耦合的結(jié)構(gòu)聲學(xué)系統(tǒng)為例，分別以波函數(shù)法與有限元法計算了其50～500 Hz頻段內(nèi)的結(jié)構(gòu)與聲學(xué)響應(yīng)。結(jié)果表明：波函數(shù)法能有效的應(yīng)用于添加自由阻尼的薄板振動以及結(jié)構(gòu)聲耦合系統(tǒng)響應(yīng)的預(yù)測與分析;相比于有限元法，其高精度、高收斂率的特點使波函數(shù)能有效解決更高頻率的聲振問題。

自由阻尼；波函數(shù)法；彎曲振動；結(jié)構(gòu)聲耦合

隨著計算機輔助工程CAE(Computer Aided Engineering)技術(shù)更廣泛的運用和對現(xiàn)在有數(shù)值分析手段的深入研究，研究者發(fā)現(xiàn)廣泛應(yīng)用的基于單元的方法如有限元法、邊界元法，由于插值誤差與污染誤差的存在，其只能應(yīng)對較低頻率問題。而基于模態(tài)重疊率假設(shè)的統(tǒng)計類方法如統(tǒng)計能量法只對高頻問題有效[1-3]。因此，20世紀末出現(xiàn)的“中頻危機”使得人們努力地探尋針對中頻問題的數(shù)值分析方法和手段。當(dāng)下代表性研究主要有：波有限元法(Wave Finite Element Method,WFEM)[4]，統(tǒng)計能量分布法(Statistical Modal Energy Distribution Analysis,SMEDA)[5]，復(fù)包絡(luò)向量(Complex Envelope Vectorization,CEV)[6]，復(fù)射線變分理論(Variational Theory of Complex Rays,VTCR)[7], 混合有限元-統(tǒng)計能量法(Hybrid Finite Element/Statistical Energy Analysis Method,Hybrid FE-SEA)[8],能量分析布法(Energy Distribution Analysis,EDA)[9],以及波函數(shù)法(Wave Based Method,WBM)[10]。其中Hybrid FE-SEA已經(jīng)應(yīng)用在一些商用軟件上，指導(dǎo)了很多工程研究，但該種方法的結(jié)果仍是時間或空間上的平均，且在仍存在一些理論待完善。而本文要討論的波函數(shù)則是完全不同現(xiàn)有基于單元法的確定性數(shù)值分析方法。

到目前，波函數(shù)法已經(jīng)成功運用到了平板彎曲振動、薄膜振動、結(jié)構(gòu)聲學(xué)耦合響應(yīng)[11-13]、多孔阻尼材料聲學(xué)耦合響應(yīng)[14]等問題。國內(nèi)Peng等[16-17]對基于薄板彎曲振動以及二維的結(jié)構(gòu)聲耦合問題也進行過相應(yīng)研究。通過對該類常見聲振問題的預(yù)測分析，波函數(shù)法展現(xiàn)出了高效高精度的特點，使其能有效的運用于中頻振動與聲學(xué)問題。而在振動特性明顯的薄板類結(jié)構(gòu)上敷設(shè)阻尼，業(yè)已成為工程中最常見的減振降噪手段，雖然基于單元法對自由阻尼結(jié)構(gòu)的也有相關(guān)研究[18]，但如上文所述，該類方法分析中頻問題往往要求付出巨大的時間成本與精度妥協(xié)。因此，研究以如波函數(shù)法等高效的數(shù)值方法分析含有自由阻尼的結(jié)構(gòu)振動與聲學(xué)響應(yīng)問題，實現(xiàn)該類系統(tǒng)的中頻聲振預(yù)測具有實際的研究意義與工程價值。

文章在完成自由阻尼薄板振動分析的基礎(chǔ)上，結(jié)合波函數(shù)法理論完成該類結(jié)構(gòu)的彎曲振動響應(yīng)預(yù)測及其與聲腔耦合系統(tǒng)的聲學(xué)響應(yīng)分析。并通過數(shù)值算例，驗證了方法的有效性，同時也再次體現(xiàn)了波函數(shù)法在針對中頻聲振問題的高精度與高收斂特性。

1 自由阻尼薄板振動響應(yīng)

對于添加均勻厚度自由阻尼的薄板如圖1所示。結(jié)構(gòu)層厚度為t1,阻尼層厚度為t2。

圖1 自由阻尼薄板截面示意圖Fig.1 Section of unconstraned damped plate

彎曲剛度為薄板類結(jié)構(gòu)振動響應(yīng)的關(guān)鍵，本次將敷設(shè)阻尼后構(gòu)建虛擬約束阻尼結(jié)構(gòu)，通過截面彎矩平衡求得帶約束層的彎曲剛度，通過消去約束層以及剪切變形的作用后得到自由阻尼薄板的彎曲復(fù)剛度B*[19]

(1)

對于含有自由阻尼的薄板，以經(jīng)典的Kirchhoff理論來描述其彎曲振動。復(fù)合薄板振動方程

(2)

對于穩(wěn)態(tài)振動響應(yīng)，其位移可表達為wt(xs,ys)=wz(xs,ys)ejωt，而在外部激勵為Qt=Qejωt時，含有自由阻尼薄板的穩(wěn)態(tài)振動響應(yīng)控制方程表示為

(3)

2 結(jié)構(gòu)聲學(xué)耦合穩(wěn)態(tài)響應(yīng)

對于含自由阻尼和薄板聲腔耦合系統(tǒng)，由于受到耦合聲壓面作用，復(fù)合薄板的振動響應(yīng)方程修正為

(4)

而聲學(xué)域仍滿足Helmholtz波動方程

(5)

(6)

3 波函數(shù)法理論與實現(xiàn)

波函數(shù)法作為基于間接Trefftz方法，與現(xiàn)有基于單元的方法不同。其將分析域變量表示為嚴格滿足控制方程的一系列波函數(shù)的疊加。再通過如Galerkin等加權(quán)余量法實現(xiàn)邊界殘差積分歸零，來求得各波函數(shù)的權(quán)系數(shù)，進而得到求解響應(yīng)。

3.1含自由阻尼薄板的彎曲振動響應(yīng)

根據(jù)波函數(shù)法理論，薄板的彎曲振動響應(yīng)表示為一系列結(jié)構(gòu)波函數(shù)的疊加

(7)

結(jié)構(gòu)波函數(shù)為

Ψs=e-j(kxsxs+kysys)

(8)

式中:Lxs,Lys分別為薄板外輪廓尺寸;j2=-1;s1=0,1,2,…,ns1;s2=0,1,2,…,ns2,ns=4(ns1+1)+4(ns2+1)波函數(shù)數(shù)量也即為模型的自由度。而對于外部點激勵定義為

(9)

對于各波函數(shù)權(quán)系數(shù)的計算，采取類似有限元中加遼金法等邊界加權(quán)余量法，將邊界誤差最小化為0，經(jīng)整理得到計算各權(quán)系數(shù)Ws的系統(tǒng)矩陣

[Ass]{Ws}={fs}

(10)

其中,

(11)

(12)

3.2結(jié)構(gòu)聲耦合系統(tǒng)

同樣，對于聲學(xué)響應(yīng)，也可表示為一系列的聲學(xué)波函數(shù)的疊加

φi≈[Φ]{ga}

(13)

式中:ga為未知的權(quán)系數(shù);Φ為滿足Helmholtz控制方程的波函數(shù)，定義為

Φ=e-j(kxax+kyay+kzaz)

(14)

由于結(jié)構(gòu)與聲腔的耦合作用，在接合面引入聲學(xué)波函數(shù)附加作用的函數(shù)

(15)

式中,Ωs為結(jié)構(gòu)與聲域耦合邊界。因此結(jié)構(gòu)位移表示式相應(yīng)變換為

(16)

(17)

其中,

(18)

(19)

(20)

4 數(shù)值驗證

對于上述基于波函數(shù)法對含有自由層結(jié)構(gòu)的聲振分析方法，以下通過兩個例子進驗證。一個為敷設(shè)自由阻尼薄板的彎曲振動響應(yīng)分析，另一個為含有自由阻尼薄板與聲腔耦合系統(tǒng)的聲學(xué)響應(yīng)。

4.1含自由阻尼矩形薄板

如圖2所示尺寸的四邊固支的薄板結(jié)構(gòu)。引入的薄板為0.5 mm的鋁板，在其表面添加1 mm厚的阻尼層，兩種材料的參數(shù)如表1所示。在點(0.2 m，0.3 m)位置處施加一個法向的單位力，R點為選取的頻響觀測點，其平面坐標為(0.44 m，0.11 m)

圖2 自由阻尼薄板Fig.2 Geometry of unconstrained plate

表1 復(fù)合板材料參數(shù)Tab.1 Material parameters of plate and damp

波函數(shù)法通過Matlab2016a平臺編程實現(xiàn)，而有限元模型則在MSC/Nastran建立。為揭示有限元的特征以及其在更高頻率的限制，建立了精細和普通單元尺寸的有限元模型加以對比分析。其中精細有限元模型包含了60 000個四邊形單元與60 000個六面體單元，而普通單元尺寸有限元模型包含2 400個四邊形單元與2 400個六面體單元。取240 Hz時薄板的法向振動響應(yīng)結(jié)果如圖3所示。波函數(shù)的計算結(jié)果與有限元符合較好。同時為了在頻域上對比兩種方法，選取了響應(yīng)點在50～500 Hz段的響應(yīng)。圖4為響應(yīng)的實部，圖5為響應(yīng)的虛部。結(jié)果同樣表明，在整個分析頻段內(nèi)，波函數(shù)法的計算結(jié)果與精細建模的結(jié)果吻合較好，而單元尺寸較大的模型則在較多頻段出現(xiàn)較大差別，特別是在頻率較高的，其偏移更加明顯。而未敷設(shè)阻尼的薄板振動響應(yīng)結(jié)果如圖6所示，薄板受其本身特征頻率影響，在很多頻率下的出現(xiàn)很大響應(yīng)，也說明了阻尼的在抑制其振動的效果十分明顯。有無阻尼材料模型的計算結(jié)果都驗證了波函數(shù)法對預(yù)測自由阻尼薄板的振動響應(yīng)的有效性。

(a) WBM

(b) FEM圖3 薄板在240 Hz法向振動位移響應(yīng)的等高線圖Fig.3 Displacement response contour of damped plate at 240 Hz

圖4 含阻尼薄板的振動位移頻響曲線(實部)Fig.4 Response plot of reference point with damp (real part)

圖5 含阻尼薄板的振動位移頻響曲線(虛部)Fig.5 Response plot of reference point with damp (imaginary part)

圖6 無阻尼薄板的振動位移響應(yīng)(實部)Fig.6 Response plot of reference point without damp (real part)

4.2含自由阻尼薄板與聲腔耦合系統(tǒng)

薄板結(jié)構(gòu)與聲腔耦合系統(tǒng)是工業(yè)上廣泛應(yīng)用的系統(tǒng)結(jié)構(gòu)，如汽車車身與內(nèi)部聲腔。而添加阻尼也成為噪聲控制的重要和基本的手段。在驗證了以波函數(shù)法實現(xiàn)含有自由阻尼的薄板彎曲振動的基礎(chǔ)上，進一步對此類結(jié)構(gòu)的聲固耦合系統(tǒng)的響應(yīng)進行分析驗證。該耦合系統(tǒng)的幾何及尺寸如圖7所示，在板上取一點(0.2 m，0.2 m，0.5 m)施加單元力，取R點(0.4 m，0.3 m，0.2 m)處的聲壓響應(yīng)。

該基于波函數(shù)法的耦合模型也在Matlab中建立和完成后處理。有限元耦合模型則在LMS.Virtual.lab中建立和計算，其中精細的有限元模型除了結(jié)構(gòu)與阻尼的單元外增加了3 000 000個聲學(xué)六面體單元，而普通單元尺寸模型則增加了120 000個聲學(xué)六面體單元。為了驗證阻尼的作用，同樣也對無阻尼耦合模型進行了計算。得到的頻響曲線如圖8、圖9所示。與薄板振動響應(yīng)的結(jié)果一致，基于波函數(shù)法的計算結(jié)果與精細的有限元模型更接近，而普通單元尺寸的有限元模型在較高頻率段的偏移明顯。添加阻尼后的模型在很多波峰處的響應(yīng)得到了很大抑制，整體的聲壓級水平也較無阻尼系統(tǒng)要小。

圖7 薄板聲腔耦合系統(tǒng)Fig.7 Geometry of coupled vibro-acoutic system

圖8 含阻尼耦合系統(tǒng)聲學(xué)響應(yīng)Fig.8 Acoustic response of reference point with damp

圖9 不含阻尼耦合系統(tǒng)聲學(xué)響應(yīng)Fig.9 Acoustic response of reference point without damp

4.3收斂性

通過以上的算例可以看出波函數(shù)法在計算含有自由阻尼的薄板聲腔等系統(tǒng)的聲振響應(yīng)的正確性，然而該方法能更好的處理中頻問題則主要是其擁有高計算精度與效率。為描述計算方法的收斂性，定義相對誤差ε

(23)

式中:R(p)為選取點的振動或聲學(xué)響應(yīng);Rref(p)為參考響應(yīng)值，此處取精細建模的有限元計算結(jié)果。

圖10、圖11分別為自由阻尼薄板算例與結(jié)構(gòu)聲耦合算例在相同硬件條件下，在100 Hz時關(guān)于計算時間的收斂曲線。雖然波函數(shù)法的計算時間包括建模與計算的時間，而有限元法的物理計算時間只包括模型計算時間(前期有限元模型劃分與生成的時間未包含)，但對于該類較簡單的結(jié)構(gòu)波函數(shù)法模型的構(gòu)建時間很短，對對比結(jié)果影響較小。從收斂曲線可以看出，在相同計算精度的情況，波函數(shù)法的計算效率高于有限元法。而對于更高頻段的聲振響應(yīng)問題，就要求單元尺寸要足夠小來捕捉短波特性。這就使得模型計算量變得十分巨大，而伴隨的還有邊界離散帶來的誤差累積，所以基于單元的方法主要用于求解低頻問題。波函數(shù)法沒有如單元類方法采用多項式擬合造成的差值誤差，同時相比之下其計算量也非常小，因此其能更好的應(yīng)用于中頻問題以及需要重復(fù)計算類問題如優(yōu)化迭代。

圖10 薄板結(jié)構(gòu)中WBM與FEM的收斂曲線Fig.10 The convergence of WBM and FEM of the plate model

圖11 耦合系統(tǒng)中WBM與FEM的收斂曲線Fig.11 The convergence of WBM and FEM of the coupled system

5 結(jié) 論

文章在完成敷設(shè)自由阻尼薄板的彎曲振動響應(yīng)的分析基礎(chǔ)上，應(yīng)用波函數(shù)法理論，推導(dǎo)了自由阻尼薄板結(jié)構(gòu)以及其與聲腔耦合的WBM模型，實現(xiàn)了波函數(shù)法對含有自由阻尼結(jié)構(gòu)的聲振預(yù)測。以兩個數(shù)值模型驗證該方法的正確性，得到了自由阻尼在中頻對薄板類結(jié)構(gòu)的減振作用與效果。通過與有限元法在頻響結(jié)果與收斂性的對比，體現(xiàn)了波函數(shù)法更能有效的處理含有自由阻尼結(jié)構(gòu)系統(tǒng)在更高頻段的聲振響應(yīng)問題，也說明了此次探究的有效性。后續(xù)將在此基礎(chǔ)上，探尋基于Middlin理論下，含有阻尼的板件類弱結(jié)構(gòu)在中頻的響應(yīng)預(yù)測方法，以及基于結(jié)構(gòu)或聲學(xué)響應(yīng)的阻尼優(yōu)化。

[1] DERAEMAEKER A, BABUSKA I, BOUILLARD P. Dispersion and pollution of the FEM solution for the helmholtz equation in one, two and three dimensions [J]. International Journal for Numerical Methods in Engineering, 1999,46(4):471-499.

[2] ZIENKIEWICZ O C, TAYLOR R L, ZHU J Z, et al. The finite element method-Vol. 1: Basic formulation and linear problems[M]. Oxford: Butterworth-Heinemann, 2005.

[3] LYON R, DEJONG R. Theory and application of statistical energy analysis [M]. 2nd ed. Oxford: Butterworth Heinemann, 1995.

[4] MENCIK J M, ICHCHOU M N. Wave finite elements in guided elastodynamics with internal fluid[J]. International Journal of Solids and Structures, 2007,44:2148-2167.

[5] MAXIT L, GUYADER J L. Extension of SEA model to subsystems with non-uniform modal energy distribution[J]. Journal of Sound and Vibration, 2003,265(2):337-358.

[6] SESTIERI A, CARCATERRA A. On the spurious solutions in complex envelope displacement analysis[J]. Journal of Sound and Vibration, 2001,240(2):293-302.

[7] LADEVEZE P, ARNAUD L, ROUCH P, et al. The variational theory of complex rays for the calculation of medium- frequency vibrations[J]. Engineering Computations, 2001, 18(1/2):193-214.

[8] LANGLEY R S, CORDIOLI J A. Hybrid deterministic-statistical analysis of vibroacoustic system domain couplings on statistical components[J]. Journal of Sound and Vibration, 2009, 321(3/4/5):893-912.

[9] MACE B R. Statistical energy analysis, energy distribution models and system modes[J]. Journal of Sound and Vibration, 2003, 264(2):391-409.

[10] DESMET W. A wave based prediction technique for coupled vibro-acoustic analysis[D]. Leuven：Katholieke Universiteit Leuven，1998.

[11] VANMAELE C, VANDEPITTE D. An efficient wave based prediction technique for plate bending vibrations[J]. Computer Methods in Applied Mechanics and Engineering, 2007,196(33/34): 3178-3189.

[12] VAN GENECHTEN B, ATAK O, BERGEN B, et al. An efficient wave based method for solving helmholtz problems in three-dimensional bounded domains[J]. Engineering Analysis with Boundary Elements, 2012,36(1):63-75.

[13] DESMET W, VAN HAL B, SAS P, et al. A computationally efficient prediction technique for the steady-state dynamic analysis of coupled vibro-acoustic systems[J]. Advances in Engineering Software, 2002,33(7/8/9/10):527-540.

[14] DECKERS E, HORLIN D E, VANDEPITTE D, et al. A wave based method for the efficient solution of the 2d poroelastic biot equations[J]. Computer Methods in Applied Mechanics and Engineering, 2012, 201/202/203/204:245-262.

[15] 何雪松,黃其柏,胡溧.WBM 法在薄板彎曲振動分析中的應(yīng)用[J].華中科技大學(xué)學(xué)報(自然科學(xué)版),2008,36(7):97-99.

HE Xuesong, HUANG Qibai, HU Li. Application of wave based method to plate bending vibration analysis [J]. Journal of Huazhong University of Science and Technology (Natural Science), 2008, 36(7):97-99.

[16] PENG Weicai, HE Zeng, LI Peng, et al. A prediction technique for dynamic analysis of flat plates in the Mid-Frequency Range[J]. Acta Mechanica Solida Sinica, 2007, 20(4): 333-341.

[17] 何锃, 彭偉才,王加強. WB法分析結(jié)構(gòu)與聲耦合問題[J]. 華中科技大學(xué)學(xué)報(自然科學(xué)版), 2007,35(8): 119-121.

HE Zeng,PENG Weicai, WANG Jiaqiang. Wave based method for coupled structural—acoustic analysis[J]. Journal of Huazhong University of Science & Technology (Natural Science),2007,35(8): 119-121.

[18] 桂洪斌，趙德有，金咸定. 自由阻尼層加筋板的穩(wěn)態(tài)簡諧響應(yīng)分析[J]. 上海交通大學(xué)學(xué)報,2002,36(11):1544-1547.

GUI Hongbin, ZHAO Deyou, JIN Xianding. Steady-state harmonic response analys is of stiffened plate with unconstrained damped layer[J]. Journal of Shanghai Jiaotong University, 2002,36(11):1544-1547.

[19] REDDY C V R，GANESAN N，RAO B V A，et al. Response of plates with unconstrained layer damping treatment to random acoustic excitation[J]. Journal of Sound and Vibration,1980,69(1): 35-43.

Responsesofunconstraineddampedvibro-acousticsystemsusingthewavebasedpredictiontechnique

XIA Xiaojun1,2, XU Zhongming1,2, LAI Shiyang3, ZHANG Zhifei1,2, HE Yansong1,2

(1. State Key Laboratory of Mechanical Transmission, Chongqing University, Chongqing 400030, China;2. School of Automotive Engineering, Chongqing University, Chongqing 400030, China;3. School of Mechanical Engineering, Chongqing Vocational Institute of Engineering, Chongqing 402260, China)

The global governing vibration equation of an unconstrained damped plate was deduced based on the Kirchhoff theory and the analysis of the complex stiffness. The methodology for predicting the vibration of the unconstrained damped plate and the acoustics of the coupled 3D vibro-acoustic system was proposed by virtue of the wave based method (WBM). Taking a four edges clamped rectangular plate coupled with a box like vibro-acoustic system as a numerical example, the out-plane displacement of the unconstrained damped plate was analysed. The responses at a selected reference point was calculated in the 50-500 Hz frequency band by using the WBM and FEM respectively. The results of the two method validate that the WBM is capable for predicting the vibration and acoustic responses of the unconstrained damped system effectively, and the WBM is more efficient to deal with vibro-acoustic problems compared with the FEM.

unconstrained damping; wave based method; bending vibration; structural-acoustic coupling

U467.4

A

10.13465/j.cnki.jvs.2017.19.024

重慶市基礎(chǔ)與前沿研究計劃項目(CSTC2015JCYJBX0075);中央高?；究蒲袠I(yè)務(wù)費(106112016CDJZR335522)

2016-05-04 修改稿收到日期：2016-08-10

夏小均 男，博士生，1988年生

徐中明 男，博士，教授，博士生導(dǎo)師，1963年生