Bearings damage detection in presence of heavy non-Gaussian noise via cyclo-stationary analysis
Radoslaw Zimroz1 , Jakub Obuchowski2 , Agnieszka Wyłomańska3
1, 2Diagnostics and Vibro-Acoustic Science Laboratory, Wroclaw University of Technology,, Wroclaw, Poland
3Hugo Steinhaus Center, Institute of Mathematics and Computer Science, Wroclaw University of Technology, Wroclaw, Poland
Vibroengineering PROCEDIA, Vol. 3, 2014, p. 88-92.
Accepted 4 August 2014; published 10 October 2014
A problem of rolling element bearings diagnostics is discussed in the paper. In the classical approach, time domain detectors (e.g. kurtosis) or frequency domain representation (envelope spectrum analysis) are used to identify a damage in bearings. However, in our case, the machine (copper ore crusher) produces randomly spaced in time heavy impulsive disturbances related to normal operation of the machine. It eliminates completely possibility of detection of damage in time domain or using analysis of the envelope spectrum of the raw vibration signal. It is well known that effectiveness of envelope analysis might be significantly improved if the signal is pre-processed before enveloping. Unfortunately, again, almost all known criteria based on maximization of impulsiveness of signal in time domain fail in this case. In the paper we propose to incorporate cyclostationarity instead of impulsiveness, namely we propose to perform the spectral coherence density analysis.
Keywords: rolling element bearing, damage detection, cyclostationary analysis, non-Gaussian noise.
A problem of rolling element bearings diagnostics is discussed in the paper. In the classical approach, time domain detectors (e.g. kurtosis) [1, 3-9] or frequency domain representation (envelope spectrum analysis ) are used to identify damage in bearings . However, for damage detection at early stages , when signal to noise ratio is poor/informative signal is weak , damage detection becomes difficult and requires advanced signal processing including decomposition, modeling, adaptive filtering, etc. [10-14]. In our case, the machine (copper ore crusher) produces randomly spaced in time heavy impulsive disturbances related to normal operation of the machine, similarly as reported in [15-16]. It makes impossible detection of damage in time domain, moreover due to very poor signal to noise ratio, even envelope spectrum analysis is not able to detect damage directly based on raw vibration. It is well known that effectiveness of envelope analysis might be significantly improved if signal is pre-processed before enveloping. Unfortunately, impulsiveness cannot be basis for signal extraction anymore. So, the commonly used kurtosis-based optimization of pre-filtering (kurtogram, spectral kurtosis) cannot be applied here. Moreover, it is known that damage signature depends on load applied to the machine. Thus we claim that these non-periodic, strong impulses are non-informative from diagnostic point of view and should be separated from informative part. After that, informative, cyclic impulses related to the damaged part of REB should be analysed with taking into account operating condition. However, first step is to detect if there is any damage. Promising approach is to search for cyclic/periodic nature of impulses [17-20]. An envelope analysis without pre-filtering do not provide acceptable results, so we focused on cyclo-stationary approach, namely spectral coherence density map [16-20].
2. Machine and experiment description
The machine discussed here is a heavy duty crusher used in a mineral processing plant. The crusher is a kind of a machine which use a metal surface to crumble materials into small fractional pieces. During this process, as well as during entering material stream into the crusher, a lot of impacts/shocks appear. They are present in vibration signal acquired from the bearing’s housing. Several crushers were investigated. Vibration from bearings housing in horizontal and vertical direction and tacho signal have been measured using Endevco accelerometers Bruel Kjaer Laser probe, NI DAQ card and Labview Signal Express. Signals were acquired for idle mode as well as for different values of load (volume and granulation of material stream) with approx 120 s duration each. Sampling frequency was set to 25 kHz.
Fig. 1. A crusher – general view (note bearings with yellow housing)
3. Preliminary analysis
As mentioned, crushers are operating under non-stationary operating conditions that results in different values of amplitudes as well as frequency structure of acquired signals (see Figs. 2, 3).
It is impossible to recognize fault based on time series or time-frequency map due to number of impulses related to crushing process other “normal” operation of the machine. As it was mentioned, due to difficulties with automatic and even manual (based on STFT) informative frequency band selection, an envelope analysis was performed without pre-filtering. Fig. 4 shows envelope and its spectrum. Based on documentation one should expect signature at fault frequency equal to 30 Hz. From Fig. 4 it is difficult to recognize the damage.
4. Application of spectral coherence density for detection of cyclic impulsive disturbance in presence of non-Gaussian impulsive noise
The vibration signal from faulty bearings, so also our considered signals here are expected to exhibit the so called second-order cyclostationarity property (also called periodically correlated property, PC). Let us mention, the periodically correlated random processes of second-order are random systems in which there exists the periodic rythm in the structure that is generally more complicated than periodicity in the mean function of given signal . We say the signal is cyclostationary with period if for each and the following properties are satisfied :
And there are no smaller values of 0 for which the above conditions hold.
The periodically correlated processes posses very interesting properties. Therefore they found many practical applications like in hydrology , meteorology , economics  and electrical engineering .
Fig. 2. Raw vibration signal a) whole observation and b) part of the signal – without and with load applied (top/bottom respectively)
Fig. 3. Time-Frequency map (STFT) of vibration signal idle mode a) and b) loaded machine
Fig. 4. Envelope and envelope spectrum of part of the signal (load was applied)
b) Envelope spectrum of raw signal
Let us mention, one of the model that is can be used to description of the periodically correlated processes is periodic autoregressive moving average model (PARMA) which is an alternative to classical ARMA models uesd to description of stationary signals, .
If the signal is considered to be second-order cyclostationary then, according to the definition of PC process, the autocorrelation function:
is periodic with respect to and then it has a Fourier series expansion:
where the Fourier coefficients are known as the cyclic autocorrelation functions.
The Fourier transform of the cyclic autocorrelation function at a given cyclic frequency then defines the cyclic power spectrum (CPS):
If one consider set of CPS’s for all , then we can get Spectral Correlation Density Map (SCDM):
In order to get a normalized version of SCDM it was proposed to use so called Spectral Coherence Map (SCM):
As shown in Fig. 5 at frequency 30 Hz and its 2nd and 3rd harmonics vertical lines appeared. They correcponds to damage signature. So, one might conclude that using SCDM representation detection of the damage is possible. Moreover there are also other excitations that are not identified at this stage. It will be a task for furhter research.
It is clear that non-identified components mask damage signature. As reported in , there is a chance to „enhance“ SCD Map, as well as to extract signal of interest using advanced blind filtering algorithms .
Fig. 5. Spectral coherence density map (note excitation at 30 Hz and its harmonics)
In the paper preliminary results regarding application of Spectral Coherence Density for detection of cyclic impulsive disturbance in presence of non-Gaussian impulsive noise are presented. The vibration signals come from the crusher – a machine used in a mineral processing plant for pre-processing of minerals. During operation of the machine a lot of impacts are generated due to falling down of material stream and and the crushing process. Thus, detection of damaged bearings is very difficult. It has been proved that known time domain indicators, time-frequency representation or even comonly used envelope spectrum fail. The Spectral Coherence Density allows to detect modulating frequency alpha and its harmonics that correspond to the bearing‘s fault frequency. Further research will be focused on extraction of the signal of interest.
- Randall R. B., Antoni J. Rolling element bearing diagnostics-A tutorial. Mechanical Systems and Signal Processing, Vol. 25, Issue 2, 2011, p. 485-520. [Search CrossRef]
- Łazarz B., Wojnar G., Czech P. Early fault detection of toothed gear in exploitation conditions. Maintenance and Reliability, Vol. 49, Issue 1, 2011. [Search CrossRef]
- Antoni J. Fast computation of the kurtogram for the detection of transient faults. Mechanical Systems and Signal Processing, Vol. 21, 2007, p. 108-124. [Search CrossRef]
- 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. [Search CrossRef]
- 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, Issue 1, 2011, p. 431-451. [Search CrossRef]
- Obuchowski J., Wylomanska A., Zimroz R. Stochastic modeling of time series with application to local damage detection in rotating machinery. Key Engineering Materials, Vol. 569-570, 2013, p. 441-448. [Search CrossRef]
- Obuchowski J., Wyłomańska A., Zimroz R. The local maxima method for enhancement of time-frequency map and its application to local damage detection in rotating machines. Mechanical Systems and Signal Processing, Vol. 46, Issue 2, p. 389-405. [Search CrossRef]
- 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, Issue 1-2, p. 138-152. [Search CrossRef]
- Lin J., Zuo M. Gearbox fault diagnosis using adaptive wavelet filter. Mechanical Systems and Signal Processing, Vol. 17, Issue 6, p. 1259-1269, 2003. [Search CrossRef]
- Czech P., Wojnar G., Burdzik R., Konieczny Ł., Warczek J. Application of the discrete wavelet transform and probabilistic neural networks in IC engine fault diagnostics. Journal of Vibroengineering Vol. 16, Issue 4, 2014, p. 1619-1639. [Search CrossRef]
- Zimroz R., Bartelmus W. Application of adaptive filtering for weak impulsive signal recovery for bearings local damage detection in complex mining mechanical systems working under condition of varying load. Diffusion and Defect Data Pt.B: Solid State Phenomena, Vol. 180, 2012, p. 250-257. [Search CrossRef]
- Makowski R. A., Zimroz R. Adaptive bearings vibration modeling for diagnosis. Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics), 2011, p. 248-259. [Search CrossRef]
- Makowski R., Zimroz R. A procedure for weighted summation of the derivatives of reflection coefficients in adaptive Schur filter with application tofault detection in rolling element bearings. Mechanical Systems and Signal Processing, Vol. 38, Issue 1, 2013, p. 65-77. [Search CrossRef]
- Makowski R. A., Zimroz R. New techniques of local damage detection in machinery based on stochastic modeling using adaptive Schur filter. Applied Acoustics, Vol. 77, 2014, p. 130-137. [Search CrossRef]
- Barszcz T., Jabłoński A. Analysis of Kurtogram performance in case of high level non-Gaussian noise. Proc. of ICSV, Kraków, Poland, 2009. [Search CrossRef]
- Barszcz T., Zimroz R., Urbanek J., Jabłoński A., Bartelmus W. Bearings fault detection in gas compressor in presence of high level of non-Gaussian impulsive noise. Key Engineering Materials, Vol. 569-570, 2013, p. 473-480. [Search CrossRef]
- Antoni J., Bonnardot F., Raad A., El Badaoui M. Cyclostationary modeling of rotating machine vibration signals. Mechanical Systems and Signal Processing, Vol. 18, Issue 6, 2004, p. 1285-1314. [Search CrossRef]
- Zimroz R., Bartelmus W. Gearbox condition estimation using cyclo-stationary properties of vibration signal. Key Engineering Materials, Vol. 413-414, 2009, p. 471-478. [Search CrossRef]
- Antoni J. Cyclostationarity by examples. Mechanical Systems and Signal Processing, Vol. 23, Issue 4, 2009, p. 987-1036. [Search CrossRef]
- Randall R., Antoni J., Chobsaard S. The relationship between spectral correlation and envelope analysis for cyclostationary machine signals. Application to Ball Bearing Diagnostics. Mechanical Systems and System Processing, Vol. 15, Issue 5, 2001, p. 945-962. [Search CrossRef]
- Hurd H. L, Miamee A. Periodically Correlated Random Sequences. Spectral Theory and Practice, Wiley & Sons, New Jersey, 2007. [Search CrossRef]
- Gladyshev E. G. Periodically correlated random sequences. Journal of Soviet Mathematics, Vol. 2, 1961, p. 385-388. [Search CrossRef]
- Vecchia A. V. Periodic autoregressive-moving average (PARMA) modeling with applications to water resources. Water Resour Bull, Vol. 21, 1965, p. 21-730. [Search CrossRef]
- Bloomfield P., Hurd H. L., Lund R. B. Periodic correlation in stratospheric ozone data. Journal of Time Series Analysis, Vol. 12, 1994, p. 127-150. [Search CrossRef]
- Broszkiewicz-Suwaj E., Makagon A., Weron R., Wylomanska A. On detecting and modeling periodic correlation in financial data. Physica A, Vol. 336, 2004, p. 196-205. [Search CrossRef]
- Wyłomańska A., Obuchowski J., Zimroz R., Hurd H. Periodic autoregressive modeling of vibration time series from planetary gearbox used in bucket wheel excavator, in Cyclostationarity: Theory and Methods Lecture Notes in Mechanical Engineering, 2014, p. 171-186. [Search CrossRef]