# Fault Detection Enhancement in Rolling Element Bearings via Peak-Based Multiscale Decomposition and Envelope Demodulation

Authors: Hua-Qing Wang; Wei Hou; Gang Tang; Hong-Fang Yuan; Qing-Liang Zhao; Xi Cao
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Vibration signals of rolling element bearings faults are usually immersed in background noise, which makes it difficult to detect the faults. Wavelet-based methods being used commonly can reduce some types of noise, but there is still plenty of room for improvement due to the insufficient sparseness of vibration signals in wavelet domain. In this work, in order to eliminate noise and enhance the weak fault detection, a new kind of peak-based approach combined with multiscale decomposition and envelope demodulation is developed. First, to preserve effective middle-low frequency signals while making high frequency noise more significant, a peak-based piecewise recombination is utilized to convert middle frequency components into low frequency ones. The newly generated signal becomes so smoother that it will have a sparser representation in wavelet domain. Then a noise threshold is applied after wavelet multiscale decomposition, followed by inverse wavelet transform and backward peak-based piecewise transform. Finally, the amplitude of fault characteristic frequency is enhanced by means of envelope demodulation. The effectiveness of the proposed method is validated by rolling bearings faults experiments. Compared with traditional wavelet-based analysis, experimental results show that fault features can be enhanced significantly and detected easily by the proposed method.


---

## Body

## 1. Introduction

A rolling bearing is one of the most widely used elements in rotating machinery. As a critical component, it carries most of the load during the running of rotating machinery. If rolling bearing fails, serious problems arise, which will in turn result in the decrease of production efficiency and large economic loss. Records show that faulty bearings contribute to about thirty percent of failures in rotating machinery. Thus, it is of great importance to study the effective fault diagnosis approaches for rolling bearings.Various methods have been developed for bearing fault diagnosis and condition monitoring, such as vibration monitoring, temperature monitoring, chemical analysis, acoustic emission monitoring, sound pressure monitoring, and laser monitoring [1]. Vibration signal analysis is one of the most efficient methods thanks to the useful information of machine’s work status carried by vibration signals [2–4]. To extract fault information, for example, partial defects of bearings, vibration signals are usually processed in time domain, frequency domain, or both [5, 6]. In time domain, some statistical parameters, for example, probability density, kurtosis, root mean square, or skewness, have been introduced into bearing defect detection [7]. However, these fault indicators are not effective in all cases, especially for weak fault signals. Therefore, analysis approaches in frequency domain are developed to detect bearing defect fault. Based on the efficiency of modern fast Fourier transform (FFT), spectral analysis of vibration signals is widely used to extract characteristic defect frequencies. In fact, bearing fault signal is amplitude modulated at its characteristic defect frequencies, so a preprocessing of demodulation should be performed before FFT applied. Traditionally, envelope demodulation is used to detect the fault frequency [5, 8]. For obvious fault features, traditional envelope demodulation has notable advantages in fault type recognition. However, it does not work well in many cases for the detection of weak fault features. To solve this problem, various approaches have been developed to denoise vibration signals and to enhance the characteristic defect frequency. Cyclostationary analysis [9], for example, spectral autocorrelation analysis [10], is such a way. In recent years, some methods based on time-frequency analysis, such as empirical mode decomposition (EMD) [11–13], local mean decomposition (LMD) [14, 15], and wavelet transform [16, 17], are also developed to enhance fault features and detect failures. Both EMD and LMD are sensitive to noise level, but they are so hard to implement properly since these algorithms contain imperfections. In addition to the mentioned signal processing approaches, some other strategies are also introduced for fault features extraction as well as pattern recognition, for example, symptom parameter waves [18] and fuzzy diagnosis method [19].As mentioned above, signal’s decomposition or sparse representation is the basis of those methods. However, because of the existing high frequency components, vibration signals are not always sparse enough after decomposition. As a popular strategy for bearing vibration signals’ processing, wavelet-based analysis is effective for feature extraction from smooth signals, signal denoising, and fault feature enhancement [20–22]. But it cannot always represent vibration signals sparsely enough, especially for those mixed by plenty of interferential signals generated by other related rolling elements. This will be illustrated further with a typical example in the next section. In related research areas, a so-called peak transform based on piecewise curve recombination has been developed. The idea is to adjust piecewise curves to improve the smoothness of the target signals or images, and it has been used in 2D image representation and coding [23, 24]. To our knowledge, there is no report in the literature by far on its applications to bearing vibration signal analysis.In order to overcome those disadvantages of wavelet-based methods for weak signals processing of bearing faults, this paper developed a new multiscale decomposition strategy with sparsity-promoting by recombination of piecewise signals under different frequencies. The new strategy improves the sparsity of original vibration signals and is of great significance in fault signals denoising and enhancement. Then the wavelet coefficients of vibration signals will be sparse enough, and a threshold method based on wavelet decomposition can eliminate most noise. Finally, the FFT-based Hilbert transform is conducted for signal demodulation and extracting fault signature of characteristic defect frequency. Fault features will be significantly enhanced and easily detected through the proposed method.The paper is organized as follows. Section2 illustrates the sparseness of vibration signals in wavelet domain, which is also the motivation of this work. The proposed method of peak-based multiscale decomposition and envelope demodulation is introduced in Section 3. Some experimental results to show the effectiveness of the method are presented in Section 4. Finally, discussions and conclusions are presented in Section 5.


## 2. Vibration Signal of Rolling Bearing Faults and Its Sparsity in Wavelet Domain

Vibration signal is always acquired by sensors related to acceleration, velocity, or displacement measurement. A series of signal processing steps are then applied to make the fault features easy to be distinguished. Basically, the key of vibration signal analysis is to find out appropriate parameters describing the running conditions of bearings clearly. The commonly used characteristic defect frequency is one of the most effective indicators for rolling bearing’s fault. If local defects exist in a bearing, the measured vibration signal will be amplitude modulated. Its carrier wave is a kind of the bearing’s inherent vibration with high frequency, while its modulating wave is a pass vibration corresponding to the local defects, whose frequency is called ball pass frequency or characteristic defect frequency. The latter is usually utilized to determine whether a fault exists, and it is also adopted by this paper.In order to reduce the potential impact of signal attenuation through complex paths, sensors are often placed as close to the bearings as possible, for example, bearing housing. However, transmission path still exists between the impact point and the vibration measurement point, which would make signals more complex and hard to be demodulated. If defects exist, a fault in one surface of bearings would strike another; then a force impulse will be generated and flowing resonances in the bearing and machine will also be excited [25]. A series of impulse responses will be generated and last continuously, which may be amplitude modulated caused by the fault passing through the load zone or the mentioned varying transmission paths. Therefore, the ideal bearing defect signal is a kind of periodical impulse-like signal with high frequency nonsmooth components. A typical vibration signal of a bearing is shown in Figure 1; the vibration waveform is obviously impacted by an existing point defect on the outer race. In practice, in addition to bearing condition information, the measured vibration signals are often mixed with lots of interfering information from other components. Meanwhile, various types of noise also make the signal more complicated. In addition to the noise introduced by measurement instruments, interfering signals can also be generated accompanied with the existence of transmission path between local impact point and vibration measurement point. In a word, bearing vibration signal is nonstationary with wide frequency band, which makes the bearing fault signature extremely weak to be detected.Figure 1
Vibration waveform of a bearing with a point defect on the outer race.Existing wavelet-based denoising methods mainly rely on the sparseness of wavelet coefficients [17]. If a signal is not sparse enough in wavelet domain, it becomes challenging for feature extraction by wavelet analysis due to the difficulties to distinguish required features from large number of wavelet coefficients. Common vibration signals from faulty bearings are often composed of impulse components. A typical vibration signal is shown in Figure 1, and its wavelet coefficients are shown in Figure 2, from which we can see that most coefficients are nonzero, that is, not sparse enough. Experiments show that this is applicable to fault vibration signals of inner race, outer race, and ball element. Thus, it is important to find out a proper way to improve the sparseness, as well as reducing noise and enhancing the weak fault features.Figure 2
Wavelet coefficients of the bearing fault signal in Figure1.


## 3. The Proposed Method

### 3.1. The Basic Principle of Peak-Based Piecewise Recombination

Peak-based piecewise recombination (PPR) or so-called peak transform is a kind of nonlinear geometric transform to make the signal smoother through transforming high frequency subsections to low frequency ones. Here, a brief description of PPR is given as follows; for more information please refer to [23].For a continuous functionf(x) that has N-point peaks (or break points), xi (i=1⋯N) is defined over [a,b]. In fact, f(x) can be viewed as a piecewise function with N+1 combined curves, where the N peaks’ positions are the connection points. Then the ith piecewise curve can be called subfunctionfi(x) over interval [a,b]. The N-point forward transform of peak-based piecewise recombination can be defined as
(1)PF[f(x)]=go(x)⊕ge(x),
where
(2)go(x)=f1(x)⊕f3(x)⊕⋯⊕f2⌊(N-1)/2⌋+1(x),ge(x)=f2(x)⊕f4(x)⊕⋯⊕f2⌊N/2⌋(x)
are the cascades of all odd- and even-numbered curve segments, respectively. Here, 2⌊(N-1)/2⌋+1 and 2⌊N/2⌋ are, respectively, the largest odd and even integers that are less than or equal to N.Geometrically speaking, as illustrated in Figure3, all piecewise curves with similar positive slopes are grouped together and cascaded to be a new curve, while all negative ones are grouped into another new curve; finally, these two new subcurves are recombined to be a new curve. It is apparent that the recombination only changes the order of the piecewise curves and is reversible. The backward transform can be done by simply recascading them according to their original orders, which is symbolized as PF-1[f(x)] according to (1).Figure 3
Schematic diagram of five-point peak-based piecewise recombination.It should be noted that, although the forward and backward transforms discussed above are defined for curve functions, they can also be defined in similar manners for discrete-time signals of rolling bearings’ vibration.


### 3.2. Wavelet Decomposition

Waveletψ(a,b)(t) is obtained by translation and dilation based on a defined single function ψ(t) as
(3)ψ(a,b)(t)=1aψ(t-ba),
where a>0 is the so-called scaling parameter and b∈R is the parameter denoting time localization, which can be continuous or discrete. ψ(t) is called “mother wavelet” to generate wavelets ψ(a,b)(t).Given a signalx(t) with finite energy, the wavelet transform with analytic wavelet ψ(t) can be viewed as the convolution of x(t) with a scaled and conjugated wavelet
(4)W(a,b)=1a∫-∞∞x(t)ψ*(t-ba)dt,
where ψ*(t) stands for the complex conjugation of ψ(t).The wavelet transformW(a,b) can be considered as a function of translation b at each scale a. From (4), we can see that wavelet transform is a kind of time-frequency analysis or a time-scaled analysis. Different from analysis either in time domain or in frequency domain, or even different from that in time-frequency domain but with fixed-length windows by short time Fourier transform (SFFT), vibration signals can be decomposed by wavelet transform with multiscale analysis through dilation and translation, so that the time-frequency features of vibration signals can be more effectively extracted.Wavelet transform is also reversible as follows:(5)x(t)=Cψ-1∬W(a,b)ψ(a,b)(t)daa2db,
where
(6)Cψ=∫-∞∞|ψ(ω)^|2|ω|dω<∞,ψ^(w)=∫ψ(t)exp⁡(-jwt)dt
which provides the possibility to reconstruct the original vibration signals.


### 3.3. Envelope Demodulation

As mentioned in the previous sections, if bearing defects exist, the measured vibration signal would be amplitude modulated at its characteristic defect frequency. The modulating wave is a pass vibration signal corresponding to local defects, and various demodulation techniques have been developed for the separation. In this work, an envelope demodulation method based on Hilbert transform is employed.For a continuous time signalx(t), the Hilbert transform x(t)^ is defined as
(7)x(t)^=1π∫-∞+∞x(τ)(t-τ)dτ.
Then, combine x(t) and x(t)^ to form a new analytic signal
(8)g(t)=x(t)+jx(t)^.
The envelope of x(t) is defined as the amplitude of g(t):
(9)A(t)=[x(t)]2+[x(t)]2.
After Hilbert transform, envelope spectra could be obtained by FFT to extract characteristic defect frequencies.


### 3.4. Enhancement Frame with Peak-Based Multiscale Decomposition and Envelope Demodulation

Wavelet transform is a kind of time-frequency analysis method for signal denoising, by adjusting the discrete detail coefficients and approximation coefficients obtained from multilevel decomposition. This scheme is called threshold processing, one of whose premises is the sparse representation of signals in wavelet domain. However, different from that of perfect bearing, the vibration signal energy of bearings with defects increases at high frequencies. Thus, vibration signals are not always sparse in wavelet domain. Piecewise curve recombination can convert energies at high frequencies to ones at low frequencies. In this way, signal energies are focused in low frequency subbands, so that the piecewise recombination can improve signal’s smoothness as well as its sparseness in wavelet domain.The present paper proposed an enhancement frame for vibration signals by the combination of peak-based multiscale decomposition and envelope demodulation. The procedure is described as follows.(1)
Peak-based piecewise recombination of vibration signal: letx(t) denote the measured vibration signal, firstly; then, transform the original signal x(t) to smooth signal PF[f(x)] through piecewise recombination according to Equation (1). The key is to determine the peaks and to recombine different segments. In the present work, a cascading strategy related to local minimum and maximum is employed. The endpoint of a segment is chosen as a peak. If two adjacent peaks are equal, one of them is defined as the peak.(2)
Wavelet decomposition of the new recombined signal: choose a wavelet basis (e.g., Morlet wavelet or Db4 wavelet); then decompose the new signalPF[f(x)] at level N:
(10)W(a,b)=1a∫-∞∞PF[x(t)]ψ*(t-ba)dt.(3)
Coefficients threshold: set a threshold value and apply it to the approximation coefficients at levelN and the detail coefficients at levels 1 to N. The hard threshold can be expressed as
(11)ys={y|y|>t0|y|<t,
where t is the selected threshold, y is the original wavelet coefficients, and ys is the wavelet coefficients after thresholding. In this work, the value which covers about 95% of the smallest wavelet coefficients is set as the threshold value. Then 5% of the largest wavelet coefficients are kept as nonzero ones to reserve the signal feature.(4)
Signal recovery: reconstruct the vibration signal based on the modified approximation coefficients at levelN and the modified detail coefficients at levels 1 to N:
(12)PF′[x(t)]=Cψ-1∬W′(a,b)ψ(a,b)(t)daa2db.(5)
Piecewise signal backward transform: put the processed signal segments back to their original positions; then a new denoised vibration signalx′(t) is generated.(6)
Envelope demodulation: according to (9), use Hilbert transform to obtain the envelope A(t) of signal x′(t).(7)
Extraction of characteristic defect frequency: get the envelope spectra via Fourier analysis and distinguish different fault types of rolling bearings.In summary, these steps are shown as in Figure4.Figure 4
Flowchart of the proposed enhancement scheme for weak fault signal.


## 3.1. The Basic Principle of Peak-Based Piecewise Recombination

Peak-based piecewise recombination (PPR) or so-called peak transform is a kind of nonlinear geometric transform to make the signal smoother through transforming high frequency subsections to low frequency ones. Here, a brief description of PPR is given as follows; for more information please refer to [23].For a continuous functionf(x) that has N-point peaks (or break points), xi (i=1⋯N) is defined over [a,b]. In fact, f(x) can be viewed as a piecewise function with N+1 combined curves, where the N peaks’ positions are the connection points. Then the ith piecewise curve can be called subfunctionfi(x) over interval [a,b]. The N-point forward transform of peak-based piecewise recombination can be defined as
(1)PF[f(x)]=go(x)⊕ge(x),
where
(2)go(x)=f1(x)⊕f3(x)⊕⋯⊕f2⌊(N-1)/2⌋+1(x),ge(x)=f2(x)⊕f4(x)⊕⋯⊕f2⌊N/2⌋(x)
are the cascades of all odd- and even-numbered curve segments, respectively. Here, 2⌊(N-1)/2⌋+1 and 2⌊N/2⌋ are, respectively, the largest odd and even integers that are less than or equal to N.Geometrically speaking, as illustrated in Figure3, all piecewise curves with similar positive slopes are grouped together and cascaded to be a new curve, while all negative ones are grouped into another new curve; finally, these two new subcurves are recombined to be a new curve. It is apparent that the recombination only changes the order of the piecewise curves and is reversible. The backward transform can be done by simply recascading them according to their original orders, which is symbolized as PF-1[f(x)] according to (1).Figure 3
Schematic diagram of five-point peak-based piecewise recombination.It should be noted that, although the forward and backward transforms discussed above are defined for curve functions, they can also be defined in similar manners for discrete-time signals of rolling bearings’ vibration.


## 3.2. Wavelet Decomposition

Waveletψ(a,b)(t) is obtained by translation and dilation based on a defined single function ψ(t) as
(3)ψ(a,b)(t)=1aψ(t-ba),
where a>0 is the so-called scaling parameter and b∈R is the parameter denoting time localization, which can be continuous or discrete. ψ(t) is called “mother wavelet” to generate wavelets ψ(a,b)(t).Given a signalx(t) with finite energy, the wavelet transform with analytic wavelet ψ(t) can be viewed as the convolution of x(t) with a scaled and conjugated wavelet
(4)W(a,b)=1a∫-∞∞x(t)ψ*(t-ba)dt,
where ψ*(t) stands for the complex conjugation of ψ(t).The wavelet transformW(a,b) can be considered as a function of translation b at each scale a. From (4), we can see that wavelet transform is a kind of time-frequency analysis or a time-scaled analysis. Different from analysis either in time domain or in frequency domain, or even different from that in time-frequency domain but with fixed-length windows by short time Fourier transform (SFFT), vibration signals can be decomposed by wavelet transform with multiscale analysis through dilation and translation, so that the time-frequency features of vibration signals can be more effectively extracted.Wavelet transform is also reversible as follows:(5)x(t)=Cψ-1∬W(a,b)ψ(a,b)(t)daa2db,
where
(6)Cψ=∫-∞∞|ψ(ω)^|2|ω|dω<∞,ψ^(w)=∫ψ(t)exp⁡(-jwt)dt
which provides the possibility to reconstruct the original vibration signals.


## 3.3. Envelope Demodulation

As mentioned in the previous sections, if bearing defects exist, the measured vibration signal would be amplitude modulated at its characteristic defect frequency. The modulating wave is a pass vibration signal corresponding to local defects, and various demodulation techniques have been developed for the separation. In this work, an envelope demodulation method based on Hilbert transform is employed.For a continuous time signalx(t), the Hilbert transform x(t)^ is defined as
(7)x(t)^=1π∫-∞+∞x(τ)(t-τ)dτ.
Then, combine x(t) and x(t)^ to form a new analytic signal
(8)g(t)=x(t)+jx(t)^.
The envelope of x(t) is defined as the amplitude of g(t):
(9)A(t)=[x(t)]2+[x(t)]2.
After Hilbert transform, envelope spectra could be obtained by FFT to extract characteristic defect frequencies.


## 3.4. Enhancement Frame with Peak-Based Multiscale Decomposition and Envelope Demodulation

Wavelet transform is a kind of time-frequency analysis method for signal denoising, by adjusting the discrete detail coefficients and approximation coefficients obtained from multilevel decomposition. This scheme is called threshold processing, one of whose premises is the sparse representation of signals in wavelet domain. However, different from that of perfect bearing, the vibration signal energy of bearings with defects increases at high frequencies. Thus, vibration signals are not always sparse in wavelet domain. Piecewise curve recombination can convert energies at high frequencies to ones at low frequencies. In this way, signal energies are focused in low frequency subbands, so that the piecewise recombination can improve signal’s smoothness as well as its sparseness in wavelet domain.The present paper proposed an enhancement frame for vibration signals by the combination of peak-based multiscale decomposition and envelope demodulation. The procedure is described as follows.(1)
Peak-based piecewise recombination of vibration signal: letx(t) denote the measured vibration signal, firstly; then, transform the original signal x(t) to smooth signal PF[f(x)] through piecewise recombination according to Equation (1). The key is to determine the peaks and to recombine different segments. In the present work, a cascading strategy related to local minimum and maximum is employed. The endpoint of a segment is chosen as a peak. If two adjacent peaks are equal, one of them is defined as the peak.(2)
Wavelet decomposition of the new recombined signal: choose a wavelet basis (e.g., Morlet wavelet or Db4 wavelet); then decompose the new signalPF[f(x)] at level N:
(10)W(a,b)=1a∫-∞∞PF[x(t)]ψ*(t-ba)dt.(3)
Coefficients threshold: set a threshold value and apply it to the approximation coefficients at levelN and the detail coefficients at levels 1 to N. The hard threshold can be expressed as
(11)ys={y|y|>t0|y|<t,
where t is the selected threshold, y is the original wavelet coefficients, and ys is the wavelet coefficients after thresholding. In this work, the value which covers about 95% of the smallest wavelet coefficients is set as the threshold value. Then 5% of the largest wavelet coefficients are kept as nonzero ones to reserve the signal feature.(4)
Signal recovery: reconstruct the vibration signal based on the modified approximation coefficients at levelN and the modified detail coefficients at levels 1 to N:
(12)PF′[x(t)]=Cψ-1∬W′(a,b)ψ(a,b)(t)daa2db.(5)
Piecewise signal backward transform: put the processed signal segments back to their original positions; then a new denoised vibration signalx′(t) is generated.(6)
Envelope demodulation: according to (9), use Hilbert transform to obtain the envelope A(t) of signal x′(t).(7)
Extraction of characteristic defect frequency: get the envelope spectra via Fourier analysis and distinguish different fault types of rolling bearings.In summary, these steps are shown as in Figure4.Figure 4
Flowchart of the proposed enhancement scheme for weak fault signal.


## 4. Experimental Results and Discussions

### 4.1. Experiment Setup

To verify the effectiveness of the proposed method, experiments are designed for fault rigs of roller element bearing, as shown in Figure5. The rig is composed of a motor, a coupling, a rotor, and a shaft with two roller bearings. Basically, the sensors are located at the position near the bearings to mitigate the effects of signal attenuation. Usually, bearing housing is one of the best locations for bearing arrangement. Therefore, here, an accelerometer is located on the bearing housing to acquire the vibration signals.Figure 5
Fault rig of a roller element bearing.In the experiments, a single point defect is introduced in inner raceway, outer raceway, and ball element of different bearings, using electron-discharge machining with fault widths of 3 mm, 7 mm, and 7 mm and depths of 5 mm, 25 mm, and 25 mm, respectively.Vibration signals are measured by an accelerometer, located at the top of the bearing house (CH1, as shown in Figure5). In all experiments, the sample frequency is 100 kHz and the shaft speed is finite, 1300 rpm. Four common conditions are studied in the present paper, including a perfect healthy bearing and three other bearings with point defect on the outer race, on the inner race bearing, and on the ball element.In theory, elements of roller bearings have their own specific rotational frequencies, which may appear in envelope spectra if defects exist. These defect frequencies are called ball pass frequencies. For a bearing with a stationary outer race, these frequencies can be given by the following formulas:(13)outerracedefectfrequency,wod=Zws2d(1-dDcos⁡α),innerracedefectfrequency,wid=Zws2(1+dDcos⁡α),rollingelementdefectfrequency,wre=wsD2d×(1-d2D2cos⁡2α),
where ws is the shaft rotation frequency in rad/s, d is the diameter of rolling element, D is the pitch diameter, Z is the number of rolling elements, and α is the contact angle.In the present experiments, outer race defect frequency, inner race defect frequency, and rolling element defect frequency are 86.32 Hz, 145.84 Hz, and 51.13 Hz, respectively.


### 4.2. Vibration Signals of Typical Bearing Faults and Their Envelop Spectra

The time-domain waveforms of bearing vibration signals obtained from accelerometers are shown in Figure6, from which we can see that, different from that of perfect bearings as shown in Figure 6(a), obvious impulsive phenomenon exists in faulty vibration signals of fault bearings as shown in Figures 6(b), 6(c), and 6(d).Vibration signals of rolling bearings: (a) a normal bearing and bearings with a point defect on (b) outer race, (c) inner race bearing, and (d) ball element.
(a)
(b)
(c)
(d)To investigate the differences of the signals and to distinguish different fault types, envelope demodulation is applied first to vibration signals of the three typical bearing faults, that is, outer race fault, inner race fault, and roller element fault. The results are shown in Figures7, 8, and 9. In Figure 7, the frequency magnitude of outer race defect is lower than the second harmonic frequency of outer race defect. In Figure 8, although the inner race defect frequency slightly sticks out, it is still possibly confused by the side-lobes with almost similar magnitudes. In Figure 9, it is almost impossible to figure out the calculated ball pass frequency of 51.13 Hz. From the figures, we can see that, simply based on envelope demodulation, it is difficult to identify the fault frequencies from those of other components, especially for weak fault signals.Figure 7
Envelope spectra of vibration signal in Figure6(b) induced by outer race fault.Figure 8
Envelope spectra of vibration signal in Figure6(c) induced by inner race fault.Figure 9
Envelope spectra of vibration signal in Figure6(d) induced by roller element fault.


### 4.3. Wavelet Coefficients Distributions of the Vibration Signals with and without Peak-Based Piecewise Recombination

Before envelope demodulation, wavelet decomposition is always applied to denoise and enhance the weak fault signals. However, vibration signals are not always sparse enough in wavelet domain, which would affect the effectiveness of the enhancement. As mentioned above, peak-based piecewise recombination can help to improve the smoothness of vibration signals and make their energies more concentrated in wavelet domain. Comparisons between wavelet decompositions with and without peak-based piecewise recombination are illustrated in this section.Taking outer race fault as an example, the time-domain vibration signal of Figure6(b) is converted to that of Figure 10 after applying piecewise recombination. The signal becomes much smoother than the original signal after piecewise recombination; that is, high frequency components have been converted to low frequency ones.Figure 10
Vibration signal induced by outer race fault in Figure6(b) after peak-based piecewise recombination.For comparisons, then wavelet transform is applied to signals in Figures6(b) and 10, respectively, where 6-level decomposition with db4 wavelet basis is employed. Both distributions of the wavelet coefficients are shown in Figure 11. Without piecewise recombination, most wavelet coefficients are nonzero, that is, not sparse enough, as depicted in Figure 11(a). However, it is improved significantly with preprocessing of peak-based piecewise recombination, as depicted in Figure 11(b). Energies of the signal in wavelet domain are concentrated in only a few low frequency intervals, and almost all detail coefficients at scales 1 to 6 are close to zero. It implies that peak-based piecewise recombination can significantly improve the sparseness of vibration signals in wavelet domain.Wavelet coefficients of the vibration signal induced by outer race fault: (a) without PPR preprocessing and (b) with PPR preprocessing.
(a)
(b)Statistically, energy distributions of wavelet coefficients obtained by the two methods are listed in Table1, from which we can see that, after employing peak-based piecewise recombination (PPR) preprocessing, energies of the signal are concentrated into only a few coefficients in wavelet domain, 99.9% in low frequency bands, compared with 49.31% without peak-based piecewise recombination.Table 1
Energy distributions of the outer race fault signal’s wavelet coefficients.
Frequency band
Low frequency
High frequency
Energies without PPR
49.31%
50.69%
Energies with PPR
99.9%
0.1%The experimental results of inner race fault and roller element fault are similar to those of the outer race fault, which are shown in Figures12 and 13.Wavelet coefficients of the vibration signal induced by inner race fault: (a) without PPR preprocessing and (b) with PPR preprocessing.
(a)
(b)Wavelet coefficients of the vibration signal induced by roller element fault: (a) without PPR preprocessing and (b) with PPR preprocessing.
(a)
(b)


### 4.4. Comparisons between Traditional Wavelet-Based Detection Method and the Proposed Strategy with Envelop Demodulation

Envelop demodulation is always combined with wavelet-based denoising and enhancement methods to analyze vibration signals. Comparisons of the results after envelop demodulation based on the results of Section4.3 are presented in this section.Envelop demodulation is applied to the signals reconstructed from the thresholding of wavelet coefficients, and a backward peak-based piecewise recombination is implemented for the proposed strategy. The demodulated envelope spectra are shown in Figures14–16. In Figure 14(a), although the fault frequency can be figured out, the magnitude of the outer race defect frequency is not the most remarkable one among all frequency components; thus, it is not always certain that a bearing defect exists in outer race. With the proposed method presented in this paper, the magnitude of outer race defect frequency becomes more prominent, where the magnitudes of outer race defect frequency have been enhanced and are more easily detected, as Figure 14(b) depicted. The proposed strategy is also more effective for inner race and roller bearing fault detection of rolling bearings, and the results are shown as in Figures 15 and 16.Envelope spectra of the vibration signal induced by outer race fault: (a) wavelet analysis method and (b) the proposed method.
(a)
(b)Envelope spectra of the vibration signal induced by inner race fault: (a) wavelet analysis method and (b) the proposed method.
(a)
(b)Envelope spectra of the vibration signal induced by roller element fault: (a) wavelet analysis method and (b) the proposed method.
(a)
(b)If defect exists in the bearing, the magnitudes of faults frequencies will increase in the envelope spectrum. Therefore, they are often used to describe the fault severity and also employed to judge whether the fault features are enhanced. In this work, magnitudes of faults frequencies are chosen as the major indictors to reflect the fault features. As recorded and listed in Table2, with the proposed method, magnitudes of faults frequencies are all enhanced and they become more easily distinguished from interfered noise for the experimental faults.Table 2
Major indictors for bearing fault features.
Defect frequency magnitude of
Traditional wavelet-based method
The proposed method
outer race
149
995.2
inner race
177.7
252.5
roller element
20.26
279.2


## 4.1. Experiment Setup

To verify the effectiveness of the proposed method, experiments are designed for fault rigs of roller element bearing, as shown in Figure5. The rig is composed of a motor, a coupling, a rotor, and a shaft with two roller bearings. Basically, the sensors are located at the position near the bearings to mitigate the effects of signal attenuation. Usually, bearing housing is one of the best locations for bearing arrangement. Therefore, here, an accelerometer is located on the bearing housing to acquire the vibration signals.Figure 5
Fault rig of a roller element bearing.In the experiments, a single point defect is introduced in inner raceway, outer raceway, and ball element of different bearings, using electron-discharge machining with fault widths of 3 mm, 7 mm, and 7 mm and depths of 5 mm, 25 mm, and 25 mm, respectively.Vibration signals are measured by an accelerometer, located at the top of the bearing house (CH1, as shown in Figure5). In all experiments, the sample frequency is 100 kHz and the shaft speed is finite, 1300 rpm. Four common conditions are studied in the present paper, including a perfect healthy bearing and three other bearings with point defect on the outer race, on the inner race bearing, and on the ball element.In theory, elements of roller bearings have their own specific rotational frequencies, which may appear in envelope spectra if defects exist. These defect frequencies are called ball pass frequencies. For a bearing with a stationary outer race, these frequencies can be given by the following formulas:(13)outerracedefectfrequency,wod=Zws2d(1-dDcos⁡α),innerracedefectfrequency,wid=Zws2(1+dDcos⁡α),rollingelementdefectfrequency,wre=wsD2d×(1-d2D2cos⁡2α),
where ws is the shaft rotation frequency in rad/s, d is the diameter of rolling element, D is the pitch diameter, Z is the number of rolling elements, and α is the contact angle.In the present experiments, outer race defect frequency, inner race defect frequency, and rolling element defect frequency are 86.32 Hz, 145.84 Hz, and 51.13 Hz, respectively.


## 4.2. Vibration Signals of Typical Bearing Faults and Their Envelop Spectra

The time-domain waveforms of bearing vibration signals obtained from accelerometers are shown in Figure6, from which we can see that, different from that of perfect bearings as shown in Figure 6(a), obvious impulsive phenomenon exists in faulty vibration signals of fault bearings as shown in Figures 6(b), 6(c), and 6(d).Vibration signals of rolling bearings: (a) a normal bearing and bearings with a point defect on (b) outer race, (c) inner race bearing, and (d) ball element.
(a)
(b)
(c)
(d)To investigate the differences of the signals and to distinguish different fault types, envelope demodulation is applied first to vibration signals of the three typical bearing faults, that is, outer race fault, inner race fault, and roller element fault. The results are shown in Figures7, 8, and 9. In Figure 7, the frequency magnitude of outer race defect is lower than the second harmonic frequency of outer race defect. In Figure 8, although the inner race defect frequency slightly sticks out, it is still possibly confused by the side-lobes with almost similar magnitudes. In Figure 9, it is almost impossible to figure out the calculated ball pass frequency of 51.13 Hz. From the figures, we can see that, simply based on envelope demodulation, it is difficult to identify the fault frequencies from those of other components, especially for weak fault signals.Figure 7
Envelope spectra of vibration signal in Figure6(b) induced by outer race fault.Figure 8
Envelope spectra of vibration signal in Figure6(c) induced by inner race fault.Figure 9
Envelope spectra of vibration signal in Figure6(d) induced by roller element fault.


## 4.3. Wavelet Coefficients Distributions of the Vibration Signals with and without Peak-Based Piecewise Recombination

Before envelope demodulation, wavelet decomposition is always applied to denoise and enhance the weak fault signals. However, vibration signals are not always sparse enough in wavelet domain, which would affect the effectiveness of the enhancement. As mentioned above, peak-based piecewise recombination can help to improve the smoothness of vibration signals and make their energies more concentrated in wavelet domain. Comparisons between wavelet decompositions with and without peak-based piecewise recombination are illustrated in this section.Taking outer race fault as an example, the time-domain vibration signal of Figure6(b) is converted to that of Figure 10 after applying piecewise recombination. The signal becomes much smoother than the original signal after piecewise recombination; that is, high frequency components have been converted to low frequency ones.Figure 10
Vibration signal induced by outer race fault in Figure6(b) after peak-based piecewise recombination.For comparisons, then wavelet transform is applied to signals in Figures6(b) and 10, respectively, where 6-level decomposition with db4 wavelet basis is employed. Both distributions of the wavelet coefficients are shown in Figure 11. Without piecewise recombination, most wavelet coefficients are nonzero, that is, not sparse enough, as depicted in Figure 11(a). However, it is improved significantly with preprocessing of peak-based piecewise recombination, as depicted in Figure 11(b). Energies of the signal in wavelet domain are concentrated in only a few low frequency intervals, and almost all detail coefficients at scales 1 to 6 are close to zero. It implies that peak-based piecewise recombination can significantly improve the sparseness of vibration signals in wavelet domain.Wavelet coefficients of the vibration signal induced by outer race fault: (a) without PPR preprocessing and (b) with PPR preprocessing.
(a)
(b)Statistically, energy distributions of wavelet coefficients obtained by the two methods are listed in Table1, from which we can see that, after employing peak-based piecewise recombination (PPR) preprocessing, energies of the signal are concentrated into only a few coefficients in wavelet domain, 99.9% in low frequency bands, compared with 49.31% without peak-based piecewise recombination.Table 1
Energy distributions of the outer race fault signal’s wavelet coefficients.
Frequency band
Low frequency
High frequency
Energies without PPR
49.31%
50.69%
Energies with PPR
99.9%
0.1%The experimental results of inner race fault and roller element fault are similar to those of the outer race fault, which are shown in Figures12 and 13.Wavelet coefficients of the vibration signal induced by inner race fault: (a) without PPR preprocessing and (b) with PPR preprocessing.
(a)
(b)Wavelet coefficients of the vibration signal induced by roller element fault: (a) without PPR preprocessing and (b) with PPR preprocessing.
(a)
(b)


## 4.4. Comparisons between Traditional Wavelet-Based Detection Method and the Proposed Strategy with Envelop Demodulation

Envelop demodulation is always combined with wavelet-based denoising and enhancement methods to analyze vibration signals. Comparisons of the results after envelop demodulation based on the results of Section4.3 are presented in this section.Envelop demodulation is applied to the signals reconstructed from the thresholding of wavelet coefficients, and a backward peak-based piecewise recombination is implemented for the proposed strategy. The demodulated envelope spectra are shown in Figures14–16. In Figure 14(a), although the fault frequency can be figured out, the magnitude of the outer race defect frequency is not the most remarkable one among all frequency components; thus, it is not always certain that a bearing defect exists in outer race. With the proposed method presented in this paper, the magnitude of outer race defect frequency becomes more prominent, where the magnitudes of outer race defect frequency have been enhanced and are more easily detected, as Figure 14(b) depicted. The proposed strategy is also more effective for inner race and roller bearing fault detection of rolling bearings, and the results are shown as in Figures 15 and 16.Envelope spectra of the vibration signal induced by outer race fault: (a) wavelet analysis method and (b) the proposed method.
(a)
(b)Envelope spectra of the vibration signal induced by inner race fault: (a) wavelet analysis method and (b) the proposed method.
(a)
(b)Envelope spectra of the vibration signal induced by roller element fault: (a) wavelet analysis method and (b) the proposed method.
(a)
(b)If defect exists in the bearing, the magnitudes of faults frequencies will increase in the envelope spectrum. Therefore, they are often used to describe the fault severity and also employed to judge whether the fault features are enhanced. In this work, magnitudes of faults frequencies are chosen as the major indictors to reflect the fault features. As recorded and listed in Table2, with the proposed method, magnitudes of faults frequencies are all enhanced and they become more easily distinguished from interfered noise for the experimental faults.Table 2
Major indictors for bearing fault features.
Defect frequency magnitude of
Traditional wavelet-based method
The proposed method
outer race
149
995.2
inner race
177.7
252.5
roller element
20.26
279.2


## 5. Conclusions

In order to improve the current wavelet-based method for fault detection of rolling bearings, the present paper developed a new way to improve vibration signals’ smoothness as well as their sparseness in wavelet domain by employing peak-based piecewise recombination and envelope demodulation. Through the proposed method, most noise mixed in vibration signals can be eliminated and the weak fault signals of rolling bearings can be enhanced and become more easily detectable. The rationality and validity of the methods are also proved by various experiments. It shows that, compared with traditional wavelet-based method, the proposed method has many advantages in bearing fault signals enhancement and detection and is very simple and easy to implement.There is still room to improve the present work. The sensitivity of the method to noise and methods to make energy more concentrated are what will be studied in future work.



---

*Source: 329458-2014-05-27.xml*