# A Hybrid Domain Degradation Feature Extraction Method for Motor Bearing Based on Distance Evaluation Technique

Authors: Baiyan Chen; Hongru Li; He Yu; Yukui Wang
Journal: International Journal of Rotating Machinery
 (2017)

---

## Abstract

The vibration signal of the motor bearing has strong nonstationary and nonlinear characteristics, and it is arduous to accurately recognize the degradation state of the motor bearing with traditional single time or frequency domain indexes. A hybrid domain feature extraction method based on distance evaluation technique (DET) is proposed to solve this problem. Firstly, the vibration signal of the motor bearing is decomposed by ensemble empirical mode decomposition (EEMD). The proper intrinsic mode function (IMF) component that is the most sensitive to the degradation of the motor bearing is selected according to the sensitive IMF selection algorithm based on the similarity evaluation. Then the distance evaluation factor of each characteristic parameter is calculated by the DET method. The differential method is used to extract sensitive characteristic parameters which compose the characteristic matrix. And then the extracted degradation characteristic matrix is used as the input of support vector machine (SVM) to identify the degradation state. Finally, It is demonstrated that the proposed hybrid domain feature extraction method has higher recognition accuracy and shorter recognition time by comparative analysis. The positive performance of the method is verified.


---

## Body

## 1. Introduction

Rolling bearing is one of the most paramount parts of the motor, which can support rotor and directly affect the operation of the motor [1, 2]. The degraded bearing can give rise to a decline in the performance of the motor and even lead to the improper work of the entire system [3]. Therefore, it is of great significance and practical value to extract degradation feature. As plenty of state information is hidden in the motor bearing vibration signal, the analysis of vibration signal is one of the main methods to extract degradation feature. The degraded bearing itself easily trigger nonconcentric fault of the motor stator and rotor. Consequently, the motor bearing vibration signals show strong nonstationary and nonlinear characteristics [4, 5]. The difficulty and complexity of extracting effective degradation feature is greatly increased. Degradation state recognition of the bearing consists of three steps, including information acquisition, feature extraction, and pattern recognition [6, 7]. Accordingly, it is the key of degradation state recognition to extract degradation feature accurately and effectively.EEMD method is especially applicable to the analysis of nonstationary and nonlinear signals, and the method has a wide range of applications in many fields [8, 9]. Žvokelj et al. [10] utilized EEMD method to process acoustic emission (AE) signal of slewing bearing, and the mixed domain fault feature of each IMF component was extracted through principal component analysis (PCA). Lei et al. [11] combined EEMD and wavelet neural network (WNN) to achieve locomotive roller bearing fault diagnosis. Jiang et al. [8] proposed an improved EEMD method with multiwavelet packet, and it was applied to rotating machinery multifault diagnosis. It is bound to be reflected in a certain IMF component when the bearing fault appears. The degradation characteristics will be more obvious after the amplification of IMF. Therefore, it can achieve a better feature extraction effect based on EEMD. Single domain feature is arduous to fully and accurately describe the different degrees of degradation in complex mechanical systems [12]. A method based on a hybrid domain feature extraction is proposed to make full use of time domain, frequency domain, and time-frequency domain characteristic information.The feature fusion method can make full use of the feature information of multiple domains, which can reflect the degradation state of the equipment more comprehensively. Yu [13] fused multiple domain features with locality preserving projections (LPP) for bearing performance degradation assessment. Chen et al. [14] utilized Local tangent space alignments (LTSA) to fuse multiple domain features for rolling bearing multifault diagnosis study. However, the fused feature does not enjoy a clear physical meaning, and it can not explain the specific performance of the feature as the degenerative state changes. The extracted feature values often appear in the following two conditions [15]. One is the characteristic parameter which has no relation with the classified target, and the other is the redundant characteristic parameter which has high correlation with other characteristic parameters. Feature selection is to select the most effective features from a set of features to reduce the dimension of the feature space [16]. Compared with feature fusion method, the advantage of using feature selection methods for dimension reduction is that the feature is not transformed and still maintains the original physical meaning [17–19]. As a classical algorithm for feature selection, DET can reduce the time of degradation state recognition and ensure the accuracy of recognition at the same time.Based on the above analysis, motor bearing hybrid domain degradation feature extraction based upon EEMD and DET method is proposed in this paper. Firstly, EEMD method is used to decompose the vibration signal of the motor bearing, and the sensitive components are selected according to the sensitive factors. Hybrid domain degradation features are constructed, and sensitive components are selected by DET method. Finally, the degradation state recognition rate obtained by SVM is used to verify the accuracy of the feature extraction.The paper is organized as follows: Section2 presents the EEMD method. And the selection of the sensitive IMF component is also given in this section. Section 3 presents a mathematical analysis in detail to select the most effective features based on DET method. Section 4 presents a favorable discussion and analysis of the experimental results by comparing the proposed method with single domain degradation feature extraction and high dimension degradation feature extraction. Finally, our conclusions are provided in Section 5.


## 2. Algorithm of EEMD

### 2.1. EEMD Method

A new noise-assisted data analysis method was proposed to overcome the problem of mode mixing in EMD [20]. It is EEMD which defines the true IMF components as the mean of an ensemble of trials. Each trial consists of the decomposition results of the signal adding a white noise of finite amplitude.The algorithm is described as follows [21].( 1 ) Add white noise sequence to the signal x(t).(1)x′t=xt+nt,where x′(t) represents the noise-added signal and n(t) indicates the added white noise.( 2 )  x′(t) is decomposed into IMFs according to the EMD method.( 3 ) Repeat steps (1) and (2) by adding a different white noise sequence.( 4 ) The average value of the corresponding IMFs is obtained.


### 2.2. Selection Algorithm Based on the Similarity Evaluation

The signal is decomposed into a group of IMFs by EEMD. A portion of the IMFs are sensitive components closely related to the degradation state, while the others are degenerate independent or noise interfering components. Consequently, it is essential to choose sensitive IMFs which are closely related to the bearing degradation before evaluating the sensitivity of each characteristic parameter. It can improve the accuracy of the degradation feature extraction by ignoring other irrelevant IMFs [22]. A sensitive IMF selection algorithm based on similarity evaluation is adopted in this paper.The specific algorithm of sensitive IMF selection algorithm [23] is as follows.( 1 ) Calculate correlation coefficients un between signal x(t) and its IMFn, where n=1,2,…,N.( 2 ) Calculate correlation coefficients βn between the IMFn of signalx(t) and signalx(t) collected under normal operating conditions.( 3 ) The correlation coefficient of bearing degradation is defined as (2)ηn=βn-un,where n=1,2,…,N.( 4 ) Define and calculate IMFs degradation sensitivity factor λn as(3)λn=ηn-min⁡ηmax⁡η-min⁡η,where n=1,2,…,N.( 5 ) All IMFs are sorted in accordance with the sensitive factor from large to small order, and the difference of the sensitivity factor of two adjacent IMFs is calculated. The index of the corresponding maximum differencenis found. Then the formern IMFs are degradation feature sensitive IMF components.IMF sensitivity factor not only considers the similarity between IMF and the signal itself but also takes into account the similarity between the IMF and the normal signal. The IMFs containing more fault information are more similar to the origin signal. The information that is not related to the fault in IMF can be evaluated through calculating the correlation coefficient between the IMF component and the normal signal. Consequently, the method is more accurate than only calculating the correlation coefficient between IMF component and origin signal in selecting the sensitive IMF components.


## 2.1. EEMD Method

A new noise-assisted data analysis method was proposed to overcome the problem of mode mixing in EMD [20]. It is EEMD which defines the true IMF components as the mean of an ensemble of trials. Each trial consists of the decomposition results of the signal adding a white noise of finite amplitude.The algorithm is described as follows [21].( 1 ) Add white noise sequence to the signal x(t).(1)x′t=xt+nt,where x′(t) represents the noise-added signal and n(t) indicates the added white noise.( 2 )  x′(t) is decomposed into IMFs according to the EMD method.( 3 ) Repeat steps (1) and (2) by adding a different white noise sequence.( 4 ) The average value of the corresponding IMFs is obtained.


## 2.2. Selection Algorithm Based on the Similarity Evaluation

The signal is decomposed into a group of IMFs by EEMD. A portion of the IMFs are sensitive components closely related to the degradation state, while the others are degenerate independent or noise interfering components. Consequently, it is essential to choose sensitive IMFs which are closely related to the bearing degradation before evaluating the sensitivity of each characteristic parameter. It can improve the accuracy of the degradation feature extraction by ignoring other irrelevant IMFs [22]. A sensitive IMF selection algorithm based on similarity evaluation is adopted in this paper.The specific algorithm of sensitive IMF selection algorithm [23] is as follows.( 1 ) Calculate correlation coefficients un between signal x(t) and its IMFn, where n=1,2,…,N.( 2 ) Calculate correlation coefficients βn between the IMFn of signalx(t) and signalx(t) collected under normal operating conditions.( 3 ) The correlation coefficient of bearing degradation is defined as (2)ηn=βn-un,where n=1,2,…,N.( 4 ) Define and calculate IMFs degradation sensitivity factor λn as(3)λn=ηn-min⁡ηmax⁡η-min⁡η,where n=1,2,…,N.( 5 ) All IMFs are sorted in accordance with the sensitive factor from large to small order, and the difference of the sensitivity factor of two adjacent IMFs is calculated. The index of the corresponding maximum differencenis found. Then the formern IMFs are degradation feature sensitive IMF components.IMF sensitivity factor not only considers the similarity between IMF and the signal itself but also takes into account the similarity between the IMF and the normal signal. The IMFs containing more fault information are more similar to the origin signal. The information that is not related to the fault in IMF can be evaluated through calculating the correlation coefficient between the IMF component and the normal signal. Consequently, the method is more accurate than only calculating the correlation coefficient between IMF component and origin signal in selecting the sensitive IMF components.


## 3. DET Method and Degradation Feature Extraction

### 3.1. Characteristic Parameters of Motor Bearing

IMF sensitive component is regarded as the research object after EEMD decomposition. 12 characteristic parameters existing in the time domain, the frequency domain, and the time-frequency domain are extracted from each IMF component, which are conducive to status identification for motor bearings. Time domain parameters include peak value, root mean square amplitude, absolute mean value, root mean square value, peak index, waveform index, and clearance factor. Frequency domain parameters include center of gravity frequency, mean square frequency, root mean square frequency, and frequency variance. Permutation entropy (PE) is determined as a characteristic parameter based on complexity, whose embedding dimensionm=6 and the time delay τ=1 [24]. This paper proposes to use hybrid features displayed in Table 1.Table 1
Hybrid domain features.
Peak value x p = x i m a x  Root mean square amplitude X r = 1 N  ∑ i = 1 N x i 1 / 2 2  Absolute mean value μ x = 1 N  ∑ i = 1 N x i  Root mean square value X r m s = 1 N  ∑ i = 1 N x i 2  Peak index C = X p X r m s  Waveform index I f = X p μ x  Clearance factor C L f = X p X r  Center of gravity frequency F C = ∑ i = 1 N x i ¯ · x i 2 π ∑ i = 1 N x i 2  Mean square frequency M S F = ∑ i = 1 N x i ¯ 2 4 π 2 ∑ i = 1 N x i 2  Root mean square frequency R M S F = M S F  Frequency variance V F = M S F - ( F C ) 2  Permutation entropy H p = - ∑ j = 1 k P j ln ⁡ ( P j ) ln ⁡ ( m ! )x i represents the vibration signal. xi¯=(xi-xi-1)·Fs; Fs represents sampling frequency.N represents sampling points. Pj is used to denote the probability distribution of each symbol sequences.


### 3.2. Calculation of Distance Evaluation Factor

DET is a kind of feature evaluation technology based on inner class distance [25]. The evaluation principle can be described as that the same kind of inner characteristic distance is the smallest, and the distance between two different classes is the biggest. Features complying with this principle are considered to be sensitive features.The specific implementation process of the evaluation method is as follows [25].Assume there are a total ofN species of degradation states of the motor bearing andMsamples are selected for each state andCparameters are extracted from each sample. An obtained characteristic matrix is {fn,m,c,n=1,2,…,N;m=1,2,…,M;c=1,2,…,C}. fn,m,c representscth characteristic ofmth sample of nth degradation state.( 1 ) The inner class average Euclidean distance of the motor bearing in the same state is calculated. (Distance is Euclidean distance without special instructions.)(4)Dn,c=1M×M-1∑l,m=1l≠mMfm,n,c-fl,n,c.The average distance between the classes ofN species is obtained.(5)Dcw=1N∑n=1NDn,c.( 2 ) The difference factor of the inner class distance in the same condition is defined by analyzing the frequency spectrum of the vibration signal.(6)νcw=max⁡Dn,cmin⁡Dn,c.( 3 ) The average value of each feature of all samples in the same class is calculated.(7)un,c=1M∑m=1Mfm,n,c.The average distance between different states of each feature is obtained.(8)Dcb=1N×N-1∑n,e=1n≠eNun,c-ue,c.( 4 ) The status of the motor bearing is examined by the variation of the characteristic parameter. The difference between class distances is defined as follows:(9)νcb=max⁡un,c-ue,cmin⁡un,c-ue,c,where n, e=1,2,…,N; n≠e.( 5 ) Weighting factors are calculated.(10)λc=1νcw/max⁡νcw+νcb/max⁡νcb.( 6 ) The ratio of all parameters average distance in different states and average distance in the same state are calculated considering the influence of weighting factors.(11)αc=λcDcbDcw.( 7 ) The results obtained are normalized.(12)βc=αcmaxαc.β c is the distance evaluation factor of characteristic parameters. It reacts to the sensitivity of each characteristic parameter directly.


### 3.3. Selection of Sensitive Degradation Features Based on Distance Evaluation Factor

It is a knotty problem how to make use of the distance evaluation factor to select the sensitive degradation features to reduce the recognition time. Currently, there is no accepted standard for the selection of thresholds [26]. Quite a few scholars have been inspired by the Wrapper method. They used the classifier to determine the threshold value of the method to screen the sensitive features [25]. However, the method has the shortcomings of high time complexity, large calculation, and poor generalization ability. In addition, a number of tests are needed [19]. The differential method is utilized to screen sensitive degradation characteristics in this paper. Compared with the method of determining the threshold value according to the recognition rate of classifier, the calculation is smaller, and the degradation features can be selected without the intervention of the classifier. Firstly, the characteristic parameters are extracted as degradation state and the operating condition (five kinds of degradation States, four operating conditions) synchronously change and the operating conditions (four conditions, the inner race degradation 0.007 inch) barely change. The distance evaluation factor of each parameter in the two cases is determined by calculating the distance evaluation factor. The difference between the distance evaluation factor of the degradation state and the operating condition at the same time is subtracted from the distance evaluation factor when the condition is changed. The difference between distance evaluation factor for degradation state and operating condition simultaneous changing and distance evaluation factor for only operating conditions changing is obtained. ForD kinds of different conditions, βcd according to the above method can be obtained. The difference between βc and βcd is χc. The characteristic parameters are more sensitive to the change of the state of the motor bearing as χc is greater. The characteristic parameters of χc>0 are selected as the sensitive parameters that exclude the influence of operating condition.


### 3.4. Support Vector Machine (SVM)

Support vector machine (SVM) is an effective computational learning method presented by Vapnik [27], which specializes for a smaller number of samples for training. The basic principle of SVM can be illustrated in a two-dimensional way as represented in Figure 1. It demonstrates the classification of a series of points for two different classes, Class 1 and Class 2. The SVM tries to calculate the distance between the boundary and the nearest data point in each class maximal. The nearest data points are regarded as support vectors.Figure 1
The optimized separating hyper plane in classification.SVM classification function output is a linear combination of intermediate nodes, and each intermediate node corresponds to a support vector. The architecture of SVM is manifested in Figure2.Figure 2
The architecture of SVM.Consider a training sample set,S=(xi,yi)∣xi∈RN,yi∈-1,1,i=1,2,…,l, where xi is an input vector and yi is a label of xi. The general form of linear discriminant function in n-dimensional space is g(x)=(w·x)+b. The data in the set can be correctly classified by the optimized hyperplane. SVM turns to the following dual optimization problem for classification problem:(13)minimize12wTw+C∑i=1nζisubject toyiw·xi+b+ζi≥1,where w andb represent undetermined coefficient.C represents penalty factor. ζi represents relaxation factor.The final result is a discrimination functionf(x) conveniently expressed as a function of the data in the lower dimensional feature space:(14)fx=sgn∑i=1nαiyiKxi,x+b.Kernel function can avoid complex operation in high dimensional feature space. The most widely adopted kernel function is the radial basis function (RBF), and it is defined as follows:(15)Kx,xi=exp⁡-x-xi22σ2.SVM is mainly used as “one to one” binary classifier and “one to others” multiclass classifier. SVM is utilized to identify the different degradation states of the motor bearings. Accordingly, multiclass classifier is determined in this paper. A multiclass classifier consists of an ocean of binary classifier. The process of the multiclass SVM is described in Figure3.Figure 3
“One to others” multiclass SVM classifier.


### 3.5. Hybrid Domain Degradation Feature Extraction Method Based on DET

The process of motor bearing degradation feature extraction based on DET is shown in Figure4.Figure 4
Motor bearing degradation feature extraction process based on DET.( 1 ) The original vibration signal is decomposed by EEMD when the degradation states and the operating conditions synchronously change, and the sensitive IMF components are selected as the analysis object according to the similarity evaluation algorithm.( 2 ) For each sensitive IMF component, 12 characteristic parameters are extracted to form a hybrid domain feature matrix.( 3 ) The distance evaluation factors of characteristic parameters are obtained according to the DET method.( 4 ) Repeat steps (1)–(3) in a special degradation state, only when the operating conditions change.( 5 ) According to the difference method, the sensitive features are selected to form the degradation feature matrix.( 6 ) Finally, the effectiveness of the proposed method is verified by SVM.


## 3.1. Characteristic Parameters of Motor Bearing

IMF sensitive component is regarded as the research object after EEMD decomposition. 12 characteristic parameters existing in the time domain, the frequency domain, and the time-frequency domain are extracted from each IMF component, which are conducive to status identification for motor bearings. Time domain parameters include peak value, root mean square amplitude, absolute mean value, root mean square value, peak index, waveform index, and clearance factor. Frequency domain parameters include center of gravity frequency, mean square frequency, root mean square frequency, and frequency variance. Permutation entropy (PE) is determined as a characteristic parameter based on complexity, whose embedding dimensionm=6 and the time delay τ=1 [24]. This paper proposes to use hybrid features displayed in Table 1.Table 1
Hybrid domain features.
Peak value x p = x i m a x  Root mean square amplitude X r = 1 N  ∑ i = 1 N x i 1 / 2 2  Absolute mean value μ x = 1 N  ∑ i = 1 N x i  Root mean square value X r m s = 1 N  ∑ i = 1 N x i 2  Peak index C = X p X r m s  Waveform index I f = X p μ x  Clearance factor C L f = X p X r  Center of gravity frequency F C = ∑ i = 1 N x i ¯ · x i 2 π ∑ i = 1 N x i 2  Mean square frequency M S F = ∑ i = 1 N x i ¯ 2 4 π 2 ∑ i = 1 N x i 2  Root mean square frequency R M S F = M S F  Frequency variance V F = M S F - ( F C ) 2  Permutation entropy H p = - ∑ j = 1 k P j ln ⁡ ( P j ) ln ⁡ ( m ! )x i represents the vibration signal. xi¯=(xi-xi-1)·Fs; Fs represents sampling frequency.N represents sampling points. Pj is used to denote the probability distribution of each symbol sequences.


## 3.2. Calculation of Distance Evaluation Factor

DET is a kind of feature evaluation technology based on inner class distance [25]. The evaluation principle can be described as that the same kind of inner characteristic distance is the smallest, and the distance between two different classes is the biggest. Features complying with this principle are considered to be sensitive features.The specific implementation process of the evaluation method is as follows [25].Assume there are a total ofN species of degradation states of the motor bearing andMsamples are selected for each state andCparameters are extracted from each sample. An obtained characteristic matrix is {fn,m,c,n=1,2,…,N;m=1,2,…,M;c=1,2,…,C}. fn,m,c representscth characteristic ofmth sample of nth degradation state.( 1 ) The inner class average Euclidean distance of the motor bearing in the same state is calculated. (Distance is Euclidean distance without special instructions.)(4)Dn,c=1M×M-1∑l,m=1l≠mMfm,n,c-fl,n,c.The average distance between the classes ofN species is obtained.(5)Dcw=1N∑n=1NDn,c.( 2 ) The difference factor of the inner class distance in the same condition is defined by analyzing the frequency spectrum of the vibration signal.(6)νcw=max⁡Dn,cmin⁡Dn,c.( 3 ) The average value of each feature of all samples in the same class is calculated.(7)un,c=1M∑m=1Mfm,n,c.The average distance between different states of each feature is obtained.(8)Dcb=1N×N-1∑n,e=1n≠eNun,c-ue,c.( 4 ) The status of the motor bearing is examined by the variation of the characteristic parameter. The difference between class distances is defined as follows:(9)νcb=max⁡un,c-ue,cmin⁡un,c-ue,c,where n, e=1,2,…,N; n≠e.( 5 ) Weighting factors are calculated.(10)λc=1νcw/max⁡νcw+νcb/max⁡νcb.( 6 ) The ratio of all parameters average distance in different states and average distance in the same state are calculated considering the influence of weighting factors.(11)αc=λcDcbDcw.( 7 ) The results obtained are normalized.(12)βc=αcmaxαc.β c is the distance evaluation factor of characteristic parameters. It reacts to the sensitivity of each characteristic parameter directly.


## 3.3. Selection of Sensitive Degradation Features Based on Distance Evaluation Factor

It is a knotty problem how to make use of the distance evaluation factor to select the sensitive degradation features to reduce the recognition time. Currently, there is no accepted standard for the selection of thresholds [26]. Quite a few scholars have been inspired by the Wrapper method. They used the classifier to determine the threshold value of the method to screen the sensitive features [25]. However, the method has the shortcomings of high time complexity, large calculation, and poor generalization ability. In addition, a number of tests are needed [19]. The differential method is utilized to screen sensitive degradation characteristics in this paper. Compared with the method of determining the threshold value according to the recognition rate of classifier, the calculation is smaller, and the degradation features can be selected without the intervention of the classifier. Firstly, the characteristic parameters are extracted as degradation state and the operating condition (five kinds of degradation States, four operating conditions) synchronously change and the operating conditions (four conditions, the inner race degradation 0.007 inch) barely change. The distance evaluation factor of each parameter in the two cases is determined by calculating the distance evaluation factor. The difference between the distance evaluation factor of the degradation state and the operating condition at the same time is subtracted from the distance evaluation factor when the condition is changed. The difference between distance evaluation factor for degradation state and operating condition simultaneous changing and distance evaluation factor for only operating conditions changing is obtained. ForD kinds of different conditions, βcd according to the above method can be obtained. The difference between βc and βcd is χc. The characteristic parameters are more sensitive to the change of the state of the motor bearing as χc is greater. The characteristic parameters of χc>0 are selected as the sensitive parameters that exclude the influence of operating condition.


## 3.4. Support Vector Machine (SVM)

Support vector machine (SVM) is an effective computational learning method presented by Vapnik [27], which specializes for a smaller number of samples for training. The basic principle of SVM can be illustrated in a two-dimensional way as represented in Figure 1. It demonstrates the classification of a series of points for two different classes, Class 1 and Class 2. The SVM tries to calculate the distance between the boundary and the nearest data point in each class maximal. The nearest data points are regarded as support vectors.Figure 1
The optimized separating hyper plane in classification.SVM classification function output is a linear combination of intermediate nodes, and each intermediate node corresponds to a support vector. The architecture of SVM is manifested in Figure2.Figure 2
The architecture of SVM.Consider a training sample set,S=(xi,yi)∣xi∈RN,yi∈-1,1,i=1,2,…,l, where xi is an input vector and yi is a label of xi. The general form of linear discriminant function in n-dimensional space is g(x)=(w·x)+b. The data in the set can be correctly classified by the optimized hyperplane. SVM turns to the following dual optimization problem for classification problem:(13)minimize12wTw+C∑i=1nζisubject toyiw·xi+b+ζi≥1,where w andb represent undetermined coefficient.C represents penalty factor. ζi represents relaxation factor.The final result is a discrimination functionf(x) conveniently expressed as a function of the data in the lower dimensional feature space:(14)fx=sgn∑i=1nαiyiKxi,x+b.Kernel function can avoid complex operation in high dimensional feature space. The most widely adopted kernel function is the radial basis function (RBF), and it is defined as follows:(15)Kx,xi=exp⁡-x-xi22σ2.SVM is mainly used as “one to one” binary classifier and “one to others” multiclass classifier. SVM is utilized to identify the different degradation states of the motor bearings. Accordingly, multiclass classifier is determined in this paper. A multiclass classifier consists of an ocean of binary classifier. The process of the multiclass SVM is described in Figure3.Figure 3
“One to others” multiclass SVM classifier.


## 3.5. Hybrid Domain Degradation Feature Extraction Method Based on DET

The process of motor bearing degradation feature extraction based on DET is shown in Figure4.Figure 4
Motor bearing degradation feature extraction process based on DET.( 1 ) The original vibration signal is decomposed by EEMD when the degradation states and the operating conditions synchronously change, and the sensitive IMF components are selected as the analysis object according to the similarity evaluation algorithm.( 2 ) For each sensitive IMF component, 12 characteristic parameters are extracted to form a hybrid domain feature matrix.( 3 ) The distance evaluation factors of characteristic parameters are obtained according to the DET method.( 4 ) Repeat steps (1)–(3) in a special degradation state, only when the operating conditions change.( 5 ) According to the difference method, the sensitive features are selected to form the degradation feature matrix.( 6 ) Finally, the effectiveness of the proposed method is verified by SVM.


## 4. Experimental Results

### 4.1. Experiment Data Sources

The data of deep groove ball bearings of 6205-2RS JEM SKF from the bearing data center of Case Western Reserve University is used as the experimental raw data in our study [28]. The test bed is manifested in Figure 5. The test rig consists of a 2 Hp (horsepower) motor, a torque converter, a dynamometer, and the control circuit. The vibration signals, including the normal, inner race fault, outer race fault, and rolling element fault signal, are collected by an acceleration sensor installed on the magnetic base shell with a sampling frequency of 12 kHz, and sampling points are 1024. In this paper, the degradation states of inner race serve as an example for analysis, and four operating conditions exist in each state. For each fault mode, data are collected under four operating conditions corresponding to motor speeds of 1730, 1750, 1772, and 1797 r/min. Inner race fault diameter is 0.007′′, 0.014′′, 0.021′′, and 0.028′′. The normal state can be seen as a special degradation state, and there are five kinds of degradation states. The vibration signal of the five degradation states and their frequency spectrum are shown in Figure 6.Figure 5
The bearing experimental system.Figure 6
Time domain and frequency spectrum of five-degradation state vibration signal of the motor bearing.


### 4.2. Degradation Feature Extraction Based on DET

First of all, it is analyzed under the condition of five kinds of degradation states, and four kinds of working conditions exist in each state.The collected signal is decomposed by EEMD method, and the sensitive components are selected based on the similarity assessment of sensitive IMF selection algorithm. The results are shown in Figure7. Each signal can be decomposed to get the first two IMF components for the sensitive IMF components. The sensitive components of the five degradation states are manifested in Figure 8. Feature extraction is performed by selecting the IMF components which contain the main degradation state information. In this paper, 12 mixed domain parameters of the IMF components are extracted as the characteristic value, including peak value, root mean square amplitude, absolute mean value, root mean square value, peak index, waveform index, clearance factor, center of gravity frequency, mean square frequency, root mean square frequency, frequency variance, and permutation entropy. Therefore, 24 characteristic values can be extracted for each signal. The characteristics of the IMF1 component values are t1~t12, and characteristics of IMF2 component values are t13~t24. Then, the degradation features are extracted through the DET method. Before feature selection, the sensitivity to class of each characteristic parameter that is the distance evaluation factor βc has been calculated. In a specific degradation state, the distance evaluation factor βcd is obtained in accordance with the same method only when four kinds of operating conditions change. The differences between βc and βcd named χc are demonstrated in Figure 9. The positive direction of y-axis indicates that it is sensitive to the change of degradation state and is not sensitive to the change of operating condition. It is contrary to the negative direction of y-axis. The sensitive characteristic values are selected by the differential method, and the degradation characteristic matrix is further formed. Distance evaluation factors and their differences are displayed in Table 2.Table 2
Distance evaluation factors and their differences betweenβc and βcd.
Feature parameter β c β c d χ c t 1 0.0315 0.0811 −0.0496 t 2 0.0409 0.0849 −0.0440 t 3 4.6864 e - 04 0.0030 −0.0026 t 4 0.0023 0.0028 - 4.6605 e - 04 t 5 0.0064 0.0321 −0.0258 t 6 0.0174 0.0281 −0.0107 t 7 0.0111 0.0951 −0.0840 t 8 0.0037 0.0927 −0.0890 t 9 0.0114 0.0082 0.0032 t 10 0.0018 0.0082 −0.0064 t 11 0.0452 0.0037 0.0415 t 12 0.0044 0.0079 −0.0035 t 13 0.0118 0.0069 0.0050 t 14 0.0038 0.0062 −0.0024 t 15 0.0765 0.0972 −0.0206 t 16 0.1938 0.1332 0.0606 t 17 0.3834 0.2453 0.1381 t 18 0.4835 0.3354 0.1482 t 19 0.3104 0.2743 0.0361 t 20 0.5367 0.2717 0.2650 t 21 0.2700 0.8912 −0.6212 t 22 1 1 0 t 23 0.0257 0.0101 0.0156 t 24 0.0257 0.0119 0.0138Figure 7
The sensitive IMF components based on similarity evaluation.
(a)
Inner race degradation diameter 0.007′′ (b)
Inner race degradation diameter 0.014′′ (c)
Inner race degradation diameter 0.021′′ (d)
Inner race degradation diameter 0.028′′Figure 8
The sensitive IMF components of the five degradation states.Figure 9
Sensitivity of each characteristic parameter.The selected feature parameters, that is,[t20,t18,t17,t16,t19,t23,t24,t9,t11,t13], can be clearly seen according to Figure 9. Degradation feature is reduced from 24 dimensions to 10 dimensions and the physical meaning of the original feature is maintained.


### 4.3. The Effectiveness of the Degradation Feature Extraction

In order to illustrate the effectiveness of the proposed method, SVM method is used to verify the effectiveness of the degradation feature extraction. The recognition accuracy of SVM output and the recognition time of SVM are used as the basis for judgment [4]. Fifty groups of samples are trained for each degenerate state, and fifty groups of samples are tested. The RBF kernel function is selected as SVM kernel function, and its classification performance is better than other types of nuclear function [29]. The penalty parameterC and kernel function parameter g are selected and the particle swarm optimization (PSO) algorithm is used to optimize the process [30]. The number of particles is set to 20 per degenerate state. Iteration maximal number is set to 200. The range ofC is set from 0.1 to 100, and the range of g is set from 0.01 to 1000. The highest recognition accuracy is taken as the optimization target. The fitness curve of parameter optimization is manifested in Figure 10 (since the iteration of the twelfth generation has been convergent, only the first fifty generations are listed). The time of PSO searching for optimized parameter is 2.852146s (operating platform for R2011b Matlab; the main computer is configured as follows: CPU Intel(R) Core(TM) i5-4590 3.30 GHz, memory 4 G). The result of particle seeking optimal parameters in two-dimensional space is that the optimal penalty parameter C=4 and the kernel function parameter g=0.125.Figure 10
The fitness curve of parameter optimization.The identification result is demonstrated in Figure11. Test sample labels 1, 2, 3, 4, and 5, respectively, represent five degradation states. The final recognition rate reaches 99.6% and the recognition time is 9.902657 s. The results indicate that most of the test data are identified correctly. The effectiveness of the proposed method is verified.Figure 11
Degradation state recognition results.


### 4.4. Comparative Analysis

In order to further verify the effect of the degradation feature extraction based on DET, we use the single domain feature extraction and high dimension feature extraction to compare the proposed feature extraction. The result of degradation state recognition rate and the recognition time is displayed in Table3.Table 3
Recognition rate and recognition time.
Method Recognition rate Recognition time The single domain feature extraction 98.4% 8.854626 s High dimension feature extraction 99.2% 19.293885 s The proposed method 99.6% 9.902657 sFirstly, single domain feature extraction is performed for sensitive IMF components. According to Table1, the root mean square frequency of IMF2 is the most sensitive. This feature is selected as the degradation feature, and the effectiveness of the feature extraction is verified by SVM. The degradation state recognition result is shown in Figure 12.Figure 12
Degradation state recognition results for single domain degradation feature.Secondly, the feature extraction of high dimension degradation is performed for sensitive IMF components, that is, screening sensitive features without the method of DET. 24 characteristics are used to form a degeneration feature matrix, and SVM is used to verify the effectiveness of the feature extraction. The recognition rate is manifested in Figure13.Figure 13
Degradation state recognition results for high dimension degradation feature.Through the comparative analysis, we can draw the following conclusions: the proposed method can extract the degradation state information of motor bearings more effectively. The redundant features of the motor bearing degradation state can be removed, and the degradation feature can be extracted more accurately.


## 4.1. Experiment Data Sources

The data of deep groove ball bearings of 6205-2RS JEM SKF from the bearing data center of Case Western Reserve University is used as the experimental raw data in our study [28]. The test bed is manifested in Figure 5. The test rig consists of a 2 Hp (horsepower) motor, a torque converter, a dynamometer, and the control circuit. The vibration signals, including the normal, inner race fault, outer race fault, and rolling element fault signal, are collected by an acceleration sensor installed on the magnetic base shell with a sampling frequency of 12 kHz, and sampling points are 1024. In this paper, the degradation states of inner race serve as an example for analysis, and four operating conditions exist in each state. For each fault mode, data are collected under four operating conditions corresponding to motor speeds of 1730, 1750, 1772, and 1797 r/min. Inner race fault diameter is 0.007′′, 0.014′′, 0.021′′, and 0.028′′. The normal state can be seen as a special degradation state, and there are five kinds of degradation states. The vibration signal of the five degradation states and their frequency spectrum are shown in Figure 6.Figure 5
The bearing experimental system.Figure 6
Time domain and frequency spectrum of five-degradation state vibration signal of the motor bearing.


## 4.2. Degradation Feature Extraction Based on DET

First of all, it is analyzed under the condition of five kinds of degradation states, and four kinds of working conditions exist in each state.The collected signal is decomposed by EEMD method, and the sensitive components are selected based on the similarity assessment of sensitive IMF selection algorithm. The results are shown in Figure7. Each signal can be decomposed to get the first two IMF components for the sensitive IMF components. The sensitive components of the five degradation states are manifested in Figure 8. Feature extraction is performed by selecting the IMF components which contain the main degradation state information. In this paper, 12 mixed domain parameters of the IMF components are extracted as the characteristic value, including peak value, root mean square amplitude, absolute mean value, root mean square value, peak index, waveform index, clearance factor, center of gravity frequency, mean square frequency, root mean square frequency, frequency variance, and permutation entropy. Therefore, 24 characteristic values can be extracted for each signal. The characteristics of the IMF1 component values are t1~t12, and characteristics of IMF2 component values are t13~t24. Then, the degradation features are extracted through the DET method. Before feature selection, the sensitivity to class of each characteristic parameter that is the distance evaluation factor βc has been calculated. In a specific degradation state, the distance evaluation factor βcd is obtained in accordance with the same method only when four kinds of operating conditions change. The differences between βc and βcd named χc are demonstrated in Figure 9. The positive direction of y-axis indicates that it is sensitive to the change of degradation state and is not sensitive to the change of operating condition. It is contrary to the negative direction of y-axis. The sensitive characteristic values are selected by the differential method, and the degradation characteristic matrix is further formed. Distance evaluation factors and their differences are displayed in Table 2.Table 2
Distance evaluation factors and their differences betweenβc and βcd.
Feature parameter β c β c d χ c t 1 0.0315 0.0811 −0.0496 t 2 0.0409 0.0849 −0.0440 t 3 4.6864 e - 04 0.0030 −0.0026 t 4 0.0023 0.0028 - 4.6605 e - 04 t 5 0.0064 0.0321 −0.0258 t 6 0.0174 0.0281 −0.0107 t 7 0.0111 0.0951 −0.0840 t 8 0.0037 0.0927 −0.0890 t 9 0.0114 0.0082 0.0032 t 10 0.0018 0.0082 −0.0064 t 11 0.0452 0.0037 0.0415 t 12 0.0044 0.0079 −0.0035 t 13 0.0118 0.0069 0.0050 t 14 0.0038 0.0062 −0.0024 t 15 0.0765 0.0972 −0.0206 t 16 0.1938 0.1332 0.0606 t 17 0.3834 0.2453 0.1381 t 18 0.4835 0.3354 0.1482 t 19 0.3104 0.2743 0.0361 t 20 0.5367 0.2717 0.2650 t 21 0.2700 0.8912 −0.6212 t 22 1 1 0 t 23 0.0257 0.0101 0.0156 t 24 0.0257 0.0119 0.0138Figure 7
The sensitive IMF components based on similarity evaluation.
(a)
Inner race degradation diameter 0.007′′ (b)
Inner race degradation diameter 0.014′′ (c)
Inner race degradation diameter 0.021′′ (d)
Inner race degradation diameter 0.028′′Figure 8
The sensitive IMF components of the five degradation states.Figure 9
Sensitivity of each characteristic parameter.The selected feature parameters, that is,[t20,t18,t17,t16,t19,t23,t24,t9,t11,t13], can be clearly seen according to Figure 9. Degradation feature is reduced from 24 dimensions to 10 dimensions and the physical meaning of the original feature is maintained.


## 4.3. The Effectiveness of the Degradation Feature Extraction

In order to illustrate the effectiveness of the proposed method, SVM method is used to verify the effectiveness of the degradation feature extraction. The recognition accuracy of SVM output and the recognition time of SVM are used as the basis for judgment [4]. Fifty groups of samples are trained for each degenerate state, and fifty groups of samples are tested. The RBF kernel function is selected as SVM kernel function, and its classification performance is better than other types of nuclear function [29]. The penalty parameterC and kernel function parameter g are selected and the particle swarm optimization (PSO) algorithm is used to optimize the process [30]. The number of particles is set to 20 per degenerate state. Iteration maximal number is set to 200. The range ofC is set from 0.1 to 100, and the range of g is set from 0.01 to 1000. The highest recognition accuracy is taken as the optimization target. The fitness curve of parameter optimization is manifested in Figure 10 (since the iteration of the twelfth generation has been convergent, only the first fifty generations are listed). The time of PSO searching for optimized parameter is 2.852146s (operating platform for R2011b Matlab; the main computer is configured as follows: CPU Intel(R) Core(TM) i5-4590 3.30 GHz, memory 4 G). The result of particle seeking optimal parameters in two-dimensional space is that the optimal penalty parameter C=4 and the kernel function parameter g=0.125.Figure 10
The fitness curve of parameter optimization.The identification result is demonstrated in Figure11. Test sample labels 1, 2, 3, 4, and 5, respectively, represent five degradation states. The final recognition rate reaches 99.6% and the recognition time is 9.902657 s. The results indicate that most of the test data are identified correctly. The effectiveness of the proposed method is verified.Figure 11
Degradation state recognition results.


## 4.4. Comparative Analysis

In order to further verify the effect of the degradation feature extraction based on DET, we use the single domain feature extraction and high dimension feature extraction to compare the proposed feature extraction. The result of degradation state recognition rate and the recognition time is displayed in Table3.Table 3
Recognition rate and recognition time.
Method Recognition rate Recognition time The single domain feature extraction 98.4% 8.854626 s High dimension feature extraction 99.2% 19.293885 s The proposed method 99.6% 9.902657 sFirstly, single domain feature extraction is performed for sensitive IMF components. According to Table1, the root mean square frequency of IMF2 is the most sensitive. This feature is selected as the degradation feature, and the effectiveness of the feature extraction is verified by SVM. The degradation state recognition result is shown in Figure 12.Figure 12
Degradation state recognition results for single domain degradation feature.Secondly, the feature extraction of high dimension degradation is performed for sensitive IMF components, that is, screening sensitive features without the method of DET. 24 characteristics are used to form a degeneration feature matrix, and SVM is used to verify the effectiveness of the feature extraction. The recognition rate is manifested in Figure13.Figure 13
Degradation state recognition results for high dimension degradation feature.Through the comparative analysis, we can draw the following conclusions: the proposed method can extract the degradation state information of motor bearings more effectively. The redundant features of the motor bearing degradation state can be removed, and the degradation feature can be extracted more accurately.


## 5. Conclusions

In this paper, a hybrid domain degradation feature extraction method based on DET is proposed to extract the degradation characteristics. Finally, the recognition rate of SVM output demonstrates the validity and accuracy of the feature extraction method. Through the analysis, we can get the following conclusions.( 1 ) The hybrid domain feature is more comprehensive than the single domain feature, which can reflect the degradation state of the bearing accurately. The proposed method can be overwhelmingly accurate to select the appropriate characteristic parameters of the motor bearing degradation state. It can achieve the degradation feature extraction more accurately.( 2 ) The DET method can remove redundant features and unrelated features effectively, which can maintain the physical meaning of the original feature parameters and shorten the recognition time.( 3 ) The sensitive feature selection is related to the quality of the feature extraction. The degradation feature is extracted through the proposed method, which is applied to the actual vibration signal of the motor bearing. The recognition rate can reach 99.6%, and the result indicates the validity and feasibility of the proposed method.



---

*Source: 2607254-2017-01-24.xml*