Planetary gearbox fault diagnosis based on pseudofault signal assisted EMD
Xianglong Ni^{1} , Jianmin Zhao^{2} , Zhendong Du^{3} , Xin Zhang^{4}
^{1, 2, 3, 4}Mechanical Engineering College, Shijiazhuang, China
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 10, 2016, p. 138143.
Received 28 October 2016; accepted 30 October 2016; published 8 December 2016
Pseudofault signal assisted empirical mode decomposition (PFSEMD) is put forward for fault detection and isolation, which is built upon EMD, envelope analysis and pseudofault signal. And this method has been proved that it is usable in fault diagnosis of simple gear pair. However, this method is still need further study that whether it can be used for fault diagnosis of complex gearbox. Moreover, planetary gearbox fault diagnosis is much more complicated than fixedaxis gearbox. This paper makes a further study of PFSEMD and focuses on its application. Combining with the seeded failure experimental data of planetary gearbox, this paper show that PFSEMD can enhance fault information and weaken other useless information in complex gearbox fault diagnosis.
Keywords: planetary gearbox, fault diagnosis, pseudofault signal, EMD.
1. Introduction
Vibration is an important monitoring method of rotating machinery, and it is great significance to accurately extract the fault information from vibration signals [1]. Signal processing methods are usually used to analyze and extract the fault characteristic frequency, and they are becoming a common technology in rotating machinery condition monitoring process. Wavelet transform (WT) [2] and empirical mode decomposition (EMD) [3] are two widely used methods of signal processing. WT method needs to choose different wavelet basis functions to determine the wavelet coefficients in order to reconstruct the signal [2]. Compared with WT, the EMD method does not need basis function in the process of characteristic frequency extraction from vibration signal, which is a kind of adaptive nonlinear and nonstationary data analysis method [4]. The EMD method decomposes nonlinear and nonstationary signal into a series of intrinsic mode function (IMF). And the IMF is a kind of basically orthogonal, adaptive and complete expression, which is confirmed by original signal rather than preset basic functions [5]. EMD is widely used in processing nonlinear and nonstationary signal [6].
Dheeraj proposed a idea of injecting pseudofault signal to the signal envelope in fault detection and isolation of rotating machinery [7]. The method is built upon EMD, envelope analysis and pseudofault signal. He proved this method is usable in fault diagnosis of simple gear pair. However, it is still need further study that whether this method can be used for fault diagnosis of complex gearbox.
As the advantages of strong loadbearing capacity, large transmission ratio, etc., planetary gearboxes are widely used in modern industry. Planetary gearboxes significantly differ from fixedaxis gearboxes because of their unique structure [8]. As a result, planetary gearbox fault diagnosis are much more complicated than fixedaxis gearbox.
This paper makes a further study of PFSEMD and focuses on its application. Combining with the seeded failure experimental data of planetary gearbox, this paper verifies the effectiveness of PFSEMD in complex gearbox fault diagnosis.
2. Pseudofault signal theory
In PFSEMD method, the frequency of masking signal is predetermined by the known system fault frequency, thus it is considered as pseudofault signal. In order to clearly describe the concept of pattern recognition, the fault signal is assumed as two closer frequency band signal components. In this article, two signal components are expressed by the sine function:
where $h$ represents the higher frequency components, $l$ on behalf of the low frequency components.
When performing EMD (The procedure for performing the EMD method is shown in Fig. 1), the first step in the sifting process is to find the extreme value point of signal and create upper and lower envelope. Moreover, the extreme value distribution depends on the signal component. Pseudofault signal is used to change the distribution of extreme value. As a result, the average envelope got by EMD should contain only low frequency components:
Fig. 1. Procedure for performing the EMD
2.1. The positive fault IMF
The positive synthesized signal ${Y}^{+}\left(t\right)$ is obtained from envelope signal ${Y}_{eop}\left(t\right)$ (containing fault information) and pseudofault signal ${Y}_{PFS}\left(t\right)$. EMD Processing ${Y}^{+}\left(t\right)$ and achieving the positive IMF ${X}^{+}\left(t\right)$:
In pseudofault signal ${Y}_{PFS}\left(t\right)$, in order to make the mean extrema rate of ${Y}^{+}\left(t\right)$ is equal to ${Y}_{PFS}\left(t\right)$ (namely ${f}_{ex}={f}_{p}$), ${B}_{P}$ should be big enough. Let $\left\{{t}_{i}^{u},{Y}^{+}\left({t}_{i}^{u}\right)\right\}$ as the set of local maxima of ${Y}^{+}\left(t\right)$, and ${Y}^{+}\left({t}_{i}^{u}\right)$ can be written as:
The ideal curve, which connect all of maxima value (namely upper envelope), can be expressed as:
where ${T}_{ex}$ is the sampling interval of extrema.
Because of sampling for the extreme value, the higher frequency components of signals will loss information as aliasing, the lower frequency components will be reconstructed by the means of minimum information loss. As the assumption ${f}_{ex}={f}_{p}$, the pseudofault signal will produce a constant DC term. Therefore, the upper envelope can be expressed as:
where ${B}_{P}^{u}$ is a constant, ${Y}_{h}^{u}\left(t\right)$ is the edition of higher frequency components after aliasing. Similarly, the lower envelope can be written as:
The average envelope can be obtained as:
In the Eq. (9), it is the lower frequency component in the first term of right hand expression, and it is average envelope for the higher frequency components after aliasing. Averaging operation further weaken the influence of higher frequency components by aliasing envelope [9]. As a result, the positive IMF ${X}^{+}\left(t\right)$ can be calculated as:
where ${\zeta}_{1}\left(t\right)$ is the error signal.
2.2. The negative fault IMF
The negative fault IMF ${X}^{}\left(t\right)$ is extracted through the similar series of operation on negative synthesized signal ${Y}^{}\left(t\right)$:
Ignoring the error terms in Eqs. (10) and (11), the output signal after PFSEMD operation can be achieved as:
Fig. 2. Procedure for performing the PFSEMD
It can be found from Eq. (12) that the pseudofault signal can help in separating signals of closer frequency bands, which is out of reach for standard EMD method. The signal dealing process of PFSEMD method are shown in Fig. 2.
3. PFSEMD method
The proposed PFSEMD method is used to extract the fault characteristics of rotating machinery, the specific steps are as follows:
(a) Calculating the failure frequency ${F}_{j}$ ($j=$ 1, 2,..., $m$) corresponding to various failure modes.
(b) Identifying the $i$th IMF that there is a peak close to ${F}_{j}$ in frequency spectrum:
where ${C}_{i}\left(F\right)$ is Fourier transform of the $i$th IMF ${C}_{i}\left(t\right)$, $\epsilon $ is a very small value and it is taken as 0.05 in this paper.
(c) Structuring a pseudofault signal, which with constant amplitude and constant frequency:
On the advice of previous research [7], this paper let ${\alpha}_{F}=$ 1.1 and ${\alpha}_{A}=$ 1.6.
(d) Performing standard EMD to get the first IMF ${X}^{+}\left(t\right)$ from ${C}_{i}^{+}\left(t\right)={C}_{i}\left(t\right)+{Y}_{PFS}\left(t\right)$ and get the second IMF ${X}^{}\left(t\right)$ from ${C}_{i}^{}\left(t\right)={C}_{i}\left(t\right){Y}_{PFS}\left(t\right)$.
(e) Obtaining fault indicator signal:
(f) Extracting fault information from fault indicator signal $X\left(t\right)$.
4. Planetary gearbox fault diagnosis
A case study is carried out for planetary gearbox fault diagnosis based on PFSEMD. In order to validate the effectiveness of the proposed method, a seeded failure experiment of planetary gearbox to obtain failure data. The experiment rig consists of a test planetary gearbox, a drive motor, a speed and torque sensor, and a magnetic powder brake, as shown in Fig. 3. There are four accelerometers fitted onto the casing of gearbox to record vibration data, as shown in Fig. 4. Schematic map of planetary gearbox structure and related configuration parameters are shown in Fig. 5. The failure is seeded on sun gear, planet gear and ring gear, respectively, as shown in Fig. 6. Fault characteristic frequencies of every gear is shown in Table 1.
Fig. 3. Planetary gearbox experiment rig
Frequency spectrum of sun gear fault is shown as Fig. 7, blue frequency spectrum is decomposed by standard EMD without pseudofault signal, red frequency spectrum is decomposed by PFSEMD. It can be found that the amplitude of red frequency spectrum is lower than blue frequency spectrum in whole frequency bands. But in fault frequency of sun gear (${f}_{S}=$ 56.97), the amplitude of red signal is significantly higher than blue signal. This phenomenon shows that fault information has been enhanced, while other useless information has been weakened.
Table 1. Fault characteristic frequencies of every gear (Hz)
Fault frequency of sun gear ${f}_{S}$

Fault frequency of planet gear ${f}_{P}$

Fault frequency of ring gear ${f}_{R}$

55.093

3.730

4.906

Fig. 4. Mounted location of accelerometers
Fig. 5. Schematic map of planetary gearbox structure
Fig. 6. Implanted wear fault: a) sun gear, b) planet gear, and c) ring gear
Fig. 7. Frequency spectrum of sun gear fault
Frequency spectrums of planet gear fault and ring gear fault are shown as Fig. 8. It cannot found the fault characteristic frequency information in the frequency spectrums, which are decomposed by standard EMD without pseudofault signal. While in the frequency spectrums which are decomposed by PFSEMD, there are two significant peaks in fault characteristic frequencies of planet gear fault and ring gear fault, the fault features are obvious.
In order to further verify the effect of the proposed method, signal to noise ratios (SNR) are analyzed. The SNR can be caudated as follow:
$S\left({f}_{0}\right)$ stand for spectrum amplitude of frequency ${f}_{0}$, $N\left({f}_{0}\right)$ is average spectrum amplitude. All SNR of every gear fault is shown in Table 2. It can be found that the SNR of PFSEMD is much greater than EMD without PFS. It shows that PFSEMD can enhance useful information in fault diagnosis.
Fig. 8. Frequency spectrum: a) planet gear fault, b) ring gear fault
a)
b)
Table 2. SNR of every gear fault
SNR of sun gear fault

SNR of planet gear fault

SNR of ring gear fault


EMD without PFS

10.2976

3.4278

2.3679

EMD with PFS

15.5547

8.4135

8.0443

5. Conclusions
This paper is meant to further investigate the theory and application in complex gearbox of PFSEMD. Pseudofault signal theory is future studied and the specific steps of PFSEMD are given. A case study of planetary gearbox fault diagnosis based on PFSEMD is analysis, which include sun gear fault, planet gear fault and ring gear fault. The frequency spectrum and SNR show that PFSEMD can enhance fault information and weaken other useless information in complex gearbox fault diagnosis.
References
 Li M., Xu J., Yang J., Yang D. Multiple manifolds analysis and its application to fault diagnosis. Mechanical Systems and Signal Processing, Vol. 23, 2009, p. 25002509. [Search CrossRef]
 Gilles J. Empirical wavelet transform. IEEE Transactions on Signal Processing, Vol. 61, 2013, p. 39994010. [Search CrossRef]
 Lei Y. G., Lin J., He Z. J., Zuo M. J. A review on empirical mode decomposition in fault diagnosis of rotating machinery. Mechanical Systems and Signal Processing, Vol. 35, 2013, p. 108126. [Search CrossRef]
 Bi N., Sun Q. Y., Huang D. R., Yang Z. H., Huang J. W. Robust image watermarking based on multiband wavelets and empirical mode decomposition. IEEE Transactions on Image Processing, Vol. 16, 2007, p. 19561966. [Search CrossRef]
 Lv Y., Yuan R., Song G. B. Multivariate empirical mode decomposition and its application to fault diagnosis of rolling bearing. Mechanical Systems and Signal Processing, Vol. 81, 2016, p. 219234. [Search CrossRef]
 Klionski D. M., Oreshko N. I., Geppener V. V. Applications of empirical mode decomposition for processing nonstationary signal. Pattern Recognition and Image Analysis, Vol. 18, 2008, p. 390399. [Search CrossRef]
 Singh D. S., Zhao Q. Pseudofault signal assisted EMD for fault detection and isolation in rotating machines. Mechanical Systems and Signal Processing, Vol. 81, 2016, p. 202218. [Search CrossRef]
 Lei Y. G., Han D., Lin J., et al. Planetary gearbox fault diagnosis using an adaptive stochastic resonance method. Mechanical Systems and Signal Processing, Vol. 38, 2013, p. 113124. [Search CrossRef]
 Yang Y., Miao C., Deng J. Research on the mechanism of masking signals technique to improve empirical mode decomposition. Journal of Information and Computational Science, Vol. 9, 2012, p. 35093516. [Search CrossRef]