Next Article in Journal
Parameter Estimation of Lunar Regolith from Lunar Penetrating Radar Data
Next Article in Special Issue
A Novel Fault Detection with Minimizing the Noise-Signal Ratio Using Reinforcement Learning
Previous Article in Journal
A Sound Source Localisation Analytical Method for Monitoring the Abnormal Night Vocalisations of Poultry
Previous Article in Special Issue
A Fiber Bragg Grating Based Torsional Vibration Sensor for Rotating Machinery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Early Fault Diagnosis for Planetary Gearbox Based Wavelet Packet Energy and Modulation Signal Bispectrum Analysis

1
School of Mechanical Engineering, Hebei University of Technology, Tianjin 300401, China
2
Centre for Efficiency and Performance Engineering, University of Huddersfield, Huddersfield HD1 3DH, UK
*
Author to whom correspondence should be addressed.
Submission received: 17 August 2018 / Revised: 28 August 2018 / Accepted: 29 August 2018 / Published: 1 September 2018
(This article belongs to the Special Issue Sensors for Fault Diagnosis and Fault Tolerance)

Abstract

:
The planetary gearbox is at the heart of most rotating machinery. The premature failure and subsequent downtime of a planetary gearbox not only seriously affects the reliability and safety of the entire rotating machinery but also results in severe accidents and economic losses in industrial applications. It is an important and challenging task to accurately detect failures in a planetary gearbox at an early stage to ensure the safety and reliability of the mechanical transmission system. In this paper, a novel method based on wavelet packet energy (WPE) and modulation signal bispectrum (MSB) analysis is proposed for planetary gearbox early fault diagnostics. First, the vibration signal is decomposed into different time-frequency subspaces using wavelet packet decomposition (WPD). The WPE is calculated in each time-frequency subspace. Secondly, the relatively high energy vectors are selected from a WPE matrix to obtain a reconstructed signal. The reconstructed signal is then subjected to MSB analysis to obtain the fault characteristic frequency for fault diagnosis of the planetary gearbox. The validity of the proposed method is carried out through analyzing the vibration signals of the test planetary gearbox in two fault cases. One fault is a chipped sun gear tooth and the other is an inner-race fault in the planet gear bearing. The results show that the proposed method is feasible and effective for early fault diagnosis in planetary gearboxes.

1. Introduction

Planetary gearboxes are widely used in wind turbines and other large machines and are an important feature of mechanical transmission systems. As planetary gearboxes usually work under heavy loading and severe environment, their failure may cause accidents and result in higher maintenance costs [1,2]. Early diagnosis of planetary gearbox contributes to the normal operation of mechanical equipment [3,4]. One of the most widely used methods to diagnose planetary gearbox failure is vibration-based analysis, because it is easily measured and contains plenty of fault information [5]. By applying advanced signal processing methods, one can extract the features of the faulty components from the measured vibration signals [6,7]. Therefore, it is important to investigate an effective and feasible method to accurately extract fault features from vibration signals to detect planetary gearbox failures.
At present, fault diagnosis methods could mainly be divided into two categories as model-based and data-driven [8]. The model-based method requires extensive system knowledge and related dynamics, which is difficult to implement in practically complex applications [9]. By contrast, the data-driven method has been widely used in rotating machinery fault diagnosis because detailed physical system dynamic knowledge is not required [10,11,12]. Moreover, many regular analysis methods for fault feature extraction have been applied to planetary gearbox fault diagnosis [13,14,15], such as envelope analysis (EA), maximum correlated kurtosis deconvolution (MCKD), Wigner–Viller distribution (WVD), empirical mode decomposition (EMD), etc. Despite the above method having been widely applied to the fault diagnosis of planetary gearboxes, all methods have limitations. For example, EA needs to choose center frequency and bandwidth for the bandpass filter in advance [13]; MCKD is restricted by its parameter selection, which will affect the denoising result [14]; WVD is subjected to cross term interference [15]; and EMD is often affected by boundary effects and modes mixing [16,17].
Different from these regular analysis methods Wavelet Packet Energy (WPE) is an alternative method formed by linear combinations of usual wavelet functions [18]. WPE can be used to reveal the energy distribution information of the vibration signal in each time-frequency subspace and highlight the energy concentrated subspaces. According to the failure mechanism of bearing, the energy level of the vibration signals related to the fault will be increased and concentrated in a certain frequency band when a fault occurs [19]. Therefore, the frequency band with high energy level in the divided time-frequency subspaces which are decomposed by WPE is more related to the fault and can be used to indicate the fault features. WPE has received much attention in the area of machinery fault diagnosis [20,21,22,23]. He proposed a novel method to probe the WPE flow characteristics of vibration signals by using the manifold learning technique [20]. Wang et al. proposed a novel method based on WPE and manifold learning to extract the weak signature and outperform the Kurtogram-based methods [21]. Gómez et al. calculated WPE as an effective fault feature to detect railway axle cracks, then classify and identify fault features through artificial neural network (ANN) [22]. Zhou et al. proposed a hybrid method that the WPE value of the first several intrinsic mode functions (IMFs)IMFs decomposed by EMD is taken as the fault feature, and then classified and identified by support vector machine (SVM)SVM [23]. The results of the above studies have indicated that WPE could effectively diagnose machine failures. However, the reconstructed signal from the relatively high energy vector of the WPE matrix still has noise and modulation components.
Recently, Modulation Signal Bispectrum (MSB) has appeared in fault diagnosis thanks to its ability to effectively reveal modulation characteristics and suppress noise. Gu et al. explored a novel analysis method for motor broken rotor bar diagnosis using the MSB-based motor current signal analysis (MCSA) [24]. Zhang attested MSB-MCSA is an effective and precise method to monitor gear wear [25]. Tian demonstrated that an MSB-based method can produce more exact and robust test results for the fault diagnosis of bearings in some scenarios [26]. However, the frequency implementation of MSB could limit MSB noise suppression performance because of spectral smearing.
In this paper, a new analysis method based on WPE and MSB for planetary gearbox fault diagnosis is proposed. First, the signal energy will be decomposed into different time-frequency subspaces by WPT, forming a WPE flowing from the low layer to the high layer. Then, the energy concentrated time-frequency subspaces are selected from WPE, and an inverse WPT is applied to the selected bands to obtain a reconstructed signal. Finally, MSB is used to suppress the noise and demodulate the modulation components from the reconstructed signal, and then extract the fault characteristic frequency for the fault diagnosis of planetary gearboxes. The results indicate that the proposed approach performs excellently in the early extraction of planetary gearbox fault characteristics.
The rest of the paper will be organized as follows. Section 2 briefly introduces the theoretical background on the WPE feature extraction. Section 3 describes the fundamental theory of MSB and its application on sideband estimation. Section 4 presents the process of the proposed method based on WPE and MSB for planetary gearbox fault diagnosis. The proposed method is verified by two planetary gearbox application cases in Section 5. The conclusions are given in Section 6.

2. Wavelet Packet Transform and Energy Feature Extraction

Let V 0 , 0 is a vector space of R n , which is the node 0 of in the parent tree. The vector space could be split into two mutually orthogonal subspaces at each level and be expressed as [27]:
V j , k = V j + 1 , 2 k U j + 1 , 2 k + 1
where j and k are the layer of the tree and index of the node in layer j respectively, and k = 0 , , 2 j 1 .
The decomposition process of WPT is depicted as follows. First, the WP function W j , k n ( t ) can be expressed as [28]:
W j , k n ( t ) = 2 j / 2 W n ( 2 j t k )
where j indicates the scale, k represents the translation parameters, and n = 0 , 1 , . When j = k = 0 , the first two wavelet packet functions Φ ( t ) and Ψ ( t ) are the scaling function and mother wavelet function respectively, and they can be expressed as follows:
W 0 , 0 0 ( t ) = Φ ( t )
W 0 , 0 1 ( t ) = Ψ ( t )
The remaining WP functions of n = 2 , 3 , can be defined by the following recursive relationships:
W 2 n ( t ) = 2 k h ( k ) W 1 , k n ( 2 t k )
W 2 n + 1 ( t ) = 2 k g ( k ) W 1 , k n ( 2 t k )
where h ( k ) = 1 / 2 φ ( t ) , φ ( 2 t k ) is the filter coefficient of low-pass filter and g ( k ) = 1 / 2 Ψ ( t ) ,   Ψ ( 2 t k ) is the filter coefficient of high-pass filter. They are orthogonal to the relationship g ( k ) = ( 1 ) k h ( 1 k ) . Here, . , . indicates the inner product operator. The wavelet packet coefficients Q j , k n can be calculated through the signal x ( t ) and the wavelet packet functions W j , k n , as [29]:
Q j , k   n = x , W j , k n ( t ) = + x ( t ) W j , k n ( t ) d t
where Q j . k n and k are the nth set of WP coefficients at the jth scale parameter and translation parameter, respectively.
As a result, at the jth layer, the signal x ( t ) is separated into 2 j packets with the order n = 1 , 2 , , 2 j . The nodes can be indexed by ( j , n ) in the binary tree structure, and the tree structure of the WPT comprises nodes with equal bandwidth for each scale parameter j. After reconstructing the coefficients at each node layer by layer in the tree structure, a waveform feature space is obtained.
For a discrete signal x ( t ) with N ( N = 2 n 0 ) points, the reconstructed signal of wavelet coefficients in node can be expressed using { R j n ( k ) , k = 1 , 2 , , 2 n 0 } . If assuming the WP node ( j , n ) is a container in the area of 2 n 0 , the waveform featured distributed in the container can be represented by R j n . Then the waveform feature (WF) space for the ( j , n ) th layer can be expressed as follows [21]:
W F j = { x ( t )   , j = 0 [ [ R j 1 ] 2 n 0 × 1 [ R j n ] 2 n 0 × 1 [ R j 2 j ] 2 n 0 × 1 n = 1 , 2 , , 2 j ] R 2 n 0 × 2 j , j ( 1 , J )
Put W F j from layer 0 to J , then the WF is expressed as follows:
WF = [ [ W F 0 ] 2 n 0 × 2 0 [ W F n ] 2 n 0 × 2 n [ W F 0 ] 2 n 0 × 2 J n = 1 , 2 , , 2 j ] T R D × 2 n 0
where D = 2 0 + 2 1 + + 2 J , which is the dimension of the W F space.

3. Modulation Signal Bispectrum

For a discrete time series x ( t ) with corresponding discrete fourier transform X ( f ) , the modulation signal bispectrum (MSB) can be expressed in the frequency domain as [24]:
B M S ( f 1 , f 2 ) = E X ( f 1 + f 2 ) X ( f 1 f 2 ) X * ( f 1 ) X * ( f 1 )
where B M S ( f 1 , f 2 ) is the bispectrum of signal x ( t ) , E < > indicates the expectation operator, f 1 is carrier frequency, f 2 is the modulating frequency, ( f 1 + f 2 ) and ( f 1 f 2 ) are the higher and lower sideband frequencies. The bispectral peak of B M S ( f 1 , f 2 ) can represent the modulation characteristics.
To quantify the sideband amplitude more accurately, the MSB is modified via removing the effect of f 1 through magnitude normalization. The MSB sideband estimator is defined as [26]:
B M S S E ( f 1 , f 2 ) = B M S ( f 1 , f 2 ) | B M S ( f 1 , 0 ) |
where B M S ( f 1 , 0 ) represents the squared power spectrum estimation at f x = 0 .
In order to obtain more robust results, the estimator is improved based on the average of several suboptimal MSB slices, and can be denoted by:
B ( f 2 ) = 1 N n = 1 N B M S S E ( f 1 n , f 2 ) ( f 2 > 0 )
where N indicates the whole number of selected f 1 suboptimal slices, the number of which relies on the significance of the peaks themselves.
To get suboptimal f 1 slices, the B M S S E ( f 1 n , f 2 ) is defined as the compound MSB slice B ( f 1 ) , calculated by averaging the main MSB peaks in the incremental direction of the f 2 :
B ( f 1   ) = 1 M 1 m = 2 M B M S S E ( f 1 , m Δ f )
where Δ f denotes the frequency resolution in the f 2 direction.

4. The Proposed Method Based on WPT and MSB

Based on the advantages of WPE and MSB, the WPE–MSB for fault diagnosis of planetary gearbox is proposed. The scheme of the proposed method is briefly illustrated in Figure 1 and depicted as follows:
Step 1:
Carry out WPT to process the vibration signal in order to separate it into different time-frequency subspaces.
Step 2:
Calculate the WPE matrix in time-frequency subspaces.
Step 3:
Select relatively high energy vectors from the WPE matrix.
Step 4:
Apply an inverse WPT to the selected bands to obtain reconstructed signal.
Step 5:
Calculate MSB of the reconstructed signal to extract the fault characteristic frequency.

5. Experimental Verification

5.1. Experimental Setup

The vibration signal was recorded based on the planetary gearbox test rig shown in Figure 2. The rated torque of the test planetary gearbox was 670   Nm , and the maximum input speed was 2800 rpm with a resulting output speed of 388 rpm . The accelerometer with a sensitivity of 28.7 mV / ms 2 was mounted on the outer housing of the ring gear along the position of the test planetary gear bearing. The collected data were sampled at 96 kHz with a resolution of 24 bits in the measurement system as shown in Figure 3. The fault modes include sun gear chipped and planetary gear bearing with inner-race fault, as illustrated in Figure 4a,b, respectively. The structural and kinematical parameters of the test planetary gearbox and planetary gear bearing are listed in Table 1 and Table 2.

5.2. The Fault Diagnosis of the Sun Gear Chipped Tooth

The measured vibration signal under the sun gear chipped tooth fault and the corresponding spectrum are shown in Figure 5a,b. It is difficult to find the characteristic frequency of the sun gear chipped due to the heavy noise.
The proposed method is used to analyze the measured vibration signal in Figure 4a. First, the WPT method is applied to extract fault features from the separated signals. The db 10 wavelet function is used for 4-layers decomposition, then 16 decomposition bands are obtained [30,31] and the frequency range of the 4-layers wavelet packet are shown in Table 3. Then the WPE matrix is calculated, and the energy distribution proportional histogram of the vibration signal of the sun gear chipped is shown in Figure 6. It is clear that the ratios of the energy distribution are higher at the node [4,2] and [4,3] than other nodes. Therefore, the above two bands are selected to reconstruct a signal. Envelope analysis is then applied to process the reconstructed signal, and the obtained normalized envelope spectrum is show in Figure 7. The results are thus not persuasive enough. Although some fault feature frequencies can be recognized, there are still many noise and modulation components.
In order to solve the above problem, MSB is then applied to suppress the noise and decompose the modulated components. The normalized analysis results of vibration signal for the sun gear chipped based on the proposed method are shown in Figure 8: Some peaks correspond exactly to sun gear rotation frequency f r s , sun gear characteristic frequency f s f , and the combination f s f ± f r s . The result demonstrates that the proposed method can get more effective and accurate fault features than the envelope analysis for sun gear fault diagnosis.

5.3. The Fault Diagnosis of the Planetary Gear Bearing

Figure 9a,b shows the waveform and spectrum of the measured vibration signal under the condition of planetary gear bearing with inner-race fault. It is difficult to identify the fault characteristic frequencies from the spectrum in Figure 9b. On the basis of the proposed method, the WPT is applied to decompose the vibration signal into a group of WP nodes with a complete binary tree form. Then the WPE matrix is calculated. The energy distribution proportional histogram of the vibration signal of the bearing inner-race fault is shown in Figure 10.
According to the WPE analysis results in Figure 10, ratios of the energy distribution at the node [4,0], [4,1], [4,2] and [4,3] are selected as the eigenvalues for the bearing inner-race fault signal. Then, the above four corresponding frequency bands are used to generate the reconstructed signal. Envelope analysis is then applied to analyze the fault signal; the corresponding normalized envelope spectrum is shown in Figure 11. It is obvious that the envelope analysis result does not present the characteristic frequencies of the bearing inner-race fault effectively owing to the strong noise interference. To deal with the problem, we took advantage of MSB to suppress noise and decompose modulations. The normalized results of the propose method is shown in Figure 12. It can be seen from the analysis results that the inner-race fault characteristic frequency f i and its harmonics can be clearly identified. In conclusion, the proposed method obtains more accurate results than the envelope analysis for the planetary gear bearing fault diagnosis.

6. Conclusions

This paper proposes a new analysis method based on WPE and MSB for planetary gearbox fault diagnosis. WPE reveals the energy distribution of the vibration signal in different time-frequency subspaces for the reconstructed signal obtained. MSB is used to enhance the reconstructed signal with relatively high energy vectors and highlight the fault features due to the effectiveness of MSB in demodulating modulation components and removing noise. Its efficiency has been evaluated on the experimental signals measured from planetary gearbox with sun gear chipped and planetary gear bearing with inner race fault. The results show that the proposed method is an effective and reliable method for early planetary gearbox faults diagnosis. However, the layers selection for the WPT and the spectral smearing of MSB would affect the diagnosis results of the proposed method. Hence, effective methods need to be studied to deal with the above issues in future works.

Author Contributions

All authors contributed to this work. J.G. analyzed the data and drafted the manuscript; F.G. and H.L. set up the experimental systems; D.Z. and Z.S. conceived and designed the experiments; D.Z., F.G. and A.D.B. contributed to the algorithm design and the results discussion. All the authors contributed to the manuscript draft and discussion.

Funding

This research is supported by the National Natural Science Foundation of China (Grant No. 51605133; 51705127), Hebei Provincial International Science and Technology Cooperation Program of China (Grant No. 17394303D).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Li, Z.; Jiang, Y.; Hu, C.; Peng, Z. Recent progress on decoupling diagnosis of hybrid failures in gear transmission systems using vibration sensor signal: A review. Measurement 2016, 90, 4–19. [Google Scholar] [CrossRef]
  2. Glowacz, A.; Glowacz, W.; Glowacz, Z.; Kozik, J. Early fault diagnosis of bearing and stator faults of the single-phase induction motor using acoustic signals. Measurement 2018, 113, 1–9. [Google Scholar] [CrossRef]
  3. Chen, X.; Feng, F.; Zhang, B. Weak Fault Feature Extraction of Rolling Bearings Based on an Improved Kurtogram. Sensors 2016, 16, 1482. [Google Scholar] [CrossRef] [PubMed]
  4. Saucedo-Dorantes, J.J.; Delgado-Prieto, M.; Ortega-Redondo, J.A.; Osornio-Rios, R.A.; Romero-Troncoso, R.D. Multiple-fault detection methodology based on vibration and current analysis applied to bearings in induction motors and gearboxes on the kinematic chain. Shock Vib. 2016, 2016, 5467643. [Google Scholar] [CrossRef]
  5. Lee, D.H.; Ahn, J.H.; Koh, B.H. Fault detection of bearing systems through EEMD and optimization algorithm. Sensors 2017, 17, 2477. [Google Scholar] [CrossRef] [PubMed]
  6. Ge, M.; Wang, J.; Zhang, F.; Bai, K.; Ren, X. A novel fault diagnosis method of rolling bearings based on AFEWT-KDEMI. Entropy 2018, 20, 445. [Google Scholar] [CrossRef]
  7. Cui, L.; Zhang, Y.; Zhang, F.; Zhang, J.; Lee, S. Vibration response mechanism of faulty outer race rolling element bearings for quantitative analysis. J. Sound Vib. 2016, 364, 67–76. [Google Scholar] [CrossRef]
  8. Benmoussa, S.; Djeziri, M.A. Remaining useful life estimation without needing for prior knowledge of the degradation features. Iet. Sci. Meas. Technol. 2017, 11, 1071–1078. [Google Scholar] [CrossRef]
  9. Liang, X.; Zuo, M.J.; Feng, Z. Dynamic modeling of gearbox faults: A review. Mech. Syst. Signal Process. 2018, 98, 852–876. [Google Scholar] [CrossRef]
  10. Pacheco, F.; Cerrada, M.; Cabrera, D.; Li, C.; Zuria, G. A statistical comparison of neuroclassifiers and feature selection methods for gearbox fault diagnosis under realistic conditions. Neurocomputing 2016, 194, 192–206. [Google Scholar] [CrossRef]
  11. Djeziri, M.A.; Ananou, B.; Ouladsine, M. Data driven and model based fault prognosis applied to a mechatronic system. In Proceedings of the 4th International Conference on Power Engineering, Energy and Electrical Drives, Istanbul, Turkey, 13–17 May 2013. [Google Scholar]
  12. Xiao, B.; Huo, M.; Yang, X.; Zhang, Y. Fault-tolerant attitude stabilization for satellites without rate sensor. IEEE. Technol. Ind. Electron. 2015, 62, 7191–7202. [Google Scholar] [CrossRef]
  13. Sui, W.; Osman, S.; Wang, W. An adaptive envelope spectrum technique for bearing fault detection. Meas. Sci. Technol. 2014, 25, 095004. [Google Scholar] [CrossRef]
  14. Feng, J.; Lei, Y.; Shan, H.; Lin, J. Early fault diagnosis of bearings using an improved spectral kurtosis by maximum correlated kurtosis deconvolution. Sensors 2015, 15, 29363–29377. [Google Scholar]
  15. Pachori, R.B.; Nishad, A. Cross-terms reduction in the Wigner–Ville distribution using tunable-Q wavelet transform. Signal Process. 2016, 120, 288–304. [Google Scholar] [CrossRef]
  16. Zhang, Y.G.; Zhang, C.H.; Sun, J.B.; Guo, J.J. Improved wind speed prediction using empirical mode decomposition. Adv. Electr. Comput. Eng. 2018, 18, 3–10. [Google Scholar] [CrossRef]
  17. Zhang, C.; Li, Z.; Chen, S.; Wang, J.; Zhang, X. Optimised ensemble empirical mode decomposition with optimised noise parameters and its application to rolling element bearing fault diagnosis. Insight 2016, 58, 494–501. [Google Scholar] [CrossRef]
  18. Li, D.; Liang, Y.; Feng, Q.; Song, G. Load monitoring of the pin-connected structure based on wavelet packet analysis using piezoceramic transducers. Measurement 2017, 122, 638–647. [Google Scholar] [CrossRef]
  19. Chen, B.; Wang, X.; Yan, W.; Tian, H.; Xiao, W.; Chen, F.; Liu, H. A RSSD fault diagnosis method for rolling bearings based on optimization of quality factors and reconstruction of sub-bands. J. Xian Jiaotong Univ. 2018, 52, 70–76. [Google Scholar]
  20. He, Q. Vibration signal classification by wavelet packet energy flow manifold learning. J. Sound Vib. 2013, 332, 1881–1894. [Google Scholar] [CrossRef]
  21. Wang, Y.; Xu, G.; Lin, L.; Jiang, K. Detection of weak transient signals based on wavelet packet transform and manifold learning for rolling element bearing fault diagnosis. Mech. Syst. Signal Process. 2015, 54–55, 259–276. [Google Scholar] [CrossRef]
  22. Gómez, M.; Corral, E.; Castejón, C.; Garciaprada, J. Effective crack detection in railway axles using vibration signals and WPT energy. Sensors 2018, 18, 1603. [Google Scholar] [CrossRef] [PubMed]
  23. Zhou, S.; Qian, S.; Chang, W.; Xiao, Y.; Cheng, Y. A novel bearing multi-fault diagnosis approach based on weighted permutation entropy and an improved SVM ensemble classifier. Sensors 2018, 18, 1934. [Google Scholar] [CrossRef] [PubMed]
  24. Gu, F.; Shao, Y.; Hu, N.; Naid, A.; Ball, A.D. Electrical motor current signal analysis using a modified bispectrum for fault diagnosis of downstream mechanical equipment. Mech. Syst. Signal Process. 2011, 25, 360–372. [Google Scholar] [CrossRef]
  25. Zhang, R.; Gu, F.; Mansaf, H.; Wang, T.; Ball, A.D. Gear wear monitoring by modulation signal bispectrum based on motor current signal analysis. Mech. Syst. Signal Process. 2017, 94, 202–213. [Google Scholar] [CrossRef]
  26. Tian, X.; Gu, J.X.; Rehab, I.; M.Abdalla, G.; Gu, F.; Ball, A.D. A robust detector for rolling element bearing condition monitoring based on the modulation signal bispectrum and its performance evaluation against the Kurtogram. Mech. Syst. Signal Process. 2018, 100, 167–187. [Google Scholar] [CrossRef]
  27. Dequeiroz, R.L.; Rao, K.R. Time-varying lapped transforms and wavelet packets. IEEE Technol. Signal Process. 1993, 41, 3293–3305. [Google Scholar] [CrossRef]
  28. Paiva, H.M.; Galvão, R.K.H. Wavelet-packet identification of dynamic systems in frequency subbands. Signal Process. 2006, 86, 2001–2008. [Google Scholar] [CrossRef]
  29. Liu, B. Selection of wavelet packet basis for rotating machinery fault diagnosis. J. Sound Vib. 2005, 284, 567–582. [Google Scholar] [CrossRef]
  30. Kim, E.Y.; Lee, Y.J.; Lee, S.K. Heath monitoring of a glass transfer robot in the mass production line of liquid crystal display using abnormal operating sounds based on wavelet packet transform and artificial neural network. J. Sound Vib. 2012, 331, 3412–3427. [Google Scholar] [CrossRef]
  31. Wang, D.; Tse, P.W.; Tsui, K.L. An enhanced Kurtogram method for fault diagnosis of rolling element bearings. Mech. Syst. Signal Process. 2013, 35, 176–199. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the proposed method.
Figure 1. Flowchart of the proposed method.
Sensors 18 02908 g001
Figure 2. The planetary gearbox test rig.
Figure 2. The planetary gearbox test rig.
Sensors 18 02908 g002
Figure 3. The measurement system.
Figure 3. The measurement system.
Sensors 18 02908 g003
Figure 4. The fault modes: (a) sun gear chipped; (b) planetary gear bearing with inner-race fault.
Figure 4. The fault modes: (a) sun gear chipped; (b) planetary gear bearing with inner-race fault.
Sensors 18 02908 g004
Figure 5. The raw vibration signal of the sun gear chipped: (a) the waveform, and (b) frequency spectrum.
Figure 5. The raw vibration signal of the sun gear chipped: (a) the waveform, and (b) frequency spectrum.
Sensors 18 02908 g005
Figure 6. Energy distribution proportional histogram of the sun gear chipped.
Figure 6. Energy distribution proportional histogram of the sun gear chipped.
Sensors 18 02908 g006
Figure 7. The envelope spectrum of the reconstruction signal.
Figure 7. The envelope spectrum of the reconstruction signal.
Sensors 18 02908 g007
Figure 8. Results of the proposed method.
Figure 8. Results of the proposed method.
Sensors 18 02908 g008
Figure 9. The raw vibration signal of the faulty bearing: (a) the waveform and (b) spectrum.
Figure 9. The raw vibration signal of the faulty bearing: (a) the waveform and (b) spectrum.
Sensors 18 02908 g009
Figure 10. Energy distribution proportional histogram of the faulty bearing.
Figure 10. Energy distribution proportional histogram of the faulty bearing.
Sensors 18 02908 g010
Figure 11. The envelope spectrum of the reconstructed signal.
Figure 11. The envelope spectrum of the reconstructed signal.
Sensors 18 02908 g011
Figure 12. Results of the proposed method.
Figure 12. Results of the proposed method.
Sensors 18 02908 g012
Table 1. Structural parameters and kinematical parameters of the planetary gearbox.
Table 1. Structural parameters and kinematical parameters of the planetary gearbox.
GearNumber of Teeth Characteristic Frequency (Hz)
Sun gear1024.10 f s f
Planet gears(number)26(3)9.77 f p f
Ring gear623.89 f r f
Table 2. Specification and characteristic frequency of the planetary gear bearing.
Table 2. Specification and characteristic frequency of the planetary gear bearing.
Bearing DesignationBall Numbers
d (mm)
Pitch Diameter D m (mm)Ball Number
z
Contact Angle
β
12547.9 0 °
6008 f o f i f b f c
49.2565.1733.604.10
Table 3. The frequency range of 4-layer wavelet packet decomposition.
Table 3. The frequency range of 4-layer wavelet packet decomposition.
Serial NumberNode SituationFrequency/Hz
1stnode [4,0]0–6000 Hz
2ndnode [4,1]6000–12,000 Hz
3rdnode [4,2]12,000–18,000 Hz
4thnode [4,3]18,000–24,000 Hz
5thnode [4,4]24,000–30,000 Hz
6thnode [4,5]30,000–36,000 Hz
7thnode [4,6]36,000–42,000 Hz
8thnode [4,7]42,000–48,000 Hz
9thnode 4,8]48,000–54,000 Hz
10thnode [4,9]54,000–60,000 Hz
11thnode [4,10]60,000–66,000 Hz
12thnode [4,11]66,000–72,000 Hz
13thnode [4,12]72,000–78,000 Hz
14thnode [4,13]78,000–84,000 Hz
15thnode [4,14]84,000–90,000 Hz
16thnode [4,15]90,000–96,000 Hz

Share and Cite

MDPI and ACS Style

Guo, J.; Shi, Z.; Li, H.; Zhen, D.; Gu, F.; Ball, A.D. Early Fault Diagnosis for Planetary Gearbox Based Wavelet Packet Energy and Modulation Signal Bispectrum Analysis. Sensors 2018, 18, 2908. https://doi.org/10.3390/s18092908

AMA Style

Guo J, Shi Z, Li H, Zhen D, Gu F, Ball AD. Early Fault Diagnosis for Planetary Gearbox Based Wavelet Packet Energy and Modulation Signal Bispectrum Analysis. Sensors. 2018; 18(9):2908. https://doi.org/10.3390/s18092908

Chicago/Turabian Style

Guo, Junchao, Zhanqun Shi, Haiyang Li, Dong Zhen, Fengshou Gu, and Andrew D. Ball. 2018. "Early Fault Diagnosis for Planetary Gearbox Based Wavelet Packet Energy and Modulation Signal Bispectrum Analysis" Sensors 18, no. 9: 2908. https://doi.org/10.3390/s18092908

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop