Research Article  Open Access
TimeFrequency Analysis Based on Improved Variational Mode Decomposition and Teager Energy Operator for Rotor System Fault Diagnosis
Abstract
A timefrequency analysis method based on improved variational mode decomposition and Teager energy operator (IVMDTEO) is proposed for fault diagnosis of turbine rotor. Variational mode decomposition (VMD) can decompose a multicomponent signal into a number of bandlimited monocomponent signals and can effectively avoid model mixing. To determine the number of monocomponents adaptively, VMD is improved using the correlation coefficient criterion. Firstly, IVMD algorithm is used to decompose a multicomponent signal into a number of monocompositions adaptively. Second, all the monocomponent signals’ instantaneous amplitude and instantaneous frequency are demodulated via TEO, respectively, because TEO has fast and high precision demodulation advantages to monocomponent signal. Finally, the timefrequency representation of original signal is exhibited by superposing the timefrequency representations of all the monocomponents. The analysis results of simulation signal and experimental turbine rotor in rising speed condition demonstrate that the proposed method has perfect multicomponent signal decomposition capacity and good timefrequency expression. It is a promising timefrequency analysis method for rotor fault diagnosis.
1. Introduction
The fault vibration signal of turbine and other largerscale rotating machinery is often characterized by multicomponent and nonstationary on account of large span, complex structure, and variable working condition. Timefrequency analysis method can provide nonstationary signal’s localization information both in time domain and in frequency domain, so it has been widely used in rotating mechanical fault diagnosis [1–10]. The key of timefrequency analysis is to find an effective method to decompose the multicomponent signal. Aiming at this problem, scholars have carried on unceasing exploration. Li and Liang used the generalized synchrosqueezing transform to detect and diagnose faults of the gearbox through analyzing the timefrequency representation [2]. Feng et al. proposed the iterative generalized synchrosqueezing transform to diagnose the lab experimental vibration signals of planetary gearboxes under nonstationary conditions [3]. A timefrequency analysis method based on the VoldKalman filter and higherorder energy separation was proposed by Feng et al., and it was used to diagnose the planetary gearbox fault under varying shaft speed [4]. The cracked rotor transient response was analyzed and monitored using HHT by Ramesh Babu et al. [5]. Liu et al. diagnosed the rotor rubimpact fault using Teager Huang transform [6], but these methods based on Empirical Mode Decomposition (EMD) had the problems of mode mixing, over envelope, less envelope, and so forth [7, 8], which affected the analysis results. Chandra and Sekhar found that Continuous Wavelet Transform (CWT) was more preferred than HHT for noisy data and it was valid to diagnose the shaft misalignment and rotor stator rubimpact faults [8]. Li et al. distinguished the severity degree of rubimpact fault using Empirical Wavelet Transform (EWT) [9]. Tang and Pang extracted the typical fault characteristics of oilwhirl using Hilbert vibration decomposition (HVD) [10].
Dragomiretskiy and Zosso proposed a new multicomponent signal decomposition method called variational mode decomposition (VMD) [11]. This method can decompose a multicomponent signal into a given number of bandlimited intrinsic mode functions using frequency domain nonrecursive iteration pattern. VMD has high precision and can effectively avoid the problem of mode mixing [12]. It has better performance than EWT to produce the classification feature vectors of power quality disturbances [13]. Wang et al. found that the filtering characteristic of VMD was wavelet packetlike expansion via numerical simulations and used it to analyze the fixed frequency rotor rubimpact fault in time domain [14]. However, how to select the appropriate parameter of VMD, especially the number of decomposed components, is very important in practice. To solve this problem, the IVMD which uses the correlation coefficient criterion to determine the number of components adaptively is applied in this paper. Combining the fast and high precision advantages of TEO to monocomponent signal, an IVMDTEO timefrequency analysis new method is proposed. The analysis results of simulation signal and experiment turbine rotor fault signal illustrate the effect of this method.
2. IVMDTEO TimeFrequency Analysis Method
2.1. Theory of VMD
For a multicomposition real valued signal , VMD assumes that is composed of a given number of subsignals (modes). Each mode is regarded as an amplitudemodulated and frequencymodulated (AMFM) signal and has mostly compact frequency around a center pulsation. These modes are called intrinsic mode functions (IMFs). The constrained variational problem is the following:
To assess the bandwidth of the modes, the following scheme is proposed. (1) Compute the associated analytic signal by means of the Hilbert transform to obtain a unilateral frequency spectrum for each mode. (2) Shift the frequency spectrum of each mode to the baseband by mixing with an exponential tuned to the estimated center frequency. (3) Estimate the bandwidth through the Gaussian smoothness of the demodulated signal, that is, the squared norm of the gradient [11].
To solve the constrained variational problem, the augmented Lagrangian is introduced and the nonconstrained variational problem is gotten bywhere α denotes the balancing parameter of the fidelity constraint. The saddle point of (2) is the optimal solution of original problem, which can be gained using alternate direction method of multipliers (ADMM) [15]. All the modes can be obtained from (3) in the frequency domain through updating each mode and its center frequency constantly:
is regarded as the Wiener filtering result of the residue , which makes the VMD algorithm much more robust to sampling and noise. The new center frequency is put at the center of gravity of the corresponding mode’s power spectrum, which can be updated by
2.2. Improved VMD
It is known from VMD algorithm that the number of decomposed components must be set in advance. However, the value is usually difficult to be accurately set for different conditions. In order to solve this problem, the correlation coefficient criterion is introduced to control the number adaptively in this paper, named IVMD.
The correlation coefficient is a quantitative measure of correlation and dependence. It shows the statistical relationships between two or more random variables or observed data values. It has been used in the chromosomes’ similarity analysis [16], the relationship analysis between cerebral venous thrombosis and oral contraceptives in adult women [17], and the mechanical vibration signal similarity analysis [18, 19]. In this paper, Pearson’s correlation coefficient is used and can be measured bywhere is the correlation coefficient of and , is the covariance between and , and and are standard deviations of and .
The correlation coefficient between the decomposed residue of VMD and original signal is computed by (5). When the absolute value is lower than the threshold (0.05 in this paper), it is thought that there is not important information in the residue and the original signal is decomposed completely and appropriately. The iterative decomposition process is finished, and the number can be determined adaptively.
The implementation steps of IVMD are as follows:(1)initial: ;(2)Do_{1}: ;(3)initial: and ;(4)Do_{2}: ;(5)for : , update and using (3) and (4), respectively;(6)update ;(7)repeat_{2} steps (4)–(6), until ;(8)repeat_{1} steps (2)–(7), until .
2.3. Teager Energy Operator
To express timefrequency representation of each mode, the demodulation analysis is necessary to each IMF. TEO is a kind of nonlinear local differential operator without involving integral transform. It can improve the estimation of the instantaneous spectrum characteristics of the vibration data under certain conditions [7]. It has been widely used in bearing fault diagnosis in recent years [20–22].
For the AMFM signal , the corresponding TEO is defined aswhere and are its first and second derivatives, respectively. In the discrete time case, the time derivatives in (6) can be approximated toTEO offers excellent time resolution and demodulation speed because only three sample points are required for the energy computation at each time instant. As a result, the instantaneous change in amplitude and frequency of signal can be quickly tracked. The AM and FM information of vibration signals can be extracted using, respectively [23],
2.4. IVMDTEO TimeFrequency Analysis Method
Combining the advantages of IVMD and TEO, a timefrequency analysis method named IVMDTEO is proposed in this paper. First of all, IVMD algorithm is used to decompose a multicomponent signal into IMFs adaptively. Then, the instantaneous amplitude and instantaneous frequency of each IMF are demodulated by TEO. Finally, the timefrequency representation of the original signal is obtained by superposing the timefrequency representations of all the monocomponents, and the fault types are diagnosed by analyzing the fault characteristic frequency components showed in the timefrequency representation diagram.
In the process of VMD, the balancing parameter α is a parameter to balance the fidelity. The bigger value is suitable to filter lowfrequency component (e.g., trend term) and the smaller value is suitable to detect impacts from a time series (or signal) [24], so the smaller value () is used in this paper.
The timefrequency analysis process of the proposed IVMDTEO method for rotor system is shown in Figure 1, and the specific implementation steps are as follows.(1)Gather original signal data from experimental rotor rig or practical rotor system.(2)Decompose the original signal (multicomposition) into a number of monocompositions ( IMFs) by IVMD method adaptively.(3)Demodulate the instantaneous amplitude and instantaneous frequency of each IMF using TEO.(4)Draw the timefrequency representation diagram of original signal by superposing all timefrequency representations of monocomponents.(5)Diagnose the rotor system fault types by observing the corresponding fault characteristic frequency components that appeared possibly in Figure 1.
3. Simulation Signal Analysis
A multicomponent simulation signal is used to verify the timefrequency analysis effect of the proposed method. is composed of two FM components , two mutation frequency cosine components , an AM cosine component , and a Gaussian noise . The components’ frequency is close to each other.
The sampling frequency is 2048 Hz, and the analysis points are 2048 points. The waveform and spectrum of are shown in Figure 2. The spectrum basically shows each component. First, is decomposed by the IVMD algorithm, and seven IMFs are extracted and shown in Figure 3. IMF1~IMF3 are corresponding to , IMF4 and IMF5 are corresponding to , and IMF6 is corresponding to . The correlation coefficient between the IMFs and is shown in Table 1. The former six IMFs are bigger than the threshold 0.05 and their timefrequency representations demodulated by TEO, respectively, are shown in Figure 4. It can be seen from Figure 4 that IVMD can effectively segregate the two FM components in , the mutation frequency at about 600th point in , and the AM component in . Although there is a little interference, all timefrequency representation of each component is clearly segregated.

(a)
(b)
For comparison, the timefrequency spectrum of simulation signal analyzed by the HHT is shown in Figure 5. It can be seen from Figure 5 that the frequency modulation phenomenon is confused. Two FM components in are extracted partly, but the low disturbance frequency is emerged. There are frequency components in the vicinity of 200 Hz, but the detailed composition cannot be expressed. The analysis result is bad.
4. Experiment Rotor System Vibration Signal Analysis
The turbine rotor faults are simulated on the Bently RK4 test rig (Figure 6). The test rig is an oilfluid supporting rotorbearing system. An electrical motor is attached to the shaft that carries two discs by a flexible coupling. A phase sensor closed to the coupling is fixed to measure the rotating speed. Two eddy current sensors, located near the rub copper thimble, are fixed by holder in the middle of test rig. A safety device is designed to ensure the experiment rotor safety. An oil film journal bearing is part of an assembly connected to the oil pump system. The bearing is equipped with two proximity eddy current sensors in the horizontal and vertical directions. Rubimpact signal and oil film instability signal are gathered by sensor 1 and sensor 2, respectively using ZonicBook/618E vibration measurement system [6]. The sampling frequency is 5120 Hz. In order to reduce the amount of analysis data, the following analysis frequency is 512 Hz. That is to say, one point in ten is selected in sequence.
4.1. RubImpact Fault Analysis
The copper thimble that has a small separation from the shaft is fixed by holder in the middle of test rig and rubs against the highspeed rotating shaft slightly. The rubimpact signal when shaft speed is 1620 r/min~1740 r/min is analyzed and its waveform and spectrum are shown in Figure 7. The waveform amplitude becomes bigger with time, which only shows the rising rotational speed process. The spectrum only shows that the continuous rotational frequency is 27~29 Hz. It is difficult to diagnose the slight rubimpact fault accurately.
(a)
(b)
Seven IMFs obtained by IVMD for slight rubimpact signal are shown in Figure 8. In apparent periodicity components, IMF5 is corresponding to the rotational frequency (), IMF3 is corresponding to , IMF2 is corresponding to , and IMF1 is corresponding to . The correlation coefficients of the former six IMFs are bigger than the threshold value. Their timefrequency representations demodulated by TEO are shown in Figure 9. It can be seen that the amplitude of rotational frequency is the maximum and frequency components are also extracted by IVMDTEO. At the same time, the high harmonics frequency fluctuates in the range of 500th to 1000th point. It illustrates that the rubimpact process is instable. All of the frequency characteristics above accord with the rubimpact fault characteristics [25], so the occurrence of rotor slight rubimpact fault can be diagnosed accurately.
In order to compare the analysis results, the timefrequency analysis result for slight rubimpact signal by HHT is shown in Figure 10. It only shows demodulated rotational frequency and cannot show the high harmonics frequency. The slight rubimpact characteristics are not extracted. The analysis effect is not as good as the effect of IVMDTEO.
4.2. Oil Film Instability Analysis
The oil film instability of the rotor system usually consists of two stages that are the nearly half speed oilwhirl and oilwhip. The monitor and control to initial oilwhirl are significant for turbine rotor to avoid great loss. The initial oilwhirl and the firstorder oilwhip spectrum characteristics of the test rig are analyzed subsequently.
In the process of oil film instability, the oilwhirl occurs at the rising speed of 1620 r/min~1860 r/min, and the waveform and spectrum of the oilwhirl are shown in Figure 11. The rising speed can be seen in the waveform. In the spectrum, the oilwhirl frequency () (12~14 Hz) and (27~31 Hz) can be seen. The value is about half value of and higher amplitude than , but all of them lack the corresponding time information.
(a)
(b)
Six IMFs analyzed by IVMD to initial oilwhirl signal are shown in Figure 12. In apparent periodicity components, IMF1 is corresponding to . IMF2 and IMF3 are corresponding to . The correlation coefficients of the former five IMFs are bigger than the threshold value. Their timefrequency representations demodulated by TEO, respectively, are shown in Figure 13. It can be seen that a little oilwhirl amplitude appears before 1000th point and the strong oilwhirl amplitude accompanied by amplitude modulation phenomenon appears after the 1000th point. The process shows that the initial oilwhirl is instability process. The timefrequency representation clearly shows the whole process of the occurrence and development of oilwhirl. A phenomenon that a frequency component whose frequency is about before the 1200th point and becomes after the 1200th point appears too. This detailed changing process suggests the high precision decomposition of IVMDTEO.
For comparison, the timefrequency spectrum of oilwhirl signal analyzed by HHT is shown in Figure 14. The frequency mixing phenomenon is obvious, which cannot clearly show the whole process of the oilwhirl.
When the rotational speed is up to 2760 r/min~3000 r/min, the rotor system test rig shows biggish vibration because of the oilwhip. Waveform and spectrum of the oilwhip signal are shown in Figure 15. The rotational frequency () (46~50 Hz) and the oilwhip frequency () (22~24 Hz) can be seen in the spectrum, but they are continuous spectrum and lack corresponding time information. The oilwhip can be diagnosed preliminary.
(a)
(b)
Eleven IMFs obtained by IVMD to oilwhip signal are shown in Figure 16. The correlation coefficients of the former ten IMFs are bigger than the threshold value. Their timefrequency representations demodulated by TEO, respectively, are shown in Figure 17. Besides the main and frequency components, the related combination frequency components , , , , , and are also extracted clearly. The oilwhip can be diagnosed accurately and reliably [25]. The analysis results also indicate the decomposition accuracy of the IVMDTEO method.
To compare the analysis results, the timefrequency spectrum of oilwhip signal analyzed by HHT is shown in Figure 18. The , , and frequency components are observed only and other frequency components cannot be extracted out. The analysis results are bad to IVMDTEO method.
Qualitative analysis results in Figures 9, 13, and 17 for rotor faults show better effectiveness of the proposed method than HHT. The quantitative description in main fault frequency amplitude’s change and detailed fault frequency amplitude’s value to every fault are shown in Table 2. For slight rubimpact fault, HHT method cannot extract fault frequency components; however, the proposed method can extract fault frequency components and their amplitudes are about 4 μm. For oil film instability fault (oilwhirl and oilwhip), most of the main fault frequency amplitudes extracted by proposed method, such as oilwhirl, , and , are improved within 12.3%～35.5% comparing with HHT except oilwhip frequency. The new extracted detail fault frequency amplitudes which cannot be extracted by HHT are also shown in Table 2. Both qualitative analysis and quantitative analysis can illustrate the advantage of IVMDTEO method.
 
Note: “—” represents no extracted fault frequency by HHT and “new” represents new extracted fault frequency by IVMDTEO. 
5. Conclusion
VMD is improved to determine the number of decomposition components adaptively using the correlation coefficient criterion. Integrating the advantages of IVMD and TEO, a timefrequency analysis method named IVMDTEO is proposed in this paper. The analysis results of experimental turbine rotor signals under variable speed conditions illustrate that the method can show the timefrequency characteristics of slight rubimpact clearly, monitor the process of oilwhirl occurrence and development exactly, and extract the complicated timefrequency components of oilwhip accurately. It enhances the precision and accuracy of rotor fault diagnosis. Comparative analysis results show that the proposed method is more effective than HHT, so this method has important engineering application value for rotor fault diagnosis.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (no. 51307058), the Natural Science Foundation of Hebei Province, China (nos. E2014502052 and E2015502013), and the Fundamental Research Funds for the Central Universities (nos. 2014MS156, 2015ZD27, and 2015XS120).
References
 Y. Qin, “Multicomponent AMFM demodulation based on energy separation and adaptive filtering,” Mechanical Systems and Signal Processing, vol. 38, no. 2, pp. 440–459, 2013. View at: Publisher Site  Google Scholar
 C. Li and M. Liang, “Timefrequency signal analysis for gearbox fault diagnosis using a generalized synchrosqueezing transform,” Mechanical Systems and Signal Processing, vol. 26, no. 1, pp. 205–217, 2012. View at: Publisher Site  Google Scholar
 Z. Feng, X. Chen, and M. Liang, “Iterative generalized synchrosqueezing transform for fault diagnosis of wind turbine planetary gearbox under nonstationary conditions,” Mechanical Systems and Signal Processing, vol. 5253, no. 1, pp. 360–375, 2015. View at: Publisher Site  Google Scholar
 Z. Feng, S. Qin, and M. Liang, “Timefrequency analysis based on voldkalman filter and higher order energy separation for fault diagnosis of wind turbine planetary gearbox under nonstationary conditions,” Renewable Energy, vol. 85, pp. 45–56, 2016. View at: Publisher Site  Google Scholar
 T. Ramesh Babu, S. Srikanth, and A. S. Sekhar, “HilbertHuang transform for detection and monitoring of crack in a transient rotor,” Mechanical Systems and Signal Processing, vol. 22, no. 4, pp. 905–914, 2008. View at: Publisher Site  Google Scholar
 S. K. Liu, G. J. Tang, and B. Pang, “Rotor local rubbing fault feature analysis based on TeagerHuang transform,” Advanced Materials Research, vol. 989–994, pp. 3244–3247, 2014. View at: Publisher Site  Google Scholar
 I. Antoniadou, G. Manson, W. J. Staszewski, T. Barszcz, and K. Worden, “A timefrequency analysis approach for condition monitoring of a wind turbine gearbox under varying load conditions,” Mechanical Systems and Signal Processing, vol. 6465, pp. 188–216, 2015. View at: Publisher Site  Google Scholar
 N. H. Chandra and A. S. Sekhar, “Fault detection in rotor bearing systems using time frequency techniques,” Mechanical Systems and Signal Processing, vol. 7273, pp. 105–133, 2016. View at: Publisher Site  Google Scholar
 Z. N. Li, M. Zhu, F. L. Chu, and Y. X. Xiao, “Mechanical fault diagnosis method based on empirical wavelet transform,” Journal of Scientific Instrument, vol. 35, no. 11, pp. 2423–2432, 2014. View at: Google Scholar
 G. Tang and B. Pang, “A mechanical fault diagnosis method based on improved Hilbert vibration decomposition,” Journal of Vibration and Shock, vol. 34, no. 3, pp. 167–182, 2015. View at: Publisher Site  Google Scholar
 K. Dragomiretskiy and D. Zosso, “Variational mode decomposition,” IEEE Transactions on Signal Processing, vol. 62, no. 3, pp. 531–544, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 G. Tang and X. Wang, “Parameter optimized variational mode decomposition method with application to incipient fault diagnosis of rolling bearing,” Journal of Xi'an Jiaotong University, vol. 49, no. 5, pp. 73–81, 2015. View at: Publisher Site  Google Scholar
 C. Aneesh, S. Kumar, P. M. Hisham, and K. P. Soman, “Performance comparison of variational mode decomposition over empirical wavelet transform for the classification of power quality disturbances using support vector machine,” Procedia Computer Science, vol. 46, pp. 372–380, 2015. View at: Publisher Site  Google Scholar
 Y. Wang, R. Markert, J. Xiang, and W. Zheng, “Research on variational mode decomposition and its application in detecting rubimpact fault of the rotor system,” Mechanical Systems and Signal Processing, vol. 60, pp. 243–251, 2015. View at: Publisher Site  Google Scholar
 D. P. Bertsekas, “Constrained optimization and Lagrange Multiplier methods,” in Computer Science and Applied Mathematics, vol. 1, Academic Press, Boston, Mass, USA, 1982. View at: Google Scholar
 C. Someswara Rao and S. Viswanadha Raju, “Similarity analysis between chromosomes of Homo sapiens and monkeys with correlation coefficient, rank correlation coefficient and cosine similarity measures,” Genomics Data, vol. 7, pp. 202–209, 2016. View at: Publisher Site  Google Scholar
 S. M. Zuurbier, S. Middeldorp, J. Stam, and J. M. Coutinho, “Sex differences in cerebral venous thrombosis: a systematic analysis of a shift over time,” International Journal of Stroke, vol. 11, no. 2, pp. 164–170, 2016. View at: Publisher Site  Google Scholar
 Z. H. Zhao, S. P. Yang, and Y. J. Shen, “Machinery fault feature extraction based on independent component analysis and correlation coefficient,” Journal of Vibration and Shock, vol. 5, no. 5, pp. 479–483, 2013. View at: Google Scholar
 Z. Qiao and Z. Pan, “SVD principle analysis and fault diagnosis for bearings based on the correlation coefficient,” Measurement Science and Technology, vol. 26, no. 8, Article ID 085014, 2015. View at: Publisher Site  Google Scholar
 M. Liang and I. S. Bozchalooi, “An energy operator approach to joint application of amplitude and frequencydemodulations for bearing fault detection,” Mechanical Systems and Signal Processing, vol. 24, no. 5, pp. 1473–1494, 2010. View at: Publisher Site  Google Scholar
 P. Henríquez Rodríguez, J. B. Alonso, M. A. Ferrer, and C. M. Travieso, “Application of the TeagerKaiser energy operator in bearing fault diagnosis,” ISA Transactions, vol. 52, no. 2, pp. 278–284, 2013. View at: Publisher Site  Google Scholar
 X. Jiang, S. Li, and C. Cheng, “A novel method for adaptive multiresonance bands detection based on VMD and using MTEO to enhance rolling element bearing fault diagnosis,” Shock and Vibration, vol. 2016, Article ID 8361289, 20 pages, 2016. View at: Publisher Site  Google Scholar
 P. Maragos, J. F. Kaiser, and T. F. Quatieri, “On amplitude and frequency demodulation using energy operators,” IEEE Transactions on Signal Processing, vol. 41, no. 4, pp. 1532–1550, 1993. View at: Publisher Site  Google Scholar
 Y. Wang and R. Markert, “Filter bank property of variational mode decomposition and its applications,” Signal Processing, vol. 120, pp. 509–521, 2016. View at: Publisher Site  Google Scholar
 H. Ma, T. Yu, Q. Han, Y. Zhang, B. Wen, and C. Xuelian, “Timefrequency features of two types of coupled rubimpact faults in rotor systems,” Journal of Sound and Vibration, vol. 321, no. 35, pp. 1109–1128, 2009. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Shangkun Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.