張彩霞，沙云東

(沈陽(yáng)航空航天大學(xué) a.能源與環(huán)境學(xué)院；b.航空航天工程學(xué)部(院)，沈陽(yáng) 110136)

?

沖擊載荷作用下約束阻尼結(jié)構(gòu)振動(dòng)和聲輻射特性研究

張彩霞a，沙云東b

(沈陽(yáng)航空航天大學(xué) a.能源與環(huán)境學(xué)院；b.航空航天工程學(xué)部(院)，沈陽(yáng) 110136)

利用部分敷設(shè)約束阻尼減小振動(dòng)和噪聲輻射對(duì)沖擊載荷作用下薄壁板的結(jié)構(gòu)振動(dòng)和聲輻射特性進(jìn)行了研究。對(duì)一四邊簡(jiǎn)支矩形板施加同峰值的三角和半周期正弦形的沖擊載荷，通過(guò)對(duì)時(shí)域振動(dòng)速度、振動(dòng)速度功率譜密度和聲輻射功率級(jí)的求解，研究?jī)煞N載荷對(duì)結(jié)構(gòu)振動(dòng)和噪聲輻射的影響。然后在板表面中心位置部分敷設(shè)約束阻尼并逐漸增加敷設(shè)面積，研究敷設(shè)面積大小對(duì)減振降噪的影響。結(jié)果表明，同峰值的半周期正弦形沖擊載荷使結(jié)構(gòu)產(chǎn)生更大的振動(dòng)及噪聲輻射；敷設(shè)約束阻尼面積越大結(jié)構(gòu)聲輻射功率級(jí)越小，但面積逐漸增大，噪聲的衰減率卻逐漸變小。

沖擊載荷；薄壁結(jié)構(gòu)；約束阻尼；減振降噪

薄壁結(jié)構(gòu)在航空航天領(lǐng)域應(yīng)用廣泛，如飛行器蒙皮、發(fā)動(dòng)機(jī)燃燒室等。它們常常受到隨時(shí)間變化的沖擊載荷作用而產(chǎn)生強(qiáng)烈的結(jié)構(gòu)振動(dòng)并輻射較大的噪聲，縮短工件壽命、影響工作人員健康。因此對(duì)沖擊載荷作用下航空薄壁結(jié)構(gòu)的減振降噪研究具有重要意義。粘彈性材料具有大阻尼特性，可以很好地抑制結(jié)構(gòu)產(chǎn)生的振動(dòng)和噪聲，因此阻尼減振降噪技術(shù)被廣泛應(yīng)用于機(jī)械結(jié)構(gòu)的振動(dòng)和噪聲控制領(lǐng)域[1-7]。約束阻尼的減振降噪效果往往要優(yōu)于自由阻尼，因此國(guó)內(nèi)外的學(xué)者更專(zhuān)注于對(duì)約束阻尼減振降噪技術(shù)的研究[3-7]。李恩奇等[8]根據(jù) Hamilton原理推導(dǎo)出了約束層阻尼梁運(yùn)動(dòng)和邊界條件方程，經(jīng) Laplace 變換后再引入狀態(tài)向量，建立系統(tǒng)的空間狀態(tài)方程。李九天等[9]基于分布參數(shù)傳遞函數(shù)的方法研究了粘彈性材料隨機(jī)性對(duì)約束層阻尼結(jié)構(gòu)動(dòng)力學(xué)特性的影響。梁超鋒等[10]由材料阻尼的損耗因子出發(fā)推導(dǎo)得到了約束阻尼彈性階段內(nèi)結(jié)構(gòu)的阻尼矩陣。Theofanis等[11]用有限元的方法計(jì)算了粘彈性?shī)A層板的阻尼，并預(yù)測(cè)了結(jié)構(gòu)的固有頻率、模態(tài)振型和模態(tài)損耗因子，最后研究了粘彈性層和約束層厚度對(duì)損耗因子和固有頻率的影響。Y.S.Chen等[12]通過(guò)實(shí)驗(yàn)對(duì)約束阻尼夾層結(jié)構(gòu)的損耗因子進(jìn)行了研究，觀察到了溫度和頻率的變化對(duì)損耗因子的影響。

本文對(duì)一四邊簡(jiǎn)支的矩形板分別施加同峰值的三角和半周期正弦形沖擊載荷，通過(guò)完全瞬態(tài)分析法求解板表面的時(shí)域振動(dòng)速度和振動(dòng)速度功率譜密度，觀察它們?cè)跁r(shí)域和頻域的振動(dòng)特性及能量分布的差異，再采用耦合有限元和邊界元方法求解它們的聲輻射功率級(jí)。然后在板表面中心位置部分敷設(shè)約束阻尼并逐漸增加敷設(shè)面積，研究敷設(shè)面積大小對(duì)減振降噪的影響。

1 結(jié)構(gòu)振動(dòng)及聲輻射理論

1.1 耦合有限元和間接邊界元分析

在物理坐標(biāo)系下，耦合有限元和間接邊界元矩陣方程為[13]

(1)

式中，Ks是結(jié)構(gòu)剛度矩陣,Ms是結(jié)構(gòu)質(zhì)量矩陣，Cs是結(jié)構(gòu)的阻尼矩陣，C是幾何耦合矩陣，H(ω)是間接邊界元影響矩陣，u是節(jié)點(diǎn)位移，μ是節(jié)點(diǎn)壓力跳動(dòng)，F(xiàn)s是結(jié)構(gòu)激勵(lì)向量，F(xiàn)A是流體激勵(lì)向量。

1.2 結(jié)構(gòu)聲輻射理論

流體是連續(xù)的介質(zhì)，因此認(rèn)為結(jié)構(gòu)表面振動(dòng)速度和結(jié)構(gòu)表面流體的振動(dòng)速度相同[14]。若結(jié)構(gòu)表面任意一點(diǎn)的聲壓為P(y,ω)，那么

(2)

式中，ρ是空氣密度，r是結(jié)構(gòu)表面任意兩點(diǎn)x和y間的距離，c是聲速。結(jié)構(gòu)表面的聲強(qiáng)可以表示為

(3)

式中，Re是雷諾數(shù)，T是時(shí)間。結(jié)構(gòu)表面聲輻射功率的表達(dá)式為

(4)

2 沖擊載荷算例

2.1 不同沖擊載荷下的響應(yīng)

研究對(duì)象為一四邊簡(jiǎn)支的矩形薄板，長(zhǎng)、寬和高分別為1 000 mm、600 mm和3 mm，材料密度為1 763 kg/m3，泊松比為0.32，楊氏模量為73.11 GPa，幾何模型如圖1所示。有限元模型采用shell181單元。對(duì)薄板進(jìn)行模態(tài)分析得到前8階振動(dòng)頻率如表1所示。

圖1 薄板幾何模型

通常采用矩形、三角形和半周期正弦形載荷模擬沖擊載荷進(jìn)行研究[15]，本文選擇后兩者。在板的中心位置分別施加峰值為100 N豎直向下的三角和半周期正弦2種形式的沖擊載荷,如圖2所示，加載時(shí)間為0.5 s。利用完全瞬態(tài)分析方法得到板表面隨時(shí)間變化的振動(dòng)速度，求振動(dòng)速度的功率譜密度，最后以板表面頻域振動(dòng)速度為邊界條件構(gòu)造模型，在距離結(jié)構(gòu)1 m處建立平面聲場(chǎng)網(wǎng)格，對(duì)薄板聲學(xué)響應(yīng)進(jìn)行求解。圖3～5為薄板中心點(diǎn)振動(dòng)及聲輻射情況，圖6為聲場(chǎng)云圖。

表1 薄板前8階模態(tài)頻率

圖2 兩種形式載荷時(shí)域曲線圖

圖3 薄板中心點(diǎn)的振動(dòng)速度時(shí)域曲線

圖4 薄板中心點(diǎn)的速度功率譜密度

圖5 薄板中心點(diǎn)的聲功率級(jí)

圖6 兩種形式載荷下聲場(chǎng)云圖

對(duì)比發(fā)現(xiàn):(1)載荷加載0.5 s到1 s過(guò)程中，三角形載荷的振動(dòng)速度可以分為正向和負(fù)向兩個(gè)階段，在峰值處發(fā)生變換，在各階段內(nèi)振動(dòng)速度變化較小。半周期正弦形載荷的振動(dòng)速度在正向到負(fù)向變化過(guò)程中出現(xiàn)中間階段,即正負(fù)向交替出現(xiàn)。自由振動(dòng)1 s后階段，半周期正弦形載荷的起始振動(dòng)速度較大、衰減較快；(2)薄板施加半周期正弦形沖擊載荷后振動(dòng)能量主要集中在特征頻率處，而施加三角形沖擊載荷后在各頻率處能量分布比較均勻，載荷形式的改變并沒(méi)有改變結(jié)構(gòu)的特征頻率；(3)半周期正弦沖擊載荷聲輻射明顯大于三角形沖擊載荷，同時(shí)對(duì)環(huán)境造成污染的噪聲主要集中在低頻率處；(4)由聲場(chǎng)云圖可以發(fā)現(xiàn)距離載荷點(diǎn)越近的場(chǎng)點(diǎn)聲壓越大，并且相同場(chǎng)點(diǎn)施加半周期正弦形沖擊載荷時(shí)聲壓較大?？梢园l(fā)現(xiàn)，施加同峰值半周期正弦形沖擊載荷后薄板的振動(dòng)更劇烈，對(duì)環(huán)境造成的噪聲污染問(wèn)題更為突出，下面主要針對(duì)半周期正弦形的沖擊載荷進(jìn)行研究。

2.2 敷設(shè)不同面積約束阻尼

在薄板的中心位置敷設(shè)不同面積的約束阻尼層，設(shè)計(jì)覆蓋面積分別為6.67%、13.33%和26.67%，如圖7所示。

有限元模型層與層之間采用節(jié)點(diǎn)共用的方法，基層和約束層分別采用向上、向下偏心的shell181單元，粘彈性層采用solid185單元，厚度分別為3 mm和15 mm。粘彈性層材料屬性為楊氏模量270 MPa，損耗因子0.1，泊松比0.49，密度1 780 kg/m3,不考慮粘彈性材料的彈性模量和損耗因子受頻率和溫度的影響，取經(jīng)驗(yàn)值。施加圖2(b)所示的半周期正弦形沖擊載荷。敷設(shè)不同面積的約束阻尼后，結(jié)構(gòu)的模態(tài)頻率如表2所示。薄板中心點(diǎn)時(shí)域振動(dòng)速度曲線和振動(dòng)速度功率譜密度如圖8和圖9所示。聲輻射功率級(jí)曲線和距薄板1 m處平面聲場(chǎng)云圖如圖10和圖11所示。

圖7 部分敷設(shè)約束阻尼幾何示意圖

Hz

圖8 中心點(diǎn)振動(dòng)速度時(shí)域曲線

圖9 中心點(diǎn)振動(dòng)速度功率譜密度

圖10 中心點(diǎn)輻射聲功率級(jí)

圖11 約束阻尼不同敷設(shè)面積聲場(chǎng)云圖

研究發(fā)現(xiàn),低頻范圍內(nèi)敷設(shè)約束阻尼面積越大結(jié)構(gòu)的模態(tài)頻率越大。敷設(shè)面積增大，結(jié)構(gòu)時(shí)域振動(dòng)速度隨之減小，但振動(dòng)速度趨勢(shì)并無(wú)變化。敷設(shè)面積從6.67%增加到26.67%時(shí)，在100 Hz頻率范圍內(nèi)減小了一階模態(tài)頻率，聲輻射功率級(jí)降低了10 dB，減振效果明顯；敷設(shè)6.67%面積的約束阻尼層后結(jié)構(gòu)聲輻射功率級(jí)從127 dB降到110 dB；當(dāng)敷設(shè)面積為6.67%時(shí)聲輻射功率級(jí)為110 dB，敷設(shè)面積增加一倍聲輻射功率級(jí)減小了5 dB；敷設(shè)面積增加三倍聲輻射功率級(jí)減小了10 dB。可以發(fā)現(xiàn)，敷設(shè)面積越大結(jié)構(gòu)聲輻射功率級(jí)越小，但面積逐漸增大，噪聲的衰減率卻逐漸變小。由圖11可以清楚地看到隨著敷設(shè)面積的增加相同場(chǎng)點(diǎn)的聲壓逐漸減小。

3 結(jié)論

通過(guò)以上研究得出如下結(jié)論：

(1)同峰值的半周期正弦形沖擊載荷與三角形載荷相比，前者在時(shí)域內(nèi)結(jié)構(gòu)振動(dòng)速度振幅較大，在頻域內(nèi)振動(dòng)能量分布更為集中，主要在特征頻率處，且結(jié)構(gòu)產(chǎn)生的聲輻射功率更強(qiáng)。

(2)沖擊載荷下薄壁結(jié)構(gòu)振動(dòng)能量和噪聲輻射主要集中在低頻范圍內(nèi)，因此針對(duì)沖擊載荷下薄壁結(jié)構(gòu)的減振降噪應(yīng)該關(guān)注對(duì)低頻范圍的研究。

(3)敷設(shè)約束阻尼面積越大結(jié)構(gòu)聲輻射功率級(jí)越小，但面積逐漸增大，噪聲的衰減率卻逐漸變小。因此，工程中敷設(shè)約束阻尼的面積應(yīng)適當(dāng)，過(guò)大的面積可能造成材料浪費(fèi)，但是得到的減振效果卻并不明顯。

[1]Xisheng Cao，Hans Peter Mlejnek.Computational prediction and redesign for viscoelastically damped structures[J].Computer Methods in Appllied Mechanics and Engineering，1995(125)：1-16.

[2]Semyung Wand，Jeawon Lee.Acoustic design sensitivity analysis and optimization for reduced exterior noise[J].AIAA Journal，2001，39(4)：1237-1245.

[3]Mohammed A，Al Ajmi,Khaled A，et al.Graded finite element modeling of constrained layer damping treatments with functionally graded[J].Journal of Intelligent Material Systems and Structures,2008(19):469-474.

[4]R Moreira，J D Rodrigues.Constrained damping layer treatments:finite element modeling[J].Journal of Vibration And Control，2004(10)：575-595.

[5]H Zheng,X M Tan,C Cai.Damping analysis of beams covered with multiple PCLD patches[J].International Journal of Mechanical Sciences,2006(48):1371-1383.

[6]Douglas Be,Yang J C S.Transverse compressional damping in the vibratory response of elastic-viscoelastic beams[J].AIAA Journal,1978(16):925-930.

[7]Marcelin J L,Trompette P,Smati A.Optimal constrained layer damping with partial coverage[J].Finite Elements in Analysis and Design,1992(12):273-280.

[8]李恩奇,唐國(guó)金,李道奎,等.局部覆蓋約束層阻尼梁動(dòng)力學(xué)問(wèn)題的解析解.振動(dòng)與沖擊,2007,26(5):85-89.

[9]李九天,唐國(guó)金,李恩奇，等.粘彈性材料參數(shù)隨機(jī)性對(duì) PCLD 板動(dòng)力學(xué)特性的影響[J].強(qiáng)度與環(huán)境,2008,35(1):21-28.

[10]梁超鋒,歐進(jìn)萍.結(jié)構(gòu)阻尼與材料阻尼的關(guān)系[J].地震工程與工程振動(dòng),2006,26(1)：49-52.

[11]Theofanis S，Plagianakos,Dimitris A,et al.High order layerwise mechanics and finite element for the damped dynamic characteristics of sandwich composite beams[J].International Journal of Solids and Structures,2004(41):6853-6871.

[12]Y S Chen,T J Hsu,S I Chen.Vibration damping characteristics of laminated steel sheet[J].Metallurgical and Materials Transactions A,1991(22A):653-656.

[13]徐張明 沈榮瀛 華宏星.利用FEM/IBEM計(jì)算流體介質(zhì)中的殼體的結(jié)構(gòu)聲耦合問(wèn)題[J].振動(dòng)工程學(xué)報(bào)，2002，15(3)：363-367.

[14]伊崗,陳花玲,陳天寧.薄板低頻聲輻射效率的研究[J].西安交通大學(xué)學(xué)報(bào),1999,33(3):108-110.

[15]張濤,劉土光,趙耀,等.加筋板彈性大撓度的沖擊響應(yīng)分析[J].爆炸與沖擊,2002,22(4):301-307.

(責(zé)任編輯：宋麗萍 英文審校：隋華)

Vibration and sound radiation analysis of plates with constrained layer damping patches under impact loading

ZHANG Cai-xiaa，SHA Yun-dongb

(a.College of Energy and Environment;b.Faculty of Aerospace Engineering,Shenyang Aerospace University,Shenyang 110136,China)

The vibration and sound radiation characteristics of thin-walled structure under impact loading were analyzed by applying some constrained layer damping patches to reduce the vibration and noise radiation.Through same volume of triangle and half-cycle sinusoidal impact loading on a simply supported rectangular plate,the vibration velocity,velocity power spectral density and acoustic radiation power levels were studied to find out the effects of these two loadings on vibration and noise radiation.Then the center of the plate was covered with some constrained layer damping patches and the laying area was increased gradually to find out the effect.The results show that the half-cycle sinusoidal impact loading resulted in greater vibration and noise radiation.The larger the constrained area,the smaller the noise radiation,while the noise attenuation rate was gradually getting smaller.

impact loading;thin-walled structure;constrained layer damping;reduction of vibration and noise radiation

2013-11-06

張彩霞(1988-)，女，遼寧遼陽(yáng)人，碩士研究生，主要研究方向：飛行器動(dòng)力學(xué)環(huán)境工程及控制技術(shù)研究，E-mail：617321997@qq.com；沙云東(1966-)，男，黑龍江阿城人，教授，主要研究方向：航空發(fā)動(dòng)機(jī)強(qiáng)度，振動(dòng)及噪聲，E-mail:Ydsha2003@vip.sina.com。

2095-1248(2015)05-0037-06

F416.5

A

10.3969/j.issn.2095-1248.2015.05.004