# Rolling Bearing Fault Diagnosis under Variable Conditions Using Hilbert-Huang Transform and Singular Value Decomposition

Authors: Hongmei Liu; Xuan Wang; Chen Lu
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

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.


---

## Body

## 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 [1] 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 [2]. Yang et al. introduced variable predictive model based class discrimination (VPMCD) into rolling bearing fault diagnosis [3]. 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 [4]. 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 [5], the Wigner-Ville distribution [6], and the wavelet transform [7] 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 [8], 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) [9] is free from the restriction of Heisenberg’s uncertainty principle. Its essence is processing the signal with Empirical Mode Decomposition (EMD) [10] 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 [11] 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) [12], Relevance Vector Machine (RVM) [13], and Extreme Learning Machine [14]. 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. Section2 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 [10], 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 ith 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 toc
(
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 signalx
(
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 [15], 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 [16]. 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 thatM 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 [8] 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 [17]. 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 Figure1, 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
)
.Figure 1
Flowchart of obtaining the feature vector.After the application of SVD toM, 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
).


## 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 [10], 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 ith 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 toc
(
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 signalx
(
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 [15], 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 [16]. 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 thatM 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 [8] 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 [17]. 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 Figure1, 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
)
.Figure 1
Flowchart of obtaining the feature vector.After the application of SVD toM, 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 [18]. 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.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 Figure3, 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.Figure 3
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 Figure4(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.
(a)
(b)
(c)
(d)After obtaining the instantaneous amplitude matrices, we can obtain the singular value vector by conducting SVD. As shown in Table1, 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.Table 1
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
24Figure 5
Scatter plots of the first three singular values extracted by HHT-SVD for datasets in Table1.


### 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 Table1, 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.Figure 6
Scatter plots of the first three singular values extracted by EMD-SVD for datasets in Table1.Afterwards, another widely used method, WPT in tandem with PCA [21], 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.Figure 7
Scatter plots of the first three PCs extracted by WPT and PCA for datasets in Table1.


### 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 Table2, 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.Table 2
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
5Table 3
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.9815Figure 8
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 Table3. 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 Figure9. 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
).Figure 9
Diagram of the time-varying condition.Figure 10
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 Figure11. 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.Figure 11
Diagram of the time-varying fault diameter.Figure 12
Diagram of instantaneous amplitude and singular valueσ
n
(
t
) corresponding to Figure 11.


## 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 Figure4(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.
(a)
(b)
(c)
(d)After obtaining the instantaneous amplitude matrices, we can obtain the singular value vector by conducting SVD. As shown in Table1, 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.Table 1
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
24Figure 5
Scatter plots of the first three singular values extracted by HHT-SVD for datasets in Table1.


## 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 Table1, 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.Figure 6
Scatter plots of the first three singular values extracted by EMD-SVD for datasets in Table1.Afterwards, another widely used method, WPT in tandem with PCA [21], 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.Figure 7
Scatter plots of the first three PCs extracted by WPT and PCA for datasets in Table1.


## 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 Table2, 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.Table 2
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
5Table 3
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.9815Figure 8
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 Table3. 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 Figure9. 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
).Figure 9
Diagram of the time-varying condition.Figure 10
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 Figure11. 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.Figure 11
Diagram of the time-varying fault diameter.Figure 12
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.



---

*Source: 765621-2014-07-16.xml*