Next Article in Journal
Performance Analyses of Energy Detection Based on Square-Law Combining in MIMO-OFDM Cognitive Radio Networks
Next Article in Special Issue
Multidimensional Data Interpretation of Vibration Signals Registered in Different Locations for System Condition Monitoring of a Three-Stage Gear Transmission Operating under Difficult Conditions
Previous Article in Journal
Gait Disorder Detection and Classification Method Using Inertia Measurement Unit for Augmented Feedback Training in Wearable Devices
Previous Article in Special Issue
An Expert System for Rotating Machine Fault Detection Using Vibration Signal Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rotating Machinery Diagnosing in Non-Stationary Conditions with Empirical Mode Decomposition-Based Wavelet Leaders Multifractal Spectra

Department of Mechanical Engineering, Kazimierz Pulaski University of Technology and Humanities in Radom, Malczewskiego 29, 26-600 Radom, Poland
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(22), 7677; https://doi.org/10.3390/s21227677
Submission received: 18 October 2021 / Revised: 10 November 2021 / Accepted: 16 November 2021 / Published: 18 November 2021

Abstract

:
Diagnosing the condition of rotating machines by non-invasive methods is based on the analysis of dynamic signals from sensors mounted on the machine—such as vibration, velocity, or acceleration sensors; torque meters; force sensors; pressure sensors; etc. The article presents a new method combining the empirical mode decomposition algorithm with wavelet leader multifractal formalism applied to diagnosing damages of rotating machines in non-stationary conditions. The development of damage causes an increase in the level of multifractality of the signal. The multifractal spectrum obtained as a result of the algorithm changes its shape. Diagnosis is based on the classification of the features of this spectrum. The method is effective in relation to faults causing impulse responses in the dynamic signal registered by the sensors. The method has been illustrated with examples of vibration signals of rotating machines recorded on a laboratory stand, as well as on real objects.

1. Introduction

When properly processed, dynamic signals recorded during the operation of rotating machines are a valuable source of information about the condition of machines. Algorithms are still being developed, which will quickly, unequivocally, and automatically determine the operational state of the machine. One of the most commonly used methods in diagnostics are vibration sensors due to the ease of assembly, the possibility of obtaining information without disassembling the machine, and the wealth of information obtained about the condition of the machine in various frequency ranges. In the case of machine operation with constant speed and load, the most effective method of analysis is the frequency spectrum [1]. However, due to the fluctuating rotational speed of the machine, often due to a variable load, the frequency spectrum becomes blurred. The real challenge is to diagnose machines in non-stationary conditions. In order to avoid resampling in the analysis of dynamic signals, time–frequency methods are used [2]. The authors present a systematic review of over 20 time–frequency methods used to detect machine damage. One of the simplest is the short-time Fourier transform STFT [3]. The improved version of the spectrogram is the Wigner–Ville analysis, the use of which in diagnosing damage to gears is described in [4,5]. The use of wavelet analysis in diagnosing damage to rotating machines can be found, among others in [6,7]. Other papers [8,9,10] present a method of detecting and diagnosing gear damage on the basis of a reference model based on the signal averaging technique. The proposed algorithm first establishes an autoregressive (AR) model of the gear vibration signal in its base state, and then diagnoses the state based on the residual signal. The use of the synchronous method (order-tracking analysis) allows to synchronize the diagnostic signal with the rotation of the diagnosed machine, as well as to eliminate the influence of background noise from other sources [11,12]. The article [13] proposes a method of informative selection of frequency bands. It utilizes the approach of non-negative matrix factorization applied to time–frequency signal representation. For the analysis of non-stationary signals, the empirical mode decomposition EMD algorithm is also used, which decomposes the signal into modes in an empirical manner [14,15]. Since the classic method can mix two or three modes, various modifications of the EMD method are developed [16,17,18].
In recent years, more and more often, efforts are being made to automate the diagnosis process, and qualified diagnosticians are to be replaced by diagnostic algorithms. So-called data-driven statistical methods can also be successfully used in diagnosing machines in non-stationary conditions. Determining the characteristics of a dynamic signal and its classification is carried out by machine learning algorithms [19,20,21,22]. More and more papers are devoted to the diagnostics of machines with the methods of deep learning [23,24,25,26,27,28,29]. Methods that allow detecting failures at an early stage, especially in non-stationary conditions that often occur in real objects, are still being searched and developed.
The fractal theory has unique advantages in dealing with transient, nonlinear, and non-stationary signals. It is a new trend in solving practical problems involving the study of short-term signals in a transient process with strong nonlinearity, based on the fractal theory in order to diagnose damage to mechanical equipment. The impact of damage to the outer race of a rolling bearing on the width of the multifractal spectrum is presented in [30]. Detrended fluctuation analysis MF-DFA is a commonly used algorithm for multifractal analysis. The method is often used to diagnose damage to rolling bearings [31,32,33] and gears [34,35,36]. In [37], the MF-DFA algorithm was used for the analysis of frictional vibrations, where the ensemble empirical mode decomposition EEMD was used to denoise the signal. The adaptive MF-DFA algorithm is presented in [38]. The disadvantage of the MF-DFA algorithm is its sensitivity to analysis parameters. These disadvantages are not present in the method of wavelet leaders (WLMF). The wavelet leaders algorithm was used to analyze the images [39] and diagnose damage to rotating machines [40].
The article proposes a novel EMD-WLMF method consisting in the multifractal analysis of the first empirical component (IMF1) of the vibration signal using WLMF algorithm. As a result, a multifractal spectrum is obtained, the parameters of which are used to distinguish the state of the machine. Since the tested damage causes impulse disturbances of the vibration signal, the signal is first filtered out using the empirical decomposition of the signal EMD. Only the first intrinsic mode function (IMF1) containing high-frequency information is subject to multifractal analysis. The method has been illustrated with examples of signals recorded on a laboratory stand as well as on real objects. The presented examples, especially the analysis of signals recorded in the vehicle, show that the combination of EMD and WLMF methods allows the detection of damage at an early stage in conditions unfavorable to diagnosis—i.e., at floating rotational speed and load.
The article is organized as follows. Section 2 describes the methods of signal analysis used in the study, i.e., the WLMF algorithm, as well as determining the first intrinsic mode function (IMF1) using the EMD method. In Section 3, the described algorithm for fault diagnosis was used on the basis of signals recorded on a laboratory stand and on real objects under operating conditions. The discussion is presented in Section 4.

2. Materials and Methods

2.1. Multifractal Formalism

Multifractal analysis, which is based on estimated signal scaling exponents, is a popular statistical tool for studying time series. The mathematical formalism is based on the increments of their values, the measure of which are the pointwise Holder exponents h of the time function x(t) at the point t0, determined by the supremum of all exponents satisfying, for the constant C > 0, the condition
| x ( t ) P n ( t t 0 ) | C | t t 0 | h
where P n ( t t 0 ) is a polynomial of the order n < h. The pointwise Holder exponent describes the regularity of the function. The greater the local regularity of the time series, the higher the singularity exponent. The multifractal spectrum D(h) represents the histogram of the Holder exponents and determines the fractal dimensions of the subsets of singularities with a given exponent value [41,42,43].
In the time–frequency approach to signal analysis, the WTMM (wavelet transform modulus maxima) method was initially used, which is based on a continuous wavelet transform (CWT). It consists in determining the chains of local, the structure function and its scaling exponents, as well as the Legendre transform. Practical implementations of such an algorithm showed a number of faults that make it impossible to conduct research for some types of real signals [44].
Another approach to the problem of estimating local scaling exponents as a method of studying the regularity of time series and their multifractality is related to the multifractal detrended fluctuation analysis (MF-DFA) method. MF-DFA enables to study the observed signals in terms of their multifractality, provides a more stable approach to multifractal formalism than the WTMM method [31,36,45,46,47].
The formalism in the time–frequency domain used in the paper allows for the estimation of multifractal parameters using wavelet leaders (WLMF), which are representatives of local Holder exponents of the signal. In the case of wavelet coefficients centered around zero values, it is difficult to guarantee numerical stability. This problem does not arise in the case of the WLMF method, based on the wavelet coefficients obtained as a result of the discrete wavelet transform (DWT). The next steps of the algorithm include the selection of coefficients called wavelet leaders, determination of the structure function with scaling exponents, and the multifractal spectrum. The algorithm shows low computational costs, numerical stability, and high versatility in terms of real signals [39,40].
The scaling exponents of the structural function do not depend on the selection of the wavelet, provided that the number of zero wavelet moments is two times greater than the largest exponent of the signal holder. The mother wavelet must be orthogonal or biorthogonal. In the work, the Daubechies (db2) wavelet mother was arbitrarily selected.
For the coefficients (2) of the discrete wavelet transform (DWT) of the function x(t) and the basic wavelet with a compact support ψ0 (t)
d x ( j , k ) = R x ( t ) 2 j ψ 0 ( 2 j t k ) d t
wavelet leaders (3), for the set of the largest coefficients d x ( j , k ) d λ in the vicinity of 3λ, are defined by the dependence on any scale
L x ( j , k ) = s u p λ 3 λ | d λ |
where j , k are integers and 3 λ = 3 λ j , k = λ j , k 1   λ j , k   λ j , k + 1 and λ = λ j , k = [ k 2 j , ( k + 1 ) 2 j ] . L x ( j , k ) consists of the largest wavelet coefficient d x ( j , k ) computed at all finer scales 2 j 2 j within a narrow time neighborhood ( k 1 ) · 2 j 2 j k < ( k + 2 ) · 2 j .
It can be shown that the Holder exponents are the scaling exponents of wavelet leaders: L x ( j , k ) ~ 2 j h . Moreover, the structure function (4) defined for wavelet leaders is described by the power dependence, the exponent of which is the multifractal scaling exponent ζ ( q ) : R R .
Z L ( q , j ) = 1 n j k = 1 n j L x ( j , k ) q = E L x ( j , k ) q ~ 2 j ζ ( q )
where q is the order of the structure function and n j is the number of intervals of the multiresolution analysis.
The function obtained by the Legendre transformation of the multifractal scaling exponent ζ ( q ) , under mild conditions of signal regularity, is the upper limit for the multifractal spectrum (5) of the tested signal
D ( h ) m i n q 0 [ 1 + q h ζ ( q ) ]
The description of dynamic properties of systems is successfully carried out on the basis of the parameters of multifractal spectra of representative time series.
The following were selected as parameters of the multifractal spectrum D ( h ) related to singularities h representing the local scaling of the measure in different places of the time series:
  • multifractality level, representing the heterogeneity of the signal under study, = h m a x h m i n , where h m a x and h m i n are the singularities corresponding to the largest and the smallest fluctuations in the time series (observed signal);
  • span of dimensions of subsets of singularities D = D ( h m a x ) D ( h m i n ) ;
  • the singularity with the greatest dimension, which is the most common singularity of the time series { h 0 : D ( h 0 ) = max D ( h ) } .
The procedure of the method is shown in Figure 1.
The effect of the algorithm is illustrated in Figure 2 for a harmonic signal and various types of noise, examined in terms of multifractality. Figure 2a shows sample waveforms of signals, Figure 2b shows their multifractal spectra with the characteristic points. The harmonic signal, three types of random signals varying in the level of multifractality, and the sinusoidal signal disturbed by the multifractal signal were analyzed. The level of multifractality can be measured by the width of the multifractal spectrum. Thus, both the harmonic signal and the white noise cannot be treated as multifractal signals. Multifractal signals are characterized by impulse disturbances, such as harmonic signals with noise.
Figure 2b shows the characteristic points of the multifractal spectrum, which were used to define the characteristics of the spectrum:
  • Spectrum width
w i d t h = = h ( q = 5 ) h ( q = 5 )
  • Spectrum asymmetry
a s y m m e t r y = D = D ( q = 5 ) D ( q = 5 )
  • The singularity exponent with the greatest fractal dimension
s h i f t = h 0 = h ( q = 0 )
Figure 2c shows the multifractal spectra determined when selecting the Daubechies mother wavelet of orders 1 to 6. As the selection of the wavelet order has little effect on the parameters of the spectrum, the Daubechies wavelet of order 2 (‘db2’) was arbitrarily selected. The selection of Daubechies wavelets of different orders should yield consistent results when used repeatedly.

2.2. Empirical Mode Decomposition

The EMD is an iterative numerical approximation algorithm designed to extract intristic mode functions IMFs from signals adaptively by cubic spline interpolation according to the local characteristic time-scale. It decomposes a signal into IMFs via iterative sifting.
Any signal can be approximated by a superposition of a series of IMFs as written by
x ( t ) = i = 1 n c i ( t ) + r n ( t ) ,
where c i ( t ) is the i-th IMF, and r n ( t ) is the residual signal which represents the slowly varying or constant trend of the signal.
Applying the Hilbert transform [14] to each IMF in Equation (9), we can construct the corresponding analytic signal. Then we can compute the instantaneous frequency via the derivative of instantaneous phase relative to time. By expressing the analytic signals in polar coordinate form, and taking the real part, we can obtain the Hilbert amplitude spectrum
T F R x ( t , f ) = R e i = 1 n { a i ( t ) e x p | j   2 π f i ( t ) d t | }
and Hilbert energy spectrum
T F R x ( t , f ) = i = 1 n a i 2 ( t ) δ [ f f i ( t ) ]
where δ (·) is the Dirac delta function.
The EMD is the core of the Hilbert–Huang transform. For a real signal x ( t ) , the EMD procedure is as follows.
1.
Initialize parameters: Set iteration index i = 1, residual signal r 0 ( t ) = x ( t ) .
2.
Extract the i-th IMF:
a.
Let j = 0, and h i j ( t ) = h i 1 ( t ) .
b.
Find the local minima and the local maxima of h i j ( t ) .
c.
Interpolate the local minima and the local maxima with cubic spline to construct the lower and the upper envelopes of h i j ( t ) .
d.
Compute the instantaneous mean of the lower and upper envelopes m i j ( t )
e.
Let h i j ( t ) = h i j ( t ) m i j ( t ) .
f.
If h i j ( t ) satisfies the stop criteria for IMF sifting, then set the i-th IMF c i ( t ) = h i j ( t ) . Otherwise, let j = j + 1, return to step 2b.
3.
Let r i ( t ) = r i 1 ( t ) c i ( t ) .
4.
If r i ( t )   satisfies the stop criteria for EMD, then set r i ( t )   as the residue, and terminate the EMD process. Otherwise, let i = i + 1 , return to Step 2.
To illustrate the method, a simulated signal was used, which shows the vibrations of a gear transmission operating at a variable rotational speed [9].
x ( t ) = m = 0 M A m [ 1 + a ˜ m ( t ) ] cos { 2 π f m ( t ) + β m + θ ˜ m ( t ) } + z ( t ) + ϕ ( t )
where m (0,1,…,M) is the meshing harmonic number, A m the amplitude at the m-th harmonic frequency f m (i.e., f m = m · N · f s ( t ) , where N is the tooth number and f s is the shaft rotation frequency), t time vector (with sampling time d t ), β m the initial phase, z ( t )   the impact-induced resonant vibration, a ˜ m ( t ) and θ ˜ m ( t ) are the modified amplitude and phase modulation functions at the m-th harmonic, respectively, ϕ ( t )   the white noise.
The following model parameters were adopted for the simulation: M = 2, A 1 = 5 , A 2 = 1.25 , N = 10, β m = 0 , dt = 0.01 s; a ˜ 1 ( t ) = 0.1 sin ( 2 π f s · t ) , a ˜ 2 ( t ) = 0.025 sin ( 4 π f s · t ) ; θ ˜ 1 ( t ) = 4.7 sin ( 2 π f s · t ) , θ ˜ 2 ( t ) = 1.2 sin ( 4 π f s · t ) . f s ( t ) is a random walk around the value 1. z ( t ) is described as a convolution of the resonant response to impact during one rotation and the pulse train z ( t ) = h ( t ) g ( t ) , where h ( t ) = 30 sin ( 40 π · t ) · e x p ( 10 t ) and g ( t ) is the sequence of pulses occurring once per revolution.
Signal-to-noise-ratio, SNR, is calculated as
S N R = 10   log 10 i = 1 n ( s i ) 2 i = 1 n ( w i ) 2
where s i is i-th amplitude of the signal and w i i-th amplitude of the noise. For a simulated signal, SNR = 4.6 dB. The main components of signal power are visible impacts.
The signal includes both quasi-harmonic components and impulse responses related to the damage of the gear. The signal presented by the Formula (12) is shown in Figure 3a, together with the following intristic mode functions from 1 to 5, which are the result of decomposition. Figure 3b shows the multifractal spectra for waveforms illustrated in Figure 3a.
The first empirical mode reveals the high-frequency impacts. This high-frequency mode increases in energy as the wear or damage progresses. The next modes show the resonance in the vibration signal and rotating frequency. The first intristic mode function will be used for further analysis. The spectrum for IMF1 is broader and more regular than the spectrum of the raw signal and the other IMF, while the IMF2 is a monofractal.

2.3. EMD-WLMF Method

Figure 4 shows a diagram of the EMD-WLMF method. From a dynamic signal using EMD analysis is extracting the first IMF. Then, using the WLMF algorithm, the signal parameters h and D(h) are determined, on the basis of which a multifractal spectrum can be drawn. Based on the characteristics of this spectrum, the operational state of the rotating machine is classified.

3. Application of the EMD-WLMF Algorithm in Diagnostics of Rotating Machines

The method presented in Section 2.3 was tested on the vibration signals recorded for the gear transmission on a laboratory stand and on a real object, which is a passenger car while driving.

3.1. Gear Transmission Vibration Signal Analysis on a Laboratory Stand

Measurements were carried out on a demonstration stand (Figure 5). The influence of assembly errors and gear teeth wear on vibrations was investigated. The electric motor enables speed control in the range of 100–3000 rpm (without load). The load is pressure regulated with an overflow valve up to 5 MPa. Vibration acceleration was measured with an accelerometer bolted to the gear bearing housing in the vertical direction. The optimal backlash was set to 0.1 mm.
Compared to previous papers [48], the experiment was extended with additional intermediate operating states and the analysis method was modified.
Measurements were carried out for the following five states:
  • fault-free (new gears, the optimal backlash, parallel shaft location);
  • new gears and misalignment by an angle up to 1/3° (two cases);
  • new gears and increased backlash +0.1 mm up to +0.3 mm (three cases);
  • worn teeth (three cases);
  • worn teeth and increased backlash +0.3 mm (two cases).
Vibration acceleration signals were recorded for a rotational speed of about 1050 rpm and a load of 12%—pressure 0.6 MPa. Each sample included a time series with a length of n = 10,000, recorded at a sampling frequency of 10 kHz.
The representative waveforms of the vibration acceleration signal for each of the five states are shown in Figure 6. The first IMF is marked on them, which is taken into account in further analysis of the signal. The signals presented in the graph represent a significant degree of damage, and their characteristic features are visible.
Twenty records for each condition were used for further analysis, each 1 s long (10,000 samples) for a total of 220 records. Each record was subjected to the EMD-WLMF analysis. Figure 7 shows the multifractal spectra for various operational states of the gears. The examples show up to two cases per state to increase the visibility of the image.
The least multifractal character is for the fault-free state. As the damage develops, the width of the spectrum increases. The spectral peak is also shifted. For the base state, the spectrum asymmetry index is close to zero. For the tested damages, the absolute value of the asymmetry index increases, and for the wear condition of the gears it has the opposite sign than for the increased backlash. For mixed states of damage, the spectrum is characterized by a large width and indices of asymmetry varying within wide limits.
Three parameters were selected for the diagnosis of the gear condition: the width of the multifractal spectrum, the asymmetry index, and the singularity of the greatest dimension. The list of these parameters is presented in the scatter plot (Figure 8a).
As the damage develops, the level of multifractality of the signal increases and the spectrum is wider. The asymmetry index for the ‘base state’ is close to zero, and not much higher for the ‘misalignment’ state. The states of increased backlash and wear are characterized by the opposite sign of the spectrum asymmetry index and an increase in its absolute value. The spectrum shift indicator also changes for different operating states.
Any known machine learning method—e.g., k-nearest neighbours (k-NN), support vector machines (SVM), or neural network (NN)—can be used to classify the damage [25]. The fault classification was performed for a newly generated test set of data (240 records) with support vector machines (SVM). The confusion matrix is shown in Figure 8b. The average true positive rate for the test set is 95.16%. The faulty recognition of the state was for ‘fault-free’ and initial ‘misalignment’. The remaining states found to be in error are ‘increased backlash’ and ‘worn teeth and increased backlash’. The tests performed on the laboratory stand brought the desired effect. The influence of damage development on the parameters of the multifractal spectrum was presented. However, the stand lacks additional disturbances that accompany the operation of real machines.

3.2. Analysis of the Transmission Vibration Signal in a Car

The tests were carried out on a Punto five-speed gearbox (Figure 9a). The experiment was carried out during road tests for various rotational speeds and loads. The experiment used the recorded signals of acceleration of vibrations of the gearbox housing and additional signals from the crankshaft position and throttle position sensor.
The vibration accelerations of the gearbox housing were measured with a Bruel & Kjaer type IEPE No. 4514 sensor. The signals were recorded with a Bruel & Kjaer PULSE type 3560E portable data recorder with a sampling frequency of 65,536 Hz.
Signals of a duration of about 1 min were recorded while driving at the most preset speed possible. Additional signals enabled signal synchronization and measurement of rotational speed and load.
The active experiment consisted in registering signals under simulated mechanical damage to the gearbox, mimicking those often diagnosed in car repair shops. Such damage includes, among others, wear of the gearbox teeth.
In order to investigate the influence of gearbox teeth wear on the vibration signal, an active experiment was performed for the following conditions:
  • gearbox in good condition
  • fifth gear drive gear teeth worn at about one-third of the circumference
  • fifth gear drive gear teeth heavily worn at about one-third of the circumference (Figure 9b)
  • gearbox after replacing worn wheels with new ones.
Figure 10a–d show the first IMF waveforms of the gearbox housing vibration signal for four operating states, and Figure 10e shows examples of changes in rotational speed during the record registration.
The EMD-WLMF algorithm was tested for 80 records, 20 records for each state. Figure 11a presents sample multifractal spectra for four states, and Figure 11b presents multifractal spectra for the same states using the only WLMF algorithm.
By comparing the multifractal spectra for the first IMF of the vibration signal and the raw of the vibration signal, it can be concluded that the introduction of the EMD enables the identification of tooth wear at an early stage. In Figure 10b, the spectrum at initial wear almost does not differ from the base state and the spectrum for a new pair of gears. The EMD-WLMF analysis showed that the spectrum for the new pair of gears (‘after repair’) is the most symmetrical. The spectrum for the ‘initial state’ is narrower, but shows a shift and asymmetry. As the wear develops, the spectrum becomes wider, shifted along the h(q) axis, and shows asymmetry. The illustration of spectrum features for all 80 records is shown in Figure 12a.
The scatter plot for more records confirms that the spectrum width increases as the wear develops. The asymmetry coefficient for all states has a large spread, while for the new pair of teeth (‘after repair’) the spectrum is the most symmetrical. As the wear of the teeth increases, the spectrum shifts to the right. The average trend of damage development was marked. For a new pair of teeth (‘after repair’), the development of the damage may follow a ‘different path’ on the scatter plot, but the increasing trend of the width indicators and spectrum shift should be analogous.
The fault classification was performed for a newly generated test set of data (120 records) with support vector machines (SVM). The confusion matrix is shown in Figure 12b. The average true positive rate for the test set is 94.17%. The faulty recognition of the state was for ‘fault-free’ and ‘initial wear’. Due to the dispersion of the data, the diagnosis can be confirmed on the basis of the average trend shown in Figure 12a.

3.3. Analysis of the Vibration Signal of the Internal Combustion Engine Head

The next example shows the analysis of the vibration signal of the engine head signal recorded under operating conditions.
The tests were carried out on the drive system of the Fiat Punto car, with a four-cylinder 1.2 spark-ignition engine. The experiment was carried out during road tests for various rotational speeds and loads. The following dynamic signals were recorded:
  • acceleration of vibrations of the cylinder head at the first cylinder in the vertical and horizontal directions
  • acceleration of vibrations of the cylinder head at the fourth cylinder in the vertical direction
and additional signals
  • from the crankshaft position sensor
  • ignition in the first cylinder
  • throttle position.
Vibration accelerations of the cylinder head were measured with Bruel & Kjaer DeltaShear type 4393 sensors, mounted with a threaded connection. The signals were recorded using a Bruel & Kjaer PULSE type 3560E portable data recorder with a sampling frequency of 65,536 Hz.
Signals of approximately 1 min duration were recorded while driving at a constant speed. Additional signals made it possible to identify engine cycles, injection and ignition timing, and valve timing.
A leak in the piston-cylinder system may be caused by burnout of the exhaust valve. During the tests, such damage was simulated by cutting the valve plug at a length of approx. 3 mm (defect 1) and 6 mm (defect 2)—Figure 13.
The process of generating vibrations in an internal combustion engine is very complex. The measured vibrations are a combination of periodic waves related to the operation of rotating elements and responses to impulse excitations related to reciprocating motion, as well as excitations caused by gas pressure. Strong transients in the vibroacoustic signal come from the work of the intake and exhaust valves, injectors, the combustion process, and piston strokes against the cylinder sleeve. Additional impulse excitation can be a response to mechanical failures not detected by the on-board diagnostic (OBD) system [49,50].
The time series of vibration accelerations after determining the first IMF for four engine operation cycles and three operating states are shown in Figure 14.
Figure 15 presents exemplary multifractal spectra for three operating states using the WLMF and EMD-WLMF algorithms. Due to the presence of transients in the signal, a large share of responses to impulse excitations—such as valve closing, ignition, piston hitting the cylinder, etc.—of the multifractal spectrum after EMD application does not differ much from the spectrum without EMD filtration. As a result of the damage development, the spectrum becomes wider and its maximum value for q = 0 is shifted to the right.
Figure 16 shows the scatter plot for the three operational states, each with 20 records. As the damage develops, the values of the ‘width’ and ‘shift’ parameters increase. The ‘asymmetry’ parameter is characterized by the largest dispersion of values, although it also increases taking into account the average values for the entire set. The average trend for damage development is marked.
The fault classification was performed for a newly generated test set of data (120 records) with support vector machines (SVM). The confusion matrix is shown in Figure 16b. The average true positive rate for the test set is 92.5%. The faulty recognition of the state was for ‘fault-free’ and ‘defect 1’. Due to the dispersion of the data, the diagnosis can be confirmed on the basis of the average trend shown in Figure 16a.

4. Discussion

The article uses the method of wavelet leaders’ multifractal analysis dedicated to the first IMF of the vibration signal of rotating machines. The presented method can be applied to dynamic signals in non-stationary conditions due to the use of the EMD and wavelet analysis methods. Previous works [32] used the EMD analysis and the MF-DFA algorithm to diagnose the damage. However, the WLMF algorithm shows lower computational costs, numerical stability, and high versatility in terms of real signals compared to MF-DFA. The method corresponds to the need to find a relatively fast algorithm for diagnosing the condition of machines in non-stationary conditions. It is a statistical method. It can be an alternative to methods based on a set of signal features. Failure classification can be done using one of the machine learning methods, such as SVM, k-NN, or a neural network based on the features of the multifractal spectrum.
The method has been tested for signals recorded on a laboratory stand and for real objects. Damage was distinguished based on three features of the multifractal spectrum, named here: width, asymmetry, and shift. The tests performed on the laboratory stand brought a good effect. Multifractal analysis of only the first IMF allowed to filter out harmonics that do not show multifractality features. The development of damage to the gear transmission was reflected in the examined features of the multifractal spectrum. However, the stand lacks additional disturbances that accompany the operation of real machines.
The next two investigated cases are the diagnosis of mechanical damage to the vehicle while driving, such as gearbox wear and burnout of the exhaust valves of the internal combustion engine. The signals were recorded in non-stationary conditions. Due to the fluctuations in rotational speed and load, which could affect the dynamics of the signal, the dispersion of the features of the multifractal spectrum was much greater than that observed in the tests on the laboratory stand. Due to the dispersion of the features of the multifractal spectrum, the average value of these features can be used to diagnose the damage. By comparing the multifractal spectra for the first IMF of the vibration signal and the raw of the vibration signal of the gearbox, it can be concluded that the introduction of the EMD enables the identification of tooth wear at an early stage. Due to the presence of transients of the engine head vibrations signal, a large share of responses to impulse excitations—such as valve closing, ignition, piston hitting the cylinder, etc.—the multifractal spectrum after EMD application does not differ much from the spectrum without EMD filtration.
Like any method, this one also has its limitations. The main limitation is the multifractality of the signal, so the signal must contain impulse excitations caused by damage, as in a toothed gear. If the damage causes an increase in the harmonic component, the change in the multifractal spectrum will be hardly noticeable and then other methods—e.g., order tracking or the AR reference model—will bring better results.

Author Contributions

Conceptualization, I.K. and A.P.; Methodology, A.P.; Software, I.K.; Validation, I.K. and A.P.; Formal analysis, A.P.; Investigation, I.K.; Writing—original draft preparation, I.K.; Writing—review and editing, I.K.; Visualization, I.K.; Supervision, A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Randall, R.B. Frequency Analysis; Bruel & Kjaer: Naerum, Denmark, 1987. [Google Scholar]
  2. Feng, Z.; Liang, M.; Chu, F. Recent advances in time–frequency analysis methods for machinery fault diagnosis: A review with application examples. Mech. Syst. Signal Process. 2013, 38, 165–205. [Google Scholar] [CrossRef]
  3. Manhertz, G.; Bereczky, A. STFT spectrogram based hybrid evaluation method for rotating machine transient vibration analysis. Mech. Syst. Signal Process. 2021, 154, 107583. [Google Scholar] [CrossRef]
  4. Baydar, N.; Ball, A.A. Comparative Study of Acoustic and Vibration Signals in Detection of Gear Failures Using Winger–Ville Distribution. Mech. Syst. Signal Process. 2001, 15, 1091–1107. [Google Scholar] [CrossRef]
  5. Wilk, A.B.; Madej, H.M.; Łazarz, B.E. Vibration Processing Techniques for Fault Detection in Gearboxes. In Proceedings of the ASME 2003 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Chicago, IL, USA, 2–6 September 2003; pp. 657–664. [Google Scholar] [CrossRef]
  6. Yan, R.; Gao, R.X.; Chen, X. Wavelets for fault diagnosis of rotary machines: A review with applications. Signal Process. 2014, 96, 1–15. [Google Scholar] [CrossRef]
  7. He, S.; Chen, J.; Zhou, Z.; Zi, Y.; Wang, Y.; Wang, X. Multifractal entropy based adaptive multiwavelet construction and its application for mechanical compound-fault diagnosis. Mech. Syst. Signal Process. 2016, 76–77, 742–758. [Google Scholar] [CrossRef]
  8. Endo, H.; Randall, R.B. Enhancement of autoregressive model based gear tooth fault detection technique by the use of minimum entropy deconvolution filter. Mech. Syst. Signal Process. 2007, 21, 906–919. [Google Scholar] [CrossRef]
  9. Wang, W.; Wong, A.K. Autoregressive Model-Based Gear Fault Diagnosis. ASME J. Vib. Acoust. 2002, 124, 172–179. [Google Scholar] [CrossRef]
  10. Komorska, I. A Vibroacoustic Diagnostic System as an Element Improving Road Transport Safety. Int. J. Occup. Saf. Ergon. 2013, 19, 371–385. [Google Scholar] [CrossRef] [PubMed]
  11. Wu, J.-D.; Bai, M.R.; Su, F.-C.; Huang, C.-W. An expert system for the diagnosis of faults in rotating machinery using adaptive order-tracking algorithm. Expert Syst. Appl. 2009, 36, 5424–5431. [Google Scholar] [CrossRef]
  12. Pawlik, P. The Use of the Acoustic Signal to Diagnose Machines Operated Under Variable Load. Arch. Acoust. 2020, 45, 263–270. [Google Scholar] [CrossRef]
  13. Wodecki, J.; Kruczek, P.; Bartkowiak, A.; Zimroz, R.; Wyłomańska, A. Novel method of informative frequency band selection for vibration signal using Nonnegative Matrix Factorization of spectrogram matrix. Mech. Syst. Signal Process. 2019, 130, 585–596. [Google Scholar] [CrossRef]
  14. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. London. Ser. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  15. Sharpley, R.; Vatchev, V. Analysis of the Intrinsic Mode Functions. Constr. Approx. 2006, 24, 17–47. [Google Scholar] [CrossRef]
  16. Moore, K.J.; Kurt, M.; Eriten, M.; McFarland, D.M.; Bergman, L.A.; Vakakis, A.F. Wavelet-bounded empirical mode decomposition for measured time series analysis. Mech. Syst. Signal Process. 2018, 99, 14–29. [Google Scholar] [CrossRef]
  17. Ge, H.; Chen, G.; Yu, H.; Chen, H.; An, F. Theoretical Analysis of Empirical Mode Decomposition. Symmetry 2018, 10, 623. [Google Scholar] [CrossRef] [Green Version]
  18. Sharma, V.; Anand, P. Extraction of weak fault transients using variational mode decomposition for fault diagnosis of gearbox under varying speed. Eng. Fail. Anal. 2020, 107, 104204. [Google Scholar] [CrossRef]
  19. Sánchez, R.; Lucero, P.; Macancela, J.; Cerrada, M.; Vásquez, R.E.; Pacheco, F. Multi-fault Diagnosis of Rotating Machinery by Using Feature Ranking Methods and SVM-based Classifiers. In Proceedings of the International Conference on Sensing, Diagnostics, Prognostics, and Control (SDPC), Shanghai, China, 16–18 August 2017; pp. 105–110. [Google Scholar] [CrossRef]
  20. Wang, Y.; Zhou, G. The Novel Successive Variational Mode Decomposition and Weighted Regularized Extreme Learning Machine for Fault Diagnosis of Automobile Gearbox. Shock. Vib. 2021, 2021, 5544031. [Google Scholar] [CrossRef]
  21. Kozłowski, E.; Mazurkiewicz, D.; Żabiński, T.; Prucnal, S.; Sęp, J. Machining sensor data management for operation-level predictive model. Expert Syst. Appl. 2020, 159, 113600. [Google Scholar] [CrossRef]
  22. Fu, Y.; Gao, Z.; Liu, Y.; Zhang, A.; Yin, X. Actuator and Sensor Fault Classification for Wind Turbine Systems Based on Fast Fourier Transform and Uncorrelated Multi-Linear Principal Component Analysis Techniques. Processes 2020, 8, 1066. [Google Scholar] [CrossRef]
  23. Zhao, R.; Yan, R.; Chen, Z.; Mao, K.; Wang, P.; Gao, R.X. Deep learning and its applications to machine health monitoring. Mech. Syst. Signal Process. 2019, 115, 213–237. [Google Scholar] [CrossRef]
  24. Qiu, G.; Gu, Y.; Cai, Q. A deep convolutional neural networks model for intelligent fault diagnosis of a gearbox under different operational conditions. Measurement 2019, 145, 94–107. [Google Scholar] [CrossRef]
  25. Zhang, S.; Zhang, S.; Wang, B.; Habetler, T.G. Deep Learning Algorithms for Bearing Fault Diagnostics—A Comprehensive Review. IEEE Access 2020, 8, 29857–29881. [Google Scholar] [CrossRef]
  26. Verstraete, D.; Ferrada, A.; Droguett, E.L.; Meruane, V.; Modarres, M. Deep Learning Enabled Fault Diagnosis Using Time-Frequency Image Analysis of Rolling Element Bearings. Shock. Vib. 2017, 2017, 5067651. [Google Scholar] [CrossRef]
  27. Barcelos, A.S.; Cardoso, A.J.M. Current-Based Bearing Fault Diagnosis Using Deep Learning Algorithms. Energies 2021, 14, 2509. [Google Scholar] [CrossRef]
  28. Hu, B.; Yang, J.; Li, J.; Li, S.; Bai, H. Intelligent Control Strategy for Transient Response of a Variable Geometry Turbocharger System Based on Deep Reinforcement Learning. Processes 2019, 7, 601. [Google Scholar] [CrossRef] [Green Version]
  29. Ren, L.; Sun, Y.; Cui, J.; Zhang, L. Bearing remaining useful life prediction based on deep autoencoder and deep neural networks. J. Manuf. Syst. 2018, 48, 71–77. [Google Scholar] [CrossRef]
  30. Zhang, S.; He, Y.; Zhang, J.M.; Zhao, Y. Multi-Fractal Based Fault Diagnosis Method of Rotating Machinery. Appl. Mech. Mater. 2012, 130–134, 571–574. [Google Scholar] [CrossRef]
  31. Liu, H.; Wang, X.; Lu, C. Rolling bearing fault diagnosis based on LCD–TEO and multifractal detrended fluctuation analysis. Mech. Syst. Signal Process. 2015, 60–61, 273–288. [Google Scholar] [CrossRef]
  32. Du, W.; Guo, Z.; Gong, X. Multifractal characterization of mechanical vibration signals through improved empirical mode decomposition-based detrended fluctuation analysis. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2017, 231, 4139–4149. [Google Scholar] [CrossRef]
  33. Wang, H.; Guo, Z.; Du, W. Diagnosis of rolling element bearing based on multifractal detrended fluctuation analyses and continuous hidden markov model. J. Mech. Sci. Technol. 2021, 35, 3313–3322. [Google Scholar] [CrossRef]
  34. Feng, Y.; Lu, B.; Zhang, D. Multifractal manifold for rotating machinery fault diagnosis based on detrended fluctuation analysis. J. Vibroeng. 2016, 18, 5153–5173. [Google Scholar] [CrossRef] [Green Version]
  35. Liu, H.; Zhang, J.; Cheng, Y.; Lu, C. Fault diagnosis of gearbox using empirical mode decomposition and multifractal detrended cross-correlation analysis. J. Sound Vib. 2016, 385, 350–371. [Google Scholar] [CrossRef]
  36. Moura, E.P.; Vieira, A.P.; Irmao, M.A.S.; Silva, A.A. Applications of detrended-fluctuation analysis to gearbox fault diagnosis. Mech. Syst. Signal Process. 2009, 23, 682–689. [Google Scholar] [CrossRef]
  37. Li, J.; Wei, H.; Fan, L. Multifractal Detrended Fluctuation Analysis of Frictional Vibration Signals in the Running-in Wear Process. Tribol Lett. 2017, 65, 50. [Google Scholar] [CrossRef]
  38. Du, W.; Kang, M.; Pecht, M. Fault Diagnosis Using Adaptive Multifractal Detrended Fluctuation Analysis. IEEE Trans. Ind. Electron. 2020, 67, 2272–2282. [Google Scholar] [CrossRef]
  39. Wendt, H.; Roux, S.G.; Jaffard, S.; Abry, P. Wavelet leaders and bootstrap for multifractal analysis of images. Signal Process. 2009, 89, 1100–1114. [Google Scholar] [CrossRef] [Green Version]
  40. Du, W.; Tao, J.; Li, Y.; Liu, C. Wavelet leaders multifractal features based fault diagnosis of rotating mechanism. Mech. Syst. Signal Process. 2014, 43, 57–75. [Google Scholar] [CrossRef]
  41. Kantelhardt, I.W. Fractal and Multifractal Time Series, Mathematics of Complexity and Dynamical Systems; Springer: New York, NY, USA, 2011; pp. 463–487. [Google Scholar]
  42. Loutridis, S.J. An algorithm for the characterization of time-series based on local regularity. Physica A 2007, 381, 383–398. [Google Scholar] [CrossRef]
  43. Jaffard, S. Multifractal formalism for functions part I: Results valid for all functions. Society for industrial and applied mathematics. SIAM J. Math. Anal. 1997, 28, 944–970. [Google Scholar] [CrossRef]
  44. Chen, A.J.; Li, Z.; Pan, J.; Chen, G.; Zi, Y.; Yuan, J.; Chen, B.; He, Z. Wavelet transform based on inner product in fault diagnosis of rotating machinery. Mech. Syst. Signal Process. 2016, 135, 70–71. [Google Scholar] [CrossRef]
  45. Lin, J.; Chen, Q. Fault diagnosis of rolling bearings based on multifractal detrended fluctuation analysis and Mahalanobis distance criterion. Mech. Syst. Signal Process. 2013, 38, 515–533. [Google Scholar] [CrossRef]
  46. Puchalski, A.; Komorska, I. A generalised entropy in multifractal time signals analysis of mechanical vibration. J. Vibroeng. 2018, 20, 1667–1675. [Google Scholar] [CrossRef]
  47. Puchalski, A.; Komorska, I. Stable Distributions and Fractal Diagnostic Models of Vibration Signals of Rotating Systems. In Advances in Condition Monitoring of Machinery in Non-Stationary Operations; CMMNO 2016; Timofiejczuk, A., Chaari, F., Zimroz, R., Bartelmus, W., Haddar, M., Eds.; Springer: Cham, Switzerland, 2018; Volume 9. [Google Scholar] [CrossRef]
  48. Puchalski, A.; Komorska, I. Data-driven monitoring of the gearbox using multifractal analysis and machine learning methods. In Proceedings of the MATEC Web of Conferences 3rd International Conference of Computational Methods in Engineering Science—CMES’18, Kazimierz Dolny, Poland, 22–24 November 2018; EDP Sciences: Paris, France, 2019; Volume 252, p. 06006. [Google Scholar] [CrossRef] [Green Version]
  49. Dąbrowski, Z.; Dziurdź, J.; Górnicka, D. Utilisation of the Coherence Analysis in Acoustic Diagnostics of Internal Combustion Engines. Arch. Acoust. 2017, 42, 475–481. [Google Scholar] [CrossRef] [Green Version]
  50. Zawisza, M. Diagnosing mechanical damages not detected by the OBD system of diesel engines. J. Vibroeng. 2017, 19, 1833–1844. [Google Scholar] [CrossRef]
Figure 1. The method scheme.
Figure 1. The method scheme.
Sensors 21 07677 g001
Figure 2. Multifractal analysis: (a) time series, (b) multifractal spectra with characteristic points, (c) multifractal spectra for Daubechies mother wavelets of orders 1 to 6.
Figure 2. Multifractal analysis: (a) time series, (b) multifractal spectra with characteristic points, (c) multifractal spectra for Daubechies mother wavelets of orders 1 to 6.
Sensors 21 07677 g002
Figure 3. Simulated signal and five intristic mode functions (a) and their multifractal spectra (b).
Figure 3. Simulated signal and five intristic mode functions (a) and their multifractal spectra (b).
Sensors 21 07677 g003aSensors 21 07677 g003b
Figure 4. Diagram of the EMD- WLMF method.
Figure 4. Diagram of the EMD- WLMF method.
Sensors 21 07677 g004
Figure 5. The test stand [48].
Figure 5. The test stand [48].
Sensors 21 07677 g005
Figure 6. Time series of the vibration signal in various operating states (a) fault-free state; (b) misalignment 1/3°; (c) increased backlash 0.2 mm; (d) worn teeth second stage; (e) worn teeth second stage, and increased backlash 0.2 mm.
Figure 6. Time series of the vibration signal in various operating states (a) fault-free state; (b) misalignment 1/3°; (c) increased backlash 0.2 mm; (d) worn teeth second stage; (e) worn teeth second stage, and increased backlash 0.2 mm.
Sensors 21 07677 g006
Figure 7. Multifractal spectra.
Figure 7. Multifractal spectra.
Sensors 21 07677 g007
Figure 8. Classification of damage (a) scatter plot of three parameters of the multifractal spectrum, (b) confusion matrix for the test set performed with SVM method.
Figure 8. Classification of damage (a) scatter plot of three parameters of the multifractal spectrum, (b) confusion matrix for the test set performed with SVM method.
Sensors 21 07677 g008
Figure 9. Test object (a) view of the gearbox, (b) pair of gear wheels with signs of wear.
Figure 9. Test object (a) view of the gearbox, (b) pair of gear wheels with signs of wear.
Sensors 21 07677 g009
Figure 10. Time series of the first IMF for (a) initial state, (b) initial wear, (c) advanced wear, (d) after repair, (e) examples of rotational speed changes during registration.
Figure 10. Time series of the first IMF for (a) initial state, (b) initial wear, (c) advanced wear, (d) after repair, (e) examples of rotational speed changes during registration.
Sensors 21 07677 g010
Figure 11. Sample multifractal spectra of (a) EMD-WLMF, (b) WLMF.
Figure 11. Sample multifractal spectra of (a) EMD-WLMF, (b) WLMF.
Sensors 21 07677 g011
Figure 12. Damage classification based on three features of the multifractal spectrum (a) and confusion matrix (b).
Figure 12. Damage classification based on three features of the multifractal spectrum (a) and confusion matrix (b).
Sensors 21 07677 g012
Figure 13. View of the damage of outlet valve.
Figure 13. View of the damage of outlet valve.
Sensors 21 07677 g013
Figure 14. Acceleration time series of the engine head vibrations—the first IMF for (a) base state, (b) defect 1, (c) defect 2.
Figure 14. Acceleration time series of the engine head vibrations—the first IMF for (a) base state, (b) defect 1, (c) defect 2.
Sensors 21 07677 g014
Figure 15. Multifractal spectrum of the vibration signal for three states of valves.
Figure 15. Multifractal spectrum of the vibration signal for three states of valves.
Sensors 21 07677 g015
Figure 16. Damage classification based on three features of the multifractal spectrum (a) and confusion matrix (b).
Figure 16. Damage classification based on three features of the multifractal spectrum (a) and confusion matrix (b).
Sensors 21 07677 g016
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Komorska, I.; Puchalski, A. Rotating Machinery Diagnosing in Non-Stationary Conditions with Empirical Mode Decomposition-Based Wavelet Leaders Multifractal Spectra. Sensors 2021, 21, 7677. https://doi.org/10.3390/s21227677

AMA Style

Komorska I, Puchalski A. Rotating Machinery Diagnosing in Non-Stationary Conditions with Empirical Mode Decomposition-Based Wavelet Leaders Multifractal Spectra. Sensors. 2021; 21(22):7677. https://doi.org/10.3390/s21227677

Chicago/Turabian Style

Komorska, Iwona, and Andrzej Puchalski. 2021. "Rotating Machinery Diagnosing in Non-Stationary Conditions with Empirical Mode Decomposition-Based Wavelet Leaders Multifractal Spectra" Sensors 21, no. 22: 7677. https://doi.org/10.3390/s21227677

APA Style

Komorska, I., & Puchalski, A. (2021). Rotating Machinery Diagnosing in Non-Stationary Conditions with Empirical Mode Decomposition-Based Wavelet Leaders Multifractal Spectra. Sensors, 21(22), 7677. https://doi.org/10.3390/s21227677

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop