# Sparse Signal Representations of Bearing Fault Signals for Exhibiting Bearing Fault Features

Authors: Wei Peng; Dong Wang; Changqing Shen; Dongni Liu
Journal: Shock and Vibration
 (2016)

---

## Abstract

Sparse signal representations attract much attention in the community of signal processing because only a few coefficients are required to represent a signal and these coefficients make the signal understandable. For bearing faults’ diagnosis, bearing faults signals collected from transducers are often overwhelmed by strong low-frequency periodic signals and heavy noises. In this paper, a joint signal processing method is proposed to extract sparse envelope coefficients, which are the sparse signal representations of bearing fault signals. Firstly, to enhance bearing fault signals, particle swarm optimization is introduced to tune the parameters of wavelet transform and the optimal wavelet transform is used for retaining one of the resonant frequency bands. Thus, sparse wavelet coefficients are obtained. Secondly, to reduce the in-band noises existing in the sparse wavelet coefficients, an adaptive morphological analysis with an iterative local maximum detection method is developed to extract sparse envelope coefficients. Simulated and real bearing fault signals are investigated to illustrate how the sparse envelope coefficients are extracted. The results show that the sparse envelope coefficients can be used to represent bearing fault features and identify different localized bearing faults.


---

## Body

## 1. Introduction

Rolling element bearings are commonly used in machines to support rotation shafts. Their failures may cause unexpected machine breakdown and lead to huge economic loss. A rolling element bearing consists of an outer race, an inner race, several rollers, and a cage. Once a defect is formed on the surface of either the outer race or the inner race, an impact is generated by each of the rollers striking the defect surface and thus it excites the resonant frequencies of the structures between bearings and transducers [1–4]. Therefore, to extract bearing fault features, envelope analysis is one of the most effective methods. To conduct envelope analysis, two steps are needed. The first step aims to use a band-pass filter to retain one of the resonant frequency bands for enhancing the signal to noise ratio of bearing fault signals because bearing fault signals are often overwhelmed by strong low-frequency vibration components and heavy noises. The second step is extracting the envelope of the signals filtered by the band-pass filter [5–7]. Moreover, if the envelope signals can be represented by a few coefficients, namely, sparse envelope coefficients [8], bearing fault signals are more understandable and easily interpreted.In this paper, a joint signal processing method for extraction of sparse envelope coefficients is proposed. Firstly, an optimal wavelet filter is tuned by particle swarm optimization (PSO). For the use of wavelet transform, the similarity between a signal and a wavelet is the most concerned and the high similarity can result in large wavelet coefficients so as to highlight hidden transients. Because the shape of a Morlet wavelet is similar with the transients caused by localized bearing faults and the Morlet wavelet has a band-pass property, which can be used to retain one of the resonant frequency bands and enhance the signal to noise ratio of bearing fault signals, the Morlet wavelet is chosen in this paper [9–14]. To automatically tune the parameters of the Morlet wavelet, two aspects including a metric and an optimization algorithm must be determined. In the past years, some metrics, such as kurtosis, entropy, smoothness index, and sparsity measurement, have been investigated for optimization of the Morlet wavelet transform [9–14]. Their comparisons show that the sparsity measurement can generate better visual inspection performance and highlight bearing fault signatures, such as bearing fault characteristic frequency and its harmonics [15]. Therefore, the metric used in this paper is the sparsity measurement. To achieve the global optimal parameters of the Morlet wavelet transform, genetic algorithm, differential evolution, and stepwise scanning have been studied [9–15]. To explore the feasibility of an easy and simple optimization algorithm, particle swarm optimization is used in this paper to tune the sparsity measurement because the core of particle swarm optimization is based on the simple physical relationship among position, velocity, and acceleration. The use of PSO is simply introduced as follows. A number of particles move in a searching space. Then, a simple mathematic algorithm searches the best position by sharing the cognitive and social influences among all particles in the searching space [16–18]. For intelligent machine fault diagnosis, PSO was used to tune the parameters of support vector machine, artificial neural network, and proximal support vector machine, respectively [19–21]. The results show that these statistical prediction models combined with PSO have good prediction accuracies for identification of different machine faults.Secondly, even though the optimal wavelet filtering is conducted on bearing fault signals for enhancement of bearing fault signatures, in-band noises still exist because the optimal wavelet filtering cannot remove the noises existing in the retained resonant frequency band. In recent years, an attracting method called morphological analysis (MA) is widely investigated due to its simplicity and effectiveness in extracting envelope signals [22–29]. MA aims to use Minkowski addition and subtraction to intersect the morphological features of bearing fault signals with a structuring element (SE). However, if the morphological features are overwhelmed by other strong vibration components and heavy noises, MA may fail to retain the morphological features of bearing fault signals. Therefore, MA can be used to postprocess bearing fault signals, if the signal to noise ratio of bearing fault signals is not high. In this paper, an adaptive MA with an iterative local maximum detection method is developed to automatically find the optimal parameter of MA and extract sparse envelope signals so as to exhibit bearing fault features.The rest of this paper is outlined as follows. Section2 introduces the fundamental algorithms used in this paper. These algorithms include wavelet transform, particle swarm optimization, and morphological analysis. In Section 3, extraction of sparse envelope coefficients from bearing fault signals is proposed. In Section 4, simulated and real bearing faults signals are analyzed by using the proposed method. Conclusions are drawn in Section 5.


## 2. Fundamental Algorithms

### 2.1. Fundamental Theory of Wavelet Transform

Wavelet transform aims to calculate the inner product between an artificial wavelet and a signal. The mathematical formula for wavelet transform is defined as follows [30, 31]:(1)Wfu,s=∫-∞+∞ft×1sψu-ts¯dt=fu∗ψs′u,(2)ψs′t=1sψ-ts¯,where s is the scale parameter and u is the translation parameter. ∗ represents the convolution operator. k- takes the complex conjugate of the signal. According to the properties of Fourier transform, (1) is rewritten as(3)Wfu,s=fu∗ψs′u=F-1f^f×sψs′^sf,where F-1 is the inverse Fourier transform and ∧ is the Fourier transform. As explained in Introduction, the complex Morlet wavelet is chosen in this paper. Its temporal waveform and the corresponding frequency spectrum are given as follows:(4)ψt=σπe-σ2t2ej2πfct,(5)ϑf=e-π2/σ2×f-fc2.From (5), it is obvious that the complex Morlet wavelet has a band-pass property and its frequency support is constrained to the frequency band fc-σ/2,fc+σ/2. Because any wavelet must satisfy the admission condition, which means that the integration of a wavelet over time must be equal to zero, the following equation should be satisfied:(6)ϑ0=e-π2/σ2×0-fc2=0.It is not difficult to verify that iffc/σ>1.3, ϑ0≈0.


### 2.2. Fundamental Theory of Particle Swarm Optimization

PSO is a population based stochastic optimization method, which optimizes a metric by iteratively moving a number of particles in a searching space, according to some simple mathematical formulas related to the positions and velocities of all particles. Each particle represents one potential solution to the optimization problem. The movements of the particles are guided by their local best positions and the best swarm position. The basic theory of PSO is described in the following [16]. Considering the physical relationship among position, velocity, and acceleration, the following basic physic principle is listed as follows:(7)xik+1=xik+vik+12aik,where xik+1 and xik mean the k+1th and the kth positions of the ith particle, respectively. vik is the kth velocity of the ith particle. aik is the kth acceleration of the ith particle. Then, the acceleration of the ith particle is divided into a cognitive acceleration, which is proportional to the distance between the current position of the ith particle and the personal best position pik of the ith particle, and a social acceleration, which is proportional to the distance between the current position of the ith particle and the global best position g(k) of the ith particle. To make these two new parts more flexible, a cognitive coefficient c1 and a social coefficient c2 are used. Consequently, (7) is reformulated as(8)xik+1=xik+vik+12c1pik-xik+12c2gk-xik.Then, in order to prevent the velocities from getting out of control, the influence of friction is considered by introducing an inertia weightω, which is smaller than 1, to (8). Besides, the constant 1/2 is replaced by two random numbers r1i and r2i, which are limited to the values between 0 and 1. Equation (8) is revised as follows:(9)xik+1=xik+ωvik+r1ic1pik-xik+r2ikc2gk-xik.At last, (9) consists of the following velocity and position update equations:(10)vik+1=ωvik+r1ic1pik-xik+r2ikc2gk-xik,xik+1=xik+vik+1.


### 2.3. Fundamental Theory of Morphological Analysis

MA aims to extract the morphological shape of a temporal signal. It uses a structuring element to intersect with the temporal signal. Letfn be the one-dimensional signal over a domain F=0,1,2,…,N-1 and let gn be the structuring element over a domain G=0,1,2,…,M-1. The basic morphology operators include the dilation operator and the erosion operator, which are related to Minkowski addition and subtraction. The equations for the erosion operator ⊖ and the dilation operator ⊕ are defined as follows [22]:(11)f⊖gn=min⁡fn+m-gm,m∈0,1,2,…,M-1,(12)f⊕gn=max⁡fn-m+gm,m∈0,1,2,…,M-1.The erosion operator reduces the wave peaks and enlarges the signal minima. On the contrary, the dilation operator increases the wave valleys and enlarges signal maxima [24]. Other morphological operators are constructed based on the combination of the above two operators. Some popular morphological operators used in bearing fault diagnosis are introduced in the following. The opening operator and the closing operator are defined as:(13)f∘gn=f⊖gs⊕gn,f·gn=f⊕gs⊖gn,where gsn is the reflection of gn. The opening operator function is smoothing the signal from the bottom by cutting wave peaks and the closing operator function is smoothing the signal from the top by filling up its wave valleys [24]. The average operator (AVG), the difference operator (DIF), the Black Top-Hat operator (BTH), and the White Top-Hat operator (WTH) are defined as [32](14)AVGf=f·g+f∘g2,DIFf=f·g-f∘g,BTHf=f·g-f,WTHf=f-f∘g.The ability of the average operator lies in flattening both the positive and negative impulsive features. On the contrary, the difference operator is used to extract both the positive and negative impulsive features. The Black Top-Hat and the White Top-Hat operators are employed to extract the negative and positive impulsive features, respectively.


## 2.1. Fundamental Theory of Wavelet Transform

Wavelet transform aims to calculate the inner product between an artificial wavelet and a signal. The mathematical formula for wavelet transform is defined as follows [30, 31]:(1)Wfu,s=∫-∞+∞ft×1sψu-ts¯dt=fu∗ψs′u,(2)ψs′t=1sψ-ts¯,where s is the scale parameter and u is the translation parameter. ∗ represents the convolution operator. k- takes the complex conjugate of the signal. According to the properties of Fourier transform, (1) is rewritten as(3)Wfu,s=fu∗ψs′u=F-1f^f×sψs′^sf,where F-1 is the inverse Fourier transform and ∧ is the Fourier transform. As explained in Introduction, the complex Morlet wavelet is chosen in this paper. Its temporal waveform and the corresponding frequency spectrum are given as follows:(4)ψt=σπe-σ2t2ej2πfct,(5)ϑf=e-π2/σ2×f-fc2.From (5), it is obvious that the complex Morlet wavelet has a band-pass property and its frequency support is constrained to the frequency band fc-σ/2,fc+σ/2. Because any wavelet must satisfy the admission condition, which means that the integration of a wavelet over time must be equal to zero, the following equation should be satisfied:(6)ϑ0=e-π2/σ2×0-fc2=0.It is not difficult to verify that iffc/σ>1.3, ϑ0≈0.


## 2.2. Fundamental Theory of Particle Swarm Optimization

PSO is a population based stochastic optimization method, which optimizes a metric by iteratively moving a number of particles in a searching space, according to some simple mathematical formulas related to the positions and velocities of all particles. Each particle represents one potential solution to the optimization problem. The movements of the particles are guided by their local best positions and the best swarm position. The basic theory of PSO is described in the following [16]. Considering the physical relationship among position, velocity, and acceleration, the following basic physic principle is listed as follows:(7)xik+1=xik+vik+12aik,where xik+1 and xik mean the k+1th and the kth positions of the ith particle, respectively. vik is the kth velocity of the ith particle. aik is the kth acceleration of the ith particle. Then, the acceleration of the ith particle is divided into a cognitive acceleration, which is proportional to the distance between the current position of the ith particle and the personal best position pik of the ith particle, and a social acceleration, which is proportional to the distance between the current position of the ith particle and the global best position g(k) of the ith particle. To make these two new parts more flexible, a cognitive coefficient c1 and a social coefficient c2 are used. Consequently, (7) is reformulated as(8)xik+1=xik+vik+12c1pik-xik+12c2gk-xik.Then, in order to prevent the velocities from getting out of control, the influence of friction is considered by introducing an inertia weightω, which is smaller than 1, to (8). Besides, the constant 1/2 is replaced by two random numbers r1i and r2i, which are limited to the values between 0 and 1. Equation (8) is revised as follows:(9)xik+1=xik+ωvik+r1ic1pik-xik+r2ikc2gk-xik.At last, (9) consists of the following velocity and position update equations:(10)vik+1=ωvik+r1ic1pik-xik+r2ikc2gk-xik,xik+1=xik+vik+1.


## 2.3. Fundamental Theory of Morphological Analysis

MA aims to extract the morphological shape of a temporal signal. It uses a structuring element to intersect with the temporal signal. Letfn be the one-dimensional signal over a domain F=0,1,2,…,N-1 and let gn be the structuring element over a domain G=0,1,2,…,M-1. The basic morphology operators include the dilation operator and the erosion operator, which are related to Minkowski addition and subtraction. The equations for the erosion operator ⊖ and the dilation operator ⊕ are defined as follows [22]:(11)f⊖gn=min⁡fn+m-gm,m∈0,1,2,…,M-1,(12)f⊕gn=max⁡fn-m+gm,m∈0,1,2,…,M-1.The erosion operator reduces the wave peaks and enlarges the signal minima. On the contrary, the dilation operator increases the wave valleys and enlarges signal maxima [24]. Other morphological operators are constructed based on the combination of the above two operators. Some popular morphological operators used in bearing fault diagnosis are introduced in the following. The opening operator and the closing operator are defined as:(13)f∘gn=f⊖gs⊕gn,f·gn=f⊕gs⊖gn,where gsn is the reflection of gn. The opening operator function is smoothing the signal from the bottom by cutting wave peaks and the closing operator function is smoothing the signal from the top by filling up its wave valleys [24]. The average operator (AVG), the difference operator (DIF), the Black Top-Hat operator (BTH), and the White Top-Hat operator (WTH) are defined as [32](14)AVGf=f·g+f∘g2,DIFf=f·g-f∘g,BTHf=f·g-f,WTHf=f-f∘g.The ability of the average operator lies in flattening both the positive and negative impulsive features. On the contrary, the difference operator is used to extract both the positive and negative impulsive features. The Black Top-Hat and the White Top-Hat operators are employed to extract the negative and positive impulsive features, respectively.


## 3. Extraction of Sparse Envelope Coefficients for Exhibiting Bearing Fault Features

It is not difficult to find that the morphological features extracted by using MA fully depend on the shape of a temporal signal. Because of the interruption from strong low-frequency periodic components and heavy noises, the morphological features of bearing fault signals are prone to be overwhelmed. Therefore, it is necessary to enhance weak bearing fault signals prior to the use of MA. As illustrated in the previous sections, the complex Morlet wavelet optimized by PSO is used to preprocess bearing fault signals and an adaptive MA is developed to postprocess bearing fault signals and to extract sparse envelope coefficients for exhibiting bearing fault features. The flowchart of the proposed method is shown in Figure1. Each step used in Figure 1 is detailed in the following paragraphs.Figure 1
The flowchart of the developed method for extracting sparse envelope coefficients.Step 1.
Load an original bearing fault vibration signal, which is collected by using a transducer attached to a bearing housing.Step 2.
The parameters of the complex Morlet wavelet are tuned by PSO. First, the parameters of particle swarm optimization must be initialized. Each particle corresponds to two-dimensional coordinates, which are used to represent the center frequency and bandwidth of the complex Morlet wavelet, respectively. The searching range for the center frequency was set to the frequency band ranging from0.1×Fs to 0.4×Fs for avoiding the interruption caused by the low-frequency vibration components. The searching range for the bandwidth is from 3×FI to 0.2×Fs. Here, to contain enough fault signatures, the minimum bandwidth should be three times larger than the inner race fault characteristic frequency. Fs is the sampling frequency and FI is the inner race fault characteristic frequency, which will be formulated later. According to the literature review [17, 18], the inertia weight, the number of particles, and the maximum iteration were empirically set to 0.9, 24, and 50, respectively, because these parameters are sufficiently large for solving an optimization problem. For all particles, the initial positions are randomly generated within the specified searching ranges. The metric, namely, sparsity measurement, is defined as [15](15)Sik=∑n=1Lbikn2∑n=1Lbikn=bikn2bikn1,1≤i≤24,1≤k≤50,where bikn2 and bikn1 are L2 norm and L1 norm, respectively. bikn represents the signal filtered by the complex Morlet wavelet transform, the parameters of which are the two-dimensional coordinates xik of the ith particle at the kth iteration. L is the length of the signal. At the kth iteration, if the Sik of the ith particle is larger than the personal largest sparsity value Sp(best) of the ith particle (recorded before the kth iteration), the personal best position of the ith particle is updated by its current position xik. Then, if the Sik of the ith particle is larger than the global largest sparsity value Sg(best) of all particles (recorded before the kth iteration), the global best position of all particles is updated by the current position xik of the ith particle. At last, repeat the above steps until the maximum iteration is reached. The optimal parameters of the complex Morlet wavelet are automatically established by PSO with the maximum sparsity. Here, the real part of the filtered signal obtained by the optimal complex Morlet wavelet is denoted by boptrealn.Step 3.
Once the weak bearing fault signal is enhanced by the optimal complex Morlet wavelet, MA is performed to get sparse envelope coefficients. As introduced in Section2.3, there are some available morphological operators. To select a proper morphological operator, their comparisons for processing a simulated signal are shown in Figure 2. The structuring element is the flat element with a length of 30 samples because the flat element is a simple and effective element to process a one-dimensional signal [22, 24, 27, 28]. From the results shown in Figure 2, the closing and opening operators can be used to extract the positive and negative envelopes of the signal, respectively. In this paper, only the closing operator is used because the positive envelope is used for further analyses. After both morphological operator and the structuring element are determined, it is necessary to determine the length of the flat element to extract the morphological features of bearing fault signals and suppress in-band noises.Figure 2
The signals obtained by using various morphological operators: (a) the signal filtered by the erosion operator; (b) the signal filtered by the dilation operator; (c) the signal filtered by opening operator; (d) the signal filtered by the closing operator; (e) the signal filtered by the AVG operator; (f) the signal filtered by the DIF operator; (g) the signal filtered by the BTH operator; and (h) the signal filtered by the WTH operator.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)Reference [22] reported a method to empirically decide the length of the flat element and indicated that the optimal length should be 0.6 to 0.7 of the bearing fault period. However, in our previous research [26], it is found that the empirical optimal length is not always effective for extracting bearing fault features. This paper develops an adaptive morphological analysis with an iterative local maximum detection method to automatically retain sparse envelope coefficients. boptrealn is processed by the closing operator with various flat element lengths and their corresponding results are denoted by mmoptrealn,l, where l is the length of the flat element used in the closing operator. The maximum length of the flat element does not exceed the desired bearing fault period T (unit: samples), which will be defined later. First, find the local maxima lmj, j=1,2,…,N of the signal mmoptrealn,l and calculate its local maxima number N. Subtract the theoretical impulsive number Z from the local maxima number N and obtain a difference (considering the influence of the negative difference, the absolute value of the difference is used in this paper). The theoretical impulsive number can be calculated as(16)Z=roundLT,where the round(∗) function is taking the element to the nearest integer.The bearing fault periodT can be represented by the outer race fault characteristic period To, the inner race fault characteristic period TI, and rolling element fault characteristic period TB. The calculations of these periods are given by(17)To=1P×fs/21-d/Dcos⁡α′,TI=1P×fs/21+d/Dcos⁡α′,TB=1D×fs/d1-d2/D2cos⁡α′,where fs is the shaft rotation frequency in Hz. d and D are diameters of the rolling element and the pitch, respectively. P is the number of rolling elements and α is the contact angle.If the difference between the local maxima numberN and the theoretical impulsive number Z is the smallest, it means that the extracted signal by using the closing operator is the best one because the signal filtered by the complex Morlet transform has a high signal to noise ratio and then the local maxima of the filtered signal are the most possible to be the local peaks of the impulses caused by localized bearing faults. Therefore, the optimal length of the flat element can be decided. However, it should be pointed out that sometimes the amplitudes of random noises may affect the number of local maxima. Therefore, the iterative local maximum selection is developed to remove the pseudo- and abnormal locations of the local maxima.Locate all local maxima and denote them byllmj, j=1,2,…,N. Then, calculate the adjacent distances between two successive locations llmj+1-llmj. There are N-1 distances disj, j=1,2,…,N-1.Assume that all the calculated distances are subject to a normal distributionNμ,o2. Here, μ and o2 are the mean and the variance, respectively. In statistics and probability theory, standard deviation o can be used to measure the diversity of samples. It is expected that all distances disj, j=1,2,…,N-1 tend to be close to the mean and have a low standard deviation. Therefore, a high standard deviation can be employed to reject the outlier (abnormal distance). In this paper, two standard deviations 2o are used. Then, it is believed that about the 95 percent of the all distances are within two standard deviations μ±2o. Moreover, only distances that are smaller than μ-2o are rejected because all distances are always positive.Find the abnormal distances, which are denoted byllmz+1-llmz. Calculate the distances llmz+2-llmz+1 and llmz-llmz-1, respectively. If llmz+2-llmz+1>llmz-llmz-1, llmz is discarded. Otherwise, llmz+1 is discarded.Repeat the iterative local maximum selection method until all distances are within two standard deviations. Finally, sparse envelope coefficients are extracted.


## 4. Case Studies

### 4.1. A Simulated Fault Signal Analyzed by the Proposed Method

In the first case study, a simulated signal is used for analyses. The simulated signal contains the impulses with an exponential decay, two sinusoidal signals, and noises. Here, the two sinusoidal signals can be regarded as two strong interruptions caused by two low-frequency components:(18)yk=∑rexp⁡-α×k-r×Fs/fm-τrFssin2πf1×k-r×Fs/fm-τrFs+0.9sin2πf2kFs+sin2πf3kFs+nnk,where α is equal to 900, nn(k) is the noise term, and fm is the modulating frequency (equal to 100 Hz). Fs is the sampling frequency set to 12000 Hz. f1 is the carrier frequency, equal to 3500 Hz. f2 and f3 are sinusoidal frequencies, equal to 70 Hz and 140 Hz, separately. 3600 samples are used. Normally distributed heavy noises with a mean of 0 and a variance of 0.2 are used as the noises. The simulated signal, the noise signal, and the mixed signal are shown in Figures 3(a), 3(b), and 3(c), respectively. From the mixed signal shown in Figure 3(c), it is hard to identify the periodic impulsive signal (the periodic intervals are equal to 10 ms).Figure 3
The simulated temporal signals: (a) the pure simulated signal consisting of impulses and two sinusoidal signals; (b) the noise signal; and (c) the pure simulated signal mixed with noises.
(a)
(b)
(c)The proposed method is employed to analyze the mixed signal. First, the optimal parameters of the complex Morlet wavelet are automatically determined by using PSO with the maximum sparsity measurement. The optimal center frequency, the optimal bandwidth, and the frequency spectrum of the complex Morlet wavelet are shown in Figure4(b). For a comparison, the frequency spectrum of the mixed signal is plotted in Figure 4(a). It is found that the optimal complex Morlet wavelet is correct to locate the simulated resonant frequency band around 3500 Hz. Then, the mixed signal is filtered by the optimal complex Morlet wavelet and its corresponding frequency spectrum is given in Figure 4(c). The real part of the filtered signal is depicted in Figure 5(a), where it is seen that the in-band noises in the resonant frequency band still exist in the filtered signal. On the other hand, to validate the correction of the filtered signal, the power spectrum of the envelope of the filtered signal is plotted in Figure 5(b), where the fundamental frequency 100 Hz and its first two harmonics are identified. The global best values by using PSO are shown in Figure 6, where it is seen that sparsity value reaches the global optimal value quickly.Figure 4
Frequency spectra: (a) the frequency spectrum of the mixed signal; (b) the frequency spectrum of the optimal complex Morlet wavelet; and (c) the frequency spectrum of the signal filtered by the optimal complex Morlet wavelet transform.
(a)
(b)
(c)Figure 5
The signals obtained by the optimal complex Morlet filtering: (a) the real part of the filtered signal; (b) the power spectrum of the envelope of the filtered signal.
(a)
(b)Figure 6
The global best values at different iterations by using PSO.At last, the morphological analysis with the closing operator is employed to extract the cyclic bearing fault characteristics. The obtained temporal signal by morphological analysis is given in Figure7(a). The absolute difference between the theoretical impulsive number and the actual impulsive number is shown in Figure 7(b), where the optimal length of 46 is found for morphological analysis. The local maximum locations of the signal shown in Figure 7(a) are displayed in Figure 7(c). It is obvious that there is a pseudolocation. The iterative local maximum selection method is applied to process these local maximum locations shown in Figure 7(c). Finally, the revised local maximum locations are plotted in Figure 7(d). Compared with the signal shown in Figure 5(a), the result plotted in Figure 7(d) can be more understandable. Moreover, it is clear to see the signal having a cyclic interval of 10 ms by picking up the interval between two successive spikes. Besides, the final signal shown in Figure 7(d) is the sparse envelope coefficients of the simulated signal mixed with noise.Figure 7
The results obtained by morphological analysis: (a) envelope extraction of the signal shown in Figure5(a) by using the closing operator; (b) the optimal length of the SE; (c) the local maximum locations of the signal shown in Figure 7(a); and (d) the revised local maximum locations of the signal shown in Figure 7(c) by using the iterative local maximum selection.
(a)
(b)
(c)
(d)


### 4.2. Experimental Fault Signals Analyzed by the Proposed Method

In this paper, the real motor bearing data picked up with a sampling frequency of 12 k Hz by an accelerometer at the drive end of the motor housing are used to validate the proposed method. A single-point defect was, respectively, introduced to the outer race and inner race of a normal bearing using electrodischarge machining. The fault diameter is 0.007 inches and the fault depth is 0.0011 inches. The motor load was 0 HP and the motor speed was 1797 rpm [33]. Outer race and inner race fault characteristic frequencies were calculated as 107 Hz and 162 Hz, respectively.The original outer race fault signal and its frequency spectrum are shown in Figures8(a) and 8(b), respectively. The optimal parameters of the complex Morlet wavelet and its frequency spectrum are shown in Figure 8(c). Then, the optimal complex Morlet wavelet is used to retain one of the resonant frequency bands. The frequency spectrum of the signal obtained by the optimal wavelet filtering is displayed in Figure 8(d). It is clear to find that one of the resonant frequency bands is kept for further morphological analysis. The real part of the filtered signal and its envelope spectrum are shown in Figures 9(a) and 9(b), respectively, which demonstrate that the optimal filtering method correctly retains the major bearing fault signatures. The global best values by using particle swarm optimization are plotted in Figure 10.Figure 8
The signals: (a) the real outer race fault signal; (b) the frequency spectrum of the original outer race fault signal; (c) the frequency spectrum of the optimal complex Morlet wavelet; and (d) the frequency spectrum of the filtered signal.
(a)
(b)
(c)
(d)Figure 9
The signals obtained by the optimal complex Morlet filtering: (a) the real part of the filtered signal; (b) the power spectrum of envelope of the filtered signal.
(a)
(b)Figure 10
The global best values at different iterations by using PSO.Finally, morphological analysis is used to modify the shape of the signal shown in Figure9(a). The extracted envelope by using the closing operator with the optimal length of the flat structuring element is displayed in Figure 11(a). The absolute differences are plotted in Figure 11(b) to find the optimal length of the flat structuring element. Here, the length of 29 is found. The local maximum locations of the envelope shown in Figure 11(a) are given in Figure 11(c). It is obvious that the noises interrupt the local maximum locations. One pseudolocation is caused by an unexpected noise. Therefore, the iterative local maximum detection method is used to remove the irrelevant local maximum location. The result is shown in Figure 11(d). After removing the interruption caused by the unexpected noise, it is clear to see the pure cyclic fault characteristics with intervals of 9.3 ms.Figure 11
The results obtained by morphological filtering: (a) envelope extraction of the signal shown in Figure9(a) by using the closing operator; (b) the optimal length of SE; (c) the local maximum locations of the signal shown in Figure 11(a); and (d) the revised local maximum locations of the signal shown in Figure 11(c) by using the iterative local maximum selection.
(a)
(b)
(c)
(d)The same procedure is applied to process the inner race fault signal shown in Figure12(a). The frequency spectrum of the inner race fault signal is plotted in Figure 12(b). The optimal parameters of the complex Morlet wavelet are given in Figure 12(c). The frequency spectrum of the filtered signal obtained by the optimal complex Morlet wavelet is displayed in Figure 12(d). Figures 13(a) and 13(b) show the real part of the filtered signal and its corresponding envelope spectrum. Obviously, the optimal complex Morlet wavelet transform reserves the inner race fault signatures. The global best values by using PSO are depicted in Figure 14.Figure 12
The signals: (a) the real inner race fault signal; (b) the frequency spectrum of the real inner race fault signal; (c) the frequency spectrum of the optimal complex Morlet wavelet; and (d) the frequency spectrum of the filtered signal.
(a)
(b)
(c)
(d)Figure 13
The signals obtained by the optimal complex Morlet filtering: (a) the real part of the signal filtered by the optimal Morlet wavelet transform; (b) the power spectrum of the envelope of the filtered signal.
(a)
(b)Figure 14
The global best values at every iteration by using particle swarm optimization.The optimal process for the length selection of the flat structuring element is shown in Figure15(b). The envelope obtained by using the closing operator with the optimal length is given in Figure 15(a). After that, the local maximum locations are extracted and the result is given in Figure 15(c). The pseudolocation is caused by an unexpected noise. In order to remove the irrelevant local maximum location shown in Figure 15(c), the iterative local maximum detection method is used. The final resulting signal is shown in Figure 15(d), where a pure cyclic fault signal is generated.Figure 15
The results obtained by morphological filtering: (a) envelope extraction of the signal shown in Figure13(a) by using the closing operator; (b) the optimal length of the flat SE; (c) the local maximum locations of the signal shown in (a); and (d) the revised local maximum locations of the signal shown in (c).
(a)
(b)
(c)
(d)


## 4.1. A Simulated Fault Signal Analyzed by the Proposed Method

In the first case study, a simulated signal is used for analyses. The simulated signal contains the impulses with an exponential decay, two sinusoidal signals, and noises. Here, the two sinusoidal signals can be regarded as two strong interruptions caused by two low-frequency components:(18)yk=∑rexp⁡-α×k-r×Fs/fm-τrFssin2πf1×k-r×Fs/fm-τrFs+0.9sin2πf2kFs+sin2πf3kFs+nnk,where α is equal to 900, nn(k) is the noise term, and fm is the modulating frequency (equal to 100 Hz). Fs is the sampling frequency set to 12000 Hz. f1 is the carrier frequency, equal to 3500 Hz. f2 and f3 are sinusoidal frequencies, equal to 70 Hz and 140 Hz, separately. 3600 samples are used. Normally distributed heavy noises with a mean of 0 and a variance of 0.2 are used as the noises. The simulated signal, the noise signal, and the mixed signal are shown in Figures 3(a), 3(b), and 3(c), respectively. From the mixed signal shown in Figure 3(c), it is hard to identify the periodic impulsive signal (the periodic intervals are equal to 10 ms).Figure 3
The simulated temporal signals: (a) the pure simulated signal consisting of impulses and two sinusoidal signals; (b) the noise signal; and (c) the pure simulated signal mixed with noises.
(a)
(b)
(c)The proposed method is employed to analyze the mixed signal. First, the optimal parameters of the complex Morlet wavelet are automatically determined by using PSO with the maximum sparsity measurement. The optimal center frequency, the optimal bandwidth, and the frequency spectrum of the complex Morlet wavelet are shown in Figure4(b). For a comparison, the frequency spectrum of the mixed signal is plotted in Figure 4(a). It is found that the optimal complex Morlet wavelet is correct to locate the simulated resonant frequency band around 3500 Hz. Then, the mixed signal is filtered by the optimal complex Morlet wavelet and its corresponding frequency spectrum is given in Figure 4(c). The real part of the filtered signal is depicted in Figure 5(a), where it is seen that the in-band noises in the resonant frequency band still exist in the filtered signal. On the other hand, to validate the correction of the filtered signal, the power spectrum of the envelope of the filtered signal is plotted in Figure 5(b), where the fundamental frequency 100 Hz and its first two harmonics are identified. The global best values by using PSO are shown in Figure 6, where it is seen that sparsity value reaches the global optimal value quickly.Figure 4
Frequency spectra: (a) the frequency spectrum of the mixed signal; (b) the frequency spectrum of the optimal complex Morlet wavelet; and (c) the frequency spectrum of the signal filtered by the optimal complex Morlet wavelet transform.
(a)
(b)
(c)Figure 5
The signals obtained by the optimal complex Morlet filtering: (a) the real part of the filtered signal; (b) the power spectrum of the envelope of the filtered signal.
(a)
(b)Figure 6
The global best values at different iterations by using PSO.At last, the morphological analysis with the closing operator is employed to extract the cyclic bearing fault characteristics. The obtained temporal signal by morphological analysis is given in Figure7(a). The absolute difference between the theoretical impulsive number and the actual impulsive number is shown in Figure 7(b), where the optimal length of 46 is found for morphological analysis. The local maximum locations of the signal shown in Figure 7(a) are displayed in Figure 7(c). It is obvious that there is a pseudolocation. The iterative local maximum selection method is applied to process these local maximum locations shown in Figure 7(c). Finally, the revised local maximum locations are plotted in Figure 7(d). Compared with the signal shown in Figure 5(a), the result plotted in Figure 7(d) can be more understandable. Moreover, it is clear to see the signal having a cyclic interval of 10 ms by picking up the interval between two successive spikes. Besides, the final signal shown in Figure 7(d) is the sparse envelope coefficients of the simulated signal mixed with noise.Figure 7
The results obtained by morphological analysis: (a) envelope extraction of the signal shown in Figure5(a) by using the closing operator; (b) the optimal length of the SE; (c) the local maximum locations of the signal shown in Figure 7(a); and (d) the revised local maximum locations of the signal shown in Figure 7(c) by using the iterative local maximum selection.
(a)
(b)
(c)
(d)


## 4.2. Experimental Fault Signals Analyzed by the Proposed Method

In this paper, the real motor bearing data picked up with a sampling frequency of 12 k Hz by an accelerometer at the drive end of the motor housing are used to validate the proposed method. A single-point defect was, respectively, introduced to the outer race and inner race of a normal bearing using electrodischarge machining. The fault diameter is 0.007 inches and the fault depth is 0.0011 inches. The motor load was 0 HP and the motor speed was 1797 rpm [33]. Outer race and inner race fault characteristic frequencies were calculated as 107 Hz and 162 Hz, respectively.The original outer race fault signal and its frequency spectrum are shown in Figures8(a) and 8(b), respectively. The optimal parameters of the complex Morlet wavelet and its frequency spectrum are shown in Figure 8(c). Then, the optimal complex Morlet wavelet is used to retain one of the resonant frequency bands. The frequency spectrum of the signal obtained by the optimal wavelet filtering is displayed in Figure 8(d). It is clear to find that one of the resonant frequency bands is kept for further morphological analysis. The real part of the filtered signal and its envelope spectrum are shown in Figures 9(a) and 9(b), respectively, which demonstrate that the optimal filtering method correctly retains the major bearing fault signatures. The global best values by using particle swarm optimization are plotted in Figure 10.Figure 8
The signals: (a) the real outer race fault signal; (b) the frequency spectrum of the original outer race fault signal; (c) the frequency spectrum of the optimal complex Morlet wavelet; and (d) the frequency spectrum of the filtered signal.
(a)
(b)
(c)
(d)Figure 9
The signals obtained by the optimal complex Morlet filtering: (a) the real part of the filtered signal; (b) the power spectrum of envelope of the filtered signal.
(a)
(b)Figure 10
The global best values at different iterations by using PSO.Finally, morphological analysis is used to modify the shape of the signal shown in Figure9(a). The extracted envelope by using the closing operator with the optimal length of the flat structuring element is displayed in Figure 11(a). The absolute differences are plotted in Figure 11(b) to find the optimal length of the flat structuring element. Here, the length of 29 is found. The local maximum locations of the envelope shown in Figure 11(a) are given in Figure 11(c). It is obvious that the noises interrupt the local maximum locations. One pseudolocation is caused by an unexpected noise. Therefore, the iterative local maximum detection method is used to remove the irrelevant local maximum location. The result is shown in Figure 11(d). After removing the interruption caused by the unexpected noise, it is clear to see the pure cyclic fault characteristics with intervals of 9.3 ms.Figure 11
The results obtained by morphological filtering: (a) envelope extraction of the signal shown in Figure9(a) by using the closing operator; (b) the optimal length of SE; (c) the local maximum locations of the signal shown in Figure 11(a); and (d) the revised local maximum locations of the signal shown in Figure 11(c) by using the iterative local maximum selection.
(a)
(b)
(c)
(d)The same procedure is applied to process the inner race fault signal shown in Figure12(a). The frequency spectrum of the inner race fault signal is plotted in Figure 12(b). The optimal parameters of the complex Morlet wavelet are given in Figure 12(c). The frequency spectrum of the filtered signal obtained by the optimal complex Morlet wavelet is displayed in Figure 12(d). Figures 13(a) and 13(b) show the real part of the filtered signal and its corresponding envelope spectrum. Obviously, the optimal complex Morlet wavelet transform reserves the inner race fault signatures. The global best values by using PSO are depicted in Figure 14.Figure 12
The signals: (a) the real inner race fault signal; (b) the frequency spectrum of the real inner race fault signal; (c) the frequency spectrum of the optimal complex Morlet wavelet; and (d) the frequency spectrum of the filtered signal.
(a)
(b)
(c)
(d)Figure 13
The signals obtained by the optimal complex Morlet filtering: (a) the real part of the signal filtered by the optimal Morlet wavelet transform; (b) the power spectrum of the envelope of the filtered signal.
(a)
(b)Figure 14
The global best values at every iteration by using particle swarm optimization.The optimal process for the length selection of the flat structuring element is shown in Figure15(b). The envelope obtained by using the closing operator with the optimal length is given in Figure 15(a). After that, the local maximum locations are extracted and the result is given in Figure 15(c). The pseudolocation is caused by an unexpected noise. In order to remove the irrelevant local maximum location shown in Figure 15(c), the iterative local maximum detection method is used. The final resulting signal is shown in Figure 15(d), where a pure cyclic fault signal is generated.Figure 15
The results obtained by morphological filtering: (a) envelope extraction of the signal shown in Figure13(a) by using the closing operator; (b) the optimal length of the flat SE; (c) the local maximum locations of the signal shown in (a); and (d) the revised local maximum locations of the signal shown in (c).
(a)
(b)
(c)
(d)


## 5. Conclusion

This paper reported a method which was used to extract sparse envelope coefficients for exhibiting bearing fault features. The proposed method consisted of two steps. Firstly, a Morlet wavelet was optimized by particle swarm optimization and then sparse wavelet coefficients were extracted from bearing fault signals. Even though sparse wavelet coefficients were useful to represent bearing fault signatures, in-band noises still existed and could not be removed by using the solely optimal wavelet transform. In this step, the sparse wavelet coefficients were not succinct enough. Secondly, to reduce in-band noises, an adaptive morphological analysis with an iterative local maximum detection method was developed to postprocess the signal obtained by the optimal wavelet filtering. The optimal flat length used in morphological analysis with the closing operator was automatically determined to retain sparse envelope coefficients and remove the pseudo- and abnormal peaks caused by unexpected noises. After the above two steps were conducted, the sparse signal representation of bearing fault signals was obtained; only a few coefficients were used to represent bearing fault features. These sparse coefficients made bearing fault signals succinct and understandable. The case studies were conducted to illustrate how the proposed method worked and the results showed that the proposed method can be used to extract the sparse representation of bearing fault signals and identify different localized bearing faults.In the future, the following works will be conducted. Firstly, more statistical metrics will be used as objective functions so as to guide optimization of wavelet transform. Secondly, besides the popular Morlet wavelet, more wavelets that are highly similar with impulses caused by bearing defects will be investigated. Thirdly, different optimization methods will be used to find optimal wavelet parameters and their comparisons will be made. Fourthly, even though wavelet transform is effective in extracting bearing fault features, in-band noises still exist. Morphological analysis with different structure elements and operators will be studied to make bearing fault features sparse and succinct. At last, based on bearing fault features extracted by the suggested steps, intelligent bearing fault diagnosis methods will be developed to automatically identify different bearing faults without requirement of expertise.



---

*Source: 1835127-2015-12-29.xml*