Abstract

Variational modal decomposition (VMD) has the end effect, which makes it difficult to efficiently obtain fault eigenvalues from rolling bearing fault signals. Inspired by the mirror extension, an improved VMD is proposed. This method combines VMD and mirror extension. The mirror extension is a basic algorithm to inhibit the end effect. A comparison is made with empirical mode decomposition (EMD) for fault diagnosis. Experiments show that the improved VMD outperforms EMD in extracting the fault eigenvalues. The performance of the new algorithm is proven to be effective in real-life mechanical fault diagnosis. Furthermore, in this article, combining with singular value decomposition (SVD), fault eigenvalues are extracted. In this way, fault classification is realized by K-nearest neighbor (KNN). Compared with EMD, the proposed approach has advantages in the recognition rate, which can accurately identify fault types.

1. Introduction

VMD is a novel method to deal with the signals. At present, VMD has shown great achievements in the processing of vibration signal, biological signal, and electrical signal. It has been used for detection, diagnosis, and prediction. In the heart sound signal, Babu et al. [1] adopted VMD to extract the Shannon entropy envelope of the heart beat pulse. The approach can accurately identify the signal features reflecting cardiac abnormalities in ECG. Li et al. [2] proposed an adaptive denoising method based on VMD, which is applied to water supply pipeline leakage location. Compared with EMD, VMD can effectively eliminate modal aliasing and is robust to noise and sampling [3, 4].

At present, many scholars use different methods to study the vibration signal of fault bearing [511]. In the field of fault diagnosis, some scholars study and apply VMD. Zhang et al. [12] constructed a new technique, which applies VMD to mechanical fault diagnosis. Through failure mechanism analysis, they established the fault model in fast Fourier transform (FFT) and envelope analysis. This approach can successfully diagnose the fault of rolling bearing. Jiang et al. [13] presented an adaptive detection method, which uses VMD to detect early defects in bearings. Yi et al. [14] employed VMD to extract bearing fault characteristics. Lv et al. [15] used VMD and multicore vector to realize mechanical fault diagnosis. Compared with the traditional fault diagnosis model, this approach had better performance. Yan et al. [16] adopted VMD to extract hybrid-domain features, identifying fault type with high accuracy. VMD has an important application value for mechanical fault diagnosis [17]. However, the drawback of the abovementioned approaches is that they have the modal aliasing in VMD, which affects the diagnostic performance.

In this paper, a new method is proposed for rolling bearing fault diagnosis. First, a mirror extension to suppress its modal aliasing is used and the improved VMD into bearing fault diagnosis is introduced [18]. Then, combine the improved VMD with SVD to achieve the effective extraction of fault eigenvalues, finally adopt KNN to complete fault classification. The contributions of this paper are presented as follows:(1)An improved VMD approach is proposed based on VMD and mirror extension, which inhibits the end effect. Improved VMD makes it effective for fault diagnosis.(2)Compared with EMD, the proposed approach has an advantage in the recognition rate, which can accurately identify the rolling bearing fault type.(3)The proposed diagnosis framework is effectively verified by detection of the rolling bearing fault. Results demonstrate that our framework has advantages in extracting the characteristic frequency, which is suitable for the detection of manufacturing systems.

The article is organized as follows. The principles of VMD, mirror extension, and KNN are introduced in Section 2. Section 3 provides a new method and the detailed diagnosis scheme of the proposed approach. Section 4 presents the experimental verification. Section 5 concludes this paper.

2. Theories

2.1. VMD

VMD is a new approach of signal processing, which is adaptive and quasi-orthogonal. The original signal is decomposed into several modal components using VMD. Each modal component has limited bandwidth, and it is near the central frequency [19]. As presented in equation (1), the constrained variational models are calculated:where are the narrow-band components and represent the central frequencies. where is the number of decompositions.

In order to find the optimal solution of the constrained variational model, the VMD method uses the quadratic penalty function with good convergence and the Lagrange multiplier operator with strong constraint ability. Therefore, a Lagrangian function can be introduced to optimize the constrained variational problem so as to minimize the narrow-band component and the central frequency obtained in equations (3) and (4). The expression of is introduced:

The formula for minimizing the narrow-band component is expressed as follows:

The formula for minimizing the central frequency narrow-band component is obtained as follows:

2.2. KNN

KNN is a nonparametric prediction algorithm. It searches for k most similar eigenvectors in the historical database to predict the future value [20]. The model has simple structure and high computational efficiency.

The KNN classification algorithm is described as follows. When the testing sample (unknown sample) is given, first search the pattern space to find the K training sample closest to the testing sample, that is, K-nearest neighbor, and then calculate the selected K-nearest neighbor. If a class has the maximum number of the nearest neighbors, the testing sample can be determined to the class. Euclidean distance is used to calculate the distance between the testing sample and all training samples. The formula is defined as follows:where X is the testing sample and Y denotes the training sample. Figure 1 is the detailed steps of KNN.

3. The Proposed Diagnosis Method

3.1. Improved VMD

The mirror extension assumes that a mirror is placed at both ends of the data. The original data sequence image in the mirror is symmetrical with respect to the mirror and the original data sequence [21]. The original data sequence image and the original data sequence in the two mirrors form a continuous and closed ring without endpoints. The internal data completely determines the upper and lower envelopes of the data, so the end effect is fundamentally avoided [22]. Above the mirror surface is the original data. Below the mirror surface is the extension data. After the operation, only the data processing results above the mirror surface are output.

The purpose of boundary extension is not to provide accurate data beyond the end, but to provide a condition that the data within the end completely determine the envelope. The mirror extension satisfies this condition, so it is an ideal extension algorithm.

When VMD analyzes and processes the fault signal, the error is caused by the influence of external factors on both ends [2327]. Taking into account the end effect in VMD, a mirror extension is adopted to suppress it. The newly proposed method has an advantage in extracting the characteristic frequency. Figure 2 is the detailed steps of the improved VMD. The main process of improved VMD can be summarized as follows:  Step 1: initialize , where are the narrow-band components, represent the central frequencies, is the number of decompositions, and represents the Lagrange multiplier.  Step 2: set and execute the whole cycle. Update and for all to reach the preset decomposition number. When , the cycle ends. The updated formula of the narrow-band component and the corresponding central frequency are as follows:  Step 3: update according to the formula:  Step 4: return to Step 2 and repeat the above process until the whole iterative process meets the constraints, and a series of narrow-band eigenmode component signals are obtained. Equation (9) is the constraint condition, where is set to 10−6:

The specific extension process is provided as follows:(1)Find all extreme points of the rolling bearing fault signal . The extreme points include the local maximum points and the local minimum points.The local maximum point sequence of the signal isThe local minimum point sequence of the signal iswhere , are the local maximum points, , are the local minimum points, , are the maximum points, and , are the minimum points.(2)According to equation (12), the symmetrical extension is employed to the extreme points of the fault signal to obtain a new extreme point sequence :(3)Correct the new extreme point sequence after extension:(4)Use the value on the left end of the fault signal as the symmetry plane, and the mirror extension is realized. Meanwhile, the new fault signal of rolling bearing is realized.

3.2. Diagnosis Scheme

Owing to the feature extraction is difficult in the incipient rolling bearing failures, this paper presents a new method combined improved VMD and KNN. The structure of the diagnosis scheme is shown in Figure 3.

The detailed steps of the diagnosis scheme are as follows:(1)Firstly, improved VMD is applied to handle the fault signal, and several narrow-band eigenmode components are realized.(2)Then, the narrow-band components are selected, which have rich fault information. These modal components are selected to reconstruct the fault signal, and the envelope spectrum is conducted for it.(3)Finally, the feature vector from modal components is extracted by using SVD. After feature extraction, there are 100 data sets. 70 of the 100 data sets are randomly employed to train the KNN classification model, and the rest are employed to test the classification performance.

4. Experimental Verification

4.1. Experimental Setting

As shown in Figure 4, the fault diagnosis platform consists of a 0.75 kW three-phase asynchronous motor, two couplings, a reducer, a magnetic powder brake, a piezoelectric accelerometer, a faulty rolling bearing, and a photoelectric speed sensor. The faulty rolling bearing is installed in the bearing pedestal of the reducer.

The signal acquisition system is used to collect the fault signals, in which the rotation speed and the sampling frequency are set at 600 r/min and 1 kHz, respectively. Table 1 is the basic parameters of rolling bearing. According to these parameters, the fault characteristic frequencies are calculated: inner race is 99 Hz and outer race is 71 Hz. The improved VMD is applied to decompose the fault signal. When the mode number K is different, their central frequency is different. The relationship between them is depicted in Figure 5.

4.2. Experimental Results and Analysis

When the value of K starts from 5, the central frequency is close [26], see Figure 5(a). This is an over decomposition phenomenon. Hence, the K value taken in the test is 5. Based on VMD experience, the balance parameter constrained by data fidelity adopts the default value of 2000, and the time step of the double rise is 0.1. Figure 6(a) illustrates the time domain diagram of the inner race fault, and its improved VMD result is demonstrated in Figure 7(a).

As illustrated in Figure 5(b), when the value of K starts from 5, the central frequency is close, which is an over decomposition phenomenon [28]. Thereby, the K value taken in the test is 5. Based on VMD experience, the balance parameter constrained by data fidelity adopts the default value of 2000, and the time step of the double rise is 0.1. Figure 6(b) shows the time domain diagram of the outer race fault. The improved VMD result is demonstrated in Figure 7(b).

As shown in Figure 6(a), the original signal of the inner race fault contains considerable background noise, which obscures the fault information. As shown in Figure 7(a), compared with the original signal, the modal components after the improved VMD eliminate the noise, which is the function of Wiener filter. There are obvious fault shock components in U3, U4, and U5. Therefore, these three modal components are selected to reconstruct the signal, and the envelope spectrum is conducted for it. The result of the envelope spectrum is depicted in Figure 8(a).

In contrast to the improved VMD, the first five components of EMD are also taken, and Figure 9(a) is the decomposition results. In Figure 9(a), IMF1, IMF2, and IMF3 contain many fault feature information, and they are reconstructed. Figure 8(b) presents the envelope spectrum of the reconstructed signal.

From Figure 6(b), the original signal time domain of the outer race fault contains large background noise, which obscures the fault information. Figure 7(b) shows that the modal components eliminate the noise, which is the function of the Wiener filter. There are obvious fault shock components in U2, U4, and U5. Hence, these three modal components are selected to reconstruct the signal, and the envelope spectrum is conducted for it. The result of the envelope spectrum is provided in Figure 10(a).

As a comparative test analysis, the first five components of EMD are taken, and Figure 9(b) is the decomposition result. As evident in Figure 9(b), IMF1, IMF2, and IMF3 contain many fault feature information, and they are reconstructed. Figure 10(b) presents the envelope spectrum of the reconstructed signal.

Figure 8(a) displays that the envelope spectrum has obvious fault shock characteristics around 99 Hz, which is basically consistent with the theoretical value. The small gap does not affect the fault identification. 198 Hz is the second frequency, and 297 Hz is the third frequency. As shown in Figure 8(b), the envelope spectrum is not ideal. 99 Hz and its octave are not demodulated. Consequently, the decomposition result of the improved VMD outperforms EMD.

Figure 10(a) shows that the envelope spectrum has obvious fault shock characteristics around 71 Hz. This is basically consistent with the theoretical value. The small gap does not affect the fault identification. 142 Hz is the second frequency, and 213 Hz is the third frequency. In Figure 10(b), the envelope spectrum is not ideal. 71 Hz and its octave are not demodulated. Therefore, the decomposition result of the improved VMD is shown to outperform EMD.

In the fault classification, K in the KNN algorithm is set to 1. The first three groups of singular values are selected as the nearest neighbor classification training sets and testing sets. For each fault type, 70 of the 100 datasets are randomly employed to train the KNN classification model, and the rest are applied to test the classification performance. The fault classification results after improved VMD + SVD are shown in Figure 11(a), and the fault classification results after EMD + SVD are provided in Figure 11(b). Figures 11(a) and 11(b) show that the performance of the improved VMD + SVD is better than EMD + SVD.

The fault recognition rates are reported in Table 2. In the classification results after the improved VMD + SVD, one inner race fault is incorrectly identified as the outer race fault, and one outer race fault is incorrectly identified as the inner race fault. In the classification results after EMD + SVD, three inner race faults are incorrectly identified as the outer race faults, and five outer race faults are incorrectly identified as the inner race faults.

5. Conclusion

A new method combining the improved VMD and KNN for the rolling bearing fault detection system is presented. First, the improved VMD decomposes the fault signal to obtain the modal component. Then, the modal components with many fault shock components are reconstructed. Finally, the envelope spectrum is applied to the reconstructed signal, and the characteristic frequency is identified successfully. Furthermore, combining with SVD, KNN realizes fault type classification. The experimental results demonstrate that the proposed approach has potential application in rolling bearing fault diagnosis.

Data Availability

Some or all data, models, or code generated or used during the study are available from the corresponding author upon request.

Conflicts of Interest

The authors have no conflicts of interest/competing interest to disclose.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 41572347).