約束阻尼結(jié)構(gòu)阻尼效果的有限元預(yù)測(cè)方法研究

伍先俊

(中國(guó)科學(xué)院 聲學(xué)研究所，北京 100190)

機(jī)械系統(tǒng)振動(dòng)消除的方法除減小源振動(dòng)外，還有隔振和阻尼消振。阻尼消振是在結(jié)構(gòu)覆蓋阻尼層來(lái)約束結(jié)構(gòu)振動(dòng)，使得結(jié)構(gòu)振動(dòng)能量轉(zhuǎn)換為阻尼材料變形的熱能，因而消耗掉振動(dòng)能量的方法。阻尼消振技術(shù)在汽車、飛機(jī)、鐵路運(yùn)輸、船舶中都有應(yīng)用，特別是結(jié)構(gòu)的共振區(qū)，阻尼材料的貼覆能夠極大減小振動(dòng)，降低噪聲。

常用約束阻尼結(jié)構(gòu)，是將黏彈性阻尼材料粘合在本體金屬板和剛度較大的約束層（通常是金屬板）之間，當(dāng)結(jié)構(gòu)彎曲變形時(shí)，本體金屬板與約束層產(chǎn)生相對(duì)滑移，黏彈性阻尼材料產(chǎn)生剪切應(yīng)變使一部分機(jī)械能損耗的結(jié)構(gòu)。阻尼消振結(jié)構(gòu)有無(wú)約束面板的阻尼消振結(jié)構(gòu)和含約束面板的阻尼消振結(jié)構(gòu)，阻尼材料的厚度通常取本體板厚度的2倍左右，選用的約束面板層彈性模量通常遠(yuǎn)大于阻尼層的彈性模量，厚度為本體金屬的1/2～1/4。

對(duì)于無(wú)約束面板的兩層結(jié)構(gòu)和含阻尼面板的三層結(jié)構(gòu)阻尼目前有一些針對(duì)梁的簡(jiǎn)單預(yù)報(bào)方法[1]，Torvik P J等[2]用復(fù)模量描述黏彈材料的本構(gòu)關(guān)系，并用模態(tài)應(yīng)變能法對(duì)夾層板進(jìn)行阻尼振動(dòng)分析；文獻(xiàn)[3]分析提出一種基于GHM模型（由格羅斯曼和哈特(Gross man&Hart,1986)、哈特和莫爾(Hart&Moore，1990)等共同創(chuàng)立的）的有限元計(jì)算方法，GHM模型實(shí)質(zhì)是測(cè)量了阻尼材料的頻域曲線后擬合出阻尼材料的傳遞函數(shù)，文獻(xiàn)[4]提出了一種基于GHM模型的有限元理論實(shí)現(xiàn)方法，目前的主流商業(yè)軟件對(duì)諧波分析的阻尼處理還是基于復(fù)模量的理論，文獻(xiàn)[5]根據(jù)復(fù)模量的理論對(duì)多層阻尼圓板的損耗因子進(jìn)行了分析，選取它們的倒數(shù)和為目標(biāo)函數(shù)進(jìn)行有限元優(yōu)化設(shè)計(jì)，不考慮自然頻率的導(dǎo)數(shù)和對(duì)于寬帶激勵(lì)不太適用。本文將對(duì)采用有限元法的阻尼因子計(jì)算法進(jìn)行研究，基于模態(tài)應(yīng)變能提出一種新的復(fù)合阻尼的快速計(jì)算方法，同時(shí)結(jié)合統(tǒng)計(jì)能量分析理論提出一種不同的結(jié)構(gòu)阻尼優(yōu)化目標(biāo)函數(shù)，并提出基于該函數(shù)的有限元優(yōu)化設(shè)計(jì)方法。通過(guò)經(jīng)典算例的仿真討論，驗(yàn)證了基于模態(tài)應(yīng)變能的阻尼計(jì)算方法的正確性。對(duì)某含約束面板的阻尼層厚度進(jìn)行了優(yōu)化設(shè)計(jì)，最后根據(jù)諧波響應(yīng)分析驗(yàn)證了基于新的阻尼效果評(píng)價(jià)函數(shù)的優(yōu)化計(jì)算的合理性。

1 復(fù)雜結(jié)構(gòu)有限元計(jì)算理論和步驟

對(duì)阻尼結(jié)構(gòu)的阻尼效果分析有諧波分析法和模態(tài)分析法，其中諧波分析法可以將阻尼參數(shù)代入諧波計(jì)算中，直接獲得有阻尼響應(yīng)，是精確的方法。但是諧波計(jì)算需要進(jìn)行整體結(jié)構(gòu)所有自由度的響應(yīng)，計(jì)算耗時(shí)長(zhǎng)，不適于優(yōu)化分析。另外一種方法是基于模態(tài)應(yīng)變能的阻尼計(jì)算方法，即先計(jì)算出各階模態(tài)，通過(guò)材料參數(shù)和模態(tài)變形能計(jì)算獲得各階模態(tài)阻尼，這種方法比諧波分析計(jì)算量小很多。下面將先介紹諧波分析矩陣阻尼模型，然后給出基于模態(tài)應(yīng)變能的阻尼計(jì)算方法。

1.1 諧波分析阻尼模型

對(duì)于諧波分析的響應(yīng)，阻尼矩陣可以寫為[6]

式中：符號(hào)含義如下：[C]結(jié)構(gòu)總阻尼矩陣；α阻尼質(zhì)量乘子，[M]質(zhì)量矩陣，β阻尼剛度乘子，g常量結(jié)構(gòu)阻尼因子，Ω激勵(lì)頻率，[K]剛度矩陣，Nm有阻尼材料個(gè)數(shù)，gj對(duì)應(yīng)于材料j的常量結(jié)構(gòu)阻尼因子，[Kj]基于材料的剛度矩陣，Ne具有阻尼輸入的單元，[Ck]單元阻尼矩陣，Ng有科式或者陀螺阻尼的單元個(gè)數(shù)，[Gl]科式或者陀螺阻尼矩陣（科式或者陀螺阻尼在一般平動(dòng)結(jié)構(gòu)中不存在，所以本文不予考慮）。

諧波響應(yīng)分析的方法有完全諧波分析、縮減諧波分析、模態(tài)疊加諧波分析和QR阻尼模態(tài)疊加分析法。在多種材料的結(jié)構(gòu)諧波響應(yīng)分析中，要求計(jì)算支持配比不同材料的βj，gj和Ck，才可以考慮不同材料阻尼對(duì)結(jié)構(gòu)響應(yīng)的影響，對(duì)以上阻尼都支持的有完全諧波分析法或者QR阻尼模態(tài)疊加分析法。

完全諧波分析法計(jì)算量比較大，QR阻尼模態(tài)疊加分析法可以減小計(jì)算量，因此本文采用后一方法。

根據(jù)動(dòng)力學(xué)方程

對(duì)于QR阻尼模態(tài)疊加分析法，[K]、[C]可以是非角對(duì)稱矩陣，QR阻尼模態(tài)疊加分析法將[K]矩陣重新排列分出對(duì)稱部分和非對(duì)稱部分，通過(guò)略去[C]矩陣和[K]矩陣的非對(duì)稱部分，求得特征值和特征向量，通過(guò)以下公式進(jìn)行坐標(biāo)變換

[Φ]為對(duì)質(zhì)量歸一化模態(tài)振型向量矩陣，{y}為模態(tài)坐標(biāo)，通過(guò)以上公式可獲得以下式子

盡管 [Φ]T[C][Φ]和 [Φ]T[Kunsys][Φ]為對(duì)角陣，但是模態(tài)坐標(biāo)數(shù)目得到減少，計(jì)算量得到降低。[Λ2]為對(duì)角特征值平方，[Kunsys]為剛度矩陣非對(duì)稱部分。

1.2 模態(tài)阻尼計(jì)算

對(duì)于具有多種材料，各材料阻尼因子不相同，可以采用基于變形能計(jì)算出一個(gè)有效模態(tài)阻尼因子，在擴(kuò)展模態(tài)時(shí)計(jì)算單元結(jié)果(采用APDL命令MXPAND,NMODE,FREQB,FREQE,Elcalc,SIGNIF，而Elcalc取1表示計(jì)算單元結(jié)果)

其中：N為材料個(gè)數(shù)，βj為材料阻尼因子，Ej為對(duì)應(yīng)材料j的模態(tài)變形能，Ej={φi}[Kj]{φi} ，{φi}為振型函數(shù)，[Kj]對(duì)應(yīng)于材料j的剛度矩陣。

1.3 阻尼評(píng)價(jià)函數(shù)

本文基于統(tǒng)計(jì)能量分析的理論，得出一種針對(duì)寬帶激勵(lì)的阻尼情形，比較適用于通用場(chǎng)合。受隨機(jī)激勵(lì)單自由度振動(dòng)物體的振速寫為[7]

其中：Sf為輸入力自功率譜，ωi為單自由振子模態(tài)頻率，M為被分析結(jié)構(gòu)質(zhì)量，ηi為單自由振子模態(tài)阻尼損耗因子ηi。

多階模態(tài)為多個(gè)振動(dòng)物體，其總的響應(yīng)為

對(duì)于多自由度系統(tǒng)，ωi為各模態(tài)頻率，ηi為模態(tài)阻尼損耗因子ηi。

該值越小阻尼效果越好，這種評(píng)價(jià)與過(guò)去的各模態(tài)阻尼平均或者倒數(shù)和都不相同[5]。

1.4 有限元計(jì)算關(guān)鍵過(guò)程

有限元分析采用單層實(shí)體殼單元代替梁?jiǎn)卧虚g層采用實(shí)體殼單元，覆蓋層采用面殼體單元。各部分單元類型、材料類型、阻尼設(shè)置的APDL命令見(jiàn)表1。

Solsh190為3-D8節(jié)點(diǎn)實(shí)體殼單元，是一種看上去象實(shí)體的殼單元，可以模擬從薄到中等厚度的殼結(jié)構(gòu)，薄殼結(jié)構(gòu)無(wú)需抽中面，直接劃分單元，只需一層Solid-shell單元，并能保證精度。

建立有限元模型，計(jì)算模態(tài)頻率，并擴(kuò)展獲得模態(tài)能量值：

2 阻尼計(jì)算方法和結(jié)果討論

2.1 對(duì)于單層梁計(jì)算結(jié)果的討論

圖1 自由阻尼約束層示意圖

對(duì)于兩層阻尼結(jié)構(gòu)，損耗因子解析計(jì)算式為

算例：模型基礎(chǔ)層采用長(zhǎng)0.5 m，截面積為H×B=7 mm×7 mm的鋁梁，阻尼材料厚度分別從0.002變化到0.01，阻尼材料彈性模量7×106，密度1 000，泊松比0.4，結(jié)構(gòu)阻尼因子1，兩端自由。（為了只比較約束方向模態(tài)阻尼，因此約束梁非彎曲方向振動(dòng)。）

表2是對(duì)一個(gè)梁結(jié)構(gòu)第1階模態(tài)采用有限元方法和解析方法計(jì)算對(duì)比，兩層阻尼梁解析方法忽略了材料泊松比、梁寬度以及阻尼層質(zhì)量的影響，在約束層厚度增大時(shí)計(jì)算誤差增大。

有限元計(jì)算阻尼層厚度大的結(jié)構(gòu)阻尼因子與解析法差別較大，因此工程中使用解析式（3）需要格外注意。

2.2 對(duì)于含約束面板三層阻尼梁計(jì)算結(jié)果討論

對(duì)于含約束的三層阻尼結(jié)構(gòu)如下

圖2 約束阻尼層示意圖

表1 各部分單元類型

對(duì)于三層梁結(jié)構(gòu)，有以下解析計(jì)算公式

對(duì)于兩端自由梁自然頻率計(jì)算如下

對(duì)于兩端自由梁，有λ1l=4.73，λ2l=7.85，…λil≈(2i+1)π/2，ρ為梁材料密度，A為梁截面積，I

模型基礎(chǔ)層采用長(zhǎng)0.5 m，截面為H×B=7 mm×7 mm的鋁梁，阻尼材料厚度0.003，阻尼材料彈性模量7×106Pa，密度1 000 kg/m3，泊松比0.4，其結(jié)構(gòu)阻尼因子為1，兩端自由，約束層為鋼，厚度0.000 5。（為了只比較約束方向模態(tài)阻尼，因此約束梁非彎曲方向振動(dòng)。）

由于約束阻尼層與波長(zhǎng)有關(guān)系，所以計(jì)算前5階模態(tài)（公式（4）計(jì)算的模態(tài)頻率為：147.08 Hz，364.29 Hz，697.31 Hz，1 136.1 Hz，1 680.2 Hz）阻尼大小對(duì)比，表3是對(duì)一個(gè)含約束面板三層阻尼梁結(jié)構(gòu)采用解析法和有限元法的計(jì)算結(jié)果對(duì)比。

有限元計(jì)算高階模態(tài)的阻尼因子與解析法差別較大，因此工程中使用解析式（4）需要格外注意。

為了驗(yàn)證有限元計(jì)算結(jié)果的正確性，采用另外一種方法獲得阻尼因子作為驗(yàn)證。即類似于通過(guò)實(shí)驗(yàn)獲得傳遞函數(shù)曲線，通過(guò)傳遞函數(shù)曲線辨識(shí)阻尼的方法。對(duì)于集中點(diǎn)力情況同時(shí)計(jì)算諧波響應(yīng)，從而獲得傳遞函數(shù)曲線。取距離梁中部0.06 m的位置計(jì)算加速度傳遞函數(shù)，傳遞函數(shù)曲線見(jiàn)圖3。圖中給出了用最小二乘法[8]與以及原始計(jì)算獲得的傳遞函數(shù)實(shí)部和虛部的對(duì)比。

圖3 計(jì)算獲得傳遞函數(shù)和最小二乘擬合傳遞函數(shù)對(duì)比

表2 解析方法和有限元計(jì)算方法對(duì)比

表3 解析方法和有限元計(jì)算方法對(duì)比

采用最小二乘法迭代法[8]識(shí)別出模態(tài)頻率、阻尼因子和解析解計(jì)算的對(duì)應(yīng)數(shù)值見(jiàn)表4，可以看出采用有限元傳遞函數(shù)辨識(shí)出的各階模態(tài)頻率與解析解符合很好，同時(shí)各階模態(tài)的阻尼與采用有限元及公式（1）計(jì)算結(jié)果符合很好，因此驗(yàn)證了采用有限元及公式（1）計(jì)算結(jié)果的正確性。

說(shuō)明模態(tài)分析得到阻尼與諧波分析得到的阻尼比是對(duì)應(yīng)的，驗(yàn)證了分析結(jié)果的正確性。

2.3 板的計(jì)算結(jié)果和優(yōu)化

具體優(yōu)化方法參考文獻(xiàn)[9]。

優(yōu)化解：hd1=0.001，hd2=0.000 69，mas=1.595 5 kg，rs=0.396。為了檢驗(yàn)?zāi)B(tài)阻尼計(jì)算結(jié)果的正確性，采用完全諧波響應(yīng)計(jì)算優(yōu)化結(jié)果與非優(yōu)化結(jié)果響應(yīng)曲線對(duì)比。

hd1=0.001，hd2=0.000 69，隨機(jī)取5個(gè)點(diǎn)激勵(lì)，對(duì)平板平均速度做平均，優(yōu)化前后響應(yīng)譜對(duì)比如下圖4所示。

表4 傳遞函數(shù)分析辨識(shí)獲得各階模態(tài)頻率和阻尼因子

對(duì)一個(gè)平板進(jìn)行優(yōu)化設(shè)計(jì)，由于平板模型沒(méi)有解析方法可以計(jì)算，本文先基于模態(tài)應(yīng)變能的阻尼計(jì)算方法進(jìn)行優(yōu)化設(shè)計(jì)，然后采用完全諧波分析法進(jìn)行驗(yàn)證。鋁板四端固支，鋁平板尺寸：0.5 m×0.5 m×0.007 m，阻尼材料彈性模量7×106Pa，密度1 000，泊松比0.4，其結(jié)構(gòu)阻尼因子為1，約束層為鋼。

初始厚度為：APDL設(shè)計(jì)變量為阻尼層厚度hd1，約束層鋼板厚度hd2；阻尼層和約束板總質(zhì)量限制條件為小于1.6 kg，即以下APDL代碼中的狀態(tài)變量mas小于1.6，以下APDL的目標(biāo)函數(shù)為rs。初始設(shè)計(jì)hd1=0.002 5 m，hd2=0.000 5 m，而初始設(shè)計(jì)目標(biāo)值rs=0.486；

對(duì)于優(yōu)化計(jì)算，關(guān)鍵優(yōu)化代碼如下，定義兩個(gè)優(yōu)化變量：

圖4 約束阻尼層平板優(yōu)化前后諧波響應(yīng)譜對(duì)比

從圖4可以看出，通過(guò)合理分配阻尼板約束層和非約束層厚度，目標(biāo)值獲得了1 dB左右的降噪效果。計(jì)算兩條曲線下的陰影面積，面積相比按分貝取對(duì)數(shù)，結(jié)果為減振1.1 dB，這與目標(biāo)函數(shù)的變化基本對(duì)應(yīng)，說(shuō)明采用簡(jiǎn)化方法的合理性。但是采用簡(jiǎn)化方法優(yōu)化計(jì)算時(shí)長(zhǎng)有大幅度降低：一次完整諧波法需要時(shí)長(zhǎng)80 s；采用模態(tài)分析法計(jì)算阻尼效果評(píng)價(jià)函數(shù)法需要時(shí)長(zhǎng)18 s，節(jié)省時(shí)間量為77.5%；預(yù)計(jì)對(duì)于復(fù)雜大模型的計(jì)算，優(yōu)化后速度可提高更多。

3 結(jié)語(yǔ)

本文進(jìn)行了一種基于有限元的阻尼因子計(jì)算方法，提出總響應(yīng)最小的模態(tài)阻尼優(yōu)化計(jì)算方法，對(duì)三層平板阻尼鋪設(shè)進(jìn)行了計(jì)算，并采用諧波響應(yīng)對(duì)優(yōu)化結(jié)果的準(zhǔn)確性進(jìn)行了檢驗(yàn)。該方法基于商業(yè)軟件ANSYS提出，具有比較好的工程可操作性，可供實(shí)際阻尼優(yōu)化設(shè)計(jì)參考。