An improved variational mode decomposition method and its application in diesel engine fault diagnosis

Gang Ren1 , Jide Jia2 , Jianmin Mei3 , Xiangyu Jia4 , Jiajia Han5 , Yu Wang6

1, 4, 5Fifth Cadet Brigade, Army Military Transportation University, Tianjin, China

2, 3Projection Equipment Support Department, Army Military Transportation University, Tianjin, China

6School of information engineering, Southwest University of Science and Technology, Sichuan, China

3Corresponding author

Journal of Vibroengineering, Vol. 20, Issue 6, 2018, p. 2363-2378. https://doi.org/10.21595/jve.2018.19694
Received 1 February 2018; received in revised form 30 May 2018; accepted 8 June 2018; published 30 September 2018

Copyright © 2018 Gang Ren, 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
Abstract.

The diesel engine is a complex mechanical device, with the characteristics of multi-source, multi moving parts, complex work. For the complex multi-component signal, it is usually necessary to decompose it into a number of single-component AM-FM signals, and each component is analyzed to extract amplitude and frequency information. VMD is essentially composed of a plurality of adaptive Wiener filter and has good noise robustness. Compared with EMD, EEMD, CEEMDAN, LMD and ITD, 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 other methods. However, it is a problem how to determine the number of decomposition layers and the penalty factor, because human factors will affect the decomposition results. In order to solve the problem, an improved adaptive genetic algorithm (IAGA) is proposed to optimize the parameters of VMD. Genetic algorithms mainly include 3 genetic operators: selection, crossover and mutation. The cross probability and mutation probability will directly affect the optimization results. In the traditional genetic algorithm, the probability of cross and mutation are fixed, and the genetic algorithm is easy to fall into the local optimal. According to the regulation of hormone regulation, the cross probability and mutation probability in evolution were improved. The permutation entropy is a new method of mutation detection, which mainly aims at the spatial characteristics of the time series itself. Therefore, the entropy of the components obtained by the VMD decomposition is used as the fitness function of the IAGA. The modal number K and penalty factor α of VMD were iteratively optimized by IAGA, and the optimal combination of parameters was obtained. Based on the proposed method, the vibration signals of the crankshaft bearing fault simulation experiment were decomposed into several components. According to the value of the permutation entropy, the fault components were selected and the energy was extracted. The fault pattern is identified by the support vector machine (SVM) successfully. The simulation analysis and the simulation experiment of the crankshaft bearing fault show that the proposed method is effective. For the diagnosis of other engines, a large number of validation experiments are needed for further research.

Keywords: improved variational mode decomposition, adaptive genetic algorithm, fuel engine, fault diagnosis.

1. Introduction

The engine is a complex mechanical device, with the characteristics of multi-source, multi moving parts, complex work. The engine has both rotational and reciprocating motion. Vibration signals are fully used in fault diagnosis because of their convenience [1, 2]. The vibration signal of engine is composed of multi-component complex signals, and its amplitude varies with time. For the complex multi-component signal, it is usually necessary to decompose it into a number of single-component AM-FM signals, and each component is analyzed to extract amplitude and frequency information.

Huang et al proposed the Empirical Mode Decomposition (EMD) [3, 4]. EMD has been widely used in mechanical fault diagnosis since it was proposed [5, 6]. EMD is fundamentally divorced from the Fourier transform. EMD analyzes the signal from itself and it is completely adaptive and unsupervised. But there are some problems for EMD, such as endpoint effect and modal aliasing [7].

Wu [8] proposed Ensemble Empirical Mode Decomposition (EEMD). Different white noises are added to the original signal for EMD, and multiple decomposition results are averaged to obtain the final Intrinsic Mode Function (IMF). The high frequency modulation information in the signal can be separated very well, and the modal aliasing of EMD is well suppressed [9]. Torres [10] proposed a Complete Ensemble Empirical Mode Decomposition with Adaptive Noise (CEEMDAN). The method reduces the reconstruction error caused by white noise while suppressing the aliasing of modes. However, the two methods share a common flaw. Their computation is very large and a lot of pseudo components will appear in the decomposition process.

Smith proposed the Local Mean Decomposition (LMD) [11]. LMD considers that a complex single-component signal is the product of its own envelope signal and an FM signal, which is called Product Function (PF). A complex signal is decomposed into several PF components with the physical meaning of instantaneous frequency. LMD avoids the problems of over envelope, incomplete envelope and negative frequency caused by Hilbert transform in EMD. However, LMD also has some problems such as frequency aliasing and endpoint effect [12]. Mark G Frei [13] proposed the Intrinsic Time-scale Decomposition (ITD). Any complex signal can be decomposed into several independent Proper Rotation (PR) components, and its instantaneous frequency is of physical significance. However, baseline definition in LMD is based on the linear transformation of the signal itself. So, from the second components, the components are different from the usual definition of IMF, and an apparent signal distortion occurs. Then the instantaneous amplitude and instantaneous frequency are distorted greatly [12].

In recent years, Konstantin Dragomiretskiy [14] 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 [15]. VMD has been fully utilized in the field of engineering [16-19].

However, it is a problem that how to determine the number of decomposition layers and the penalty factor, because human factors will affect the decomposition results. Genetic algorithm is a global optimization random probability search algorithm that simulates the inheritance and variation of chromosomes in a population based on natural selection and natural genetic mechanism in the biological world. Genetic algorithm is often used in parameter optimization, and has fast convergence speed and wide adaptability. Genetic algorithms mainly include 3 genetic operators: selection, crossover and mutation. The cross probability and mutation probability will directly affect the optimization results. In the traditional genetic algorithm, the probability of cross and mutation are fixed, and the genetic algorithm is easy to fall into the local optimal [20]. M. Srinvas [21] proposed an adaptive genetic algorithm, which could adaptively adjust the cross probability and mutation probability in the genetic algorithm. As a result, the convergence rate is improved, but there is an evolutionary stagnation at the early stage of evolution. In order to solve the problem, an improved adaptive genetic algorithm (IAGA) was proposed to optimize the modal number K and penalty factor α of VMD in this paper. According to the law of hormone regulation, the cross probability and mutation probability in evolution were improved. Compared with the traditional adaptive genetic algorithm, the convergence speed is faster and the local optimum is avoided. Based on the proposed method, the vibration signal of diesel engine crankshaft was analyzed. The result shows that the proposed method can quickly converge to the global optimal value and have high stability.

2. Variational mode decomposition

The VMD algorithm defines the intrinsic mode function as a non-stationary AM-FM signal. The intrinsic mode is considered as follows:

(1)
u k t = A k t c o s φ k t ,

where the phase φkt shall satisfy the following condition: φkt' 0; the envelope line Akt should satisfy the following condition: Akt 0; the instantaneous frequency ωkt should satisfy the following condition: ωkt=φkt'. Akt and ωkt change slowly, and φkt changes more rapidly.

The Hilbert transform is performed for each modal function ukt, and exponential correction is applied to obtain K modal functions. Then the frequency spectrum of the modal function is corrected to the estimated central frequency, and the bandwidth of the modal component is calculated by using Gauss smoothing. The variational constraint problem can be defined as follows:

(2)
m i n u k , ω k k = 1 K a t δ t + j π t * u k t e j ω k t 2 2 ,       s t k = 1 K u k t = f t ,

where uk is the modal component, ωk 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, so that the optimal solution of the variational constraint problem can be solved. The expression for the modal component ukn+1 is:

(3)
u k n + 1 = a r g m i n u k X α a t δ t + j π t * u k t e j ω k t 2 2 + f t i u i t + λ t 2 2 2 ,

where α is the penalty factor, and λ is the Lagrange multiplier.

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

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

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 is:

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

The stopping condition of the iteration is:

(6)
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:

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

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

3. Improved adaptive genetic algorithm

Genetic algorithm is a global optimization random probability search algorithm that simulates the inheritance and variation of chromosomes in a population based on natural selection and natural genetic mechanism in the biological world. Genetic algorithms mainly include 3 genetic operators: selection, crossover and mutation. The cross probability and mutation probability will directly affect the optimization results. When the fitness of individuals is lower than the average fitness, the probability of crossover and mutation should be improved. If the fitness of individuals is higher than the average fitness, the probability of crossover and mutation should be reduced. The cross probability and mutation probability in the traditional genetic algorithm are fixed. As a result, the genetic algorithm is easy to fall into the local optimal.

M. Srinvas [21] proposed an adaptive genetic algorithm to adaptively adjust the cross probability and mutation probability in the genetic algorithm:

(8)
p c = k 1 f m a x - f ' f m a x - f m i n ,       f ' f a v g ' , k 2 ,           f ' < f a v g ' .
(9)
p m = k 3 f m a x - f f m a x - f m i n ,       f f a v g , k 4 ,         f < f a v g .

where Pc is the cross probability; Pm is the mutation probability; f' is the larger fitness of the two individuals involved in the cross; f is the individual’s fitness value for a mutation operation; fmax is the maximum fitness value in the population; fmin is the smallest fitness value in the population; k1, k2, k3 and k4 are constants in the range of (0,1). According to Eqs. (12) and (13), the convergence rate of the adaptive genetic algorithm is improved. However, there is an evolutionary stagnation at the early stage of evolution for the adaptive genetic algorithm.

In order to solve the above problems, according to the law of hormone regulation, the following improvements are made to optimize the cross probability and mutation probability in evolution:

(10)
p c = P c 0 1 + γ f a v g k c f a v g - f m i n k c + f a v g k c ,
(11)
p m = P m 0 1 + β f a v g k m f a v g - f m i n k m + f a v g k m ,

where Pc0 is the initial cross probability; Pm0 is the initial mutation probability; γ, β, kc and km are coefficient factors. Through the above improvement, the individual is more diverse, and the problems of premature convergence and slow evolution are effectively solved.

4. Proposed method

According to the VMD algorithm introduced in the chapter 2, VMD needs a given decomposition mode number K during the decomposition process, and the value of K has a great influence on the decomposition results. In addition, the penalty factor α in the algorithm also has a great influence on the decomposition results. If the modal number K and the penalty factor α are set artificially, the decomposition results will be subject to subjective influence. It is a prerequisite for VMD that how to set these two parameters.

The IAGA is a global optimization random probability search algorithm, which avoids the problem of premature convergence and slow evolution. Therefore, the IAGA is used to optimize the modal number K and penalty factor α of VMD in this paper, and the optimal input parameters are obtained.

In the process of optimization, IAGA needs to define a fitness function. The permutation entropy is a new method of mutation detection, which mainly aims at the spatial characteristics of time series [22, 23]. The permutation entropy is very simple in theory and has good noise robustness. Besides, the permutation entropy also has a high resolution, and the output results are very intuitive. The permutation entropy reflects the random degree of the signal. In other words, the smaller the permutation entropy, the more regular the time series is; the larger the entropy, the more random the time series is. In phase space, a set of time series Xi is reconstructed to get a set of symbol sequence, and the probability of each symbol sequence appears as P1, P2, …, Pk. Then the permutation entropy of the time series Xi can be defined as:

(12)
H p m = - j = 1 k P j l n P j ,

when the crankshaft of a diesel engine fails, it is often accompanied by the violent vibration and abnormal noise of the body. If the fault features are obvious, the signal has periodic fluctuation and the permutation entropy is relatively small. If there is a lot of noise in the signal, the randomness of the signal is large, and the permutation entropy is relatively large. Therefore, the fitness function of IAGA is the permutation entropy of the components obtained by VMD. When the signal is decomposed by VMD, the minimum permutation entropy is called the local minimum permutation entropy value, and the corresponding component is the local optimal component. Then the local minimal permutation entropy is used as the fitness value of the whole evolutionary process to search the optimal combination of parameters. The signal is decomposed by the optimized VMD, and the component with the minimum permutation entropy is selected as the fault component to extract the feature.

5. Experimental results

5.1. Simulation

The simulation signal is set up to verify the decomposition effect of the VMD optimized by IAGA. The engine structure is complex, and the vibration excitation source is multitudinous. Besides, the engine often operates in transient conditions, and the source vibration signal is modulated by a plurality of components. Therefore, the simulation signals of the engine must satisfy the requirements of multi- component, non-stationary and AM-AF signals. In this paper, Gauss noise is added to the simulation signal, which is more in line with the actual vibration signal characteristics of the diesel engine. The simulation signal is composed of the following:

(13)
a = e - 8 ( t - 0.5 ) 2 ,
(14)
s 1 = a s i n ( 2 π 100 t + c o s 5 π t ,
(15)
s 2 = a s i n 2 π 150 t + 50 t 2 ,
(16)
s 3 = a s i n 2 π 250 t + 80 t 3 ,
(17)
s = s 1 + s 2 + s 3 + 0.1 r a n d n s i z e s 1 .

The signal sampling frequency is 1500 Hz, and the sampling number is 1024. The time domain waveform of the simulation signal components are shown in Fig. 2. As shown in Fig. 2, the simulation signal consists of three non-stationary AM-FM signals. The time-domain and frequency-domain waveforms of the noisy simulation signals are shown in Fig. 3. The composition and variation of the signal are not clear from Fig. 3. We don't know what the signal contains.

Fig. 1. Detailed experimental scheme

Detailed experimental scheme

Fig. 2. The source components of simulation signal

The source components of simulation signal

Fig. 3. The simulation signal and its spectrum

The simulation signal and its spectrum

a)

The simulation signal and its spectrum

b)

In order to further analyze the signal, the simulation signal is decomposed by the VMD optimized by IAGA. In the IAGA, the population size is 8, and the evolution subalgebra is 100. In the process of IAGA optimization of modal number K and penalty factor α, the curve of the best fitness values in each generation is shown in Fig. 4.

Fig. 4. The curve of the best fitness values in each generation

The curve of the best fitness values in each generation

As shown in Fig. 4, the best fitness value of the proposed method appears in the 36th generation, and the best fitness value of the traditional genetic algorithm appears in the 50th generation. In comparison, the proposed method converges faster and inhibits the precocious phenomenon. After the end of evolution, the optimal parameters obtained by the proposed method are K,α= 3,646. Therefore, the modal number K is set to 3, and the penalty factor α is set to 646. The simulation signal is decomposed by VMD, and the decomposition results are shown in Fig. 5.

As shown in Fig. 5, the time domain waveform of the component of VMD decomposition is in good agreement with the time domain waveform of the signal source signal. That is to say, the proposed method has achieved good decomposition effect. In addition, in order to verify the advantages of the proposed method, EMD, EEMD, CEEMDAN, LMD and ITD are used to decompose the simulation signal, as shown in Fig. 6.

The energy of the first three components accounts for more than 90 % of the total energy of the simulation signal, so the first three components were selected for analysis. By comparing the signal source component in Fig. 2 with the signal decomposition component in Fig. 5 and Fig. 6, we can tentatively judge that the effect of the proposed method is the best. In order to further compare the authenticity of the components obtained by different methods, correlation coefficient is used as evaluation index. The correlation coefficients between the components decomposed by different methods and the simulation signal sources are compared, as shown in Table 1.

Fig. 5. The decomposition results of the proposed method

The decomposition results of the proposed method The decomposition results of the proposed method The decomposition results of the proposed method

Table 1. Correlation between decomposition components and simulation signal source components

Decomposition methods
Component 1
Component 2
Component 3
VMD
0.99
0.99
0.98
EMD
0.85
0.73
0.86
EEMD
0.57
0.76
0.86
CEEMDAN
0.84
0.4
0
LMD
0.1
0.03
0.01
ITD
0.86
0.6
0.05

According to Table 1, it is easy to see that, compared with other methods, the components obtained by the proposed method are more closely related to the simulation signal source components. Thus, the components obtained by the proposed method are closer to the real signal, and the decomposition effect is the best.

5.2. Experiment condition

The structure of diesel engine is complex, and the working environment is abominable. As a result, it is prone to malfunction. The crankshaft bearing is located inside the engine, so it is difficult to diagnose the fault. In this paper, vibration signals are collected from the vibration sensors on the experimental stand, as shown in Fig. 7. The basic parameters of the vibration sensor are shown in Table 2. The engine on the experimental stand is Cummins 6BT diesel engine, and its parameters are shown in Table 3.

Table 2. Vibration sensor parameters

Model
Sensitivity
Frequency range (±3 dB)
Range
Resolution
Temperature
range
Weight
Output connector
603C01
100 mV/g
0.5 Hz-10 KHz
±50 g
350 μg
–54-121℃
51 g
Top

The fourth crankshaft bearings of Cummins EQ6BT diesel engine are set with different clearance (0.10 mm, 0.26 mm, 0.40 mm, 0.55 mm) to simulate the normal, slight, moderate and heavy wear of the crankshaft bearing. Vibration signals are collected on the left side of the fourth main bearings on the surface of the engine block. The sampling frequency is 20000 Hz and the sampling points are 4096 points.

Testing temperature is important when acquiring vibration signals. In the experiment, the temperature of cooling water is measured to reflect the internal temperature of diesel engine. The temperature is controlled at 60-70°C.

Table 3. Basic parameters of the engine

Engine type
6BT5.9-G2
Fuel type
Diesel oil
Type
Inline 6 cylinders
Rated power (KW)
118
Compression ratio
17.5:1
Ignition sequence
153624
Rated speed (RPM)
2600
Continuous
Power (KW)
86
Maximum
torque (N·m)
558
Radius (mm)
×Distance (mm)
102×120
Maximum torque
speed (r/min)
1600

Fig. 6. The decomposition results of other methods

The decomposition results of other methods

a) EMD

The decomposition results of other methods

b) EEMD

The decomposition results of other methods

c) CEEMDAN

The decomposition results of other methods

d) LMD

The decomposition results of other methods

e) ITD

Fig. 7. Measuring position of vibration sensor

 Measuring position of vibration sensor

5.3. Data acquired

The acquisition system is composed of collector, computer, sensor and connecting circuit, as shown in Fig. 8. The acquisition system set the speed of the engine to 1800 r/min.

Fig. 8. Vibration signal acquisition system

 Vibration signal acquisition system

The vibration signals of the engine under different wear conditions are collected, as shown in Fig. 9. In Fig. 9, there is a large amount of background noise in the vibration signals of different wear conditions of the crankshaft bearing. The shock component of the signal is not obvious, and the fault feature cannot be extracted effectively.

5.4. Experimental data processing

As an example, the slight wear signal of the crankshaft is analyzed by the method proposed. The vibration signal is decomposed by the VMD optimized by IAGA. In the IAGA, the population size is 8, and the evolution subalgebra is 100. In the process of IAGA optimization of modal number K and penalty factor α, the curve of the best fitness values in each generation is shown in Fig. 10.

As shown in Fig. 10, the best fitness value of the proposed method appears in the 19th generation, and the best fitness value of the traditional genetic algorithm appears in the 55th generation. In comparison, the proposed method converges faster and inhibits the precocious phenomenon. After the end of evolution, the optimal parameters obtained by the proposed method are K,α= 5,873. Therefore, the modal number K is set to 5, and the penalty factor α is set to 873. The slight wear vibration signal is decomposed by VMD, and the decomposition results are shown in Fig. 11.

As shown in Fig. 11, the slight wear vibration signal is decomposed into 5 components by VMD. In the same way, the vibration signals of normal wear, moderate wear and heavy wear are decomposed into 5 components. In order to select the appropriate component as the fault component, the permutation entropy of each component under different wear states is calculated respectively, as shown in Table 4.

Fig. 9. Vibration signals in different wear conditions

 Vibration signals in different wear conditions

a) Normal wear

 Vibration signals in different wear conditions

b) Slight wear

 Vibration signals in different wear conditions

c) Moderate wear

 Vibration signals in different wear conditions

d) Heavy wear

Fig. 10. The curve of the best fitness values in each generation

The curve of the best fitness values in each generation

As can be seen from Table 4, for different wear conditions, the entropy permutation value of the component u1 is all minimal. That is, the randomness of the signal is minimal, and the shock component is the strongest. Therefore, the component u1 is selected as the fault component. The spectrums of the component u1 under different wear conditions are shown in Fig. 12.

Table 4. The decomposition component permutation entropy of vibration signals under different wear conditions

Wear condition
u 1
u 2
u 3
u 4
u 5
Normal wear
0.38
0.42
0.56
0.61
0.64
Slight wear
0.36
0.42
0.57
0.69
0.67
Moderate wear
0.33
0.41
0.53
0.55
0.64
Heavy wear
0.30
0.41
0.56
0.62
0.65

Fig. 11. The decomposition components of vibration signal under slight wear condition

The decomposition components of vibration signal under slight wear condition

Fig. 12. The decomposition components of vibration signal under slight wear condition

The decomposition components of vibration signal under slight wear condition

a) Normal wear

The decomposition components of vibration signal under slight wear condition

b) Slight wear

The decomposition components of vibration signal under slight wear condition

c) Moderate wear

The decomposition components of vibration signal under slight wear condition

d) Heavy wear

As can be seen from Fig. 12, the more serious the fault is, the greater the energy of the component u1. In order to further analyze the fault components, the energy of the component u1 under different wear conditions is extracted, as shown in Table 5.

As can be seen from Table 5, for the components obtained by VMD, the change of energy is not consistent with the wear fault of the crankshaft bearing. But for the fault component u1, with the increase of the degree of wear, its energy is also increasing. Therefore, the component u1 can be used as the characteristic component of the fault diagnosis of the crankshaft bearing wear, which is in agreement with the results of the previous permutation entropy. The wear fault diagnosis of the crankshaft bearing can be realized by monitoring the change of the energy of component u1.

Table 5. Basic parameters of the engine

Wear state
u 1
u 2
u 3
u 4
u 5
Normal wear
39.7
102.4
90.9
48.0
125.0
Slight wear
69.1
147.8
49.0
26.8
52.4
Moderate wear
141.4
193.8
85.9
136.1
97.9
Heavy wear
366.5
151.2
76.6
60.3
74.9

In order to accurately analyze the wear failure of diesel engine crankshaft bearing, 40 sets of components u1 energy values are extracted as training samples for support vector machine (SVM), as shown in Table 6. In addition, other 40 sets of components u1 energy values are extracted as test samples for SVM, as shown in Table 7. Specifically, 1, 2, 3 and 4 represents the normal, slight, moderate, and heavy wear of the crankshaft bearing respectively. The results of SVM’s training and test are shown in Fig. 13.

Table 6. Training samples

Feature parameters
Fault type
Feature parameters
Fault type
Feature parameters
Fault type
Feature parameters
Fault type
38.7
1
76.4
2
153.8
3
350.0
4
45.6
1
71.2
2
153.4
3
374.6
4
40.1
1
71.2
2
223.4
3
357.2
4
41.4
1
71.6
2
215.3
3
349.8
4
40.6
1
73.4
2
222.8
3
343.9
4
40.1
1
77.0
2
223.1
3
402.0
4
38.3
1
76.5
2
222.3
3
357.8
4
40.9
1
73.6
2
222.4
3
338.8
4
40.9
1
73.4
2
218.1
3
350.6
4
39.4
1
75.1
2
211.3
3
367.1
4

Table 7. Test samples

Feature parameters
Fault type
Feature parameters
Fault type
Feature parameters
Fault type
Feature parameters
Fault type
43.1
1
70.3
2
150.4
3
388.8
4
45.1
1
75.7
2
165.6
3
346.3
4
39.6
1
71.4
2
203.5
3
353.2
4
39.8
1
70.2
2
201.8
3
347.5
4
39.5
1
69.6
2
216.3
3
407.2
4
40.0
1
74.5
2
163.9
3
371.6
4
99.7
1
76.5
2
212.9
3
356.1
4
39.6
1
74.0
2
206.3
3
345.9
4
38.3
1
71.0
2
155.9
3
350.6
4
44.0
1
74.9
2
156.8
3
358.0
4

As can be seen from Fig. 13, the recognition rate of the test set sample reached 97.5 %, and fault types can be effectively identified by SVM. Therefore, the proposed method has achieved good diagnostic results. For the diagnosis of other diesel engines, a large number of confirmatory experiments are needed to make further research.

Fig. 13. The training and test results of SVM

The training and test results of SVM

a) Accuracy rate 100 %

The training and test results of SVM

b) Accuracy rate 97,5 %

6. Conclusions

Aiming at the problem that VMD’s modal number K and penalty factor α are difficult to be determined, an improved VMD algorithm based on adaptive genetic algorithm is proposed, which reduces the interference of human factors to decomposition results. The simulation signal and the fault diagnosis of diesel engine crankshaft are studied. The conclusions are as follows:

(1) In this paper, the crossover and mutation operators of traditional adaptive genetic algorithm are optimized, so that individuals maintain good diversity. Then the precocious phenomenon is avoided, and the problem of slow evolution is solved. VMD is optimized by IAGA, and the optimal parameter combination is obtained by the evolution of the population.

(2) As is known to all, mechanical fault diagnosis of diesel engine is difficult. According to the proposed method, the vibration signals under different wear conditions are studied, and the characteristic component of the wear fault of the diesel engine crankshaft bearing and its energy are effectively extracted. The wear fault of the diesel engine crankshaft bearing can be effectively identified by SVM, which provides a theoretical basis for the condition monitoring and fault diagnosis of the diesel engine crankshaft bearing wear.

References

  1. Rauber T. W., Boldt F. D. A., et al. Heterogeneous feature models and feature selection applied to bearing fault diagnosis. IEEE Transactions on Industrial Electronics, Vol. 62, Issue 1, 2015, p. 637-646. [Publisher]
  2. Xie Z. J., Song B. Y., Zhang Y., et al. Application of an improved wavelet threshold denoising method for vibration signal processing. Advanced Materials Research, Vol. 889, Issue 890, 2014, p. 799-806. [Publisher]
  3. Huang N. E., Shen Z., Long S. R. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proceedings Mathematical Physical and Engineering Sciences, Vol. 454, Issue 1971, 1998, p. 903-995. [Publisher]
  4. Huang N. E., Shen Z., Long S. R. A new view of non- linear water waves: the Hilbert spectrum. Annual Review of Fluid Mechanics, Vol. 31, Issue 1, 2003, p. 417-457. [Publisher]
  5. Yu D., Cheng J., Yang Y. Application of EMD method and Hilbert spectrum to the fault diagnosis of roller bearings. Mechanical Systems and Signal Processing, Vol. 19, Issue 2, 2005, p. 259-270. [Publisher]
  6. Cheng J., Yu D., Yu Y. The application of energy operator demodulation approach based on EMD in machinery fault diagnosis. Mechanical Systems and Signal Processing, Vol. 21, Issue 2, 2007, p. 668-677. [Publisher]
  7. Wang T., Zhang M., Yu Q., et al. Comparing the application of EMD and EEMD on time-frequency analysis of seismic signal. Journal of Applied Geophysics, Vol. 83, Issue 6, 2012, p. 29-34. [Publisher]
  8. Wu Z. H., Huang N. E. Ensemble empirical mode decomposition: a noise-assisted data analysis method. Advances in Adaptive Data Analysis, Vol. 1, Issue 1, 2009, p. 1-41. [Publisher]
  9. Wang H., Chen J., Dong G. Feature extraction of rolling bearing’s early weak fault based on EEMD and tunable Q-factor wavelet transform. Mechanical Systems and Signal Processing, Vol. 48, Issue 2, 2014, p. 103-119. [Publisher]
  10. Torres M. E., Colominas M. A., Schlotthauer G., et al. A complete ensemble empirical mode decomposition with adaptive noise. IEEE International Conference on Acoustics, Speech and Signal Processing, 2011, p. 4144-4147. [CrossRef]
  11. Smith J. S. The local mean decomposition and its application to EEG perception data. Journal of the Royal Society Interface, Vol. 2, Issue 5, 2005, p. 443-454. [Publisher]
  12. Cheng J. S., Zheng J. D., Yu Y. A nonstationary signal analysis approach-the local characteristic-scale decomposition method. Journal of Vibration Engineering, Vol. 25, Issue 2, 2012, p. 215-220. [CrossRef]
  13. Frei Mark G. Intrinsic time-scale decomposition: time-frequency-energy analysis and real-time filtering of non-stationary signals. Proceedings Mathematical Physical and Engineering Sciences, Vol. 463, Issue 2078, 2007, p. 321-342. [Publisher]
  14. Dragomiretskiy K., Zosso D. Variational mode decomposition. IEEE Transactions on Signal Processing, Vol. 62, Issue 3, 2013, p. 531-544. [Publisher]
  15. 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]
  16. An X. L., Tang Y. J. Denoising of hydropower unit vibration signal based on variational mode decomposition and approximate entropy. Transactions of the Institute of Measurement and Control, Vol. 38, Issue 3, 2016, p. 282-292. [Publisher]
  17. Li Z. P., Chen J. L., Zi Y. Y. Independence-oriented VMD to identify fault feature for wheel set bearing fault diagnosis of high speed locomotive. Mechanical Systems and Signal Processing, Vol. 85, 2017, p. 512-529. [Publisher]
  18. E. J. W., Bao Y. L., Ye J. M. Crude oil price analysis and forecasting based on variational mode decomposition and independent component analysis. Physica A, Vol. 484, 2017, p. 412-427. [Publisher]
  19. An X. L., Tang Y. J. Application of variational mode decomposition energy distribution to bearing fault diagnosis in a wind turbine. Transactions of the Institute of Measurement and Control, Vol. 39, Issue 7, 2017, p. 1000-1006. [Publisher]
  20. Pająk M. Genetic-fuzzy system of power units maintenance schedules generation. Journal of Intelligent and Fuzzy Systems, Vol. 28, Issue 4, 2015, p. 1577-1589. [CrossRef]
  21. Srinivas M., Patnaik L. M. Adaptive probabilities of crossover and mutation in genetic algorithms. IEEE Transactions on Systems Man and Cybernetics, Vol. 24, Issue 4, 2002, p. 656-667. [Publisher]
  22. Bandt C., Pompe B. Permutation entropy: a natural complexity measure for time series. Physical Review Letters, Vol. 88, Issue 17, 2002, p. 174102. [Publisher]
  23. Yan R., Liu Y., Gao R. X. Permutation entropy: A nonlinear statistical measure for status characterization of rotary machines. Mechanical Systems and Signal Processing, Vol. 29, Issue 5, 2012, p. 474-484. [Publisher]