1. Introduction
Rolling bearings are the key component of rotating machinery. Their vibration signals are generally very complex, and especially when faults occur, the vibration signals are complicated by strong background noise, and obvious non-linear and non-stationary features. The time-frequency analysis method is to describe the changed frequency spectrum of the signal with time, establishing a distribution spectrum that can simultaneously express the energy or intensity of the signal in time and frequency [
1]. They can extract the contained fault features from the signal. The time-frequency analysis method has been widely applied in the field of mechanical fault diagnosis because it can provide local information about signals in the time domain and frequency domain [
2,
3,
4,
5,
6,
7]. The time-frequency analysis methods include Wigner-Ville distribution, wavelet transform, S transform, spectral kurtosis, Gabor transform, sparse decomposition, empirical mode decomposition (EMD), local mean decomposition (LMD), local feature scale decomposition (LCD), ensemble empirical mode decomposition (EEMD) and so on [
8,
9,
10,
11,
12]. Research on wavelet transforms in fault diagnosis is mainly based on combining wavelet analysis with other signal analysis methods, diagnosis methods. EMD is an adaptive signal time-frequency processing method proposed by Huang. It can directly decompose the signal and adaptively obtain the basis function, but it suffers from endpoint effects and mode mixing problems, which limit its application in actual engineering problems [
13,
14,
15,
16,
17,
18]. The EEMD method can effectively solve the mode mixing phenomenon, but this method requires a large amount of computation time.
Empirical wavelet transform (EWT) is a novel signal processing method proposed by Gilles [
19]. This method combines the good characteristics of the wavelet transform and EMD methods, and can extract a series of amplitude modulated-frequency modulated (AM-FM) signals from the given signal. The EWT is used for weak fault diagnosis and identification. A compact supporting Fourier spectrum can be obtained by using the AM-FM components. The adaptive wavelets capable of AM-FM components are constructed using differential methods [
20]. The differential mode is used to segment the Fourier spectrum, and some filters are used for each obtained supporting feature. The EWT has good adaptability and a reliable theoretical mathematical deduction basis, and has been effectively applied to engineering practice, but in actual engineering applications, high frequency interference components often exist in the signal spectrum due to strong background noise, which makes the signal spectrum interval of EWT decomposition not be accurate enough, and the mode components after decomposition will not contain accurate fault feature information. In recent years, some experts and scholars have deeply studied and improved the EWT method, and many improved EWT methods have been proposed. Amezquita-Sanchez and Adeli [
21] proposed a new adaptive multiple signal classification EWT method in order to accurately represent time-frequency results. Kedadouche et al. [
22] proposed a new signal analysis method based on EWT and operational modal analysis, which is used to decompose the signal into multiple components to extract the related frequency. Zheng et al. [
23] proposed an adaptive parameterless EWT to fulfill an adaptive separation of the Fourier spectrum in EWT. Hu et al. [
24] proposed an enhanced EWT method based on the advantages of the waveform in order to eliminate the disadvantage of basic EWT in spectrum segmentation. Thirumala et al. [
25] proposed a generalized empirical wavelet transform to estimate the time-varying PQ indices for accurate interpretation of disturbances. Merainani et al. [
26] used the EWT method to develop an appropriate wavelet filter bank for the early detection and condition monitoring of faults. Shi et al. [
27] proposed an enhanced EWT by developing a feasible and efficient spectrum segmentation method in order to improve the accuracy of the EWT results. Hu et al. [
28] proposed an enhanced empirical wavelet transform algorithm to extract bearing fault features. Hu et al. [
29] proposed multi-taper empirical wavelet transform to separate the fault features from the masking components for a wind turbine planetary gearbox under nonstationary conditions. Wang et al. [
30] proposed a sparsity guided EWT method, which is used to automatically build Fourier segments for fault diagnosis. Song et al. [
31] proposed an improved empirical wavelet transform with adaptive empirical mode segmentation and the merging of redundant empirical modes for fault diagnosis of roller bearings. Merainani et al. [
32] proposed a new self-adaptive time-frequency analysis, called Hilbert empirical wavelet transform, for obtaining the instantaneous amplitude matrices of the vibration signals. Bhattacharyya et al. [
33] enhanced the existing empirical wavelet transform by using Fourier-Besse series expansion in order to obtain an improved TF representation of non-stationary signals.
Although these improved EWT methods can improve the quality of the decomposed mode components, the spectrum segmentation mode based on the frequency domain extremum point of the EWT method does not change. In addition, when motor bearing faults occur, the fault feature is usually weak and easily influenced by interference signals. The signal has a strong non-linear and non-stationary multi-component amplitude modulation feature. Because the fault features are reflected in the complex edge band information, when the EWT method is directly used to decompose the vibration signal, it will occur that many mode components will be close to the component and it is easy to lose the fault feature signal. Therefore, an enhanced empirical wavelet transform (MSCEWT) based on maximum-minimum length curve method is proposed in this paper. The maximum-minimum length curve method transforms the spectrum of the original vibration signal to the scale space in order to obtain a set of minimum length curves, and find the maximum length curve value in the set of the minimum length curve values for obtaining the number of the spectrum decomposition intervals. The MSCEWT method is used to decompose the vibration signal to obtain a series of IMFs. The Hilbert transform is the used to process these IMF components in order to obtain the fault feature frequencies of motor bearings. The actual motor bearing vibration signals are used to verify the effectiveness of the MSCEWT method for fault diagnosis.
The rest of the paper is organized as follows: in
Section 2, the empirical wavelet transform method is introduced. The empirical wavelet transform method based on scale space is introduced and analyzed in
Section 3. In the
Section 4, an enhanced EWT (MSCEWT) based on maximum-minimum length curve method is proposed to realize fault diagnosis of motor bearings. Experimental comparison and analyses are discussed in the
Section 5. Finally, our conclusions are drawn in
Section 6.
3. The Maximum-Minimum Length Curve Method
The EWT is an adaptive signal processing method. Based on the adaptive Fourier spectrum segmentation of the signal, a bank of wavelet filters is constructed in order to decompose the signal into a set of amplitude modulation and frequency modulation terms. Because the conventional LocalMaxmin method needs to determine the number of segmentations in advance, it is affected by human subjective factors in the process of implementation, which will result in the inaccurate interval decomposition and cannot achieve the purpose of adaptive spectrum segmentation. Gilles proposed the adaptive and adaptivereg methods to segment the spectrum, but one needs to preset the initial boundary vector when the EWT method segments the interval. Next, Gilles also proposed a new method of Fourier spectrum segmentation based on scale-space representation and K-Means clustering algorithm, which uses the scale space to achieve adaptive segmentation of the spectrum without presetting interval numbers.
Let a function
be defined in the interval
and the Gauss kernel function
is defined. The scale space representation is described as follows:
where
denotes the convolution product and denotes scale parameter.
In practical applications, the Fourier transform
of the original function and the sampled Gaussian kernel function are used to describe the scale space representation:
When M is large enough, the error of Gaussian approximation is negligible. In the text, is set, where . In order to meet the error that is small enough, the set .
The scale space representation can be interpreted as the global trend of signals under different scale parameters, that is, all the scale space curves under different scale parameters are calculated, and then the threshold of the scale parameter is determined. For the scale space representation under the initial scale parameter ( = 0.5), the minimum value among the local maximum values in is obtained. When there is t = 0, the number of local minimum value obtained is obtained and regarded as . These local minimum values are used to construct the scale space curve . In the scale space representation plane , the calculated at t = 0 is recorded as the number of minimum length curves . With the increasing of the scale t, the number of the local minimum values in is decreasing until the number of minimums is reduced to 0. At this time, the value of the minimum length curve is calculated. The local minimum value of at t = 0 is used to represent the starting of scale space curve that the first length of scale space curve is 1. The scale space curve is sequentially accumulated under different scales () in order to obtain the first minimum length curve . Then the position of the local minimum value is determined at t = 0 in turn. From the starting points of these positions, the scale space curves are sequentially accumulated under different scales () in order to obtain the minimum length curves . Therefore, a proportional threshold T can be set from the minimum length curve , that is to say all minimum length curves corresponding to the proportional threshold T is set in the minimum length curves . That is, all obtained minimum length curves are greater than the threshold T, which are found out. The problem is transformed into a clustering problem, so that the clusters of the set are divided into two clusters (meaningful/non-meaningful minima). Here, the method of k-Means is used for clustering.
The calculation flow of the minimum length curve is shown in
Figure 2.
The calculation steps of the minimum length curve are described as follows:
Step 1. Initialize the scale space plane, including the first column number i in the scale space plane (the local minimum point of original function), and the number ic of the local minimum value of the original function.
Step 2. Determine whether the first column position in the scale space plane is 1. If the first column position is 1, Step 3 is executed. Otherwise execute Step 6.
Step 3. When the first position in the scale space plane exists the local minimum value of the original function, the scale space curves are sequentially accumulated under different scales in order to obtain the minimum length curves . Otherwise execute Step 4.
Step 4. When the last position in the scale space plane exists the local minimum value of the original function, the scale space curves are sequentially accumulated under different scales in order to obtain the minimum length curves . Otherwise execute Step 5.
Step 5. When the local minimum value of the original function in the scale space plane does not exist the first position or the last position, the scale space curves are sequentially accumulated under different scales in order to obtain the minimum length curves . Otherwise execute Step 7.
Step 6. Calculate the first column number i = i + 1 in the scale space plane.
Step 7. If the end condition i ≤ size (plane, 1) is satisfied, the minimum length curve L(ic) is output. Otherwise go to Step 6.
4. An Enhanced EWT (MSCEWT) Based on Maximum-Minimum Length Curve
4.1. The Idea of the MSCEWT Method
As a novel signal processing method, the EWT method still has a problem to extract physical frequency components from actual complex mechanical signals. In the process of Fourier spectrum segmentation of the signal, when a fine frequency interval is set, many narrowband modes are extracted, and several modes could display the same modulation information, which will result in unnecessary redundancy. In order to avoid this problem of the EWT method, an improved scale space representation (maximum-minimum length curve) method is proposed in this paper. That is an enhanced empirical wavelet transform (MSCEWT) method based on maximum-minimum length curve, which is used to decompose mechanical vibration signals in order to realize mechanical fault diagnosis. In the MSCEWT, the local minimum value in is used to adaptively segment the spectrum in order to determine the boundaries, and the corresponding position of the maximum value of all the minimum length curves in the scale space representation is determined. Then the signal is decomposed into as a series of single packets with fault information, and the power spectral density is applied to the single component in order to extract the modulation information. The MSCEWT method can avoid the over decomposition phenomenon of the mode at utmost.
4.2. The Flow and Steps of the MSCEWT Method
Firstly, the Fourier transform
of the original function and the sampled Gaussian kernel function
are convoluted in order to obtain the scale space representation. Secondly, in the calculation process, with the increasing of the scale
t, the scale space curve
is obtained. Then, the position of the local minimum is determined in turn at
t = 0. From the starting points of these positions, the scale space curves are sequentially accumulated under different scales
(
) in order to obtain the minimum length curves
. The final difference between scale space representation method and maximum-minimum length curve is that the problem is converted to calculate all the maximum thresholds of the minimum length curve
in the scale space representation set
, then the boundaries are determined by detecting the local minimum value. The number
N of the obtained maximum threshold by calculating is defined as the number of spectrum segmentation. The flow of the MSCEWT method is shown in
Figure 3. The specific steps of the MSCEWT method are described as follows:
Step 1. The original signal is transformed to obtain Fourier spectrum by FFT.
Step 2. Execute convolution operations between and the Gaussian kernel function to obtain a scale space representation, denoted as .
Step 3. Calculatethe point of the local minimum value of when there is t = 0. From the starting points of these positions, the scale space curves are sequentially accumulated under different scales () in order to obtain the minimum length curves .
Step 4. Calculate all maximum of the minimum length curve in the set .
Step 5. Calculate the set of the maximum values of the minimum length curve, and the number of the maximum values is denoted as N. These points are used as the segmentation points of the spectrum segmentation interval.
Step 6. Establish a band-pass filter based on wavelet transform on the spectrum segmentation interval.
Step 7. Reconstruct the signal in order to obtain the IMFs component.
Step 8. The power spectrum analysis is performed on the IMFs component in order to extract the fault feature frequency.
4.3. Experimental Environment and Results
4.3.1. Experimental Environment and Data
The experiment vibration data comes from Bearing Data Center of Case Western Reserve University [
34]. The experimental platform is shown in
Figure 4. The test stand consists of a 2 hp motor, a torque transducer/encoder, a dynamometer, and control electronics. The experiment uses 6205-2RS JEM SKF deep groove ball bearings. The motor is coupled to the dynamometer and torque sensor through self-aligned coupling. The test bearings support the motor shaft. Single point faults were introduced to the test bearings using electro-discharge machining with fault diameters of 7, 14, 21, 28 and 40 mils (1 mil = 0.001 inches). The data source are the accelerometers on the motor housing of the motor drive end. The vibration signals have a rotational speed of 1797 r/min. Vibration data was collected using accelerometers, which were attached to the housing with magnetic bases. Accelerometers were placed at the 12 o’clock position at both the drive end and fan end of the motor housing. In order to quantify this effect, experiments were conducted for both fan and drive end bearings with outer raceway faults located at 3 o’clock (directly in the load zone), at 6 o’clock (orthogonal to the load zone), and at 12 o’clock. Vibration signals were collected using a 16 channel DAT recorder, and were post processed in a Matlab environment. The vibration signal sampling frequency of the motor bearing is 12,800 Hz, and the duration of each vibration signal is 10 s. Divide the original vibration signals into multiple samples, each with a data sample length of 2048 points. According to the theoretical calculation formula of the fault characteristic frequency, the fault frequency of the outer race, inner race and rolling element of the motor bearing is obtained [
35]. The experiment environment is described: Intel Core I5 2450, 4 GB RAM, Win 7 and Matlab 2014b. The fault diameter (14 mils) is selected in here, and the calculated results are shown in
Table 1.
4.3.2. Experimental Results
The inner race fault vibration signal of motor bearing is selected to verify the effectiveness of the improved scale space representation method. The spectrum segmentation and power spectrum of the inner race vibration signal of the motor bearing are shown in
Figure 5 and
Figure 6, respectively. The power spectrum results of scale space representation (SSR) and modified scale space representation based on maximum-minimum length curve (MSSR) are compared in
Table 2.
As can be seen from
Figure 5 and
Figure 6 and
Table 2, the scale space method is used to segment the inner race vibration signal of motor bearing, and the number of spectrum segmentation intervals is N = 8. According to the power spectrum of the IMFs component, the maximum frequency value of the extracted IMF3, IMF5, IMF6 and IMF7 is 29.30 Hz, the fault feature frequency value of the extracted IMF1 and IMF8 is 164.06 Hz. The modified scale space representation based on maximum-minimum length curve (MSSR) is used to spectrally segment the inner race vibration signal of motor bearing, and the number of spectrum segmentation intervals is N = 5.
According to the power spectrum of the IMFs component, the maximum frequency value of the extracted IMF3 and IMF4 is 29.30 Hz, the maximum frequency value of the extracted IMF1 is 357.42 Hz, the maximum frequency value of the extracted IMF2 is 117.19 Hz, and the maximum frequency value of the extracted IMF5 is 164.06 Hz. 29.30 Hz is the rotation frequency, 117.19 Hz is the quadruple frequency of rotation frequency, 357.42 Hz is the 12 times frequency of rotation frequency, and 164.06 Hz is the fault feature frequency. Therefore, the modified scale space representation based on maximum-minimum length curve reduces the obtained IMFs components of vibration signal from 8 to 5, which reduces the number of adaptive spectrum segmentation and the redundancy of the same frequency, effectively reduces the same modulation of the obtained IMFs and accurately extracts the rotation frequency and the fault feature frequency.