# Planetary Gearbox Fault Diagnosis Using Envelope Manifold Demodulation

Authors: Weigang Wen; Robert X. Gao; Weidong Cheng
Journal: Shock and Vibration
 (2016)

---

## Abstract

The important issue in planetary gear fault diagnosis is to extract the dependable fault characteristics from the noisy vibration signal of planetary gearbox. To address this critical problem, an envelope manifold demodulation method is proposed for planetary gear fault detection in the paper. This method combines complex wavelet, manifold learning, and frequency spectrogram to implement planetary gear fault characteristic extraction. The vibration signal of planetary gear is demodulated by wavelet enveloping. The envelope energy is adopted as an indicator to select meshing frequency band. Manifold learning is utilized to reduce the effect of noise within meshing frequency band. The fault characteristic frequency of the planetary gear is shown by spectrogram. The planetary gearbox model and test rig are established and experiments with planet gear faults are conducted for verification. All results of experiment analysis demonstrate its effectiveness and reliability.


---

## Body

## 1. Introduction

Planetary gear systems are widely used in industrial machinery because they have advantages of large transmission ratios, strong load-bearing capacity, and high transmission efficiency [1]. In practical industrial fields, the planetary gear system whose structure is more complicated than fixed-shaft gearbox is subject to be damaged under heavy load and tough working conditions. Accurate fault detection of planetary gearbox is important to reduce unscheduled machine downtime and avoid catastrophic accidents [2]. Vibration signal analysis which was used extensively in various machine condition monitoring techniques has become one of the most important methods applied for planetary gearbox fault diagnostics.The planetary gear system is composed of sun gear, planet gears, ring gear, planet carrier, shaft, and bearings. Usually, the ring gear is stationary, the sun gear rotates around its own center, and the planet carrier rotates under very low speeds because of the large transmission ratio. Multiple planet gears, which mesh simultaneously with the sun gear and ring gear, not only rotate around their own center but also revolve around the center of the sun gear [3]. Such behaviors cause the vibration signals of the planetary gearbox with different meshing phases and varied propagation paths [4]. The unique behavior mode of planetary gearbox presents a challenge to fault diagnosis of planetary gear system. Studies on the fault diagnosis of planetary gearboxes are gaining more attention of the researchers [5].The fault diagnosis of planetary gearboxes comes from various methodologies. Modeling methods are useful in condition monitoring. Tooth pitting and crack on the sun gear of planetary gearboxes were modeled and the influence of faults on the gear mesh stiffness was analyzed in [6]. An approach based on physical models to detect damage of a planetary gearbox in a helicopter transmission system was proposed in [7]. The spectral structure of vibration signal of planetary gearbox in different health conditions was investigated and vibration signal models of gear damage were proposed in [8]. But the most popular fault diagnosis of planetary gearbox is signal processing methods, for example, time-domain methods, frequency-domain methods, time-frequency-domain methods, and other signal processing methods. Time-domain signal processing methods, such as statistical indicators and time synchronous averaging, are relatively easy and direct compared to other methods. A technique for estimating the time-domain average of the tooth meshing vibration of individual planet gears and of the sun gear was demonstrated in [9]. Two special diagnostic parameters in time-domain for detecting faults of planetary gearbox were studied in [10]. In frequency-domain, the Fourier series of the vibration data from a planetary gearbox was explored and the source of the asymmetry phenomena in the frequency spectra was analyzed in [1]. The additional sidebands in frequency spectra produced by planet carrier torque modulations masking potentially the sidebands of fault in planetary gearboxes were investigated and a frequency-domain method to eliminate the effects of structural-path-caused amplitude changes was presented in [11]. The time-frequency-domain methods, such as Wigner-Ville distribution and wavelets, are developed to diagnose planetary gearboxes. The Wigner-Ville distribution of the simulated signals of fault modes in planetary gearboxes was studied in [12]. A method based on adaptive Morlet wavelets and singular value decomposition was applied to extract impulse features of a planetary gearbox in [13]. Reference [14] illustrated the spectral kurtosis for detecting a tooth crack on the ring in the planetary gearbox of a wind turbine. The cyclostationary properties of planetary gearbox vibration signals were investigated and a fault diagnosis based on spectral coherence map was developed in [15]. These literatures demonstrated the different effects of signal processing methods. The time-domain and frequency-domain methods just depict one aspect of the signal. The characteristic in time-domain and frequency-domain cannot be drawn simultaneously. Time-frequency-domain methods are generally much more effective than the other two. But it is better to combine them with other methods to implement demodulating and denoising.The vibration signals of mechanical systems, especially for ones with fault, often show mutation, nonlinearity, and nonstationarity [16]. With the development of nonlinear dynamic theories, a number of nonlinear parameters and methods have been applied to machine condition monitoring and fault diagnosis [17–19]. Manifold learning especially is adopted to implement nonlinear dimensionality reduction for revealing the nonlinear intrinsic characteristics of vibration signal. The signal processing technique and manifold learning method are combined to strengthen the fault features of the rotating machine and reduce the impact of noise [20]. The manifold learning algorithm on the time-frequency distributions or the time-scale distributions generated manifold signature of fault to enhance nonlinear fault-induced components and suppress in-band noise in [21–23].In planetary gearbox, the vibration signals from multiple planet gears may couple with each other leading the faulty components to be neutralized and cancelled [24]. Especially in practical industrial applications, the vibration signal is contaminated by background noise. Here, in order to extract the fault feature of planetary gear from the noisy vibration signal, the methodology implementing decoupling, demodulating, and denoising is investigated. Because of the demodulation capability of wavelet envelope and the denoising ability of manifold learning [20, 25], these methods are integrated to achieve demodulating and denoising of the planetary gear vibration signal simultaneously in this paper.The rest of the paper is organized as follows. In the second section, the wavelet envelope manifold demodulation method is introduced for fault diagnosis of planetary gearbox. Then in the third section, the vibration models and simulation experiments of the planetary gearbox in faulty condition are studied. In the fourth section, the practical vibration signals from planetary gearbox rig are analyzed to testify the effectiveness of the proposed method. The conclusions are finally presented in the fifth section.


## 2. Principle and Methodology

### 2.1. Framework of Methodology

In fault diagnosis of planetary gears, the accelerometers are fixed on the gearbox casing to collect vibration signals. The relative locations of sun-planet and planet-ring meshing with respect to the accelerometers are time-varying, which result in amplitude modulation (AM) effect on the meshing vibration signal. Time variation of meshing stiffness or the errors of the planetary gear also causes frequency modulation (FM) of the meshing vibration signal. The gear damage generates both AM and FM of the meshing vibration signal. These effects increase the complexity of the signal. Furthermore, there exists the random background noise in every frequency band for practical fault diagnosis of planetary gear in industrial fields. The noise can blight the signal at meshing frequency. So, extracting the meshing signal, reducing the noise effect in meshing frequency band, and demodulating the meshing vibration signal are the critical issues for fault diagnosis of planetary gearbox.The enveloping analysis of resonance is a powerful tool to achieve demodulation of the signal of rotating machine. In the traditional way, a filter with resonance frequency and specific pass band is utilized to preprocess the vibration signal, and then the envelope of the resonance is drawn out by Hilbert transform. Here, a synthetic methodology that combines the complex wavelet transform, manifold learning method, and spectral analysis is proposed to demodulate and denoise the meshing vibration signal of planetary gears. The methodology is shown in Figure1. At first, the continuous wavelet transform (CWT) with complex Morlet wavelet base is utilized to demodulate the vibration signal over scales to obtain the wavelet envelopes. Then, the meshing vibration sensitive indicator is built to select the specific scale band which corresponds to meshing frequency band. The first two steps both extract the envelope and filter out the noise not within the meshing frequency band. Thirdly, the manifold learning algorithm is introduced to extract the envelope manifold from the wavelet envelopes at the selected scales to reduce the in-band noise effect and reveal the envelope structure of the meshing vibration signal. Finally, the fault characteristic frequency is discovered by frequency spectrogram analysis.Figure 1
The process diagram of the methodology.


### 2.2. Complex Wavelet Enveloping

It is well known that the wavelet transform can decompose a signal onto a time-scale plane [25]. The wavelet transform of a signal is equal to a set of band-pass filters with certain frequencies relevant to the scales. Continuous wavelet transform (CWT) is defined as the convolution of signal x(t) with the complex conjugate of a set of wavelets. Here, complex wavelet φ(t) is applied. The complex wavelet has the analytic property. It is defined as(1)φt=φRt+jφIt=φRt+jHφRt,where φRt and φI(t) represent the real and imaginary parts of the complex wavelet and H[] is the Hilbert transform. Usually, the complex Morlet wavelet is used as the base wavelet for the vibration signal of rotating machine. The complex Morlet wavelet is defined as(2)φt=1πfbe-t2/fbej2πfct,where fb is the bandwidth parameter and fc is the wavelet center frequency. Here, the parameter optimization of the Morlet wavelet is not the research focus of this paper. The parameters are set as fb=1 and fc=2 by which the satisfied results can be achieved. The base wavelet with specific scale s is used to extract the component of frequency f in the signal. The relationship can be written as(3)f=fcsΔ,where Δ is the sampling period. From formula (3), the meshing frequency of the planetary gears can be estimated according to the special indicator at the scale. The wavelet transform of the signal based on complex wavelet can be written as(4)WCs,τ=WRs,τ+jWIs,τ=WRs,τ+jHWRs,τ.According to the property of Hilbert transform, the envelope of signal at scales can be computed from the modulus of the wavelet coefficients as(5)EWs,τ=WCs,τ=WR2s,τ+WI2s,τ.Through complex wavelet transformation and scale band selection, the band-pass filtering and enveloping can be achieved simultaneously.


### 2.3. Meshing Scale Band Selection

In the enveloping technique for fault diagnosis of rotating machine, the critical issue is to estimate the resonance frequency of the structure. In current literatures different indicators have been applied for decision of resonance modes. The spectral kurtosis (SK) is an effective measure to estimate the resonance frequency of impulsive components in vibration signal [26–29]. The smoothness index (SI) defined as the ratio of the arithmetic mean to the geometric mean is also utilized to measure the impulsiveness of the power spectrum of wavelet coefficient at each scale in [30]. After complex wavelet filtering and enveloping, the energy of vibration signal envelope concentrated within meshing frequency band more than the other band. The mean envelope energy can be adopted as the indicator for selection of meshing scale band. For the wavelet envelope of the signal x(t), the mean energy (ME) of wavelet envelope at scale s is defined as(6)MEs=∑τEWs,τEW∗s,τN,where EW∗s,τ is the complex conjugate of EWs,τ and N is number of sample points. When the ME, SK, and SI in frequency-domain are considered as indicators for selecting the appropriate scale for the meshing frequency in the data simulations and practical experiments, it is revealed that the ME is most suitable to select scale corresponding to the meshing frequency.The scale corresponding to the global maximum of the indicator is defined as central scalesC, which is utilized to evaluate the planetary gear meshing central frequency. The wavelet envelopes in meshing scale band, which construct a high-dimensional data space, contain the fault information in the meshing frequency band. The meshing scale band can be divided by the left limit sL and right limit sR around the central scale sC. The left limit sL and right limit sR can be determined by the indicator value reducing to 1/2 of the maximum value of the indicator at the central scale sC. So the scale band for the meshing frequency band can be defined as sB=[sL,…,sC,…,sR].The fault-induced impacts modulating the meshing frequency form the manifold structure embedded in the wavelet envelopes at this scale band. There is also another background noise within meshing frequency band. The in-band noise will corrupt the embedded manifold structure. Here, manifold learning is applied to retain the envelope manifold of the fault impacts while leaving out the in-band noise. And then-dimension (n=R-L+1, where R is the index of sR and L is the index of sL) data space can be reduced tom-dimension (m≪n) by mean of the nonlinear dimensionality reduction for posttreatment.


### 2.4. Manifold Learning

There are several algorithms of manifold learning. The techniques of Isomap and locally linear embedding (LLE) were firstly proposed. Isomap is a global method because its embedding is based upon the geodesic distances between all pairs of points, while LLE is a local method because its embedding is based upon the relationship of each data point to its neighboring points [31, 32]. Laplacian eigenmap, Hessian eigenmap, local tangent alignment, diffusion map, and several other techniques were introduced.A local approach is computationally more efficient than a global approach because the former involves sparse matrices and is more generally applicable to different types of manifolds. Laplacian eigenmap, which is a typical local approach, is applied to extract the manifold structure from the wavelet envelopes at planetary gear meshing frequency band because of its efficiency and robustness to the noise [33].The methodology of Laplacian eigenmap algorithm consists of three steps: (1) constructK nearest neighbor graph on the set of sample points; (2) approximate the local manifold geometry within the neighborhood of each sample point; (3) minimize a global error function to obtain the global embedding [34].Here, the number of nearest neighbor is set asK=40, the reduced dimension is set as m=2, and the first dimension data are taken as the signal for the next step of frequency spectrogram analysis for fault characteristic frequency extraction because of its minimum reconstruction error.


## 2.1. Framework of Methodology

In fault diagnosis of planetary gears, the accelerometers are fixed on the gearbox casing to collect vibration signals. The relative locations of sun-planet and planet-ring meshing with respect to the accelerometers are time-varying, which result in amplitude modulation (AM) effect on the meshing vibration signal. Time variation of meshing stiffness or the errors of the planetary gear also causes frequency modulation (FM) of the meshing vibration signal. The gear damage generates both AM and FM of the meshing vibration signal. These effects increase the complexity of the signal. Furthermore, there exists the random background noise in every frequency band for practical fault diagnosis of planetary gear in industrial fields. The noise can blight the signal at meshing frequency. So, extracting the meshing signal, reducing the noise effect in meshing frequency band, and demodulating the meshing vibration signal are the critical issues for fault diagnosis of planetary gearbox.The enveloping analysis of resonance is a powerful tool to achieve demodulation of the signal of rotating machine. In the traditional way, a filter with resonance frequency and specific pass band is utilized to preprocess the vibration signal, and then the envelope of the resonance is drawn out by Hilbert transform. Here, a synthetic methodology that combines the complex wavelet transform, manifold learning method, and spectral analysis is proposed to demodulate and denoise the meshing vibration signal of planetary gears. The methodology is shown in Figure1. At first, the continuous wavelet transform (CWT) with complex Morlet wavelet base is utilized to demodulate the vibration signal over scales to obtain the wavelet envelopes. Then, the meshing vibration sensitive indicator is built to select the specific scale band which corresponds to meshing frequency band. The first two steps both extract the envelope and filter out the noise not within the meshing frequency band. Thirdly, the manifold learning algorithm is introduced to extract the envelope manifold from the wavelet envelopes at the selected scales to reduce the in-band noise effect and reveal the envelope structure of the meshing vibration signal. Finally, the fault characteristic frequency is discovered by frequency spectrogram analysis.Figure 1
The process diagram of the methodology.


## 2.2. Complex Wavelet Enveloping

It is well known that the wavelet transform can decompose a signal onto a time-scale plane [25]. The wavelet transform of a signal is equal to a set of band-pass filters with certain frequencies relevant to the scales. Continuous wavelet transform (CWT) is defined as the convolution of signal x(t) with the complex conjugate of a set of wavelets. Here, complex wavelet φ(t) is applied. The complex wavelet has the analytic property. It is defined as(1)φt=φRt+jφIt=φRt+jHφRt,where φRt and φI(t) represent the real and imaginary parts of the complex wavelet and H[] is the Hilbert transform. Usually, the complex Morlet wavelet is used as the base wavelet for the vibration signal of rotating machine. The complex Morlet wavelet is defined as(2)φt=1πfbe-t2/fbej2πfct,where fb is the bandwidth parameter and fc is the wavelet center frequency. Here, the parameter optimization of the Morlet wavelet is not the research focus of this paper. The parameters are set as fb=1 and fc=2 by which the satisfied results can be achieved. The base wavelet with specific scale s is used to extract the component of frequency f in the signal. The relationship can be written as(3)f=fcsΔ,where Δ is the sampling period. From formula (3), the meshing frequency of the planetary gears can be estimated according to the special indicator at the scale. The wavelet transform of the signal based on complex wavelet can be written as(4)WCs,τ=WRs,τ+jWIs,τ=WRs,τ+jHWRs,τ.According to the property of Hilbert transform, the envelope of signal at scales can be computed from the modulus of the wavelet coefficients as(5)EWs,τ=WCs,τ=WR2s,τ+WI2s,τ.Through complex wavelet transformation and scale band selection, the band-pass filtering and enveloping can be achieved simultaneously.


## 2.3. Meshing Scale Band Selection

In the enveloping technique for fault diagnosis of rotating machine, the critical issue is to estimate the resonance frequency of the structure. In current literatures different indicators have been applied for decision of resonance modes. The spectral kurtosis (SK) is an effective measure to estimate the resonance frequency of impulsive components in vibration signal [26–29]. The smoothness index (SI) defined as the ratio of the arithmetic mean to the geometric mean is also utilized to measure the impulsiveness of the power spectrum of wavelet coefficient at each scale in [30]. After complex wavelet filtering and enveloping, the energy of vibration signal envelope concentrated within meshing frequency band more than the other band. The mean envelope energy can be adopted as the indicator for selection of meshing scale band. For the wavelet envelope of the signal x(t), the mean energy (ME) of wavelet envelope at scale s is defined as(6)MEs=∑τEWs,τEW∗s,τN,where EW∗s,τ is the complex conjugate of EWs,τ and N is number of sample points. When the ME, SK, and SI in frequency-domain are considered as indicators for selecting the appropriate scale for the meshing frequency in the data simulations and practical experiments, it is revealed that the ME is most suitable to select scale corresponding to the meshing frequency.The scale corresponding to the global maximum of the indicator is defined as central scalesC, which is utilized to evaluate the planetary gear meshing central frequency. The wavelet envelopes in meshing scale band, which construct a high-dimensional data space, contain the fault information in the meshing frequency band. The meshing scale band can be divided by the left limit sL and right limit sR around the central scale sC. The left limit sL and right limit sR can be determined by the indicator value reducing to 1/2 of the maximum value of the indicator at the central scale sC. So the scale band for the meshing frequency band can be defined as sB=[sL,…,sC,…,sR].The fault-induced impacts modulating the meshing frequency form the manifold structure embedded in the wavelet envelopes at this scale band. There is also another background noise within meshing frequency band. The in-band noise will corrupt the embedded manifold structure. Here, manifold learning is applied to retain the envelope manifold of the fault impacts while leaving out the in-band noise. And then-dimension (n=R-L+1, where R is the index of sR and L is the index of sL) data space can be reduced tom-dimension (m≪n) by mean of the nonlinear dimensionality reduction for posttreatment.


## 2.4. Manifold Learning

There are several algorithms of manifold learning. The techniques of Isomap and locally linear embedding (LLE) were firstly proposed. Isomap is a global method because its embedding is based upon the geodesic distances between all pairs of points, while LLE is a local method because its embedding is based upon the relationship of each data point to its neighboring points [31, 32]. Laplacian eigenmap, Hessian eigenmap, local tangent alignment, diffusion map, and several other techniques were introduced.A local approach is computationally more efficient than a global approach because the former involves sparse matrices and is more generally applicable to different types of manifolds. Laplacian eigenmap, which is a typical local approach, is applied to extract the manifold structure from the wavelet envelopes at planetary gear meshing frequency band because of its efficiency and robustness to the noise [33].The methodology of Laplacian eigenmap algorithm consists of three steps: (1) constructK nearest neighbor graph on the set of sample points; (2) approximate the local manifold geometry within the neighborhood of each sample point; (3) minimize a global error function to obtain the global embedding [34].Here, the number of nearest neighbor is set asK=40, the reduced dimension is set as m=2, and the first dimension data are taken as the signal for the next step of frequency spectrogram analysis for fault characteristic frequency extraction because of its minimum reconstruction error.


## 3. Simulated Signal Analysis

### 3.1. Simulation of Faulty Planet Gear

The planet gears mesh with ring gear and sun gear simultaneously. When there are faults on a planet gear, the meshing location of the faulty tooth and mating gear changes with the planet carrier rotation. Thus this changes the transfer path of the vibration signal of faulty tooth to the sensor and generates an extra AM effect on the meshing vibration signal. The meshing vibration amplitudes sampled by the sensor are modulated by the faulty tooth and the varied propagations synchronously [8]. Here, the faulty planet gear is analyzed. The faulty planet gear vibration signal can be written as(7)xt=1-Ccos⁡2πfcarriert1+Acos⁡2πfplanett+∅cos⁡2πfmesht+Bsin⁡2πfplanett+φ+θ,where xt is vibration signal in time-domain, fcarrier is planet carrier rotation frequency which represents propagation modulation, fplanet is planet gear rotation frequency which represents faulty planet gear modulation, fmesh is meshing frequency, and ∅, φ, and θ are the initial phases of AM, FM, and meshing, respectively. A, B, and C are the modulation amplitudes. After amplitude demodulation, the envelope of the meshing vibration signal is(8)yt=1-Ccos⁡2πfcarriert+Acos⁡2πfplanett+∅+12ACcos⁡2πfplanet+fcarriert+∅+12ACcos⁡2πfplanet-fcarriert+∅.So the envelope amplitude spectrum is(9)Yf=δf+Cδf-fcarrier+Aδf-fplanet+12ACδf-fplanet+fcarrier+12ACδf-fplanet-fcarrier.The simulation parameters of the planetary gearbox are shown in Table1. For comparison, the time-domain diagram and frequency spectrogram are adopted to analyze the vibration signals at first. Then the envelope manifold demodulation is utilized to illustrate the advantage of this method.Table 1
The simulation parameters of faulty planet gear.
Parameters
Values
AM magnitudeA
1
AM magnitudeC
0.8
FM magnitudeB
0.5
Meshing frequency (Hz)
437.5
Planet carrier rotation frequency (Hz)
5.468
Planet gear rotation frequency (Hz)
15.191The simulated vibration signal in time-domain is shown in Figure2(a). The frequency spectrogram of the signal is shown in Figure 2(b). The meshing frequency of 437.5 Hz and the sideband can be seen in the spectrogram. It is not very clear to identify the planet carrier rotation frequency and the planet gear rotation frequency.Figure 2
The time-domain diagram and frequency spectrogram of simulated signal.
(a)
The time-domain diagram
(b)
The frequency spectrogramIn envelope manifold demodulation, a series of scales of wavelet ranging from 1 to 30 are chosen to cover the meshing frequency of all the planetary gears. An increment of 0.1 for scale is employed in the data experiments. The wavelet envelopes are shown in Figure3. The lighter color in the graph illustrates greater wavelet envelopes. So the area circled by the ellipse corresponds to the greatest wavelet envelopes. The meshing frequency computed by formula (3) is the central frequency of the circled area.Figure 3
The wavelet envelope of simulated signal.The mean energy (ME), smoothness index (SI), and spectral kurtosis (SK) of wavelet envelope are utilized to choose the proper indictor of meshing frequency for comparison. The ME, SI, and SK as indicators are shown in Figures4(a), 4(b), and 4(c). The central scale, respectively, is 22.8, 22.2, and 29.0 according to the ME, SI, and SK. The corresponding frequency, respectively, is 438.5 Hz, 450.5 Hz, and 344.8 Hz. The result from ME is closest to the meshing frequency band which is circled by the ellipse in Figure 3. The indicator ME is most accurate for selection of meshing frequency band.Figure 4
The indicators versus scale curve of simulated signal.
(a)
The ME versus scale curve
(b)
The SI versus scale curve
(c)
The SK versus scale curveThe wavelet envelope at the central scale is shown in Figure5(a) and the power spectrum of the wavelet envelope is shown in Figure 5(b). The characteristic frequencies of the power spectrum are 4.89 Hz, 9.78 Hz, 14.66 Hz, and 19.55 Hz corresponding to fcarrier, fplanet-fcarrier, fplanet, and fplanet+fcarrier. The characteristic frequency of faulty planet gear can be seen clearly on the spectrogram of the wavelet envelope because there is no in-band noise interference. The wavelet envelope is proved to be adept at demodulation of planetary gear vibration signal.Figure 5
The wavelet envelope and spectrogram at the central scale of simulated signal.
(a)
The wavelet envelope
(b)
The envelope frequency spectrogram


### 3.2. Simulation of Contaminated Signal

In order to simulate the contamination of the in-band noise, the random background noise is added to the simulated signal of the faulty planet gear model to generate the contaminated signal up to SNR= −10 dB. Meanwhile, the meshing frequency is switched from 437.5 Hz to 492.0 Hz to testify the accuracy of selection of meshing frequency by ME.For comparison, the amplitude demodulation methods based on empirical mode decomposition (EMD) and Hilbert transform are used to reveal the characteristic frequency of the planetary gears at first. It is well known that these methods are capable of demodulating the signal when it is not contaminated seriously by noise [35, 36]. But when random background noise is added to the signal, the demodulation ability of these methods diminishes dramatically. When the signal is contaminated by noise with SNR= −6 dB, the EMD analysis results are shown in Figures 6(a) and 6(b). The fault characteristic frequency is not very clear in the spectrogram.Figure 6
The EMD analysis result of contaminated signal with SNR = −6 dB.
(a)
An intrinsic mode function (IMF) of EMD decomposition
(b)
The frequency spectrogram of the IMFWhen the signal is contaminated by noise with SNR = −8 dB, the envelope and the frequency spectrogram of Hilbert transform after filtering are shown in Figures7(a) and 7(b). The fault characteristic frequency cannot be seen in the envelope spectrum.Figure 7
The Hilbert analysis result of contaminated signal with SNR = −8 dB.
(a)
The envelope in time-domain
(b)
The envelope frequency spectrogramIn these methods, the EMD can be seen as a process of band filter, and the band-pass filter is applied prior to the Hilbert transform. The noise out of meshing frequency band is filtered out. But the in-band noise distorts the original vibration signal causing the invisibility of faulty characteristic frequency.When the signal is contaminated by noise with SNR = −10 dB, the wavelet envelopes can be shown in Figure8(a). There is so much noise pollution that the original signal cannot be recognized. The ME as indicator is shown in Figure 8(b). The central scale 20.3 according to the ME is the central frequency of 492.6 Hz. The results show that the indicator of ME is accurate and robust to the noise up to SNR = −10 dB. After manifold learning, the wavelet envelope and the power spectrum are shown in Figures 8(c) and 8(d). All the characteristic frequencies can be seen in the envelope manifold spectrogram. These data experiments fully demonstrate the superiority of the envelope manifold demodulation method.Figure 8
The envelope manifold result of contaminated signal with SNR = −10 dB.
(a)
The wavelet envelope on time-scale plane
(b)
The ME versus scale curve
(c)
The envelope manifold in time-domain
(d)
The envelope manifold frequency spectrogram


## 3.1. Simulation of Faulty Planet Gear

The planet gears mesh with ring gear and sun gear simultaneously. When there are faults on a planet gear, the meshing location of the faulty tooth and mating gear changes with the planet carrier rotation. Thus this changes the transfer path of the vibration signal of faulty tooth to the sensor and generates an extra AM effect on the meshing vibration signal. The meshing vibration amplitudes sampled by the sensor are modulated by the faulty tooth and the varied propagations synchronously [8]. Here, the faulty planet gear is analyzed. The faulty planet gear vibration signal can be written as(7)xt=1-Ccos⁡2πfcarriert1+Acos⁡2πfplanett+∅cos⁡2πfmesht+Bsin⁡2πfplanett+φ+θ,where xt is vibration signal in time-domain, fcarrier is planet carrier rotation frequency which represents propagation modulation, fplanet is planet gear rotation frequency which represents faulty planet gear modulation, fmesh is meshing frequency, and ∅, φ, and θ are the initial phases of AM, FM, and meshing, respectively. A, B, and C are the modulation amplitudes. After amplitude demodulation, the envelope of the meshing vibration signal is(8)yt=1-Ccos⁡2πfcarriert+Acos⁡2πfplanett+∅+12ACcos⁡2πfplanet+fcarriert+∅+12ACcos⁡2πfplanet-fcarriert+∅.So the envelope amplitude spectrum is(9)Yf=δf+Cδf-fcarrier+Aδf-fplanet+12ACδf-fplanet+fcarrier+12ACδf-fplanet-fcarrier.The simulation parameters of the planetary gearbox are shown in Table1. For comparison, the time-domain diagram and frequency spectrogram are adopted to analyze the vibration signals at first. Then the envelope manifold demodulation is utilized to illustrate the advantage of this method.Table 1
The simulation parameters of faulty planet gear.
Parameters
Values
AM magnitudeA
1
AM magnitudeC
0.8
FM magnitudeB
0.5
Meshing frequency (Hz)
437.5
Planet carrier rotation frequency (Hz)
5.468
Planet gear rotation frequency (Hz)
15.191The simulated vibration signal in time-domain is shown in Figure2(a). The frequency spectrogram of the signal is shown in Figure 2(b). The meshing frequency of 437.5 Hz and the sideband can be seen in the spectrogram. It is not very clear to identify the planet carrier rotation frequency and the planet gear rotation frequency.Figure 2
The time-domain diagram and frequency spectrogram of simulated signal.
(a)
The time-domain diagram
(b)
The frequency spectrogramIn envelope manifold demodulation, a series of scales of wavelet ranging from 1 to 30 are chosen to cover the meshing frequency of all the planetary gears. An increment of 0.1 for scale is employed in the data experiments. The wavelet envelopes are shown in Figure3. The lighter color in the graph illustrates greater wavelet envelopes. So the area circled by the ellipse corresponds to the greatest wavelet envelopes. The meshing frequency computed by formula (3) is the central frequency of the circled area.Figure 3
The wavelet envelope of simulated signal.The mean energy (ME), smoothness index (SI), and spectral kurtosis (SK) of wavelet envelope are utilized to choose the proper indictor of meshing frequency for comparison. The ME, SI, and SK as indicators are shown in Figures4(a), 4(b), and 4(c). The central scale, respectively, is 22.8, 22.2, and 29.0 according to the ME, SI, and SK. The corresponding frequency, respectively, is 438.5 Hz, 450.5 Hz, and 344.8 Hz. The result from ME is closest to the meshing frequency band which is circled by the ellipse in Figure 3. The indicator ME is most accurate for selection of meshing frequency band.Figure 4
The indicators versus scale curve of simulated signal.
(a)
The ME versus scale curve
(b)
The SI versus scale curve
(c)
The SK versus scale curveThe wavelet envelope at the central scale is shown in Figure5(a) and the power spectrum of the wavelet envelope is shown in Figure 5(b). The characteristic frequencies of the power spectrum are 4.89 Hz, 9.78 Hz, 14.66 Hz, and 19.55 Hz corresponding to fcarrier, fplanet-fcarrier, fplanet, and fplanet+fcarrier. The characteristic frequency of faulty planet gear can be seen clearly on the spectrogram of the wavelet envelope because there is no in-band noise interference. The wavelet envelope is proved to be adept at demodulation of planetary gear vibration signal.Figure 5
The wavelet envelope and spectrogram at the central scale of simulated signal.
(a)
The wavelet envelope
(b)
The envelope frequency spectrogram


## 3.2. Simulation of Contaminated Signal

In order to simulate the contamination of the in-band noise, the random background noise is added to the simulated signal of the faulty planet gear model to generate the contaminated signal up to SNR= −10 dB. Meanwhile, the meshing frequency is switched from 437.5 Hz to 492.0 Hz to testify the accuracy of selection of meshing frequency by ME.For comparison, the amplitude demodulation methods based on empirical mode decomposition (EMD) and Hilbert transform are used to reveal the characteristic frequency of the planetary gears at first. It is well known that these methods are capable of demodulating the signal when it is not contaminated seriously by noise [35, 36]. But when random background noise is added to the signal, the demodulation ability of these methods diminishes dramatically. When the signal is contaminated by noise with SNR= −6 dB, the EMD analysis results are shown in Figures 6(a) and 6(b). The fault characteristic frequency is not very clear in the spectrogram.Figure 6
The EMD analysis result of contaminated signal with SNR = −6 dB.
(a)
An intrinsic mode function (IMF) of EMD decomposition
(b)
The frequency spectrogram of the IMFWhen the signal is contaminated by noise with SNR = −8 dB, the envelope and the frequency spectrogram of Hilbert transform after filtering are shown in Figures7(a) and 7(b). The fault characteristic frequency cannot be seen in the envelope spectrum.Figure 7
The Hilbert analysis result of contaminated signal with SNR = −8 dB.
(a)
The envelope in time-domain
(b)
The envelope frequency spectrogramIn these methods, the EMD can be seen as a process of band filter, and the band-pass filter is applied prior to the Hilbert transform. The noise out of meshing frequency band is filtered out. But the in-band noise distorts the original vibration signal causing the invisibility of faulty characteristic frequency.When the signal is contaminated by noise with SNR = −10 dB, the wavelet envelopes can be shown in Figure8(a). There is so much noise pollution that the original signal cannot be recognized. The ME as indicator is shown in Figure 8(b). The central scale 20.3 according to the ME is the central frequency of 492.6 Hz. The results show that the indicator of ME is accurate and robust to the noise up to SNR = −10 dB. After manifold learning, the wavelet envelope and the power spectrum are shown in Figures 8(c) and 8(d). All the characteristic frequencies can be seen in the envelope manifold spectrogram. These data experiments fully demonstrate the superiority of the envelope manifold demodulation method.Figure 8
The envelope manifold result of contaminated signal with SNR = −10 dB.
(a)
The wavelet envelope on time-scale plane
(b)
The ME versus scale curve
(c)
The envelope manifold in time-domain
(d)
The envelope manifold frequency spectrogram


## 4. Experimental Signal Analysis

### 4.1. Experiment Design

In this section, the cases of practical experiments are studied to verify the effectiveness of the envelope manifold demodulation methodology. The planetary gearbox test rig is shown in Figure9. The planetary gear set in the gearbox is shown in Figure 10. The motor drives the input shaft of planetary gearbox through coupler. The input shaft is coupled to the sun gear, the ring gear is standstill, and the planet carrier is coupled to the output shaft which drives the other load gearboxes and load motor. The rotation speed of the planetary gearbox input shaft is measured by the tachometer. The vibration signals are collected by the accelerometer mounted on the planetary gearbox. The experimental planet gears are shown in Figure 11. The parameters of planetary gearbox are listed in Table 2.Table 2
The planetary gearbox parameters.
Parameters
Values
Planet gear teeth
36
Sun gear teeth
28
Ring gear teeth
100
Number of planet gears
3Figure 9
The planetary gearbox rig.Figure 10
The planetary gear set.Figure 11
The experimental planet gears.
(a)
(b)
(c)Three experiments are designed for comparison. The first case is the experiment of the normal planet gear, the second is the man-made tooth missing planet gear that simulates the local fault, and the third is the teeth surface worn planet gear that simulates the distributed fault.During experiments, the rotation speed of drive motor is set to 1500 rpm. The working parameters of the planetary gearbox in the experiments can be computed which are listed in Table3. The collected vibration signals in practical experiments are very complicated because of the abundant components in every frequency band. These practical vibration signals are analyzed with envelope manifold demodulation method to extract the fault features properly as follows.Table 3
The working parameters of experiments.
Parameters
Value (Hz)
Sun gear frequency (fs)
25
Planet gear meshing frequency (fm)
546.88
Planet carrier rotation frequency (fca)
5.47
Planet gear pass frequency (fpp)
21.88
Planet gear rotation frequency (fpr)
9.72
Faulty planet gear frequency (fpf)
15.19


### 4.2. Experiment Analysis

#### 4.2.1. Normal Planet Gear

In order to demonstrate the denoising effect of envelope manifold demodulation at first, the wavelet envelope at the central scale and the envelope manifold are studied for comparison. The central scale is 18.2 according to ME indicator. The corresponding meshing frequency, respectively, is 549.5 Hz. But the central scale, respectively, is 8.0 or 1.8 according to SI or SK. The corresponding frequency, respectively, is 1250 Hz or 5555.6 Hz. So the ME indicator is most suitable for identifying the planet gear meshing frequency in the practical experiments. The wavelet envelope at the central scale is shown in Figure12(a) and the wavelet envelope frequency spectrogram is shown in Figure 12(b). After manifold learning, the envelope manifold and the frequency spectrogram are shown in Figure 13. The in-band noise causes the distortion of wavelet envelope. The characteristic frequency of envelope manifold is clearer than the wavelet envelope by comparing the frequency spectrograms in these figures.Figure 12
The wavelet envelope at the central scale of normal planet gear.
(a)
The wavelet envelope in time-domain
(b)
The wavelet envelope frequency spectrogramFigure 13
The envelope manifold of normal planet gear.
(a)
The envelope manifold in time-domain
(b)
The envelope manifold frequency spectrogramThe meshing vibration signal of planetary gearbox is modulated by different meshing phases, varied propagation paths, and even the manufacture error of the gear pairs. So in the envelope manifold spectrogram of the normal planet gear shown in Figure13(b), the characteristic frequencies of fpp, fpr, fpp+fca, fpp+2fca, 2fpp, 2fpp+fca, 2fpp+2fca,… can be recognized. The characteristic frequencies are related to fpp, fpr, and fca but not associated with the faulty planet gear frequency. From spectrogram, we can find the characteristic frequencies caused by the manufacture error, that is, the unbalance of the planet carrier and multiple planet gears. Their amplitudes are not big. So there is no fault on the planet gears.


#### 4.2.2. Tooth Missing Planet Gear

The envelope manifold of the tooth missing planet gear is shown in Figure14. The characteristic frequencies of fpf, fpr, fpp+fpf, fpp+fca, fpp+2fca,… can be seen in the frequency spectrogram. The modulating frequency of faulty planet gear fpf can be extracted from the demodulated spectrum. fpf, fpp, and fpr jointly modulated the meshing frequency of the planet gear. fpf and fpp+fpf become the mainly modulating frequencies besides those characteristic frequencies of normal planet gear. The faulty planet gear frequency and the harmonic frequencies can be distinguished in the frequency spectrogram.Figure 14
The envelope manifold of tooth missing planet gear.
(a)
The envelope manifold in time-domain
(b)
The envelope manifold frequency spectrogram


#### 4.2.3. Teeth Surface Worn Planet Gear

The envelope manifold of the teeth surface worn planet gear is shown in Figure15. The characteristic frequencies of fca, fpp+fpf, fpp+2fca, fpf+fca, fpf+fpr,… can be viewed in the frequency spectrogram. Because there is just one teeth surface worn planet gear in the gearbox, the unbalance of the planet carrier increases dramatically. The effect of planet carrier characteristic frequency is enhanced. When the teeth surface worn planet gear meshes with mating gears, it can cause the faulty planet gear characteristic frequency. And so the frequencies of fpf+fpr and fpf+fca appear in the power spectrogram. The teeth surface worn planet gear can be identified by this way.Figure 15
The envelope manifold of teeth surface worn planet gear.
(a)
The envelope manifold in time-domain
(b)
The envelope manifold frequency spectrogramThe planetary gearbox rig sets and experiment conditions are different from each other. All the experiment results are based on the special rig set and experiment condition, but the experiments validate the effectiveness of the envelope manifold demodulation methodology.


## 4.1. Experiment Design

In this section, the cases of practical experiments are studied to verify the effectiveness of the envelope manifold demodulation methodology. The planetary gearbox test rig is shown in Figure9. The planetary gear set in the gearbox is shown in Figure 10. The motor drives the input shaft of planetary gearbox through coupler. The input shaft is coupled to the sun gear, the ring gear is standstill, and the planet carrier is coupled to the output shaft which drives the other load gearboxes and load motor. The rotation speed of the planetary gearbox input shaft is measured by the tachometer. The vibration signals are collected by the accelerometer mounted on the planetary gearbox. The experimental planet gears are shown in Figure 11. The parameters of planetary gearbox are listed in Table 2.Table 2
The planetary gearbox parameters.
Parameters
Values
Planet gear teeth
36
Sun gear teeth
28
Ring gear teeth
100
Number of planet gears
3Figure 9
The planetary gearbox rig.Figure 10
The planetary gear set.Figure 11
The experimental planet gears.
(a)
(b)
(c)Three experiments are designed for comparison. The first case is the experiment of the normal planet gear, the second is the man-made tooth missing planet gear that simulates the local fault, and the third is the teeth surface worn planet gear that simulates the distributed fault.During experiments, the rotation speed of drive motor is set to 1500 rpm. The working parameters of the planetary gearbox in the experiments can be computed which are listed in Table3. The collected vibration signals in practical experiments are very complicated because of the abundant components in every frequency band. These practical vibration signals are analyzed with envelope manifold demodulation method to extract the fault features properly as follows.Table 3
The working parameters of experiments.
Parameters
Value (Hz)
Sun gear frequency (fs)
25
Planet gear meshing frequency (fm)
546.88
Planet carrier rotation frequency (fca)
5.47
Planet gear pass frequency (fpp)
21.88
Planet gear rotation frequency (fpr)
9.72
Faulty planet gear frequency (fpf)
15.19


## 4.2. Experiment Analysis

### 4.2.1. Normal Planet Gear

In order to demonstrate the denoising effect of envelope manifold demodulation at first, the wavelet envelope at the central scale and the envelope manifold are studied for comparison. The central scale is 18.2 according to ME indicator. The corresponding meshing frequency, respectively, is 549.5 Hz. But the central scale, respectively, is 8.0 or 1.8 according to SI or SK. The corresponding frequency, respectively, is 1250 Hz or 5555.6 Hz. So the ME indicator is most suitable for identifying the planet gear meshing frequency in the practical experiments. The wavelet envelope at the central scale is shown in Figure12(a) and the wavelet envelope frequency spectrogram is shown in Figure 12(b). After manifold learning, the envelope manifold and the frequency spectrogram are shown in Figure 13. The in-band noise causes the distortion of wavelet envelope. The characteristic frequency of envelope manifold is clearer than the wavelet envelope by comparing the frequency spectrograms in these figures.Figure 12
The wavelet envelope at the central scale of normal planet gear.
(a)
The wavelet envelope in time-domain
(b)
The wavelet envelope frequency spectrogramFigure 13
The envelope manifold of normal planet gear.
(a)
The envelope manifold in time-domain
(b)
The envelope manifold frequency spectrogramThe meshing vibration signal of planetary gearbox is modulated by different meshing phases, varied propagation paths, and even the manufacture error of the gear pairs. So in the envelope manifold spectrogram of the normal planet gear shown in Figure13(b), the characteristic frequencies of fpp, fpr, fpp+fca, fpp+2fca, 2fpp, 2fpp+fca, 2fpp+2fca,… can be recognized. The characteristic frequencies are related to fpp, fpr, and fca but not associated with the faulty planet gear frequency. From spectrogram, we can find the characteristic frequencies caused by the manufacture error, that is, the unbalance of the planet carrier and multiple planet gears. Their amplitudes are not big. So there is no fault on the planet gears.


### 4.2.2. Tooth Missing Planet Gear

The envelope manifold of the tooth missing planet gear is shown in Figure14. The characteristic frequencies of fpf, fpr, fpp+fpf, fpp+fca, fpp+2fca,… can be seen in the frequency spectrogram. The modulating frequency of faulty planet gear fpf can be extracted from the demodulated spectrum. fpf, fpp, and fpr jointly modulated the meshing frequency of the planet gear. fpf and fpp+fpf become the mainly modulating frequencies besides those characteristic frequencies of normal planet gear. The faulty planet gear frequency and the harmonic frequencies can be distinguished in the frequency spectrogram.Figure 14
The envelope manifold of tooth missing planet gear.
(a)
The envelope manifold in time-domain
(b)
The envelope manifold frequency spectrogram


### 4.2.3. Teeth Surface Worn Planet Gear

The envelope manifold of the teeth surface worn planet gear is shown in Figure15. The characteristic frequencies of fca, fpp+fpf, fpp+2fca, fpf+fca, fpf+fpr,… can be viewed in the frequency spectrogram. Because there is just one teeth surface worn planet gear in the gearbox, the unbalance of the planet carrier increases dramatically. The effect of planet carrier characteristic frequency is enhanced. When the teeth surface worn planet gear meshes with mating gears, it can cause the faulty planet gear characteristic frequency. And so the frequencies of fpf+fpr and fpf+fca appear in the power spectrogram. The teeth surface worn planet gear can be identified by this way.Figure 15
The envelope manifold of teeth surface worn planet gear.
(a)
The envelope manifold in time-domain
(b)
The envelope manifold frequency spectrogramThe planetary gearbox rig sets and experiment conditions are different from each other. All the experiment results are based on the special rig set and experiment condition, but the experiments validate the effectiveness of the envelope manifold demodulation methodology.


## 4.2.1. Normal Planet Gear

In order to demonstrate the denoising effect of envelope manifold demodulation at first, the wavelet envelope at the central scale and the envelope manifold are studied for comparison. The central scale is 18.2 according to ME indicator. The corresponding meshing frequency, respectively, is 549.5 Hz. But the central scale, respectively, is 8.0 or 1.8 according to SI or SK. The corresponding frequency, respectively, is 1250 Hz or 5555.6 Hz. So the ME indicator is most suitable for identifying the planet gear meshing frequency in the practical experiments. The wavelet envelope at the central scale is shown in Figure12(a) and the wavelet envelope frequency spectrogram is shown in Figure 12(b). After manifold learning, the envelope manifold and the frequency spectrogram are shown in Figure 13. The in-band noise causes the distortion of wavelet envelope. The characteristic frequency of envelope manifold is clearer than the wavelet envelope by comparing the frequency spectrograms in these figures.Figure 12
The wavelet envelope at the central scale of normal planet gear.
(a)
The wavelet envelope in time-domain
(b)
The wavelet envelope frequency spectrogramFigure 13
The envelope manifold of normal planet gear.
(a)
The envelope manifold in time-domain
(b)
The envelope manifold frequency spectrogramThe meshing vibration signal of planetary gearbox is modulated by different meshing phases, varied propagation paths, and even the manufacture error of the gear pairs. So in the envelope manifold spectrogram of the normal planet gear shown in Figure13(b), the characteristic frequencies of fpp, fpr, fpp+fca, fpp+2fca, 2fpp, 2fpp+fca, 2fpp+2fca,… can be recognized. The characteristic frequencies are related to fpp, fpr, and fca but not associated with the faulty planet gear frequency. From spectrogram, we can find the characteristic frequencies caused by the manufacture error, that is, the unbalance of the planet carrier and multiple planet gears. Their amplitudes are not big. So there is no fault on the planet gears.


## 4.2.2. Tooth Missing Planet Gear

The envelope manifold of the tooth missing planet gear is shown in Figure14. The characteristic frequencies of fpf, fpr, fpp+fpf, fpp+fca, fpp+2fca,… can be seen in the frequency spectrogram. The modulating frequency of faulty planet gear fpf can be extracted from the demodulated spectrum. fpf, fpp, and fpr jointly modulated the meshing frequency of the planet gear. fpf and fpp+fpf become the mainly modulating frequencies besides those characteristic frequencies of normal planet gear. The faulty planet gear frequency and the harmonic frequencies can be distinguished in the frequency spectrogram.Figure 14
The envelope manifold of tooth missing planet gear.
(a)
The envelope manifold in time-domain
(b)
The envelope manifold frequency spectrogram


## 4.2.3. Teeth Surface Worn Planet Gear

The envelope manifold of the teeth surface worn planet gear is shown in Figure15. The characteristic frequencies of fca, fpp+fpf, fpp+2fca, fpf+fca, fpf+fpr,… can be viewed in the frequency spectrogram. Because there is just one teeth surface worn planet gear in the gearbox, the unbalance of the planet carrier increases dramatically. The effect of planet carrier characteristic frequency is enhanced. When the teeth surface worn planet gear meshes with mating gears, it can cause the faulty planet gear characteristic frequency. And so the frequencies of fpf+fpr and fpf+fca appear in the power spectrogram. The teeth surface worn planet gear can be identified by this way.Figure 15
The envelope manifold of teeth surface worn planet gear.
(a)
The envelope manifold in time-domain
(b)
The envelope manifold frequency spectrogramThe planetary gearbox rig sets and experiment conditions are different from each other. All the experiment results are based on the special rig set and experiment condition, but the experiments validate the effectiveness of the envelope manifold demodulation methodology.


## 5. Conclusions

In this paper, the envelope manifold demodulation methodology was introduced to implement the fault diagnosis of planetary gearbox. The method consisted of the complex wavelet transform, manifold learning, and frequency spectrogram analysis method. The continuous wavelet transform with complex Morlet wavelet was utilized to extract the envelope for demodulation of the planetary gears meshing vibration signal at first. The scale band selection was actualized through mean envelope energy (ME) indicator to locate the meshing frequency and remove the noisy signal not within the meshing frequency band. Then, the manifold learning method is applied to extract the intrinsic manifold of the envelope and reduce the effect of in-band noise. At last, the frequency spectrogram is used to show the characteristic frequency of the faulty planetary gears meshing envelope. By this method, the effect of noise among the vibration signals was immensely weakened, the envelope manifold structure was clearly revealed, and the fault characteristics of planetary gearbox were conveniently extracted. The result of simulated signal and experimental dataset all verified the effectiveness of the method. This method revealed the intrinsic characteristics of the planetary gearbox vibration signal. It was observed that this method outperformed other approaches in various conditions. The further work of us will be research on the exact mechanism of the planetary gearbox vibrations based on the results of this paper.



---

*Source: 3952325-2015-12-28.xml*