Vibration Diagnostics Methods
Vibration Diagnostics Methods
ABSTRACT
This paper discusses the marine diesel engine vibration diagnostics methods, according to the international
standards. The analysis of vibration diagnostics methods applied to piston engines is being shown. The
analysis of root-mean-square value of vibration velocity (RMS) for various classes of mechanisms is being
shown. The main defects of various diesel units and the frequencies of the harmonics that indicates their
occurrence have been presented. A method for eliminating the “leakage effect” in the discrete spectrum
DFT is offered. The proposed method is based on solving a system of complex equations. The article
considers the method for diagnosing a turbocharger, based on the spectral analysis of vibroacoustic signals
of a compressor.
1.0 INTRODUCTION
The structure of marine diesel engines combines various mechanisms and units associated with different
functions and operations, such as:
- reciprocating motion mechanisms (the crank mechanism, pistons);
- rotary mechanisms (gear, belt and chain drives, oil and water pumps, crankshafts and camshafts);
- turbocharger;
- high pressure fuel injection elements (high pressure fuel pumps, valves and injectors);
- gas distribution mechanism (inlet and exhaust valve drives);
- roller bearings and connecting couplings;
- other components (for instance, generator as a part of a diesel-generator set).
Normal functioning of a marine diesel engine is ensured by a complex of different operational mechanisms.
Each one of these mechanisms generates its own specific vibrations, which greatly complicate the task of
vibration diagnosis of marine diesel malfunctions, despite relatively easy malfunction diagnosis based on
analyzing vibrations of rotary type machinery. The fact that the vibration signal at some location results
from a sum of very heterogeneous signals from the different nearby mechanisms complicates the task even
more. In this case, it is logical to simplify the task of diagnosis by analysing the individual components and
mechanisms [1].
STO-MP-AVT-306 9-1
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
fn = 1500 / 60 = 25 Hz
fcyl = 25 × 6 × 0,5 = 75 Hz
There is an evaluation method of engine technical condition based on analyzing vibration frequencies of
non-rotating parts (ISO 10816[2]). According to the requirements of [2], technical condition diagnosis is
dependent on general vibration class for each mechanism. In order to analyze this, certain mechanisms RMS
vibrospeed values between 10-1000 Hz are being used. For a continuous signal V(t), RMS is defined as
V
T
~RMS = 1
T 0
2
(t )dt , where T is the sampling period, which should be longer (at least 10 times) than any
basic period of any analyzed signal fundamental frequencies contained in V(t).
In case of a discretely recorded vibration signal Vi (i = 1,2,…N) and availability of N values of vibrospeed
Vi, the RMS value is being defined as
N
VRMS = 1
N i =1
Vi 2 .
Certain filters have to be applied for choosing between necessary frequency range of 10 - 1000 Hz when
calculating RMS for a time signal. VRMS calculations are done easier using spectrum amplitude values sj
k2
VRMS = 1
2 j = k1
s 2j ,
where k1 and k2 are harmonics indexes between 10 Hz and 1000 Hz, respectively.
It is possible to estimate overall technical condition of the mechanism by calculating the RMS value using
the Table 1, where vibrospeed RMS limits for normal (A, B) and abnormal (C, D) conditions for
mechanisms of distinct classes are specified. MAN D20 diesel corresponds to a class III in case of its
installation on the rigid basis and to a class IV in case of an installation on flexible support.
Table 2-1: ISO 10816: zones* vibrospeed RMS for mechanisms of various classes [2].
Vibration level higher than the normal indicates existing problems in the mechanism, like cylinder
9-2 STO-MP-AVT-306
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
unbalance, breaking or misalignment of couplings, fasteners weakening, mounting problems, etc. Since each
unit of a diesel engine is unique in its vibration, the most effective method of diagnosis is an over-time
monitoring of measurement results. Such monitoring can be accomplished by trending (scatter plan of RMS
values of vibration velocity in time), see Fig. 2-1.
Table 2-1 shows that the vibration level is determined by the vibration RMS and describes the general state
of the mechanism. Specific future defects in individual units can be predicted when analyzing the vibration
spectrum, as shown in Table 3-1.
STO-MP-AVT-306 9-3
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
Table 3 - 1. Continued
Parallel shear or fracture of coupling cause
Vibration couplings fn mass unbalance of rotating rotor due to the
in the axial direction axis displacement of the mass centers.
Axial vibration with frequency fn appears
Vibration depends on a condition of working
surfaces of sliding bearings, gaps between
shaft and bearings as well as on properties of
Increased vibration of FCR =(0,42÷0,48) × fn lubricating oil. Under certain conditions
crankshaft bearings in the radial direction oscillations of the shaft in the radial direction
may arise. Appearance in the spectrum of
harmonics FCR may characterize defect inserts
or shaft
- defect in an outer ring BPFO: Defects or damage to the rolling contact
z d surfaces (outer / inner ring, separator or
f0 = f n 1 − cos rolling elements).
2 D More often vibration acceleration is analyzed
- defect of the inner ring BPFI: (or vibration velocity in some cases).
Vibration of rolling bearings z d z - number of rolling elements;
(generators etc.) fi = f n 1 + cos D - diameter of a circle passing through the
2 D center of the rolling elements (average
- separator defect FTF: diameter of the separator);
fn d d - the diameter of the rolling element;
fс = 1 − cos β - contact angle of the rolling element.
2 D
- defect of the rolling elements
BSF:
D d
2
fr = f n 1 − cos
d D
- BPFO: .
f 0 0,45 z RPM
Vibration of rolling bearings z - number of rolling elements;
(simplified formula, through - BPFI: RPM - revolutions per minute of crankshaft
RPM) fi 0,55 0,6 z RPM
- BSF:
f r 3,5 RPM
Deviation from axial symmetry of the rotor
and the heterogeneity of the material.
Turbocharger rotor unbalance fTUR Causes forced oscillation with a speed of the
turbocharger due to rotation of the
unbalanced centrifugal forces.
fTUR × nb.comp
Vibration of rotor blades of the Compressor blade vibration is always present
air compressor (gas turbine) in the vibration spectrum and does not
nb.comp - number of the
turbocharger indicate a defect *.
compressor wheel blades of
Defect diagnosed by increased levels of blade
turbocharger
harmonics over time.
Defect formed during the assembling or
Defect turbocharger mounting fn during operation.
in the vertical and transverse The general level of vibration and harmonic
directions spectrum rises at engine speed in the vertical
and transverse directions.
* Presence of a bladed compressor harmonica in a spectrum of vibration allows carrying out the
turbocharger express analysis [7].
9-4 STO-MP-AVT-306
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
where Am , Ai with i 1...9 are the amplitudes of the m-th and i-th harmonics in the spectrum.
The references [4, 6, 7] – from which Table 3-2 was extracted – propose a number of criteria for the
condition diagnosis of reciprocating engines. We think that those criteria and table, directly or in a modified
form, can also be used for marine diesel engines.
Rotating details of a
shaft, the coupling Imbalance K ((11)9 )
Bearing of
reciprocating
Shaft of the engine and
engine from K ((1293))
generator Misalignment
the drive side
or flywheel
Coupling Increased clearance, slackening, decreasing of K ((1399))
rigidity
Two facts emerge from the analysis of the fore mentioned proposed reciprocating engines diagnostic
method:
• The use of the squares of multiple harmonics quotient normalizes levels of diagnostic coefficients in
the range [0 to 1];
• Identical coefficients calculated from sensor data installed at various locations contain information
about the corresponding units and mechanisms condition.
STO-MP-AVT-306 9-5
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
Unfortunately, this technique [6, 7] could not be extended yet to the most important units of modern diesels -
the high pressure fuel equipment (nozzles and high pressure fuel pumps). It is known that the more frequent
operational failures of ship diesels - from 50% to 80% - correspond to the high pressure fuel equipment [8,
9].
The proposed methodology in [4, 6, 7] has to be checked for high-speed diesels. The equations for
calculation of diagnostic coefficients have to be specified and modified if necessary. Thus it is possible to
recognize that in a spectrum of diesel vibration the most informative harmonics are 0,5fn, fn, fcyl and multiples
of them up to 20-th harmonic.
As an example, Fig. 4-1 shows the amplitude spectra of the same sinusoidal signal with an integer (a) and a
non-integer (b) number of samples per one signal period. Let the frequency of a signal be represented by γ =
M/T, where T is the period of the signal; M = n + σ, where n is an integer and 0< σ <1, then the maximum
distortions of the amplitude, frequency and phase of the central harmonica and leakage of power into the
neighbouring ones will be observed at σ = 0,5 See [12].
a) σ = 0
b) σ ≠ 0
Figure 4 - 1: DFT leakage effect
Thus, when analyzing the parameters of the original spectrum signal, i.e. the central harmonic, the resulting
amplitude, frequency and phase will be distorted in the case of a non-integer number of signal samples per its
period. In practice, for discrete recording of signals, an ADC with a selected and fixed sampling rate is used.
It is understandable that the number of samples per period will never be an integer and the value of σ will
change from 0 to 1 depending on the natural frequency of the measured signal, and the accuracy of
estimating the signal parameters along the central harmonic will change.
9-6 STO-MP-AVT-306
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
The most common solution for leakage effect reduction is based on window transform methods. The essence
of the method is simple: to reduce the number of discontinuities at the edges in order to reduce leakage it is
necessary to reduce the amplitude of the signal near the edges. This scaling is carried out during the
implementation of the multiplication by the window with the special form
s wj = s j W ( j ) ,
As a result of applying window functions, the spectrum of the original signal is changed and its RMS
decreases by RMS Coeff times, as shown in Table 2. The decrease in RMS when using the Hanning window
is 0,707 / 0,433 = 1,633. Thus, the dependence of the amplitude, frequency and phase of the fundamental
harmonic in the spectrum from the value of σ decreases. This means that the fundamental harmonic swj can
be used to approximate the signal parameters with a certain constant error, which can be taken into account.
Table 4-1. Window functions and RMS coefficient used to reduce the effect of "leakage" [12]
2πn
Hemming window RMS Coeff = 1,414 (n) = 0,53836− 0,46164 cos
N −1
2πn
Hanning window RMS Coeff = 1,633 (n) = 0,5 1 − cos
N −1
I 0 1 − ( 2 nN−−N1+1 )
2
Kaiser window RMS Coeff = 1,61 ( n) =
I 0 ( )
Blackman-Harris RMS Coeff = 1,585
(n) = 0,42 − 0,5 cos(2πn /(N − 1)) + 0,8 cos(4πn /(N − 1))
window
More precisely, we can eliminate the “leakage” effect by a numerical method based on the processing of the
complex DFT results. In [12], a suggestion was made that the frequency m, the phase φ, and the amplitude A
of the original signal from the values of two maximum harmonics in the spectrum k k +1 should be
X ,X
specified. For this it is proposed to solve numerically the system of complex equations. To do so, the system
of complex equations is proposed to be solved numerically:
E (m, )k / E (m, ) k +1 = X k / X k +1
(1)
Arg( E (m, ) k ) = Arg( X k )
STO-MP-AVT-306 9-7
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
The harmonic coefficients can be represented in the form X k = (Ak / 2)E(m, )k , where E (m, )k is a
complex function independent of the amplitude, but dependent on the frequency and phase:
e 2 j ( m − k ) − 1 e −2 j ( m + k ) − 1
E (m, ) k = e j 2 j ( m − k )
+ e − j −2 j ( m + k )
e N
−1 e N
−1
The system of equations (1) must be solved in the case where the harmonics to the left and right of the
central one are not equal to zero (in practice it is more than a given small value δ):
X k −1 , X k +1 .
If X k −1 = 0, X k +1 = 0 , then the leakage effect is absent and the frequency, amplitude and phase of the central
harmonica correspond to parameters of the measured initial signal See Fig. 4-2d.
When solving the system (1) for the situation of strong leakage effects (σ ~ 0,5 See Fig. 4-2c), only five full
iterations were required to provide a specified error of less than 0.5% in frequency and phase. For a
sinusoidal signal, the amplitude and frequency are recovered to the value specified in the original signal with
accuracy to 5 decimal places. In this case, the amplitude of the central harmonic in the spectrum after the
DFT before the recovery procedure was with an error of 35% (!) See Fig. 4-2c [10, 11]
An error in estimating the frequency of the original signal with respect to the frequency of the central
harmonic can also be significant. It depends on the frequency of the ADC and the frequency of the original
signal. As the frequency of the ADC increases, the error in estimating the frequency will decrease.
The solution of the system (1) is not associated with additional memory as is the case for the fast Fourier
transform (FFT). Despite the iterative numerical solution for system (1), such procedure only very slightly
increases the overall computation time, and make it possible to obtain not only the spectrum of the signal,
but also the restored value of the fundamental frequency, amplitude and phase of the measured signal, when
it is close to sinusoidal.
The solution of system (1) is not connected with allocating additional memory for storing volumetric data
sets and computed coefficients, as in case of fast Fourier transform (FFT). In this regard, the algorithm can
be programmed on a modern DSP controller which implements the FFT.
9-8 STO-MP-AVT-306
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
Despite the iterative numerical solution (1), such a recovery procedure increases the total calculation time
very slightly and it allows recovery of not only the signal spectrum, but also the restored value of the
fundamental frequency and amplitude and phase of the measured signal, if it is close to sinusoidal.
This method was investigated in case of noise in the original signal (with a white noise of 5% and 10% of the
amplitude of the sinusoid). Fig. 4-2 shows the solution of the system (1) for a sinusoid with an amplitude of
0,8 and for the cases a) ϭ = 0,1, b) ϭ = 0,3 c) ϭ = 0,5 and c) ϭ = 0. The central green line in each figure 4-2 a,
b, c, d is the main harmonic of a sinusoid with amplitude of 0,8 with the restored amplitude, frequency and
phase, being a result of solving the system of equations (1).
For all the cases, not more than 5 complete iterations were required to ensure a given accuracy. As a result of
the solution of the system (1), the phase and frequency of the signal with the addition of white noise to 10%,
are restored to the initial value with an error of not more than 0.5%.
High efficiency of MAN ME and MAN MC diesel engines (with a real specific effective flow rate of 160-
170 g/kWh) is provided by the high charge air pressure, in particular. When the efficiency of the
turbocharger decreases, the efficiency of the diesel engine drops abruptly, the level of emission of carbon
oxides and soot increases.
During operation of marine diesel engines, the exhaust manifolds become clogged with products of
incomplete combustion. As a result, the throughput of the exhaust manifolds and the nature of the internal
flow of gases before the blades of the turbocharger impeller may vary. In this case, the appearance of
pulsations is possible which leads to rotor oscillation. The increased level of rotor oscillation creates
additional loads on the turbocharger bearings and reduces their life. In the event of microdefects in the
bearings of the turbocharger, the vibration level of the rotor increases even further that may lead to a severe
accident.
Constant operational monitoring of the vibration level of the turbocharger rotor can prevent an emergency
situation.
The experiments on diesel engines in laboratory and in sea conditions have revealed that the turbocharger
compressor blades generate oscillations which are always present in the overall vibration spectrum,
regardless of the technical condition of the turbocharger. The spectral analysis of the turbocharger vibration
has shown that the compressor blades generate a vibroacoustic signal with a frequency equal to the speed of
the turbocharger rotor multiplied by the number of air blades [11, 12],
υb = nb × RPMtur / 60,
where υb - blade frequency of the turbocharger compressor, Hz; nb - the number of compressor air blades,
RPMtur - the speed of the turbocharger rotor min-1.
To determine the blade frequency of the turbocharger compressor and the subsequent calculation of the
turbocharger speed, the amplitude spectrum of vibroacoustic signals was used. The recording was made
opposite the compressor air filter (see Figure 5-1) through the use of a broadband industrial microphone with
STO-MP-AVT-306 9-9
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
In the case of recording vibration of the turbocharger of diesel engine 6L80MCE with the
TURBOCHARGER VTR 564-31 (20 compressor blades), the expected blade frequency of the compressor
was calculated on the basis of the turbocharger rotor speed rating at the nominal conditions:
Due to the fact that the operating mode of the diesel engine was at a lower load, the expected speed of the
turbocharger rotor shall be less than the nominal one. Thus, the value of the blade frequency calculated for
the nominal conditions can be used as the upper limit for determining the actual operational value.
Figure 5-2 shows the vibration spectrum of the TURBOCHARGER VTR 564-31 recorded at a load close to
the nominal one. It can be seen from Figure 5-2 that the harmonic closest to 3 kHz has a frequency of 2948
Hz. The nearest harmonic on the left has a frequency of 1474 Hz and is a subharmonic with a frequency
equal to half of the blade frequency υb / 2. This leftmost subharmonic in the spectrum can be considered as
the left boundary when determining of the harmonic corresponding to the blade frequency of the
turbocharger compressor.
9-10 STO-MP-AVT-306
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
RPMtur = 60 × υb / nb,
RPMtur = 60 × 2948 Hz / 20 = 8844 RPM,
The regular tachometer of the turbocharger showed a rotation speed of 8800 RPM, which in comparison with
the value determined by the spectrum gives a relative error of 0,5%.
It is necessary to take into account the industrial accuracy class of the standard tachometer (division scale of
200 RPM).
Spectral analysis of a vibroacoustic signal recorded at a frequency of 44,1 kHz makes it possible to analyze
harmonics in steps of less than 1 Hz at a recorded signal frequency up to 20 kHz [12]. The blade frequency
of the turbocharger compressor is significantly lower. Thus, as a result of the spectral analysis of the
vibroacoustic signal of the turbocharger compressor, an error in determining the frequency less than 1 RPM
can be reasonably obtained. Such accuracy is much higher than the accuracy of the standard tachometers,
which makes it possible to use the blade frequency of the turbocharger compressor in accurate calculations
of the main rotational speed of the turbocharger and the subsequent estimation of the diesel engine power.
After determining the compressor blade frequency and the main speed of the turbocharger (RPMtur), we can
STO-MP-AVT-306 9-11
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
υ turbocharger = υb / nb
Obviously, if there is a significant increase in the amplitude Δ of the harmonic at the main speed of the
turbocharger rotor, this demonstrates an increased vibration of the rotor. Fig. 5 - 2 shows a slight increase in
the amplitude of the fundamental harmonic Δ, which characterizes the permissible vibration level of the
turbocharger rotor.
Preliminary experiments on MAN MC diesel engines have shown that an increase in the amplitude of the
harmonic at the main frequency υ turbocharger in 2-3 times regarding the average level of the amplitude
spectrum characterizes the dangerous vibration level of the turbocharger rotor. The average level of the
harmonic amplitudes was estimated in the frequency range
To better quantify the limits of vibration level of the turbocharger rotor, further research is required. It may
be noted that the spectrum analysis of vibroacoustic signals of the turbocharger compressor can be made
quickly under operating conditions.
6.0 CONCLUSIONS
The methods of vibrodiagnostics of marine diesel engines with turbocharging considered in the article can be
helpful for practical use. The “leakage effect” method improves the reliability of diagnostic findings.
A vibroacoustic method for determining the speed of a turbocharger rotor and estimating the level of the
oscillation amplitude at the main rotational speed can be used as a basis for the express diagnostics of the
turbocharger under operating conditions.
9-12 STO-MP-AVT-306
Vibration Diagnostics Methods of Marine Diesel Engines with Turbocharger
7.0 REFERENCES
[1] Zigelman E.B., Skvortzov D.F., Loshinin I.A. Study of Possibility for vibrodiagnostics of medium
diesel generators. Izvestiya vuzov, 2013, N6, pp.-42-48.
[2] ISO 10816 Series. Mechanical vibration -- Evaluation of machine vibration by measurements on non-
rotating parts.
[5] Suri Ganeriwala, Phd. Review of Techniques for Bearings & Gearbox Diagnostics. IMAC Conference
- Feb. 3, 2010, Jacksonville FL.
[7] Naumenko A. P. Real-time condition monitoring of reciprocating machines. – In: The 6th Intern.
Conf. on Condition Monitoring and Machinery Failure Prevention Technologies. – Dublin (Ireland):
2009, pp. 1202–1213.
[8] Varbanets R.A. Turbocharged Marine diesel engine frequency parameters monitoring / R.A. Varbanets,
Y.M. Kucherenko, A.І. Halavan // Bulletin of the Astrakhan State Technical University. Series: Marine
equipment and technology. "Astrakhan. - 2013. - № 1.- pp 103-110.9. Varbanets R., 9.
[9] Varbanets R.A., Karianskiy A. Marine diesel engine performance analyze of // Journal of Polish
CIMAC. Energetic Aspects. Gdansk: Faculty of Ocean Engineering and Ship Technology Gdansk
University of Technology, 2012, vol. 7, no. 1, pp. 269–275.
[10] Varbanets R.A. The “leakage” elimination methods analysis of the spectrum in the diagnosis system of
marine diesel engines turbo system / R.A. Varbanets, Y.N. Kucherenko, A.I. Golovanov,
N. Alexandrovskaya // Artificial intelligence. - 2013. - № 4 (62). - pp. 289-295.
[11] Varbanets R.A. Turbocharged Marine diesel engine frequency parameters monitoring / R.A.
Varbanets, Y.M. Kucherenko, A.І. Halavan // Bulletin of the Astrakhan State Technical University.
Series: Marine equipment and technology. Astrakhan. - 2013. - № 1.- pp. 103-110.9.
[12] Applied Time Series Analysis, by Robert K. Otnes and Loren Enochson.New York: Wiley, 1978, pp.
428.
STO-MP-AVT-306 9-13