# Combined Failure Diagnosis of Slewing Bearings Based on MCKD-CEEMD-ApEn

Authors: Fengtao Wang; Chenxi Liu; Wensheng Su; Zhigang Xue; Qingkai Han; Hongkun Li
Journal: Shock and Vibration
 (2018)

---

## Abstract

Large-size and heavy-load slewing bearings, which are mainly used in heavy equipment, comprise a subgroup of rolling bearings. Owing to the complexity of the structures and working conditions, it is quite challenging to effectively diagnose the combined failure and extract fault features of slewing bearings. In this study, a method was proposed to denoise and classify the combined failure of slewing bearings. First, after removing the mean, the vibration signals were denoised by maximum correlated kurtosis deconvolution. The signals were then decomposed into several intrinsic mode functions (IMFs) by complementary ensemble empirical mode decomposition (CEEMD). Appropriate IMFs were selected based on the correlation coefficient and kurtosis. The approximate entropy values of the selected IMFs were regarded as the characteristic vectors and then inputted into the support vector machine (SVM) based on multiclass classification for training. The practical combined failure signals of the 3 conditions were finally recognized and classified using SVMs. The study also compared the proposed method with 5 other methods to demonstrate the superiority and effectiveness of the proposed method.


---

## Body

## 1. Introduction

Large-size and heavy-load slewing bearings are widely used in major and heavy equipment. Faults more easily occur in slewing bearings, compared with other machine parts under practical working conditions [1]. Owing to the strong background noise interference, extracting combined fault features from vibration signals has long been the focus of current studies. A common approach is to diagnose the faults of time-domain vibration signals after the denoising process [2–5], aiming to increase the signal-to-noise ratio (SNR) and reduce the noise interference.The maximum correlated kurtosis deconvolution (MCKD) algorithm was first proposed by McDonald et al. [6]. This method can increase the correlated kurtosis values of the signals as well as the SNR of the weak fault signals by deconvolution. Meanwhile, the submerged successive pulses can be given emphasis. Compared with the minimum entropy deconvolution [7] algorithm, the MCKD algorithm can also eliminate strong background noise interference to achieve predenoising. The MCKD algorithm is currently applied in many fault diagnosis studies on large-size rolling bearings. Jia et al. [8] combined MCKD with improved spectrum kurtosis to diagnose vibration signals and extract the fault features of rolling bearings in wind turbines and hot-strip rolling mills. Zhao and Li [9] used MCKD and empirical mode decomposition (EMD) method to extract the early-stage fault features of rolling bearings in wind turbines from strong background noise. The present study focused on the denoising effect of the MCKD algorithm and the time-domain analysis of the vibration signals. Thus, the study of frequency spectrum was not demonstrated.Complementary ensemble empirical mode decomposition (CEEMD) [10] was proposed based on ensemble empirical mode decomposition (EEMD) [11]. In CEEMD, a pair of opposite white noises is added to the original signals first. The recombined signals are then decomposed by EMD to obtain a series of intrinsic mode functions (IMFs). Compared with EEMD, CEEMD can reduce the reconstruction error of the added white noise as well as mode mixing under the same decomposition effects. As a signal decomposition method, CEEMD is typically combined with other techniques to diagnose faults in machines [12–14]. In our study, CEEMD was used to denoise signals and extract characteristic parameters for further fault classification.Currently, some research achievements on the combination of signal decomposition methods and many kinds of entropy have been published to apply on the fault diagnosis and classification of rolling bearings. Entropy is a nonlinearity index for evaluating the irregularity of signals. As a category of significant phase-space indexes, many kinds of entropy were used as characteristic parameters in fault diagnosis and classification of rolling bearings [15–18]. Approximate entropy (ApEn), proposed by Pincus [19], can quantize vibration signals and achieve feature extraction, which exhibits strong anti-interference. The more complex the time series, the bigger the ApEn value. ApEn has been considered in recent years as an effective characteristic parameter in condition monitoring and fault diagnosis [20–22]. Caesarendra et al. [23] analyzed four nonlinear features, including ApEn, largest Lyapunov exponent (LLE), and correlation dimension (CD), to provide more superior descriptive information about slewing bearings than time-domain features. Although most studies focused on a single artificial fault, this study focused on practical combined faults under different conditions. This study used ApEn as the characteristic parameter to diagnose the faults and extract the features of large-size and heavy-load slewing bearings under various working conditions. For comparison, sample entropy (SE), LLE, and CD were also taken as the characteristic parameters to justify the benefits of ApEn in this study.Support vector machine (SVM) is a kind of learning algorithm proposed by Vapnik [24] in 1998. In a classification process, SVM separates a set of labeled training data with a hyperplane that maximizes its distance to the data. So, the SVM can work in combination with kernel functions to compute a nonlinear mapping to the features space. The hyperplane corresponds to a nonlinear decision boundary in the input space [25]. SVM has strong function fitting advantages and it is suited for resolving small-sample and nonlinear model identification. Compared with other intelligent classifiers, this one requires a smaller sample size to seek for the optimal solution based on finite samples, which can be taken as a convex optimization problem. The existence and acquisition of the globally optimal solution can be guaranteed. With sound robustness, the empirical risk and confidence risk are taken into consideration in SVM for resolving small-sample problems. The optimal solution based on finite information rather than infinite samples is acquired. In terms of popularization, SVM is more appropriate to be used in fault diagnosis of small-sample conditions. Xiong et al. [26] developed a scheme using SVMs for diagnosing bearing conditions, where a 97.42% accuracy is achieved. Zhang et al. [27] proposed a fault diagnosis scheme for locomotive roller bearings using an SVM classifier. Compared with some methods based on other neural networks, SVM has both a simple structure and an improved generalization capability.The original vibration signals were denoised by removing the mean and the MCKD algorithm as well as conducting decomposition by CEEMD to obtain IMFs. The selected IMFs, which could represent the features of the original signals, were extracted ApEn values. The values consisted of characteristic vectors and were then inputted into the SVMs to train in order to identify the faults of the slewing bearings. All the vibration data of the slewing bearings applied in this study resulted from the practical working conditions. The remaining sections of this paper are arranged as follows. The theories supporting the proposed method are presented in Section2. The practical signal verification conducted on the slewing bearing of a practical port crane is demonstrated in Section 3. The last section provides the summary and conclusions of the study.


## 2. Theoretical Framework

### 2.1. Maximum Correlated Kurtosis Deconvolution

MCKD uses correlated kurtosis as its index and emphasizes periodic pulse components submerged in the background noise by deconvolution computation. Iteration ends when the correlated kurtosis values of the original signals reach the peak. The algorithm achieves denoising for the impact failure of rolling bearings and highlights shock pulses. The vibration signal when faults occur in a rolling bearing can be described as follows:(1)yn=hn∗xn+en.The MCKD algorithm aims to determine the FIR filter w(l), which can achieve the maximal correlated kurtosis value of the original signal. The input signal x(n) can be obtained using the output signal y(n) using the following:(2)xn=wn∗yn.The correlated kurtosis of the MCKD algorithm is defined as(3)CKMT=∑n=1N∏m=0Myn-mT2∑n=1Ny2nM+1,where T represents the period of the impulse signal and M denotes the shift order. The higher the value of M, the greater the impulse sequences of deconvolution, which indicates that the fault-detecting ability is stronger. The optimized objective function of MCKD can be described as follows:(4)MCKDMT=maxwl⁡∑n=1N∏m=0Myn-mT2∑n=1Nyn2M+1,where l=1,2,…,L. The optimized filter w(l) should satisfy the equation(5)ddwlCKmT=0.The results of w(l) can be described in the form of a matrix as follows:(6)w=y22B2X0X0T-1∑m=0MXmTAm,where(7)Xr=x1-rx2-r⋯xN-r0x1-r⋯xN-1-R⋮⋮⋮⋮00⋯xN-L-r+1L×N,r=0T2T⋯mT,Am=y1-mT-1y12y1-T2⋯y1-mT2⋮yN-mT-1yN2yN-T2⋯yN-mT2N×1,B=y1y1-T⋯y1-MT⋮yNyN-T⋯yN-MTN×1.Finally, the coefficient inputs are incorporated into (2) to obtain x(n). This study used MCKD for signal denoising. Different parameters affected the result of denoising. In general, the shift order is selected from 1 to 7 [28]. Although a higher value can obtain more successive pulses, the period of the impulse signal should be estimated more accurately; in addition, the computational cost can increase. The practical fault frequency also varies from the theoretical result. With all factors considered, the shift order selected is 4. The filter size is another important parameter. This study determined the sizes as the kurtosis values of the vibration signals from different working conditions reaching their maximum.


### 2.2. Complementary Ensemble Empirical Mode Decomposition

As a signal decomposition method, CEEMD adds a pair of Gaussian white noises with opposite signs, which can not only ensure the continuity of signals but also completely neutralize the added noise. First,I power white noise is added in the original signal x(n), i=1,2,…,I. The numbers of white noises are added in pairs to obtain 2I signals:(8)x+x-=111-1xωi.Each signal xi(n) in the set is decomposed by EMD; a group of IMFki(n) can then be obtained. The first one can be described as IMF1n=1/2I∑i=12IIMFkin and the first residual error can be described as r1(n)=x(n)-IMF1(n). As k=2,3,…,K, the kth residual error can be calculated as follows:(9)rkn=rk-1n-IMFkn.The decomposition function becomes rkn+ξkEk(ωi(n)):(10)IMFk+1n=12I∑i=12IE1rk+ξkEkωin,where ξk represents SNR and Ej(·) denotes the jth decomposition function after EMD. The aforementioned steps are repeated until the selection is ended. The signal can ultimately be described as follows:(11)xn=∑k=1kIMFnk+Rn.We referred to the study by Huang and Wu [29] to decrease the error caused by the residual noise below 1%, and the amplitude of white noise was set to 1/10th of the original signal.


### 2.3. Approximate Entropy

ApEn can be used to quantify the vibration signals on the basis of the complexity of its information, which can be applied in feature extraction. As for the time seriesx1,xx,…,x(N), which contains N data, the ApEn algorithm predetermines the mode dimension m to reconstruct the phase space. The elements in the time series are extracted to form an m-D vector:(12)Xi=xi,xi+1,…,xi+m-1.The maximum difference between X(i) and X(j) is defined as the distance:(13)dXi,Xj=maxk=0,1,…,m-1⁡xi+k-xj+k.On the basis of the predetermined threshold value r, the number of distances lower than r accumulates:(14)Cimr=1N-m+1NumdXi,Xj<r,where i,j=1,2,…,N-m+1, i≠j. Cim(r) represents the degree of correlation between X(i) and X(j). Φm(r) denotes the degree of autocorrelation of Xi:(15)Φmr=1N-m+1∑i=1N-m+1ln⁡Cimr.The mode dimension is expanded from m to m+1. Φm+1(r) can be obtained after the repetition of the aforementioned steps. The ApEn value of the time series can be obtained based on Φm(r) and Φm+1(r):(16)ApEnm,r=limN→∞⁡Φmr-Φm+1r.This study sets ApEn as the feature parameter. To obtain improved statistical property and reduce error, the threshold value was set to 1/10th of the standard deviation [19].


### 2.4. Support Vector Machines

In order to resolve the 3-classification problem in this study, one-against-one method is adopted, in which SVMs are constructed for all possible pairs of classes including (I) normal condition–early fault condition (SVM1), (II) normal condition–severe fault condition (SVM2), and (III) early fault condition–severe fault condition (SVM3). The decision function that separates theith from the jth class for the features vector is defined as follows:(17)dijx≜wijTx+bij,where i,j=1,2,…,N, wij∈Rn is a weight vector, xi∈Rn is the ith features sample, and bij is a bias term. The number of times that the features vector x is assigned to the ith class is computed by(18)dix=∑j=1j≠iNsigndijx.The following decision function is considered:(19)dx=argmaxi=1,2,…,N⁡dix.The highest number of votes that a feature x vector may have is N-1. For the design of the SVMs, the Gaussian RBF kernel is considered:(20)Kxi,xj=exp⁡-xi-xj22σ2,where xi-xj denotes the Euclidean norm and σ is a free parameter related to the dispersion of the support vectors. RBF functions have been documented extensively in the literature, making them a benchmark for classification applications based on SVMs [25].


### 2.5. Fault Diagnosis Based on MCKD-CEEMD-ApEn

Figure1 demonstrates the specific implementation steps of the proposed method.Figure 1
Flow chart of the proposed method.Step 1.
Remove the mean of the original signals and predenoise the signals by MCKD.Step 2.
Decompose the denoised signals by CEEMD and select appropriate IMFs based on kurtosis and the correlation coefficient.Step 3.
Calculate the ApEn values of each selected IMF and compose the characteristic vectors.Step 4.
Input the characteristic vectors into the SVM for training and verify the fault classification rate and validity of the proposed method.


## 2.1. Maximum Correlated Kurtosis Deconvolution

MCKD uses correlated kurtosis as its index and emphasizes periodic pulse components submerged in the background noise by deconvolution computation. Iteration ends when the correlated kurtosis values of the original signals reach the peak. The algorithm achieves denoising for the impact failure of rolling bearings and highlights shock pulses. The vibration signal when faults occur in a rolling bearing can be described as follows:(1)yn=hn∗xn+en.The MCKD algorithm aims to determine the FIR filter w(l), which can achieve the maximal correlated kurtosis value of the original signal. The input signal x(n) can be obtained using the output signal y(n) using the following:(2)xn=wn∗yn.The correlated kurtosis of the MCKD algorithm is defined as(3)CKMT=∑n=1N∏m=0Myn-mT2∑n=1Ny2nM+1,where T represents the period of the impulse signal and M denotes the shift order. The higher the value of M, the greater the impulse sequences of deconvolution, which indicates that the fault-detecting ability is stronger. The optimized objective function of MCKD can be described as follows:(4)MCKDMT=maxwl⁡∑n=1N∏m=0Myn-mT2∑n=1Nyn2M+1,where l=1,2,…,L. The optimized filter w(l) should satisfy the equation(5)ddwlCKmT=0.The results of w(l) can be described in the form of a matrix as follows:(6)w=y22B2X0X0T-1∑m=0MXmTAm,where(7)Xr=x1-rx2-r⋯xN-r0x1-r⋯xN-1-R⋮⋮⋮⋮00⋯xN-L-r+1L×N,r=0T2T⋯mT,Am=y1-mT-1y12y1-T2⋯y1-mT2⋮yN-mT-1yN2yN-T2⋯yN-mT2N×1,B=y1y1-T⋯y1-MT⋮yNyN-T⋯yN-MTN×1.Finally, the coefficient inputs are incorporated into (2) to obtain x(n). This study used MCKD for signal denoising. Different parameters affected the result of denoising. In general, the shift order is selected from 1 to 7 [28]. Although a higher value can obtain more successive pulses, the period of the impulse signal should be estimated more accurately; in addition, the computational cost can increase. The practical fault frequency also varies from the theoretical result. With all factors considered, the shift order selected is 4. The filter size is another important parameter. This study determined the sizes as the kurtosis values of the vibration signals from different working conditions reaching their maximum.


## 2.2. Complementary Ensemble Empirical Mode Decomposition

As a signal decomposition method, CEEMD adds a pair of Gaussian white noises with opposite signs, which can not only ensure the continuity of signals but also completely neutralize the added noise. First,I power white noise is added in the original signal x(n), i=1,2,…,I. The numbers of white noises are added in pairs to obtain 2I signals:(8)x+x-=111-1xωi.Each signal xi(n) in the set is decomposed by EMD; a group of IMFki(n) can then be obtained. The first one can be described as IMF1n=1/2I∑i=12IIMFkin and the first residual error can be described as r1(n)=x(n)-IMF1(n). As k=2,3,…,K, the kth residual error can be calculated as follows:(9)rkn=rk-1n-IMFkn.The decomposition function becomes rkn+ξkEk(ωi(n)):(10)IMFk+1n=12I∑i=12IE1rk+ξkEkωin,where ξk represents SNR and Ej(·) denotes the jth decomposition function after EMD. The aforementioned steps are repeated until the selection is ended. The signal can ultimately be described as follows:(11)xn=∑k=1kIMFnk+Rn.We referred to the study by Huang and Wu [29] to decrease the error caused by the residual noise below 1%, and the amplitude of white noise was set to 1/10th of the original signal.


## 2.3. Approximate Entropy

ApEn can be used to quantify the vibration signals on the basis of the complexity of its information, which can be applied in feature extraction. As for the time seriesx1,xx,…,x(N), which contains N data, the ApEn algorithm predetermines the mode dimension m to reconstruct the phase space. The elements in the time series are extracted to form an m-D vector:(12)Xi=xi,xi+1,…,xi+m-1.The maximum difference between X(i) and X(j) is defined as the distance:(13)dXi,Xj=maxk=0,1,…,m-1⁡xi+k-xj+k.On the basis of the predetermined threshold value r, the number of distances lower than r accumulates:(14)Cimr=1N-m+1NumdXi,Xj<r,where i,j=1,2,…,N-m+1, i≠j. Cim(r) represents the degree of correlation between X(i) and X(j). Φm(r) denotes the degree of autocorrelation of Xi:(15)Φmr=1N-m+1∑i=1N-m+1ln⁡Cimr.The mode dimension is expanded from m to m+1. Φm+1(r) can be obtained after the repetition of the aforementioned steps. The ApEn value of the time series can be obtained based on Φm(r) and Φm+1(r):(16)ApEnm,r=limN→∞⁡Φmr-Φm+1r.This study sets ApEn as the feature parameter. To obtain improved statistical property and reduce error, the threshold value was set to 1/10th of the standard deviation [19].


## 2.4. Support Vector Machines

In order to resolve the 3-classification problem in this study, one-against-one method is adopted, in which SVMs are constructed for all possible pairs of classes including (I) normal condition–early fault condition (SVM1), (II) normal condition–severe fault condition (SVM2), and (III) early fault condition–severe fault condition (SVM3). The decision function that separates theith from the jth class for the features vector is defined as follows:(17)dijx≜wijTx+bij,where i,j=1,2,…,N, wij∈Rn is a weight vector, xi∈Rn is the ith features sample, and bij is a bias term. The number of times that the features vector x is assigned to the ith class is computed by(18)dix=∑j=1j≠iNsigndijx.The following decision function is considered:(19)dx=argmaxi=1,2,…,N⁡dix.The highest number of votes that a feature x vector may have is N-1. For the design of the SVMs, the Gaussian RBF kernel is considered:(20)Kxi,xj=exp⁡-xi-xj22σ2,where xi-xj denotes the Euclidean norm and σ is a free parameter related to the dispersion of the support vectors. RBF functions have been documented extensively in the literature, making them a benchmark for classification applications based on SVMs [25].


## 2.5. Fault Diagnosis Based on MCKD-CEEMD-ApEn

Figure1 demonstrates the specific implementation steps of the proposed method.Figure 1
Flow chart of the proposed method.Step 1.
Remove the mean of the original signals and predenoise the signals by MCKD.Step 2.
Decompose the denoised signals by CEEMD and select appropriate IMFs based on kurtosis and the correlation coefficient.Step 3.
Calculate the ApEn values of each selected IMF and compose the characteristic vectors.Step 4.
Input the characteristic vectors into the SVM for training and verify the fault classification rate and validity of the proposed method.


## 3. Experimental Evaluation

All the time-domain signals used in this study are selected from the large-size, heavy-load slewing bearing of a practical port crane. The type of the slewing bearing type was 132.45.2800.03. The slewing bearing began to serve from Jul. 2013 and its vibration data were collected from Feb. 2015 to Dec. 2017 about every 2 months. No artificial defect was manufactured on the slewing bearing for fault simulation. The complicated practical working condition and equipment led to combined faults occurring on the slewing bearing. In order to reduce human interference as much as possible, the outer ring, inner ring, rolling elements, and gear teeth of the slewing bearing did not, on the condition of operational safety, undergo any maintenance except for a disassembly check every 6 to 7 months and regular lubrication. Except for the periods of data collection, the slewing bearing maintained normal operation process. For analysis convenience, during the periods of data collection, the rotating speed of the inner race was 1.5 rpm and the rotation direction was counterclockwise. The container carried by the crane remained the same each time. Each acquisition duration of the vibration signals collected using PCB acceleration sensors was longer than a rotation cycle. The sensors were placed on the radial surface at a 90° angle relative to one another around the ring (Figure2). As shown in Figure 3, the vibration signals under 3 different conditions were collected: those under the normal condition, early fault condition, and severe fault condition. The sampling frequency was 512 Hz. To decrease the calculated amount and for ease of comparison, a one-rotation-cycle (40 s) vibration signal was acquired each time. A single fault rarely occurs in practical working slewing bearings, unlike the theoretical analysis as well as the experimental simulation. In addition, according to the disassembly check records (Figure 4), combined faults at varying degrees were observed in the slewing bearing in the early fault and severe fault stages. The working conditions of the slewing bearing were divided into 3 stages based on the deterioration degree of the combined failure. In the normal working stage (from Feb. 2015 to Jun. 2016), the whole bearing was complete, and all components were of normal size. No obvious pitting or abrasion occurred; the lubrication was sound, and the noise was weak. In the early fault stage (from Jul. 2016 to Aug. 2017), large-scale pitting or abrasion occurred in the outer race and the inner race of the slewing bearing. In addition, some rolling elements, retainers, and gear teeth were damaged. Several metal particles could be found in the lubricating grease. In the severe fault stage (from Sep. 2017 to Dec. 2017), abrasion occurred in most rolling elements of the three rows. More serious damage, abrasion, and cracks occurred in the outer gear teeth, some of which were broken. The effects were apparently accompanied by strong noise. Many metal impurities could be found in the lubricating grease. Thus, the combined failure could be determined in the early and severe fault stages.Figure 2
Arrangement of sensors.Figure 3
Time-domain signals of various conditions (the data were acquired on Oct. 15, 2015; Oct. 25, 2016; and Oct. 26, 2017, resp.).
(a)
Normal condition (b)
Early fault condition (c)
Severe fault conditionFigure 4
Faults of various conditions.To verify the superiority and validity of the proposed method, 3 control methods were conducted in this study as follows:(1)
Decompose the signals by CEEMD directly after removing the mean without predenoising by MCKD.(2)
Decompose the signals by EEMD rather than CEEMD after MCKD.(3)
Use SE/LLE/CD, rather than ApEn, as the characteristic indexes.The flow in Figure5 shows the comparison among the 4 methods.Figure 5
Methods comparison.First, the direct-current components of the 3 original signals in Figure2 were eliminated by removing the mean. The signals were predenoised using the MCKD algorithm. Figure 6 shows the time-domain signals after MCKD. Part of the background noise was eliminated, and the effects became more apparent.Figure 6
Time-domain signals of various conditions after MCKD (the data were acquired on Oct. 15, 2015; Oct. 25, 2016; and Oct. 26, 2017, resp.).
(a)
Normal condition (b)
Early fault condition (c)
Severe fault conditionThe signals were then decomposed by CEEMD, and the results are shown in Figures7(a), 8(a), and 9(a). For comparison, Figures 7(b), 8(b), and 9(b) show the EEMD results of the corresponding signals.Figure 7
Signal decomposition of Stage (A).
(a)
CEEMD (b)
EEMDFigure 8
Signal decomposition of Stage (B).
(a)
CEEMD (b)
EEMDFigure 9
Signal decomposition of Stage (C).
(a)
CEEMD (b)
EEMDThe time-domain signals are decomposed into a series of IMFs by CEEMD and EEMD. Although some of them contain plenty of failure feature components, others contain stronger background noise. Therefore, the IMFs should be selected before signal reconstruction, aiming to increase the SNR. In addition, the noise components are not to be considered. Kurtosis is a characteristic parameter that is sensitive to impact components. The value of kurtosis increases with the number of impact components in the signal. Correlation coefficient reflects the degree of correlation between the decomposed signals and the corresponding original signals. Therefore, on the condition of maintaining the basic features of the original signals, in order to extract fault features of slewing bearing from strong background noise, kurtosis and correlation coefficient are used to select IMFs to extract the fault features of the slewing bearing from strong background noise. The decomposition results of the vibration signals from the 3 stages are presented in Figures10, 11, and 12 with different colors, representing MCKD-CEEMD, MCKD-EEMD, and CEEMD, respectively.Figure 10
Characteristic parameters of Stage (A).
(a)
Kurtosis (b)
Correlation coefficientFigure 11
Characteristic parameters of Stage (B).
(a)
Kurtosis (b)
Correlation coefficientFigure 12
Characteristic parameters of Stage (C).
(a)
Kurtosis (b)
Correlation coefficientThe cumulative percentages of kurtosis and the correlation coefficient about the firstn IMFs of different signals comprising all of the corresponding features are shown in Tables 1 and 2. After the calculation of the first 7 IMFs, all cumulative percentages of the correlation coefficient values are not lower than 85%, and the cumulative percentages are all not lower than 50% (except for the result of MCKD-EEMD in the normal stage). Therefore, the first 7 IMFs, which contain the majority of the principal characteristics and achieve denoising to a certain extent, are selected to represent the original signals of the corresponding stages.Table 1
Ratio of parameter: kurtosis.
Methods IMFn Stages 1 2 3 4 5 6 7 Ratio (%)  MCKD-CEEMD 5.16 12.1 19.5 27.5 35.0 42.5 50.8 (A) 5.28 13.2 21.2 29.9 38.8 48.7 57.7 (B) 11.9 24.1 40.2 53.1 58.8 64.8 69.5 (C)  MCKD-EEMD 4.36 10.3 16.6 23.5 29.5 36.3 45.8 (A) 4.74 11.7 19.0 27.5 36.1 43.2 51.9 (B) 11.0 22.5 39.0 51.9 56.6 61.5 66.0 (C)  CEEMD 4.93 11.4 18.2 28.4 36.1 45.3 52.2 (A) 4.27 10.1 35.7 43.4 49.4 56.0 65.5 (B) 12.7 27.1 43.3 49.2 53.7 59.8 64.9 (C)Table 2
Ratio of parameter: correlation coefficient.
Methods IMFn Stages 1 2 3 4 5 6 7 Ratio (%)  MCKD-CEEMD 46.0 66.2 83.1 89.6 94.5 96.9 98.4 (A) 31.2 52.2 67.9 80.2 85.7 87.8 93.3 (B) 46.2 66.8 80.2 87.4 91.6 95.8 98.7 (C)  MCKD-EEMD 44.2 64.0 80.5 87.2 92.5 95.4 97.2 (A) 28.8 47.4 62.6 74.7 81.1 84.2 89.8 (B) 44.0 64.0 77.7 85.1 90.2 95.0 98.2 (C)  CEEMD 19.4 31.3 40.6 50.9 62.5 73.7 85.2 (A) 7.82 14.0 23.6 41.5 59.1 78.4 87.9 (B) 13.5 23.8 34.4 48.8 64.2 78.1 88.7 (C)The entropy features ApEn, SE, LLE, and CD values of the first 7 IMFs from different working stages, obtained by CEEMD or EEMD, are extracted. Table3 shows the ApEn values of the 3 stages obtained using the proposed method. According to this idea, 25 one-rotation-cycle time-domain segments of each stage are selected to calculate ApEn, SE, LLE, and CD values, obtaining 3150 (75 signals × 7 IMFs × ApEn/SE/LLE/CD values) values.Table 3
Results of ApEn of the proposed method.
Stages IMFn 1 2 3 4 5 6 7 (A) 0.447 0.251 0.172 0.089 0.058 0.046 0.032 (B) 0.598 0.393 0.288 0.202 0.150 0.113 0.085 (C) 0.771 0.450 0.350 0.173 0.179 0.171 0.134The 3150 values are inputted into SVM for training. In addition, 3150 values of 75 additional signals from 3 different working stages were calculated using the same methods. The overall data collection dates of the 150 segments are shown in Table4. The latter 3150 values are used to verify the classification results of the SVM, which are listed in Table 5, where labels 1, 2, and 3 represent the normal condition, early fault condition, and severe fault condition, respectively. Compared with the 3 control methods, the proposed methods obtain the highest classification rate. Besides, ApEn could be taken as a better characteristic parameter than LLE and CD in fault diagnosis and classification of slewing bearings. With the practical working environment of the slewing bearing and the complexity of the port crane considered, the noise interference in vibration signals is much stronger than the noise interference collected from simulant test platforms. Therefore, although the 100% classification rate is rather difficult to obtain, the proposed method and the results can also provide a great value for reference in practical engineering.Table 4
Data collection dates of the 150 segments.
Stages 1 2 3 4 5 6 7 8 9 Notes (A) 150506 150625 150906 151015 151120 160122 160315 160420 160604 Dates 3 3 2 3 3 3 2 3 3 Segments 1 2 3 3 3 3 2 3 3 3 Segments 2  (B) 160702 160915 161025 161210 170222 170411 170522 170618 - Dates 2 3 3 4 4 3 3 3 - Segments 1 2 2 4 4 4 4 3 2 - Segments 2  (C) 170920 171026 171110 171228 - - - - - Dates 6 6 6 7 - - - - - Segments 1 5 7 8 5 - - - - - Segments 2 7 5   1 segments for feature extraction; 752 segments for condition classification.Table 5
Classification results.
Labels Label 1 Label 2 Label 3 Classification rate (%) Samples 25 25 25 MCKD-CEEMD-ApEn 23 2  93.3 2 23 1 24  MCKD-CEEMD-LLE 23 1 1 89.3 1 22 2 3 22  MCKD-CEEMD-CD 24  1 92.0 1 22 2 2  23  MCKD-CEEMD-SE 23 2  88.0 3 21 1 1 2 22  MCKD-EEMD-ApEn 22 3  89.3 3 22 2 23  CEEMD-ApEn 20 5  85.3 3 21 1 2 23


## 4. Conclusions

This study proposed a combined failure diagnosis method based on MCKD-CEEMD-ApEn for large-size and heavy-load slewing bearings. The original signals were denoised by removing the mean, MCKD algorithm, and CEEMD to obtain the first 7 IMFs, which could maintain the features of original signals and remove part of the strong background noise. In this method, ApEn was adopted as the classification feature and SVM was used to verify the reliable classification rate. Owing to the greater interest drawn by the severe fault stage compared with the early fault stage of the slewing bearings in practical engineering, the methods considered the features of the 3 conditions and conducted method comparison and classification, which could be widely used in a wide range of applications in fault diagnosis and could provide a method to distinguish various working conditions of slewing bearings.



---

*Source: 6321785-2018-04-23.xml*