Basic Vibration Signal Processing For Bearing Fault Detection
Basic Vibration Signal Processing For Bearing Fault Detection
1, FEBRUARY 2003
149
AbstractFaculty in the College of Engineering at the University of Alabama developed a multidisciplinary course in applied
spectral analysis that was first offered in 1996. The course is
aimed at juniors majoring in electrical, mechanical, industrial,
or aerospace engineering. No background in signal processing
or Fourier analysis is assumed; the requisite fundamentals are
covered early in the course and followed by a series of laboratories
in which the fundamental concepts are applied. In this paper, a
laboratory module on fault detection in rolling element bearings is
presented. This module is one of two laboratory modules focusing
on machine condition monitoring applications that were developed
for this course. Background on the basic operational characteristics of rolling element bearings is presented, and formulas given
for the calculation of the characteristic fault frequencies. The
shortcomings of conventional vibration spectral analysis for the
detection of bearing faults is examined in the context of a synthetic
vibration signal that students generate in MATLAB. This signal
shares several key features of vibration signatures measured on
bearing housings. Envelope analysis and the connection between
bearing fault signatures and amplitude modulation/demodulation
is explained. Finally, a graphically driven software utility (a set
of MATLAB m-files) is introduced. This software allows students
to explore envelope analysis using measured data or the synthetic
signal that they generated. The software utility and the material
presented in this paper constitute an instructional module on
bearing fault detection that can be used as a stand-alone tutorial
or incorporated into a course.
Index TermsBearings, envelope analysis, fault detection, kurtosis, spectral analysis.
I. INTRODUCTION
150
TABLE I
TOPICS COVERED IN THE FUNDAMENTALS MODULE
Fig. 1.
very high speeds, and robust in operation and require little maintenance. Angular-contact ball bearings have an angle between
Fig. 2.
151
the inner and outer rings, as shown in Fig. 2, which enables them
to support both radial and axial loads. This type of bearing is
used in the bearing fault simulator in the laboratory module, but
the basic ideas and equations developed for these bearings can
be readily adapted to roller bearings.
B. Bearing Failure Mechanisms
There are a number of mechanisms that can lead to bearing
failure, including mechanical damage, crack damage, wear
damage, lubricant deficiency, and corrosion. [3] Abusive
handling can induce nicks and dents, which are especially
harmful when located in areas tracked by the rolling elements.
When the smooth rolling contact surfaces are marred, higher
stress conditions imposed on the surface reduce bearing life
significantly. Permanent indentation created by rolling element
overload is called brinelling. A crack in a bearing component
may begin as a manufacturing-related defect or be induced by
operating stress via overload or cyclic loading.
Wear results in a gradual deterioration of the bearing
components, leading to a loss of dimensioning and associated
problems. When lubrication conditions become inadequate, the
increased friction results in metal-to-metal contact. Operating
forces can cause enlarged plastic deformation by tearing the
locally friction-welded regions from the metal matrix. Lighter
adhesive damage is often called scuffing or scoring; whereas,
more intensive damage is referred to as seizing or galling.
Abrasive wear occurs when hard particles become entrained
between the contact surfaces.
As discussed in connection with wear damage, above, insufficient lubrication accelerates the adhesive wearing progress.
Poor lubrication increases bearing component temperatures,
which speeds up the deterioration processes. Bearings that
operate in an environment of high humidity may be subjected
to surface oxidation and produce rust particles and pits. These
particles can produce rapid wear via abrasion. The pits function
as stress raisers, providing sites for crack initiation.
and
(3)
152
where
is the ball diameter,
the pitch diameter,
the
number of balls, and the shaft rotation rate in hertz. These
formulas are theoretical, and discrepancies arise when bearings
carry significant thrust loads or there is slippage. Software and
software modules (sold as options with vibration condition
monitoring packages) that incorporate databases of parameters
for a wide variety of bearings are commercially available.
III. BASIC BEARING FAULT DETECTION TECHNIQUES
If vibration signatures measured on bearings were as uncomplicated as that in Fig. 3, detection of bearing flaws would be
an easy matter. In practice, however, vibrations measured on a
bearing are dominated by high-level imbalance and misalignment components and include random vibrations associated
with friction and other sources. Imbalance vibration occurs at
the shaft rate of rotation (referred to as the 1X); misalignment
shows up at the fundamental (1X) and its harmonics. The
spectral components associated with the ringing pulse sequence
are not integer harmonics of the fundamental, and would not
be mistaken for them, but are of relatively small amplitude.
In a conventional spectrum, these components are lost in the
spectral noise floor generated by random vibrations and leakage
from the high level harmonics.
A. Synthetic Signal
In the laboratory module, students construct a synthetic
signal that shares some of the characteristic features of the vibration signatures measured on bearings. This signal is synthesized
as
a) A five-term approximation to a 0.05 amplitude 6-Hz
square wave [see Fig. 5(a)]. A square wave contains only
odd harmonics, thus this approximation contains only the
fundamental (1X), 3X, 5X, 7X, and 9X harmonics. An
imbalance generates a high 1X level, and a misalignment
peaks, but the absence of
generates
even harmonics in the square wave does not affect the
instructive value of the synthetic waveform. The peak
amplitude of the waveform in Fig. 5(a) is 0.059; the rms
amplitude is 0.049.
b) Random noise with a Gaussian amplitude distribution
[see Fig. 5(b)]. Random noise is present in all measurements and arises from a variety of sources (background
vibrations, rubs, electrical noise). Random noise with
a Gaussian amplitude distribution is readily generated
using the randn command in MATLAB. Here, the noise
was scaled so that in Fig. 5(b) the maximum amplitude
is 0.5 and the rms amplitude is 0.11.
c) A ringing pulse sequence with a repetition rate of 40 Hz
(6.67X) and a 0.1 peak amplitude [see Fig. 5(c)]. A pulse
sequence generated by appending multiple records of a
single measured 4096 point pulse (the response of a table
to an impact) is used. This pulse sequence models the
resonant response of a highly damped bearing housing as
ball rolls over a race fault.
The individual signal components are shown in Fig. 5 along
with their sum the composite signal. A MATLAB m-file, syn-
Fig. 6. Spectra of the (a) Square wave approximation. (b) Random noise.
(c) Ringing pulse sequence. (d) Composite waveform in Fig. 5.
153
Fig. 8. Application of the Hilbert transform to two time waveforms. a) Modulated sine wave. b) Pulse sequence from Fig. 1(c).
154
Fig. 9.
(1)
This signal can be viewed as the real part of the complex signal
Fig. 10. Two signals with different kurtosis coefficients. The kurtosis
coefficient of the signal in (a) is larger than that of the signal in (b).
The kurtosis coefficient is a metric that is also used in the detection of bearing faults. It is calculated from the time domain
data. [10][13] The variance, , of a zero mean (no dc component) signal is the mean square value [i.e., the square of the
root-mean-square (rms) value]. AC coupled accelerometer signals are zero mean signals. A rectified waveform has a nonzero
average or mean value. In this case, the variance is given by
(3)
where is the number of data points in the data sequence, and
is the average value. The coefficient of kurtosis is a fourth-order
statistic normalized by the square of the variance
D. Envelope Analysis
The phrase envelope analysis typically refers to the
sequence of procedures indicated in Fig. 9. The purpose of
bandpass filtering is to reject the low-frequency high-amplitude
signals associated with imbalance and misalignment and to
eliminate random noise outside the passband. [8], [9] Analyzers and data loggers that offer envelope analysis generally
offer user selectable bandpass settings, such as 12.5 kHz,
2.55 kHz, 510 kHz, 1020 kHz, and 2040 kHz. After
bandpass filtering and rectification (here, using the Hilbert
transform), the final step in the envelope analysis process is
calculation of the spectrum of the rectified bandpassed signal.
(4)
Large values of indicate large excursions from the RMS value.
The physical meaning of kurtosis may be best explained by illustration. Fig. 10(a) shows a signal with kurtosis value that is
large, relative to that of the signal in Fig. 10(b). This statistical
moment is sensitive to the pulses induced by bearing defects and
has proven effective in the rolling element bearing diagnostics.
IV. INSTRUCTIONAL SOFTWARE UTILITY
In this section, the use of the software utility developed for
this module is described. The utility is invoked from within
Fig. 11.
Fig. 12.
155
156
V. CONCLUSION
S. A. McInerny received the B.S. degree in mechanical engineering from California State University, Long Beach (CSULB), in 1979, and the M.S. degree in
engineering and the Ph.D. degree in mechanical engineering from the University of California, Los Angeles, in 1984 and 1987, respectively.
She was a Member of the Technical Staff at The Aerospace Corporationfrom
1987 to 1993 and an Associate Professor of Mechanical Engineering at CSULB
from 1988 to 1993. She conducts applied research in fault detection and diagnostics in mechanical systems and has spent the last several summers working
in the propulsion diagnostics group at the Naval Air Warfare Center, Patuxent
River, MD, on bearing and gear diagnostics.