基于實測應(yīng)變的時域模態(tài)參與因子識別研究

胡偉鋼，劉志明，王 曦，林浩博

(北京交通大學(xué) 機(jī)械與電子控制工程學(xué)院，北京 100044)

結(jié)構(gòu)疲勞強(qiáng)度決定了機(jī)械系統(tǒng)的正常運轉(zhuǎn)及使用壽命，且關(guān)系到飛機(jī)、鐵路列車、汽車等交通工具的運行安全，在結(jié)構(gòu)強(qiáng)度設(shè)計和運營評估中，了解結(jié)構(gòu)的整體應(yīng)變分布場和關(guān)鍵位置結(jié)構(gòu)響應(yīng)是十分必要的[1]。

結(jié)構(gòu)響應(yīng)一般分為非結(jié)構(gòu)振動響應(yīng)和結(jié)構(gòu)振動響應(yīng)[2]。鐵道車輛部件對非結(jié)構(gòu)振動響應(yīng)測試一般采用直接測試法和準(zhǔn)靜態(tài)標(biāo)定法[3-4]，對結(jié)構(gòu)振動響應(yīng)測試一般采用直接測試法[5-6]。目前，對于復(fù)雜結(jié)構(gòu)振動響應(yīng)測量一般是定性給出傳感器布置方案[7-8]，在布置傳感器過程中會有大量的經(jīng)驗性和不確定性。

本文提出一種通過優(yōu)化結(jié)構(gòu)表面可測部位測點的方法，逆求結(jié)構(gòu)動態(tài)響應(yīng)模態(tài)參與因子，獲得結(jié)構(gòu)整體應(yīng)變響應(yīng)。該方法利用應(yīng)變片質(zhì)量輕,在結(jié)構(gòu)測試時不會改變被測結(jié)構(gòu)動態(tài)參數(shù)的特點，通過在結(jié)構(gòu)表面上粘貼合適數(shù)量的應(yīng)變片，基于D優(yōu)化理論優(yōu)化最佳測點位置和方向，根據(jù)模態(tài)疊加基本原理，得到整體結(jié)構(gòu)響應(yīng)特征。該方法為復(fù)雜結(jié)構(gòu)無法實測部位結(jié)構(gòu)響應(yīng)提供了一種間接測試方法。

1 模態(tài)參與因子估計基本理論

1.1 結(jié)構(gòu)振動

對于一個結(jié)構(gòu)可離散為N自由度的系統(tǒng)，其平衡運動方程為[9]

( 1 )

其中n×1的x(t)為結(jié)構(gòu)振動響應(yīng)。結(jié)構(gòu)振動系統(tǒng)的振動響應(yīng)可以表示為位移模態(tài)的線性組合。

( 2 )

式中：Φ=[φ1φ2…φn]，為n×n結(jié)構(gòu)模態(tài)矩陣；qi(t)是模態(tài)振型φi的模態(tài)參與因子；q(t)=[q1q2…qn]為結(jié)構(gòu)振型Φ模態(tài)參與因子向量。

在模態(tài)坐標(biāo)下，結(jié)構(gòu)運動方程可以表示為[10]

( 3 )

對于一個連續(xù)的結(jié)構(gòu)體，在時域上結(jié)構(gòu)變形與應(yīng)變存在如下的線性微分關(guān)系

ε(t)=Dx(t)

( 4 )

式中：D是微分算子。通過微分運算對式( 2 )進(jìn)行計算，對特定時間點模態(tài)的線性疊加為

( 5 )

結(jié)構(gòu)模態(tài)的微分即為結(jié)構(gòu)應(yīng)變模態(tài)，根據(jù)式( 4 )、式( 5 )可以得到

( 6 )

式中：Ψε=[φ1φ2…φn]為結(jié)構(gòu)模態(tài)應(yīng)變矩陣。式( 6 )表明對于在線彈性范圍內(nèi)的小變形結(jié)構(gòu)上任意點的時間域應(yīng)變響應(yīng)可以表示為模態(tài)應(yīng)變的線性疊加，即模態(tài)疊加。

1.2 模態(tài)參與估計

( 7 )

( 8 )

( 9 )

2 應(yīng)變片測點布置方法

2.1 測試候選區(qū)

定義在結(jié)構(gòu)表面上適宜安裝應(yīng)變片的所有區(qū)域為測試候選區(qū)，在測試候選區(qū)內(nèi)可粘貼應(yīng)變片測點為測試候選點，測試候選點與位置和方向有關(guān)，通過優(yōu)化得到參與模態(tài)參與因子估計的測點組合為優(yōu)化測點組。結(jié)合文獻(xiàn)[3，8]所述的傳感器布置方案及結(jié)構(gòu)的幾何特征，給出適宜布置應(yīng)變片測試候選區(qū)的選擇原則：

(1)測試候選區(qū)應(yīng)充分反映結(jié)構(gòu)彈性模態(tài)變形的整體特征，特別是結(jié)構(gòu)的主要模態(tài)振動方向的變形。

(2)測試候選區(qū)必須保證本區(qū)域結(jié)構(gòu)為彈性變形且應(yīng)變梯度變化不應(yīng)過大，這些區(qū)域包括載荷激勵區(qū)、結(jié)構(gòu)幾何應(yīng)力集中區(qū)、結(jié)構(gòu)缺陷區(qū)等。

(3)測試候選區(qū)應(yīng)避開應(yīng)變片安裝空間位置無法滿足區(qū)域和傳感器信號傳輸困難的區(qū)域，這些位置包括結(jié)構(gòu)裝配區(qū)域，如箱體內(nèi)部、結(jié)構(gòu)設(shè)計孔等。

在建立由實體單元離散的有限元模型時，利用薄膜單元在結(jié)構(gòu)振動中只傳遞面內(nèi)力，沒有彎曲剛度的特點，在結(jié)構(gòu)表面測試候選區(qū)域增加一層薄膜單元，通過薄膜單元應(yīng)變信息來表達(dá)測試候選區(qū)的應(yīng)變片信息。由于節(jié)點應(yīng)變信息存在不唯一性，因此，選擇單元應(yīng)變信息進(jìn)行分析計算。

測點方向θ是測點敏感度的變量，單元在不同角度θ(0°≤θ<180°)得到各自的應(yīng)變張量，應(yīng)變張量由xyz坐標(biāo)系轉(zhuǎn)化到任意x′y′z′坐標(biāo)系[12]

式中：Tm是x′y′z′坐標(biāo)系對xyz坐標(biāo)系的轉(zhuǎn)換矩陣。薄膜單元的局部坐標(biāo)系內(nèi)，z軸方向始終垂直于單元面，其轉(zhuǎn)換矩陣為

2.2 D優(yōu)化設(shè)計

由優(yōu)化矩陣[ψε]opt可得優(yōu)化模態(tài)參與因子估計為

3 懸臂梁模態(tài)參與因子識別

3.1 激振懸臂梁有限元模型

將懸臂梁作為工程對象進(jìn)行模態(tài)參與因子識別，搭建懸臂梁試驗激振臺，如圖1所示，懸臂梁方形端后半部分固定，另一端開孔為激振桿裝配孔，激振桿與懸臂梁之間為螺栓連接，激振器通過激振桿傳遞激振力。

圖1 試驗懸臂梁

按照2.1節(jié)所述方法選擇測試候選區(qū)，如圖1中紅色區(qū)域。對懸臂梁模型進(jìn)行模態(tài)分析，表1為懸臂梁前3階模態(tài)頻率。

表1 試驗懸臂梁前3階模態(tài)頻率

選擇6個測點估計懸臂梁前3階模態(tài)參與因子，按照2.2所述運用D優(yōu)化法得到測試候選區(qū)內(nèi)測點最佳位置和方向，圖1中給出由測點G#1～G#6組成的優(yōu)化測點組，預(yù)留候選區(qū)內(nèi)部測點G#7和外部測點G#8驗證模態(tài)參與因子估計。

圖2所示為懸臂梁優(yōu)化測點組G#1～G#6和驗證測點G#7與G#8的前3階模態(tài)應(yīng)變矩陣。

圖2 懸臂梁各測點前3階模態(tài)應(yīng)變矩陣

3.2 模態(tài)參與因子識別及驗證

對激振器提供25～300 Hz的等幅值掃頻信號，同步采集各應(yīng)變片信號，采樣頻率為2 000 Hz，得到各測點的應(yīng)變響應(yīng)，根據(jù)式(13)，采用最小二乘法估計各階模態(tài)參與因子，如圖3所示。

圖3 前3階模態(tài)參與因子估計

根據(jù)上文的模態(tài)參與因子和模態(tài)應(yīng)變矩陣，由式( 6 )模態(tài)疊加原理可得各測點應(yīng)變響應(yīng)時域波形。圖4(a)和圖5(a)所示分別為驗證測點G#7和G#8的試驗值和反求值時域波形。為了清晰表達(dá)時域波形的重合度，將試驗信號和反求信號建立散點圖，圖4(b)和圖5(b)所示分別為驗證測點G#7和G#8的試驗值和反求值的散點圖。

圖4 測點G#7應(yīng)力的試驗值和反求值

圖5 測點G#8應(yīng)力的試驗值和反求值

表2所示為各測點試驗值和反求值的散點擬合相關(guān)性及RMS誤差，各測點線性擬合斜率接近1，相關(guān)性接近1，均方根誤差均較小，說明試驗應(yīng)力信號和反求應(yīng)力信號波形重合度高。驗證測點G#7和G#8的散點圖擬合斜率接近1，擬合相關(guān)性高，均方根誤差小，反求波形和實際測試波形重合度高，說明運用模態(tài)參與因子估計的方法估計準(zhǔn)確度高。

表2 各測點擬合相關(guān)性及RMS誤差

4 箱體模態(tài)參與因子識別

4.1 箱體有限元模型

齒輪箱作為高速列車傳動系統(tǒng)重要組成部分，對于保障高鐵列車的動力傳輸具有重要作用，其主要由大小齒輪、箱體和聯(lián)軸節(jié)構(gòu)成，其中，小齒輪通過聯(lián)軸節(jié)與牽引電機(jī)相連，大齒輪直接壓裝在輪軸上，齒輪箱箱體的一端通過抱軸承懸掛在車軸上，另一端通過吊桿吊掛在構(gòu)架橫懸臂梁上[15]。圖6所示為齒輪箱箱體有限元模型，在箱體安裝軸承的中心建立參考點RP1、RP2、RP3和RP4，吊餅區(qū)中心建立參考點RP5，箱體和軸承的接觸面與對應(yīng)參考點之間建立耦合關(guān)系，參考點RP5與吊餅面建立耦合關(guān)系，約束參考點RP5的Y和Y轉(zhuǎn)動方向，箱體通過抱軸方式懸掛在車軸上，因此，約束參考點RP1和RP2的X、Y、Z、X轉(zhuǎn)動和Y轉(zhuǎn)動方向。

圖6 齒輪箱箱體模型

對箱體進(jìn)行模態(tài)分析，圖7所示為第1階模態(tài)應(yīng)變響應(yīng)，表3為前11階模態(tài)頻率。

圖7 第1階約束模態(tài)響應(yīng)

模態(tài)階次模態(tài)頻率/Hz1659.721023.631265.241389.151826.261900.772220.382594.092631.1102813.3112890.1

按照2.1節(jié)所述方法選擇適宜的測試候選區(qū)，如圖8所示，去掉所有不適合粘貼應(yīng)變片位置(包括箱蓋、吊桿、銘牌、軸承蓋等附近)。選擇22個測點估計箱體的前11階模態(tài)參與因子，按照2.2節(jié)所述D優(yōu)化得到測試候選區(qū)內(nèi)測點最佳位置和方向，圖8給出了由優(yōu)化得到測點G&1～G&22的優(yōu)化測點組。

圖8 箱體測試候選區(qū)及優(yōu)化測點組測點

圖9所示為箱體優(yōu)化測點組測點G&1～G&22前11階模態(tài)應(yīng)變矩陣。

圖9 箱體各測點前11階模態(tài)應(yīng)變矩陣

4.2 箱體模態(tài)參與因子識別

高速列車齒輪傳動系統(tǒng)通常采用斜齒輪單級傳動，輪齒嚙合產(chǎn)生軸向力、徑向力和垂向力，這三個方向的力由軸承傳遞給箱體。根據(jù)文獻(xiàn)[15]齒輪箱箱體軸承載荷分布關(guān)系和文獻(xiàn)[16]中運營工況下齒輪箱箱體載荷值，參考點RP3的縱向載荷幅值、垂向載荷幅值分別2 194、4 227 N，參考點RP4的縱向載荷幅值、垂向載荷幅值、橫向載荷幅值分別2 194、4 227、3 077 N。對RP3和RP4按照上述對應(yīng)載荷幅值加載動態(tài)掃頻載荷，掃頻范圍為500～1 500 Hz，掃頻時間為1 s，運用瞬態(tài)動力學(xué)仿真模擬，提取測點G&1～G&22應(yīng)變響應(yīng)。

實際測試和模擬仿真一般存在誤差，如式(14)所示，對應(yīng)變測點G&1～G&22理論值加15%的白噪聲幅值誤差進(jìn)行模態(tài)參與因子估計，如圖10所示。

εTest(t)=εth(t)+Δε(t)

(14)

式中：εTest(t)為實際測量應(yīng)變響應(yīng)；εth(t)為理論應(yīng)變響應(yīng)；Δε(t)為誤差應(yīng)變。

圖10 測點G&1理論應(yīng)變響應(yīng)波形和15%白噪聲幅值誤差波形

根據(jù)式(13)，運用最小二乘法估計各階模態(tài)參與因子，如圖11所示(GU1為仿真波形，MPF1為識別波形)。表4給出15%應(yīng)變幅值誤差的各階模態(tài)參與因子理論值和反求值擬合相關(guān)性及RMS誤差。

圖11 第1階模態(tài)參與因子波形及散點圖

模態(tài)階次模態(tài)參與量/%斜率k擬合優(yōu)度r2RMS誤差/%180.1600.99980.99961.03211.8000.99990.99593.6646.9500.99960.99812.1831.0101.00040.953711.94100.0700.99860.96809.2850.0041.00120.923622.0760.0031.00100.942418.9570.0021.00000.904224.75801.00420.3822130.2590.0011.00090.97008.801101.01520.0920171.50

由表4可知，第1、2、4階模態(tài)參與量之和占總模態(tài)的98.91%，模態(tài)參與因子的識別波形與理論波形之間RMS誤差較小，第3、10階模態(tài)參與量之和占總模態(tài)1.08%，模態(tài)參與因子的識別信號與理論信號的RMS誤差不超過12%；第5、6、7、8、9和11階模態(tài)的模態(tài)參與量之和占總模態(tài)的0.01%，對結(jié)構(gòu)振動貢獻(xiàn)量較少，即使RMS誤差較大，對結(jié)構(gòu)應(yīng)變響應(yīng)的能量貢獻(xiàn)較小，可以忽略。

箱體的第1、2、4階模態(tài)可以表達(dá)結(jié)構(gòu)的整體振動特征，通過本文提出的方法可以準(zhǔn)確識別結(jié)構(gòu)振動的主要模態(tài)參與因子，說明運用有限測點測試可以識別箱體整體結(jié)構(gòu)特征。

5 結(jié)論

本文提出了一種基于優(yōu)化結(jié)構(gòu)表面應(yīng)變測點實現(xiàn)結(jié)構(gòu)振動模態(tài)參與因子的時域識別方法，結(jié)合模態(tài)應(yīng)變響應(yīng)矩陣，可以獲得結(jié)構(gòu)整體響應(yīng)。文中將懸臂梁板作為研究對象，選擇6個應(yīng)變測點對25～300 Hz的掃頻信號進(jìn)行前3階模態(tài)參與因子時域波形估計，由模態(tài)疊加反求的兩驗證測點應(yīng)力的RMS誤差值較小，說明該方法模態(tài)參與因子估計準(zhǔn)確度高。

對高速列車齒輪箱箱體進(jìn)行小齒輪軸承座位置500～1 500 Hz掃頻仿真，選擇22個測點的仿真應(yīng)變響應(yīng)并對每個測點加15%的白噪聲幅值誤差，根據(jù)模態(tài)疊加基本原理，結(jié)合模態(tài)應(yīng)變響應(yīng)，運用最小二乘法估計前11階模態(tài)參與因子時域波形，其中，模態(tài)參與量總和為98.9%的第1、2、4階模態(tài)參與因子估計的RMS誤差最大為3.66%，而其他階模態(tài)參與因子雖然誤差較大，但模態(tài)參與量較小，對結(jié)構(gòu)振動響應(yīng)誤差較小，箱體的模態(tài)參與因子識別效果好。

參考文獻(xiàn)：

[1]王萌,李強(qiáng),孫守光.基于應(yīng)力響應(yīng)的多頻率激勵載荷識別研究[J].鐵道學(xué)報,2015,37(2):27-33.

WANG Meng,LI Qiang,SUN Shouguang.Study of Multi-frequency Exciting Load Identification Based on Dynamic Stress Response[J].Journal of the China Railway Society,2015,37(2):27-33.

[2]姚起杭,姚軍.工程結(jié)構(gòu)的振動疲勞問題[J].應(yīng)用力學(xué)學(xué)報,2006,23(1):12-15.

YAO Qihang,YAO Jun.Vibration Fatigue in Engineering Structure[J].Chinese Journal of Applied Mechanics,2006,23(1):12-15.

[3]ZHU N,SUN S,LI Q,et al.Theoretical Research and Experimental Validation of Elastic Dynamic Load Spectra on Bogie Frame of High-speed Train[J].Chinese Journal of Mechanical Engineering,2016,29(3):498-506.

[4]鄒驊,李強(qiáng),孫守光.基于載荷標(biāo)定的城際列車轉(zhuǎn)向架載荷及應(yīng)力分布特征研究[J].鐵道學(xué)報,2016,38(10):27-33.

ZOU Hua,LI Qiang,SUN Shouguang.Study on Intercity Train Load Spectrum Distribution Estimation and Calibration Methods Based on Load Demarcation[J].Journal of the China Railway Society,2016,38(10):27-33.

[5]FU D,WANG W,DONG L.Analysis on the Fatigue Cracks in the Bogie Frame[J].Engineering Failure Analysis,2015,58:307-319.

[6]王文靜,惠曉龍,馬紀(jì)軍.高速列車設(shè)備艙支架疲勞裂紋機(jī)理研究[J].機(jī)械工程學(xué)報,2015,51(6):142-147.

WANG Wenjing,HUI Xiaolong,MA Jijun.Fatigue Crack Mechanism Research on High Speed Train Equipment Cabin Frame[J].Journal of Mechanical Engineering,2015,51(6):142-147.

[7]廖芳,高衛(wèi)民,王承,等.基于模態(tài)擴(kuò)展的變速器箱體振動識別及輻射噪聲優(yōu)化[J].同濟(jì)大學(xué)學(xué)報:自然科學(xué)版,2012,40(11):1698-1703.

LIAO Fang,GAO Weimin,WAN Cheng,et al.Vibration Identification of Gearbox Housing by Modal Expansion and Radiation Sound Optimization[J].Journal of Tongji University:Natural Science,2012,40(11):1698-1703.

[8]張方,秦遠(yuǎn)田.工程結(jié)構(gòu)動載荷識別方法[M].北京:國防工業(yè)出版社,2011.

[9]克雷格.結(jié)構(gòu)動力學(xué)[M].常玲,李振邦,譯.北京:人民交通出版社,1996.

[10]RAO S S.Mechanical Vibrations[M].5th ed.New Jersey:Prentice Hall,2010.

[11]CHATTERJEE S,HADI A S.Sensitivity Analysis in Linear Regression[M].New Jersey:John Wiley,1988.

[12]MITCHELL T J.An Algorithm for the Construction of D-optimal Experimental Designs[J].Technometrics,2000,16(2):203-210.

[13]GUPTA D K,DHINGRA A K.Input Load Identification from Optimally Placed Strain Gages Using D-optimal Design and Model Reduction[J].Mechanical Systems & Signal Processing,2013,40(2):556-570.

[14]GALIL Z,KIEFER J.Time and Space-saving Computer Methods,Related to Mitchell’s DETMAX,for Finding D-Optimum Designs[J].Technometrics,1980,22(3):301-313.

[15]黃冠華.高速列車齒輪傳動系統(tǒng)動態(tài)特性研究[D].成都:西南交通大學(xué),2015.

[16]STREET T,LONDON E C.GMRT2100:Railway Group Standard-requirements for Rail Vehicle Structures PartⅡ[S].2012.