Published: 24 September 2018

Informative frequency band identification method using bi-frequency map clustering for fault detection in rotating machines

Jacek Wodecki1
Piotr Kruczek2
Justyna Hebda-Sobkowicz3
Agnieszka Wylomanska4
Radoslaw Zimroz5
Konstantinos Gryllias6
1, 3, 5Diagnostics and Vibro-Acoustic Science Laboratory, Wroclaw University of Science and Technology, Na Grobli 15, 50-421, Wroclaw, Poland
2, 4KGHM Cuprum Ltd, Research and Development Centre, Sikorskiego 2-8, 53-659 Wroclaw, Poland
6Department of Mechanical Engineering, KU Leuven, Celestijnenlaan 300 – Box 2420, 3001 Leuven, Belgium
6Core Lab Dynamics of Mechanical and Mechatronic Systems, Flanders Make, Belgium
Corresponding Author:
Jacek Wodecki
Views 74
Reads 38
Downloads 1315


In presented work the problem of local damage detection in rolling element bearings is addressed. Usually such issues require the usage of the techniques of decomposition, separation etc. In such real industrial cases main difficulty lies in relatively low signal-to-noise ratio as well as unpredictable distribution of damage-related information in the frequency domain, hence the typical methods cannot be used. In this paper such industrial scenario is addressed and a simple yet effective approach to underlying component extraction will be discussed. Proposed method analyzes Cyclic Spectral Coherence map as starting data representation, and Expectation-Maximization is used as analytical tool to determine the informative frequency band (IFB) for impulsive component localization in the carrier frequency spectrum. Finally, based on identified IFB, the bandpass filter is constructed to extract the impulsive component from the input signal.

1. Introduction

he area of machine diagnostics is often based on the analysis of vibration signals. In such signals, damage-related information is located in an Informative Frequency Band (IFB) [1]. The width, the location and the consistency of the band depends on numerous factors, i.e. the type of damage, the operational conditions, the kinematic structure and the damage location in the machine. In many cases those factors and information about the IFB are not known. Regarding such constraints, it is understandable that research on data-driven methods is a highly pursued direction. Various methods for IFB detection have been proposed over the last years, such as the usage of spectral kurtosis [2], selectors [3], the protrugram [4], and many others. However, it is very important that not only IFB-related techniques have been successfully developed over recent years. Examples of other directions and approaches include the selective time analysis in time-frequency domain [5], the multichannel information fusion [6], the time-varying filtration [7], or the heavy-tailed distributions modeling [8-10]

In the presented work the authors attempt to determine the IFB based on a bi-frequency map clustering. In this case the modulation spectra vectors are considered as points in the multidimensional feature space and each carrier frequency bin holds a single point. An Expectation-Maximization technique is used to cluster such data, which allows to obtain a precisely bounded IFB. The boundaries of the appropriate cluster are used to design a bandpass filter allowing the extraction of the damage-related impulsive component from the signal.

2. Methodology

The proposed procedure consists of three main steps. Firstly, the bi-frequency representation of the vibration signal is calculated. In particular the Cyclic Spectral Coherence is used. Secondly the Expectation-Maximization technique is used to cluster individual modulation spectra (for each carrier frequency bin). Finally, the boundaries of the cluster of the CSC map, that can be interpreted as containing the informative frequency band, are used to construct a bandpass filter for impulsive component extraction.

All these techniques are well-known, however their combination is proposed as a successful approach in condition monitoring, focusing specifically on cases, where an impulsive damage-related component investigation is needed but on the other hand a meaningful determination of the IFB is difficult using simpler methods. In this section the definitions and the main properties of the proposed methodology are briefly introduced. Their fusion is illustrated and described afterwards.

2.1. Cyclic spectral coherence

The Cyclic Spectral Coherence is a function, which depends on the modulation and the carrier frequency. It has been introduced by Antoni in 2007 [11] and it has been proven as a powerful tool to indicate the cyclostationary property of signals. A detailed description of this method can be found in [12].

2.2. Expectation-Maximization clustering algorithm

The Expectation-Maximization (EM) algorithm is an optimization method for the estimation of parameters for given measurement data X [13, 14]. EM is especially targeted for separating mixed Gaussian distributions (or any other distribution) in the considered feature space. It is constructed using two main steps which are iterated alternately until achieving the convergence. Clustering in multiple dimensions is the most classic application of this method.

Fig. 1Raw input signal

Raw input signal

It is very important to be aware of the limitations of the EM methodology. It only approximates the maximum likelihood estimate and does not guarantee finding the theoretical value. If the function of likelihood is not concave, EM does not guarantee to find the global maximum of the function. In practice, it should not be executed only once. Usually one should repeat the EM several times with random initialization and select the one presenting the more effective result. In this paper EM has been run 10 times.

3. Application to industrial data

In this article, a vibration signal measured on a rolling bearing is investigated. The bearing is a part of a drive pulley installed in a belt conveyor driving station. The signal is presented in Fig. 1. The sampling frequency has been selected equal to 19.2 kHz. The signal contains a single wideband impulsive component corresponding to local damage of the outer race. The interested readers may find more detailed information regarding this machine to the Ref. [15].

Fig. 2Top: Cyclic Spectral Coherence map of the input signal; Middle: Layout of obtained clusters along the modulating frequency dimension; Bottom: Bandpass filter designed based on obtained clustering results

Top: Cyclic Spectral Coherence map of the input signal; Middle: Layout of obtained clusters along the modulating frequency dimension; Bottom: Bandpass filter designed based on obtained clustering results

The Cyclic Spectral Coherence map of the signal is presented in the top panel of Fig. 2. In the middle panel there is a visual representation of the clustering results, with the carrier frequency axis aligned with the CSC map. EM correctly focuses on the essential part of the carrier frequency band that contains the information about the damage component. The lower and the higher cutoff frequencies of this cluster have been identified equal to 1204 Hz and 5005 Hz respectively. Based on those values, an equiripple bandpass Finite Impulse Response (FIR) filter has been designed using a filter design tool available in Matlab software (see Fig. 2 bottom panel).

Fig. 3Extracted damage-related impulsive component

Extracted damage-related impulsive component

The obtained filter has been used to filter the input vibration signal. As a result, the impulsive component has been extracted allowing for a clear and precise observation of damage-related impulsiveness, which is present in the signal.

4. Conclusions

In this paper the authors presented a new method for local damage detection using bi-frequency map clustering for informative frequency band identification, applied to a real vibration signal emitted by a faulty bearing. Instead of using clustering of two-dimensional data representation for recognizing the location of individual impulses in time (in case of using time-frequency maps as a base domain) or in modulating frequency domain (for bi-frequency maps), the Expectation-Maximization algorithm is used in order to identify the Informative Frequency Band on the Cyclic Spectral Coherence map, which is then used to construct the bandpass filter for the processing of the original signal. The results demonstrate the effectiveness of the method to successfully identify the impulsive component of the signal related to the machine fault.


  • Obuchowski J., Wyłomańska A., Zimroz R. Selection of informative frequency band in local damage detection in rotating machinery. Mechanical Systems and Signal Processing, Vol. 48, 2014, p. 138-152.
  • Antoni J. Cyclostationarity by examples. Mechanical Systems and Signal Processing, Vol. 23, 2009, p. 987-1036.
  • Wyłomańska A., Zimroz R., Janczura J., Obuchowski J. Impulsive noise cancellation method for copper ore crusher vibration signals enhancement. IEEE Transactions on Industrial Electronics, Vol. 63, 2016, p. 5612-5621.
  • Barszcz T., Jabłoński A. A novel method for the optimal band selection for vibration signal demodulation and comparison with the kurtogram. Mechanical Systems and Signal Processing, Vol. 25, 2011, p. 431-451.
  • Wodecki J., Zdunek R., Wyłomańska A., Zimroz R. Local fault detection of rolling element bearing components by spectrogram clustering with semi-binary nmf. Diagnostyka, Vol. 18, 2017, p. 3-8.
  • 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.
  • Kruczek P., Obuchowski J., Wylomanska A., Zimroz R. Cyclic sources extraction from complex multiple-component vibration signal via periodically time varying filter. Applied Acoustics, Vol. 126, 2017, p. 170-181.
  • Wyłomańska A., Żak G., Kruczek P., Zimroz R. Application of tempered stable distribution for selection of optimal frequency band in gearbox local damage detection. Applied Acoustics, Vol. 128, 2017, p. 14-22.
  • Żak G., Wyłomańska A., Zimroz R. Data-driven vibration signal filtering procedure based on the α-stable distribution. Journal of Vibroengineering, Vol. 18, Issue 2, 2016, p. 826-837.
  • Żak G., Wyłomańska A., Zimroz R. Application of alpha-stable distribution approach for local damage detection in rotating machines. Journal of Vibroengineering, Vol. 17, Issue 6, 2015, p. 2987-3002.Antoni J. Cyclic spectral analysis of rolling-element bearing signals: facts and fictions. Journal of Sound and Vibration, Vol. 304, 2007, p. 497-529.
  • Wodecki J., Barszcz T., Zimroz R. Informative frequency band identification for automatic extraction of impulsive components in vibration data from rotating machinery. Vibroengineering Procedia, Vol. 13, 2017, p. 109-114.
  • Sundberg R. Maximum likelihood theory for incomplete data from an exponential family. Scandinavian Journal of Statistics, Vol. 1, 1974, p. 49-58.
  • Dempster A. P., Laird N. M., Rubin D. B. Maximum likelihood from incomplete data via the em algorithm. Journal of the royal statistical society. Series B (Methodological), Vol. 39, 1977, p. 1-38.
  • Wodecki J., Michalak A., Zimroz R. Optimal filter design with progressive genetic algorithm for local damage detection in rolling bearings. Mechanical Systems and Signal Processing, Vol. 102, 2018, p. 102-116.

Cited by

A Diagnostics of Conveyor Belt Splices
Tomasz Kozłowski | Jacek Wodecki | Radosław Zimroz | Ryszard Błażej | Monika Hardygóra

About this article

08 September 2018
15 September 2018
24 September 2018
Fault diagnosis based on vibration signal analysis
local damage detection
cyclic spectral coherence

The work of J. Wodecki and J. Hebda-Sobkowicz was supported by the statutory grant.