Published: 10 October 2014

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
Corresponding Author:
Radoslaw Zimroz
Views 32
Reads 15
Downloads 1129

Abstract

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.

1. Introduction

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 [20]) are used to identify damage in bearings [1]. However, for damage detection at early stages [2], when signal to noise ratio is poor/informative signal is weak [11], 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. 1A crusher – general view (note bearings with yellow housing)

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 [21]. We say the signal x(t) is cyclostationary with period T if for each s and t the following properties are satisfied [22]:

1
mt=Ext=mt+T,Rx(t,τ)=E{x(t)x(t-τ)*}=Rx(t+T,τ).

And there are no smaller values of T> 0 for which the above conditions hold.

The periodically correlated processes posses very interesting properties. Therefore they found many practical applications like in hydrology [23], meteorology [24], economics [25] and electrical engineering [26].

Fig. 2Raw vibration signal a) whole observation and b) part of the signal – without and with load applied (top/bottom respectively)

Raw vibration signal a) whole observation and b) part of the signal – without and with load applied (top/bottom respectively)

a)

Raw vibration signal a) whole observation and b) part of the signal – without and with load applied (top/bottom respectively)

b)

Fig. 3Time-Frequency map (STFT) of vibration signal idle mode a) and b) loaded machine

Time-Frequency map (STFT) of vibration signal idle mode a) and b) loaded machine

a)

Time-Frequency map (STFT) of vibration signal idle mode a) and b) loaded machine

b)

Fig. 4Envelope and envelope spectrum of part of the signal (load was applied)

Envelope and envelope spectrum of part of the signal (load was applied)

a) Envelope

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, [27].

If the signal x(t) is considered to be second-order cyclostationary then, according to the definition of PC process, the autocorrelation function:

2
Rx(t,τ)=E{x(t)x(t-τ)*},

is periodic with respect to t and then it has a Fourier series expansion:

3
Rx(t,τ)=αiΑRxαi(τ)ej2παit,

where the Fourier coefficients Rxαi(τ) are known as the cyclic autocorrelation functions.

The Fourier transform of the cyclic autocorrelation function at a given cyclic frequency αi then defines the cyclic power spectrum (CPS):

4
Sxαi(f)=F{Rxαi(τ)}.

If one consider set of CPS’s for all αiA, then we can get Spectral Correlation Density Map (SCDM):

5
SCxα(f)=αiΑSxαi(f)δ(α-αi).

In order to get a normalized version of SCDM it was proposed to use so called Spectral Coherence Map (SCM):

6
γxα=SCxα(f)SCx0(f+α/2)SCx0(f-α/2) .

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 [16], there is a chance to „enhance“ SCD Map, as well as to extract signal of interest using advanced blind filtering algorithms [19].

Fig. 5Spectral coherence density map (note excitation at 30 Hz and its harmonics)

Spectral coherence density map (note excitation at 30 Hz and its harmonics)

5. Conclusions

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.

References

  • Randall R. B., Antoni J. Rolling element bearing diagnostics-A tutorial. Mechanical Systems and Signal Processing, Vol. 25, Issue 2, 2011, p. 485-520.
  • Łazarz B., Wojnar G., Czech P. Early fault detection of toothed gear in exploitation conditions. Maintenance and Reliability, Vol. 49, Issue 1, 2011.
  • Antoni J. Fast computation of the kurtogram for the detection of transient faults. Mechanical Systems and Signal Processing, Vol. 21, 2007, p. 108-124.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • Lin J., Zuo M. Gearbox fault diagnosis using adaptive wavelet filter. Mechanical Systems and Signal Processing, Vol. 17, Issue 6, p. 1259-1269, 2003.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • Zimroz R., Bartelmus W. Gearbox condition estimation using cyclo-stationary properties of vibration signal. Key Engineering Materials, Vol. 413-414, 2009, p. 471-478.
  • Antoni J. Cyclostationarity by examples. Mechanical Systems and Signal Processing, Vol. 23, Issue 4, 2009, p. 987-1036.
  • 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.
  • Hurd H. L, Miamee A. Periodically Correlated Random Sequences. Spectral Theory and Practice, Wiley & Sons, New Jersey, 2007.
  • Gladyshev E. G. Periodically correlated random sequences. Journal of Soviet Mathematics, Vol. 2, 1961, p. 385-388.
  • Vecchia A. V. Periodic autoregressive-moving average (PARMA) modeling with applications to water resources. Water Resour Bull, Vol. 21, 1965, p. 21-730.
  • 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.
  • 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.
  • 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.

About this article

Received
Accepted
04 August 2014
Published
10 October 2014
Keywords
rolling element bearing
damage detection
cyclostationary analysis
non-Gaussian noise