郝 勇 王經(jīng)磊 劉 鉑 孫浩誠 韓蕓徽 陳 豐

(河北建筑工程學(xué)院，河北 張家口 075000)

0 前 言

隨著我國經(jīng)濟的快速發(fā)展，一些高精密電子廠房，如芯片廠、面板廠等，以及一些大型橋梁在不斷地建設(shè).這些建筑不僅要滿足承載力要求，還要防止外界振動對結(jié)構(gòu)產(chǎn)生不利影響.為獲取結(jié)構(gòu)動力特性的相關(guān)參數(shù)，防止外界振動對結(jié)構(gòu)安全及設(shè)備正常運行產(chǎn)生不利影響，需要對結(jié)構(gòu)進行模態(tài)測試及相關(guān)數(shù)據(jù)分析.

模態(tài)分析是以識別振動系統(tǒng)模態(tài)參數(shù)為目的的一種分析方法.它把振動系統(tǒng)中相互耦合的運動方程通過模態(tài)坐標變換，使方程組解耦，進而求出各階模態(tài)參數(shù).由于模態(tài)分析技術(shù)在振動控制方面存在巨大的應(yīng)用價值，自上世紀中期開始，模態(tài)分析技術(shù)迅速發(fā)展，如今在線性模態(tài)分析方面已發(fā)展至成熟階段.目前國內(nèi)外已研制出多款模態(tài)分析軟件，這些軟件為解決振動問題發(fā)揮了巨大作用.

1 理論基礎(chǔ)

為引入模態(tài)參數(shù)，以n個自由度的無阻尼系統(tǒng)實模態(tài)分析為研究對象，闡述相關(guān)參數(shù)的物理意義[1].

(1)

1.1 自由振動

(2)

設(shè)特解x=φejωt

其中φ為自由響應(yīng)的位移幅值列向量，將特解代入(1-2)式得(K-ω2M)φ=0

(3)

式(1-3)為n階齊次線性方程組，式中K-ω2M為系數(shù)矩陣，ω2為方程組特征值，φ為方程組特征向量，要使該方程組有非零解，則系數(shù)矩陣的行列式|K-ω2M|=0,設(shè)該方程無重根，可解得此方程有n個不等正根ω0i(i=1,2,…,n)，ω0i(i=1,2,…,n)即為結(jié)構(gòu)各階模態(tài)的自振頻率，將ω0i(i=1,2,…,n)代入式(1-3)，可得n個線性無關(guān)的特征向量φi(i=1,2,…,n)，φi為第i階頻率對應(yīng)的振型.將φi(i=1,2,…,n)組合成一個n×n的方陣，即φ=[φ1,φ2,…,φn]，φ被稱為模態(tài)矩陣.

把其中一個特征值ω0i與其對應(yīng)的特征向量φi代入式(1-3)，并左乘另外一個特征向量的轉(zhuǎn)置φjT得φjT(K-ω0i2M)φi=0

(4)

然后再將第j個特征值ω0j和其對應(yīng)的特征向量φj代入式(1-3)得

(K-ω0j2M)φj=0

(5)

將式(1-5)進行轉(zhuǎn)置(K、M為對稱矩陣)，右乘φi得

φjT(K-ω0j2M)φi=0

(6)

(1-4)與(1-6)相減得(w0j2-ω0i2)φjTMφi=0

(7)

當i≠j時，ω0i2≠ω0j2，則φjTMφi=0.

當i=j時，令mi=φiTMφi，定義mi為模態(tài)質(zhì)量，因為M為正定矩陣，所以mi≥0.

當i≠j時，因φjTMφi=0，故φjTKφi=0.

當i=j時，令ki=φiTKφi，定義ki為模態(tài)剛度，因為K為正定或半正定矩陣，所以ki≥0.

通過以上理論推導(dǎo)可知，各階模態(tài)的運動是相互獨立的，彼此之間不發(fā)生能量交換，且每階模態(tài)的能量守恒.

根據(jù)模態(tài)展開定理，可得物理坐標系的自由響應(yīng)函數(shù)

(8)

式中yi(i=1,2,…,n)為向量x對應(yīng)的模態(tài)坐標.

令初始條件y0=φ-1x0，的模態(tài)坐標自由響應(yīng)函數(shù)yi=yisin(ω0it+θi)

(9)

將式(1-9)代入式(1-8)可得物理坐標系的自由響應(yīng)函數(shù)

1.2 頻響函數(shù)

無阻尼振動受簡諧激勵，則系統(tǒng)的運動為與激勵同頻的簡諧運動，令f(t)=Fejωt穩(wěn)態(tài)位移響應(yīng)為x(t)=Xejωt,將x(t)與f(t)代入式(1-1)可得(K-ω2M)X=F，令X=H(ω)F，所以H(ω)=(K-ω2M)-1，H(ω)被稱為頻響函數(shù)或傳遞函數(shù).

1.2.1 頻響函數(shù)的模態(tài)展開式

(10)

設(shè)穩(wěn)態(tài)位移響應(yīng)函數(shù)y=Uejet，代入式(1-10)，整理得

(11)

頻響函數(shù)是模態(tài)分析中頻域法參數(shù)識別的基礎(chǔ).

2 模態(tài)參數(shù)測試流程

模態(tài)測試系統(tǒng)包含三個部分：①激振系統(tǒng)；②傳感系統(tǒng)；③分析系統(tǒng).其中激振系統(tǒng)主要包括激振器、功率放大器和信號源.傳感系統(tǒng)主要包括傳感器、連接線等.分析系統(tǒng)主要包括數(shù)據(jù)接收裝置、機箱及相關(guān)軟件等.分析儀接受傳感器產(chǎn)生的信號，并進行數(shù)字分析與計算.測試相關(guān)流程見圖1.

圖1 模態(tài)測試流程

以錘擊法為例，利用LMS測試系統(tǒng)[2]對模態(tài)測試流程進行詳細敘述.

在測試之前首先應(yīng)設(shè)置力錘通道和加速度傳感器通道，然后進行錘擊示波，定義示波的帶寬及頻率分辨率.用力錘(如圖2)敲擊被測結(jié)構(gòu)，錘擊力度應(yīng)盡量保持一致，通過錘擊示波使系統(tǒng)確定一個合適的量程范圍.

接下來進行錘擊測試的設(shè)置，包括觸發(fā)級設(shè)置，帶寬、窗函數(shù)及錘擊點的選取.

用力錘錘擊被測結(jié)構(gòu)，直到在電腦屏幕上看到脈沖，選取應(yīng)用建議值，軟件可自動進行觸發(fā)級設(shè)置.然后繼續(xù)錘擊被測結(jié)構(gòu)，觀察頻率帶寬范圍內(nèi)激勵點譜的分布來檢查所定義帶寬是否合適.錘頭的選擇應(yīng)根據(jù)頻率帶寬來確定.對窗函數(shù)的選擇仍然需要用力錘錘擊被測結(jié)構(gòu)，通過錘擊，軟件可以自動選擇窗函數(shù)，也可以手動設(shè)置.激勵點的選取應(yīng)對被測結(jié)構(gòu)的不同位置進行錘擊，通過驅(qū)動點設(shè)置菜單可以將所測驅(qū)動點進行顯示和比較，選取能夠有效激起各階模態(tài)的激勵點作為最佳激勵點.以上設(shè)置完成后，可進行正式的模態(tài)測試.

圖2 力錘 圖3 三軸加速度傳感器

圖4 數(shù)字式頻率分析儀

模態(tài)測試中，要在結(jié)構(gòu)的X向、Y向和Z向分別選取一個最佳激勵點.將三軸加速度傳感器(如圖3)依次放置在測點上.用力錘依次錘擊三個最佳激勵點，即可記錄每個測點上的激勵響應(yīng).力錘和傳感器將激勵和響應(yīng)的時域信號傳輸?shù)椒治鰞x(如圖4)，通過分析儀和電腦軟件的相關(guān)處理與計算可獲得結(jié)構(gòu)的非模態(tài)參數(shù)，如頻響函數(shù)或脈沖響應(yīng)函數(shù).對頻響函數(shù)或脈沖響應(yīng)函數(shù)進一步分析得到結(jié)構(gòu)振動特性.

3 模態(tài)分析的應(yīng)用

隨著振動理論的發(fā)展，基于承載力和位移的設(shè)計指標已不再滿足結(jié)構(gòu)工程的設(shè)計要求，尤其是大型復(fù)雜結(jié)構(gòu)和對振動較敏感的結(jié)構(gòu).為了解結(jié)構(gòu)的振動特性，模態(tài)分析技術(shù)在解決工程振動問題中得到了廣泛的應(yīng)用，如對結(jié)構(gòu)模型進行模態(tài)分析，推算原型的動力特性，或通過模態(tài)分析診斷結(jié)構(gòu)缺陷等.

3.1 模態(tài)分析在汽輪機基礎(chǔ)中的應(yīng)用

汽輪機的正常運行會對其基礎(chǔ)產(chǎn)生持續(xù)振動，通過LMS Test Lab系統(tǒng)對汽輪機基礎(chǔ)的自振特性進行測試，可獲得結(jié)構(gòu)的各階自振頻率、振型及阻尼比等相關(guān)參數(shù).

汽輪機是一個連續(xù)的結(jié)構(gòu)，故其階數(shù)有無窮多階，在實際工程中，一般只取前50至前80階自振頻率及對應(yīng)的振型進行研究.通過對相關(guān)文獻的對比研究(見表1)發(fā)現(xiàn)，汽輪機的工作頻率一般在50Hz左右，為防止結(jié)構(gòu)發(fā)生共振，基礎(chǔ)各階振型所對應(yīng)的自振頻率在50Hz左右的范圍內(nèi)分布不應(yīng)太密，越稀疏越有利于汽輪機的安全運行.在低階振型范圍內(nèi)，汽輪機基礎(chǔ)主要以整體平動為主，在高階振型的范圍內(nèi)，汽輪機基礎(chǔ)主要以局部彎曲、扭轉(zhuǎn)為主.基礎(chǔ)整體平動對汽輪機的安全運行影響較小，而基礎(chǔ)局部的扭轉(zhuǎn)、彎曲對汽輪機的安全運行影響較大.

表1 汽輪機基礎(chǔ)自振特性的對比研究

除對汽輪機基礎(chǔ)進行模態(tài)分析，獲得各階自振頻率、阻尼比及振型外，還應(yīng)對擾力點的振動線位移進行計算.擾力點的位置一般為汽輪機軸承所在位置，擾力值一般由汽輪機生產(chǎn)廠家提供.在LMS Virtual.Lab軟件中應(yīng)建立汽輪機基礎(chǔ)模型，在擾力點位置施加X、Y、Z三個方向的擾力.通過軟件計算可知各測點的振動線位移響應(yīng)，檢測各測點的振動線位移是否符合《動規(guī)》要求.對于振動線位移響應(yīng)較大的測點，應(yīng)對相關(guān)測點加強振動控制，如增加減震裝置或加強局部剛度等.

3.2 模態(tài)分析在結(jié)構(gòu)損傷檢測和損傷預(yù)警中的應(yīng)用

利用模態(tài)分析對結(jié)構(gòu)損傷檢測和損傷預(yù)警有多種方法，對相關(guān)文獻進行研究(見表2)，常見方法有剛度法、模態(tài)置信因子法、模態(tài)柔度法等.其原理都是由結(jié)構(gòu)損傷引起的結(jié)構(gòu)質(zhì)量和剛度發(fā)生改變，這兩種因素的改變會導(dǎo)致結(jié)構(gòu)的模態(tài)發(fā)生變化.對比損傷前與損傷后某一模態(tài)參數(shù)的變化或?qū)p傷后某一模態(tài)參數(shù)進行研究，可確定結(jié)構(gòu)的損傷位置和損傷程度.

表2 模態(tài)分析在結(jié)構(gòu)損傷檢測及損傷預(yù)警中的研究

4 結(jié) 論

通過對模態(tài)分析基本理論的研究，得出各階模態(tài)正交的理論依據(jù)及頻響函數(shù)的推導(dǎo)過程，為模態(tài)分析中參數(shù)識別法的研究奠定基礎(chǔ).

對于現(xiàn)場模態(tài)測試，由于受多種環(huán)境因素的影響，一般多采用錘擊法測試模態(tài).把所測數(shù)據(jù)進行模態(tài)分析，其分析結(jié)果可為結(jié)構(gòu)設(shè)計、振動控制和結(jié)構(gòu)損傷及損傷預(yù)警等工程問題提供幫助.