基于ACTRAN的結(jié)構(gòu)振動(dòng)聲輻射故障檢測(cè)

張春良1，姚明鏡2，岳 夏1，朱厚耀1，齊觀壇1

（1.廣州大學(xué)機(jī)械與電氣工程學(xué)院，廣東廣州 510006；2.成都理工大學(xué)工程技術(shù)學(xué)院，四川樂山 614000）

文章對(duì)隨機(jī)激勵(lì)下結(jié)構(gòu)振動(dòng)聲輻射問題進(jìn)行研究，利用ACTRAN軟件對(duì)簡(jiǎn)單結(jié)構(gòu)體的振動(dòng)聲輻射問題進(jìn)行了系統(tǒng)的分析，運(yùn)用有限元－無限元方法對(duì)不同監(jiān)測(cè)點(diǎn)進(jìn)行定義加載，分析了聲輻射聲壓與振動(dòng)幅值的關(guān)系，同時(shí)，利用實(shí)驗(yàn)進(jìn)行分析驗(yàn)證.結(jié)果證明：隨著激勵(lì)頻率的增大，結(jié)構(gòu)體振動(dòng)幅值與聲輻射聲壓值變化呈正相關(guān)關(guān)系，即激振頻率越大，振動(dòng)幅值與聲壓值的變化也變大.同時(shí)，有裂紋的結(jié)構(gòu)體裂紋處的振動(dòng)幅值與聲壓值比正常結(jié)構(gòu)體變化非常明顯，從而為運(yùn)用振動(dòng)與聲學(xué)進(jìn)行機(jī)械設(shè)備非接觸故障診斷的研究提供了新的方法.關(guān)鍵詞：聲固耦合；振動(dòng)響應(yīng)；聲輻射響應(yīng)

聲是由振動(dòng)引起的，根據(jù)振動(dòng)系統(tǒng)的大小和形狀，可以抽象出不同的聲源模型，例如，點(diǎn)聲源、線聲源等.可是，工程上有些振動(dòng)系統(tǒng)不能理想化，需要對(duì)系統(tǒng)本身進(jìn)行模擬［1］.

ACTRAN采用有限元－無限元法，能夠分析復(fù)雜結(jié)構(gòu)振動(dòng)的聲輻射，包括聲場(chǎng)和結(jié)構(gòu)的耦合與非耦合［2］.研究結(jié)構(gòu)振動(dòng)的聲輻射，就是研究結(jié)構(gòu)振動(dòng)和聲之間的相互關(guān)系.結(jié)構(gòu)振動(dòng)會(huì)產(chǎn)生聲，聲也會(huì)引起結(jié)構(gòu)振動(dòng)［3－4］.在實(shí)際問題中，結(jié)構(gòu)振動(dòng)會(huì)帶動(dòng)其周圍流體介質(zhì)一起振動(dòng)而產(chǎn)生聲波，聲波也會(huì)產(chǎn)生壓力作用在結(jié)構(gòu)上［5］.結(jié)構(gòu)中受激勵(lì)的部件可以直接向周圍空氣輻射聲波，也可以沿著結(jié)構(gòu)傳播引起其他部件振動(dòng)［6－8］.因此，利用ACTRAN軟件研究結(jié)構(gòu)振動(dòng)輻射聲場(chǎng)問題［9－10］，對(duì)于解決噪聲問題及大型機(jī)械設(shè)備中利用聲音進(jìn)行定位及故障診斷問題都有很大的幫助.

1 聲與結(jié)構(gòu)的耦合方程

結(jié)構(gòu)振動(dòng)的有限元?jiǎng)恿W(xué)方程，不失一般性忽略結(jié)構(gòu)阻尼，可以表示為

其中，x為位移，Ms為質(zhì)量矩陣，Ks為剛度矩陣，F(xiàn)s為結(jié)構(gòu)載荷，F(xiàn)as為聲場(chǎng)對(duì)結(jié)構(gòu)的作用力.

其中，Ns為結(jié)構(gòu)單元形函數(shù)，為聲壓，n為方向余弦.根據(jù)聲學(xué)流體單元的插值關(guān)系，有＝Nape.其中Na為聲學(xué)流體單元的形函數(shù).

聲學(xué)有限元方程為

為了推導(dǎo)過程的簡(jiǎn)明，不失一般性這里不考慮流體速度、吸聲邊界條件，但是假設(shè)存在聲源.

把式（3）的右端項(xiàng)表示為

其中，F(xiàn)a是聲載荷，F(xiàn)sa是結(jié)構(gòu)振動(dòng)產(chǎn)生的.對(duì)式（3）忽略阻尼項(xiàng)后進(jìn)行變換，得

由上式的變換，得

這就是聲與結(jié)構(gòu)的耦合方程，矩陣A是它們的耦合矩陣.從上式可見結(jié)構(gòu)在載荷作用下振動(dòng)會(huì)產(chǎn)生聲波，聲源輻射聲波會(huì)引起結(jié)構(gòu)振動(dòng).

2 基于ACTRAN的振動(dòng)聲學(xué)仿真分析

ACTRAN軟件是由FFT公司開發(fā)的新一代計(jì)算聲學(xué)的工具，具備豐富的單元庫(kù)、材料庫(kù)、邊界條件、解決問題方案、求解器和前后處理接口.本文通過正方體結(jié)構(gòu)對(duì)振動(dòng)聲輻射進(jìn)行分析.

2.1 結(jié)構(gòu)有限元的模型建立及參數(shù)設(shè)置

首先采用Hypermesh中建立好正方體的模型，正方體大小為70 mm×70 mm（圖1），然后導(dǎo)入ACTRAN中進(jìn)行參數(shù)設(shè)置.圖中小正方體材料屬性為鋼材料，大正方體材料屬性為空氣.分別對(duì)結(jié)構(gòu)體進(jìn)行屬性定義，然后加載激勵(lì)進(jìn)行分析.對(duì)于有裂紋的結(jié)構(gòu)體，裂紋長(zhǎng)度為20 mm，深1 mm，寬2 mm，結(jié)構(gòu)體其他部分與正常結(jié)構(gòu)體相同.

圖1 正方體結(jié)構(gòu)模型Fig.1 Model for the cube structure

2.2 正方體結(jié)構(gòu)振動(dòng)與聲輻射仿真分析

響應(yīng)分析求解方法有2種：直接法和模態(tài)法.考慮到計(jì)算精度問題，本次分析采用直接頻率響應(yīng).對(duì)結(jié)構(gòu)體的正常與故障狀態(tài)進(jìn)行分析，加載激勵(lì)為100 Hz～1 500 Hz，每次間隔100 Hz；激振點(diǎn)為正方體的一個(gè)面的中心，激振力為18 N，方向?yàn)椋?，0，1］.其中400 Hz與500 Hz正常狀態(tài)和故障（有裂紋）狀態(tài)的振動(dòng)模態(tài)見圖2～3.

由圖2可知，激振頻率在不斷增大的過程中，結(jié)構(gòu)體各個(gè)位置的振動(dòng)位移幅值變化越大，且離激振點(diǎn)的位置越近，位移幅值變化越明顯，激振點(diǎn)處的位移幅值從1.90E－06 m變化到3.00E－06 m；同時(shí)從正方體結(jié)構(gòu)可以得到，當(dāng)在一個(gè)面的中心位置進(jìn)行激振時(shí)，其它面上的中心位置振動(dòng)位移幅值要比其它各處變化明顯.

圖2 正常狀態(tài)振動(dòng)模態(tài)圖Fig.2 Vibration mode diagram for the normal state

由圖3可知，有裂紋的結(jié)構(gòu)體其振動(dòng)位移變化同樣滿足正相關(guān)的關(guān)系.對(duì)于有裂紋的結(jié)構(gòu)體，裂紋處的能量變化與正常結(jié)構(gòu)體相比較，變化很明顯.

為了研究聲輻射模態(tài)，對(duì)結(jié)構(gòu)體自身及周圍空氣的聲壓分布進(jìn)行了分析，圖4～5分別為結(jié)構(gòu)體400 Hz與500 Hz自身正常與故障狀態(tài)下的聲輻射模態(tài).

圖3 故障狀態(tài)振動(dòng)模態(tài)圖Fig.3 Vibration mode diagram for the fault state

圖4 正常狀態(tài)結(jié)構(gòu)體自身聲輻射模態(tài)Fig.4 The acoustic radiation modes of the structure for the normal state

圖5 故障狀態(tài)結(jié)構(gòu)體自身聲輻射模態(tài)Fig.5 The acoustic radiation modes of the structure for the fault state

圖4～5分析可知，隨著激振頻率的增大，結(jié)構(gòu)體自身的聲輻射值也隨著增大.完好正方體非激振面的聲壓由－4.52 Pa變化到－8.50 Pa；裂紋結(jié)構(gòu)體其裂紋處的聲壓值從－16.40 Pa也變化到－20.30 Pa.同時(shí)，裂紋對(duì)聲輻射值影響巨大.400 Hz時(shí)完好結(jié)構(gòu)體裂紋處的聲壓值為－4.52 Pa，而裂紋結(jié)構(gòu)體裂紋處為－16.40 Pa；500 Hz時(shí)完好結(jié)構(gòu)體裂紋處的聲壓值為－8.50 Pa，而裂紋結(jié)構(gòu)體裂紋處為－20.30 Pa.2種情況下變化均大于10 Pa.因此，可以利用這種變化關(guān)系來判斷結(jié)構(gòu)體是否發(fā)生了故障.同時(shí)，聲壓變化與振動(dòng)模態(tài)變化一致，呈正相關(guān)關(guān)系，其總體變化幅頻見圖6.

圖6 聲壓變化曲線圖Fig.6 The curve of pressure change

為進(jìn)一步對(duì)仿真結(jié)果進(jìn)行實(shí)驗(yàn)驗(yàn)證，本文對(duì)結(jié)構(gòu)體周圍空氣中的聲輻射模態(tài)進(jìn)行了仿真，其中，400 Hz與500 Hz的空氣聲輻射模態(tài)仿真結(jié)果見圖7～8（豎直面為激振面正對(duì)的面，水平面為有裂紋面正對(duì)的面）.

對(duì)圖7～8分析可知，隨著激振頻率的增大，結(jié)構(gòu)體周圍空氣中的聲壓值也在不斷增大；空氣中聲場(chǎng)的聲壓是圍繞著結(jié)構(gòu)體向外進(jìn)行輻射，且有裂紋的結(jié)構(gòu)體周圍空氣中的聲壓值與結(jié)構(gòu)體自身的聲壓值變化一致，這樣就可以利用傳聲器陣列采集結(jié)構(gòu)體周圍空氣中的聲壓信息來判斷結(jié)構(gòu)體的故障位置，為利用聲學(xué)信息進(jìn)行故障診斷提供了理論依據(jù).

2.3 正方體結(jié)構(gòu)振動(dòng)與聲輻射實(shí)驗(yàn)分析

圖8 故障狀態(tài)下結(jié)構(gòu)體周圍空氣聲輻射模態(tài)Fig.8 The acoustic radiation modes of the structure for the fault state

為驗(yàn)證上面的仿真結(jié)果，實(shí)驗(yàn)室利用半消音室對(duì)正方體結(jié)構(gòu)進(jìn)行了實(shí)驗(yàn)分析.實(shí)驗(yàn)利用傳聲器陣列架測(cè)量正常與有裂紋正方體在激振情況下的聲輻射聲壓，實(shí)驗(yàn)同樣是對(duì)與仿真同尺寸正常結(jié)構(gòu)與故障結(jié)構(gòu)進(jìn)行激勵(lì).實(shí)驗(yàn)裝置見圖9.裝置包括傳聲器陣列架、SDG5162型兩通道函數(shù)發(fā)生器、EPA-104型信號(hào)放大器、NIPXI-1033型采集箱及電腦.實(shí)驗(yàn)中將2個(gè)壓電陶瓷片用導(dǎo)電膠粘貼在模型上面，然后利用函數(shù)發(fā)生器發(fā)出信號(hào)，利用放大器對(duì)信號(hào)進(jìn)行放大后對(duì)模型進(jìn)行激振，同時(shí)利用壓電陶瓷片進(jìn)行拾振，將經(jīng)過激振的模型振動(dòng)信號(hào)及聲壓信號(hào)用采集箱及電腦進(jìn)行采集拾取，采樣頻率為10 kHz，這樣便可以得到實(shí)驗(yàn)的數(shù)據(jù).

圖10～11為400 Hz和500 Hz時(shí)利用傅里葉變化對(duì)信號(hào)進(jìn)行處理的正常結(jié)構(gòu)體和故障結(jié)構(gòu)體在不同頻率下的聲壓云圖.對(duì)圖10～11進(jìn)行對(duì)比分析可得：在正常情況下的400 Hz激振后圖中尖峰最高處獲得數(shù)據(jù)最大為0.009 7 Pa；在500 Hz激振頻率下其尖峰最高處的聲壓值為0.012 5 Pa，故障結(jié)構(gòu)體中間部位的聲壓幅值變化很大，圖11中尖峰最高處的數(shù)據(jù)分別為0.069 0 Pa和0.080 5Pa，由此可知隨著激振頻率的增大，結(jié)構(gòu)體周圍空氣中的聲壓值也在不斷增大，同時(shí)有故障的結(jié)構(gòu)體其聲壓值變化比正常情況下要大7倍左右，與仿真所得到的測(cè)量結(jié)果的趨勢(shì)是一致的.

圖9 實(shí)驗(yàn)裝置圖Fig.9 Experimental device

圖10 正常結(jié)構(gòu)體聲壓云圖Fig.10 The acoustic image of the structure for the normal state

圖11 故障結(jié)構(gòu)聲壓云圖Fig.11 The acoustic image of the structure for the fault state

3 結(jié) 論

本文對(duì)隨機(jī)激勵(lì)下正方體結(jié)構(gòu)振動(dòng)聲輻射問題進(jìn)行了研究.通過仿真與實(shí)驗(yàn)分析得到了以下結(jié)論：隨激振頻率的增大，本文結(jié)構(gòu)體自身的聲壓及振動(dòng)位移幅值在300 Hz到1 100 Hz內(nèi)不斷增大.同時(shí)有裂紋（故障）的結(jié)構(gòu)體故障位置處的聲壓及振動(dòng)幅值與正常結(jié)構(gòu)體相比的變化十分明顯；而且裂紋對(duì)應(yīng)的空氣中聲場(chǎng)分布的聲壓變化也比較大.因此，可以依據(jù)聲壓云圖的異常變化對(duì)塊狀結(jié)構(gòu)的缺陷進(jìn)行初步診斷，也可以為其它利用聲學(xué)信息進(jìn)行故障診斷的研究提供參考.

［1］ 王治國(guó).工程聲學(xué)有限元分析理論與應(yīng)用［M］.北京：國(guó)防工業(yè)出版社，2007.

WANG Z G.Engineering acoustics finite element analysis theory and its application［M］.Beijing：National Defense Indus-try Press，2007.

［2］ 曾旭，左曙光.基于有限元和邊界元的輪胎振動(dòng)聲輻射仿真計(jì)算［J］.佳木斯大學(xué)學(xué)報(bào)：自然科學(xué)版，2012，30（4）：485-489.

ZENG X，ZUO SG.Simulation calculation of tire vibration sound radiation based on the finite elementmethod and boundary elementmethod［J］.JJiamusi Univ：Nat Sci Edi，2012，30（4）：485-489.

［3］ 李凱.基于聲強(qiáng)可視化的船舶結(jié)構(gòu)聲振能量特性研究［D］.大連：大連理工大學(xué)，2011.

LIK.Study on the acoustic energy of ship structure based on sound intensity visualizationmethod［D］.Dalian：Dalian University of Technology，2011.

［4］ 劉紅光，陸森林.車廂壁面振動(dòng)對(duì)其內(nèi)部聲場(chǎng)的影響度分析與阻尼降噪［J］.農(nóng)業(yè)工程學(xué)報(bào)，2002，18（2）：62-64.

LIU H G，LU SL.The influence of the vibration of the compartmentwall on the internal acoustic field and the damping and noise reduction［J］.Transact Chin Soc Agr Engin，2002，18（2）：62-64.

［5］ 謝昌林.機(jī)械結(jié)構(gòu)振動(dòng)噪聲的仿真［D］.西安：西北工業(yè)大學(xué)，2004.

XIE C L.Simulation of mechanical structure vibration and noise［D］.Xi'an：North Western Poly Technical University，2004.

［6］ 劉寶山，趙國(guó)忠，李瑞勇.基于虛擬激勵(lì)法的結(jié)構(gòu)隨機(jī)振動(dòng)聲輻射分析［J］.固體力學(xué)學(xué)報(bào)，2011，32（6）：581-586.

LIU B S，ZHAO G Z，LIR Y.Analysis of the acoustic radiation of random vibration based on the virtual excitationmethod［J］.Chin JSolid Mechan，2011，32（6）：581-586.

［7］ 周健.板結(jié)構(gòu)的聲輻射與隔聲性能研究［D］.武漢：華中科技大學(xué)，2002.

ZHOU J.The acoustic radiation and sound insulation of the structure of slab［D］.Wuhan：Huazhong University of Science and Technology，2002.

［8］ ELLIOTT S J，JOHNSON M E.Radiation modes and the active control of sound power［J］.JAcoust Soc Amer，1993，94（4）：2194-2204.

［9］ TAO J，GEH，QIU X.A new rule of vibration sampling for predicting acoustical radiation from rectangular plates［J］.Appl Acoust，2006，67（8）：756-770.

［10］喬思茂，包善斐，于駿一.單點(diǎn)激振試驗(yàn)中激振點(diǎn)位置的選擇［J］.吉林工業(yè)大學(xué)學(xué)報(bào)，1992，65（1）：102-104.

QIAO X M，BAO SF，YU JY.The choice of exciting point in single point excited vibration test［J］.J Jilin Indust Univ，1992，65（1）：102-104.

Fault detection of structural vibro-acoustical radiation based on ACTRAN

ZHANG Chun-liang1，YAO M ing-jing2，YUEXia1，ZHU Hou-yao1，Q IGuan-tan1

（1.School of Mechanical and Electrical Engineering，Guangzhou University，Guangzhou 510006，China；
2.The Engineering＆Technical College，Chengdu University of Technology，Leshan 614000，China）

This paper first introduces the dynamic equation of fluid-solid coupling，followed by numerical simulation using ACTRAN software on the vibration of the cable structure acoustic radiation problem.Finite element method is used to define differentmonitoring loading infinity and analyze the relationship between the sound radiation and vibration amplitude.Results show that vibration and sound radiation has a proportional relationship，namely when the vibration is severe，the sound pressure of the correspondingmechanical parts is high；when the vibration is small，the corresponding sound pressure is low.Furthermore，a cracked structure presents high vibration amplitude and the sound pressure value is higher than the healthy structure.This study provides the theoretical foundation for the intelligent fault diagnosis ofmechanical equipment research using acoustical signals.

sound-solid coupling；vibro-acoustic；acoustic radiation

TH 172

A

【責(zé)任編輯：陳 鋼】

1671-4229（2015）03-0053-06

2015－04－06；

2015－04－21

國(guó)家自然科學(xué)基金資助項(xiàng)目（51275099，51305086）；廣東省自然科學(xué)基金資助項(xiàng)目（S2012010009505）；廣州市羊城學(xué)者首席科學(xué)家基金資助項(xiàng)目（12A006S）

張春良（1964－），男，教授，博導(dǎo)，博士.E-mail：nhzcl＠163.com