# Gear Fault Diagnosis in Variable Speed Condition Based on Multiscale Chirplet Path Pursuit and Linear Canonical Transform

Authors: Xu Shuiqing; Zhang Ke; Chai Yi; He Yigang; Feng Li
Journal: Complexity
 (2018)

---

## Abstract

The vibration signals analysis is a very effective and reliable method for detecting the gear failures. Because the vibration signals acquired from the gear in the variable speed condition often contain more useful fault information, the analysis of the gear vibration signals during the variable speed condition has been a hot research topic. In this paper, a method based on the multiscale chirplet path pursuit (MSCPP) and the linear canonical transform (LCT) has been applied to diagnose the gear fault in the variable speed condition for the first time. First, by using the MSCPP method to estimate the instantaneous meshing frequency, the suitable signal segment approximation to the acceleration or deceleration process can be selected. Then, because the LCT is a novel and efficient nonstationary signals analysis tool, the optimal LCT spectrum of the selected signal has been attainted to diagnose the gear faults based on the properties of the LCT. In addition, the simulations and the experimental evaluation are provided to verify the effectiveness of the proposed method.


---

## Body

## 1. Introduction

The gear is an important device and has a wide range of applications in industry. However, owing to abrasion and other reasons, the gears may have many kinds of faults, such as pitting, chipping, and the serious crack [1–5]. When the gear causes a local fault, the amplitude and phase of the vibration signal of the gear are modulated [2, 3, 6, 7]. So the vibration signals obtained from the gear can reflect the gear’s state very well. The vibration signals analysis also has been the extremely effective and reliable method for detecting the gear failures [8–11]. Meanwhile, on the one hand, the gear vibration signals obtained from the variable speed condition, that is, the acceleration or deceleration processes, often can contain more fault information compared to the stationary processes, which can detect the gear faults earlier. On the other hand, the gear vibration signals obtained from the acceleration or deceleration processes are nonstationary signals, which have low signal-to-noise ratio (SNR) in practice [11–14], which make it difficult to obtain the gear vibration signals’ fault features. Hence, the diagnosis of gear faults by analyzing the gear vibration signals obtained from the variable speed condition has been a hot research topic recently.A number of failure diagnosis methods have been used to diagnose the gear faults in variable speed condition, for example, the traditional time frequency analysis methods, self-adaptive signal processing methods, and data driven methods [4–12, 15–25]. The traditional time frequency analysis methods, such as the short time Fourier transform [14], the Wigner-Ville distribution [10, 12], and the wavelet transform [2], will result in spectral aliasing, cross term interference, and low resolution because the vibration signals obtained from the fault gears in practice are nonstationary with low SNR [4]. The self-adaptive signal processing methods, for example, the empirical mode decomposition [5] and the local mean decomposition [17], will lead to over envelope, end effects, and distorted components, respectively. In addition, the data driven methods, such as the Elman neural network [24] and support vector machine [22], require a lot of data for training and classification. To solve these problems, some other data driven methods have been introduced, such as the self-organizing maps (SOM) method [26–28]. In [26], the self-organizing feature map neural network has been used to diagnose the faults of the wind turbine’s converter. In [27], the SOM and minimum quantization error (MQE) method have been selected to achieve degradation assessment and fault localization. In [28], an intelligent approach based on the SOM method has been presented to machine component health prognostics. Although the SOM method does not need a lot of data, it still needs some a priori knowledge. Therefore, the diagnosis of gear faults in acceleration and deceleration processes remains an open research field and new signal analysis tools for nonstationary vibration signals are needed.The linear canonical transform (LCT) is a generalization of the Fourier transform (FT) and the fractional Fourier transform (FRFT), which has four-parameter family of linear integral transform [29, 30]. It performs an affine mapping of the time frequency distribution of the signal because it has additional degrees of freedom [31, 32]. So, the LCT is more flexible and suitable for processing nonstationary signals, particularly in the linear frequency modulated (LFM) signals [33–37]. Simultaneously, the gear vibration signals obtained from the variable speed condition can be approximated to the multicomponent LFM signals within a short period of time. From the above analysis, because the vibration signals attained from the variable speed condition are nonstationary and the LCT has advantages in nonstationary signals processing, it is therefore worthwhile to explore the diagnosis of gear faults by using the LCT.However, in variable speed condition, the gear vibration signals contain many different components, which are not approximated to the LFM signals. The vibration signals acquired from the acceleration or deceleration processes of the gears are often a short period of the gear vibration signals in variable speed condition. In order to apply the LCT method to diagnosis of the gear faults in variable speed condition, the vibration signals acquired from the acceleration or deceleration processes should be selected at first. In practice, it is difficult to only obtain the gear vibration signals during acceleration or deceleration processes directly. Nevertheless, the shaft rotational frequency (SRF) is time varying in the acceleration or deceleration processes, which can be seen as an indicator to select the acceleration or deceleration processes [1]. Hence, we can obtain the acceleration or deceleration time by estimating the instantaneous frequency (IF) of the gear vibration signals in variable speed condition. In addition, the optimal LCT parameters of the LFM signals are also determined by gradient of the IF [30]. For these reasons, we need to estimate the IF of the gear vibration signals in variable speed condition at first.In this paper, the multiscale chirplet path pursuit (MSCPP) has been used in the estimation of the IF of the gear vibration signals, which is a widely used and efficient method for IF estimation [38]. The MSCPP method divides the time length of the analysis signal into a series of dynamic time support areas in binary form, finding atom with the largest correlation in each dynamic support region from the defined chirplet atom library, and the selected chirplet atoms are connected one by one by using the principle of the best connection [38, 39]. Thus, it can adaptively find the instantaneous frequency trend line with the largest correlation in the signal [39]. Compared with other instantaneous frequency methods, the MSCPP method can choose shorter chirplet atoms flexibly and can effectively suppress noise interference. In addition, the MSCPP method does not have any a priori knowledge. So the instantaneous frequency of the vibration signals obtained from the variable speed condition can be estimated by the MSCPP method. Based on this, the vibration signals obtained from the acceleration or deceleration processes in a short period of time can be picked out.In this paper, the MSCPP method and the LCT have been applied to analyze the vibration signals to diagnose the gear faults in variable speed condition for the first time. First, by using the MSCPP, the IF of the vibration signals obtained from the variable speed condition can be estimated. Then, based on the IF of the vibration signals, we can choose the suitable signal segment, which can be seen as the acceleration or deceleration processes. In addition, according to the gradient of the instantaneous frequency of the selected signal segment, the optimal LCT parameters can be obtained based on the properties of the LCT [30]. Finally, the optimal LCT spectrum of the selected signal is presented, which can be used to diagnose the gear fault. The block diagram of the proposed method is presented in Figure 1 vividly.Figure 1
The block diagram of the proposed method.The remaining parts of this article are organized as follows. In Section2, the preliminaries of the MSCPP and the LCT are introduced. In Section 3, the application of the MSCPP and the LCT to the simulated vibration signal is presented, which shows that the proposed method is effective. The experimental evaluations have been provided in Section 4. Section 5 concludes this paper.


## 2. Preliminaries

### 2.1. The Multiscale Chirplet Path Pursuit Method

The MSCPP method was first introduced in [38]; in this method, if the IF of chirplet atoms is linear, they can be used to adaptively piecewise approximate the nonstationary signals. For example, the nonstationary signal in the following is considered.(1)yt=Acos⁡wt+φ+nt,where A, φ, n(t), and w(t) represent the amplitude, the initial phase, the noise, and the successive and derivable instantaneous phase of the signal, respectively. In this algorithm, the chirplet atoms are obtained from the chirplet dictionary that can be written as a set of functions as follows [38].(2)fm,n,It=I-1/2eimt2/2+ntLIt,where m is the slope coefficient and n is the offset coefficient. According to (2), we can easily have that the IF of a chirplet is linear and equal to mt+n. Based on the sampling theorem, we know that mt+n should be less than a half of the sampling rate fs. I is the dyadic time span, which has the binary scale division and can be defined as I={Ij,k∣j=0,1,…,log2N-1,k=0,1,…,2j-1}. In this equation, Ij,k=[k2-jT,(k+1)2-jT]. T and N are the total sampling time and the number of samples, respectively. j is the scale coefficients and j=0,1,…,log2N-1, k=0,1,…,2j-1. LI(t) is a rectangular window function, which is 1 when t∈I and 0 when t∉I. I-1/2 is the normalization factor that makes fm,n,I(t)2=1.In each time interval, by computing the maximum projection coefficientCIj,k, where CIj,k equals CIj,k=max⁡x,fm,n,Ij,k(t) (· stands for the interior product) and includes the amplitude and initial phase information of the signal component in the time interval Ij,k, the chirplet atom which has the highest correlation to signal y(t) can be obtained from the chirplet dictionary. Suppose SIj,k is the signal component expressed by CIj,k in the time interval Ij,k, then SIj,k can be written in the following [38]:(3)SIj,kt=CIj,ke-imt2+nt+∠CIj,kIj,k.


### 2.2. The Linear Canonical Transform

The LCT of a signalf(t) with parameter A=(a,b,c,d), represented as LfA(u), is defined by [29](4)
L
f
A
u
=
∫
-
∞
∞
f
t
K
A
u
,
t
d
t
b
≠
0
d
e
j
c
d
u
2
/
2
f
d
u
b
=
0
,where a,b,c,d are real numbers satisfying ad-bc=1. In addition, the kernel KA(u,t) is given by(5)KAu,t=CAeja/2bt2-ut/b+d/2bu2and CA=1/j2πb. According to (4), it can easily be obtained that the LCT has three degrees of freedom and when the LCT parameters are certain, the kernel KA(u,t) of the LCT is a LFM signal; thus the LCT of a LFM signal with suitable parameters can be a dirac function [31]. Therefore, the LCT is more flexible and suitable for processing nonstationary signals, especially for the LFM signals.


## 2.1. The Multiscale Chirplet Path Pursuit Method

The MSCPP method was first introduced in [38]; in this method, if the IF of chirplet atoms is linear, they can be used to adaptively piecewise approximate the nonstationary signals. For example, the nonstationary signal in the following is considered.(1)yt=Acos⁡wt+φ+nt,where A, φ, n(t), and w(t) represent the amplitude, the initial phase, the noise, and the successive and derivable instantaneous phase of the signal, respectively. In this algorithm, the chirplet atoms are obtained from the chirplet dictionary that can be written as a set of functions as follows [38].(2)fm,n,It=I-1/2eimt2/2+ntLIt,where m is the slope coefficient and n is the offset coefficient. According to (2), we can easily have that the IF of a chirplet is linear and equal to mt+n. Based on the sampling theorem, we know that mt+n should be less than a half of the sampling rate fs. I is the dyadic time span, which has the binary scale division and can be defined as I={Ij,k∣j=0,1,…,log2N-1,k=0,1,…,2j-1}. In this equation, Ij,k=[k2-jT,(k+1)2-jT]. T and N are the total sampling time and the number of samples, respectively. j is the scale coefficients and j=0,1,…,log2N-1, k=0,1,…,2j-1. LI(t) is a rectangular window function, which is 1 when t∈I and 0 when t∉I. I-1/2 is the normalization factor that makes fm,n,I(t)2=1.In each time interval, by computing the maximum projection coefficientCIj,k, where CIj,k equals CIj,k=max⁡x,fm,n,Ij,k(t) (· stands for the interior product) and includes the amplitude and initial phase information of the signal component in the time interval Ij,k, the chirplet atom which has the highest correlation to signal y(t) can be obtained from the chirplet dictionary. Suppose SIj,k is the signal component expressed by CIj,k in the time interval Ij,k, then SIj,k can be written in the following [38]:(3)SIj,kt=CIj,ke-imt2+nt+∠CIj,kIj,k.


## 2.2. The Linear Canonical Transform

The LCT of a signalf(t) with parameter A=(a,b,c,d), represented as LfA(u), is defined by [29](4)
L
f
A
u
=
∫
-
∞
∞
f
t
K
A
u
,
t
d
t
b
≠
0
d
e
j
c
d
u
2
/
2
f
d
u
b
=
0
,where a,b,c,d are real numbers satisfying ad-bc=1. In addition, the kernel KA(u,t) is given by(5)KAu,t=CAeja/2bt2-ut/b+d/2bu2and CA=1/j2πb. According to (4), it can easily be obtained that the LCT has three degrees of freedom and when the LCT parameters are certain, the kernel KA(u,t) of the LCT is a LFM signal; thus the LCT of a LFM signal with suitable parameters can be a dirac function [31]. Therefore, the LCT is more flexible and suitable for processing nonstationary signals, especially for the LFM signals.


## 3. The Proposed Method of the Simulated Vibration Signal

When the gear causes a partial failure, the amplitude and phase of the vibration signal of the gear are modulated, which are periodic with the gear’s rotation frequency [1]. The vibration signal x(t) attained from a pair of meshing gears with tooth fault could be represented as [1](6)xt=∑m=0MXm1+amtcos⁡2πmkfst+ϕm=∑m=0MXm1+∑n=0NAmncos⁡2πnfst+σmn×cos⁡2πmkfst+ϕm,where am(t) is the amplitude modulating function and M represents the amount of tooth meshing harmonics. Xm and ϕm indicate the mth meshing harmonic’s amplitude and phase, separately [1]. kfs is the meshing frequency, k expresses the gear teeth’s amount, and fs stands for the SRF [1]. N represents sidebands’ amount around the tooth meshing harmonics. Amn and σmn are the amplitudes and phases at the nth sidebands of the amplitude modulated signals around the mth meshing harmonic, separately [2, 4].The SRFfs in (6) is time varying and the vibration signal can be approximated to the LFM signal when the instantaneous rotational speed is ascending and descending straight in the acceleration and deceleration, respectively. To verify the proposed algorithm is effective for the vibration signals; next, a signal is considered to simulate the vibration signal obtained from the acceleration process as follows based on (6):(7)xt=1+cos⁡2πfmtcos⁡2π·10·fmt+φt,where fm stands for the SRF in the simulated vibration signal, which is chosen as fm=2t+30 when 0≤t<1.5 and equal to 33 when 1.5≤t≤2. The sampling frequency is 2000 Hz and the number of the sample points is set as 4096. The time interval is set as 0 s to 2 s. Moreover, because the vibration signals attained from the gear in variable speed condition are usually with low SNR, φ(t) is selected as the Gaussian white noise with the SNR is −4 dB. Based on the above, Figure 2 depicts the signal shown in (7) without the noise being depicted in it. The simulated vibration signal based on (7) is presented in Figure 3. Figure 4 shows the FT spectrum of the simulated vibration signal. According to Figure 3, because the signal contains noise, any useful information can not be obtained directly. Analogously, since the simulated vibration signal is nonstationary and has noise, this leads to the FT spectrum in Figure 4 that has spectral aliasing and blurring.Figure 2
The simulated gear vibration signal without noise.Figure 3
The simulated gear vibration signal with noise.Figure 4
The FT of the simulated gear vibration signal with noise.Now, to obtain the fault features from the simulated gear vibration signal and diagnose the gear faults, the MSCPP and the LCT are applied. According to the above analysis, we first use the MSCPP method to estimate the instantaneous frequency of the simulated gear vibration signal. Because the sampling frequency is 2000 Hz, we let the search range and the search step lengths of the frequency slope be −100 to 100 Hz/s and 1 Hz/s, respectively. The search range and the search step lengths of offset coefficients are 0 to 100 Hz and 1 Hz, separately. The points of any dyadic time span is 32. Therefore, we can obtain thatjmax=log2(N/32), which is equal to 7. Based on these, by applying the MSCPP algorithm, the IF of the simulated gear vibration signal is obtained in Figure 5. In Figure 5, the blue line and the red line represent the actual and estimated instantaneous frequency of the simulated gear vibration signal, respectively. The red line is coincided with the blue line very well in Figure 5, which shows that the MSCPP algorithm can be used to estimate IF of the simulated gear vibration signal efficiently. Based on Figure 6, we can know that the time interval from 1 s to 1.5 s can be selected to approach the acceleration process of the gear in a short period of time. Next, based on the gradient of the signal segment from 1 s to 1.5 s and the properties of the LCT [30], the optimal LCT spectrum with parameters (0.93,0.42,0.34,1.89) of the selected signal segment is obtained in Figure 6. It is shown that Figure 6 has three peaks, which are (595,1066), (661,1915), and (727,869.6), respectively. The highest peak is (661,1915). The LCT frequency distances of each of the two adjacent peaks are 66 Hz. It is shown that the frequency of (661,1915) is almost ten times of the SRF associated with the LCT. This means that the peaks (595,1066), (661,1915), and (727,869.6) are 9, 10, and 11 times of the SRF associated with the LCT, separately. From the results, we can know that it is very consistent with the analysis of (7).Figure 5
Estimated instantaneous frequency of the simulated gear vibration signal.Figure 6
The LCT of the simulated gear vibration signal with noise.Moreover, since the results obtained in Figure6 show that the range of any two adjacent peaks are equal, hence the signal presented in Figure 3 could be considered as the fault gear vibration signal based on the description in [1, 4]. Thus, the proposed method is an effective method to diagnose the gear fault according to this simulation.


## 4. Experimental Evaluation

In Section3, the validity of the MSCPP and the LCT in gear faults diagnosis has been verified by the signal presented in (7). In order to further verify the correctness of the proposed algorithm, in this part, the gear faults of experiment setup presented in Figure 7 have been used. In this gear experiment equipment, the output gear has 75 teeth and the input pinion gear has 55 teeth. The sample frequency and the sample time of the signals obtained from this experiment equipment are 4000 Hz and 2 s, respectively.Figure 7
Platform of experiment.Firstly, the vibration signal attained from a normal gear in variable speed condition is presented in Figure8. The signal presented in Figure 8 is nonstationary since the gear signals with low SNR are obtained from variable speed condition. On this account, Figure 9 depicts the FT of the signal shown in Figure 8 which has spectrum aliasing, and it cannot be used to diagnose the gear fault. Next, by applying the method we presented in this paper, we used the MSCPP method to estimate the instantaneous frequency of the gear vibration signal at first. In the MSCPP method, the seeking scope of frequency offset coefficients and slope coefficient are [0, 2000 Hz] and [−100, 100 Hz/s], respectively. The search resolutions of these are chosen as 1 Hz/s and 1 Hz, separately. The samples are 4096. Based on these, we can obtain the maximum scales of 7. Thus, the IF of the gear vibration signal is obtained in Figure 10. From Figure 10, it is easy to know that the time interval from 1 s to 1.4 s can be selected to approach the acceleration process of the gear in a short period of time. Subsequently, by utilizing the advantages of the LCT and the properties of the selected signal segment, the optimal LCT spectrum with parameters (-0.17,1.53,-0.98,2.93) of the selected signal segment is shown in Figure 11. We can easily know that Figure 11 has only one peak. That is to say, the signal presented in Figure 8 has not amplitude modulated. Hence, the gear can be seen as a normal gear.Figure 8
The vibration signal obtained from the normal gear.Figure 9
The FT of the vibration signal obtained from the normal gear.Figure 10
Estimated instantaneous frequency of the signal presented in Figure13.Figure 11
The LCT of the selected signal segment.Then, Figure12 shows a gear vibration signal, which is attained from the gear which has a broken tooth in variable speed condition. The FT of the signal presented in Figure 12 is presented in Figure 13, which shows that it also has spectrum aliasing and cannot be used to diagnose the gear fault. Now, by applying the proposed method, the search range of frequency offset coefficients and slope coefficient are also set as [0, 2000 Hz] and [−100, 100 Hz/s]. The search resolutions of these are also chosen as 1 Hz/s and 1 Hz, separately. The samples are 4096. Based on these, the estimated instantaneous frequency of the gear vibration signal can been obtained in Figure 14 based on the MSCPP method. From Figure 14, the time interval from 1 s to 1.5 s can be selected to approach the acceleration process of the gear in a short period of time. Subsequently, according to the above analysis, the optimal LCT spectrum with parameters (-0.51,0.87,-1.47,0.56) of the selected gear vibration signal segment can be obtained in Figure 15. We can know that Figure 15 has three peaks, that is, (613,2403), (636,3512), and (659,2524). Since (636,3512) is the highest and there are 55 teeth in the pinion gear, thus the input SRF associated with the LCT can be expressed as 636 Hz/55 = 11.56 Hz [4], which is equal to half of the frequency distance. That is to say, the meshing frequency is twofold of the SRF. Based on the above analysis, the gear should be considered to have fault happened.Figure 12
The vibration signal obtained from the fault gear.Figure 13
The FT of the vibration signal obtained from the fault gear.Figure 14
Estimated instantaneous frequency of the signal presented in Figure17.Figure 15
The LCT of the selected signal segment from the fault gear.Figure 16
Another vibration signal obtained from the fault gear.Figure 17
The FT of the vibration signal presented in Figure16.In addition, in order to further verify the correctness of the proposed algorithm, Figure16 presents another gear vibration signal, which is attained from the fault gear in different variable speed condition. The FT of the signal presented in Figure 16 is presented in Figure 17, which also shows that it cannot be used to diagnose the gear fault. Now, by applying the proposed method, similar to the above cases, the largest scale is also 7. Then, the estimated IF of the gear vibration signal has been obtained in Figure 18 based on the MSCPP method. From Figure 18, the time interval from 1 s to 1.4 s can be selected to approach the deceleration process of the gear in a short period of time. Then, the optimal LCT spectrum with parameters (1.53,0.57,0.42,0.81) of the selected gear vibration signal segment can be obtained in Figure 19. It is shown that Figure 15 also has three peaks (591,2484), (613,3543), and (635,2103). Since (613,3543) is the highest and there are 55 teeth in the pinion gear, thus the input SRF associated with the LCT can be expressed as 613 Hz/55 = 11.14 Hz. That is to say, the meshing frequency is twofold of the SRF in the LCT domain. Thus, it can be said that the gear also have fault that occurred.Figure 18
Estimated instantaneous frequency of the signal presented in Figure16.Figure 19
The LCT of the selected signal segment presented in Figure16.Based on the normal case and two fault cases, it is shown that the method proposed in this can diagnose the gear faults in variable speed condition. However, the MSCPP and LCT method only can show that the gear faults happened, and our research directions will be the diagnosis of the types and the severity of the gear faults.


## 5. Conclusion

In this paper, a method based on the MSCPP and the LCT has been applied to diagnose the gear faults in the variable speed condition for the first time. Firstly, the preliminaries of the MSCPP and the LCT have been presented. Then, the proposed method of the simulated gear vibration signals have been showed. At last, in order to further verify the correctness of the proposed algorithm, the diagnosis of actual gear vibration signals also has been presented. It is indicated that the proposed method can diagnose the gear faults availably. In the future, the diagnosis of early gear faults, intermittent gear faults, and multiple gear faults also will be our future research directions.



---

*Source: 3904598-2018-03-11.xml*