王超峰, 林建輝
(西南交通大學牽引動力國家重點實驗室,四川成都 610031)
高速列車在行駛過程中,一般受到軌道不平順、軌道變形、風載荷以及車輪踏面等造成的激勵。而且在兩車交會或是在進出隧道時,可能會在車體表面形成峰值很高的瞬態壓力波[1]。這些激勵都會對車門構件帶來很大的安全影響。所以,對高速列車車門進行模態分析是結構優化設計與結構損傷檢測基礎,具有十分重要的意義。
高速列車運行過程中,在車門的響應存在的非平穩的隨機信號、有局部斷點的信號及一些不能用傅里葉變換來分析的信號,可以利用小波變換將信號分解到不同尺度,以達到模態自動解耦的效果。
模態分析實質上是一種坐標的線性變換,將振動系統以物理坐標和物理參數所描述的、互相耦合的運動方程組,能夠變為一組彼此獨立的方程(每個獨立方程只含一個獨立的模態坐標)[2-3]。
經過離散化處理后的多自由度系統,它在物理坐標系統中的運動微分方程實際上為結構動力學模態分析中的特征值問題。一個n自由度系統的運動微分方程為:

對該方程進行變換可得:

此為一個廣義特征值問題。對于n自由度系統,求解該方程便可確定ωn和{φn},即特征解,其中特征值ω1, ω2,…, ωn代表系統對應的n個固有頻率(正交模態)。當結構振動時,它在任何時候的運動都是正交模態的線性組合。
對多自由度黏性阻尼振動系統自由振動響應作小波變換:

由于小波變換為線性變換,可以交換上式積分和求和的順序,同時將Ai(t)。
在t=τ附近的泰勒級數展開[4]并忽略無窮小量,得到:

不同的尺度下的小波變換相當于用一組帶通濾波器對信號進行帶通濾波處理[5],選擇不同的尺度參數,使小波基的頻窗中心近似等于某階模態頻率,這樣就實現了模態濾波和系統解耦。所以,解耦后的各獨立模態的小波變換可表示為:

根據上式,多自由度系統黏性阻尼系統自由振動響應信號經過小波變換后[6],由其模及相位可以得到:

系統第i階振型在某一時刻t=τt每一個自由度的自由振動響應的小波變換系數,從而可以得到該多自由度系統的第i階振型[7-8]為:

對i取不同的值即可得到不同階的振型。
現對某一型列車車門進行運行模態分析。測試系統包括美國National Instruments(NI)公司的數字采集器2個,型號為NI16255,精度級為±0.1%;Lance(朗斯)公司的加速度計12個,型號為LC0711A,精度級為±0.5%;采用可視化編程的LabVIEW對數據進行顯示。通過便攜式筆記本電腦,實時觀測數據采集情況。利用手持式GPS測定客車的運行速度。
在車門處布置12個加速度傳感器,具體布置如圖1所示。

圖1 車門測點布置
測點5測試結果經過隨機減量法處理后的連續小波變換譜圖如圖2所示。

圖2 小波譜圖
從圖2可以看出,譜圖上的4條脊線分別對應車門的前四階頻率。而且,隨著時間推移,小波脊線上的系數的能量逐漸減小,這與信號是自由振動衰減響應的情況是一致的。通過讀取相應的脊線在所對應頻域的位置,即可得到車門的前四階模態頻率為:42.8 Hz,76.8 Hz,112.1 Hz,162.2 Hz。
車門前兩階歸一化振型分別如圖3、圖4所示。

圖3 第一階振型

圖4 第二階振型
可以看出,基于連續小波變換的模態識別方法很好的解耦了車門的模態參數。
通過以上理論分析和實驗驗證,可以利用小波變換,可以對高速列車車門進行模態參數識別。但該方法在強噪聲背景下的識別效率需要進一步的研究。此外,現階段的模態識別,均假設結構的振動特性是線性的,而且假設外部激勵是平穩的。所以,這種方法在結構非線性和外部激勵多樣化的模態解耦效果仍有待考慮研究與改進。
[1] 雷明準,張豐利,王建楠,陳劍.基于有限元的車門模態分析與優化研究[J].設計計算研究,2008(12):4-7.
[2] 朱宏平.結構損傷檢測的智能方法[M].北京:人民交通出版社,2009:28-30.
[3] 李宏男.結構健康檢測[M].大連:大連理工大學出版社,2005:246-260.
[4] 伊廷華,李宏男,王國新.基于小波變換的結構模態參數識別[J].振動工程學報,2003,19(1):51-56.
[5] Robert X Gao,Ruqiang Yan.Wavelets theory for manufacturing[M].New York:Springer,2011:205-207.
[6] 劉寧.系統模態參數辨識的連續小波方法研究[D].天津:天津大學,2000.
[7] 羅光坤,張令彌.基于Morlet小波變換的模態參數識別研究[J].振動與沖擊,2007,26(7):135-138.
[8] 朱宏平,翁順.運用小波分析方法進行結構模態參數識別[J].振動與沖擊,2007, 26(4):1-4.