Abstract

Adaptive sparsest narrow-band decomposition (ASNBD) method is proposed based on matching pursuit (MP) and empirical mode decomposition (EMD). ASNBD obtains the local narrow-band (LNB) components during the optimization process. Firstly, an optimal filter is designed. The parameter vector in the filter is obtained during optimization. The optimized objective function is a regulated singular local linear operator so that each obtained component is limited to be a LNB signal. Afterward, a component is generated by filtering the original signal with the optimized filter. Compared with MP, ASNBD is superior in both the physical meaning and the adaptivity. Drawbacks in EMD such as end effect and mode mixing are reduced in the proposed method because the application of interpolation function is not required. To achieve the fault diagnosis of roller bearings, raw signals are decomposed by ASNBD at first. Then, appropriate features of the decomposed results are chosen by applying distance evaluation technique (DET). Afterward, different faults are recognized by utilizing maximum margin classification based on flexible convex hulls (MMC-FCH). Comparisons between EMD and ASNBD show that the proposed method performs better in the antinoise performance, accuracy, orthogonality, and extracting the fault features of roller bearings.

1. Introduction

Analysis methods for vibration signals play an important role in mechanical fault diagnosis. Techniques like matching pursuit (MP) [13] and empirical mode decomposition (EMD) [4, 5] have been widely used for analyzing vibration signals of mechanical components recently.

EMD method, which was proposed by Huang [4], is a suitable method for nonstationary signal analysis. As the extremas of an original signal is calculated, the signal is divided into several intrinsic mode functions (IMFs) with physical meaning when EMD is applied. EMD has been proved to be effective in the fields such as fault diagnosis [6], condition monitoring [7], and image identification [8]. Nevertheless, some disadvantages including mode mixture [9, 10] and end effect [11, 12] are existed in EMD as interpolation function is involved in the calculation process.

MP was brought forward by Mallat [1] in 1993. Various dictionaries are utilized to gain the sparsest decomposition results of an original signal. MP has been used in many areas since it was proposed [1317]. Nevertheless, the highly redundant dictionary must be selected beforehand, and the physical meaning of the obtained decomposition results are not explicit enough when MP is used; thus, the application range of MP is limited [3].

An adaptive time-frequency approach for nonstationary signal analysis, namely, adaptive sparsest narrow-band decomposition (ASNBD) is proposed based on the basic idea of MP and EMD [18]. ASNBD decomposes an original signal into several local narrow-band (LNB) signals [19]. An optimal filter must be determined at first in ASNBD method. The filter parameters are defined during the optimization process. The optimized objective function is a regulated operator [19] so that each obtained component is limited to be a LNB signal. Afterward, a component is generated by utilizing the optimized filter. The obtained components are constrained to be LNB signals as a regulated differential operator is applied as the objective function.

The similarity between ASNBD and EMD is that both the two methods can generate a series of meaningful component. However, to gain the decomposition results, an optimization method is applied in ASNBD instead of calculating the envelopes of the extreme points in EMD. Compared to other improved EMD methods, as the calculation of extreme points is inevitable [10], the accompanied disadvantages such as mode mixing and end effect are alleviated.

An optimization problem has to be solved to gain the sparsest representation in both MP and ASNBD. Nevertheless, the construction of a highly redundant dictionary is no more needed in ASNBD before the method is applied. Meanwhile, the generated components are LNB signals with physical meaning. Therefore, compared to MP, ASNBD is better at either the adaptivity or the physical meaning. In short, ASNBD is inspired by the basic theory of EMD and MP and proposed by avoiding their disadvantages and learning their advantages.

The maximum margin classification based on flexible convex hulls (MMC-FCH), which was inspired by SVM, was proposed by Zeng [20, 21]. The class regions of the sample sets are approximated by convex hulls of the training samples in SVM. The convex hull of SVM is an unrealistic approximation of the class region of the samples when the samples are not enough. MMC-FCH was proposed aiming at solving the problems accompanied with the application of convex hull by extending the flexible convex hulls. As MMC-FCH was used to recognize rolling element faults and showed a satisfactory performance in reference [20], it is chosen to fulfill the fault diagnosis in the paper.

Firstly, the raw signals are decomposed by using ASNBD and the features of the components are generated. Then, the dimensions of the features are reduced by distance evaluation technique (DET) [22, 23] to improve classification accuracy. Finally, MMC-FCH is applied for classifying the fault type of the bearing.

The paper is organized as follows. ASNBD is proposed in section 2. In section 3, descriptions about feature selection are given. In section 4, the MMC-FCH algorithm is given out. In section 5, the procedure of the proposed method is stated, and the effectiveness of is verified by experimental and simulation signals. The conclusions are given in the end.

2. ASNBD

2.1. Vibration Model

Vibration impulses appear when a fault occurred on the roller bearing [24], and the collected vibration signals are modulated by the corresponding impulses. Therefore, the signal model can be described as [25] is the impulse signal, is the fault signal, denotes the periodic modulation, and is the noise disturbance. Besides, discrete time and continuous time are denoted as “n” and “t” in this paper, respectively.

It can be inferred from Fourier transform (FT) that is consist of several cosine signals. The vibration signal is a relative narrow-band signal. Hence, the model of the vibration signal is as follows: is the modulated signal and is a narrow-band signal. The definition of LNB signal is stated in the following section.

2.2. Definitions of LNB and INBC

For a given signal, if any neighborhood range of the data points can be regarded as a narrow-band signal, the given signal is a local narrow-band (LNB) signal [19].

If a generated component satisfies the conditions of LNB signal, then it can be regarded as an intrinsic narrow-band component (INBC).

2.3. Singular Local Linear Operator (SLLO)

For an operator T, if there is a neighborhood range , when , and 0; otherwise, T is a SLLO [19].

A singular local linear operator T proposed in reference [19] is used by the ASNBD method.

2.4. ASNBD Procedure

For an original signal , the ASNBD algorithm is as follows:(1)Set and , .(2)Construct the optimal problem :T is given in equation (4). Minimizing indicates that is a LNB signal. is the differentiation of . is applied for regulating the residue. is the weight parameter, which is set to 17 in the experiments of this paper. It is worth mentioning that in the optimal problem , all the data in is optimized, which needs a big amount of calculation.(3)Set .(4)If , stop. If not go to step 2 and set .

To reduce the calculated amount in step (2), the optimization problem is replaced with the optimization procedure as below, and the genetic algorithm (GA) is applied to solve the optimal problem. The approach for solving the nonlinear optimization problem is as follows. Meanwhile, Figure 1 shows the calculation process of ASNBD.(1)The FFT of is calculated as .(2)The filter is as follows (see also Figure 2):(3)The optimal problem as below is solved by GA:The optimal is gained during the optimization process. equals the inverse fast Fourier transformation (IFFT) of .(4)Set to be .

2.5. Comparison between ASNBD and EMD

To verify the validity of ASNBD, a simulation signal is analyzed: is consist of an AM-FM signal and an attenuating signal. Figures 3(a)3(c) are the waveforms of , the AM-FM component, and the attenuating component, respectively.

is decomposed by EMD and ASNBD (mirror extension [26] is used to reduce end effect in this paper). Figures 4 and 5 show the analysis results. The INBCs are basically the same with the real signals, the results gained by EMD have obviously mixed modes.

To compare the decomposing results, the parameters such as the index of orthogonality , the correlation coefficient , the rate of energy error , and time T are considered, which are defined in equations (9)–(11). and are used to evaluate the effectiveness of ASNBD. and T are applied to compare the orthogonality and the computational efficiency of ASNBD and EMD, respectively.

Table 1 shows the calculated parameters. of the components generated by EMD is much bigger than that obtained by ASNBD. And the INBCs are more close to the actual signals compared to EMD. Thus, the results obtained by ASNBD are more accurate than EMD. Meanwhile, the IO of ASNBD is smaller which shows the orthogonality is better. Nevertheless, ASNBD still spends more time, which is one drawback in ASNBD.

For comparing the antinoise performance, x2(n) composed of an intermittent noised signal with different signal-to-noise ratios (SNRs) is applied:

SNR is obtained as follows:

Signal x2(n) (SNR equals 5, 10, or 15 dB) and the components are shown in Figure 6. The results generated by EMD and ASNBD are given in Figures 79. It can be seen from the figures that there exists severe mode mixing of the components obtained by EMD (labeled by the red line) because of the noise interference. ASNBD can restrain mode mixing effectively. The INBC1 generated by ASNBD is in excellent accordance with the AM-FM component of .

In Table 2, the calculated parameters are given. The parameters indicate that ASNBD method excels EMD method both in possessing better orthogonality and the decomposition accuracy when the SNR varies.

Therefore, at least for the signals shown in Figure 6, ASNBD excels EMD both in orthogonality and accuracy, which means that ASNBD is more effective than EMD when utilized to the simulation signals with different SNRs.

The simulation analysis results indicate that ASNBD excels EMD in the accuracy, the orthogonality, and the antinoise performance when applied to the simulation signals above.

3. Feature Selection

After every raw data set is decomposed into several components, different features of the data sets can be obtained by calculating various statistical parameters. For providing enough feature information for selection, some common time and frequency domain statistical parameters shown in Table 3 are applied to extract the fault feature.

It is noted that not all features could positively affect the identification of bearing faults. Thus, useful features which are usually called as salient features have to be chosen. DET [22, 23] is applied to get salient features, and the procedure of DET is illustrated as follows:(1)The mean distance is calculated. is the jth feature of the kth sample with the ith condition. C, J, and Ni are the number of conditions, features, and samples:then the mean distance of all conditions is obtained:(2)The mean value of the features is obtained:then the average distance of various conditions is gained:(3)The effectiveness factor is obtained:

In the end, the features with larger are selected to be the salient features.

4. MMC-FCH

Inspired by the basic theory of SVM and flexible convex hull, MMC-FCH method was proposed by Zeng [20, 21]. The main idea of MMC-FCH is replacing the convex hull in SVM with the flexible convex hull, while the optimized hyper-plane is obtained by the same approach as SVM.

A convex hull in SVM is considered to be a convex combination of a data set as all the combination parameters amount to one and are nonnegative. For , the convex hull is defined as follows:

Normally, the tightest convex approximation of the class region of the data is gained by using the convex hull. Nevertheless, for the classes with complicated convex forms, the convex hull is not accurate enough as the class region always extends beyond it. It can be inferred that convex hull should be made looser to fit more complicated situations. Therefore, reference [20] used a new model which is called flexible convex hull instead of the convex hull applied in SVM. The improvement of the flexible convex hull is that it uses a flexibility factor to construct the bounds of the coefficients. Flexible convex hull is defined as follows:

For a prescribed , a data point expands along with . is the average of X. Then the expanded is obtained by calculation to construct the flexible convex hull . The bounds of coefficients can be different by applying the flexible convex hull. It can be easily inferred that flexible convex hull becomes convex hull when ; thus, convex hull is a special situation of the flexible convex hull. If , then the flexible convex hull becomes looser than the convex hull. If the flexible factor is set appropriately, it is possible that a nice approximation of the class region can be obtained.

The main difference between MMC-FCH and SVM is that MMC-FCH uses flexible convex hull instead of convex hull to get a more satisfactory matching of the class region. The method of obtaining the optimal hyper-plane of the two methods is similar. More details of MMC-FCH can be found in references [20, 21]. Results in reference [20] have already shown that MMC-FCH performed better than SVM when applied to the recognition of roller bearing fault type. Hence, no more comparison is made between MMC-FCH and SVM in this paper as it is unnecessary.

5. The Proposed Method and Experimental Analysis

5.1. The Proposed Method

By applying ASNBD, DET, and MMC-FCH, an approach for mechanical fault diagnosis is proposed. Figure 10 shows the flow chart of the proposed method. The main calculation procedure is as follows:(1)Decompose the experimental signals by using ASNBD, suppose every data set is decomposed into M INBCs(2)Calculate the features (shown in Table 3) of the training signals and the INBCs(3)The normalized effectiveness factors are calculated; then salient features are selected by applying DET(4)Classify the fault type by training and testing with MMC-FCH using the salient features

5.2. Experimental Analysis

The experimental data are offered by Case Western Reserve University [2731]. The experimental system is shown in Figure 11. The type of all the tested bearings is 6205-2RS JEM SKF. Several single-point faults were separately set to the tested roller bearings to generate three fault types including inner race, outer race, and ball fault (the location of outer race faults is 6 o’clock). The fault diameters were set to be 0.021 in., 0.014 in., and 0.007 in., which implied that there were three categories of the degrees of the roller bearing faults: severe, moderate, and slight. The vibration data sets which were collected at the motor housing’s drive end are used. The sampling frequency, rotating speed, and the number of the points of each sample were 12 kHz, 1797 rpm, and 4000, respectively. The details about these bearing data sets are given in Table 4. Data sets in Table 4 cover ten conditions and each condition contains 30 samples.

5.2.1. Experiment 1

In experiment 1, the ten classes of fault type are considered to be classified by MMC-FCH. Firstly, to compare the performance of EMD and ASNBD on extracting useful features from the original vibration signals, all the 300 data sets are decomposed into eight components by ASNBD and EMD, respectively. The examples of the vibration data samples from each class are shown in Figure 12. The decomposition results of a data sample from class 2 by applying ASNBD and EMD are given in Figures 13 and 14. Then the features from Table 3 of the original signals and the components are calculated. DET is used for choosing salient features, the normalized effectiveness factors of the features which are larger than 0.5 are chosen, and the results are shown in Figures 15 and 16. To assess the proposed method when the number of training data samples varies, 1, 2, 3, 4, 5, 10, 15, 20, and 25 signals of every class were randomly chosen to construct the training samples as well as the rest samples are utilized for testing. The whole test process was repeated for 20 times so that the influence of randomness on the results of classifications can be alleviated. The final classification accuracy in Tables 5 and 6 are the mean values of the 20 results. The flexible factor is set to 1.82 in this experiment.

There are 20 salient features when ASNBD is applied and 15 salient features when EMD is applied, which indicates that ASNBD can extract more effective features than EMD at least in this experiment. Results in Table 5 show that ASNBD and EMD both can generate satisfying classification accuracy. The accuracy is above 90% even when there is only one training sample. However, the classification accuracies of ASNBD are always higher than EMD in this experiment. Thus, ASNBD is more effective in the analysis of experiment 1 when the fault locations and fault severity varies.

5.2.2. Experiment 2

The severity of the faults is ignored and only the locations of the faults are identified in this experiment. Thus, only 4 faults consisting of ball fault, regular condition, inner fault, and outer fault are considered. As there are 30 samples from the regular condition and 90 samples from every other condition, 30 samples are randomly chosen from the 90 samples of inner fault or outer fault. The sample partitions of this experiment are just like those of experiment 1. As all the data sets are decomposed into eight components by ASNBD and EMD, there are still 198 statistical parameters for testing the effectiveness of both ASNBD and EMD. In this experiment, the normalized effectiveness factors of the features which are bigger than 0.8 are chosen. Salient features of the components obtained by EMD and ASNBD in experiment 2 are shown in Figures 17 and 18. The flexible factor is set to 2.59 in this experiment.

There are 15 salient features when ASNBD is applied and 12 salient features when EMD is applied; thus, ASNBD still extracts more effective features than EMD in experiment 2. Results in Table 6 show that the accuracies are above 90% when the training samples are more than 5, while the classification accuracy drops rapidly when the number of the samples is very few (less than or equal to 5). This is because that the classification has become more difficult when the severity of the faults are ignored, since the samples with the same fault locations still have intrinsic differences when the fault severity varies. Even so, ASNBD still outperforms EMD.

Experiment 1 is conducted for analyzing the effectiveness of the proposed method when applying to recognize rolling bearing with different fault types, and experiment 2 is conducted for recognizing the fault locations. Experiments 1 and 2 show that ASNBD and EMD both show a satisfying performance on identifying different fault types of roller bearings when combined with MMC-FCH. Moreover, ASNBD is superior to EMD, either the samples can be easily classified (experiment 1) or the recognition task is relatively difficult (experiment 2). However, the mean time of the identification process of every data set shown in Tables 5 and 6 shows that ASNBD still costs more calculation than EMD when combined with MMC-FCH.

6. Conclusions

ASNBD is proposed by applying the idea of EMD and MP for nonstationary signal analysis. The theoretical analysis indicates that the proposed method excels MP in both the physical meaning and the adaptivity. And problems like end effect and mode mixing of EMD are reduced because the application of interpolation function is not required. The simulation analysis results show that ASNBD excels EMD in the accuracy of the decomposition results, the orthogonality and the antinoise performance. Furthermore, ASNBD is applied to classify the roller bearing fault. The vibration signals of roller bearings with fault are decomposed by ASNBD at first. Then DET is utilized to choose salient features of the generated components. To fulfill the fault diagnosis, the new classification method MMC-FCH is introduced to identify the fault types of roller bearings automatically. The experimental analysis indicates that ASNBD not only shows a satisfying performance on identifying different fault types of roller bearings when combined with MMC-FCH, but also is superior to EMD on extracting accurate fault features.

Data Availability

Previously reported .mat data were used to support this study and are available at DOI: 10.1016/j.ymssp.2015.04.021. These prior studies (and data sets) are cited at relevant places within the text as references [23].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Key Research and Development Program of China (2018YFF0212902), National Natural Science Foundation of China (51805161, 51475159, and 51575178), Hunan Provincial Key Research and Development Program (2018GK2044), and Hunan Provincial Natural Science Foundation of China (2018JJ3187, 2017JJ1015, 2017JJ2086, and 2018JJ4084).