Application of principal component analysis of time-frequency representation for gearbox fault detection
Jacek Wodecki^{1} , Justyna Hebda Sobkowicz^{2} , Agnieszka Wyłomanska^{3} , Radoslaw Zimroz^{4} , Konstantinos Gryllias^{5}
^{1, 2, 4}Diagnostics and Vibro-Acoustic Science Laboratory, Wroclaw University of Science and Technology, Na Grobli 15, 50-421, Wroclaw, Poland
^{3}KGHM Cuprum Ltd, Research and Development Centre, Sikorskiego 2-8, 53-659, Wroclaw, Poland
^{5}Department of Mechanical Engineering, KU Leuven, Celestijnenlaan 300 – box 2420, 3001, Leuven, Belgium
^{5}Core Lab Dynamics of Mechanical and Mechatronic Systems, Flanders Make, Belgium
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 19, 2018, p. 82-85.
https://doi.org/10.21595/vp.2018.20209
Received 5 September 2018; accepted 11 September 2018; published 24 September 2018
Dimensionality reduction methods are very useful and effective tools in the field of data analytics, used either independently or as a pre-processing step in the frames of a complex algorithm. In this paper a simple yet powerful technique for local damage detection in heavy-duty industrial machinery is presented with particular focus on gearboxes. It assumes that the cyclic component present in the vibration signal carrying information about the damage, can be extracted from relevant frequency bands of the signal. Although this assumption is usually a starting point for selective filtration in the notion of Informative Frequency Band (IFB) identification, in this case the frequency bands are not addressed directly. The authors propose to apply Principal Component Analysis (PCA) as a dimensionality reduction method on the time-frequency representation of the input data in such a way, that the dimension of frequency is reduced. In this way, the variance maximized in the first principal component is expected to capture the cyclic information which is related to the damage present in the machine.
Keywords: local damage detection, vibration, principal component analysis, time-frequency analysis.
1. Introduction
The topic of fault detection in rotating machines is still an open problem in the field of machine diagnostics. Several reviews on damage detection in bearings and gears can be found in literature [1-3]. The typical methods take advantage of Higher-Order Statistics [4], the Wavelet Transform [5], the Time-Frequency domain analysis [6, 3], the Bi-frequency analysis [7], etc. The Vibration signals capture on a machinery system are often a mixture of several source signals. For instance, a signal acquired on a bearing operating in a belt conveyor driving station might be contaminated with vibrations from a neighboring gearbox or by vibrations caused by other damage. Another example is the analysis of a multi-source signal that is acquired on a gearbox with two faults, each of different nature [8]. In this paper PCA is used for local damage detection in a two-stage gearbox operating in a belt conveyor driving station. The PCA is applied on the Short-Time Fourier Transform (STFT) of the measured vibration signal in order to reduce its dimensionality, which is a common method for this purpose [9, 10]. The authors propose to focus on the analysis of the first principal component produced by the PCA algorithm. The method allows to discover the waveform that represents the fault component within the overall vibration data.
2. Methodology
Firstly, the signal must be transformed into a time-frequency domain. For such transformation, the spectrogram has been chosen. The Short-time Fourier transform is given for discrete data $x\left[0\right],x\left[1\right],...,x[N-1]$ is given by the formula [11]:
where $0\le k\le N-1$ is the frequency bin, $n$ is the time point and $w[\cdot ]$ is the window of length $L$. One can observe that, in the STFT for each time point the Fourier transform is calculated using the FFT algorithm. Furthermore, the spectrogram is equal to the squared absolute value of the STFT:
In the next step the dimensionality of the obtained spectrogram matrix is reduced using principal component analysis with respect to frequency, so that the output contains vectors in time domain.
Fig. 1. Raw input signal
Principal component analysis is one of the most common and widespread methods for multivariate linear data analysis. It serves for investigating data structure, data mining, data smoothing and approximation as well as for exploring data dimensionality. The method permits to build new features, called Principal Components (PCs), which may serve for visualization of the data [9, 10, 12, 13].
Finally, the first PC is expected to capture the time-domain waveform of the damage component hidden in the signal.
3. Application to industrial data
The vector of observation contains a vibration signal of a two-stage gearbox in a belt conveyor driving station, commonly used in mining industry for material transportation (see Fig. 1). The parameters of the data acquisition are selected as follows: duration 2.5 s, sampling frequency 16384 Hz and the expected fault frequency is equal to 16.5 Hz. Preliminary observations of the raw signal allow to observe a clearly visible amplitude modulation that is not related to the local damage under investigation but is related to a misalignment of the neighboring shaft. In order to decompose the process into more informative sub-processes the STFT has been applied. Hence the time-frequency spectrogram matrix, presented in Fig. 2, is obtained.
Fig. 2. Spectrogram of the input signal
Fig. 3. Fourier spectra of raw signal and first three components
In the spectrogram one can observe three main frequency bands: a first one containing a low frequency with high energy responsible for the shape of the signal set at 0-1.5 kHz, a second one which contains two Informing Frequency Bands (IFB) placed at 2.5-3.5 kHz and 4-4.5 kHz and a third non-informative high frequency band above 5.3 kHz. PCA is applied on the spectrogram matrix and the first principal component is investigated and plotted in Fig. 3 along with the original time series. The extracted principal component contains by definition the most relevant information content of the data. The remaining PCs have been omitted at this work.
The first principal component maximizes the variance, namely it takes into account as much of the variability in the data as possible and therefore the PC1 is treated as the component which is responsible for the fault description. The spectral signatures of the observed PC1 have been investigated in Figure 3. As it was expected, the amplitude spectrum graph of PC1 gives very precise information about the fundamental frequency of the damage, equal to 16.5 Hz.
4. Conclusions
In this article, the authors proposed a simple yet robust approach to the detection of periodically impulsive behaviors in the vibration signal. This behavior is associated with a component with information about the fault. The analyzed data have been acquired from a two-stage gearbox of a complex mechanical system working in mining environment. In order to detect the frequency of the component which has been damaged, principal component analysis has been applied, using a time-frequency representation of the vibration data as the basis for the analysis. As a result, a clear and precise periodical component has been obtained, that can indicate the fault frequency, which is confirmed by the spectrum of the result component. Such information can be further used in identifying the faulty component in the process of the machine maintenance.
Acknowledgements
The work of J. Wodecki and J. Hebda-Sobkowicz was supported by the statutory grant.
References
- Randall R. B., Antoni J. Rolling element bearing diagnostics – a tutorial. Mechanical Systems and Signal Processing, Vol. 25, 2011, p. 485-520. [Publisher]
- Samuel P. D., Pines D. J. A review of vibration-based techniques for helicopter transmission diagnostics. Journal of Sound and Vibration, Vol. 282, 2005, p. 475-508. [Publisher]
- Feng Z., Liang M., Chu F. Recent advances in time-frequency analysis methods for machinery fault diagnosis: A review with application examples. Mechanical Systems and Signal Processing, Vol. 38, 2013, p. 165-205. [Publisher]
- Antoni J., Randall R. The spectral kurtosis: application to the vibratory surveillance and diagnostics of rotating machines. Mechanical Systems and Signal Processing, Vol. 20, 2006, p. 308-331. [Publisher]
- Lin J., Zuo M. Gearbox fault diagnosis using adaptive wavelet filter. Mechanical Systems and Signal Processing, Vol. 17, 2003, p. 1259-1269. [Publisher]
- Burdzik R., Konieczny L., Folęga P. Structural health monitoring of rotating machines in manufacturing processes by vibration methods. Advanced Materials Research, Vol. 1036, 2014, p. 642-647. [Publisher]
- Borghesani P., Pennacchi P., Chatterton S. The relationship between kurtosis-and envelope-based indexes for the diagnostic of rolling element bearings. Mechanical Systems and Signal Processing, Vol. 43, 2014, p. 25-43. [Publisher]
- Żak G., Obuchowski J., Wyłomańska A., Zimroz R. Novel 2d representation of vibration for local damage detection. Mining Science, Vol. 21, 2014, p. 105-113. [Search CrossRef]
- Bartkowiak A., Zimroz R. Dimensionality reduction via variables selection–linear and nonlinear approaches with application to vibration-based condition monitoring of planetary gearbox. Applied Acoustics, Vol. 77, 2014, p. 169-177. [Publisher]
- Bartkowiak A., Zimroz R. Data dimension reduction and visualization with application to multi-dimensional gearbox diagnostics data: comparison of several methods. Solid State Phenomena, Vol. 180, 2012, p. 177-184. [Publisher]
- Allen J. Short term spectral analysis, synthesis, and modification by discrete Fourier transform. IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. 25, 1977, p. 235-238. [Publisher]
- Moore B. Principal component analysis in linear systems: Controllability, observability, and model reduction. IEEE Transactions on Automatic Control, Vol. 26, 1981, p. 17-32. [Publisher]
- Wodecki J., Stefaniak P., Obuchowski J., Wyłomańska A., Zimroz R. Combination of principal component analysis and time-frequency representations of multichannel vibration data for gearbox fault detection. Journal of Vibroengineering, Vol. 18, 2016, p. 2167-2175. [Publisher]