約束阻尼板結(jié)構(gòu)振動(dòng)聲輻射優(yōu)化

鄭 玲, 祝喬飛

(重慶大學(xué) 機(jī)械傳動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，重慶 400044)

盡管主動(dòng)振動(dòng)噪聲控制能有效抑制結(jié)構(gòu)的低頻噪聲，但低噪聲設(shè)計(jì)以其優(yōu)異的性能和可靠性而備受關(guān)注。Kikuchi[1]首先將拓?fù)鋬?yōu)化方法應(yīng)用于結(jié)構(gòu)低噪聲設(shè)計(jì)，通過(guò)材料布局優(yōu)化，實(shí)現(xiàn)了結(jié)構(gòu)低噪聲設(shè)計(jì)。Luo等[2]采用拓?fù)鋬?yōu)化方法，通過(guò)對(duì)封閉結(jié)構(gòu)的優(yōu)化，大大降低了封閉結(jié)構(gòu)內(nèi)部的聲輻射。Baz等[3]采用有限元法，建立了板-聲腔耦合動(dòng)力學(xué)模型，以聲腔模態(tài)頻率峰值聲壓最小化為目標(biāo)函數(shù)，對(duì)板的材料布局進(jìn)行了拓?fù)鋬?yōu)化設(shè)計(jì)，實(shí)驗(yàn)驗(yàn)證了理論分析結(jié)果的正確性。對(duì)結(jié)構(gòu)表面進(jìn)行阻尼處理，能有效降低結(jié)構(gòu)的振動(dòng)與聲輻射，廣泛應(yīng)用于航空、航天、汽車(chē)及船舶等領(lǐng)域。阻尼結(jié)構(gòu)可分為自由阻尼和約束層阻尼，相比自由阻尼，約束層阻尼消耗的能量更多，減振作用更大。

如何根據(jù)結(jié)構(gòu)的動(dòng)力學(xué)特性要求，確定阻尼材料的最優(yōu)布局，是阻尼結(jié)構(gòu)設(shè)計(jì)的一個(gè)關(guān)鍵問(wèn)題。郭中澤等[4]采用漸進(jìn)結(jié)構(gòu)拓?fù)鋬?yōu)化方法，以結(jié)構(gòu)模態(tài)損耗因子最大化為目標(biāo)，研究了阻尼結(jié)構(gòu)設(shè)計(jì)中阻尼材料的最優(yōu)分布問(wèn)題。楊德慶等[5]采用均勻化拓?fù)鋬?yōu)化方法，以薄板模態(tài)頻率最大化為目標(biāo)，研究了薄板的最優(yōu)材料布局問(wèn)題。李以農(nóng)等[6]采用漸進(jìn)結(jié)構(gòu)拓?fù)鋬?yōu)化算法，以結(jié)構(gòu)模態(tài)阻尼比最大化為目標(biāo)函數(shù)，研究了約束層阻尼梁的阻尼材料最優(yōu)分布問(wèn)題。鄭玲等[7]建立了約束阻尼板的有限元模型，推導(dǎo)了模態(tài)阻尼比對(duì)單元胞位置的靈敏度，提出了一種基于進(jìn)化方式的約束阻尼優(yōu)化布局方法。李超等[8]研究了圓柱殼體上阻尼材料的最優(yōu)分布問(wèn)題。王明旭等[9]將變密度法引入阻尼材料的布局設(shè)計(jì)，推導(dǎo)了靈敏度表達(dá)式，提出了模態(tài)振型的跟蹤方法。呂毅寧等[10]給出了附加自由阻尼復(fù)雜結(jié)構(gòu)的有限元建模方法，以結(jié)構(gòu)的動(dòng)態(tài)柔度最小化為目標(biāo)函數(shù)，計(jì)算了薄板和車(chē)身地板附加自由阻尼的拓?fù)鋬?yōu)化布局。以上文獻(xiàn)均以結(jié)構(gòu)振動(dòng)損耗因子、固有頻率或動(dòng)柔度為目標(biāo)函數(shù)，采用拓?fù)鋬?yōu)化方法，研究了阻尼材料的最優(yōu)布局。

在某些特殊環(huán)境下，結(jié)構(gòu)振動(dòng)導(dǎo)致的聲輻射更為人們所關(guān)注，如潛艇的結(jié)構(gòu)靜音設(shè)計(jì)、飛機(jī)或汽車(chē)艙內(nèi)噪聲的控制等。目前，以聲壓或聲功率等聲學(xué)性能為目標(biāo)函數(shù)，采用拓?fù)鋬?yōu)化方法，對(duì)結(jié)構(gòu)進(jìn)行優(yōu)化已有大量研究。Sigmund等[11]將拓?fù)鋬?yōu)化方法應(yīng)用于室內(nèi)聲學(xué)設(shè)計(jì)和噪聲屏蔽墻設(shè)計(jì)，對(duì)比了優(yōu)化前后的材料布局與聲壓均方值變化。Sun等[12]建立了蜂窩夾層梁的有限元模型，以特定頻率輻射聲功率最小化為目標(biāo)函數(shù)，對(duì)蜂窩的形狀及尺寸進(jìn)行了優(yōu)化。但是，采用拓?fù)鋬?yōu)化方法，以輻射聲功率最小化為目標(biāo)，對(duì)約束阻尼材料布局進(jìn)行優(yōu)化的研究尚不多見(jiàn)。

本文首先建立四邊簡(jiǎn)支約束層阻尼板有限元模型，用Rayleigh積分推導(dǎo)結(jié)構(gòu)輻射聲功率，并進(jìn)行靈敏度分析，以聲功率最小化為目標(biāo)函數(shù)，采用漸進(jìn)優(yōu)化算法，對(duì)約束阻尼材料布局進(jìn)行優(yōu)化。

1 約束層阻尼板振動(dòng)和聲輻射計(jì)算

1.1 約束層阻尼板有限元模型

圖1 約束層阻尼板單元

約束阻尼單元模型如圖1所示。每個(gè)單元均由基板層、粘彈性(VEM)層和約束層組成，下標(biāo)分別用p、v、c表示。每個(gè)節(jié)點(diǎn)有7個(gè)自由度，分別為約束層中性面內(nèi)x方向的位移uci、約束層中性面內(nèi)y方向的位移vci、基板層中性面內(nèi)x方向的位移upi、基板層中性面內(nèi)y方向的位移vpi、單元中性面的橫向位移wi以及繞x軸方向的轉(zhuǎn)角θxi和y軸方向的轉(zhuǎn)角θyi。

節(jié)點(diǎn)位移向量為：

Δi={ucivciupivpiwiθxiθyi}T

(i=1,2,3,4)

(1)

單元位移向量為：

(2)

根據(jù)單元位移模式，可得約束層阻尼單元的形函數(shù)[13]：

N={[Nuc];[Nvc];[Nup];

[Nvp];[Nw];[Nwx];[Nwy]}

(3)

其中，Nuc、Nvc、Nup、Nvp、Nw、Nwx、Nwy分別為對(duì)應(yīng)單元節(jié)點(diǎn)自由度uc、vc、up、vp、w、θx、θy的形函數(shù)。

約束阻尼單元各層的動(dòng)能和應(yīng)變勢(shì)能如下[14]：

(4)

(5)

(Nvc)TNvc+(Nw)TNw]dxdyΔ(e)=

(6)

(7)

(Nvv)TNvv+(Nw)TNw]dxdyΔ(e)=

(8)

(9)

(10)

其中，T和P分別表示結(jié)構(gòu)單元的動(dòng)能和應(yīng)變勢(shì)能。運(yùn)用Hamilton原理，可導(dǎo)出約束層阻尼單元的質(zhì)量矩陣和剛度矩陣：

m=mp+mc+mv

(11)

k=kp+kc+kv+kβv

(12)

其中，mp、mc、mv分別為基板、約束層以及粘彈性層的質(zhì)量矩陣，kp、kc、kv、kβv分別為基板、約束層、粘彈性層位移及橫向剪切變形的剛度矩陣。將單元質(zhì)量矩陣和剛度矩陣進(jìn)行組集，獲得約束層阻尼板結(jié)構(gòu)的總體質(zhì)量矩陣和剛度矩陣：

(13)

(14)

因此，約束層阻尼板的有限元?jiǎng)恿W(xué)模型為：

(15)

這里，粘彈性層的剪切模量為G=Gv(1+ηi)，η為損耗因子，[K]為復(fù)剛度矩陣。

1.2 約束層阻尼板的振動(dòng)模型

假設(shè)約束層阻尼板的外激勵(lì)為簡(jiǎn)諧激勵(lì)f(t)=Fejωt，位移響應(yīng)滿(mǎn)足x(t)=Xejωt。將以上兩式代入(15)式，消去ejωt，得：

(-ω2[M]+[K]){X}={F}

(16)

因此，約束層阻尼板節(jié)點(diǎn)位移響應(yīng)為：

{X}=[H(ω)]{F}

(17)

節(jié)點(diǎn)速度響應(yīng)為：

{V}=iω{X}=iω[H(ω)]{F}

(18)

其中：H(ω)=(-ω2[M]+[K])-1為位移傳遞函數(shù)。

1.3 約束層阻尼板聲輻射計(jì)算

將約束層阻尼板鑲嵌于無(wú)限大剛性障板上，運(yùn)用簡(jiǎn)化后的Helmholtz公式，即Rayleigh積分計(jì)算其聲輻射。振動(dòng)輻射聲功率可由表面聲壓p(P)和表面法向速度vn(P)表示：

(19)

(20)

將(20)式代入(19)式：

(21)

考慮到格林函數(shù)的對(duì)稱(chēng)性，式(21)可表示為：

(22)

將式(22)寫(xiě)成矩陣形式，其中{V}為各單元上法向速度構(gòu)成的向量。

W(ω)={V}T[R]{V}

(23)

在結(jié)構(gòu)均勻離散情況下，矩陣R的形式為[15]：

(24)

在一個(gè)頻帶內(nèi)的平均聲功率表示為：

(25)

離散后得：

(26)

其中，Wave表示頻帶內(nèi)的平均輻射聲功率，N表示頻帶內(nèi)激勵(lì)力頻率的總數(shù)，ωi表示第i個(gè)激勵(lì)力頻率。

2 優(yōu)化模型

以簡(jiǎn)諧激勵(lì)作用下約束阻尼板輻射的聲功率最小化為目標(biāo)函數(shù)，約束阻尼材料體積分?jǐn)?shù)為約束條件，建立拓?fù)鋬?yōu)化模型：

(27)

其中，X={β1β2…βn}T為設(shè)計(jì)變量，表示約束阻尼單元的存在狀態(tài)。當(dāng)βi=1時(shí)，表明第i個(gè)單元位置有約束阻尼單元，當(dāng)βi=0時(shí)，表明第i個(gè)單元位置無(wú)約束阻尼單元，n為約束阻尼單元總數(shù)，W為約束阻尼板輻射聲功率，V為約束阻尼材料體積，V*是體積約束。

3 靈敏度分析

約束阻尼板的速度響應(yīng)為：

V=iωX

(28)

速度對(duì)設(shè)計(jì)變量β的偏導(dǎo)數(shù)為：

(29)

將式(15)兩端對(duì)設(shè)計(jì)變量h求偏導(dǎo)數(shù)：

(30)

將上式整理并代入式(29)，可得速度對(duì)設(shè)計(jì)變量的靈敏度[16]：

(31)

結(jié)合式(23)，可求出約束阻尼板輻射聲功率對(duì)設(shè)計(jì)變量的靈敏度：

(32)

這里的設(shè)計(jì)變量β是約束阻尼單元的存在狀態(tài)，根據(jù)式(31)，結(jié)合設(shè)計(jì)變量的定義，速度對(duì)設(shè)計(jì)變量的靈敏度可用約束阻尼單元?jiǎng)h除前后的速度變化量來(lái)表示。同樣，根據(jù)式(32)，聲功率對(duì)設(shè)計(jì)變量的靈敏度也可以用約束單元?jiǎng)h除前后的聲功率變化值來(lái)表示。因此，速度對(duì)設(shè)計(jì)變量的靈敏度改寫(xiě)為：

ΔV=-iω(K-ω2M)-1(ΔK-ω2ΔM)X

(33)

聲功率對(duì)設(shè)計(jì)變量的靈敏度改寫(xiě)為：

ΔW=ΔVHRV+VHRΔV

(34)

4 獨(dú)立網(wǎng)格濾波

在連續(xù)結(jié)構(gòu)拓?fù)鋬?yōu)化過(guò)程中，常出現(xiàn)棋盤(pán)格式和網(wǎng)格依賴(lài)性等數(shù)值不穩(wěn)定現(xiàn)象。采用獨(dú)立網(wǎng)格濾波技術(shù)，以有效抑制棋盤(pán)格式和網(wǎng)格依賴(lài)性。獨(dú)立網(wǎng)格濾波技術(shù)屬于靈敏度過(guò)濾法，是一種局部意義上的約束方法，它以過(guò)濾半徑范圍內(nèi)單元的靈敏度的加權(quán)平均值代替原來(lái)的靈敏度值。獨(dú)立網(wǎng)格濾波方程為:

(35)

采用獨(dú)立網(wǎng)格濾波技術(shù)，最重要的是選擇合適的濾波半徑，對(duì)阻尼單元的靈敏度信息進(jìn)行濾波操作，它不僅能抑制棋盤(pán)格式，而且能消除網(wǎng)格依賴(lài)性。

5 優(yōu)化流程

本文采用漸進(jìn)優(yōu)化算法，求解優(yōu)化模型。漸進(jìn)優(yōu)化算法的基本思想是通過(guò)逐步刪除設(shè)計(jì)域中的無(wú)效或低效單元，搜尋最佳的材料布局，使目標(biāo)函數(shù)最小化或最大化。基于約束阻尼板輻射聲功率的拓?fù)鋬?yōu)化模型，可設(shè)計(jì)使約束層阻尼板輻射聲功率最小化的拓?fù)鋬?yōu)化計(jì)算流程，見(jiàn)圖2所示。

圖2 約束阻尼板聲功率優(yōu)化設(shè)計(jì)流程圖

優(yōu)化步驟如下：

① 建立約束阻尼板的有限元模型和邊界元模型；

② 給定約束阻尼材料的體積約束量以及初始刪除量RR0和其進(jìn)化量ER；

③ 對(duì)激勵(lì)力作用下的結(jié)構(gòu)進(jìn)行聲輻射分析，獲得輻射聲功率；

④ 按式(34)計(jì)算各阻尼單元的聲功率靈敏度；

⑤ 設(shè)定適當(dāng)?shù)臑V波半徑，采用獨(dú)立網(wǎng)格濾波技術(shù)對(duì)各單元的聲功率靈敏度進(jìn)行濾波；

⑥ 對(duì)濾波后的單元聲功率靈敏度進(jìn)行排序；

⑦ 刪除RR0個(gè)單元的聲功率靈敏度值最小的單元；

⑧ 判定是否達(dá)到約束條件。如果不滿(mǎn)足約束條件，更新RR0=RR0+ER，并重復(fù)步驟③-⑧；如果滿(mǎn)足約束條件，則結(jié)束迭代，輸出結(jié)果。

6 算例分析

將四邊簡(jiǎn)支的約束阻尼板鑲嵌于無(wú)限大剛性障板上，計(jì)算約束層阻尼板的振動(dòng)聲輻射。約束層阻尼板的參數(shù)如下：

約束阻尼板長(zhǎng)L=0.4 m，寬B=0.3 m。基層：Ep=7.0e10 Pa，ρp=2 800 kg/m3，μp=0.3，hp=0.000 8 m；粘彈性層：Ev=1.2e7 Pa，ρv=1 200 kg/m3，μv=0.495，hv=0.000 05 m，η=0.5；約束層：Ec=7.0e10 Pa，ρc=2 700 kg/m3，μc=0.3，hc=0.000 13 m。有限元網(wǎng)格單元數(shù)16×12，邊界元網(wǎng)格單元數(shù)16×12，在(0.20,0.15,0)處即中心點(diǎn)受到正弦激勵(lì)力：F=sinωt=sin(2πft)。空氣密度ρ0=1.225 kg/m3，空氣中聲速c0=343 m/s。

6.1 一階峰值聲功率的拓?fù)鋬?yōu)化

隨著阻尼材料的刪除，約束層阻尼板結(jié)構(gòu)的聲功率曲線(xiàn)會(huì)發(fā)生頻移，因此，針對(duì)約束層阻尼板一階峰值聲功率的拓?fù)鋬?yōu)化，必須在每次刪除單元后，計(jì)算輻射聲功率時(shí)，重新計(jì)算固有頻率，進(jìn)行頻率跟蹤。

優(yōu)化過(guò)程中，每次刪除4個(gè)單元，濾波半徑取為4.5，將192個(gè)單元全部刪除須進(jìn)行48次迭代。圖3是輻射聲功率隨刪除率變化的迭代過(guò)程圖。

圖3 聲功率隨阻尼材料刪除率的變化圖

由圖3可知，在阻尼材料刪除過(guò)程中，約束阻尼板的輻射聲功率有一定程度波動(dòng)，因?yàn)樵趧h除過(guò)程中被刪除的單元靈敏度可能為正也可能為負(fù)，從而導(dǎo)致刪除后約束阻尼結(jié)構(gòu)的輻射聲功率有波動(dòng)。其中，約束層阻尼板的聲功率有幾個(gè)低點(diǎn)，當(dāng)刪除率為12.5%和35.4%時(shí)，輻射聲功率均為112.1 dB，刪除率為62.5%時(shí)，輻射聲功率為112.2 dB，刪除率為81.2%時(shí)，輻射聲功率為113.3 dB。選取刪除率為35.4%、62.5%和81.2%，分析一階峰值聲功率，并計(jì)算最優(yōu)的拓?fù)錁?gòu)型圖。

約束阻尼最優(yōu)拓?fù)錁?gòu)型如圖4～圖6，黑色表示約束阻尼材料，白色為無(wú)約束阻尼材料，從動(dòng)力學(xué)特性的角度來(lái)分析，粘貼約束阻尼材料的區(qū)域均是一階模態(tài)應(yīng)變能較大的區(qū)域，表明這些單元均是有效單元。基板、全覆蓋板及不同阻尼材料用量下的輻射聲功率對(duì)比見(jiàn)圖7～圖8所示，可以看出，隨著約束阻尼材料用量的增加，聲功率曲線(xiàn)向后發(fā)生頻移，此外，在第一階峰值處，35.4%、62.5%和81.2%刪除率時(shí)，其輻射聲功率均比基板和全覆蓋板有所下降，由于拓?fù)鋬?yōu)化保留了有效阻尼單元，刪除了無(wú)效阻尼單元，在阻尼材料用量減小的情況下，獲得了比全覆蓋更低的輻射聲功率。

圖4 35.4%刪除率下阻尼材料的最優(yōu)拓?fù)鋱D

圖7 一階峰值優(yōu)化的輻射聲功率對(duì)比圖

圖10 頻帶優(yōu)化的輻射聲功率對(duì)比圖

圖11 頻段內(nèi)輻射聲功率對(duì)比圖

6.2 頻帶輻射聲功率的拓?fù)鋬?yōu)化

選取頻率區(qū)間10 Hz-60 Hz，這是基板第一階峰值所在的頻率范圍。阻尼材料刪除率取為50%，步長(zhǎng)為1 Hz，每次刪除4個(gè)單元，濾波半徑取為4.5。圖9是約束阻尼材料最優(yōu)拓?fù)錁?gòu)型，圖10是基板、全覆蓋以及最優(yōu)拓?fù)錁?gòu)型下的聲功率對(duì)比圖。

圖11表明，優(yōu)化后的約束阻尼板，其輻射平均聲功率比基板和全覆蓋板更低；在目標(biāo)頻率區(qū)間，優(yōu)化結(jié)構(gòu)的平均聲功率為99.7 dB，基板為111.4 dB，全覆蓋為102.1 dB，優(yōu)化后的約束阻尼板，其平均聲功率分別下降11.7 dB和2.4 dB。

7 結(jié) 論

本文建立了約束層阻尼板的有限元模型，用Rayleigh積分推導(dǎo)了約束層阻尼板的輻射聲功率及其靈敏度表達(dá)式，以輻射聲功率最小化為目標(biāo)，采用漸進(jìn)優(yōu)化算法對(duì)約束阻尼材料布局進(jìn)行了拓?fù)鋬?yōu)化。研究結(jié)果表明：本文提出的理論模型和拓?fù)鋬?yōu)化算法是正確和有效的，無(wú)論是針對(duì)一階峰值或是某一頻帶，它均能以較少的約束阻尼材料，實(shí)現(xiàn)結(jié)構(gòu)振動(dòng)聲輻射最小化的目標(biāo)。約束層阻尼板振動(dòng)聲輻射的優(yōu)化為結(jié)構(gòu)低噪聲設(shè)計(jì)提供了重要的理論參考和技術(shù)手段。

參 考 文 獻(xiàn)

[1]Bendsoe M P, Kikuchi N. Generating optimal topologies in structural design using a homogenization method[J]. Computer Methods in Applied Mechanics and Engineering, 1988, 71(2): 197-224.

[2]Luo J, Gea H C. Optimal stiffener design for interior sound reduction using a topology optimization based approach[J], Journal of Vibration and Acoustics, 2003, 125:267-273.

[3]Akl W, Baz A. Topology optimization of a plate coupled with acoustic cavity[J]. International Journal of Solids and Structures, 2009, 46: 2060-2074.

[4]郭中澤,陳裕澤,侯 強(qiáng),等. 阻尼材料布局優(yōu)化研究[J]. 兵工學(xué)報(bào), 2007, 28(5): 638-640.

GUO Zhong-ze, CHEN Yu-ze, HOU Qiang, et al. Damping material optimal placement in damping structure design[J]. Acta Armamentarii, 2007, 28(5): 638-640.

[5]楊德慶,柳擁軍,金咸定. 薄板減振降噪的拓?fù)鋬?yōu)化設(shè)計(jì)方法[J]. 船舶力學(xué), 2003, 7(5): 91-96.

YANG De-qing, LIU Yong-jun, JIN Xian-ding, et al. Structural topology optimal design to reduce vibration and noise of thin plate[J]. Journal of Ship Mechanics, 2003, 7(5): 91-96.

[6]李以農(nóng),謝熔爐,王 宜,等. 約束阻尼結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)的進(jìn)化算法[J]. 重慶大學(xué)學(xué)報(bào), 2006, 10(2): 120-125.

LI Yi-nong, XIE Rong-lu, WANG Yi, et al. Topology optimization for constrained layer damping material in structures using ESO method[J]. Journal of Chongqing University(Natural Science Edition), 2006, 10(2): 120-125.

[7]鄭 玲,謝熔爐,王 宜,等. 基于優(yōu)化準(zhǔn)則的約束阻尼材料優(yōu)化配置[J]. 振動(dòng)與沖擊, 2010, 29(11): 156-179.

ZHENG Ling, XIE Rong-lu, WANG Yi, et al. Optimal placement of constrained damping material in structures based on optimality criteria[J]. Journal of Vibration and Shock, 2010, 29(11): 156-179.

[8]李 超,李以農(nóng),施 磊,等. 圓柱殼體阻尼材料布局拓?fù)鋬?yōu)化研究[J]. 振動(dòng)與沖擊, 2012, 31(4): 48-52.

LI Chao, LI Yi-nong, SHI Lei, et al. Topological optimization for placement of damping material on cylindrical shells[J]. Journal of Vibration and Shock, 2012, 31(4): 48-52.

[9]王明旭,陳國(guó)平. 基于變密度方法約束阻尼層動(dòng)力學(xué)性能優(yōu)化[J]. 南京航空航天大學(xué)學(xué)報(bào), 2010, 42(3): 283-287.

WANG Ming-xu, CHEN Guo-ping. Dynamics performance optimization of constrained damping layer using variable density method[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2010, 42(3): 283-287.

[10]呂毅寧,呂振華,趙 波,等. 附加自由阻尼結(jié)構(gòu)系統(tǒng)的有限元分析和拓?fù)鋬?yōu)化設(shè)計(jì)方法研究[J]. 計(jì)算力學(xué)學(xué)報(bào), 2012, 29(2).

Lü Yi-ning, Lü Zhen-hua, ZHAO Bo, et al. A finite element analysis and topology optimization method for structures with free damping[J]. Chinese Journal of Computational Mechanics, 2012, 29(2).

[11]Duhring M B, Sigmund O. Acoustic design by topology optimization[J]. Journal of Sound and Vibration, 2008, 317: 557-575.

[12]Denli H, Sun J Q. Structural-acoustic optimization of sandwich structures with cellular cores for minimum sound radiation[J]. Journal of Sound and Vibration, 2007, 301: 93-105.

[13]劉天雄,華宏星,陳兆能,等. 約束層阻尼板的有限元建模研究[J]. 機(jī)械工程學(xué)報(bào), 2002, 38(4): 108-114.

LIU Tian-xiong, HUA Hong-xing, CHEN Zhao-neng, et al. Study on the model of finite element of constrained layer damping plate[J]. Chinese Journal of Mechanical Engineering, 2002, 38(4): 108-114.

[14]龔培康. 汽車(chē)拖拉機(jī)有限元法基礎(chǔ)[M]. 北京:機(jī)械工業(yè)出版社, 1995.

[15]劉寶山,趙國(guó)忠,顧元憲. 復(fù)合材料層合板結(jié)構(gòu)振動(dòng)聲輻射優(yōu)化[J]. 振動(dòng)與沖擊, 2008, 27(12): 31-35.

LIU Bao-shan, ZHAO Guo-zhong, GU Yuan-xian, et al. Optimization of acoustic radiation caused by structural vibration of composite laminated plates[J]. Journal of Vibration and Shock, 2008, 27(12): 31-35.

[16]袁國(guó)清. 基于聲輻射模態(tài)模型求解聲學(xué)靈敏度[D]. 鎮(zhèn)江:江蘇大學(xué), 2009.