Jeas 0912 776
Jeas 0912 776
Jeas 0912 776
9, SEPTEMBER 2012
ISSN 1819-6608
www.arpnjournals.com
ABSTRACT
Three Phase Induction Motors are widely used for industrial and domestic applications. There are various faults
that occur in induction motors like stator inter-turn fault, bearings faults and eccentricity fault. Out of these faults, the
rotor broken bar fault is very specific in squirrel cage induction machines. This paper deals with the detection of broken
bars in three phase squirrel cage induction motor using finite element model of the induction machine. Finite element
method is more precise than the winding function approach method, as it is based on the actual geometry of the machine.
This paper uses a CAD package called Magnet for the Static 2D and Transient 2D analysis. The various machine
parameters like flux density, flux function, magnetic energy, etc are calculated using this CAD package and their values are
compared under healthy and faulty conditions.
Keywords: rotor broken bars, finite element method, induction motor, magnet, flux density, flux function, magnetic energy.
INTRODUCTION
Three phase induction motors have wide
applications in various industries. It is recognized that a
variety of faults can occur in these motors during normal
operation such as rotor fault (broken bars or end ring),
stator inter-turn fault, eccentricity fault and bearing fault.
Hence early detection and diagnosis of such faults are very
essential for the protection of induction motors against
failures and permanent damages. A sudden motor failure
may reduce productivity and may be catastrophic in an
industrial system if undetected. In recent years, the
problem of failures in large machines has become more
significant. The desire to improve the reliability of the
industrial drive systems has led to researches and
developments in various countries to evaluate the cause
and consequences of various fault conditions.
The interior faults of induction motor accounts
for more than 70% in proportion of induction motor
failures. Interior faults include stator and rotor faults of
induction motors. Rotor faults are related to broken bars.
Rotor failures are caused by a combination of various
stresses that act on the rotor and these stresses can be
identified as electromagnetic, thermal, dynamic,
environmental and mechanical. Therefore these leads to
low-frequency torque harmonics, which increases noise
and vibration.
The transient analysis was done using coupled
electric circuit with 2D finite element electromagnetic
field analysis. The designed geometric dimension is
modeled in the finite element domain and transient
performances are predicted at the starting of motor with no
load [1]. An MCSA technique is adopted for the diagnosis
of rotor breakages in squirrel cage induction motor and
finite method to calculate the parameters and modeled
using state space modeling approach. It computes the
characteristic frequency components which are indicative
of rotor bar and connector breakages and developed torque
profiles [2].
1170
ISSN 1819-6608
www.arpnjournals.com
numerical approach which uses 2D, nonlinear, time
stepping finite element method for excitation from a
constant voltage source [12].
The effect of pole pair and rotor slot numbers was
presented under different harmonics in healthy and
eccentric conditions. The simulation technique was not as
accurate as Finite Element Method [13]. A dynamic model
for IM under inter turn insulation failure fault was derived
using reference frame theory. Finite element analysis is
used for parameter determination of the machine in
healthy and faulty condition [14]. An MCSA technique to
diagnose the faults in the three phase induction motor
drives was focused [15]. The use of Partial Relative
Indexes (PRI) is proposed as a new fault indicator to
ameliorate the reliability of fault detection task and uses
MCSA method [16].
The above literature uses vibration monitoring
techniques, MCSA, and Thermal Monitoring for the
detection of faults. In this paper Finite element analysis is
adopted to perform Static and Transient 2D analysis to
detect the faults of three phase squirrel cage induction
motor. The analysis is carried out with a CAD package
called Magnet.
FINITE ELEMENT METHOD
Finite element analysis (FEA) is a computer
based numerical technique for calculating the parameters
of electromagnetic devices. It can be used to calculate the
flux density, flux linkages, inductance, torque; induced
emf etc., in the finite element method, the large
electromagnetic device is broken down into many small
elements. The behavior of an individual element can be
described with a relatively simple set of equations. The
equations describing the behavior of the individual
elements are joined into an extremely large set of
equations that describe the whole device. The computer
can solve this large set of simultaneous equations. From
the solution, the computer extracts the behavior of the
individual elements. Finite element methods (FEM) of
analysis have emerged in the past decades as the useful
numerical methods for magnetic field analysis of electrical
machines.
22KW
415V
70.87A
Rated frequency
50Hz
Rated speed
1458 rpm
Number of poles
36
28
Slip
0.028
Efficiency
0.85
Power factor
0.88
(1)
1171
ISSN 1819-6608
www.arpnjournals.com
global energy into its element components and to
minimize one triangle at a time. Then appropriate solution
technique is used to solve the equations and obtain the
necessary parameters like energy, flux function, current,
torque etc.
STATIC ANALYSIS
In this section, the simulation results for the
Static Analysis of Three Phase Induction Motor for
healthy motor, faulty motor with 2, 4, 6 and 8 broken bars
under no load, half load and full load are presented.
(d)
(a)
(a)
(b)
(b)
(c)
(c)
(e)
1172
ISSN 1819-6608
www.arpnjournals.com
5. The flux density and flux function values under various
loads are tabulated in Table-3 and Table-4.
(d)
(a)
(e)
Figure-3. Magnetic field and flux distribution under full
load (a) Healthy (b) 2 broken bar (c) 4 broken bar
(d) 6 broken bar (e) 8 broken bar.
(b)
Healthy
Stored
magnetic
energy (Joules)
0.0005013
Percentage
Change
(%)
-
2 broken bar
0.0004984
0.58
4 broken bar
0.0004939
1.48
6 broken bar
0.0004871
2.83
8 broken bar
0.0004783
4.58
Healthy
0.0273506
2 broken bar
0.0270617
1.05
4 broken bar
0.0266133
2.69
6 broken bar
0.0259341
5.17
8 broken bar
0.0250502
8.41
Healthy
0.0501441
2 broken bar
0.0500637
0.16
4 broken bar
0.0497823
0.72
6 broken bar
0.0490010
2.28
8 broken bar
0.0484359
3.40
Condition
No
load
Half
load
Full
load
(c)
(d)
(e)
Figure-4. Field and flux distribution under no load
(a) Healthy (b) 2 broken bar (c) 4 broken bar
(d) 6 broken bar (e) 8 broken bar.
1173
ISSN 1819-6608
www.arpnjournals.com
and under faulty condition with two-broken bars is 0.0025
Wb, with four-broken bars is 0.0032 Wb, with six-broken
bars is 0.0034 Wb and with eight-broken bars is 0.0035
Wb.
Table-3. Flux function.
Condition
(a)
Healthy
2 broken bar
4 broken bar
6 broken bar
8 broken bar
Healthy
2 broken bar
4 broken bar
6 broken bar
8 broken bar
Healthy
2 broken bar
4 broken bar
6 broken bar
8 broken bar
No
load
(b)
Half
load
Full
load
(c)
Flux
function
(Wb)
0.0023
0.0025
0.0032
0.0034
0.0035
0.0115
0.0131
0.0141
0.0157
0.0191
0.0229
0.0348
0.0428
0.0432
0.0435
Percentage
change (%)
8.69
39.13
47.82
52.17
13.91
22.60
36.52
66.08
51.96
86.89
88.64
89.95
(d)
No
load
Half
load
(e)
Figure-5. Field and flux distribution under half load (a)
Healthy (b) 2 broken bar (c) 4 broken bar (d) 6 broken bar
(e) 8broken bar.
It is observed that the amplitude of the flux
function under no load for healthy condition is 0.023 Wb
Full
load
Healthy
2 broken bar
4 broken bar
6 broken bar
8 broken bar
Healthy
2 broken bar
4 broken bar
6 broken bar
8 broken bar
Healthy
2 broken bar
4 broken bar
6 broken bar
8 broken bar
Flux density
(Wb/m)
0.1229
0.1250
0.1298
0.1363
0.1401
0.8532
0.9061
0.9548
1.0043
1.0425
1.2278
1.3745
1.4830
1.4861
1.5032
Percentage
change (%)
1.70
5.61
10.90
13.99
6.20
11.90
12.89
22.18
11.94
21.03
20.78
22.43
1174
ISSN 1819-6608
www.arpnjournals.com
From the above analysis, it is clear that the
amplitude of flux function and flux density increases as
the number of broken bars increases. There is a drastic
increase in flux function and flux density when the load
increases from no load to full load. Also when the number
of broken bars increased from two broken bars to eight
broken bars, the flux function and flux density increases
and respective percentage change is obtained.
TRANSIENT ANALYSIS
In this section, the simulation results for the
Transient Analysis of Three Phase Induction Motor for
healthy motor, faulty motor with 2, 4, 6 and 8 broken bars
under no load, and full load are presented.
Transient model of induction motor
The electrical model representation of three phase
squirrel cage star connected induction motor is shown in
Figure-6. The power source is considered as a voltage
source connected with the series resistances and
inductances of the stator windings in each phase.
(2)
V2 = Vmsin (wt-2/3)
(3)
V3 = Vmsin (wt+2/3)
(4)
(a)
(b)
(c)
(d)
1175
ISSN 1819-6608
www.arpnjournals.com
(e)
Figure-8. Instantaneous magnetic energy under no load
(a) Healthy (b) 2 broken bar (c) 4 broken bar
(d) 6 broken bar (e) 8broken bar.
Figure-10. Instantaneous magnetic energy- graphical
representation.
Table-5. Instantaneous magnetic energy.
Instantaneous
magnetic
energy (Joules)
Percentage
change (%)
Healthy
1.58446
2 broken bar
1.57958
0.307
4 broken bar
1.28950
18.61
6 broken bar
1.21112
23.56
8 broken bar
1.16965
26.17
Healthy
194.2153
2 broken bar
157.6726
18.81
4 broken bar
125.8326
35.20
6 broken bar
107.2033
44.80
8 broken bar
91.4319
52.92
Condition
(a)
(b)
No load
(c)
(d)
Full
load
(a)
(b)
1176
ISSN 1819-6608
www.arpnjournals.com
the current obtained for healthy is 57.23 A, for 2 broken
bars is 69.97 A, for 4 broken bars is 84.30 A, for 6 broken
bars is 98.35 A and for the 8 broken bars is 106.90 A.
Table-6. Stator phase current.
Condition
(c)
(d)
No
load
(e)
Figure-11. Stator phase current under no load (a) Healthy
(b) 2 broken bar (c) 4 broken bar (d) 6 broken bar
(e) 8broken bar.
Full
load
Healthy
2 broken bar
4 broken bar
6 broken bar
8 broken bar
Healthy
2 broken bar
4 broken bar
6 broken bar
8 broken bar
Stator phase
current (A)
10.40
14.13
14.87
15.07
15.79
57.23
69.97
84.30
98.34
106.9
Percentage
change (%)
35.86
42.96
44.89
51.79
22.25
47.29
71.83
86.78
(a)
(b)
(C)
(d)
(a)
(b)
(c )
(d)
(e)
Figure-2. Stator phase current under full load (a) Healthy
(b) 2 broken bar (c) 4 broken bar (d) 6 broken ba
(e) 8 broken bar.
The stator phase current values under various
loads are tabulated in Table-6. Under no load, the current
obtained for healthy is 10.40 A, for 2 broken bars is 14.13
A, for 4 broken bars is 14.87 A, for 6 broken bars is 15.07
A and for the 8 broken bars is 15.79 A. Under full load,
(e)
Figure-13. Flux linkage under no load (a) Healthy (b) 2
broken bar (c) 4 broken bar (d) 6 broken bar
(e) 8 broken bar.
1177
ISSN 1819-6608
www.arpnjournals.com
Hence it is observed that the value of flux linkage
increases when the number of broken bars increases.
(a)
(b)
(c)
(d)
Healthy
1.4562
Percentage
change
(%)
-
2 broken bar
2.0683
42.03
4 broken bar
2.5302
73.75
6 broken bar
2.7809
90.96
8 broken bar
2.8211
93.73
Healthy
123.684
2 broken bar
161.468
30.54
4 broken bar
187.245
51.38
6 broken bar
207.047
67.39
8 broken bar
225.006
81.92
Magnetic
torque (Nm)
Condition
No load
(d)
Figure-14. Flux linkage under full load (a) Healthy
(b) 2 broken bar (c) 4 broken bar (d) 6 broken bar
(e) 8 broken bar
It can be observed that under no load, the flux
linkage obtained for healthy is 0.01104 Wb, for 2 broken
bars is 0.01327 Wb, for 4 broken bars is 0.01791 Wb, for
6 broken bars is 0.01822 and for 8 broken bars is 0.02205
Wb. Under full load, the flux linkage obtained for healthy
is 0.16668 Wb, for 2 broken bars is 0.19386 Wb, for 4
broken bars is 0.22034 Wb, for 6 broken bars is 0.25851
Wb and for 8 broken bars is 0.31012 Wb.
Full
load
Healthy
0.01104
Percentage
change
(%)
-
2 broken bar
0.01327
20.19
4 broken bar
0.01791
62.22
6 broken bar
0.01822
65.04
8 broken bar
0.02205
99.72
Flux Linkage
(Wb)
Condition
No
load
Full
load
Healthy
0.16668
2 broken bar
0.19386
16.30
4 broken bar
0.22034
32.19
6 broken bar
0.25851
55.09
8 broken bar
0.31012
86.05
1178
ISSN 1819-6608
www.arpnjournals.com
CONCLUSIONS
A three phase squirrel cage induction motor is
modeled on the basis of finite element method. The
simulation result was obtained for broken bar faults.
Comparisons are made with the healthy motor condition
and the result was tabulated. It was found that the faults
due to the broken bars saturate the magnetic field
distribution on the rotor tooth adjacent to the bars that
were broken. In the static analysis, the stored magnetic
energy decreased when the number of broken bars
increased. The Flux Function and Flux Density was
increased when the number of broken bars was increased.
In the transient analysis, it was found that the flow of
current in the stator phases and the flux linkage produced
in the motor was increased, whereas, the stored magnetic
energy was decreased when the number of broken bars
was increased. Also the magnetic torque was increased
when the number of broken bars was increased. The
simulated results are verified theoretically for the values
like flux function and flux density.
REFERENCES
[1] Balamurugan S., Arumugam R., Paramasivam S. and
Malaiappan M. 2004. Transient Analysis of induction
Motor Using Finite Element Analysis. IEEE Industrial
Electronics Society, 30th annual conference. pp. 15261529.
[2] John.F. Bangura, Richard J. Povinelli, Nabeel A. O.
Demerdash and Ronald H. Brown. 2003. Diagnostics
of Eccentricities and Bar/End-Ring Connector
Breakages in Polyphase Induction Motors through a
combination of Time-Series Data Mining and Time
Stepping Coupled FE - State Space Technique. IEEE
Trans. on Industrial Applications. 39(4): 1005-1013.
[3] Bentounsi A. and Nicolas A. 1998. On Line Diagnosis
of Defaults on Squirrel Cage Motor Using FEM. IEEE
Trans. Magn. 34(5): 3511-3574.
[4] Bianchi N., Bolognani S. and Comelato G. 1999.
Finite Element Analysis of Three Phase Induction
Motors: Comparison of Two Different Approaches.
IEEE Trans. on Energy Conversion. 14(4): 15231528.
[5] Elkasabgy N. M. and Eastham A. R. 1992. Detection
of Broken Bars in the Cage Rotor on an Induction
Machine. IEEE Trans. Ind. Appl. 28(1): 165-171.
[6] Genadi Y. Sizov and Ahmed Sayed-Ahmed nabeel.
2009. Analysis and Diagnostics of Adjacent and
Nonadjacent Broken-Rotor-Bar Faults in SquirrelCage Induction Machines. IEEE Trans. on Industrial
Electronics. 56(11): 4627-4641.
1179