Design of Combined Scalar Indicators For Fault Detection of Rotating Machines
Design of Combined Scalar Indicators For Fault Detection of Rotating Machines
12(10), 974-987
RESEARCH ARTICLE
DESIGN OF COMBINED SCALAR INDICATORS FOR FAULT DETECTION OF ROTATING
MACHINES
the bearing is greased or when the gears are loaded, and that Kurtosis increases according to the size of the defect. In
parallel with their work, in this paper, we show that Kurtosis is not able to detect the severity of a bearing defect
compared to another defect on the bearing, but as in their case, it is sensitive to shocks, especially to gear defects.
Other criteria are also evaluated in this work. The authors, Antoni and Sidahmed briefly presented different methods
for fault diagnosis in acoustic and vibration signals, applicable to stationary or non-stationary signals. In the case of
non-stationary signals, especially cyclostationary signals where the speed is variable, a direct analysis of the signal
would not give a correct result, because the number of points is not the same in each period of the signal. Specific
signal processing techniques are required to complement or improve the presented methods as presented by Abboud,
2015. To analyze signals in this context, one approach is to resample the signal (angular resampling), apply
synchronous averaging to eliminate first-order components of the signal, apply filtering (Antoni and Randall, 2006),
and finally apply spectral correlation to detect second-order faults. Lejeune et al. 1997 also present work on first-
and second-order cyclostationarities applied to gear vibration signals. In the following sections, we present our
approach for the work.
For this work, the signals delivered by the sensor of channel 3 are used. Table 1 below shows the tests carried out
and the typology of the defects.
975
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
Methods:-
Method 1: Fault detection by traditional scalar indicators
In this work, we use the temporal analysis of vibration signals, which is a method of analyzing stationary signals in
the time domain. The temporal analysis makes it possible to highlight the presence of faults in machines. In this
analysis, the presence of faults and, if necessary, their severity in rotating machines are detected. The temporal
methods are based on the statistical analysis of the collected signal and use well-known scalar indicators to follow
the evolution of a quantity derived from the power or peak amplitude of the signal. These indicators include, among
others, the RMS (effective value), the Kurtosis, the crest factor, the peak-to-peak and the K factor presented below.
The crest factor has the advantage of detecting defects before the effective value. This comes from the fact that for a
bearing without defects, the ratio remains approximately constant and increases when degradation appears, whereas
the crest value increases while the effective value remains approximately constant.
976
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
The Kurtosis:-
Kurtosis is an indicator that characterizes shocks in mechanical organs. It is a statistical quantity that allows us to
analyze the “pointy” or “flat” character of a distribution and therefore to observe the shape of the signal. Kurtosis is
defined as the ratio of the mean value of the signal raised to the power of 4 to the square of its energy (see equation
3):
N (X t −x )4
i=1
N
Kurt = N (s i −x )2 (3)
( i=1 )2
N
A bearing in good condition generates a vibration signal with a Kurtosis close to 3. For a degraded bearing,
presenting flaking, indentations or significant mechanical error or irrespective of tolerances, the shape of the
distribution of the signal amplitude is modified and the Kurtosis is greater than 3.
The K Factor:-
K factor of a signal is defined as the product of the peak value by the effective value. It is an indicator specific to
bearings. It is expressed by equation (4):
N X(i)2
i=1
K = Acc crete ∗ Acc eff = Sup|X i | × (4)
N
The interpretation of the crest factor is done through its evolution as the bearing deteriorates. The value of K
increases with bearing wear.
I1 = a ∗ RMS + b ∗ FC + c (7)
With a, b et c are constants to be determined.
For all indicators, we can write the system of equations in the following matrix form (equation 8):
I = H. X (8)
Where I : Column matrix of combined indicators, with:
𝐈𝟏
𝐈𝟐
𝐈
𝐈 = 𝟑
⋮
⋮
𝐈𝐍
H : Matrix of scalar indicators, with:
977
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
𝐑𝐌𝐒 𝐅𝐂 𝟏
𝐑𝐌𝐒 𝐊𝐔𝐑𝐓 𝟏
𝐊 𝐅𝐂 𝟏
𝐇 =
⋮ ⋮ 𝟏
⋮ ⋮ 𝟏
𝐏𝐏 𝐊𝐔𝐑𝐓 𝟏
X : Column matrix of constants a, b and c, with:
𝐚
𝐗 = 𝐛
𝐜
The coefficients a, b and c are determined using equation (9). Since the matrix [H] is not square, there is no exact
solution to this system of equations. The pseudo-inverse of Gauss (see Boukar et al., 2014 and Bernard G-B., 1967)
or the least square minimization of multilinear regression is given by equation (10), it allows to find a better solution
for a, b and c.
X = H∗ . I (9)
With: H ∗ = ( H t . H )−1 . H t the Pseudo-inverse of Gauss.
This gives equation (10):
A
B =
c
−1 I1
RMS FC 1
RMS KURT 1 I2
RMS RMS K … … PP RMS RMS K … … PP
K FC 1 I3
FC KURT FC… … KURT X X FC KURT FC… … KURT X
⋮
(10)
⋮ ⋮ 1
1 1 1… … 1 1 1 1… … 1
⋮ ⋮ 1 ⋮
PP KURT 1 IN
The coefficients of the matrix [H] are the scalar indicators such as RMS, Kurtosis, FC… determined using the
vibration measurements carried out on the experimental test bench in Figure 1. The measurement results were
analyzed using MATLAB®. For detecting defects with the combined scalar indicators I, initial values were imposed
on these indicators depending on the severity of the defects. Equations (11) and (12) make it possible to decide
whether or not a defect is present:
If Icalculed ≤ Iinitial , no fault (11)
If Icalculed > Iinitial , there is a fault (12)
In our case, we calculate several R² between the scalar indicators (2-parameter correlations). The best R² is obtained
between the peak value and the crest factor as presented in Figure 2. The values of the six scalar indicators studied
(Kurtosis, RMS, FC, P, PP, and K) are close to the regression line in Figure 2(a) presented below.
978
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
Fig. 2c-R² between Kurtosis and K Fig. 2d-R² between Kurtosis and P
Figure 1 above shows the results of some values of R² allowing to obtain the values of a, b and c. Among these
results, we retain the best R² = 0, 99719 obtained between the peak-to-peak value of the acceleration and the crest
factor. Which gives the constants a = - 0, 0099356 and b = 0, 44901 et c = - 2.0601.
To obtain these results of R², a, b and c, the simulations were carried out with the signals from the tests of six (6)
fault configurations among the twelve (12) presented in Table 1. The six signals chosen to calculate these
parameters are those from the following tests:
Test 1: Healthy bearing (Rs);
979
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
Indeed, there are fifteen (15) combined scalar indicators obtained by linear combination of the six traditional scalar
indicators. These constants will therefore be used to calculate the remaining fourteen (14) combined scalar
indicators in order to determine the best combined scalar indicator among the fifteen. The results are presented in the
following section.
K = Acc(peak)*Acc(eff)
2.5
2
Rs
1.5 Rd_bext
Rd_bint
1 Rd_cage
0.5
980
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
14
12 Peak factor
10
8
6
4
2
0
981
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
RMS
12
10
8
6
4
2
0
Kurtosis
12
10
8
6
4
2
0
982
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
Peak value
120
100
80
60
40
20
0
983
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
As mentioned above, these constants are used to calculate the other fourteen (14) combined scalar indicators. To
calculate these values, we arbitrarily imposed in the simulation an initial value for each combined scalar indicator
(see equation 10). This initial value is set according to the severity of the defect for the six (6) tests used to calculate
R², a, b and c.
Thus, for:
Test 1 (Rs) : I_initial = 0,1 ;
Test 2 (Rs_Es) : I_initial = 0,45 ;
Test 3 (Rs_Ed) : I_initial = 0,6 ;
Test 7 (Rd_bint) : I_initial = 0,5 ;
Test 8 (Rd_bint_Es) : I_initial = 0,8 ;
Test 9 (Rd_bint_Ed) : I_initial = 1,4.
To know whether there is a defect or not, the decision is given by equations (11) and (12) above.
To better evaluate the method, we proceeded to the tests in the following manner: a test is carried out with a healthy
bearing or a defective bearing, then a healthy mesh is added, then a defective gear so as to have a single defect or a
multiple defect. Of the fifteen (15) combinations of the combined scalar indicators, fourteen (14) do not give the
expected results. In the 14 calculated indicators, at least one test out of twelve is not verified. Only one combined
indicator gives satisfactory results in the twelve configurations of defects tested. This is the following combined
scalar indicator:
The results obtained by this indicator (see Fig. 4c) show that when a healthy gear is added, it is detected, because the
value of the indicator I increases compared to the case without gear. When a defective gear is added, its value
increases further. This shows the sensitivity of this indicator to shocks, and therefore its sensitivity to gear defects.
In the context of the detection of bearing defects, the result is also satisfactory. Indeed, the initial value of the
indicator chosen for test 1 (Healthy bearing) is almost the same as the calculated one, this means that it is able to
recognize when there is no bearing defect (Iinitial = 0,1 et Icalculed = 0,1246).
From all of the above, we logically retain for the temporal analysis, the monitoring indicator I obtained by linear
combination of RMS and FC.
Figure 4c below shows the results of this indicator. We can see that the values of the indicator calculated for the
tests with defects are higher than their initial values. When there is no defect, its two values are close. Also, as we
can see in this figure, the defects of low amplitude are less highlighted, for example tests 4 (Rd_bext) and 7
(Rd_bint).
Rd_cage_Ed Rd_cage_Ed
Rd_cage_Es Rd_cage_Es
Rd_cage Rd_cage
Rd_bint_Ed Rd_bint_Ed
Rd_bint_Es Rd_bint_Es
Rd_bint Rd_bint
Rd_bext_Ed Rd_bext_Ed
Rd_bext_Es Rd_bext_Es
Rd_bext Rd_bext
Rs_Ed Rs_Ed
Rs_Es Rs_Es
Rs Rs
0 10 20 30 40 50
0 1 2 3 4
Ic1=a*CC+b*FC+c Ic initinial Ic2=a*FC+b*C+c Ic initinial
Fig. 4a-I1 = f (CC, FC) Fig. 4b -I2 = f (FC, C)
984
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
Rd_cage_Ed
Rd_cage_Ed
Rd_cage_Es
Rd_cage_Es
Rd_cage
Rd_cage
Rd_bint_Ed
Rd_bint_Ed
Rd_bint_Es
Rd_bint_Es
Rd_bint
Rd_bint
Rd_bext_Ed Rd_bext_Ed
Rd_bext_Es Rd_bext_Es
Rd_bext Rd_bext
Rs_Ed Rs_Ed
Rs_Es Rs_Es
Rs Rs
0 2 4 6 8 10
0 1 2 3 4
Ic3=a*RMS+b*FC+c Ic initinial Ic4=a*K+b*FC+c Ic initinial
Fig. 4c-I3 = f (RMS, FC) Fig. 4d-I4 = f (K, FC)
Rd_cage_Ed Rd_cage_Ed
Rd_cage_Es Rd_cage_Es
Rd_cage Rd_cage
Rd_bint_Ed Rd_bint_Ed
Rd_bint_Es Rd_bint_Es
Rd_bint Rd_bint
Rd_bext_Ed Rd_bext_Ed
Rd_bext_Es Rd_bext_Es
Rd_bext Rd_bext
Rs_Ed Rs_Ed
Rs_Es Rs_Es
Rs Rs
985
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
Rd_cage_Ed Rd_cage_Ed
Rd_cage_Es Rd_cage_Es
Rd_cage Rd_cage
Rd_bint_Ed Rd_bint_Ed
Rd_bint_Es Rd_bint_Es
Rd_bint Rd_bint
Rd_bext_Ed Rd_bext_Ed
Rd_bext_Es Rd_bext_Es
Rd_bext Rd_bext
Rs_Ed Rs_Ed
Rs_Es Rs_Es
Rs Rs
Conclusion:-
In this article, we proposed a new indicator for monitoring rotating machines by vibration analysis, based on the
linear combination of scalar indicators commonly used in maintenance by time analysis. Fifteen different indicators
are studied and evaluated on twelve tests, with different fault configurations performed on a vibration test bench.
Only one combined indicator out of the fifteen studied gives satisfactory results on all the tests performed, namely
the indicator obtained by linear combination of RMS and FC. This is what we now retain for the detection of faults
in the time domain. However, to definitively validate this indicator, full-scale tests would have to be carried out,
because we identified a misalignment of shafts on the experimental setup using spectral analysis, which is not
presented here in this work.
References:-
1. Zezz R. Chaib, S. Meziani, I. Verzea, Surveillance des roulements par analyse vibratoire, Sciences
&Technologie B – N°21, Juin (2004), pp. 23-27.
986
ISSN: 2320-5407 Int. J. Adv. Res. 12(10), 974-987
2. Wärsilä, Lubricating oils for WÄRTSILÄ® 32, and 32GD engines, DATA & SPECIFICATIONS, Technical
Services, Issue 7, 02 March 2016
3. R. Younes, N. Ouelaa, N. Hamzaoui, Optimisation d’indicateurs de défauts combinés d’engrenages et de
roulements par la perception sonore, 22ème Congrès Français de Mécanique Lyon, 24 au 28 Août 2015
4. N. Hamzaoui, R. Younes, N. Ouelaa, Indicators of real gear transmission defects using sound perception,
NOVEM 2015 Noise and vibration, Dubrovnik- Croatia, April 13-15 2015
5. A. Djebala, N. Ouelaa et N. Hamzaoui, « Optimisation de l’analyse multirésolution en ondelettes des signaux de
choc. Application aux signaux engendrés par des roulements défectueux », Mécanique & Industries 8, 379-389,
© AFM, EDP Sciences 2007, DOI : 10.1051/meca:2007060
6. A. Djebala, N. Ouelaa, C. Benchaabane and D. F. Laefer, Application of the wavelet Multirésolution Analysis
and Hilbert transform for the prediction of gear tooth defects, Meccanica, DOI: 10.1007/s11012-012-9538-1,
2012
7. A. Djebala, N. Ouelaa, N. Hamzaoui and L. Chaabi, Detecting mechanical failures inducing periodical shocks
by wavelet multirésolution analysis. Application to rolling bearings faults diagnosis, ISSN 1392-1207.
MECHANIKA, 2006. Nr. 2(58)
8. A. Djebala, Habilitation universitaire en Génie Mécanique, Université 08 Mai 1945 Guelma, Mai 2012
9. M. El Badaoui, F. Guillet, N. Nejjar, P. Martini et J. Danière, « Diagnostic d'un train d'engrenages par analyse
cepstrale synchrone », Seizième Colloque GRETSI, 15-19 Septembre 1997, Grenoble
10. J. Antoni, M. Sidahmed, Contrôle et diagnostic à partir des signaux acoustiques et vibratoires, Acoustique &
Techniques n° 38, Spécial congrès Surveillance 5
11. D. Abboud, Vibration-based condition monitoring of rotating machines in non-stationary regime, Thesis INSA
Lyon, 2015
12. J. Antoni, R. B. Randall, The spectral kurtosis: application to the vibratory surveillance and diagnostics of
rotating machines, Mechanical Systems and Signal Processing vol. 20, pp. 308–331, 2006
13. G. Lejeune, J. L. Lacoume, P. Marchand, M. Durnerin, N. Martin, J. Liénard, A. Silvent CEPHAG,
Cyclostationnarités d’ordre 1 et 2 : application à des signaux vibratoires d’engrenages, Seizième Colloque
GRETSI, 15-19 Septembre 1997, Grenoble
14. A. Boukar, M. Pillet, E. Pairel, and A. Sergent, “Correction of machining operations with the Total Inertial
Steering” , Int. J. Metrol. Qual. Eng. 5, 21–27 (2014), 103 (2014), © EDP Sciences 2014, DOI:
10.1051/ijmqe/2014003
15. Bernard Germain-Bonne, « Méthodes de calcul de pseudo-inverses de matrices », Thèse de l’université de Lille,
1967.
987