Misfire fault diagnosis of diesel engine based on VMD and XWT

Jide Jia1 , Gang Ren2 , Jianmin Mei3

1School of Mechanical and Manufacturing Engineering, Xiamen Institute of Technology, Xiamen, 361024, China

1, 3Projection Equipment Support Department, Army Transportation University, Tianjin, China

2Automobile NCO School, Army Military Transportation University, Bengbu, 233011, China

2Corresponding author

Vibroengineering PROCEDIA, Vol. 24, 2019, p. 17-22. https://doi.org/10.21595/vp.2019.20767
Received 25 April 2019; accepted 4 May 2019; published 7 June 2019

Copyright © 2019 Jide Jia, 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 References
Cite this article
Views 44
Reads 18
Downloads 923
Abstract.

For diesel engine misfire fault diagnosis under strong noise, a new method based on variational mode decomposition (VMD) and cross wavelet transform (XWT) is proposed. Firstly, the vibration signal of cylinder head is processed by stages according to the working cycle of diesel engine, and then the vibration signal of each working cycle is re-sampled according to the same angle. Then the vibration signal of cylinder head is decomposed by VMD to denoise adaptively and reconstruct the signal. Next, XWT is used to analyze the time-frequency correlation of any two continuous working cycle signals, and the non-periodic components and random noise in the vibration signals are further eliminated to extract the combustion characteristics of diesel engines. Finally, the misfire fault of diesel engine is diagnosed by calculating the energy proportion of each cylinder in time-frequency space. The effectiveness of the proposed method is verified by simulation and experiment.

Keywords: variational modal decomposition, cross wavelet transform, feature extraction, fault diagnosis, diesel engine.

1. Introduction

There are many factors affecting the combustion of diesel engine and the working process is complex, which results in the multi-component and non-stationary characteristics of the signal. It is difficult to monitor and diagnose the combustion state of diesel engine. Domestic and foreign scholars have done a lot of research on diesel engine misfire fault, and put forward a variety of diagnostic methods, such as cylinder pressure diagnostic method [1], instantaneous speed diagnostic method [2], exhaust noise diagnostic method, and cylinder head vibration and sound diagnostic method [3]. However, no matter what method is used, feature extraction is one of the key steps in the diagnosis process. Whether the extracted feature information can characterize the combustion characteristics of diesel engine is directly related to the accuracy of diagnosis results.

Empirical mode decomposition (EMD) is used to extract eigenvalues, which can reconstruct signals and denoise. However, EMD has the problems of endpoint effect and modal aliasing. Although a series of improvements have been made, such as EEMD, CEEMDAN, endpoint effect is still unavoidable, as well as the problems of large amount of computation and poor anti-noise ability. Wavelet transform (WT) is used to extract eigenvalues, which has been widely used in the field of diesel engine fault diagnosis. However, WT can only analyze a single time series signal, and it is difficult to analyze the interaction and time-frequency correlation between two time series.

In recent years, Konstantin Dragomiretskiy [4] proposed a new variable adaptive decomposition method, Variational Mode Decomposition (VMD). The method is essentially composed of a plurality of adaptive Wiener filter and has good noise robustness. Compared with EMD, VMD has strong mathematical theory basis. At the same time, VMD rejects the method of recursive screening stripping. So VMD can effectively alleviate or avoid a series of problems which appear in EMD and has high operation efficiency of decomposition [5]. VMD has been fully utilized in the field of engineering.

Cross-wavelet transform (XWT) is a new signal analysis technology based on traditional wavelet transform. Through the combination of wavelet transform and cross-spectrum analysis, the relationship between two time series in time-frequency domain can be studied from the perspective of multi-time scale. At the same time, the different performance of signal and noise in time-frequency correlation is utilized to achieve the purpose of denoising and signal enhancement [6].

After a lot of research, a new diagnosis method for diesel engine misfire fault based on VMD and XWT is presented in this paper. The proposed method is applied to diesel engine misfire fault diagnosis. The fault feature is successfully extracted, and diesel engine misfire fault can be diagnosed effectively.

2. Feature extraction based on VMD and XWT

2.1. VMD

Detailed algorithms for VMD can be found in reference [4]. The variational constraint problem can be defined as follows:

(1)
m i n u k , w k k = 1 K a t { [ ( δ ( t ) + j π t ] * u k ( t ) } e - j ω k t 2 2     s . t .     k = 1 K f t = u k t ,

where uk is the modal component, wk is the central frequency for the modal component, δ(t) is the unit pulse function, and * is the convolution symbol.

In the VMD algorithm, the secondary penalty factor and the Lagrangian multiplication operator are used. Then, the alternating direction method is introduced. ukn+1, ωkn+1, and λ n+1 are constantly updated.

The expression for the modal component ukn+1 in frequency domain is:

(2)
u ^ k n + 1 ( ω ) = f ^ ω - i u ^ i ω + λ ^ ω 2 1 1 + 2 α ( ω - ω k ) 2

where ωk  is the center of the modal component power spectrum. The Wiener filter is introduced, which makes the VMD algorithm have better noise robustness.

Similarly, the expression for the central frequency ωkn+1:

(3)
ω k n + 1 ω = 0 ω u ^ k ω 2 d ω 0 u ^ k ω 2 d ω .

The stopping condition of the iteration is:

(4)
k = 1 K u ^ k n + 1 - u ^ k n 2 2 u ^ k n 2 2 < e .

The VMD algorithm is a linear transformation, so the signal can be reconstructed. The reconstructed signal can be represented as:

(5)
f ^ t = k = 1 K u ^ k ,

where u^k is the final modal component, after the iteration is stopped.

2.2. XWT

The continuous wavelet transform of the signal x(t) is:

(6)
W x a , b = - x t ψ a , b * t d t = x t , ψ a , b t ,
(7)
ψ a , b t = a - 1 2 ψ t - b a ,

where a is the scale factor and b is the translation factor; ψa,b(t) is a wavelet basis function, and * denotes complex conjugation.

Wavelet transform expands one-dimensional signal in time-frequency domain and can analyze any details of signal change. For any two signals xt and yt, the XWT of two signals established on the basis of wavelet transform can be expressed as [7]:

(8)
C x , y a , b = W x a , b W y * a , b .

The coefficients of XWT represent the correlation of two signals with wavelet coefficients at the same scale. XWT can analyze the correlation between the two signals in both time and frequency domains. It reflects the correlation between the two signals in the time-frequency domain. The larger the correlation, the closer the correlation between the two signals [6].

Morlet wavelet is chosen as the wavelet basis function in XWT analysis, which has good localization properties in both time and frequency domains. The expressions in time domain and frequency domain are as follows:

(9)
ψ t = π - 1 4 e x p j ω 0 t e x p - t 2 2 ,
(10)
ψ ^ ( a ω ) = π 1 4 e x p - a ω - ω 0 2 2 ,

where the value of ω0 is 6 to ensure good local performance in time-frequency distribution.

2.3. Proposed method

Firstly, the vibration signal of cylinder head is re-sampled in segments and equal angles. Then the VMD is used to decompose and reconstruct the cylinder head vibration signal, and XWT is used to analyze the time-frequency correlation. Finally, according to the energy of each cylinder in the time-frequency correlation spectrum, the misfire diagnosis of diesel engine is carried out.

3. Experimental results

3.1. Simulation

To verify the proposed method, a multi-component simulation signal consisting of four AM-FM components and noise is established as follows:

(11)
S = s 1 + s 2 + s 3 + s 4 + n t ,
(12)
s 1 = e x p - a t s i n 2 π f 1 k 1 t 1 ,
(13)
s 2 = e x p - a t s i n 2 π f 2 k 2 t 2 ,
(14)
s 3 = e x p - a t s i n 2 π f 3 k 3 t 3 ,
(15)
s 4 = e x p - a t s i n 2 π f 4 k 4 t 4 ,

where t=modkdt,1/fm; k is a natural number; dt= 5×10-5; fm= 10 Hz; a= 1800 Hz; k1= 1:250; k2= 1:250; k3= 1:250; k4= 1:150; f1= 1 kHz; f2= 3 kHz; f3= 5 kHz; f4= 7 kHz. The signal sampling frequency is 1500Hz, and the sampling number is 2000 Hz.

WT, VMD and VMD+XWT are used to analyze the time-frequency of the simulation signal, and the time-frequency analysis results are shown in Fig. 1.

Fig. 1. The time-frequency diagram of simulation signal

 The time-frequency diagram of simulation signal

a)

 The time-frequency diagram of simulation signal

b)

 The time-frequency diagram of simulation signal

c)

The change of contour nephogram corresponds to the change of signal energy in time-frequency space. The light-colored area corresponds to the higher intensity energy value, while the dark-colored area corresponds to the lower intensity energy value. Other time-frequency diagrams in this paper show the same. From Fig. 1(a), it can be seen that WT has insufficient noise attenuation for simulation signal. From Fig. 1(b), it can be seen that VMD has obvious attenuation for simulation signal noise. From Fig. 1(c), it can be seen that the method based on VMD and XWT has a great attenuation to the simulation signal noise. By combining VMD and XWT, the main components of the simulation signal are extracted completely, and the interference noise is basically eliminated.

3.2. Experiment condition

The vibration signal of engine cylinder head contains abundant information, which can effectively reflect the change of combustion state of diesel engine. The cylinder head vibration signal can be used to monitor and diagnose the combustion state of diesel engine. The measurement method is reliable and can realize on-line monitoring and diagnosis. The test object is HJ493 diesel engine.

A 601A01 ICP vibration sensor is installed on the cylinder head between the second cylinder and the third cylinder of the diesel engine. QSY8104 external card type high-pressure tubing pressure sensor is clamped on the first cylinder high-pressure tubing of diesel engine. The speed of the engine is measured by the oil pressure pulse signal of the high-pressure tubing and the position of the first cylinder stop point is measured. The signal sampling frequency is 12.8 kHz and the data sample length is 12288. Five working conditions of diesel engine are simulated by controlling the fuel supply of each cylinder injector: normal combustion state of diesel engine, poor combustion of the first cylinder, poor combustion of the second cylinder, poor combustion of the third cylinder and poor combustion of the fourth cylinder. Vibration signals of cylinder head surface were collected at 700 r/min of diesel engine speed under normal and abnormal conditions.

3.3. Experimental data processing

According to the proposed method, the vibration signals of cylinder head in five states are processed respectively. The effects of aperiodic component and random noise are effectively attenuated, and the combustion state characteristics are enhanced. The results are shown in Fig. 2.

Fig. 2. Three-dimensional time-frequency diagram of vibration signals in different combustion states

 Three-dimensional time-frequency diagram of vibration signals in different combustion states

a)

 Three-dimensional time-frequency diagram of vibration signals in different combustion states

b)

 Three-dimensional time-frequency diagram of vibration signals in different combustion states

c)

 Three-dimensional time-frequency diagram of vibration signals in different combustion states

d)

 Three-dimensional time-frequency diagram of vibration signals in different combustion states

e)

As can be seen from Fig. 2, when the diesel engine is in normal combustion state, the combustion state of the four cylinders is basically the same. The energy of each cylinder is different, but the difference is not significant. When a cylinder is misfired, the combustion balance of the whole diesel engine is destroyed, and the energy of the misfired cylinder decreases most obviously. Therefore, the cylinder energy is selected as the characteristic parameter of diesel engine combustion state change, and the proportion of cylinder energy is used as the quantitative evaluation criterion for diagnosing diesel engine misfire. When the energy proportion of a cylinder is less than or equal to 10 %, it indicates that the cylinder has a misfire fault. The results of quantitative diagnosis of combustion state of each cylinder are shown in Table 1.

Table 1. Energy proportion of cylinders in different combustion states

Cylinder number
1
2
3
4
Normal combustion
28 %
17 %
22 %
32 %
First cylinder misfire
10 %
20 %
33 %
37 %
Second cylinder misfire
45 %
6 %
32 %
16 %
Third cylinder misfire
70 %
14 %
4 %
11 %
Fourth cylinder misfire
48 %
12 %
35 %
5 %

As can be seen from Table 1, when the first cylinder is misfired, the energy of the first cylinder is obviously less than that of other cylinders, accounting for only 10 % of the total energy. When the second cylinder is misfired, the energy of the second cylinder is also obviously less than that of other cylinders, accounting for only 6 % of the total energy. When the third cylinder is misfired, the energy of the third cylinder is also obviously less than that of other cylinders, accounting for only 4 % of the total energy. When the fourth cylinder is misfired, the energy of the fourth cylinder is also obviously less than that of other cylinders, accounting for only 5 % of the total energy.

4. Conclusions

Both VMD and XWT have certain ability of noise reduction in signal analysis. However, the analysis of cylinder head vibration signal under strong background noise of diesel engine cannot achieve the desired effect. Using the method based on VMD and XWT to process the vibration signal of diesel engine cylinder head can effectively eliminate the influence of interference noise and enhance combustion characteristics, which lays a foundation for diesel engine combustion quality monitoring and misfire diagnosis. The proposed method is applied to the misfire diagnosis of diesel engines. The results show that the method can eliminate the non-periodic component and random interference in the cylinder head vibration signal, and significantly enhance the fault characteristics. The proposed method can effectively identify the misfire fault of diesel engine.

References

  1. Payri F., Luján Jm, Martín J., Abbad A. Digital signal processing of in-cylinder pressure for combustion diagnosis of internal combustion engines. Mechanical Systems and Signal Processing, Vol. 24, Issue 6, 2010, p. 1767-1784. [Publisher]
  2. Li Z., Yan X., Peng Z. Intelligent fault diagnosis method for marine diesel engines using instantaneous angular speed. Journal of Mechanical Science and Technology, Vol. 26, Issue 8, 2012, p. 2413-2423. [Publisher]
  3. Liu J., Li X., Zhang X., et al. Misfire diagnosis of diesel engine based on rough set and neural network. Procedia Engineering, Vol. 16, Issue 16, 2011, p. 224-229. [Search CrossRef]
  4. Dragomiretskiy K., Zosso D. Variational mode decomposition. IEEE Transactions on Signal Processing, Vol. 62, Issue 3, 2013, p. 531-544. [Publisher]
  5. Zhao C., Feng Z. P. Application of multi-domain sparse features for fault identification of planetary gearbox. Measurement, Vol. 104, 2017, p. 169-179. [Publisher]
  6. Jide J. Combustion feature enhancement in diesel engine based on cross-wavelet transform and teager energy operator. Automotive Engineering, Vol. 36, Issue 2, 2014, p. 164-167. [Search CrossRef]
  7. Torrence C., Compo G. P. A practical guide to wavelet analysis. Bulletin of the American Meteorological Society, Vol. 79, Issue 1, 1998, p. 61-78. [Publisher]

Cited By

Vibroengineering PROCEDIA
Jide Jia, Gang Ren
2021