# Fault Feature Extraction and Diagnosis of Gearbox Based on EEMD and Deep Briefs Network

Authors: Kai Chen; Xin-Cong Zhou; Jun-Qiang Fang; Peng-fei Zheng; Jun Wang
Journal: International Journal of Rotating Machinery
 (2017)

---

## Abstract

A gear transmission system is a complex nonstationary and nonlinear time-varying coupling system. When faults occur on gear system, it is difficult to extract the fault feature. In this paper, a novel fault diagnosis method based on ensemble empirical mode decomposition (EEMD) and Deep Briefs Network (DBN) is proposed to treat the vibration signals measured from gearbox. The original data is decomposed into a set of intrinsic mode functions (IMFs) using EEMD, and then main IMFs were chosen for reconstructed signal to suppress abnormal interference from noise. The reconstructed signals were regarded as input of DBN to identify gearbox working states and fault types. To verify the effectiveness of the EEMD-DBN in detecting the faults, a series of gear fault simulate experiments at different states were carried out. Results showed that the proposed method which coupled EEMD and DBN can improve the accuracy of gear fault identification and it is capable of applying to fault diagnosis in practical application.


---

## Body

## 1. Introduction

Gearbox is an indispensable part in the modern industry, especially for lots of large equipment [1]. As a device of transmission power, it is often prone to failure because of its complex structure and the poor working conditions. Failure of the gearbox not only leads to shutting down or threatening personal safety but also causes considerable economic losses. Therefore, it is very important for engineers and researchers to monitor the gear conditions to prevent this kind of malfunction of the plants.The internal structure of the gearbox system is interrelated and coupled inside [2]. During the operation process, a lot of factors need to be considered when evaluating the performance status of the equipment, such as vibration, noise temperature, the debris contaminants in the oil and grease, torque of the power input and output, and stress distribution on the tooth surface. There are many methods for gearbox fault diagnosis, including vibration analysis [3], noise analysis [4], and oil analysis. Being online, real-time, and nondamage, vibration signal analysis shows many advantages which make it widely used for gearbox fault diagnosis.The most common signal processing methods include time-domain analysis and time-frequency analysis [5] such as amplitude spectrum analysis, order analysis [6], cepstrum analysis [7, 8], envelope spectrum analysis [6], Hilbert transform demodulation analysis [9], wavelet analysis, autocorrelation analysis [10], and empirical mode decomposition (EMD) [11]. Among the available vibration analysis methods, EMD is an effective signal analysis method which is suited for dealing with nonlinear and nonstationary signals [12]. It consists in local and fully data-driven separation of a signal in fast and slow oscillations. However, EMD experiences some problems, such as the presence of oscillations of very disparate amplitude in a mode or the presence of very similar oscillations in different modes, named as “mode mixing.” To overcome shortcomings, the ensemble empirical mode decomposition (EEMD) was proposed [13]. EEMD performs the EMD over an ensemble of the signal plus Gaussian white noise. Considering that vibration signals of the gearbox are nonlinear and nonstationary, the EEMD algorithm is suitable for analyzing and judging the gearbox signals [14].There are several failure modes in gear breakdown; some researchers utilized intelligent pattern recognition techniques in fault detection and diagnosis. SVM was used in machine condition monitoring and fault diagnosis [15]. RBF neural network is applied to improve the fault recognition rate in the fault diagnosis of gearbox [16]. Artificial neural network expert system was proposed to detect and localize defects in rolling element bearings [17]. However, the conventional neural network is of slow convergence speed and the possibilities of finding optimal solution are small, which have limited the neural network’s application in the practical projects. But since the Deep Learning was introduced by Hinton et al. in 2006 [18], a new method in machine learning techniques started to be applied in visual data classification [19], decoding analysis [20], language information retrieval [21], and classifying the faults [22]. Especially in 2016, AlphaGo, an outstanding product of the depth learning algorithm, defeated a famous chess player with the score of 4 : 1 [23]. Its excellent performance in the game revealed the power of the learning algorithm. Therefore, it is necessary to introduce such an advanced algorithm into the fault diagnosis of gearbox. In this paper, a new method coupled with EEMD and DBN will be presented to diagnose the gear fault.


## 2. Description of Fault Diagnosis Method

### 2.1. EEMD Algorithm

EMD is a time-frequency signal analysis method of nonlinear signals, which can decompose the data adaptively and obtain a series of IMFs [24]. These IMFs reflect the characteristics of the signal itself and are distributed from high frequency to low frequency. EMD algorithm is very suitable for analyzing the vibration signal of gearbox, because the vibration signal of gearbox is nonstationary.However, in practical application, it is found that there is a problem of mode mixing in the EMD decomposition process [13]; that is, a single IMF contains different frequency components or the same frequency components are decomposed into different IMFs.To alleviate the mode mixing problem occurring in EMD, an ensemble empirical mode decomposition (EEMD) is presented [25]. The essential of EEMD algorithm is to decompose the original signal which adds Gaussian noise by using the EMD method repeatedly, and the original vibration data is decomposed into a series of IMFs with different scales and continuous characteristics because of the characteristic of frequency uniform distribution of Gaussian white noise, which can suppress the appearance of modal mixing.The EEMD algorithm can be summarized as follows:(1) Adding a white noise N(t) to the original signal x(t) can get a new signal X(t):(1)Xt=xt+Nt.(2) The signal X(t) is decomposed by EMD to obtain a set of IMFs:(2)Xt=∑j=1ncjt+rnt.(3) Repeat steps (1) and (2) m times; get m groups of IMF by adding different amplitude of white noise Ni(t), i = 1~m:(3)Xit=xt+NitXit=∑j=1ncijt+rint.(4) Preserve the mean cj(t)(i=1,2,…,m) of each of the m IMFs as the final IMFs:(4)cjt=1m∑i=1mcijt.To verify the validity of EEMD, a simulated signal as shown in formula (5) is decomposed by using EMD and EEMD, respectively; the simulated signal and decomposition results are shown in Figure 1.(5)y1=sin⁡2π30ty2=sin⁡2π10tt>0.2y3=sin⁡2π50tt<0.2y=y1+y2+y3.Figure 1
Comparison results of EMD and EEMD.
(a)
Simulated signal
(b)
EMD results
(c)
EEMD resultsIt can be seen from Figure1 that EEMD can solve the mode mixing phenomenon which is better than EMD. Moreover, IMF2, IMF3, and IMF4 decomposed by EEMD are basically the same as the three components in the original simulated signal.


### 2.2. Deep Belief Networks

#### 2.2.1. Architecture of Restricted Boltzmann Machine

Restricted Boltzmann machine (RBM), a generated random neural network, was proposed by Hinton and Sejnowski in 1986 [26]. There are two components in RBM: the visible layer (v) (input layer) in the bottom and the hidden layer (h) above. All visible units are connected to all hidden units, but there are no visible-visible or hidden-hidden connections. The weight between the visible layer and the hidden layer is denoted by W. The random variables of unit (v,h) take values ∈(0,1), and the joint probability distribution P(v,h) under the model follows Boltzmann distribution [27]. The architecture of RBM is shown in Figure 2.Figure 2
Architecture of restricted Boltzmann machine.The learning process is as follows: when inputting the visible layerv, the hidden layer h is obtained by p(h∣v), and vice versa. In other words, the visible layer is determined by p(h∣v) and the hidden layer. Then the parameter is adjusted until difference between the visible layer and derived visible layer is minimized. It means that the hidden layer obtained in the middle process is an alternative representation of the visible layer.The derivation process of the RBM weight using contrastive divergence (CD) is described as follows [28]:(1)
Randomly choose a set of data in training sample as inputv(1) (superscript denotes train time); determine learning rate L and node number m of the hidden layer.(2)
Randomly initializeW, b, and c.(3)
Update state of hidden variable according to formula (4): (6)Phjk=1∣vk=11+exp⁡-∑iWijvik-bj,wherej=1,2,3,…,m is the node number of the hidden layer.(4)
Reconstitutev(k+1) according to h(k) and formula (5):(7)Pvik+1=1∣hk=11+exp⁡-∑jWijhjk-ci,wherei=1,2,3,…,n is the node number of the input layer.(5)
Then calculateh(k+1) according to v(k+1) and formula (4).(6)
Update weight and bias according to formula (6):(8)W=W+LPhjk=1∣vkvkT-Phjk+1=1∣vk+1vk+1Tc=c+Lvk-vk+1b=b+LPhjk=1∣vk-Phjk+1=1∣vk+1.(7)
Iterate steps (4), (5), and (6) K times; fulfill parameter update of RBM.


#### 2.2.2. Deep Belief Networks

A DBN is formed by stacking a number of the RBMs layer by layer as shown in Figure3. Each layer of the RBM is independent, and the layers are interconnected by weight. Deep belief network structure is shown in Figure 3.Figure 3
Deep belief network structure.A usable DBN network model can be obtained via pretraining and fine-tuning. Note that this learning procedure, so-called pretraining, is unsupervised. The training procedure is described as follows [29]:(1)
Each RBM layer is trained individually, insuring that the output of each layer contains input features as much as possible. The training process is as follows: the input data is mapped to the output data by weight; then the output data try to reconstruct the input data. The weight of the network is updated according to the difference after reconstruction. Repeat this process until the difference between input data and output data is very small; this is the RBM learning process presented in the previous section.(2)
In the previous step, the training only optimizes the mappings between the input layer and the output layer rather than the entire DBN. Therefore, it is necessary to set the last DBN layer as softmax classifier. The output feature vector of RBM is input feature vector of softmax classifier and it is trained under monitoring. All parameters of the network are fine-tuned according to the error of output and input to each RBM layer spreading from top to bottom. The process is illuminated as a flowchart in Figure4.Figure 4
Flowchart of the DBN training process.


## 2.1. EEMD Algorithm

EMD is a time-frequency signal analysis method of nonlinear signals, which can decompose the data adaptively and obtain a series of IMFs [24]. These IMFs reflect the characteristics of the signal itself and are distributed from high frequency to low frequency. EMD algorithm is very suitable for analyzing the vibration signal of gearbox, because the vibration signal of gearbox is nonstationary.However, in practical application, it is found that there is a problem of mode mixing in the EMD decomposition process [13]; that is, a single IMF contains different frequency components or the same frequency components are decomposed into different IMFs.To alleviate the mode mixing problem occurring in EMD, an ensemble empirical mode decomposition (EEMD) is presented [25]. The essential of EEMD algorithm is to decompose the original signal which adds Gaussian noise by using the EMD method repeatedly, and the original vibration data is decomposed into a series of IMFs with different scales and continuous characteristics because of the characteristic of frequency uniform distribution of Gaussian white noise, which can suppress the appearance of modal mixing.The EEMD algorithm can be summarized as follows:(1) Adding a white noise N(t) to the original signal x(t) can get a new signal X(t):(1)Xt=xt+Nt.(2) The signal X(t) is decomposed by EMD to obtain a set of IMFs:(2)Xt=∑j=1ncjt+rnt.(3) Repeat steps (1) and (2) m times; get m groups of IMF by adding different amplitude of white noise Ni(t), i = 1~m:(3)Xit=xt+NitXit=∑j=1ncijt+rint.(4) Preserve the mean cj(t)(i=1,2,…,m) of each of the m IMFs as the final IMFs:(4)cjt=1m∑i=1mcijt.To verify the validity of EEMD, a simulated signal as shown in formula (5) is decomposed by using EMD and EEMD, respectively; the simulated signal and decomposition results are shown in Figure 1.(5)y1=sin⁡2π30ty2=sin⁡2π10tt>0.2y3=sin⁡2π50tt<0.2y=y1+y2+y3.Figure 1
Comparison results of EMD and EEMD.
(a)
Simulated signal
(b)
EMD results
(c)
EEMD resultsIt can be seen from Figure1 that EEMD can solve the mode mixing phenomenon which is better than EMD. Moreover, IMF2, IMF3, and IMF4 decomposed by EEMD are basically the same as the three components in the original simulated signal.


## 2.2. Deep Belief Networks

### 2.2.1. Architecture of Restricted Boltzmann Machine

Restricted Boltzmann machine (RBM), a generated random neural network, was proposed by Hinton and Sejnowski in 1986 [26]. There are two components in RBM: the visible layer (v) (input layer) in the bottom and the hidden layer (h) above. All visible units are connected to all hidden units, but there are no visible-visible or hidden-hidden connections. The weight between the visible layer and the hidden layer is denoted by W. The random variables of unit (v,h) take values ∈(0,1), and the joint probability distribution P(v,h) under the model follows Boltzmann distribution [27]. The architecture of RBM is shown in Figure 2.Figure 2
Architecture of restricted Boltzmann machine.The learning process is as follows: when inputting the visible layerv, the hidden layer h is obtained by p(h∣v), and vice versa. In other words, the visible layer is determined by p(h∣v) and the hidden layer. Then the parameter is adjusted until difference between the visible layer and derived visible layer is minimized. It means that the hidden layer obtained in the middle process is an alternative representation of the visible layer.The derivation process of the RBM weight using contrastive divergence (CD) is described as follows [28]:(1)
Randomly choose a set of data in training sample as inputv(1) (superscript denotes train time); determine learning rate L and node number m of the hidden layer.(2)
Randomly initializeW, b, and c.(3)
Update state of hidden variable according to formula (4): (6)Phjk=1∣vk=11+exp⁡-∑iWijvik-bj,wherej=1,2,3,…,m is the node number of the hidden layer.(4)
Reconstitutev(k+1) according to h(k) and formula (5):(7)Pvik+1=1∣hk=11+exp⁡-∑jWijhjk-ci,wherei=1,2,3,…,n is the node number of the input layer.(5)
Then calculateh(k+1) according to v(k+1) and formula (4).(6)
Update weight and bias according to formula (6):(8)W=W+LPhjk=1∣vkvkT-Phjk+1=1∣vk+1vk+1Tc=c+Lvk-vk+1b=b+LPhjk=1∣vk-Phjk+1=1∣vk+1.(7)
Iterate steps (4), (5), and (6) K times; fulfill parameter update of RBM.


### 2.2.2. Deep Belief Networks

A DBN is formed by stacking a number of the RBMs layer by layer as shown in Figure3. Each layer of the RBM is independent, and the layers are interconnected by weight. Deep belief network structure is shown in Figure 3.Figure 3
Deep belief network structure.A usable DBN network model can be obtained via pretraining and fine-tuning. Note that this learning procedure, so-called pretraining, is unsupervised. The training procedure is described as follows [29]:(1)
Each RBM layer is trained individually, insuring that the output of each layer contains input features as much as possible. The training process is as follows: the input data is mapped to the output data by weight; then the output data try to reconstruct the input data. The weight of the network is updated according to the difference after reconstruction. Repeat this process until the difference between input data and output data is very small; this is the RBM learning process presented in the previous section.(2)
In the previous step, the training only optimizes the mappings between the input layer and the output layer rather than the entire DBN. Therefore, it is necessary to set the last DBN layer as softmax classifier. The output feature vector of RBM is input feature vector of softmax classifier and it is trained under monitoring. All parameters of the network are fine-tuned according to the error of output and input to each RBM layer spreading from top to bottom. The process is illuminated as a flowchart in Figure4.Figure 4
Flowchart of the DBN training process.


## 2.2.1. Architecture of Restricted Boltzmann Machine

Restricted Boltzmann machine (RBM), a generated random neural network, was proposed by Hinton and Sejnowski in 1986 [26]. There are two components in RBM: the visible layer (v) (input layer) in the bottom and the hidden layer (h) above. All visible units are connected to all hidden units, but there are no visible-visible or hidden-hidden connections. The weight between the visible layer and the hidden layer is denoted by W. The random variables of unit (v,h) take values ∈(0,1), and the joint probability distribution P(v,h) under the model follows Boltzmann distribution [27]. The architecture of RBM is shown in Figure 2.Figure 2
Architecture of restricted Boltzmann machine.The learning process is as follows: when inputting the visible layerv, the hidden layer h is obtained by p(h∣v), and vice versa. In other words, the visible layer is determined by p(h∣v) and the hidden layer. Then the parameter is adjusted until difference between the visible layer and derived visible layer is minimized. It means that the hidden layer obtained in the middle process is an alternative representation of the visible layer.The derivation process of the RBM weight using contrastive divergence (CD) is described as follows [28]:(1)
Randomly choose a set of data in training sample as inputv(1) (superscript denotes train time); determine learning rate L and node number m of the hidden layer.(2)
Randomly initializeW, b, and c.(3)
Update state of hidden variable according to formula (4): (6)Phjk=1∣vk=11+exp⁡-∑iWijvik-bj,wherej=1,2,3,…,m is the node number of the hidden layer.(4)
Reconstitutev(k+1) according to h(k) and formula (5):(7)Pvik+1=1∣hk=11+exp⁡-∑jWijhjk-ci,wherei=1,2,3,…,n is the node number of the input layer.(5)
Then calculateh(k+1) according to v(k+1) and formula (4).(6)
Update weight and bias according to formula (6):(8)W=W+LPhjk=1∣vkvkT-Phjk+1=1∣vk+1vk+1Tc=c+Lvk-vk+1b=b+LPhjk=1∣vk-Phjk+1=1∣vk+1.(7)
Iterate steps (4), (5), and (6) K times; fulfill parameter update of RBM.


## 2.2.2. Deep Belief Networks

A DBN is formed by stacking a number of the RBMs layer by layer as shown in Figure3. Each layer of the RBM is independent, and the layers are interconnected by weight. Deep belief network structure is shown in Figure 3.Figure 3
Deep belief network structure.A usable DBN network model can be obtained via pretraining and fine-tuning. Note that this learning procedure, so-called pretraining, is unsupervised. The training procedure is described as follows [29]:(1)
Each RBM layer is trained individually, insuring that the output of each layer contains input features as much as possible. The training process is as follows: the input data is mapped to the output data by weight; then the output data try to reconstruct the input data. The weight of the network is updated according to the difference after reconstruction. Repeat this process until the difference between input data and output data is very small; this is the RBM learning process presented in the previous section.(2)
In the previous step, the training only optimizes the mappings between the input layer and the output layer rather than the entire DBN. Therefore, it is necessary to set the last DBN layer as softmax classifier. The output feature vector of RBM is input feature vector of softmax classifier and it is trained under monitoring. All parameters of the network are fine-tuned according to the error of output and input to each RBM layer spreading from top to bottom. The process is illuminated as a flowchart in Figure4.Figure 4
Flowchart of the DBN training process.


## 3. Experimental Setup and Fault Data

### 3.1. Experimental Setup

To evaluate the performance of the fault diagnosis method, experimental analysis of a gear transmission has been carried out in an experimental setup. As shown in Figure5, the experimental setup is composed of electrical motor, timing belt pulley, coupling, gearbox, and magnetic powder brake. An electrical motor with frequency converter, whose speed is up to 1,500 r/min, is the drive. The number of teeth of the input gear is 55, while that of the output gear is 75; thus the speed ratio is 55/75. A magnetic powder brake is used as the external load. The sensors are the piezoelectric accelerometers (CA-YD-186).Figure 5
Experimental setup and sensor arranged.The amplitude of the vibration signal is attenuated during transmission, and the amplitude attenuation of the high frequency component is much faster than the low frequency component. The location near the bearing, in which attenuation and distortion of the vibration signal were minimum, is the best sensor measuring point. In the experiments, the acceleration sensor is used to measure the vibration signal of the gearbox. And the sensor was arranged in the gearbox near each bearing seat on vertical direction, which is closer to the bearing than the horizontal direction in this experimental setup. The sensor measuring point was shown in Figure5.


### 3.2. Fault Data Collection

In order to collect fault data of gear fault type, a series of gears were processed to simulate common fault gears. As shown in Figure6, three typical faults of gear were simulated, which are tooth breaking, root crack, and pitting corrosion.Figure 6
Simulation of gear fault type.
(a)
Normal
(b)
Tooth breaking
(c)
Pitting corrosion
(d)
Tooth wearExperiment was designed to collect the mass of data under various faults and work conditions, which involve four gear statuses (normal, break, wear, and pit), two types of motor speed (1500 r/min and 1200 r/min), and three kinds of load (0 A, 0.1 A, and 0.2 A). The sampling frequency is 10240 Hz. There are 131072 sampling points in each condition; then the data is divided into 64 parts, each of which is regarded as a sample with 2048 data points. Thus, 384 sets of data were generated in each gear status, and 64 sets of them were selected as test samples by random drawing, while the other 320 sets were chosen as training samples. In total, 1280 training samples and 256 test samples were extracted from the vibration data.


## 3.1. Experimental Setup

To evaluate the performance of the fault diagnosis method, experimental analysis of a gear transmission has been carried out in an experimental setup. As shown in Figure5, the experimental setup is composed of electrical motor, timing belt pulley, coupling, gearbox, and magnetic powder brake. An electrical motor with frequency converter, whose speed is up to 1,500 r/min, is the drive. The number of teeth of the input gear is 55, while that of the output gear is 75; thus the speed ratio is 55/75. A magnetic powder brake is used as the external load. The sensors are the piezoelectric accelerometers (CA-YD-186).Figure 5
Experimental setup and sensor arranged.The amplitude of the vibration signal is attenuated during transmission, and the amplitude attenuation of the high frequency component is much faster than the low frequency component. The location near the bearing, in which attenuation and distortion of the vibration signal were minimum, is the best sensor measuring point. In the experiments, the acceleration sensor is used to measure the vibration signal of the gearbox. And the sensor was arranged in the gearbox near each bearing seat on vertical direction, which is closer to the bearing than the horizontal direction in this experimental setup. The sensor measuring point was shown in Figure5.


## 3.2. Fault Data Collection

In order to collect fault data of gear fault type, a series of gears were processed to simulate common fault gears. As shown in Figure6, three typical faults of gear were simulated, which are tooth breaking, root crack, and pitting corrosion.Figure 6
Simulation of gear fault type.
(a)
Normal
(b)
Tooth breaking
(c)
Pitting corrosion
(d)
Tooth wearExperiment was designed to collect the mass of data under various faults and work conditions, which involve four gear statuses (normal, break, wear, and pit), two types of motor speed (1500 r/min and 1200 r/min), and three kinds of load (0 A, 0.1 A, and 0.2 A). The sampling frequency is 10240 Hz. There are 131072 sampling points in each condition; then the data is divided into 64 parts, each of which is regarded as a sample with 2048 data points. Thus, 384 sets of data were generated in each gear status, and 64 sets of them were selected as test samples by random drawing, while the other 320 sets were chosen as training samples. In total, 1280 training samples and 256 test samples were extracted from the vibration data.


## 4. Application of Proposed Diagnosis Method

### 4.1. EEMD-DBN Diagnosis Method

The original vibration signals were decomposed into several IMFs components with EEMD method. The decomposed IMFs components represent the components of different frequencies in the original vibration signal. EEMD results in one sample of the different fault signals are indicated in Figure7. There are 11 IMFs components in each EEMD result, and each decomposed fault signal is shown from high frequency to low frequency.Figure 7
EEMD results of fault signal.
(a)
Normal
(b)
Break
(c)
Wear
(d)
PitIn real working condition, the collected vibration signal generally contains the noise background, so the pseudocomponent influences the EEMD decomposition. The frequency components contained in these pseudocomponents have the possibility of coinciding with the IMF characteristic frequency band and they should be eliminated. The process of gearbox fault diagnosis based on EEMD noise reduction and DBN is illuminated in Figure8.Figure 8
Procedure of fault diagnosis based on EEMD and DBN.As the natural frequency of the gear is high, the noise in the signal is decomposed into the low frequency band in EEMD. However, the information contains the gearbox working status and fault characteristics which are decomposed into high frequency. Therefore, the high frequency parts of EEMD results are what is mainly studied. Front 4th-order IMFs decomposed by original data were chosen to reconstitute signals. After EEMD noise reduction, DBN is used to classify the four states of the gearbox.


### 4.2. DBN Training Results

Theoretically, the length of the input data should be 2048 points, but the type of input data is the frequency domain of reconstituted signal, which was symmetry, so the first half was taken. This choice will not affect the classification results. Moreover, it can reduce the training time. Therefore, the size of the input data is 1024. The parameters of DBN are shown in Table1.Table 1
Parameters of DBN.
Network layer (including input layer)
4
Type of input data
Frequency domain
Length of the input data
1024
Node number of input layer
1024
Node number of RBM1
500
Node number of RBM2
100
Node number of output layer
4
Classifier of output layer
Softmax
Learning rate
0.1
Maximum iteration
100A comparative analysis was carried out between coupled EEMD with DBN and DBN without EEMD to verify the availability of the coupled method.To increase the separation among the feature clusters in order to gain more diagnostic accuracy and reduce the feature dimensionality for effective computation, the principle component analysis (PCA) is successively applied to the obtained feature sets. The distribution of the features in the new feature space can be observable through the visualization.The first three principal components of first layer and second RBM layer are extracted, and the visual graphic is shown in Figures9 and 10. Figure 9 is the result of DBN training without EEMD method, and Figure 10 is the result of DBN training with EEMD noise reduction.Figure 9
Features distribution of DBN result.
(a)
First RBM layer
(b)
Second RBM layerFigure 10
Features distribution of EEMD-DBN result.
(a)
First RBM layer
(b)
Second RBM layer


### 4.3. DBN Test Results

In order to estimate classification effect of the two training data samples on DBN, the data test was carried out. Test data was processed the same as the training data (if the training data is processed by EEMD, so is the test data; if the training data was unprocessed, so is the test data). The test of the DBN network was carried out 20 times, and the accuracy rate of each time is shown in Figure11.Figure 11
Accuracy of test result.It can be observed from Figure11 that the average recognition rate with EEMD processing is 99.25%, while the average recognition rate without processing is 96.25%. The recognition rate was enhanced almost 3% after EEMD processing. Moreover, recognition rate is more stable when the data is processed. It can be said that DBN methods effectively assist in improving the accuracy of gear fault diagnosis.


## 4.1. EEMD-DBN Diagnosis Method

The original vibration signals were decomposed into several IMFs components with EEMD method. The decomposed IMFs components represent the components of different frequencies in the original vibration signal. EEMD results in one sample of the different fault signals are indicated in Figure7. There are 11 IMFs components in each EEMD result, and each decomposed fault signal is shown from high frequency to low frequency.Figure 7
EEMD results of fault signal.
(a)
Normal
(b)
Break
(c)
Wear
(d)
PitIn real working condition, the collected vibration signal generally contains the noise background, so the pseudocomponent influences the EEMD decomposition. The frequency components contained in these pseudocomponents have the possibility of coinciding with the IMF characteristic frequency band and they should be eliminated. The process of gearbox fault diagnosis based on EEMD noise reduction and DBN is illuminated in Figure8.Figure 8
Procedure of fault diagnosis based on EEMD and DBN.As the natural frequency of the gear is high, the noise in the signal is decomposed into the low frequency band in EEMD. However, the information contains the gearbox working status and fault characteristics which are decomposed into high frequency. Therefore, the high frequency parts of EEMD results are what is mainly studied. Front 4th-order IMFs decomposed by original data were chosen to reconstitute signals. After EEMD noise reduction, DBN is used to classify the four states of the gearbox.


## 4.2. DBN Training Results

Theoretically, the length of the input data should be 2048 points, but the type of input data is the frequency domain of reconstituted signal, which was symmetry, so the first half was taken. This choice will not affect the classification results. Moreover, it can reduce the training time. Therefore, the size of the input data is 1024. The parameters of DBN are shown in Table1.Table 1
Parameters of DBN.
Network layer (including input layer)
4
Type of input data
Frequency domain
Length of the input data
1024
Node number of input layer
1024
Node number of RBM1
500
Node number of RBM2
100
Node number of output layer
4
Classifier of output layer
Softmax
Learning rate
0.1
Maximum iteration
100A comparative analysis was carried out between coupled EEMD with DBN and DBN without EEMD to verify the availability of the coupled method.To increase the separation among the feature clusters in order to gain more diagnostic accuracy and reduce the feature dimensionality for effective computation, the principle component analysis (PCA) is successively applied to the obtained feature sets. The distribution of the features in the new feature space can be observable through the visualization.The first three principal components of first layer and second RBM layer are extracted, and the visual graphic is shown in Figures9 and 10. Figure 9 is the result of DBN training without EEMD method, and Figure 10 is the result of DBN training with EEMD noise reduction.Figure 9
Features distribution of DBN result.
(a)
First RBM layer
(b)
Second RBM layerFigure 10
Features distribution of EEMD-DBN result.
(a)
First RBM layer
(b)
Second RBM layer


## 4.3. DBN Test Results

In order to estimate classification effect of the two training data samples on DBN, the data test was carried out. Test data was processed the same as the training data (if the training data is processed by EEMD, so is the test data; if the training data was unprocessed, so is the test data). The test of the DBN network was carried out 20 times, and the accuracy rate of each time is shown in Figure11.Figure 11
Accuracy of test result.It can be observed from Figure11 that the average recognition rate with EEMD processing is 99.25%, while the average recognition rate without processing is 96.25%. The recognition rate was enhanced almost 3% after EEMD processing. Moreover, recognition rate is more stable when the data is processed. It can be said that DBN methods effectively assist in improving the accuracy of gear fault diagnosis.


## 5. Conclusion

This paper has presented the new approach to diagnose gearbox fault using vibration signals. The original vibration signals were decomposed into several IMFs components with EEMD method. The decomposed IMFs components represent the different frequency components of the original vibration signal. To deal with the instabilities occurring in vibration signal, only main IMFs were chosen to reconstitute signals for eliminating noise effects in signal. The fault features were extracted by DBN training.Finally, the features are inputted to DBN test, so that their performance is appraised. The results show that the proposed approach which coupled EEMD and DBN can improve the accuracy of gear fault identification and it is capable of applying to fault diagnosis in practical application.



---

*Source: 9602650-2017-06-11.xml*