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

(1.重慶大學(xué) 機(jī)械傳動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,重慶 400030；2.重慶大學(xué) 汽車工程學(xué)院,重慶 400030;3.重慶工程職業(yè)技術(shù)學(xué)院 機(jī)械工程學(xué)院,重慶 402260)

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

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

(1.重慶大學(xué) 機(jī)械傳動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,重慶 400030；2.重慶大學(xué) 汽車工程學(xué)院,重慶 400030;3.重慶工程職業(yè)技術(shù)學(xué)院 機(jī)械工程學(xué)院,重慶 402260)

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

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

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

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

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

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

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

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

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

(1)

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

(2)

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

(3)

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

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

(4)

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

(5)

(6)

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

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

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

根據(jù)波函數(shù)法理論，薄板的彎曲振動(dòng)響應(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ù)量也即為模型的自由度。而對(duì)于外部點(diǎn)激勵(lì)定義為

(9)

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

[Ass]{Ws}={fs}

(10)

其中,

(11)

(12)

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

同樣，對(duì)于聲學(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ù)值驗(yàn)證

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

4.1含自由阻尼矩形薄板

如圖2所示尺寸的四邊固支的薄板結(jié)構(gòu)。引入的薄板為0.5 mm的鋁板，在其表面添加1 mm厚的阻尼層，兩種材料的參數(shù)如表1所示。在點(diǎn)(0.2 m，0.3 m)位置處施加一個(gè)法向的單位力，R點(diǎn)為選取的頻響觀測(cè)點(diǎn)，其平面坐標(biāo)為(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ù)法通過(guò)Matlab2016a平臺(tái)編程實(shí)現(xiàn)，而有限元模型則在MSC/Nastran建立。為揭示有限元的特征以及其在更高頻率的限制，建立了精細(xì)和普通單元尺寸的有限元模型加以對(duì)比分析。其中精細(xì)有限元模型包含了60 000個(gè)四邊形單元與60 000個(gè)六面體單元，而普通單元尺寸有限元模型包含2 400個(gè)四邊形單元與2 400個(gè)六面體單元。取240 Hz時(shí)薄板的法向振動(dòng)響應(yīng)結(jié)果如圖3所示。波函數(shù)的計(jì)算結(jié)果與有限元符合較好。同時(shí)為了在頻域上對(duì)比兩種方法，選取了響應(yīng)點(diǎn)在50～500 Hz段的響應(yīng)。圖4為響應(yīng)的實(shí)部，圖5為響應(yīng)的虛部。結(jié)果同樣表明，在整個(gè)分析頻段內(nèi)，波函數(shù)法的計(jì)算結(jié)果與精細(xì)建模的結(jié)果吻合較好，而單元尺寸較大的模型則在較多頻段出現(xiàn)較大差別，特別是在頻率較高的，其偏移更加明顯。而未敷設(shè)阻尼的薄板振動(dòng)響應(yīng)結(jié)果如圖6所示，薄板受其本身特征頻率影響，在很多頻率下的出現(xiàn)很大響應(yīng)，也說(shuō)明了阻尼的在抑制其振動(dòng)的效果十分明顯。有無(wú)阻尼材料模型的計(jì)算結(jié)果都驗(yàn)證了波函數(shù)法對(duì)預(yù)測(cè)自由阻尼薄板的振動(dòng)響應(yīng)的有效性。

(a) WBM

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

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

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

圖6 無(wú)阻尼薄板的振動(dòng)位移響應(yīng)(實(shí)部)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)部聲腔。而添加阻尼也成為噪聲控制的重要和基本的手段。在驗(yàn)證了以波函數(shù)法實(shí)現(xiàn)含有自由阻尼的薄板彎曲振動(dòng)的基礎(chǔ)上，進(jìn)一步對(duì)此類結(jié)構(gòu)的聲固耦合系統(tǒng)的響應(yīng)進(jìn)行分析驗(yàn)證。該耦合系統(tǒng)的幾何及尺寸如圖7所示，在板上取一點(diǎn)(0.2 m，0.2 m，0.5 m)施加單元力，取R點(diǎn)(0.4 m，0.3 m，0.2 m)處的聲壓響應(yīng)。

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

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

(23)

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

圖10、圖11分別為自由阻尼薄板算例與結(jié)構(gòu)聲耦合算例在相同硬件條件下，在100 Hz時(shí)關(guān)于計(jì)算時(shí)間的收斂曲線。雖然波函數(shù)法的計(jì)算時(shí)間包括建模與計(jì)算的時(shí)間，而有限元法的物理計(jì)算時(shí)間只包括模型計(jì)算時(shí)間(前期有限元模型劃分與生成的時(shí)間未包含)，但對(duì)于該類較簡(jiǎn)單的結(jié)構(gòu)波函數(shù)法模型的構(gòu)建時(shí)間很短，對(duì)對(duì)比結(jié)果影響較小。從收斂曲線可以看出，在相同計(jì)算精度的情況，波函數(shù)法的計(jì)算效率高于有限元法。而對(duì)于更高頻段的聲振響應(yīng)問(wèn)題，就要求單元尺寸要足夠小來(lái)捕捉短波特性。這就使得模型計(jì)算量變得十分巨大，而伴隨的還有邊界離散帶來(lái)的誤差累積，所以基于單元的方法主要用于求解低頻問(wèn)題。波函數(shù)法沒(méi)有如單元類方法采用多項(xiàng)式擬合造成的差值誤差，同時(shí)相比之下其計(jì)算量也非常小，因此其能更好的應(yīng)用于中頻問(wèn)題以及需要重復(fù)計(jì)算類問(wèn)題如優(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è)自由阻尼薄板的彎曲振動(dòng)響應(yīng)的分析基礎(chǔ)上，應(yīng)用波函數(shù)法理論，推導(dǎo)了自由阻尼薄板結(jié)構(gòu)以及其與聲腔耦合的WBM模型，實(shí)現(xiàn)了波函數(shù)法對(duì)含有自由阻尼結(jié)構(gòu)的聲振預(yù)測(cè)。以兩個(gè)數(shù)值模型驗(yàn)證該方法的正確性，得到了自由阻尼在中頻對(duì)薄板類結(jié)構(gòu)的減振作用與效果。通過(guò)與有限元法在頻響結(jié)果與收斂性的對(duì)比，體現(xiàn)了波函數(shù)法更能有效的處理含有自由阻尼結(jié)構(gòu)系統(tǒng)在更高頻段的聲振響應(yīng)問(wèn)題，也說(shuō)明了此次探究的有效性。后續(xù)將在此基礎(chǔ)上，探尋基于Middlin理論下，含有阻尼的板件類弱結(jié)構(gòu)在中頻的響應(yīng)預(yù)測(cè)方法，以及基于結(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 法在薄板彎曲振動(dòng)分析中的應(yīng)用[J].華中科技大學(xué)學(xué)報(bào)(自然科學(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] 何锃, 彭偉才,王加強(qiáng). WB法分析結(jié)構(gòu)與聲耦合問(wèn)題[J]. 華中科技大學(xué)學(xué)報(bào)(自然科學(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)簡(jiǎn)諧響應(yīng)分析[J]. 上海交通大學(xué)學(xué)報(bào),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ǔ)與前沿研究計(jì)劃項(xiàng)目(CSTC2015JCYJBX0075);中央高?；究蒲袠I(yè)務(wù)費(fèi)(106112016CDJZR335522)

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

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

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