基于小波降噪和短時傅里葉變換的主軸突加不平衡非平穩信號分析

胡振邦, 許睦旬, 姜歌東,2, 張東升

(1.西安交通大學 機械工程學院，西安 710049; 2.西安交通大學 機械制造系統工程國家重點實驗室，西安 710049)

現代制造業向著高速高效的方向不斷發展，裝備制造起著至關重要的決定性作用，電主軸是高速加工工藝中最為主要的核心部件，其工作性能將直接影響到裝備制造的水平[1]，由此可見，通過各種技術手段提高電主軸工作性能指標具有極其重要的意義[2]。由于材質、加工誤差、安裝誤差等多種因素導致的主軸不平衡是電主軸最為常見的故障，通過改善電主軸不平衡狀態能夠有效提高電主軸加工精度[3-4]，因此，主軸動平衡技術越來越受到人們的關注。目前，主軸在線動平衡技術大多采用試重的方式用來確定主軸失衡量與振動響應之間的關系[5-6]，對主軸在線添加試重可視為對主軸突加不平衡，由此可以改變主軸振動響應狀態，造成主軸振動信號變化，此過程的振動信號具有極其明顯的非平穩性，并且由于電主軸結構復雜，多種噪聲信號混雜其中，很難從測量信號中直接判斷主軸運行狀態，傳統的傅里葉變換難以對這種含有大量噪聲的非平穩信號進行分析[7]。基于此種情況，本文提出一種采用小波降噪和短時傅里葉變換相結合對主軸突加不平衡非平穩信號進行分析處理的方法，通過小波降噪對原始測量信號按目標頻帶進行分解，提取轉頻信號所在的頻帶進行信號重構，再利用短時傅里葉變換對重構信號進行分析，提取轉頻信號的振動幅值信息。

1 小波降噪技術

1.1 諧波小波基本原理

小波變換是一種具有明確的時域、頻域函數的復變換，諧波小波最早由Newland[8]提出，該方法相較于其他小波算法具有算法簡單、鎖相功能好、適用范圍廣等優點[9-10]，并且能夠適用于非平穩信號的處理，諧波小波時域對應的函數實部為偶函數，虛部為奇函數，即偶函數和奇函數都是零相濾波器，由此可以判斷諧波小波具有“鎖定”信號相位的能力，這種特點對于主軸振動信號的特征提取具有十分重要的作用。

1.2 諧波小波濾波與信號重構

在實際應用中，為了得到頻率為2πm≤ω≤2πn范圍內的信號，需要對諧波小波進行伸縮和平移，此時諧波小波的函數如下所示[11]：

(1)

其中m，n決定了諧波小波變換的尺度j。

對于實際信號x(t)相對尺度j的諧波小波變換如下所示：

(2)

Wx(m,n,τ)=X(ω)Wm,n(ω)

(3)

信號x(t)在尺度j下的諧波小波變換過程如圖1所示。

圖1 諧波小波變換過程

利用諧波小波在頻域的優良特性，在不同分解層和同一層上的不同頻段的局部頻譜進行細化分析，并對信號進行重構，得到轉頻信號所在頻帶范圍的還原信號。

1.3 濾波效果分析

主軸運行過程中的信號極為復雜，以轉頻信號為主同時伴有多種倍頻信號，突加不平穩時的信號中轉頻信號會發生非線性變化，假定主軸轉速為1 500 r/min，采樣頻率為2 048 Hz，采樣時間為1 s，以如式(4)形式的信號作為仿真信號：

x(t)=x1(t)+x2(t)+x3(t)

(4)

其中x1(t)=3*e-2t*sin(2π*25*t)為變化的轉頻信號，x2(t)=2*sin(2π*150*t)為6倍頻諧波信號，x2(t)=4*sin(2π*250*t)為10倍頻諧波信號，仿真信號x(t)的時域波形及頻譜特性如圖2所示。

圖2 仿真信號時域波形及頻譜特性圖

可見，由于信號中存在大量的倍頻信號，已經無法從時域波形中將轉頻信號進行準確提取，從頻譜特性圖中可以提取到轉頻信號的幅值，然而利用傅里葉變換得到的轉頻信號幅值是定值，而該幅值實際上是隨時間變化的，因此頻譜特性圖無法準確反映出信號的幅值變化過程，對該信號進行諧波小波變換并對轉頻及倍頻信號進行重構后所得的時域波形及頻譜分析結果如圖3所示。

由圖3可以看出，利用諧波小波可以清晰的將轉頻信號從原有信號中分離出來，重構的轉頻信號具有非線性特征，與仿真信號的基本特征相吻合，實現了對于高次諧波信號的濾除，同時采用諧波小波技術可以將不同頻域范圍內的倍頻信號根據需要進行分離重構，可見傳統的傅里葉變換方法僅能夠對平穩信號的頻譜進行分析，而對于非平穩信號的分析能力還有待提高。

圖3 信號分離并重構后的信號

2 短時傅里葉變換

目前在工程應用領域，對于信號的頻譜分析一般都采用快速傅里葉變換的方法，該方法的應用最為廣泛，技術也相對比較成熟，但是快速傅里葉變換方法無法滿足非平穩信號的分析，對非平穩過程進行分析時存在無法避免的缺陷[12]，由于快速傅里葉變換只能給出時域或頻域內的統計平均結果，而不能提供信號時域和頻域局部化的信息，因此采用短時傅里葉變換的方法來克服快速傅里葉變換的不足之處，進而為主軸振動特性分析提供準確可靠的信息[13]。短時傅里葉變換的基本運算公式如下：

(5)

其中gω,τ(t)=g(t-τ)e-jωt為積分核函數，e-jωt為限頻函數，g(t-τ)為限時函數，Gx(ω,τ)則是反映x(t)在τ時刻頻率為ω信號頻率成分。

圖4 短時傅里葉變換基本原理示意圖

由此可以看出，短時傅里葉變換是通過固定的時域窗函數g(t-τ)對目標函數x(t)進行截取，對于截取的局部時域信號進行傅里葉變換，通過改變時窗的起始時刻t來調整時域窗函數g(t-τ)的位置，進而得到以t為變化過程的一系列傅里葉變換結果，由這些結果所構成的集合即是短時傅里葉變換的時頻聯合分析結果。

以仿真信號x(t)中1 500 r/min的轉頻信號x1(t)=3*e-2t*sin(2π*25*t)為例，利用短時傅里葉變換計算該振動信號幅值，選用矩形窗函數，窗寬256，移動步長為1，分析結果如圖5所示。

圖5 短時傅里葉變換計算幅值

圖6 振動測試試驗臺

從圖5可以看出，采用短時傅里葉變換能夠對非平穩信號的幅值進行計算，計算值與真實值非常接近，基本滿足應用需要，由此可以知道利用固定時窗函數進行短時傅里葉變換，時域內和頻域內都有較好的局部性，適合對非平穩信號進行分析。

3 應用實例

為驗證本文提出的方法的可行性，采用現有電主軸在線動平衡實驗平臺對方法進行驗證，實驗現場如圖6所示。

該實驗平臺采用150SD40Q7型電主軸，最高設計轉速可達30 000 r/min，最大功率為7 kW，利用在線式動平衡裝置[14]在主軸平穩運行過程中添加不平衡量，利用振動加速度傳感器對主軸振動信號進行測量并換算為振動速度進行處理，振動加速度的測量采用PCB公司的352C33型ICP振動加速度傳感器，信號采集卡采用NI公司出品的4472多路同步信號采集卡，分辨率為24位，實驗用系統采樣頻率為10 kHz，采樣時間為4 s，電主軸測試轉速設定為3 000 r/min。

測試信號的原始時域波形圖如圖7(a)所示，經過傅里葉變換后的頻譜分析結果如圖7(b)所示，可見，在主軸突加不平衡的測試信號中已無法直觀分辨出轉頻信號，同時，由于存在大量的高倍頻信號和多種噪聲的干擾，導致不能夠通過時域波形和頻譜分析結果看出轉頻信號的幅值變化情況，無法在動平衡過程中準確判定主軸振動變化情況。

圖7 突加不平衡測試信號時域波形頻譜

圖8 重構的轉頻信號時域波形及頻譜

利用諧波小波分析對測量信號進行濾波處理，并對轉頻信號進行重構后的信號時域波形如圖8所示，從時域波形可以看出轉頻信號是隨時間變化的，而通過傅里葉變換無法準確表達出轉頻信號的幅值變化情況。

圖9 短時傅里葉變換計算的轉頻信號幅值

選取短時傅里葉變換的時窗寬度為2048，移動步長為1，得到如圖9所示的幅值計算結果。

從實驗結果可以看出，從第1秒開始主軸轉頻振動幅值隨時間推移而不斷的減小并呈現出非平穩變化趨勢，，當到達第3秒時振動幅值逐漸趨于平緩，通過短時傅里葉變換可以準確計算出振動幅值，取得了良好的實驗結果。

4 結 論

主軸平穩運行過程中，突加不平衡將產生非平穩振動信號，并且信號中含有大量的倍頻信號和噪聲信號干擾，無法通過傳統的傅里葉變換對信號進行準確分析，本文對這類問題進行了研究，得到如下結論：

(1)本文提出一種基于小波降噪和短時傅里葉變換相結合對主軸非平穩振動信號進行分析處理的方法，該方法利用諧波小波技術可以將轉頻信號利用時頻分析的方法從原始信號中分離出來并進行信號的重構，實現了突加不平衡非平穩信號的消噪，還原了轉頻信號的變化過程。

(2)采用短時傅里葉變換對降噪后的轉頻信號進行分析計算，獲取轉頻信號振動幅值，能夠直觀、清晰地反映主軸振動幅值變化情況。

(3)該信號處理方法可為在線式動平衡裝置在動平衡調整過程中準確監測振動變化情況提供預警參考，避免不平衡調整量添加過量而導致主軸振動的驟增變化，可以防止由此而發生的主軸嚴重運行故障。

在實際應用過程中，應根據現場需要適當調整短時傅里葉變換的窗函數及移動步長等參數，合理控制算法的計算量，提高應用效率。

參 考 文 獻

[1]吳玉厚.數控機床電主軸單元技術[M].北京:機械工業出版社,2006.

[2]Rodrigues D J, Champneys A R, Friswell M I, et al. Automatic two-plane balancing for rigid rotors[J]. International Journal of Non-Linear Mechanics,2008,43(6): 527-541.

[3]LIAO Yu-he, LANG Gen-feng, WU Fang-ji,et al.An improvement to holospectrum based field balancing method by reselection of balancing object[J]. Journal of Vibration and Acoustics, Transactions of the ASME, 2009,131(3): 0310051-0310057.

[4]ZHANG Yun, ZHANG Dong-sheng, XU Mu-xun, et al. Research on accuracy improvement of holospectrum technique in machine tool spindle[J]. Advances in Communication Technology, 2011, 4:598-603.

[5]Zhou S Y, Dyer Stephen W, Shin K K, et al. Extended influence coefficient method for rotor active balancing during acceleration[J]. Journal of Dynamic Systems, Measurement and Control, Transactions of the ASME, 2004,126(1): 219-223.

[6]章 云,梅雪松,鄒冬林,等.應用模態分析及傅里葉變換的柔性轉子無試重動平衡方法[J].振動與沖擊,2008, 31(11):7-10.

ZHANG Yun，MEI Xue-song, ZOU Dong-lin. A field balancing method of a flexible rotors based on modal analysis and fourier transformation[J]. Journal of Vibration and shock, 2008,31(11): 7-10.

[7]張 鵬,孟 光,吳峰崎.小波分析降噪方法在碰磨和突加不平衡信號特征提取中的應用[J].汽輪機技術,2005, 47(6):441-443.

ZHANG Peng, MENG Guang, WU Feng-qi. The application of wavelet denoising method in the fault feature extraction of friction and sudden unbalance signal[J]. Turbine Technology,2005,47(6): 441-443.

[8]Newland D E. Harmonic wavelet analysis[J]. Proceeding of Royal Society of London, 1993, 443(10):203-205.

[9]孫 偉,熊邦書,黃建萍,等.小波包降噪與LMD相結合的滾動軸承故障診斷方法[J].振動與沖擊,2012,31(18):153-156.

SUN Wei，XIONG Bang-shu，HUANG Jian-ping, et al. Fault diagnosis of a rolling bearing using wavelet packet de-noising and LMD[J]. Journal of Vibration and shock, 2012,31(18):153-156.

[10]李 輝,鄭海起,唐力偉.基于改進雙樹復小波變換的軸承多故障診斷[J].振動、測試與診斷,2013,33(1):53-59.

LI Hui, ZHENG Hai-qi, TANG Li-wei. Bearing multi-faults diagnosis based on improved dual-tree complex wavelet transform[J]. Journal of Vibration, Measurement & Diagnosis,2013,33(1): 53-59.

[11]胡廣書.現代信號處理教程[M].北京:清華大學出版社,2004.

[12]趙鳳展,楊仁剛.基于短時傅里葉變換的電壓暫降擾動檢測[J].中國電機工程學報, 2007,27(10):28-34.

ZHAO Feng-zhan, YANG Ren-gang.Voltage sag disturbance detection based on short time fourier transform[J]. Proceedings of the CSEE, 2007,27(10),28-34.

[13]胡曉依,何慶復,王華勝,等.基于STFT的振動信號解調方法及其在軸承故障檢測中的應用[J].振動與沖擊,2008,27(2):82-86.

HU Xiao-yi, HE Qing-fu, WANG Hua-sheng, et al. Vibration signal demodulation method based on STFT and its application in rolling bearing fault detections [J].Journal of Vibration and shock, 2008,27(2):82-86.

[14]章 云,梅雪松,胡振邦,等.注液式高速切削主軸動平衡裝置設計及其性能研究[J].西安交通大學學報,2013, 47(3):13-17.

ZHANG Yun, MEI Xue-song, HU Zhen-bang, et al. Design and performance analysis of hydrojet-typed balancing device for high-speed machine tool spindle[J].Journal of Xi’an Jiaotong University,2013,47(3):13-17.