Fault diagnosis of gearbox based on the nonlinear output frequency response functions and PNN

Xingliu Zhang1 , Zili Wang2 , Chen Lu3

1, 2, 3School of Reliability and Systems Engineering, Beihang University, Beijing, 100191, China

1, 2, 3Science and Technology on Reliability and Environmental Engineering Laboratory, Beijing, 100191, China

2, 3Corresponding authors

Vibroengineering PROCEDIA, Vol. 19, 2018, p. 71-75. https://doi.org/10.21595/vp.2018.20170
Received 14 August 2018; accepted 27 August 2018; published 24 September 2018

Copyright © 2018 Xingliu Zhang, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Table of Contents Download PDF Acknowledgements References
Cite this article
Views 48
Reads 24
Downloads 775
CrossRef Citations 0
Abstract.

In this study, the fault diagnosis of nonlinear gearbox is investigated by using the frequency domain approach. The Nonlinear Output Frequency Response Functions (NOFRFs) of the nonlinear gearbox are evaluated and be used to diagnosis the abnormal conditions of the gears and the bearings. A case study is used to illustrate the NOFRFs based feature extraction method, then a probabilistic neural network (PNN) is mentioned to classify the failure modes. Finally, the results show that this method can complete the fault diagnosis of gearbox perfectly. The study provides a basis of the application of the NOFRFs and PNN in the fault diagnosis or condition monitoring of complex nonlinear systems.

Keywords: fault diagnosis, the NOFRFs, PNN, gearbox.

1. Introduction

Gearbox plays an important part in modern industry and the its advantage makes it widely used in automobile, airplane, wind power system [1] etc. In general, the inside of the gearbox is made up of bearings, shafts and gears. Different from simple components, the failure and failure modes are various. Thus, the importance of fault diagnosis for gearbox is obviously.

In general, the frequency domain fault detection approaches are preferred in practice [2]. However, most of the frequency domain methods are conducted based on the linear theory, where the nonlinearities of the system are ignored or linearized.

Recently, a new concept known as the Nonlinear Output Frequency Response Functions (NOFRFs) has been proposed based on the Volterra series representation of a general class of nonlinear systems [3], which is expecting to solve these issues in the fault detection and fault diagnosis of nonlinear rotor-bearing systems. The NOFRFs is the natural extension of the Frequency Response Function (FRF) of linear systems, making it possible to analyze the nonlinear systems in the frequency domain in a similar way to linear systems. The NOFRFs based fault detection and fault diagnosis of nonlinear systems has also been widely studied, i.e. Cao et al. [4] combined the NOFRFs of the nonlinear system and the evidence theory for the fault detection and diagnosis of nonlinear complex systems. Peng et al. [5] used the NOFRFs to detect cracks in beams, showing that the NOFRFs are a sensitive indicator of the presence of cracks under an appropriate input excitation. In Zhao et al. [6], the nonlinear transmissibility derived by the NOFRFs were applied in the condition monitoring of structure damages.

In this study, the NOFRFs concept is applied in the feature extraction of nonlinear gearbox systems, and the results indicate that the higher order NOFRFs are more sensitive than the output frequency response to nonlinear faults which can be used to detect and quantify the nonlinear faults of the system. Then, the PNN is mentioned to classify the failure modes. These two methods work together to complete the fault diagnosis of the gearbox.

2. The NOFRFs based analysis of nonlinear systems

2.1. The frequency domain representation of nonlinear systems

In general, the nonlinear system that stable at zero equilibrium can be represented by using Volterra series as:

(1)
y t = n = 1 N y n t = n = 1 N - - h n τ 1 , , τ n i = 1 n u t - τ i d τ i ,

where yt and ut are the system output and input, respectively; hnτ1,,τn is the nth order Volterra kernel of the nonlinear system.

According to the Volterra series-based representation Eq. (1), in the frequency domain, the output spectrum of the system can be written as:

(2)
Y j ω = n = 1 N Y n j ω , Y n j ω = 1 / n 1 / 2 2 π n - 1 ω 1 + + ω n = ω H n ω 1 , , ω n × i = 1 n U j ω i d σ ω ,

where ω is the output frequency and:

(3)
H n ω 1 ,   , ω n = - - h n τ 1 ,   , τ n × e - j ω 1 τ 1 + + ω n τ n d τ 1 d τ n ,

is the nth order GFRF of the nonlinear system.

The GFRFs Eq. (3) of the nonlinear system can be evaluated by using the recursive algorithm. However, it can be seen that the GFRFs Eq. (3) is a multi-dimensional function which is difficult to compute and analysis, especially for higher order cases. To solve this problem, a new concept known as the NOFRFs of the nonlinear system was proposed [3], which is introduced as below.

2.2. The NOFRFs of nonlinear systems

The NOFRFs of the nonlinear system is defined as:

(4)
G n j ω = Y n j ω U n j ω ,         n = 1 , , N ,

where Ynjω is the nth order output spectrum of the system, and:

(5)
U n j ω = 1 / n 1 / 2 2 π n - 1 ω 1 + + ω n = ω i = 1 n U j ω i d σ ω = F F T u n t .

Is the nth input, such that the output spectrum Eq. (2) can be rewritten as:

(6)
Y j ω = n = 1 N Y n j ω = n = 1 N G n j ω U n j ω .

The new representation Eq. (6) indicate that the nonlinear system which can be represented by a Volterra series can be analyzed, similar to the FRF method of linear systems, by suing the NOFRFs. The approach is illustrated in Fig. 1.

It is worth noticing that, when the input of the nonlinear system is harmonic as:

(7)
u t = A c o s ω F t + β ,

where A is the magnitude of the input signal, ωF and β are the input frequency and the phase, respectively, the NOFRFs based representation of the system can be written as:

(8)
Y j ω F = G 1 j ω F U j ω F + G 3 j ω F U 3 j ω F + G 5 j ω F U 5 j ω F + ,
Y j 2 ω F = G 2 j 2 ω F U j 2 ω F + G 4 j 2 ω F U 4 j 2 ω F + G 6 j 2 ω F U 6 j 2 ω F + ,

where the NOFRFs under a harmonic input is independent to the input magnitude A.

In practice, the NOFRFs of the nonlinear system can be estimated by using the LS, where the truncation order N is assumed to be sufficient and a set of input with magnitudes are required. In this study, the NOFRFs under the harmonic input is applied for the feature extraction of the nonlinear gearbox system instead of the traditional frequency response detection method.

Fig. 1. The NOFRFs based analysis of a nonlinear system

 The NOFRFs based analysis of a nonlinear system

3. Fault classification based on PNN classifier

The PNN, first proposed by Dr. D. F. Speeht in 1989, is a branch of the radial basis Network, belonging to the feedforward Network. It has the following advantages: simple learning process, fast training speed; Classification is more accurate, fault tolerance and so on. In essence, it belongs to a supervised network classifier based on Bayesian minimum risk criteria.

Fig. 2. Basic structure of PNN

 Basic structure of PNN

PNN generally have the following four layers: input layer, mode layer, summation layer, and output layer. In some materials, the pattern layer is also called the hidden layer, and the summation layer is called the competition layer. The input layer is responsible for transmitting the feature vector to the network, and the number of input layers is the number of sample features. The mode layer is connected to the input layer by a connection weight. Calculate the degree of matching between the input feature vector and each mode in the training set, that is, the similarity, and send the distance to the Gaussian function to obtain the output of the mode layer. The number of neurons in the pattern layer is the number of input sample vectors, that is, how many samples there are, and how many neurons there are. The summation layer is responsible for connecting the pattern layer units of each class. The number of neurons in this layer is the number of categories of samples. The output layer is responsible for outputting the highest score in the summation layer.

4. Case study

The data sets used in this paper come from the 2009 PHM Conference Data Analysis Competition.

Table 1. Failure modes chosen of the samples

Case
Gear
Bearing
32T
96T
48T
80T
IS: IS
ID: IS
OS: IS
IS: OS
ID: OS
OS: OS
1
Intact
Intact
Intact
Intact
Intact
Intact
Intact
Intact
Intact
Intact
2
Chipped
Intact
Eccentric
Intact
Intact
Intact
Intact
Intact
Intact
Intact
3
Intact
Intact
Eccentric
Intact
Intact
Intact
Intact
Intact
Intact
Intact
4
Intact
Intact
Eccentric
Broken
Ball
Intact
Intact
Intact
Intact
Intact
5
Chipped
Intact
Eccentric
Broken
Inner
Ball
Outer
Intact
Intact
Intact
6
Intact
Intact
Intact
Broken
Inner
Ball
Outer
Intact
Intact
Intact

4.1. Evaluation of the NOFRFs and Feature extraction of the vibration signals

The NOFRFs representation Eq. (6) of the system output spectrum under the new input αut can be rewritten as:

(9)
Y j ω = α G 1 j ω U j ω + α 2 G 2 j ω U 2 j ω + + α N G N j ω U N j ω .

Such that the NOFRFs G1jω, G2jω, etc. can be evaluated by the followings.

Step 1: Test the system by using a set of different input signals αiut, i= 1, …, N- with N-N;

Step 2: Compute the output spectrum of the system output under the N- test inputs, denote as Yijω, i= 1, …, N-;

Step 3: Rewritten the NOFRFs representation into a matrix form as:

(10)
Y j ω = A U j ω G j ω .

According to the N- tests results, where:

Y j ω = Y 1 j ω ,   ,   Y N - j ω T ,
G j ω = G 1 j ω ,   ,   G N j ω T ,
A U j ω = α 1 U j ω α 1 N U N j ω     α N - U j ω α N - N U N j ω .

Step 4: Evaluate the NOFRFs by using the LS method as:

(11)
G j ω = A U j ω T A U j ω - 1 A U j ω T Y j ω .

In this study, the truncation order of the NOFRFs representation is assumed to be N= 5, and five different inputs are applied to evaluated the NOFRFs.

4.2. Fault classification based on PNN classifier

After completing the work of NOFRF, the PNN model needs to be trained to classify the various patterns. First, the first 200 groups of each pattern were used as training data, and the last 100 groups were used as test data. Then the PNN classifier is obtained by using the training data. Finally, the trained PNN classifier is used to identify the test data and obtain the final fault diagnosis results. The results of the test data are shown in Table 2, indicating the excellent classification ability of PNN.

As is shown in Table 2, the method based on the NOFRFs and PNN has a good diagnostic result that the accuracies are not lower than 0.90 even when the gearbox under different cases.

Table 2. Classification results based on the NOFRFs and PNN

Variety
Load
30 Hz
35 Hz
40 Hz
45 Hz
50 Hz
Input side data
High1
1.00
1.00
0.90
1.00
0.98
High2
1.00
1.00
1.00
1.00
1.00
Low1
1.00
0.95
0.95
1.00
1.00
Low2
1.00
1.00
1.00
0.97
1.00
Output side data
High1
1.00
0.91
0.91
1.00
1.00
High2
1.00
1.00
1.00
1.00
1.00
Low1
1.00
0.97
0.96
1.00
0.97
Low2
1.00
0.99
1.00
1.00
1.00

5. Conclusions

The NOFRFs of nonlinear systems are the natural extension of the FRF for linear systems, which can be used for nonlinear systems analyses. In this study, the feature extraction of nonlinear gearbox system was completed by using NOFRFs concept, and then the fault mode was classified by using PNN classifier. The diagnostic results are obtained by study the vibration signals of gearbox. The results show that the way presented in this study can diagnose the gearbox failure perfectly.

Acknowledgements

This research was supported by the National Natural Science Foundation of China (Grant Nos. 51605014, 51105019 and 51575021), the Aviation Science Fund (Grant No. 20163351018), the Technology Foundation Program of National Defense (Grant No. Z132013B002), and the Fundamental Research Funds for the Central Universities (Grant No. YWF-18-BJ-Y-159).

References

  1. Ewins D. J. Control of vibration and resonance in aero engines and rotating machinery – an overview. International Journal of Pressure Vessels and Piping, Vol. 87, Issue 9, 2010, p. 504-510. [Publisher]
  2. Lei Y., Lin J., He Z., Zuo M. J. A review on empirical mode decomposition in fault diagnosis of rotating machinery. Mechanical Systems and Signal Processing, Vol. 35, Issue 1, 2013, p. 108-126. [Publisher]
  3. Lang Z. Q., Billings S. A. Energy transfer properties of non-linear systems in the frequency domain. International Journal of Control, Vol. 78, Issue 5, 2005, p. 345-362. [Publisher]
  4. Cao J., Chen L., Zhang J., Cao W. Fault diagnosis of complex system based on nonlinear frequency spectrum fusion. Measurement, Vol. 46, Issue 1, 2013, p. 125-131. [Publisher]
  5. Peng Z. K., Lang Z. Q., Billings S. A. Crack detection using nonlinear output frequency response functions. Journal of Sound and Vibration, Vol. 301, Issue 3, 2007, p. 777-788. [Publisher]
  6. Zhao X. Y., Lang Z. Q., Park G., Farrar C. R., Todd M. D., Mao Z., Worden K. A new transmissibility analysis method for detection and location of damage via nonlinear features in MDOF structural systems. IEEE/ASME Transactions on Mechatronics, Vol. 20, Issue 4, 2015, p. 1933-1947. [Publisher]