Engine fault diagnosis based on multiresolution analysis and dual microphones
Rong Zeng^{1} , Ruili Zeng^{2} , Jianmin Mei^{3} , Lei Ding^{4}
^{1, 4}Postgraduate Training Brigade, Army Military Transportation University, Tianjin, China
^{2, 3}Military Vehicle Department, Army Military Transportation University, Tianjin, China
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 18, 2018, p. 4651.
https://doi.org/10.21595/vp.2018.19931
Received 21 April 2018; accepted 28 April 2018; published 22 May 2018
JVE Conferences
In order to diagnose the abnormal sound fault of the diesel without disassemble, dual microphones are used to collect noise signal of the engine in this paper, and the sound pressure signal is decomposed with wavelet decomposition technique. By comparing different wavelet decomposition techniques, it determines the sound feature parameter for fault diagnosis, and establishes a BP neural network for pattern recognition. The test result shows that diagnosing engine fault with dual microphones is accuracy, and it can be applied to diagnose the abnormal sound fault of engine.
Keywords: diesel engine fault diagnosis, dual microphone, multiresolution analysis, BP neural network.
1. Introduction
Because of the low signaltonoise ratio, the feature of acoustic signals is more difficult to be extracted. For the problem of low signaltonoise ratio, many researchers have carried out some research, such as using wavelet analysis, applying wavelet transform to the processing of acoustic signals for rolling bearing failure [1], and improving the signaltonoise ratio of acoustic signals through multiscale decomposition [2]; using independent component analysis and autocorrelation analysis of noise reduction [3]; the use of morphological component analysis method for generalized soft threshold decrease noising [4]; and the use of MEL frequency inverse spectral coefficient algorithm to extract fault features [5]. These noise reduction methods achieve better results under the condition that the sound source environment is single.
However, these existing studies are only the types of faults with obvious diagnostic characteristics. In practice, most of the engine faults are gradual processes from light to heavy (for example, slight oil leakage to complete oil shutoff), while it's useful that minor faults are diagnosed.
For acoustic signals, analyzing the frequency domain is more effective than the time domain because the more important characteristics are mainly reflected in the power spectrum, and the spectrum is not easily affected by external interference, and the acoustic characteristics are relatively stable.
This article mainly analyzes the low frequency signal of the engine. Based on multiresolution analysis and band energy analysis, combined with the strong ability of neural network in nonlinear mapping processing, this paper presents a method of engine fault diagnosis using dualmicrophone measurement combined with wavelet multiresolution analysis.
2. Signal processing and pattern recognition method
2.1. Wavelet multiresolution analysis
According to the characteristics of engine noise signal, this paper uses wavelet multiresolution analysis method to deal with engine noise signal. The timefrequency resolution of the wavelet transform varies with the time decomposed scale. For the lowfrequency component of the signal, the time resolution is higher, and the frequency resolution is lower. For the highfrequency component of the signal, the time resolution is lower, and the frequency resolution is higher. S. Mallat proposed the concept of multiresolution analysis when constructing orthogonal wavelet bases, and was first applied to wavelet analysis. Wavelet analysis is also mature in the field of fault diagnosis at the present stage [68]. The multiresolution analysis wavelet decomposition tree is shown in Fig. 1.
Fig. 1. Fourlayer multiresolution analysis tree structure
($X$ is the signal to be analyzed; L1L4 is the lowfrequency part of 14 layers in sequence; H1H4 is the highfrequency part of 14 layers in sequence)
If the highest frequency is $f$, the frequency range for the highfrequency H1 signal is [$f$/2, $f$], and the lowfrequency L1 frequency range is [0, $f$/2]. When the second layer is decomposed, the high frequency portion H1 obtained by the decomposition of the original first layer is not decomposed, only the low frequency portion L1 is decomposed, and so on, the signal is decomposed layer by layer. After the signal $X$ is decomposed by the fourlayer multiresolution shown in the figure, the signal coefficients ($i=$ 1, 2, 3, 4) can be obtained, and the corresponding frequency ranges are [0, $f$/16), [$f$/16, $f$/8), [$f$/8, f/4), [$f$/4, $f$/2), [$f$/2, $f$]. That is, signal $X$ has the following equation:
Since the energy of the timedomain calculation signal is equal to the energy value calculated in the frequency domain, each band energy is ${E}_{i}$, and the feature vector is constructed:
If $E=\sum {E}_{i}$, then normalized eigenvector is:
Fig. 2. 4layer wavelet decomposition results
a) Engine duty cycle signal
b) After 4layert wavelet decomposition
As shown in Fig. 2, the amplitude of the frequency band coefficient of the acoustic signal of the engine after multiresolution analysis is obviously different, and five waveform segments with different amplitudes can be clearly seen.
2.2. Band energy accumulation feature extraction
According to the wavelet decomposition results, characteristic parameters are obtained by calculating the band energy accumulation. The waveletband energy accumulating method mainly analyzes the timefrequency transform of the signal and decomposes it into different frequency bands, performs statistical analysis on the total energy in each frequency band, judges the change of the fault signal with the change of the energy in the band, and uses these characteristic changes to establish the signal. The mapping relationship between frequency band energy and engine abnormal noise status.
The preliminary analysis of the acoustic signal of the engine shows that the Daubechies family of wavelet basis functions db4 or db5 and the number of decomposition layers 3, 4, and 5 (energy bands 4, 5, and 6) can achieve better results. Comprehensive analysis of normal conditions, normal oil leakage, slight oil leakage, moderate oil leakage, and complete oil breakage (misfire failure), as shown in Fig. 4. The length shown in each band is the proportion of the corresponding band to the total band energy. The abscissa is the percentage of the total energy in each frequency band, and the ordinate is the four fault conditions after wavelet decomposition. From the figure, it can be found that when the number of decomposition layers is 3, the energy characteristics of the frequency band of slight oil leakage failure and moderate oil leakage failure are completely indistinguishable, and when the number of decomposition layers is 5, the stability of the band energy characteristics is poor. From the comparison of (a)(f) in Fig. 4, when the number of decomposition layers is 4 and the wavelet basis function is db4, the fault condition is better differentiated and the stability of the same state feature is better. In this paper, the fault eigenvectors are established by the wavelet basis function db4 on the 10 frequency bands obtained by performing 4layer wavelet decomposition on the two microphone acoustic signals.
Fig. 3. Analysis diagram of different wavelet decomposition layers and basis function bands: (s1: Complete oil failure, s2: Moderate oil leakage, s3: Light oil leakage, s4: Normal operating conditions)
a) Wavelet decomposition layer4(db5)
b) Wavelet decomposition layer3(db5)
c) Wavelet decomposition layer5(db5)
d) Wavelet decomposition layer3(db4)
e) Wavelet decomposition layer4(db4)
f) Wavelet decomposition layer4(db5)
2.3. Constructing BP neural network
BP neural network has powerful nonlinear mapping, selflearning and selfadaption, generalization and faulttolerant capabilities, mainly composed of input layer, hidden layer and output layer. If the sample output goes to $S=[{s}_{1},{s}_{2},{s}_{3},{s}_{4}]$, and the training result vector is $R=[{r}_{1},{r}_{2},{r}_{3},{r}_{4}]$, then the square error between the network training results and the target output of a single sample is:
Thus, training system error mean:
where $k$ is the number of samples, we can obtain the weight correction amount $\u2206{V}_{ij}$ and $\u2206{W}_{jk}$ of the hidden layer and the output layer by the value of $D$, and then obtain the network weight correction relationship:
The structure of the BP algorithm is shown in Fig. 4, where $X$(1), $X$(2), ..., $X\left(n\right)$ are the input values of the network, and $Y$(1), $Y$(2), ..., $Y\left(n\right)$ are the predicted values of the network. The input values of the network propagate from the input layer through the hidden layer to the output layer, and the input responses of the network are obtained by each neuron in the output layer. Then, in order to reduce the direction of the target output and the actual error, each connection weight value is corrected layer by layer from the output layer through each intermediate layer, and finally returns to the input layer. Within a certain range, with the back propagation of the error, the accuracy of the feedback of the network test also increases. Through the two processes of forward propagation and backward propagation, the algorithm advances and the error is constantly corrected, and the output value approaches the ideal value.
Fig. 4. BP neural network structure
The eigenvectors composed of the normalized frequency band coefficients are the input values of the BP neural network, and the network outputs correspond to the four working condition setting codes. The failure type (V) is set as the normal oil (1000), slight oil leakage (0100), moderate oil leakage (0010) and complete oil cut (0001), which supply according to the severity of “oil leakage”.
3. Fault diagnosis experiment
3.1. Data collection and processing
The experiment used a 40 PH type microphone from GRAS, and its operating frequency was 1020 kHz. Integrated error intensity measurement mechanism, combined with experimental conditions, take the microphone measurement position from the noise source about 0.8 m, the specific location shown in Fig. 6, 1 for the single sound pressure measurement microphone, numbered 1 and 2 The microphone is a microphone set for measuring double sound pressure. The experimental sampling frequency is 25.6 kHz. We take the Weichai WD615 sixcylinder engine as the experimental object, set the highpressure tubing oil leakage fault, and use the NI data acquisition system to collect acoustic signals. The engine speed is 710 r•min^{1}.
Fig. 5. Microphone measurement position
By adjusting the amount of fuel injected into the cylinder by the high pressure fuel pipe of the sixth cylinder to simulate the oil leakage fault, which included the high pressure fuel pipe, the normal working condition, slight leaking condition, moderate leaking condition and completely broken condition. Fixed sound signal acquisition device (microphone), two microphones simultaneously measure the acoustic signal of the engine under four operating conditions, the sampling time for each working condition is 5 s, 5 sets for each sample.
The wavelet base function db4 is used to decompose the acoustic signal into four layers to obtain five frequency band coefficient values, namely the wavelet coefficient S1 of the high frequency signal and the low frequency wavelet coefficients S2, S3, S4, S5. After summing the wavelet coefficients of the signals of the five frequency bands obtained by decomposing the acoustic signal, they are summed and further normalized to obtain ${e}_{1},{e}_{2},{e}_{3},{e}_{4},{e}_{5}$ as a characteristic parameter of the signal. Therefore, for a single microphone feature vector, there are:
For dual microphone feature vectors:
Pattern recognition of the energy band characteristic parameters of four different working conditions, 500 working cycles (feature vectors) collected by the microphone under each working condition are divided into sample data and test data, among which 400 groups are randomly selected for BP neural network samples Training, another 100 groups as test samples.
3.2. BP network structure training and results
When the BP neural network is used for experimental analysis, since the experiment includes four conditions of normal engine oil supply, slight oil leakage, moderate oil leakage, and complete oil cutoff, the type of sample to be classified is 4, and the number of nodes corresponding to the output layer is also It is determined as 4 accordingly. According to the Kolmogorov theorem, the number of neurons in the hidden layer m and the number of neurons in the input layer $n$ and the optimal range of the output layer nodes $q$ are: where $a$ is in the range of 1 to 10. After repeated experiments, the number of neurons in the hidden layer is set to 9 to obtain better results. The neuron transfer function selects the tangent Sshaped transfer function, the training function selects the LM algorithm, the learning function selects the gradient descent momentum learning function, and the network performance function is the mean squared error function. After many experiments, the learning rate is set to 0.001 for better training. The effect is that the maximum number of iterations is 600.
Under the same conditions, using BP to extend into the network and fuzzy clustering for pattern recognition, the results are shown in Table 2. From the Table, it can be found that with the same algorithm, dual microphones are superior to single microphones in terms of diagnostic accuracy. After many experiments, the average training time of the BP neural network is about 2 seconds, and the fuzzy clustering training time is about 10 seconds. Therefore, the BP neural network is selected as the fault identification method.
Table 1. Comprehensive analysis of network training
Number of microphones

Pattern recognition method

Average accuracy/%

Wrong distribution area

Training time/s

Dual

BP neural network

92

Mild5, Moderate3

2.767

Fuzzy clustering

91

Mild3, Moderate5, Severe1

3.417


Single

BP neural network

73

normal3, Mild11, Moderate12, Severe1

9.113

Fuzzy clustering

66

normal12, Mild12, Moderate9, Severe1

1.653

4. Conclusions
By selecting the wavelet basis function and the number of decomposition layers, using two microphones to collect engine acoustic signals, and adopting the method of wavelet multiresolution analysis combined with BP neural network to diagnose diesel engine faults, good results are obtained, as follows:
(1) For engine acoustic signals, the wavelet basis function db4 can well decompose the lowfrequency part of the acoustic signal and can extract obvious fault features.
(2) Using dual microphones to collect sound pressure signals, combined with BP neural network to diagnose abnormal engine noise, compared with a single microphone, can greatly improve the diagnostic accuracy and diagnostic efficiency, provided a reference for engine abnormal noise diagnosis.
References
 Wang Xiaolong, Tang Guiji. A new method for early fault diagnosis of rolling bearings based on continuous wavelet transform. Propulsion Technology, Vol. 37, Issue 8, 2016, p. 14311437. [Search CrossRef]
 Pan Yi, Liu Xiong, Yan Zijian. Intelligent fault diagnosis of rolling bearings based on wavelet packet band energy characteristics and BP Neural network. Automation and Instrumentation, Vol. 5, 2015, p. 8284. [Search CrossRef]
 Yao J., Xiang Y., Qian S., et al. Noise source identification of diesel engine based on variational mode decomposition and robust independent component analysis. Applied Acoustics, Vol. 116, 2017, p. 184194. [Publisher]
 Guo Yingying, Zhang Chunliang, Yue Xia, et al. Application of morphological component analysis of acoustic signals in bearing fault diagnosis. Mechanical and Electrical Engineering Technology, Vol. 45, Issue 4, 2016, p. 5561. [Search CrossRef]
 Kemalkar A. K., Bairagi V. K. Engine fault diagnosis using sound analysis. International Conference on Automatic Control and Dynamic Optimization Techniques (ICACDOT), 2016, p. 943946. [Publisher]
 Ma S., Wang H. Based on wavelet packet and the LSSVM ship with warship engine fault diagnosis. Control and Decision Conference (CCDC), Vol. 2016, 2016, p. 26872691. [Search CrossRef]
 Qu J., Zhang Z., Gong T. A novel intelligent method for mechanical fault diagnosis based on dualtree complex wavelet packet transform and multiple classifier fusion. Neurocomputing, Vol. 171, 2016, p. 837853. [Publisher]
 Liu Qing, Zhu Bing An. Fault diagnosis of gearbox based on wavelet analysis. Mechanical Design and Manufacture, Vol. 11, 2016, p. 245248. [Search CrossRef]