Engr Moiz Assignment Video
Engr Moiz Assignment Video
Methods of ConditionMonitoring
and Fault Diagnosis for Induction Motors
0.V. Thorsen, M. Dalva
Abstract
The study gives a synopsis over condition monitoring methods both as a diagnostic tool and as a techniquef o r
failure identification in high-voltage induction motors. Methods of fault detection are reviewed and analyzed:
mechanical (vibrations, shockpulses, acoirstic, speedfluctuations), electro-mechanical (currents, surges, par-
tial discharges, leakagejluxes) and temperature-, oil particle-, gas analysis, and performance- and visual me-
thods. The study indicates how different c p e s of faults generate different frequency components in vibration
spectra, in motor supply current, in torque and in magnetic field, and relationships between fault types and
spectra have been developed. Supply current waveform analysis has been applied to demonstrate the detec-
tion of broken rotor bars and bearing faults.
Here are: fect parts. a crack on the innedouter race or on the roll-
h number of balls or rollers. ing elements in the bearing. A crack on the stationary
race gives pulses of the same strength, a crack on the ro-
fi contact angle of the balls on the races. tating race gives impulses that will vary with changes in
jnl, relativespeed(in s-I) between innerandouterraces,
ball diameter, the rolling-element load [ 7 ] .
bearing pitch diameter.
3.4 Acoustic Emission Monitoring
There will alsobevibrationsof very high frequencies.
The test machine has the following bearing data: Acoustic emission monitoring works with ultrason-
Shaft end bearing: ic as well as audible frequencies. Contact between roll-
b = 9, dh = 9.52 mm, d, = 53. I mm. ing elements with and without cracks generates waves
Fan end bearing: that propagate through the machine with the speed of
b = 9, d,, = 7.94 mm, dp= 45.25 mm. sound. The waves have little energy, but their high fre-
quencies can be detected by piezoelectric or capacitive
Nominal rotor speed: I 435/60 = 23.92 s-'. transducers. The method is particularly promising for
Calculating the frequencies for the bearings for test the detecting of bearing faults.
motor:
- Shaft end bearing: j ,= 88.3 Hz, J = 126.9 Hz. 3.5 Speed Fluctuations Monitoring
fh = 129. I H Z , fc.= 9.8 HZ.
Speed fluctuations monitoring can detect defects by
- Fan end bearing: ji, = 88.8 Hz, f;= 126.5 Hz. measuring fluctuations in the rotational period of the
f h = 132.3 Hz, = 9.9 Hz. motor. The method is particularly useful at rotor faults.
Unhtrltrncetl rotor can be detected by the increas- vibrations, air-pap eccentricity, rotor asymmetry, dam-
ing of vibrations with rotational frequency (n/60).or aged bearings/couplings and misalignment.
twice. three or four times that frequency. Misalignment The current in a normal rotor bar will fluctuate si-
can be detected by the increase of vibrations of twice nusoidally with slip frequency, and give a contribution
the rotational frequency of the shaft 171. Since magnet- to thedeveloped torque that will vary sinusoidally with
ic forces are proportional to the square of the flux den- twice the slip frequency. A defective bar will not con-
sity, vibrations will come into existence, also at normal tribute to the shaft torque. The shaft will consequently
run, with the frequency 25, and slot harmonic compo- slow twice per slip cycle. Summed up for the whole
nents as shown in the following ( , f ,is the supply fre- rotor. one can say that the resulting torque consists of
quency). two components. one constant and one that varies with
In the core vibmtion specrriim the rotor vibration twice slip frequency. Provided all rotor bars are in
slot harmonics will be:f; [ ( Z ? / p () 1 -s) f v ' ] . v ' = 0,2,4, order, the mean ofthe variable part ofthe torque is zero.
6, ...,and sidebands+2sfi. Broken rotorburs will mod- However. in most cases the motorhas a loadof variable
ulate stator core vibrations with f2sfi. torque, and the measuring instruments must be capable
Statorasymmetr_v(winding fault or unbalanced sup- of distinguishing between load tluctuations and fluctu-
ply/single phasing) produces vibration signals contain- ations of twice slip frequency indicating rotor faults.
ing 25, for asymmetric rotor also twice slip frequency The selective measurement therefore requires a refer-
side band ( f 2 s f i ) on the shaft frequency [ 7 ] . ence signal at slip frequency. [8] describes principles
Static eccentricity gives frequency components of of instruments capable of detecting faults on a single
twice supply frequency (2J ) in addition toslot harmon- rotor bar of a squirrel cage motor by speed fluctuation
ic components:fi [ ( v Z l ) (I - s ) / p f v ' ] , v' = 0. 2,4, ..., measurements.
v = 1,2,3, ...; Z? is the number of rotor slots andp is the Determination of slip and speed can be performed
number of pole pairs. indirectly by isolating the fundamental slot harmonic
Dynamic eccentricity also produces twice slip fre- frequency component from the supply current waveform
quency side bands on the shaft frequency in addition to [9]. With siniisoidal mmf the frequency of the slot har-
the slot harmonic components: fi [ ( v & -t nd)(1 - s)/p monic component for the stator current can be ex-
f v ' ] , n d = I , 2, ... [4,7]. pressed:
An alternative signal processing method for vibra-
tional analysis of bearing monitoring is the statistically Z2(1 -S,fi/Pffi.
treatment of the measuring data from accelerometers. A To take stator mmf time h a h o n k s into considera-
bearing that runs satisfactorily gives a Gaussian distri- tion, f i must be multiplied by h; the slot harmonic fre-
bution of the frequency components. Defects in the bear- quency for the stator current will then be:
ing disturb this distribution. This method gives an abso-
lute measure of the bearing condition, one does not need Zz(1 - s ) f i / p * h f i .
to know the original condition. h being the harmonic of the stator mmf time harmonics.
h = l , 3 , 5 , 7 .....
3.3 Shock-Pulse Monitoring The motor speed is given by:
n = 60fi ( I - s)/p.
Shock-pulse monitoring (SPM) applies piezo-
electric transducers fixed to the bearings to detect shock where the slip s is found from the measurement of the
waves caused by impact between moving parts and de- frequence components.
ndis an integer (eccentricity order number 0 for static ec- b) n 50Hz (1-2s) dB Diff. (1+2s) Diff.
centricity and 1,2,3,4, ... for dynamic eccentricity) [7]. mid dB .50Hz dB .50Hz dB dB
1435 -0.2 45.7 -48.3 48.1 54.3 -55 54.8
3.6 Current Monitoring 1454 -1.5 46.9 -51.1 49.6 53.1 -56.1 54.6
1471 -3 48.1 -55.4 52.4 51.9 -57.1 54.1
Current monitoring deals with the analysis of the C) n 50Hz (1-2s) dB Diff. (1+2s) Diff.
motor supply current because this is influenced to agreat m i d dB .50Hz dB .50Hz dB dB
extent by the motor condition. Current analysis may re-
I435 -0.4 45.7 -38.2 37.8 54.3 -42.7 42.3
veal damaged rotor bars as well as mechanical problems I454 -1.8 46.9 -43.1 41.3 53.1 -44.8 43
associated with rotating machines (MCSA - Motor Cur- 1471 -3.1 48.1 -49.8 46.7 51.9 -49.5 46.4
rent Signature Analysis) [ 101.
Faults in rotor bars will introduce different frequen- d) n 50Hz (1-2s) dB Diff. (1+2s) Diff.
cy componentsin the supply current to the stator. Because min-' dB .50Hz dB .50Hz dB dB
of asymmetry in the rotor circuit as a result of damaged I435 -0.1 45.7 -34 33.9 54.3 -38.6 38.5
rotor bars the rotor current will produce positive- and 1454 - 1.6 46.9 -39.3 37.7 53.1 -39.5 37.9
negative-sequencerotormmf at frequencies sol. Here a 1 1471 -3.1 48.1 -45.2 42.1 51.9 -43.6 40.5
is theangular frequency of the supply fundamental.
One assumes as a starting point that the mmf is sinu- Tab. 1. Overview of the results by measurements on the test
motor (Diff.: Difference)
soidally distributed. Positive rotor mmf experiences an a) Sound motor
angularspeedsw, + wl(1 -s)= a1relative tostator, while b) One broken rotor bar
angular speed of negative sequence rotor mmf relative ro c) Two broken rotor bars
stator is soI- a,( 1 - s) = (2s - I)al,which gives an in- d) Three broken rotor bars
duced stator current of the angular frequency (2s - I )al.
If we considered the v-th space harmonic we will
have:sol +.vwI(1 -s), v = 1,3,5,7, I1,13, ... [ I I ] being In asimilar way one can describe the next components.
the order of space harmonics. This gives the following One can conclude by sayingthat damaged rotor bars
frequencies for different values of v: WI, (2s - I)wI, and generally rotor asymmetry will produce the follow-
(3- 2s)Wl, (4s - 3)a1,( 5 - 4s)ol, (6s - 5)a1, (7 - ~ S ) O I , ing frequency components in the srator current:
( 8 ~ - 7 ) ~...
l , [7, 121.
( 2 s - l ) 0 1 , ( 4 ~ - 1 ) ~ 1 , ( 2 sl)wl,(6s-
+ I)o1,(4s+I)w1.
Between the angular frequencies of the positive- and
negative-sequencestator currents components there will These are sidebands f2 s f , , 2 4sfi, f6sfi, ... around
hence be an angular frequency difference of 2 ~ 0 1that , the supply frequency fi [lo]. For the test motor this
again leads to rorque pulsations of the angular frequen- means: 45.7 Hz, 54.3 Hz,41.3 Hz,58.7 Hz,37 Hz,
cy 2 s a l and their accompanying speed fluctuations 63 Hz, ... for 1 435 min-I.
2 ~ 0 The
1 . speed fluctuations relative to stator is: 2 s w 1 + The magnitudes of the above described frequency
0l(1 -s) = (S + I)WI and 2~01- W I ( 1 -s) = (3s- I)Ol. components are small compared with the main supply
The stator flux with its electrical angular frequency wI frequency component and are decreasing with the har-
therefore produces additional rotor currents with angu- monic order, and therefore difficult to detect. Particular-
lar frequencies swI and 3so1. 3sal produces positive- ly at low load or no load, these frequency components
and negative rotor mmf that relative to stator has an an- are close to the mains frequency, and that also leads to
gular speed: 3swl - al( -
1 - s) = (4s I)wI, and 3s01 measuring complications.
+ OI(1 - s) = (2s + l ) o l which induce stator currents of Tab. 1presents an overview of the results by meas-
the angular frequencies (4s- l ) w l and (2s + 1)ol. If we urements on the mentioned test motor. First measure-
considered the v-th space harmonic we will have: ments (a) are performed on a sound motor, thereafter on
a motor with one (b), two (c) and three broken bars (d),
3SOl 2 V O l ( 1 -S).
respectively. For each case the measurements are carried
New torque pulsations get the angular frequency out with a load of 100 % (1 435 min-'), 75 % (1 454 min-I)
4sw1,and associated angular speed fluctuations referred and 50 % ( 147 1 mid), respectively,Fig. 1to Fig. 4 show
to the stator: 4 s a 1+ al(1- s) = (3s + I)wI, and 4 s a 1 the current spectra.
- 01 ( 1 - s) = (5s - l)ol.The stator flux produces new [ 131 applies a method for statistical signal process-
rotor currents of the angular frequency 5 ~ 0 1 . 5 ~pro- 0 1 ing for detecting rotor faults on-line.
duces positive- and negative-sequence rotor mmf that Statistical analysis can be used for signal process-
again will cause stator currents with angular frequencies ing, often as a supplement to the direct method.
(6s- l)wl and(4s+ l)al. If weconsidered the v-thspace A current component with a frequency lower than
harmonic we will have: the main supply frequency by twice the slip frequency
will cause an analog ammeter in the supply lead to
fluctuate slowly. A spectrum analyzer will provide
I
i
gT
01
5 I5 25 35 45 55 65 75 85 Hz 105 5 15 25 35 45 55 65 75 85 Hz 105
f--’ f-
Fig. 1. Current spectrum for motor with one broken rotor bar Fig. 3. Current spectrum for motor with three broken rotor bars
(black)together with that of a sound motor (grey), 100 % load (black) together with that of a sound motor (grey), 100 % load
mT
01
21
5 15 25 35 45 55 65 75 85 Hz 105
f -*
Fig. 2. Current spectrum for motor with two broken rotor bars
(black)together with that of a sound motor (grey), 100 % load
5 15 25 35 45
f -.
55 65 75 85 Hz 105
Fig. 4. Current spectrum for motor with three broken rotor bars
(black) together with that of a sound motor (grey), 50 % load
f i [ I ~ m ( l - s ) / p ] , m=1,2,3,... [lo].
This can be the same frequencies as given by certain 205 225 245 265 285 305 325 345 365 Hz 405
unsymmetries and faults in rotor bars. f b
For a two-poles motor the following frequencies are Fig. 6.Current spectrum for motor with an outer race fault on
given: sfi, (2- s)fl, (2s - I )fi, (3 - 2s)fi,(3s - 2)fi, both bearings and cage fault on one of them (black), together
(4- 3s)fi,(4s- 3)fi,(5- 4s)fi, form = 1,2,3. with current spectrum for sound motor (grey), 100 % load
For a four-poles motor the following frequencies are a) Frequency range 5 Hz ..,205 Hz
given: ( I -s)fi/2, (3-s)f,/2,sf!,(2-s)fi,(1 - 3s)J/2, b) Frequency range 205 Hz ..,405Hz
(5 - 3s)fi/2,for rn = I , 2,3.
A defect in rolling element bearings causes the air
gap of the motor to vary, and generates stator currents at fan end 138.8 Hz,38.8Hz,227.6 Hz,127.6Hz etc. For
frequencies: an inner race defect the frequencies will be 176.9 Hz,
76.9 Hz, 303.8 Hz,203.8 Hz, ..., 176.5 Hz, 76.5 Hz,
ffikmJ.1, m = 1,2,3,... 1141. 303 Hz,203 Hz etc. For a ball defect the frequencies will
f v is one of the vibration frequencies for the bearings as be 179.1 Hz,79.1 Hz,308.2Hz,208.2 Hz, ..., 182.3Hz,
given before, calculated from the bearings dimensions. 82.3 Hz,314.6Hz. 214.6Hz,..., and for a cage defect
To the test motor several bearing faultswere imposed. 59.8Hz,40.2Hz,69.6Hz,30.4Hz ,...,59.9Hz,40.1Hz,
An outer race defect on the shaft end of the motor gener- 69.8Hz,30.2 Hz etc.
ates currents of frequencies: (50 2 88.3)Hz = 138.3 HZ Fig. 6 presents current spectra with fault on outer
and 38.3 Hz,226.6Hz and 126.6 Hz, ..., Fig. 5, on the race for both bearings, and cage defect for one of them,
together with current spectra for sound motor, 100%
load.
Air-gap dynamic eccentricity, caused by bearing
faults, produces anomalies in the air-gap flux density
and thereby harmonic current components as shown.
tI
3.7 Air-Gap Torque Monitoring
I T
$1 ber of PD pulses (above a specified magnitude) per sec-
ond measured over a selected time. Integrated quantities
are current I, qiiadrntic rate D and discharge power P . I
is the sum of the absolute values of the apparent charg-
es during a certain time intervai divided by this time
interval. D is the sum of the squares of the apparent
charges, during a certain time interval, divided by this
+ & time interval and Pis the average power fed into the ter-
Fig. 7. The abc model for partial discharges in a dielectric minals of the test object due to PD. Alternating voltages
with a void related to PD are given by their peak values divided by
a) Test object b) Equivalent circuit 42. Inception voltage U,is the lowest voltage at which
partial discharges are observed in the test arrangement,
when the voltage applied to the test object is gradually
age. The little internal movement of charge is nearly a increased from a lower value at which no (or specified
Dirac pulse. that is a time-current graph of nearly zero low values of) such discharges are observed. PD extinc-
time but a defined area. The voltage drop UAacross the tion voltage U,is the lowest voltage at which PD are ob-
void caused by the Dirac pulse, causes a voltage drop in served in the test arrangement when the voltage applied
the outer circuit that can be expressed as a voltage divi- to the object is gradually decreased from a higher value
sion between c b and C,, since one can suppose that dur- at which such discharge are observed.
ing the extreme short time of this process, there will be
- Alternative t o the abc Model
no current in the outer circuit. However, the voltage drop
over C,is detectable across the terminals oftheelectrodes The abc model has been successfully used for many
as UA = u ,cb/(c b + ca). years despite many shortcomings. Strictly speaking, the
Immediately after the Dirac pulse is over, a current concept of void capacitance is not correct, and the other
will start to flow in the outer circuit. This current is also capacitances involved cannot be calculated [23,24].By
of short duration, but lasts much longer than the Dirac using classical electrostatic theory one can find expres-
pulse. This current represents a charge 9 which can be sions for apparent charge and current due to PD for a die-
expressed as: q = uA(C, + CbCc)/fCh+ C,) or nearly: lectric with given defect geometries. It will lead too far
q = u A C ~ The
. charge q is defined apparent charge. AC- to describe the method more precisely, but one of the
cording to [20] the apparent charge of a partial discharge conclusions is that a defect (a void or another occlusion)
is that charge which, if injected instantaneously between which causes a PD with a certain charge in one system
the terminals of the test object, would momentarily will produce a PD magnitude of half that value in a
change the voltage between its terminals by the same system of double size but with the same defect.
amount as the partial discharge itself. The apparent
charge is expressed in pC, and is not equivalent to the - Reasons and Locationsfor PD
charge qlnrat the PD site, that is, in the above description, The insulation of a modern high-voltage rotating AC
q is different from q,",.Let the voltage across the test ob- machine stator winding is mainly epoxy/mica (former-
ject be: u = Usin(wt). Then the voltage across the ly asphalt/mica) groundwall insulation covered with a
voidpriortothePDis: uI = [ c b / ( c b + C c ) ] f l Usin(wt). semi-conductive layer. and stress grading coating. The
These two voltages are drawn in Fig. 8. The voltage at bars are normally pressed against the bottom of the slot
which the PD starts is U,(inception voltage, see below), by wedges at the top of it. By casting or extruding insu-
and the rest voltage just after the PD is Urns,. lation materials it is in practice nearly impossible to
Referring to Fig. 8, when uI at the instant tl has completely avoid small occlusions of air or gas (voids,
reached U,, the breakdown starts, and the voltage across cavities) in the solid. Insulation systems build up of
the void drops to Umt,the differencebeing: UA= Ui- URyl. layers of different or equal materials will often have
Immediately thereafter the voltage across the void will
build up according to the equation u2 = U I - UA(that is,
according to the sloping part of the nearly saw-tooth part
of the graph (parallel to ul). Normally u2 will reach U,,
and another PD will occur (t2), and the voltage across the '1 2
void will increase according to u 3= u I - 2 UA. Now it may u3
well happen that 1.43 will reach U,, and hence another sin-
gle PD starts. However, in the example of Fig. 8 it does +
not do so until it has become negative, and the single PD t
processes starts over again with negative polarity ( t 3 ) .
The whole PD process is characterized by a series of sin-
gle PD in the first part of the positive sine wave, and a
similar series in the first part of the negative part of the
sine wave. In accordance with the particular cases, there Fig. 8. Test voltage and voltage across a void during PD
may be more or less single PD in each series. in a dielectric
small slips between the layers, or there may be small par- Impedance or filter
0 -
ties where the insulation does not adhere perfectly to the T I
metallic electrodes, or a pre-insulated bar of a stator Coupling -L
winding may not be side-packed properly in the slot, or +Test object capacitanceT
the wedge may have become loose. As a result of poor 1 Connecting lead I
side-packing and/or loose wedges the electromagnetic
forces under load will cause small movements that may
abrade the field-equalizingpaint in the slot section of the
coil. Further one can often observe PD at the end wind-
ings, mainly because of dirt or oil contaminations. PD
may also occur around a conductor insulated with a gas Fig. 9. Example of a test circuit; the dotted Z, is an alterna-
(corona). tive position [20]
We may sum up here that PD may take place:
in voids in the dielectric itself in the slot part,
pressor for noise from the mains. The so-called Phase
in voids between the copper conductors and the Resolved Partial Discharge Analysis (PRDA) applies
groundwall insulation, very fast responding computer techniques so that it has
been made possible to analyse the PD as to the number
in voids between the surface of the stator insulation
of discharges, the charge in each of them, and, if supplied
and the core iron (slot discharge), and
with other sensors in addition to the RFCT, the location
in voids in the dielectric itself in the end winding, or of the PD in the motor and the type of fault that have led
surface discharges there due to contamination or pure to PD. In addition to the principle that noise pulses from
design as to clearances or stress grading (altogether the mains will be suppressed by cables or condensers,
called end winding discharges). one applies the principle that PD pulses and noise puls-
es have given different distances to travel to the sensors.
Consequences of PD
PD will deteriorate the insulation in three ways [ 2 2 ] : - Calibruting PD Measuring Instruments
Calibration is a major concern with PD. To date it has
Bombardment of the insulation around the discharge not been possible to establish an accepted and reliable
spots from ions and electrons accelerated in the dis- method
of absolute measurement so that one for instance
charge paths. could design a standard transportable calibrating instru-
Chemical reactions in materials adjacent to the PD ment. That means among other implications,that it is dif-
spots caused by the PD, particularly organic materi- ficult to establish a sort of standard recommended limits
als are subject to these reactions. for the “severity of the PD’. The calibration that has been
possible to be established, is only valid for a particular
Radiation from PD (PD are extremely fast phenom-
machine with particular instrumentation. However.
ena, pulse rise time and duration are in the order of
much research is currently going on in this field, and even
ns); ultraviolet light may have enough energy to de- at the present level the PD monitoring has proved to be a
teriorate organic compounds.
very useful. not too expensive, and possible to be run by
PD will erode insulating materials such as polyester non-experts after a short training period. PD quantities
and epoxy resins. Mica on ;he other hand is very resist- that may directly be recorded and displayed by a PD in-
ant against attack from PD. The compound of epoxy strument arepuise-repetition rate np(pulse count rate) in
/mica is therefore relatively long living as far as PD is the unit “pulses persecond”,npparentchargeq in the unit
concerned, so that it will be necessary to watch the trend pC (or pulse magnitude in the unit mV) andpulse-phase
of the deterioration of the insulation in order to predict angle in the unit degree. The three quantities may be pre-
the condition in the future. sented in the form of a three-dimensional graph that dis-
tinguishes between positive and negative pulses. A less
- Methods of Detecting PD complicate variant presents pulse count rate and magni-
Until some few years ago detecting and measuring tude with polarity in a two-dimensional graph.
of PD was performed off-line, and Fig. 9 shows a basic For many years the insulation of high-voltage ma-
test circuit [20,25]. The impedance (or filter) prevents chines was based on mica flakes and bitumen- or shellac
discharge pulses from being bypassed through the high resins. As the voltage and power ratings increased, one
test-voltage supply, it also will attenuate disturbance experienced localized premature failures. Traditional
(noise) from the supply. The recording instrument is a methods for monitoring the bulk insulation condition, e.
very fast responding oscilloscope or a recorder. Much g. the “Loss Tangent Variation”, showed to be insuffi-
effort is placed on designing such test circuits for dif- cient to detect localized premature failures in a reliable
ferent test object and different noise condition and other way. And when the more modem thermoset insulation,
demands. based on mica chipping and glass fibres, was introduced,
For measurements of PD in high-voltage machines the need for detecting localized faults was still there. The
on-line, high (radio) frequency current transformers localized faults (attack on one coil or conductor) on this
(RFCT) are clamped on the machines supply cables. For modem insulation can appear unexpectedly and devel-
machines supplied by a bus duct system, the sensoring op relatively rapidly, so one has to follow the situation
is more complicate. The cable itself serves as noise sup- with not more than months interval.
Temperature measurements with thermocouples Visual monitoring means either inspection with the
and thermistors have a long tradition. Bearing tempera- bare eye or closed circuit television cameras, or a by pho-
tures give useful information on possible friction prob- tography or video tape to record trends.
lems, and are surveyed on a routine basis. Coolant bulk
outlet temperature is often monitored, particularly when
the machine is stressed beyond its rated data. 4 Application of Monitoring Methods
Winding temperatures on the hottest spots are difi-
cult to detect because of insulation problems, especially Use of monitoring methods are continual increasing
in high-voltage motors. May be a fibre-optic system is for medium to large induction motors. Both off-and on-
an answer to that problem. line systems are in use.