MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2014/765621 765621 Research Article Rolling Bearing Fault Diagnosis under Variable Conditions Using Hilbert-Huang Transform and Singular Value Decomposition Liu Hongmei 1 Wang Xuan 1 Lu Chen 1,2 Wong Pak-Kin 1 School of Reliability and Systems Engineering Beihang University Beijing 100191 China buaa.edu.cn 2 Science & Technology on Reliability & Environmental Engineering Laboratory Beijing 100191 China 2014 1672014 2014 11 04 2014 08 06 2014 25 06 2014 16 7 2014 2014 Copyright © 2014 Hongmei Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Fault diagnosis precision for rolling bearings under variable conditions has always been unsatisfactory. To solve this problem, a fault diagnosis method combining Hilbert-Huang transform (HHT), singular value decomposition (SVD), and Elman neural network is proposed in this paper. The method includes three steps. First, instantaneous amplitude matrices were obtained by using HHT from rolling bearing signals. Second, the singular value vector was acquired by applying SVD to the instantaneous amplitude matrices, thus reducing the dimension of the instantaneous amplitude matrix and obtaining the fault feature insensitive to working condition variation. Finally, an Elman neural network was applied to the rolling bearing fault diagnosis under variable working conditions according to the extracted feature vector. The experimental results show that the proposed method can effectively classify rolling bearing fault modes with high precision under different operating conditions. Moreover, the performance of the proposed HHT-SVD-Elman method has an advantage over that of EMD-SVD or WPT-PCA for feature extraction and Support Vector Machine (SVM) or Extreme Learning Machine (ELM) for classification.

1. Introduction

In modern rapidly developing industries, the rolling bearing, a vital component of most rotary machines, faces an increasingly complex working environment. For the purpose of satisfying the production requirement, rolling bearing fault diagnosis is essential to ensure that rolling bearings work in a good state. In essence, rolling bearing fault diagnosis is approximately divided into several steps: signal processing to extract the feature vector, pattern recognition, and classification according to the extracted fault feature vector. Therefore, it is extremely important that the fault feature vector correctly and efficiently represents the fault characteristics hidden in the original signal, which also directly affects the classification results and diagnostic accuracy. In rolling bearing fault diagnosis, the variation of working conditions frequently leads to the change of most feature vital parameters  and even the entire diagnostic method. Adjusting the feature parameter or the diagnostic method according to the working condition is particularly labor intensive and time consuming. Thus, it is essential to find a rolling bearing diagnostic method that is applicable to different working conditions.

Recently, researchers have paid more attention to the study of the rolling bearing fault diagnosis under variable working conditions. Based on a smart bearing, Shao presented a fault diagnosis system including several sensing devices to track the time-varying parameter . Yang et al. introduced variable predictive model based class discrimination (VPMCD) into rolling bearing fault diagnosis . In order to achieve the fault diagnosis under variable conditions, Yang et al. employed the feature vector obtained from different working conditions as the input of the classifier. He et al. addressed manifold learning on generated time-frequency distributions for machine fault signature analysis . These previous studies have a pivotal function in fault diagnosis under variable conditions. However, some of these studies require additional sensing devices that increase the system cost; some conduct feature extraction under as many working conditions as possible to ensure the robustness; some decrease the algorithm robustness because of complex calculation. How to certify the robustness and low costing of the fault diagnosis under variable working conditions still needs further study.

At present, there exist many conventional methods employed to process the vibration signal: the Fourier transform, the short-time Fourier transform, the Wigner-Ville distribution, and the wavelet transform. When processing the signal, the Fourier transform is only applicable to linear system and stationary signal. Otherwise, the result will lose its physical meaning. However, the vibration signal of rollingbearings is characterized by nonstationarity, which makes the Fourier transform unsuitable for extracting the feature of the nonstationary signal. As time-frequency analysis methods, the short-time Fourier transform , the Wigner-Ville distribution , and the wavelet transform  can process the nonstationary signal. However, they all have respective limitation. Due to Heisenberg’s uncertainty principle, the finest time location and the best frequency resolution cannot be reached simultaneously , which restricts the application of the short-time Fourier transform. Because of the presence of cross-term interference, the Wigner-Ville distribution is limited. Meanwhile, although the wavelet transform has demonstrated its superiority in processing the nonstationary signal, the energy leakage issue is inevitable. Besides, the base function needs a preselection according to the structure of the signal, and the optimization problem for base function preselection still requires a widely accepted resolution. By contrast, the Hilbert-Huang transform (HHT)  is free from the restriction of Heisenberg’s uncertainty principle. Its essence is processing the signal with Empirical Mode Decomposition (EMD)  to get the intrinsic mode function (IMF) and applying the Hilbert transform to each IMF for Hilbert spectrum analysis. Based on its own time scale characteristics, the EMD method decomposes signals without any basis function preselection. In theory, the EMD method can be applied to any type of signal decomposition; therefore, it has an extremely distinct advantage in handling nonstationary and nonlinear data. It is able to decompose a complicated signal into a finite number of IMFs, and the decomposed IMFs contain the local characteristic of the original signal at different time scales.

Under variable working conditions, the rolling bearing fault characteristic frequency changes over time, thereby necessitating time-frequency analysis. Compared with other time-frequency methods, Hilbert-Huang transform (HHT) has a distinct advantage in handling nonstationary and nonlinear data. Consequently, HHT is implemented in this research for time-frequency analysis. Nevertheless, the result obtained from HHT under variable working conditions is tremendous and complex; thus, dimensionality reduction is required for the purpose of enhancing classification accuracy. In addition, it is extremely complicated to extract the corresponding amplitude at the time-varying characteristic frequency. In comparison, searching a fault feature insensitive to the working condition variation is preferred. Theoretically, the PCA and the singular value decomposition (SVD) can both be employed for dimensionality reduction. Moreover, the singular value has considerable stability  and changes negligibly when matrix elements change. Thus, the singular value can maintain satisfactory stability under variable working conditions, thereby enhancing the robustness of the fault diagnosis under variable conditions.

Thus, HHT in tandem with SVD is proposed for rolling bearing fault feature extraction under variable working conditions in this paper. To begin with, time-frequency analysis is implemented on the vibration signal by utilizing HHT, and the instantaneous amplitude and instantaneous frequency are obtained. Afterwards, for the purpose of dimensionality reduction and obtaining the feature insensitive to working condition variation, SVD is employed to the instantaneous amplitude to obtain the feature vector. After the feature vector is obtained, some intelligent classification techniques can be implemented for intelligent rolling bearing fault classification, such as Support Vector Machine (SVM) , Relevance Vector Machine (RVM) , and Extreme Learning Machine . Similar to other binary classification algorithms, SVM cannot solve multiclassification efficiently. RVM is restricted by the long training time. As a feedforward network, ELM is characterized by fast learning and considerable generalization performance. However, the model accuracy of the feedforward network is lower than that of the feedback network. Consequently, Elman neural network, a typical feedback neural network, is employed for rolling bearing fault identification and classification in this paper.

This paper is organized as follows. Section 2 introduces the principle of the diagnostic method based on HHT and SVD. Section 3 presents fault identification and classification scheme. The experimental results are reported in Section 4. Conclusions are drawn in Section 5.

2. Feature Extraction Method Based on HHT and SVD

Fault diagnosis for rolling bearings is comprised of fault feature extraction and fault classification. Fault features extraction is the most important procedure during the fault diagnosis, and the fault is then classified according to the extracted feature.

2.1. Time-Frequency Signal Decomposition Based on HHT

As a self-adaptive data decomposition method, HHT is based on partial time scales and therefore can deal with both nonlinear and nonstationary process. HHT could transfer any complex signal into the linear superposition of several intrinsic mode function (IMF) components. In HHT, the analytic signal can be constructed by Hilbert transform, in advance; the instantaneous amplitude and instantaneous frequency can be obtained. Additionally, the Hilbert spectrum can be obtained after further processing.

HHT , which combines Empirical Mode Decomposition with Hilbert transform, can be applied to determine the instantaneous amplitude and the instantaneous frequency.

( 1) Empirical Mode Decomposition. After the implementation of the EMD, the original signal x ( t ) can be decomposed as follows: (1) x ( t ) = i = 1 n c i ( t ) + r n ( t ) , where c i ( t ) is the i th IMF and r n ( t ) is the residue. As the constant sequence or monotonic sequence, r n ( t ) represents the average trend of signal. The IMF component c i ( i = 1 , , n ) represents the signal from high frequency to low frequency in different frequency bands. The frequency component of the signal in each frequency band changes according to the original signal.

From the perspective of the basis function, the basis functions are different during the application of EMD. It differs from the Fourier transform (taking a series of sine waves or cosine waves with constant amplitude and frequency as the basis function) and the wavelet transform (taking the preset basis function). Consequently, EMD is beneficial to enhance the signal decomposition efficiency and to process the nonstationary data.

( 2) Hilbert Transform. In essence, Hilbert transform is the convolution of signal with 1 / π t , which makes it possible to calculate instantaneous parameters for transient signal and complex signal.

When the IMF is set as equal to c ( t ) , the analytic signal h ( t ) of c ( t ) is as follows: (2) h ( t ) = c ( t ) + j c ~ ( t ) = a ( t ) e j φ ( t ) , where c ~ ( t ) is the outcome function after applying Hilbert transform to c ( t ) and is expressed as (3) c ~ ( t ) = 1 π - + c ( t ) t - τ d τ .

a ( t ) is the amplitude function and the envelope signal of the IMF component c ( t ) is expressed as (4) a ( t ) = c 2 ( t ) + c ~ 2 ( t ) .

φ ( t ) is the phase function: (5) φ ( t ) = ta n - 1 ( c ~ ( t ) c ( t ) ) .

We calculate the derivation of the phase function to obtain the instantaneous frequency, as follows: (6) ω ( t ) = d φ ( t ) d t .

Accordingly, the original signal x ( t ) can be expressed by a time-dependent function containing instantaneous amplitude and instantaneous frequency, as follows: (7) H ( w , t ) = Re i = 1 n a i ( t ) exp ( j w i ( t ) d t ) .

In this way, we are equipped with the ability to conduct the observation and analysis intuitively from the Hilbert spectrum, which is a three-dimensional plot with time axis, frequency axis, and amplitude axis. Because the Hilbert spectrum can provide the amplitude contribution at each frequency and each time , it is possible to handle the rolling bearing fault diagnosis under time-varying conditions by applying HHT.

The instantaneous amplitude of the IMF is its envelope signal, which is similar to the periodic signal with large amplitude, long duration, and unchangeable repetition rate, as well as no low-frequency interference . The envelope signal is beneficial in identifying faults, even weak ones. Thus, we can calculate the amplitude of each analytic signal h ( t ) to obtain a ( t ) for further analysis. However, it is complicated to calculate the accurate fault characteristic frequency and to extract the corresponding amplitude from a ( t ) under variable working conditions. Thus, rolling bearing fault identification is difficult under variable operating conditions. To address this issue, we further propose SVD in tandem with HHT for rolling bearing fault feature extraction.

2.2. SVD Based on HHT

SVD is a matrix transformation algorithm, which is defined in the following process. We assume that M is a m × n matrix. In SVD, any matrix M R m * n can be decomposed into three matrices, which are U R m * m , R m * n , and V R n * n  as follows: (8) M = U V * , where U and V are orthogonal and is a diagonal matrix with r = rank ( M ) principal positive diagonal entries: (9) = ( σ 1 0 0 0 0 σ 2 0 0 0 0 σ r - 1 0 0 0 0 σ r ) .

σ i are considered as the singular values of M . The singular values are listed in descending order automatically by SVD function, as follows: (10) σ 1 ( M ) σ 2 ( M ) σ r ( M ) .

The singular value has considerable stability and changes negligibly when matrix elements change . In this paper, SVD is used to determine the intrinsic characteristic of instantaneous amplitude matrices obtained from HHT. Singular values can express the feature of the original matrix in the form of several values, which is conducive to compressing the scale of the feature vector.

As shown in Figure 1, the original vibration signal is firstly decomposed by EMD, and IMFs are then obtained. The analytic signal h i is obtained after each major IMF is processed by the Hilbert transform (see formula (2)). Subsequently, the instantaneous amplitude a i is obtained after calculating the amplitude function of each h i (see formula (4)). The feature matrix M is constructed by a i ( i = 1 , , k ). To a certain degree, the number of IMFs is decided by the characteristic, length, and complexity of the target signal. For the purpose of convenient analysis, for a signal with given length, normally the signal needs to be divided into a number of subsignals with fixed length. Generally, due to the difference in the subsignals, the numbers of obtained IMFs will normally be different. To all subsignals, if the minimum number of obtained IMFs is k 1 , the feature matrix’s dimension k should be no larger than k 1 . In addition, the first few IMFs have the highest frequency and largest energy, while the last few IMFs keep relatively moderate. That is, the first few IMFs comprise most fault information. Consequently, the prior k IMFs are then extracted for analysis. From the perspective of practice, k is usually set as 3~5. In this paper, according to the decomposition results of raw vibration signal, we set k = 5 , so the feature matrix M can be expressed as follows: (11) M = ( a 1 a 2 a k ) .

Flowchart of obtaining the feature vector.

After the application of SVD to M , the singular value vector as well as the fault feature vector σ ( M ) can be obtained as follows: (12) σ ( M ) = [ σ M 1 , σ M 2 , , σ M k ] .

The fault is then identified and classified according to the obtained fault feature vector σ ( M ) .

3. Fault Identification and Classification Based on Elman Neural Network

In intelligent fault diagnosis, neural networks are always used to classify the fault mode. In this study, the Elman neural network was used for fault classification.

The Elman neural network includes input, output, hidden, and context layers. When the output of the hidden layer neurons goes through the context layer, it associates with the input of the hidden layer . This type of association is sensitive to the historical state of data; internal feedback increases the capacity of the network to deal with dynamic information, overcoming the drawback of the feed-forward network, which does not have dynamic characteristics. Therefore, compared with other neural networks, the Elman neural network is more applicable to classify the fault mode. The structure of the Elman neural network is shown in Figure 2.

Elman neural network.

The transfer function of the hidden layer in the Elman neural network is a sigmoid function. The activation function of the output layer neuron is a linear function. The nonlinear state equations of the network are (13) y ( k ) = g [ w 3 x ( k ) ] , where y represents the output vector of the network, x denotes the output vector of the hidden layer, w 3 is the connection matrix from the context layer to the output layer, and g ( · ) denotes the activity function of the output neurons. Consider the following: (14) x ( k ) = f [ w 1 x c ( k ) + w 2 u ( k - 1 ) ] , x c ( k ) = x ( k - 1 ) . u is the input vector of the network, w 2 is the connection matrix from the input layer to the hidden layer, w 1 represents the linked weights from the context layer to the hidden layer, x c represents the feedback state vector, and f ( · ) represents the transfer function of the hidden layer.

4. Experimental Result and Discussion

To validate the effectiveness of the proposed fault diagnosis method, we used the vibration data from Bearing Data Center of Case Western Reserve University. The 6205-2RS deep-groove ball bearing is employed in the experiment. As shown in Figure 3, the test bed consists of 1.5 kW motor, torque sensor/encoder, dynamometer, and electrical control devices. The motor drives the input shaft and the output shaft drives the load. The vibration data are sampled under different working conditions: for motor loads of 0 to 3 horsepower (HP) (motor speeds of 1797 to 1720 RPM). In addition, different faults such as inner-race fault, outer-race fault, and rolling element fault were injected in the experiment. The data of different fault states are acquired with the sampling frequency of 12 kHz.

Test bed.

4.1. Performance for Different Working Conditions

Before the implementation of SVD, the original vibration signal is processed by HHT to obtain the Hilbert spectrum and the instantaneous amplitude. The Hilbert spectrum of different states is shown in Figure 4(a)~Figure 4(d). The Hilbert spectrum of the normal state in Figure 4(a) indicates that the energy distribution is uniform, whereas the peak appears in the frequency band at approximately 1,000 Hz. The maximum amplitude is approximately 0.09. The Hilbert spectrum in Figure 4(b) indicates that the dominant energy distributes in the frequency band at approximately 3,000 Hz, and the time-frequency distribution is complicated under the inner-race fault mode. The maximum amplitude is approximately 0.9. For the outer race fault mode, the Hilbert spectrum in Figure 4(c) shows that the dominant energy distribution is similar to that of the inner-race fault mode. Nevertheless, the maximum amplitude is approximately 3, which indicates that the vibration energy is relatively high. The Hilbert spectrum of the rolling element fault illustrates that energy rarely distributes in the frequency band from 2,000 Hz to 3,000 Hz. Moreover, in Figure 4(d), the vibration energy is low and its maximum amplitude is only approximately 0.35. The time-frequency distribution varies under different states. Therefore, the feature matrix can be constructed for fault identification.

Hilbert spectra of (a) normal, (b) inner-race fault, (c) outer-race fault, and (d) rolling element fault.

After obtaining the instantaneous amplitude matrices, we can obtain the singular value vector by conducting SVD. As shown in Table 1, four datasets sampled under normal state, inner-race fault state, outer-race fault state, and rolling element fault state were obtained. Each dataset is divided into four subdatasets by different working conditions, and each subdataset contains 24 groups of data. Therefore, the whole four datasets contain 384 samples, and the feature extraction method based on HHT and SVD is implemented on each sample. For visualization, the scatter plots of the first three singular values extracted by HHT and SVD for datasets in Table 1 are shown in Figure 5. As shown in Figure 5, even under different working conditions, the singular value vectors extracted under the same failure mode retain a high degree of coincidence. Namely, the feature vector extracted by HHT in tandem with SVD is insensitive to the working condition variation. In addition, the singular value vectors corresponding to the three fault modes and the normal state have satisfactory mode separability. Due to the good separability of different fault states, the singular value can be feasibly used as the input of the neural network classifier.

Relevant information for dataset.

Label Groups
Status Working condition
0 HP (1797 RPM) 1 HP (1772 RPM) 2 HP (1750 RPM) 3 HP (1720 RPM)
1 Normal 24 24 24 24
2 Inner-race fault 24 24 24 24
3 Outer-race fault 24 24 24 24
4 Rolling element fault 24 24 24 24

Scatter plots of the first three singular values extracted by HHT-SVD for datasets in Table 1.

4.2. Comparison between the Proposed Method with EMD and WPT in Tandem with PCA for Feature Extraction

In this subsection, another two widely applied methods, EMD and WPT (wavelet packet transform) in tandem with PCA, are also implemented on the same datasets in Table 1, and their performance is compared with that of HHT in tandem with SVD for rolling bearing feature extraction under variable working conditions.

To begin with, EMD in tandem with SVD [19, 20] is exploited, SVD is applied to the first 5 IMFs, and then the feature vector is extracted. The scatter plots of the first three singular values extracted by EMD and SVD for datasets in Table 1 are shown in Figure 6. As shown in Figure 6, there exist major fluctuations for the corresponding feature vector, especially for inner-race fault and outer-race fault; that is, the feature vector obtained by EMD and SVD is affected greatly by working condition variation. The clustering interval between inner-race fault and outer-race fault is relatively small, which may lead to the misclassification. Consequently, the feature extraction method based on EMD and SVD is unable to work effectively and stably under variable working conditions.

Scatter plots of the first three singular values extracted by EMD-SVD for datasets in Table 1.

Afterwards, another widely used method, WPT in tandem with PCA , is employed for feature extraction. The signal is decomposed by WPT into 3 levels, and then the signal can be decomposed into 8 frequency bands. The energy of each frequency band is then calculated to obtain wavelet packet energy feature. Subsequently, PCA is used to extract principal components from the wavelet packet energy feature for dimensionality reduction. The scatter plots of the first three PCs (principal components) extracted by WPT and PCA for datasets in Table 1 are illustrated in Figure 7. As shown in Figure 7, there exist major fluctuations for all states. Namely, the feature vector obtained by WPT and PCA is sensitive to the working condition variation. In addition, the clustering intervals between three fault states are relatively small, which may lead to the incorrect classification. Consequently, WPT in tandem with PCA is also unsuited for rolling bearing fault diagnosis under variable working conditions. By contrast, the feature vector extracted by HHT in tandem with SVD keeps considerable stability even under different working conditions, and the mode separability for different fault states is satisfactory. For rolling bearing feature extraction under variable working conditions, the proposed method performs better than EMD in tandem with SVD and WPT in tandem with PCA.

Scatter plots of the first three PCs extracted by WPT and PCA for datasets in Table 1.

4.3. Classification Performance Using Elman Neural Network

Through the proposed feature extraction based on HHT and SVD, singular value vectors are obtained as the feature vectors. The faults status is then identified and classified by using the Elman neural network. The input of the Elman neural network is the singular value feature, and the output is set as (1 0 0 0), (0 1 0 0), (0 0 1 0), and (0 0 0 1), which represent normal, inner-race fault, outer-race fault, and rolling element fault in sequence. As shown in Table 2, the datasets are divided into a training dataset with 20 groups of data (the neural network training set is not listed in this paper because of space limitations) and a testing dataset with 364 groups of data. The partial results of network neural recognition are reported in Table 3. In addition, another two classification techniques, Support Vector Machine (SVM) and Extreme Learning Machine (ELM), are applied to the same training data and testing data in Table 2. Afterwards, the corresponding classification performance of Elman neural networks, SVM, and ELM is illustrated in Figure 8.

Training data and testing data statistics.

Label Status Working condition
Testing data Training data
0 HP 1 HP 2 HP 3 HP 3 HP
1 Normal 24 24 24 19 5
2 Inner-race fault 24 24 24 19 5
3 Outer-race fault 24 24 24 19 5
4 Rolling element fault 24 24 24 19 5

Elman neural network recognition results.

Sequence State Operating condition Target output Actual output of network
1 Normal (1750 r/min, 2 HP) (1 0 0 0) 0.993781 3.45 E - 08 0.043759 0.001039

2 Inner-race fault (1720 r/min, 3 HP) (0 1 0 0) 0.0072 0.9910 0.0072 0.0070
3 (1750 r/min, 2 HP) (0 1 0 0) 0.0067 0.9885 0.0078 0.0079
4 (1772 r/min, 1 HP) (0 1 0 0) 2.94 E - 06 0.999569 5.87 E - 05 0.078958
5 (1797 r/min, 0 HP) (0 1 0 0) 0.0067 0.9883 0.0078 0.0080

6 Outer-race fault (1720 r/min, 3 HP) (0 0 1 0) 0.000413 0.100543 0.990864 6 E - 11
7 (1750 r/min, 2 HP) (0 0 1 0) 7.241 e - 05 0.0155 0.9845 0.0128
8 (1772 r/min, 1 HP) (0 0 1 0) 0.000264 0.007905 0.997584 0.000573
9 (1797 r/min, 0 HP) (0 0 1 0) 5.880 e - 05 0.0109 0.9900 0.0112

10 Rolling element fault (1720 r/min, 3 HP) (0 0 0 1) 0.0051 0.000405 0.0109 0.9866
11 (1750 r/min, 2 HP) (0 0 0 1) 0.0051 0.000236 0.0117 0.9898
12 (1772 r/min, 1 HP) (0 0 0 1) 0.0064 0.000245 0.0100 0.9899
13 (1797 r/min, 0 HP) (0 0 0 1) 0.0100 0.000583 0.0072 0.9815

Classification performance of different classification techniques.

To verify the robustness of the proposed method, the Elman neural network is trained under a harsh situation: small training dataset (5 groups of data for each fault status) is used, and only data under a single working condition (3HP for this paper) is contained. Even so, the actual output of the neural network significantly agrees with the target output under variable conditions, as shown in Table 3. To a certain degree, parameters do not need to change under variable working conditions. As shown in Figure 8, the classification accuracy of Elman neural network is 100%. In comparison, some samples of inner-race fault are classified incorrectly into the outer-race fault by SVM. The same misclassification follows in ELM. The classification accuracy of SVM and ELM is 97.5% and 98.6%, respectively. The performance of Elman neural network outperforms that of SVM or ELM for classification. Consequently, the proposed method combining HHT-SVD with Elman neural network can be effectively used in rolling bearing fault diagnosis under variable working conditions.

4.4. Case Study for Rolling Bearing Fault Diagnosis under Time-Varying Conditions

To further verify the effectiveness of the rolling bearing fault diagnosis method proposed in this paper, this method is also used to detect the burst fault under the time-varying condition, which is exhibited in Figure 9. The time-varying data is acquired by connecting data segments under different working conditions and different fault modes. The initial state of the rolling bearing was normal, and an inner-race fault was injected at 0.5 s. The initial motor load was 3HP, which changed 0HP at 0.3 s and 1HP at 0.7 s. The instantaneous amplitude matrices are analyzed by conducting SVD for every 0.05 s. Figure 10 shows the corresponding instantaneous amplitude and singular value σ n ( t ) . The function σ n ( t ) remains relatively stable in one state. With the alteration of the motor load, the change of σ n ( t ) is not obvious. After an unexpected fault is injected, the value of σ 1 ( t ) changes rapidly. Moreover, even when the motor load changes to 1HP, the fault identification could still be achieved by observing the value of σ n ( t ) .

Diagram of the time-varying condition.

Diagram of instantaneous amplitude and singular value σ n ( t ) corresponding to Figure 9.

After the detection of a burst fault, this method is used to measure the fault severity under time-varying condition in advance, which is illustrated in Figure 11. To realize fault injection of different severities, grooves with different diameters were engraved on different positions of bearings (inner race, outer race, and rolling element), and fault diameter is used to represent fault severity. Here, the data is acquired by connecting data segments under different working conditions and different fault severities. The unit of fault diameter in this paper is inch. The initial motor load was 3HP, which changed to 0HP at 0.3 s and 1HP at 0.7 s. In addition, the initial fault diameter of inner-race fault was 0.007 inches, which changed to 0.014 inches at 0.2 s, 0.021 inches at 0.4 s, and 0.028 inches at 0.6 s. Similarly, the instantaneous amplitude matrices are analyzed by conducting SVD for every 0.05 s. Figure 12 shows the corresponding instantaneous amplitude and singular value σ n ( t ) , from which we can figure that the function σ n ( t ) remains relatively stable under a fault severity. With the alteration of the motor load, the change of σ n ( t ) is not obvious. The function σ n ( t ) tends to increase simultaneously with the fault severity. To this end, after an early slight fault is diagnosed under variable working conditions, the proposed method can be employed to observe whether the fault severity is deepened in advance.

Diagram of the time-varying fault diameter.

Diagram of instantaneous amplitude and singular value σ n ( t ) corresponding to Figure 11.

5. Conclusion

Rolling bearing fault diagnosis under time-varying working condition has always been a thorny issue in the field of bearing fault diagnosis. This paper presents a new fault diagnosis method combining HHT, SVD, and Elman neural network. HHT can generate the time domain signal and the frequency domain signal at the same time; thus, it supplies abundant information for the fault diagnosis under variable conditions. Since the time-frequency feature vector generated by the HHT is indeed redundant, the singular value decomposition is proposed in this paper for reducing dimensionality and obtaining feature insensitive to working condition variation, which is beneficial to enhance the classification accuracy under variable working conditions.

First, the IMFs were obtained by applying EMD. Second, the instantaneous amplitude matrices were calculated by applying Hilbert transform to each IMF. The matrix was then decomposed by SVD. This method finds out the singular value commonality for different operating conditions, which achieves rolling bearing fault diagnosis under time-varying conditions. Finally, the Elman neural network was used for fault identification and classification according to the extracted feature vector. The experimental results demonstrated the effectiveness of the singular value for different operating conditions.

Taking the existing EMD-SVD and WPT-PCA for comparison, we demonstrated the robustness of the proposed feature extraction based on HHT-SVD. Finally, by comparison to SVM and ELM, the effectiveness of the proposed method using the Elman neural network for classification was further verified. Based on the preceding discussion, the superiority of the proposed method in handling fault diagnosis under the time-varying condition can be verified. Without the need for human intervention and additional cost, this approach is an excellent automatic method that can be promoted for fault diagnosis under variable conditions.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research is supported by the National Natural Science Foundation of China (Grant nos. 51105019, 61074083, and 50705005) and the Technology Foundation Program of National Defense (Grant no. Z132013B002).

Allotta B. Papi R. Pugi L. Toni P. Violi A. G. Experimental campaign on a servo-actuated pantograph 1 Proceedings of the IEEE/ASME International Conference on Advanced Intelligent Mechatronics July 2001 Como, Italy 237 242 2-s2.0-0034866229 Shao Y. Ge L. Fang J. Fault diagnosis system based on smart bearing Proceedings of the International Conference on Control, Automation and Systems (ICCAS '08) July 2008 1084 1089 10.1109/ICCAS.2008.4694313 2-s2.0-58149096378 Yang Y. Wang H. Cheng J. Zhang K. A fault diagnosis approach for roller bearing based on VPMCD under variable speed condition Measurement 2013 46 8 2306 2312 10.1016/j.measurement.2013.03.023 2-s2.0-84879095999 He Q. Liu Y. Long Q. Wang J. Time-frequency manifold as a signature for machine health diagnosis IEEE Transactions on Instrumentation and Measurement 2012 61 5 1218 1230 10.1109/TIM.2012.2183402 2-s2.0-84862825547 Tang X. G. Guo Y. Zheng Y. C. Application of rolling element bearing envelope analysis based on short time Fourier transition and independent components analysis Journal of Mechanical Strength 2012 34 1 1 3 2-s2.0-84863392809 Shi L. Zhang Y. Mi W. Application of Wigner-Ville-distribution-based spectral kurtosis algorithm to fault diagnosis of rolling bearing Journal of Vibration, Measurement and Diagnosis 2011 31 1 27 31 2-s2.0-79953833442 Kankar P. K. Sharma S. C. Harsha S. P. Rolling element bearing fault diagnosis using wavelet transform Neurocomputing 2011 74 10 1638 1645 10.1016/j.neucom.2011.01.021 2-s2.0-79954415606 Feng Z. Liang M. Chu F. Recent advances in time-frequency analysis methods for machinery fault diagnosis: a review with application examples Mechanical Systems and Signal Processing 2013 38 1 165 205 10.1016/j.ymssp.2013.01.017 2-s2.0-84876940227 Huang N. E. Shen S. S. Hilbert-Huang Transform and Its Applications 2005 Singapore World Scientific MR2188500 Huang N. E. Shen Z. Long S. R. Wu M. C. Shih H. H. Zheng Q. Yen N. Tung C. C. Liu H. H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis Proceedings of the Royal Society of London Series A: Mathematical, Physical and Engineering Sciences 1998 454 1971 903 995 10.1098/rspa.1998.0193 MR1631591 2-s2.0-5444236478 Safonov M. G. Stability margins of diagonally perturbed multivariable feedback systems IEE Proceedings D Control Theory and Applications 1982 129 6 251 256 10.1049/ip-d.1982.0054 MR685110 2-s2.0-0020207749 Yang J. Zhang Y. Zhu Y. Intelligent fault diagnosis of rolling element bearing based on SVMs and fractal dimension Mechanical Systems and Signal Processing 2007 21 5 2012 2024 10.1016/j.ymssp.2006.10.005 2-s2.0-34047275789 Lei L. Zhang Q. Relevance Vector Machine based bearing fault diagnosis Proceedings of the IEEE International Conference on Machine Learning and Cybernetics August 2006 Dalian, China 3492 3496 10.1109/ICMLC.2006.258539 2-s2.0-33947214583 Huang G. B. Zhu Q. Y. Siew C. K. Extreme learning machine: a new learning scheme of feedforward neural networks Proceedings of the IEEE International Joint Conference on Neural Networks July 2004 985 990 10.1109/IJCNN.2004.1380068 2-s2.0-10944272650 Li H. Zhang Y. Bearing localized fault detection based on Hilbert-Huang transformation Proceedings of the 4th International Conference on Fuzzy Systems and Knowledge Discovery (FSKD '07) August 2007 138 142 10.1109/FSKD.2007.203 2-s2.0-44049098517 Yang Y. Yu D. Cheng J. A fault diagnosis approach for roller bearing based on IMF envelope spectrum and SVM Measurement 2007 40 9-10 943 950 10.1016/j.measurement.2006.10.010 2-s2.0-34848858238 Zhan R. Wu H. Feature extraction & application of engineering non-stationary signals based on EMD-AR model and SVD Proceedings of the IEEE International Conference on Computer Design and Applications (ICCDA '10) June 2010 V4-528 V4-532 10.1109/ICCDA.2010.5540704 2-s2.0-77955933942 Qi Z. Liu M. Wang H. Process modeling method based on an improved Elman Neural Network Proceedings of the 7th World Congress on Intelligent Control and Automation (WCICA '08) June 2008 Chongqing, China 8188 8192 10.1109/WCICA.2008.4594210 2-s2.0-52149104065 Cheng J. Yu D. Tang J. Yang Y. Application of SVM and SVD technique based on EMD to the fault diagnosis of the rotating machinery Shock and Vibration 2009 16 1 89 98 10.3233/SAV-2009-0457 2-s2.0-57649230902 Zhu Z. Sun Y. Application of EMD and SVD in fault identification Proceedings of the 8th International Power Engineering Conference (IPEC '07) April 2007 1247 1250 2-s2.0-51349142631 Shao R. Hu W. Wang Y. Qi X. The fault feature extraction and classification of gear using principal component analysis and kernel principal component analysis based on the wavelet packet transform Measurement 2014 54 118 132 10.1016/j.measurement.2014.04.016