Abstract

In view of the incipient fault characteristics are difficult to be extracted from the raw bearing fault signals, an incipient bearing fault diagnosis method based on parameter-adaptive variational mode decomposition (VMD) is proposed. The beetle antennae search (BAS) algorithm is adopted to seek for the optimal combination of the VMD parameters. The reciprocals of the calculated kurtosis values of intrinsic mode functions (IMFs) decomposed via VMD are employed as a fitness function in the searching process. The optimal mode number and the quadratic penalty term of VMD are adaptively set after the search. Afterwards, a vibration signal is decomposed into a set of IMFs using the parameter-adaptive VMD, and the IMF with the maximal kurtosis value is selected as the sensitive one. The selected IMF is further analyzed by Hilbert envelope demodulation. The resulting envelope spectrum can show the significant fault impulse characteristics which are highly helpful to diagnose incipient bearing faults. The kurtosis and the proportion of fault energy are introduced as the input vector of the extreme learning machine (ELM). Comparisons have been conducted via ELM to evaluate the performance by using EMD and the fixed-parameter VMD. The experimental results demonstrate that the proposed method is more effective in extracting the incipient bearing fault characteristics.

1. Introduction

Bearing is one of the most critical components in rotating machinery. The presence of defects in the bearing may lead to noise, vibration, or even system breakdown. For this reason, the bearing fault diagnosis has received considerable attention in the past decades.

It has become a widely acceptable method to assess the bearing performance by studying the vibration signals. As a result, this method has been adopted in many fields, such as bearing fault diagnosis [13], performance evaluation [4, 5], and equipment health monitoring [68].

However, the bearing vibration signals are ordinarily nonstationary and may be polluted by noise interference [6, 9]. Additionally, the measured fault characteristic frequencies usually could not exactly match the theoretically values due to geometry imperfection, speed variation, or load variation [10], and so on. All of these factors perplex the extraction of fault characteristics from the raw signals.

Bearing fault diagnosis is mainly applied in two fields which are the quality control in bearing manufactories and the condition monitoring of the bearing in service, respectively. This paper is based on the former application and focuses on the fault diagnosis techniques. In this case, the bearing faults are mostly characterized as incipient faults which mean that the faults are slight or tiny, and the characteristics are not obvious. It gets difficult to accurately recognize and diagnose the incipient faults. At the same time, the measurements suffer little negative influence introduced by other mechanical transmission parts such as gears. Bearing manufactories usually complete bearing vibration measurement on the dedicated vibration tester. The measurement conditions are that the outer ring of the bearing is stationary, the inner ring rotates together with the shaft at the specified speed, and the specified load is applied at the same time.

Empirical mode decomposition (EMD) [11, 12] is a well-known adaptive time-frequency processing method and can recursively decompose signals without the prior knowledge. A large number of scholarly articles covering a variety of applications have been published. However, EMD suffers from some shortcomings, such as the limited mathematical understanding, mode mixing, pseudocomponents, and end effects [9].

Variational mode decomposition (VMD) is a novel and self-adaptive signal decomposition algorithm proposed by Dragomiretskiy and Zosso in 2014 [13]. In comparison with the recursive EMD method, VMD can decompose a signal into an ensemble of band-limited intrinsic mode functions (IMFs) synchronously. The literature [1315] compared the performance of VMD and EMD and deduced that VMD outperformed EMD with regard to the noise robustness and the effectiveness of feature extraction.

However, a significant drawback of VMD is that the parameters, namely, the mode number M and the quadratic penalty term α, need to be set in advance. Unreasonable preset parameters may induce the loss of effective modes or the mixing of different components [16] and affect the subsequent feature extraction results. M should be specified based on the number of different frequency components contained in the raw signals, while α, namely, the mode frequency bandwidth control parameter, is determined based on each central frequency of the mode. They are correlated with each other. A large number of modes may result in redundant VMD information. Correspondingly, a small number of modes may lead to mode mixing in the VMD results. α is related to the performance of suppressing noise interference. As α decreases, the bandwidth of the mode tends to be wide, and the VMD results may involve more background noise. But, the VMD results are likely to be distorted if the bandwidth is too narrow [14]. Since the two parameters need to be specified in advance without any prior knowledge about the raw signals, it is difficult to guarantee the accuracy of the VMD results. Accordingly, seeking for the optimum of M and α that matches with the analyzed signals is the key to the VMD method.

In order to overcome the critical drawback, some studies have been conducted. Lian et al. [17] chose M according to a series of indicators including permutation entropy, extreme value in the frequency domain, kurtosis criterion, and energy loss coefficient. Liu et al. [18] estimated M by using the minimum redundancy maximum relevance. The above two methods in [17, 18] both optimized the value of M but neglected the influence of α on decomposition results. Based on the grasshopper optimization algorithm, Zhang et al. [19] proposed a method to obtain the optimum of M and α. Wang et al. [20] used multi-objective particle swarm optimization (MOPSO) to obtain the optimum of M and α, in which the symbol dynamic entropy and the power spectral entropy were selected as the objective function of MOPSO. In [16], an improved adaptive genetic algorithm was proposed to optimize the two parameters of VMD. The performance of several decomposition methods for the simulation signal was compared, and authors considered that the decomposition effect of VMD preceded those of other methods. In [19], the effective value range of α was preassigned in the interval of [1000, 10,000]. However, with separate optimization algorithm, our experimental results described in the following as well as the literature [16] indicate that the optimal α may lie out of the interval stated in [19], respectively.

Compared with the above key parameters, the noise tolerance and the convergence tolerance level have little influence on the decomposition results; thus, the default values in the original VMD method are usually adopted [19].

In this paper, a novel optimization algorithm called beetle antennae search (BAS) is employed to estimate the optimal M and α. BAS is a nature-inspired algorithm developed by Jiang and Li [21]. In [21], the global optimization performance of BAS was benchmarked on the Michalewicz function and the Goldstein-Price function, in which the numerical results validated the efficacy of this algorithm.

2. BAS Principles

BAS is a nature-inspired algorithm to solve the optimization problems, which mimics the detecting and searching behaviors of long-horn beetles. A beetle wobbles its two antennae to detect the odour while preying or finding mates, i.e., the beetle explores nearby area randomly using a pair of antennae. When the antenna in one side receives a higher concentration of odour, the beetle would turn to the direction towards the same side; otherwise, it would turn to the other side. Searching behavior of beetles may be formulated in a way which is associated with an objective function to be optimized.

The position of the beetle can be expressed as a vector at t-th time instant (). is denoted as a fitness function which indicates the concentration of odour at position . The maximal value of corresponds to the source point of the odour, namely, the optimal solution of the function.

The random direction of searching behavior can be modeled as follows:where presents a random function and c denotes the dimensions of the position. The barycenter coordinates of both right-hand and left-hand antennae can be then described in the following equation:where and stand for the barycenter coordinate of the right-hand antenna and that of the left-hand antenna, respectively, and s is the sensing length. The value of s should be large enough to cover an appropriate searching area for fear of falling into local minimum points at the beginning and then attenuate with the iterations.

Considering the searching behavior, the iterative model of detecting behavior can be formulated as follows:where represents a sign function and is the step size of searching. In consideration of convergence speed, follows a decreasing function of t. The initialization of should be equivalent to the sensing length s. As examples, the sensing length and the step size may be updated as follows:

As aforementioned, the BAS algorithm can achieve efficient optimization without the acquisition of functional expression or the gradient of function. Compared with the particle swarm optimization algorithm [22, 23], the BAS algorithm only employs one individual, i.e., one beetle; thus, the computation is reduced.

3. Parameter-Adaptive VMD Method Based on BAS

In this section, a parameter-adaptive VMD method based on BAS for incipient bearing fault diagnosis is introduced. The basic idea of the proposed method is to seek for the optimal M and α by using the BAS algorithm.

Kurtosis criterion [2426] has been developed based on the resonance demodulation technique. As a dimensionless index, the kurtosis can be used as a measure of impact components in the signals. The impact impulses caused by defect will increase the kurtosis obviously. Dong [27] introduced the spectral Lp/Lq norm and deduced that the kurtosis could be explained as a special case of the spectral Lp/Lq norm. Dong drew a conclusion that the spectral Lp/Lq norm could be used for characterizing the repetitive impulses in the bearing vibration signals.

The bearing faults, in view of the application described in this paper, are mostly characterized by surface defects such as pitting and spot, which account for more than 90% of all faults. In our experiments, the kurtosis is effective for the diagnosis of bearing faults and can be used to identify the conditions of the bearings. In consideration of this advantage, the reciprocal of kurtosis values calculated from modes may be adopted as the fitness function to optimize VMD parameters. The optimization objective is to search for the mode with the maximal kurtosis, which indicates the most obvious fault information is contained in the screening mode.

After obtaining the optimal parameters of VMD, the vibration signals are decomposed into a set of IMFs via VMD, and the mode with the maximal kurtosis is selected as the sensitive one. The selected mode may be further analyzed by Hilbert envelope demodulation to estimate whether the bearing is defective.

The following steps detail the procedures of the proposed fault diagnosis method:Step 1: calculate the fault characteristic frequencies according to the geometry of the bearing and its rotational speed.Step 2: in order to optimize M and α, initialize the vector . In this paper, M is preset as an integer in the interval of [3, 8], and the initial values of and are equal to 3 and 2000, respectively.Step 3: with each , decompose the vibration signals into a set of IMFs via VMD. Compute the kurtosis values of all IMFs, and specify the reciprocal of the maximal kurtosis as the value of .Step 4: model the normalized random searching direction according to equation (1). Calculate and with the aid of equation (2), where and .Step 5: decompose the vibration signals via VMD by using and , separately. And then calculate the respective and .Step 6: update by solving equation (3).Step 7: reiterate Steps 4–6 through 100 loops to acquire the optimal M and α.Step 8: decompose the vibration signals into a set of IMF components via VMD by using the optimal M and α, and choose the IMF with the maximal kurtosis.

The flowchart of the proposed method is illustrated in Figure 1.

The proposed method is applied to the low-noise deep groove ball bearing 6203. The geometry of bearing 6203 is listed in Table 1. More details about the calculation of the fault characteristic frequencies can be found in [28]. The measurement conditions are in accordance with those described in Section 1.

4. Applications on Bearing Fault Diagnosis

The experimental setup is a self-developed bearing vibration tester, as shown in Figure 2. Vibration signals are sampled via two vibration velocity sensors. A data acquisition card PCI-6143 is used to sample the conditioned vibration signals. The sampling rate and the sampling numbers are set to 45 KHz and 8192, respectively. The spindle speed is 1800 r/min.

4.1. Incipient Inner Raceway Fault

Figure 3 shows a time-domain vibration signal of bearing 6203 with the incipient inner raceway fault. As shown, it is difficult to find periodic impact impulses in the time-domain waveform.

The above vibration signal is processed by the proposed method. and in the VMD method are both set to the default values, i.e., and . The optimal M and α acquired by BAS are 8 and 380, respectively. Decomposition results of the signal via VMD with the optimal parameters are presented in Figure 4. The kurtosis values (K) of IMFs are calculated and listed in Table 2. The maximal kurtosis appears in IMF3, namely, the IMF3 has the most observable impact composition.

As shown in Figure 4, despite IMF3 is still polluted by a little noise, significant time-domain impulses can be found. For further evaluating the fault characteristic frequency of the inner raceway, IMF3 is demodulated by Hilbert envelope analysis. The consequent envelope spectrum is shown in Figure 5.

It can be observed that after processed by the proposed method, 5 spectrum peaks are highlighted in the envelope spectrum of IMF3. These 5 spectrum peaks correspond to the fault characteristic frequency fi (148.3 Hz) and its harmonics. The evident features suggest the inner raceway of the bearing is defective.

For the sake of comparison, the signal shown in Figure 3 is decomposed by means of EMD. The results are illustrated in Figure 6. As shown, faint time-domain impulses can be found in mode4 and mode3. The kurtosis values of the modes are calculated and listed in Table 3.

Similarly, mode4 with the maximal kurtosis value is chosen and demodulated by Hilbert envelope analysis. The corresponding envelope spectrum is shown in Figure 7.

In order to further validate the effectiveness of the proposed method in extracting the bearing fault characteristics, the fixed-parameter VMD method is adopted to process the signal shown in Figure 3. Drawing on the experiences in the literature [15], M is set to 4, while α, , and are all set to the default values in original VMD, i.e., , , and . Likewise, decomposition results of the same signal are illustrated in Figure 8, and the kurtosis values of all IMFs are listed in Table 4.

Following the same ideas, IMF2 with the maximal kurtosis value is chosen and demodulated by Hilbert envelope analysis. The result is presented in Figure 9.

Adopting the proposed method, the maximal kurtosis value of the decomposed IMF listed in Table 2 is 11.06, which is much greater than maximal kurtosis values listed in Tables 3 and 4. The comparison implies that the sensitive component decomposed by the proposed method shows more obvious fault characteristics. As shown in Figure 4, IMF3 selected by the proposed method demonstrates distinct impact impulses in time-domain waveform. In contrast, mode4 illustrated in Figure 6 and IMF2 illustrated in Figure 8 both exhibit more noise and less obvious impact impulses. Comparing Figures 5, 7, and 9, it can be observed that, in Figure 5, the fault characteristic frequency of inner raceway fi and its harmonics are highlighted in the envelope spectrum. A fundamental tone at 148.3 Hz can be clearly found in both Figure 7 and 9, but the higher harmonics like the third and the fourth components are smeared or completely obscured by noise in two envelope spectra.

The above experiments demonstrate that the proposed parameter-adaptive VMD outperforms EMD and the fixed-parameter VMD with regard to noise robustness and effectiveness in extracting weak fault characteristics contained in bearing vibration signals. Meanwhile, the experiments also indicate the VMD parameters have appreciably impacts on decomposition results. The performance of VMD with inappropriate parameters may be inferior to that of EMD.

4.2. Incipient Outer Raceway Fault

Figure 10 illustrates a time-domain vibration signal with the incipient outer raceway fault. The optimal M and α sought by BAS are 5 and 130, respectively. Decomposition results of the signal via VMD with the optimal parameters are displayed in Figure 11. The kurtosis values of the resulting IMFs are calculated and listed in Table 5.

As shown in Figure 11, distinct time-domain impulses are exhibited in IMF2, and meanwhile, the maximal kurtosis value appears in IMF2. In view of the aforementioned analysis, IMF2 is separately selected for further envelope analysis via Hilbert transform. The result is displayed in Figure 12.

The fault characteristic frequency of outer raceway fo (93.4 Hz) and its harmonics have been clearly demonstrated in Figure 12. The distinct spectrum peaks indicate the outer raceway of the bearing is defective.

For comparison, the signal shown in Figure 10 is processed via EMD. Decomposed mode4 with the maximal kurtosis value 10.15 and its envelope spectrum are separately displayed in Figures 13 and 14.

The signal shown in Figure 10 is decomposed by fixed-parameter VMD. The key parameters of VMD are the same as those mentioned above, i.e., M and α are set to 4 and 2000, respectively.

IMF2 has the maximal kurtosis value 3.13 in the decomposition results. The time-domain waveform of IMF2 and its corresponding envelope spectrum are separately illustrated in Figures 15 and 16.

Comparing the above experimental results, the following conclusions may be drawn: both parameter-adaptive VMD and EMD can demonstrate clear impact impulses in proper decomposition components, as shown in Figures 11 and 13. The corresponding envelope spectra shown in Figures 12 and 14 highlight the fault characteristic frequency of outer raceway fo and its harmonics, which clearly indicate the outer raceway of the bearing is defective. By contrast, the results decomposed by parameter-adaptive VMD exhibit more prominent and clearer fault characteristic information, while sensitive mode4 decomposed by EMD involves more noise interference.

It can be observed that IMF2 decomposed by fixed-parameter VMD shows little sign of impact impulses but strong noise interference. And it is hard to distinguish the fault characteristic components in the corresponding envelope spectrum displayed in Figure 16. The analysis results indicate that inappropriate parameters of VMD may not identify the fault characteristic frequency.

4.3. Incipient Rolling Element Fault

Figure 17 illustrates a time-domain vibration signal with the incipient rolling element fault. The optimal M and α sought by BAS are 6 and 170, respectively. Using parameter-adaptive VMD, the vibration signal is processed to acquire the IMF components, and the results are displayed in Figure 18. The kurtosis values of the resulting IMFs are calculated and listed in Table 6.

The envelope spectrum of IMF2 with maximal kurtosis value is displayed in Figure 19. As shown, the raw vibration signal is contaminated by heavy noise, and the weak fault characteristics are buried by other interference frequency components. Nevertheless, satisfactory results can be acquired by the proposed method. As illustrated in Figure 19, the fault characteristic frequency of rolling element fb (126.3 Hz) along with its second and third harmonics can be easily recognized.

The signal shown in Figure 17 is processed by EMD, in which decomposed mode4 has the maximal kurtosis value 10.73. Mode4 and its envelope spectrum are displayed in Figures 20 and 21, respectively.

Comparing Figure 21 with Figure 19, it can be noted that the fault characteristic frequency fb acquired by EMD may be roughly identified, but the envelope spectrum of mode4 is mixed with a large amount of noise. It may have something to do with mode mixing effects in the EMD method.

Figure 22 illustrates the time-domain waveform of IMF2 with the maximal kurtosis value 3.06, which is decomposed by fixed-parameter VMD. The envelope spectrum of IMF2 is illustrated in Figure 23.

It is difficult to identify the fault characteristic components in the envelope spectrum shown in Figure 23. The existence of serious mode mixing leads to the selected IMF containing heavy noise and other frequency components. The strong interference masks the useful fault characteristic information and results in the failure of bearing incipient fault diagnosis.

For further comparison, an extreme learning machine (ELM) [2931] classifier is adopted to estimate the proposed method. Considering that the energy at the fault characteristic frequencies will increase if the bearing is faulty, the proportion of fault energy (PFE) is firstly introduced. As aforementioned, the measured fault characteristic frequencies may be slightly different from the theoretical ones in the envelope spectrum. Therefore, the calculated fault characteristic frequency f is updated to a narrow band frequency range , where Δf is set to 5 Hz. The fault energy E1 can be calculated by the following equation:

The frequency range [0.5f, 1.5f] is chosen to calculate the signal energy E1a to avoid the influences derived from other frequency bands.

The ratio of is calculated to evaluate the proportion of the 1st harmonic fault energy. Similarly, the 2nd harmonic fault energy E2 is calculated. Considering that the deviation between the theoretical 2nd harmonic and the measured one will increase further, the 2nd harmonic 2f is updated to , where Δf is still set to 5 Hz.

The frequency range [1.5f, 2.5f] is chosen to calculate the signal energy E2a:

The ratio of is calculated to evaluate the proportion of the 2nd harmonic fault energy. The deviation between the theoretical 3rd harmonic and the measured one will continue to increase, which may introduce a significant error. Therefore, the proportion of the 3rd harmonic fault energy is no longer calculated. The PFEs of the envelope spectra above are calculated and listed in Table 7.

As displayed in Table 7, the PFEs obtained by the proposed method are higher than those obtained by the other two methods, which is beneficial for incipient fault diagnosis of bearing.

In this study, K, R1i + R2i, R1o + R2o, and R1b + R2b are calculated and chosen as the input vector of ELM. Four states of vibration signals, separately corresponding to the normal state, inner raceway fault, outer raceway fault, and rolling element fault, are assigned as the output vector of ELM. There are four groups of signals corresponding to the four states. Each of four groups involves 20 training signals and 20 testing signals. The hidden layer nodes of ELM are set to 80. The testing results of ELM are listed in Table 8.

As shown in Table 8, the fault diagnosis accuracy of the proposed method is higher than that of the other two methods. Compared with the proposed method and EMD, the testing accuracy of the fixed-parameter VMD is much lower, which indicates that M and α will greatly affect the VMD results.

5. Conclusions

The VMD method can completely decompose the raw vibration signals into a set of IMFs from low frequency to high frequency. The mode number M and the quadratic penalty term α of VMD will have considerable influence on decomposition results. The key to the VMD method lies in seeking for the optimal combination of M and α that matches with the analyzed signals.

In this paper, the BAS optimization algorithm is employed to adaptively estimate the optimal M and α. The optimized parameters can guarantee the availability of VMD.

The proposed parameter-adaptive VMD method has been applied in the field tests for low-noise deep groove ball bearing 6203. Comparisons have also been conducted to evaluate the performances by using the proposed method, EMD, and the fixed-parameter VMD. Three cases studies demonstrate that the proposed method outperforms EMD and the fixed-parameter VMD in suppressing the noise interference and highlighting the weak fault characteristic frequency information contaminated by heavy noise.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.