Next Article in Journal
Benchmarking Particle Filter Algorithms for Efficient Velodyne-Based Vehicle Localization
Next Article in Special Issue
Automatic Fault Detection and Isolation Method for Roller Bearing Using Hybrid-GA and Sequential Fuzzy Inference
Previous Article in Journal
Decision Tree Pattern Recognition Model for Radio Frequency Interference Suppression in NQR Experiments
Previous Article in Special Issue
A Low-Delay Lightweight Recurrent Neural Network (LLRNN) for Rotating Machinery Fault Diagnosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Fault Feature Recognition Method for Time-Varying Signals and Its Application to Planetary Gearbox Fault Diagnosis under Variable Speed Conditions

1
The Key Laboratory of Metallurgical Equipment and Control of Education Ministry, Wuhan University of Science and Technology, Wuhan 430081, China
2
Hubei Key Laboratory of Mechanical Transmission and Manufacturing Engineering, Wuhan University of Science and Technology, Wuhan 430081, China
*
Authors to whom correspondence should be addressed.
Sensors 2019, 19(14), 3154; https://doi.org/10.3390/s19143154
Submission received: 10 June 2019 / Revised: 9 July 2019 / Accepted: 15 July 2019 / Published: 17 July 2019
(This article belongs to the Special Issue Sensors for Fault Diagnosis)

Abstract

:
The existing time-frequency analysis (TFA) methods mainly highlight the time-frequency ridges of the interested components by optimizing the time-frequency plane to facilitate the extraction of the relevant components. Generalized demodulation (GD), order tracking (OT), and other methods are generally used in conjunction with the TFA methods to realize the transition from a time-varying signal to a stationary signal, and finally identify the fault feature through a time-frequency plane. Generally, it is necessary to clarify the accuracy of the estimated components such as the rotational frequency or the fault characteristic frequency (FCF) during the operation of the GD or OT methods. Unfortunately, it is not only difficult to extract and locate rotational frequency or FCF, but also complicated in the whole estimation process. In this paper, a simple yet readable method is proposed to reveal the fault feature of time-varying signals. First, the method only needs to extract an arbitrary instantaneous frequency (IF). This is different from the GD method which needs to estimate and locate all phase functions. Then, it converts all variable frequency curves into corresponding lines parallel to the frequency axis based on the extracted IF to determine the proportional relationship between the components. Finally, to further improve the readability of the final results, we reduce the dimension of the transformed time-frequency representation to generate a two-dimensional (2D) energy-frequency map with high resolution and the same proportion. Subsequently, the performance is validated by simulated and experimental data.

1. Introduction

With the increasingly serious energy crisis and environmental pollution, as a green renewable energy, wind energy has attracted more and more attention throughout the world. As the main transmission component of wind turbines, planetary gearboxes are prone to failure due to long-term impact and alternating loads [1,2]. Therefore, accurate identification of gearbox fault characteristics is the prerequisite for fault diagnosis and health management of mechanical equipment [3,4].
Under the condition of variable speed, the signal of planetary gearbox has typical aperiodicity and non-stationarity. Undoubtedly, directly obtaining the fault characteristic frequencies of non-stationary signals by using the traditional Fourier spectrum analysis such fast Fourier transform is almost impossible [5,6]. Due to the aperiodic characteristics of non-stationary signals, the energy of the characteristic frequency is dispersed on the frequency axis, which is difficult to highlight the energy of a certain characteristic frequency on spectrum.
Time-frequency analysis (TFA) has always been a powerful tool for dealing with a non-stationary signal [7,8,9,10]. In recent decades, many TFA methods based on the short-time Fourier transform (STFT) and the continuous wavelet transform (CWT) have been proposed by scholars and applied in the field of mechanical fault diagnosis [11,12,13,14,15]. However, although the linear time-frequency analysis methods represented by the STFT and the CWT can be used to analyze non-stationary signals, time-frequency distribution with perfect time and frequency resolution cannot be obtained simultaneously, which is limited by the Heisenberg uncertainty principle [16,17]. Similarly, although the bilinear time-frequency analysis method represented by the Wigner–Ville distribution has effectively improved the resolution of time-frequency distribution, the analysis results of this method are prone to interference by cross terms and not suitable for analyzing multicomponent signals [17,18]. Moreover, due to the typical modulation characteristics of time-varying signals, energy divergence and mode aliasing are easy to occur when using traditional TFAs, which directly blur the time-frequency representation (TFR) and make it difficult to extract the accurate component information [19].
In order to achieve the concentrated time-frequency expression, I. Daubechies proposed the synchrosqueezed transform (SST), which is based on time-frequency coefficient rearrangement technology and has been applied in many fields, such as medicine, geophysics, and audio [20]. The SST is a special redistribution algorithm, which corresponds to a reallocation operator to improve the time-frequency resolution. In SST, the time-frequency coefficient is reallocated only on the frequency axis, which not only simplifies the distribution process but also keeps contact with the original parameters. Thakur et al. proposed the SST based on the STFT and this method had been successfully applied to ECG signals [21]. Presently, the SST has been widely used in the fault diagnosis of key parts of mechanical equipment [22,23,24,25]. However, the frequency reallocation operator in SST cannot provide an unbiased estimation for real-time accuracy, which makes it only sharpen the slowly changing signals and cannot meet the requirement of the high time-frequency resolution. For this reason, Pham et al. put forward the high-order synchrosqueezing transform (HSST) [26]. The HSST redefines the reallocation operator based on the high-order amplitude and phase approximation by using the Taylor expansion, which further improves the accuracy of the extracted instantaneous frequency (IF). Theoretically, the HSST is based on sharpening fast changing signals, however, the ability of HSST is related to the degree of the Taylor expansion. With the improvement of the order of the Taylor expansion, the corresponding expression of the algorithm will become complex. To make matters worse, the computation numbers of the STFT increase in the form of a square of the order, which greatly increases the calculation cost.
Recently, the proposed multi-synchrosqueezing transform (MSST) method has a better performance in terms of further sharpening the TFR than the HSST [27]. The MSST method introduces a progressive idea into the SST, and therefore makes it possible to estimate IF more accurately through repeated estimation of the reallocation operator, since this method only optimizes the reallocation operator and the whole process only needs to calculate the STFT once. With the increase of iterations, the MSST method can obtain time-frequency representation with higher resolution than the HSST. Although the HSST and MSST methods achieve further energy concentration, it is still difficult to intuitively obtain interesting information from the TFR in the background noise, for the same theoretical basis.
Generalized demodulation (GD) [28,29,30] transforms the IF of any time-varying signal into a constant frequency, and each component signal is extracted by means of bandpass filtering, which largely makes up for the deficiency of time-frequency analysis methods based on the STFT and the CWT in signal processing. However, GD needs to estimate the phase function of all interested components, which is very complicated in practical use. Similarly, order tracking (OT) [31,32] is a powerful tool commonly used in non-stationary signal processing. First, this approach estimates the shaft frequency accurately, then, realizes the signal stabilization, and finally, presents in the form of an order spectrum. However, the angle domain resampling method is adopted in OT, which inevitably leads to the problems of missing values and interpolation. In addition, due to the lack of noise reduction operation, the background noise is still retained in the order spectrum.
In this paper, a new method of fault feature recognition for time-varying signals is proposed. The main purpose of this method is to simplify the analysis process of GD and present more explicit and readable results. For this reason, we focus on the following:
  • Can the process of estimating and locating all phase functions in GD with only one instantaneous frequency be simplified and achieve similar smooth transformation?
  • Can the process be further simplified? For example, can the instantaneous frequency that needs to be estimated be arbitrary?
  • Can more intuitive and readable results on the basis of TFR be presented?
In order to realize this idea, we mainly improve the following algorithms:
  • To obtain stable high resolution TFR, an improved multi-synchrosqueezing transform (IMSST) algorithm based on MSST is proposed. According to the intermediate frequency extracted in the first step, the time-varying full frequency can be directly converted into stable full frequency.
  • To transform the instantaneous frequency ridges into a series of lines parallel to the frequency axis, we improve the instantaneous frequency estimation operator based on the MSST algorithm, so that we achieve the result which is similar to GD using only an arbitrary extracted IF.
  • To present more intuitive and readable results, we propose a simple data dimension reduction method, which generates a more readable two-dimensional (2D) energy-frequency diagram.
This paper contributes to the following three major aspects:
  • A three-step model is used to enhance readability of the final results.
  • The proposed method is validated in multicomponent and planetary gearbox simulation signals in the form of increasing signal complexity.
  • The proposed method is applied to diagnose the planetary gearbox fault under a time-varying condition and directly recognize the fault type from the 2D energy-frequency map without using any other method.

2. Theoretical Description

According to the previous ideas, the proposed method is mainly divided into three steps: firstly, extract an arbitrary instantaneous frequency; then, obtain steady and high resolution TFR from time-varying signals; finally, produce a 2D energy-frequency map. Therefore, in this section, we will introduce the relevant theories of these three steps in turn.

2.1. Extraction Algorithm of IF

At present, the peak searching algorithm is the most common algorithm in IF extraction. It obtains the IF of interested components by searching the peak of time-frequency ridge energy. However, due to the interference of modulation phenomenon and strong noise in measured signals, it is generally difficult to extract high accuracy time-frequency ridges directly.
In Ref. [33], G. Yu proposed a general linear chirplet transform (GLCT). He used a time-varying demodulation operator to produce a rotational effect on the time-frequency plane, so as to solve the problem of energy divergence caused by modulation. The GLCT constructs discrete time-varying demodulation operators by defining rotating degree α which is written as follows:
{ α = arctan ( T f s c ) ,   α ( π 2 , π 2 ) α = π 2 + π N + 1 , π 2 + 2 π N + 1 , , π 2 + N π N + 1 , ( N = 1 , 2 , )
where, T is sampling time and f s is sampling frequency. N denotes the number of discrete time-varying demodulation operator c . From Equation (1), the rotating degree α is equally divided into N parts. However, considering the difference of the instantaneous change trend of general time-varying signals, it is difficult to completely represent the change trend of the time-frequency curve only through the division of equal intervals.
In this paper, we abandon the above estimation form and reconstruct the time-varying demodulation operator directly by estimating the trend of the instantaneous frequency. In this way, the energy concentration of the instantaneous frequency curve is maintained to the greatest extent.
An IF extraction algorithm based on time-varying demodulation and the STFT (TDSTFT) is proposed, which combines the idea of time-varying demodulation and the peak searching algorithm. It repeatedly estimates the time-varying demodulation operator through the changing trend of IF to optimize the TFR, and then extracts ridge line using the peak searching algorithm. The following is a detailed theoretical introduction.
An AM-FM signal is defined as:
s ( t ) = A ( t ) e i Φ ( t ) = A ( t ) e i φ ( t ) d t
where, A(t) and Φ(t) are their instantaneous amplitude and instantaneous phase, respectively. And φ ( t ) = Φ ( t ) is regarded as the ideal IF of the analytical signal. In a short time u , the IF of signal s ( t ) can be approximately regarded as a linear equation, which is expressed by the first-order Taylor expansion:
φ ( u ) = φ ( t ) + φ ( t ) ( u t )
Then, the corresponding STFT is expanded into the following form:
| S ( φ ( t ) | = | + g ( u t ) s ( u ) e i φ ( t ) u d u | = | + g ( u t ) A ( u ) e i [ φ ( t ) u + φ ( t ) ( u t ) 2 / 2 ] e i φ ( t ) u d u | = | + e i φ ( t ) ( u t ) 2 / 2 g ( u t ) A ( u ) e i φ ( t ) u e i φ ( t ) u d u | | + g ( u t ) A ( u ) d u |
where, g ( t ) is a window function. Here we choose the Gauss window, whose expression is written as:
g ( t ) = 1 2 π σ e 1 2 σ 2 t 2
According to Equation (4), the signal is modulated and the energy of STFT diverges due to the presence of e i φ ( t ) ( u t ) 2 / 2 . Therefore, we introduce the time-varying demodulation operator e i c ˜ ( t ) ( u t ) 2 / 2 to eliminate the influence of e i φ ( t ) ( u t ) 2 / 2 , which makes the TFR energy more concentrated. After introducing the time-varying demodulation operator, the STFT is rewritten as:
T D S T F T = + g ( u t ) s ( u ) e i φ ( t ) u e i c ˜ ( t ) ( u t ) 2 / 2 d u
Because the c ˜ ( t ) is unknown, here, we adopt the approximation strategy and achieve the accuracy allowed by the algorithm through the change trend of the approximate IF step by step. Considering the uniformity of energy concentration, this paper directly uses the change trend of instantaneous frequency to express c ˜ ( t ) .
c ˜ ( t ) = φ ( t )
The specific process of the whole extraction algorithm (Algorithm 1) is as follows:
Algorithm 1 TDSTFT
Input: s ( t )
Output: c ˜ ( t ) ,   T D S T F T
1:Initialize: c ˜ ( t ) = 0 ; the largest number of iterations: n = 1 , 00 ; convergence threshold ε .
2:Calculate:
     T D S T F T = + g ( u t ) s ( u ) e i φ ( t ) u e i c ˜ ( t ) ( u t ) 2 / 2 d u ;
     L t = l e n g t h ( s ( t ) ) ;
     P = p e a k s e a r c h i n g ( T D S T F T ) ;
     I ( t ) = p o l y f i t ( P ) ;
     I ( t ) = p o l y f i t ( P ) ;
     I ( t ) = d i f f ( I ( t ) ) ;
3:Iteration:
    for j = 1 : n ;
      I i ( t ) = I ( t ) ;
      c ˜ ( t ) = I ( t ) ;
      T D S T F T = + g ( u t ) s ( u ) e i φ ( t ) u e i c ˜ ( t ) ( u t ) 2 / 2 d u ;
      P = p e a k s e a r c h i n g ( T D S T F T ) ;
      I ( t ) = p o l y f i t ( P ) ;
      I ( t ) = d i f f ( I ( t ) ) ;
      ξ = 1 L t 0 L t | I ( t ) I i ( t ) | ;
     if ξ < ε ;
      breake;
     end
    end
The IF of a component can be obtained more accurately after the above steps. The following is an example to illustrate the precision of the IF extraction algorithm.
We use a numerical example to illustrate the effectiveness of this method and the simulated signal is s ( t ) = sin [ 2 π ( 10 t + 10 t 2 19 9 t 3 + 9 50 t 4 1 200 t 5 ) ] + n ( t ) , where n ( t ) is −2 dB while Gaussian noise and the sampling frequency is set to 200 Hz. Figure 1a,c shows the TFRs by STFT and TDSTFT after four iterations, respectively. The energy of TFR, as shown in Figure 1c, is more concentrated. Figure 1b,d shows the extracted IFs obtained from Figure 1a,c, respectively, using the peak searching algorithm (where, the blue line is the extracted IF and the red line is the real value). Through the comparative analysis, it is illustrated that the IF extracted from graph d in Figure 1 is obviously closer to the true value, which also proves the effectiveness of the proposed method.

2.2. The TFRs by IMSST

According to the expression of time-varying signals in Equation (2), a standard STFT is written as follows:
V f g ( t , η ) = s ( τ ) g ( τ t ) e i η ( τ t ) d τ
where, g ( t ) is the complex common of g ( t ) .
Then, the SST based on the STFT is expressed as:
T s ( t , ω ) = 1 g ( 0 ) { η , | V f g ( t , η ) | > γ } ( t , η ) ( t , ω ) V f g ( t , η ) δ ( ω ω ^ f ( t , η ) ) d η
where, γ is a very small threshold to prevent excessive error when V f g ( t , η ) = 0 . And δ function is defined as:
δ ( ω ω f ) = { 1 , i f   ω = ω f 0 , f   ω ω f
In Equation (9), ω f ( t , η ) represents the IF estimation of signals on t and η . Since the phase after STFT is not affected by the window function g ( t ) , we use V f g ( t , η ) to estimate the IF:
ω f ( t , η ) = t arg [ V f g ( t , η ) ] = t V f g ( t , η ) i V f g ( t , η )
where, arg [ A ] denotes the angle of the complex A . Considering that ω f ( t , η ) is a complex number, we only take the real part as the estimation of IF.
ω ^ f ( t , η ) = Re { ω f ( t , η ) } = Re { t V f g ( t , η ) i V f g ( t , η ) }
Although SST cannot deal with strong time-varying signals well, TFR with high concentration can also be obtained by continuously utilizing SST concentrating energy. This is the core idea of MSST method.
T s [ 2 ] ( t , ω ) = T s [ 1 ] ( t , η ) δ ( ω ω ^ f ( t , η ) ) d η T s [ 3 ] ( t , ω ) = T s [ 2 ] ( t , η ) δ ( ω ω ^ f ( t , η ) ) d η T s [ N ] ( t , ω ) = T s [ N 1 ] ( t , η ) δ ( ω ω ^ f ( t , η ) ) d η
when N = 2 ,
T s [ 2 ] ( t , ω ) = T s [ 1 ] ( t , ξ ) δ ( ξ ω ^ f ( t , ξ ) ) d ξ = V f g ( t , η ) δ ( ω ω ^ f ( t , η ) ) d η δ ( ω ω ^ f ( t , ξ ) ) d ξ = V f g ( t , η ) δ ( ω ω ^ f ( t , η ) ) δ ( ω ω ^ f ( t , ξ ) ) d ξ d η = V f g ( t , η ) δ ( ω ω ^ f ( t , ω ^ f ( t , η ) ) ) d η
After the second iteration, the original time-frequency plane is changed from ( t , η ) to ( t , ω ^ f ( t , ω ^ f ( t , η ) ) ) .
In a short time, the phase Φ ( τ ) of the original signal s ( τ ) is expanded as:
Φ ( τ ) = Φ ( t ) + Φ ( t ) ( τ t ) + 1 2 Φ ( t ) ( τ t ) 2
The purpose of (15) is to simplify the STFT coefficient V f g ( t , η ) and V f g ( t , η ) is rewritten as follows:
V f g ( t , η ) = A ( τ ) e i [ Φ ( t ) + Φ ( t ) ( τ t ) + 1 2 Φ ( t ) ( τ t ) 2 ] g ( τ t ) e i η ( τ t ) d τ = A ( τ ) e i Φ ( t ) 1 1 i Φ ( t ) e ( η Φ ( t ) ) 2 2 ( 1 i Φ ( t ) )
Similarly, the IF estimation is expressed as:
ω ^ f ( t , η ) = Re { t V f g ( t , η ) i V f g ( t , η ) } = Φ ( t ) + ( Φ ( t ) ) 2 1 + ( Φ ( t ) ) 2 ( η Φ ( t ) )
When N = 2 ,
ω ^ f ( t , ω ^ f ( t , η ) ) = Φ ( t ) + ( Φ ( t ) ) 2 1 + ( Φ ( t ) ) 2 ( ω ^ f ( t , η ) Φ ( t ) ) = Φ ( t ) + [ ( Φ ( t ) ) 2 1 + ( Φ ( t ) ) 2 ] 2 ( η Φ ( t ) )
When N = 3 ,
ω ^ f [ 3 ] ( t , η ) = ω ^ f ( t , ω ^ f ( t , ω ^ f ( t , η ) ) ) = Φ ( t ) + [ ( Φ ( t ) ) 2 1 + ( Φ ( t ) ) 2 ] 2 ( ω ^ f ( t , η ) Φ ( t ) ) = Φ ( t ) + [ ( Φ ( t ) ) 2 1 + ( Φ ( t ) ) 2 ] 2 [ ( Φ ( t ) ) 2 1 + ( Φ ( t ) ) 2 ( η Φ ( t ) ) ] = Φ ( t ) + [ ( Φ ( t ) ) 2 1 + ( Φ ( t ) ) 2 ] 3 ( η Φ ( t ) )
In order to facilitate the writing, here, we use ω ^ f [ 3 ] ( t , η ) instead of ω ^ f ( t , ω ^ f ( t , ω ^ f ( t , η ) ) ) . Similarly, the estimation of IF ω ^ f ( t , ω ^ f ( t , ω ^ f ( t , ω ^ f ( t , ) ) ) ) after the Nth iteration is written as ω ^ f [ N ] ( t , η ) . According to the law, it is not difficult to deduce the estimation of IF after the Nth iteration.
ω ^ f [ N ] ( t , η ) = Φ ( t ) + [ ( Φ ( t ) ) 2 1 + ( Φ ( t ) ) 2 ] N ( η Φ ( t ) )
Therefore, MSST is expressed as:
T s [ N ] ( t , ω ) = 1 g ( 0 ) { η , | V f g ( t , η ) | > γ } ( t , η ) ( t , ω ^ f [ N ] ( t , η ) ) V f g ( t , η ) δ ( ω ω ^ f [ N ] ( t , η ) ) d η
To illustrate the advantages of MSST in sharpening time-frequency domain, we take a simple time-varying signal, s ( t ) = 10 cos [ 2 π ( 500 t 2 e 2 t sin ( 20 π t ) ) ] as an example, and the sampling frequency of this signal is set as 4096 Hz. Figure 2 shows its TFRs obtained by SST, 2-SST, 4-SST and MSST ( N = 6 ), respectively. And their resolution is quantified by Rényi entropy shown in Table 1. In general, the smaller the value of Renyi entropy [34], the higher the resolution obtained. It proved that MSST method has the smallest value and produces more energy-intensive TFRs as shown in Table 1.
In order to obtain stable and high resolution TFR, we directly improve MSST by referring to the idea of demodulation, so that the final TFR is a series of lines parallel to the frequency axis.
Suppose we obtained an arbitrary instantaneous frequency f I F ( t ) using TDSTFT and then the proportional coefficient at any time can be expressed as:
K = f I F ( i ) f I F ( t )
where, f I F ( i ) represents the frequency value at the time i . Generally, fIF(i) is chosen among the four frequencies: initial frequency, cutoff frequency, maximum frequency, and minimum frequency, i.e., f I F ( i ) = f I F ( 0 ) , f I F ( i ) = f I F ( e n d ) , f I F ( i ) = f I F max ( t ) , f I F ( i ) = f I F min ( t ) .
According to the proportional coefficient K , the transformed IF is written as:
ω ˜ f [ N ] ( t , ω ) = ω ^ f [ N ] ( t , K × ω )
According to (22), the IFs are redistributed on the time-frequency plane to obtain the target time-frequency representation.
Finally, the expression of IMSST is written as follows:
I T s [ N ] ( t , ω ) = 1 g ( 0 ) { η , | V f g ( t , η ) | > γ } ( t , η ) ( t , ω ˜ f [ N ] ( t , η ) ) V f g ( t , η ) δ ( ω ω ˜ f [ N ] ( t , η ) ) d η

2.3. Two-Dimensional Energy-Frequency Map Obtained by IMSST

In the fault diagnosis of time-varying signals, scholars not only expect to retain as much original information as possible, but also hope to be able to produce familiar analysis graphs to directly judge, such as the Fourier spectrum, order spectrum, and so on. We can obtain high resolution time-frequency plane using IMSST. In order to further clarify the proportional relationship between signal components, it is necessary to eliminate the influence of time and transform the TFR into a 2D frequency-energy map.
The specific conversion is as follows:
I T s r e ( ω j ) = j = 1 M | I T s [ N ] ( t j , ω ) | M
where, I T s r e ( ω j ) represents the amplitude (energy) at frequency ω j and M is the number of samples, i.e., the length of the original time series. By (25), the influence of t is eliminated and the corresponding 2D energy-frequency map is obtained.
The flowchart of the proposed method is shown in Figure 3, which is also summarized as follows:
  • Apply acceleration sensor to obtain the vibration signal s ( t ) .
  • Extract an arbitrary IF f I F ( t ) from s ( t ) using TDSTFT.
  • Calculate the proportionality coefficient K according to f I F ( t ) .
  • Obtain the high resolution of TFR T s [ N ] ( t , ω ) using MSST.
  • Transforming instantaneous frequency curves to a series of lines which are parallel to the frequency axis according to the proportionality coefficient K using IMSST, and then obtaining reconstructed TFR I T s [ N ] ( t , ω ) .
  • Eliminate the influence of time t to obtain I T s r e ( ω j ) which is a result of dimensionality reduction comparing with I T s [ N ] ( t , ω ) , and then derive the 2D energy-frequency map.
  • Identify the fault pattern using the 2D energy-frequency map.

3. Simulation Validation

In this section, we describe a multicomponent signal and a planetary gear simulation signal which are employed to illustrate the effectiveness of the proposed method.

3.1. Multicomponent Signal

To better illustrate the principle of this method, a simple time-varying multicomponent signal is established in this section, which consists of four sinusoidal functions with proportional frequencies.
{ s ( t ) = k = 1 4 A k cos ( 2 π 0 t k × f ( t ) d t ) + n ( t ) f ( t ) = 12 sin ( 0.35 t + π 4 ) + 3 sin ( 2 t 4 π 5 ) + 1.62 sin ( 3.5 t 7 π 3 ) + 50
where, f ( t ) is the IF and n ( t ) is the white Gaussian noise. In this signal, the sampling frequency is 1000 Hz and the sampling number is 6000. Similarly, a SNR of −5 dB is added to form a synthetic signal, as shown in Figure 4a. The Fourier spectrum of the synthetic signal is shown in Figure 4b. For the influence of variable speed, the energy of Fourier spectrum is divergent and blurry, which means the useful information is hard to find.
For comparison, the TFR results of the signal obtained using STFT, MSST, and IMSST are presented in Figure 4c–e. Due to the influence of the Heisenberg uncertainty principle [17], the TFR of STFT has the characteristics of a large frequency bandwidth and a low frequency resolution. Especially in noise background, the whole time-frequency plane is basically blurred, which corresponds to the result in Figure 4c. Figure 4d,e are TFRs obtained using MSST and IMSST, respectively. In comparison to Figure 4c, Figure 4d and e shows sharper time-frequency ridges. Moreover, Figure 4e transforms the time-varying instantaneous frequency ridge into a line parallel to the frequency axis on the basis of retaining the high resolution of Figure 4d, which makes the proportional relationship between the components clearer. Figure 4f is a 2D energy-frequency map derived from the TFR of Figure 4e which shows the four proportional frequencies: 20.33 Hz, 40.67 Hz, 61 Hz, and 81.67 Hz, which correspond to the integral times of f ( 0 ) = 20.32 Hz basically. There is no doubt that Figure 4f is more readable than Figure 4e.
For further comparison, we draw the final results of the OT method and the proposed method in Figure 5. It should be noted that the result of the OT method is obtained through three steps: firstly, obtain an IF using the peak searching algorithm on the TFR by STFT. Secondly, apply the angle domain resampling technology to obtain the stable angle domain signal. Thirdly, draw the Fourier spectrum of the angle domain signal. As compared with the results of OT and the proposed method, we see that the latter not only highlights the feature components, but also reduces the influence of noise, which indicates that the proposed method produces clearer results.

3.2. The Simulated Planetary Gear Signal under Time-Varying Rotating Speed

To evaluate the effectiveness of the proposed method for analyzing planetary gear vibration signals, the simulated planetary gear signals are generated according to [35]:
s ( t ) = [ 1 cos ( 2 π f s ( r ) ( t ) d t ) ] [ 1 + A 1 cos ( 2 π f s s ( t ) d t ) + θ 1 ] × cos [ 2 π f m ( t ) d t + A 2 sin ( f s s ( t ) d t + θ 2 ) + θ 3 ] + n ( t )
where, f s ( r ) ( t ) is the instantaneous rotational frequency of sun gear and f s s ( t ) is the instantaneous fault characteristic frequency (IFCF). f m ( t ) denotes the instantaneous meshing frequency of the simulated signal. A 1 and A 2 are expressed as modulation factors. θ 1 , θ 2 and θ 3 are initial phases. n ( t ) corresponds to the while Gaussian noise. All parameters are shown in Table 2.
Considering the variable speed, the instantaneous rotation frequency f s ( r ) ( t ) is written as:
f s ( r ) ( t ) = 20 t 2 + 20 t + 27
Figure 6a,b shows the time domain waveforms and the Fourier spectrum of the simulated signals, respectively. The TFR obtained using STFT is shown in Figure 6c, where we can roughly see that the energy is concentrated between 400 and 600 Hz. However, it is difficult to further distinguish the meshing frequency and its sideband frequency. Figure 6d is the extracted IF by TDSTFT, which is basically close to the true value (where the blue line is the extracted IF and the red line is the real IF). Figure 6e,f are TFRs obtained using MSST and IMSST, respectively. As comparing with Figure 6d, they are clearer to find the meshing frequency and its sideband frequency.
To better illustrate the advantage of the method, we respectively generate the 2D energy-frequency maps by SST and 4-SST base on the same algorithmic theory in Section 2.2 and Section 2.3, in Figure 7. Due to the aliasing of each component, it is difficult to capture useful peaks from the graph. This is mainly because two methods cannot provide enough clear TFRs, and therefore it results in blurred 2D energy-frequency maps. Figure 8 is the 2D energy-frequency map using IMSST ( N = 20 ) of a planetary gear simulated signal, which clearly shows the meshing frequency and its corresponding sideband. In comparison to Figure 6b, the phenomenon of frequency divergence has been greatly improved. Furthermore, the proportional relationship among the components exhibits more clearly in the energy-frequency map as comparing with those in Figure 7. Moreover, the results perfectly realize the transformation from non-stationary signal to stationary signal, which is also convenient for fault diagnosis.
For further comparison, we also draw the results of the OT method and the proposed method in Figure 9 which shows that the proposed method is better for highlighting fault characteristic frequency and eliminating background noise.
To further testify the robustness of the proposed method, more simulated signals under different noise levels are analyzed. Figure 10 presents the TFRs using IMSST and the corresponding 2D energy-frequency maps of the simulated signals with SNR of 0 dB, −5 dB, and −10 dB, respectively. With the enhancement of noise, it becomes more and more difficult to directly recognize the interest components from TFR using IMSST. Correspondingly, the peak value of characteristic frequencies becomes less prominent in the 2D energy-frequency maps. Therefore, when the noise SNR of a simulated signal is greater than −5 dB, the proposed method performs well with respect to fault feature recognition. If the noise level is too high, the effectiveness of the proposed method is not ideal. It is worth noting that, when the signal noise is very strong, it can be combined with the existing singular spectrum analysis (SSA), convex optimization, time-varying filter, and other denoising methods to reduce the noise level of the signal.

4. Experimental Validation

As we all know, the simulated signal is usually a vibration model based on the meshing characteristics of gearbox which is a simplified AM-FM model. The real acquisition signal contains not only the coupling vibration information of key components such as gears, bearings, and so on, but also the other interference information caused by background noise and vibration of other components of the mechanical equipment. It is more complex than the simulation signal. Therefore, it is necessary to verify the proposed method by experiment.

4.1. Experimental Rig

In this section, the proposed method is applied to detect the gear damage. The experiment is carried out on the test rig shown in Figure 11. The experimental rig is composed of a driving motor, planetary gear box, helical gearbox, magnetic powder brake, and other components. The internal structure of the test rig is shown in Figure 12. The test fault is located on the sun gear of planetary gearbox, as shown in Figure 13. A PCB accelerometer is mounted on the top of the planetary gearbox to measure the vibration signal.
In this experiment, the rotational speed of the motor decreases linearly from 3000 r/min with a sampling frequency of 2560 Hz and a sampling number of 7680. The applied torque by the magnetic powder brake is 0. The parameters of the planetary gearbox are shown in Table 3, and the corresponding fault characteristic frequencies are calculated with reference to [36], as shown in Table 4.

4.2. Result Analysis

Figure 14a,b are the time domain waveform and the Fourier spectrum of the fault vibration signals of the sun gear, respectively. As shown in Figure 14a, the impact decreases with a decrease of rotational speed. Due to the influence of variable speed, it is still difficult to identify the fault type directly from the Fourier spectrum. Figure 14c shows the TFR using STFT, where we can roughly see the changing trend of the vibration signal. f I F ( t ) is the IF extracted using TDSTFT (see Figure 14d), which shows that the speed of the vibration signal basically decreases in a linear way, which is also in line with our Introduction.
According to the extracted IF, select f I F ( 0 ) as the molecule of the proportional coefficient K, and then we obtain K = f I F ( 0 ) f I F ( t ) . At this time, the variable speed signal will be converted to a constant frequency according to the proportional relationship. Then the rotation frequency of the sun gear in the energy-frequency map is f s ( r ) = 23.5 Hz, and the other frequencies are calculated according to Table 4.
Generally, if the three planetary gears meshing with the sun gear are identical, the FCF of the sun gear is expressed as f s s = N f m / Z s , where N is the number of planetary gears and Z s is the number of sun gear teeth. And the peak value of sideband usually appears at the frequency of f m ± k f s ( r ) + n f s s . On the basis of the equation f s s = N ( f s ( r ) f c ) , the above frequency can also be written as f m ± k f c + ( n ± k / N ) f s s . However, manufacturing error is inevitable in the actual gearbox, which means the three planetary gears meshing with the sun gear cannot be exactly identical, leading to the difference of the impact sequence between the fault teeth of the sun gear and three planetary gears. If these impacts are considered as three different sequences, the FCF of the sun gear can be written as f s s / N . Correspondingly, the peak value of the sideband appears at frequency f m ± k f c + ( ( n ± k ) / N ) f s s . Figure 15a is the final 2D energy-frequency map, where the rotational frequency f s ( r ) of the sun gear; the rotational frequency f c of the planetary frame; and the meshing frequency, f m and 2 f m , 3 f m are easily found. Figure 15b is an enlarged view of Figure 13 near the frequency fm and shows that the corresponding peak values are related to one-third times the characteristic frequency of the sun gear fault characteristic frequency, which could be caused by the manufacturing or installation errors of the three planetary gears. However, this does not prevent us from judging that the fault of the planetary gearbox is the broken teeth of sun gear.

5. Conclusions

In this paper, a novel fault diagnosis method for time-varying signals was proposed, and the performance of this method has been verified using simulation signals and a planetary gearbox experiment. Its advantages are as follows:
Firstly, this method does not need to estimate the instantaneous shaft frequency. It only needs to estimate an arbitrary instantaneous frequency with time-varying characteristics, such as gear meshing frequency, which is easier than extracting the instantaneous shaft frequency.
Secondly, this method can smooth the time-varying frequency ridges and achieve the effect similar to GD and OT. Moreover, this method omitted the steps of filtering and resampling, so that it can retain the complete information of the vibration signal.
Finally, this method can generate a 2D energy-frequency map, which is similar to the spectrum map. In addition, in comparison to direct analysis on time-frequency plane, this form is more intuitive and friendly.

Author Contributions

B.P. and Y.L. conceived method; B.P. and C.Y. designed the experiments; B.P. and Y.M. performed the experiments; B.P. and C.Y. analyzed the data; M.G., Y.L. and C.Y. contributed reagents/materials/analysis tools; B.P. and Y.M. wrote the paper.

Funding

This work was supported by the National Natural Science Foundation of China under Grants No. 51875416 and 51805382.

Conflicts of Interest

The authors declare that no conflict of interests regarding the publication of this article.

References

  1. Chen, B.; Zhang, Z.; Hua, X.; Basu, B.; Nielsen, S.R. Identification of aerodynamic damping in wind turbines using time-frequency analysis. Mech. Syst. Signal Process. 2017, 91, 198–214. [Google Scholar] [CrossRef] [Green Version]
  2. Yi, C.; Lv, Y.; Dang, Z.; Xiao, H.; Yu, X. Quaternion singular spectrum analysis using convex optimization and its application to fault diagnosis of rolling bearing. Measurement 2017, 103, 321–332. [Google Scholar] [CrossRef]
  3. Lv, Y.; Yuan, R.; Wang, T.; Li, H.; Song, G. Health degradation monitoring and early fault diagnosis of a rolling bearing based on CEEMDAN and improved MMSE. Materials 2018, 11, 1009. [Google Scholar] [CrossRef] [PubMed]
  4. Yi, C.; Lv, Y.; Ge, M.; Xiao, H.; Yu, X. Tensor singular spectrum decomposition algorithm based on permutation entropy for rolling bearing fault diagnosis. Entropy 2017, 19, 139. [Google Scholar] [CrossRef]
  5. Qiao, W.; Lu, D. A Survey on Wind Turbine Condition Monitoring and Fault Diagnosis—Part II: Signals and Signal Processing Methods. IEEE Trans. Ind. Electron. 2015, 62, 6546–6557. [Google Scholar] [CrossRef]
  6. Tu, X.; Hu, Y.; Li, F.; Abbas, S.; Liu, Z.; Bao, W. Demodulated High-Order Synchrosqueezing Transform With Application to Machine Fault Diagnosis. IEEE Trans. Ind. Electron. 2019, 66, 3071–3081. [Google Scholar] [CrossRef]
  7. Feng, Z.; Chu, F.; Zuo, M.J. Time–frequency analysis of time-varying modulated signals based on improved energy separation by iterative generalized demodulation. J. Sound Vib. 2011, 330, 1225–1243. [Google Scholar] [CrossRef]
  8. Huang, H.; Baddour, N.; Liang, M. Bearing fault diagnosis under unknown time-varying rotational speed conditions via multiple time-frequency curve extraction. J. Sound Vib. 2018, 414, 43–60. [Google Scholar] [CrossRef]
  9. Chen, S.; Yang, Y.; Wei, K.; Dong, X.; Peng, Z.; Zhang, W. Time-varying frequency-modulated component extraction based on parameterized demodulation and singular value decomposition. IEEE Trans. Instrum. Meas. 2016, 65, 276–285. [Google Scholar] [CrossRef]
  10. Liu, D.; Cheng, W.; Wen, W. Rolling Bearing Fault Diagnosis via ConceFT-Based Time-Frequency Reconfiguration Order Spectrum Analysis. IEEE Access 2018, 6, 67131–67143. [Google Scholar] [CrossRef]
  11. Kwok, H.K.; Jones, D.L. Improved instantaneous frequency estimation using an adaptive short-time Fourier transform. IEEE Trans. Signal Process. 2000, 48, 2964–2972. [Google Scholar] [CrossRef]
  12. Ta, R.; Li, Y.L.; Wang, Y. Short-time fractional Fourier transform and its applications. IEEE Trans. Signal Process. 2010, 58, 2568–2580. [Google Scholar] [CrossRef]
  13. Hyvärinen, A.; Ramkumar, P.; Parkkonen, L.; Hari, R. Independent component analysis of short-time Fourier transforms for spontaneous EEG/MEG analysis. NeuroImage 2010, 49, 257–271. [Google Scholar] [CrossRef] [PubMed]
  14. Liu, W.; Cao, S.; Chen, Y. Seismic Time-Frequency Analysis via Empirical Wavelet Transform. IEEE Geosci. Remote Sens. Lett. 2016, 13, 28–32. [Google Scholar] [CrossRef]
  15. Jaganathan, K.; Eldar, Y.C.; Hassibi, B. Recovering signals from the short-time Fourier transform magnitude. In Proceedings of the 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Brisbane, QLD, Australia, 19–24 April 2015; pp. 3277–3281. [Google Scholar]
  16. Yu, G.; Wang, Z.; Zhao, P.; Zhen, L. Local maximum synchrosqueezing transform: an energy-concentrated time-frequency analysis tool. Mech. Syst. Sig. Process. 2019, 117, 537–552. [Google Scholar] [CrossRef]
  17. Kalyani, A.; Pachori, R.B. Cross-Terms Free Time-Frequency Representation Using Empirical Wavelet Transform and Wigner-Ville Distribution. Ph.D. Thesis, Discipline of Electrical Engineering, IIT Indore, West Bengal, India, 2018. [Google Scholar]
  18. 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]
  19. Canal, M.R. Comparison of Wavelet and Short Time Fourier Transform Methods in the Analysis of EMG Signals. J. Med. Syst. 2010, 34, 91–94. [Google Scholar] [CrossRef]
  20. Daubechies, I.; Lu, J.; Wu, H.T. Synchrosqueezed wavelet transforms: An empirical mode decomposition-like tool. Appl. Comput. Harm. Anal. 2011, 30, 243–261. [Google Scholar] [CrossRef] [Green Version]
  21. Thakur, G.; Wu, H.T. Synchrosqueezing-based recovery of instantaneous frequency from nonuniform samples. SIAM J. Math. Anal. 2011, 43, 2078–2095. [Google Scholar] [CrossRef]
  22. Li, C.; Liang, M. Time–frequency signal analysis for gearbox fault diagnosis using a generalized synchrosqueezing transform. Mech. Syst. Signal Process. 2012, 26, 205–217. [Google Scholar] [CrossRef]
  23. Feng, Z.; Chen, X.; Liang, M. Iterative generalized synchrosqueezing transform for fault diagnosis of wind turbine planetary gearbox under nonstationary conditions. Mech. Syst. Signal Process. 2015, 52, 360–375. [Google Scholar] [CrossRef]
  24. Fourer, D.; Auger, F.; Czarnecki, K.; Meignen, S.; Flandrin, P. Chirp Rate and Instantaneous Frequency Estimation: Application to Recursive Vertical Synchrosqueezing. IEEE Signal Process. Lett. 2017, 24, 1724–1728. [Google Scholar] [CrossRef]
  25. Guo, Y.; Chen, X.; Wang, S.; Sun, R.; Zhao, Z. Wind Turbine Diagnosis under Variable Speed Conditions Using a Single Sensor Based on the Synchrosqueezing Transform Method. Sensors 2017, 17, 1149. [Google Scholar] [Green Version]
  26. Pham, D.H.; Meignen, S. High-Order Synchrosqueezing Transform for Multicomponent Signals Analysis-With an Application to Gravitational-Wave Signal. IEEE Trans. Signal Process. 2017, 65, 3168–3178. [Google Scholar] [CrossRef]
  27. Yu, G.; Wang, Z.; Zhao, P. Multi-synchrosqueezing Transform. IEEE Trans. Ind. Electron. 2018, 66, 5441–5455. [Google Scholar] [CrossRef]
  28. Li, C.; Liang, M. A generalized synchrosqueezing transform for enhancing signal time–frequency representation. Signal Process. 2012, 92, 2264–2274. [Google Scholar] [CrossRef]
  29. Shi, J.; Liang, M.; Necsulescu, D.S.; Guan, Y. Generalized stepwise demodulation transform and synchrosqueezing for time–frequency analysis and bearing fault diagnosis. J. Sound Vib. 2016, 368, 202–222. [Google Scholar] [CrossRef]
  30. Li, C.; Sanchez, V.; Zurita, G.; Lozada, M.C.; Cabrera, D. Rolling element bearing defect detection using the generalized synchrosqueezing transform guided by time–frequency ridge enhancement. ISA Trans. 2016, 60, 274–284. [Google Scholar] [CrossRef]
  31. Sapena-Bano, A.; Burriel-Valencia, J.; Pineda-Sanchez, M.; Puche-Panadero, R.; Riera-Guasp, M. The harmonic order tracking analysis method for the fault diagnosis in induction motors under time-varying conditions. IEEE Trans. Energy Convers. 2017, 32, 244–256. [Google Scholar] [CrossRef]
  32. Hu, Y.; Tu, X.; Li, F.; Meng, G. Joint high-order synchrosqueezing transform and multi-taper empirical wavelet transform for fault diagnosis of wind turbine planetary gearbox under nonstationary conditions. Sensors 2018, 18, 150. [Google Scholar] [CrossRef]
  33. Yu, G.; Zhou, Y. General linear chirplet transform. Mech. Syst. Signal Process. 2016, 70, 958–973. [Google Scholar] [CrossRef]
  34. Sheu, Y.L.; Hsu, L.Y.; Chou, P.T.; Wu, H.T. Entropy-based time-varying window width selection for nonlinear-type time–frequency analysis. Int. J. Data Sci. Anal. 2017, 3, 231–245. [Google Scholar] [CrossRef]
  35. Guan, Y.; Liang, M.; Necsulescu, D.S. A velocity synchrosqueezing transform for fault diagnosis of planetary gearboxes under nonstationary conditions. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2017, 231, 2868–2884. [Google Scholar] [CrossRef]
  36. Feng, Z.; Zuo, M.J. Vibration signal models for fault diagnosis of planetary gearboxes. J. Sound Vib. 2012, 331, 4919–4939. [Google Scholar] [CrossRef]
Figure 1. (a,b) are the TFR and extracted instantaneous frequency (IF) obtained using short-time Fourier transform (STFT); (c,d) are the TFR and extracted IF obtained using time-varying demodulation and the STFT (TDSTFT).
Figure 1. (a,b) are the TFR and extracted instantaneous frequency (IF) obtained using short-time Fourier transform (STFT); (c,d) are the TFR and extracted IF obtained using time-varying demodulation and the STFT (TDSTFT).
Sensors 19 03154 g001
Figure 2. TFRs provided using (a) STFT; (b) synchrosqueezed transform (SST); (c) 4-SST; and (d) multi-synchrosqueeezing transform (MSST) ( N = 6 ).
Figure 2. TFRs provided using (a) STFT; (b) synchrosqueezed transform (SST); (c) 4-SST; and (d) multi-synchrosqueeezing transform (MSST) ( N = 6 ).
Sensors 19 03154 g002
Figure 3. Flowchart of the proposed method.
Figure 3. Flowchart of the proposed method.
Sensors 19 03154 g003
Figure 4. The result of a multicomponent simulated signal: (a,b) are the waveform and the Fourier spectrum; (ce) are the TFRs obtained using STFT, MSST, and IMSST, respectively; and (f) is the 2D energy-frequency map.
Figure 4. The result of a multicomponent simulated signal: (a,b) are the waveform and the Fourier spectrum; (ce) are the TFRs obtained using STFT, MSST, and IMSST, respectively; and (f) is the 2D energy-frequency map.
Sensors 19 03154 g004
Figure 5. The result analysis of a multicomponent simulated signal: (a) order spectrum and (b) 2D energy-frequency map.
Figure 5. The result analysis of a multicomponent simulated signal: (a) order spectrum and (b) 2D energy-frequency map.
Sensors 19 03154 g005
Figure 6. The result of simulated signal: (a,b) are the waveform and the Fourier spectrum; (d) is the extracted IF by TDSTFT, where the blue line is the extracted IF and the red line is the real IF; (c,e,f) are the TFRs obtained using STFT, MSST, and IMSST, respectively.
Figure 6. The result of simulated signal: (a,b) are the waveform and the Fourier spectrum; (d) is the extracted IF by TDSTFT, where the blue line is the extracted IF and the red line is the real IF; (c,e,f) are the TFRs obtained using STFT, MSST, and IMSST, respectively.
Sensors 19 03154 g006
Figure 7. (a,b) are the 2D energy-frequency maps by SST and 4-SST based on the same algorithmic theory in Section 2.2 and Section 2.3, respectively.
Figure 7. (a,b) are the 2D energy-frequency maps by SST and 4-SST based on the same algorithmic theory in Section 2.2 and Section 2.3, respectively.
Sensors 19 03154 g007
Figure 8. 2D energy-frequency map by proposed method.
Figure 8. 2D energy-frequency map by proposed method.
Sensors 19 03154 g008
Figure 9. The result analysis of a multicomponent simulated signal: (a) order spectrum and (b) 2D energy-frequency map.
Figure 9. The result analysis of a multicomponent simulated signal: (a) order spectrum and (b) 2D energy-frequency map.
Sensors 19 03154 g009
Figure 10. The result analysis of simulated signals under different noise levels: (a,c,e) are TFRs using IMSST, respectively; (b,d,f) are the 2D energy-frequency maps, respectively.
Figure 10. The result analysis of simulated signals under different noise levels: (a,c,e) are TFRs using IMSST, respectively; (b,d,f) are the 2D energy-frequency maps, respectively.
Sensors 19 03154 g010
Figure 11. The experimental rig.
Figure 11. The experimental rig.
Sensors 19 03154 g011
Figure 12. The structure of experimental rig.
Figure 12. The structure of experimental rig.
Sensors 19 03154 g012
Figure 13. The damage of the sun gear.
Figure 13. The damage of the sun gear.
Sensors 19 03154 g013
Figure 14. Vibration signal acquired from the planetary gearbox. (a) waveform; (b) Fourier spectrum, (c) TFR by STFT, and (d) the extracted IF by TDSTFT.
Figure 14. Vibration signal acquired from the planetary gearbox. (a) waveform; (b) Fourier spectrum, (c) TFR by STFT, and (d) the extracted IF by TDSTFT.
Sensors 19 03154 g014
Figure 15. (a) 2D energy-frequency map of the vibration signal and (b) an enlarged view.
Figure 15. (a) 2D energy-frequency map of the vibration signal and (b) an enlarged view.
Sensors 19 03154 g015
Table 1. Renyi entropy values of TFRs.
Table 1. Renyi entropy values of TFRs.
TFASTFTSST4-SST MSST   ( N = 6 )
Renyi entropy19.116215.351113.591112.0610
Table 2. Parameters of simulated signal.
Table 2. Parameters of simulated signal.
f s ( r ) ( t ) f m ( t ) A 1 A 2 θ 1 θ 2 θ 3 n ( t )
( 10 / 3 ) f s ( r ) ( t ) ( 50 / 3 ) f s ( r ) ( t ) 10.05000−1 dB
Table 3. Configuration of the planetary gearbox.
Table 3. Configuration of the planetary gearbox.
GearSun GearPlanet GearRing Gear
Number of Gear teeth1827 (3)72
Table 4. Fault characteristic frequency of planetary gearbox.
Table 4. Fault characteristic frequency of planetary gearbox.
Gear Meshing FrequencyAbsolute Rotating FrequencyFault Characteristic Frequency
Sun GearPlanet CarrierSun GearPlanet GearRing Gear
f m = 14.4 f s ( r ) f s ( r ) f c = 0.2 f s ( r ) f s s = 2.4 f s ( r ) 0.53 f s ( r ) 0.2 f s ( r )

Share and Cite

MDPI and ACS Style

Lv, Y.; Pan, B.; Yi, C.; Ma, Y. A Novel Fault Feature Recognition Method for Time-Varying Signals and Its Application to Planetary Gearbox Fault Diagnosis under Variable Speed Conditions. Sensors 2019, 19, 3154. https://doi.org/10.3390/s19143154

AMA Style

Lv Y, Pan B, Yi C, Ma Y. A Novel Fault Feature Recognition Method for Time-Varying Signals and Its Application to Planetary Gearbox Fault Diagnosis under Variable Speed Conditions. Sensors. 2019; 19(14):3154. https://doi.org/10.3390/s19143154

Chicago/Turabian Style

Lv, Yong, Bingqi Pan, Cancan Yi, and Yubo Ma. 2019. "A Novel Fault Feature Recognition Method for Time-Varying Signals and Its Application to Planetary Gearbox Fault Diagnosis under Variable Speed Conditions" Sensors 19, no. 14: 3154. https://doi.org/10.3390/s19143154

APA Style

Lv, Y., Pan, B., Yi, C., & Ma, Y. (2019). A Novel Fault Feature Recognition Method for Time-Varying Signals and Its Application to Planetary Gearbox Fault Diagnosis under Variable Speed Conditions. Sensors, 19(14), 3154. https://doi.org/10.3390/s19143154

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