Massachusetts Institute of Technology

Department of Electrical Engineering and Computer Science

6.685 Electric Machinery
Class Notes 10: Induction Machine Control and Simulation
c 2003 James L. Kirtley Jr.

Introduction

The inherent attributes of induction machines make them very attractive for drive applications.
They are rugged, economical to build and have no sliding contacts to wear. The difficulty with
using induction machines in servomechanisms and variable speed drives is that they are hard to
control, since their torque-speed relationship is complex and nonlinear. With, however, modern
power electronics to serve as frequency changers and digital electronics to do the required arithmetic,
induction machines are seeing increasing use in drive applications.
In this chapter we develop models for control of induction motors. The derivation is quite brief
for it relies on what we have already done for synchronous machines. In this chapter, however, we
will stay in ordinary variables, skipping the per-unit normalization.

Volts/Hz Control

Remembering that induction machines generally tend to operate at relatively low per unit slip, we
might conclude that one way of building an adjustable speed drive would be to supply an induction
motor with adjustable stator frequency. And this is, indeed, possible. One thing to remember is
that flux is inversely proportional to frequency, so that to maintain constant flux one must make
stator voltage proportional to frequency (hence the name constant volts/Hz). However, voltage
supplies are always limited, so that at some frequency it is necessary to switch to constant voltage
control. The analogy to DC machines is fairly direct here: below some base speed, the machine
is controlled in constant flux (volts/Hz) mode, while above the base speed, flux is inversely
proportional to speed. It is easy to see that the maximum torque varies inversely to the square of
flux, or therefore to the square of frequency.
To get a first-order picture of how an induction machine works at adjustable speed, start with
the simplified equivalent network that describes the machine, as shown in Figure 1
In Chapter 8 of these notes it is shown that torque can be calculated by finding the power
dissipated in the virtual resistance R 2 /s and dividing by electrical speed. For a three phase machine,
and assuming we are dealing with RMS magnitudes:
p
R2
Te = 3 |I2 |2

s
where is the electrical frequency and p is the number of pole pairs. It is straightforward to find
I2 using network techniques. As an example, Figure 2 shows a series of torque/speed curves for
an induction machine operated with a wide range of input frequencies, both below and above its
base frequency. The parameters of this machine are:

X1

Ia Ra

X2 I2

<

Xm <> R2

>
< s

Figure 1: Equivalent Circuit

Number of Phases
Number of Pole Pairs
RMS Terminal Voltage (line-line)
Frequency (Hz)
Stator Resistance R1
Rotor Resistance R2
Stator Leakage X1
Rotor Leakage X2
Magnetizing Reactance Xm

3
3
230
60
.06
.055
.34
.33
10.6

A strategy for operating the machine is to make terminal voltage magnitude proportional to frequency for input frequencies less than the Base Frequency, in this case 60 Hz, and to hold voltage
constant for frequencies above the Base Frequency.
Induction Motor Torque

250

Nm

200

150

100

50

0
0

50

100

150
Speed, RPM

200

250

Figure 2: Induction Machine Torque-Speed Curves

For high frequencies the torque production falls fairly rapidly with frequency (as it turns out,
2

it is roughly proportional to the inverse of the square of frequency). It also falls with very low
frequency because of the effects of terminal resistance. We will look at this next.

2.1

Idealized Model: No Stator Resistance

Ia X1

X2 I2


<

+
Xm <> R2

V
<> s


Figure 3: Idealized Circuit: Ignore Armature Resistance

Ignore, for the moment, R1 . An equivalent circuit is shown in Figure 3. It is fairly easy to show
that, from the rotor, the combination of source, armature leakage and magnetizing branch can be
replaced by its equivalent circuit, as shown in in Figure 4.

X10
Ia


+

X2 I2

<>
<> R2
< s

V0



Figure 4: Idealized Equivalent

In the circuit of Figure 4, the parameters are:
Xm
Xm + X 1
= Xm ||X1

V0 = V
X0

If the machine is operated at variable frequency , but the reactance is established at frequency
B , current is:
I=

V0
j(X10 + X2 ) B +

R2
s

and then torque is

|V 0 |2 Rs2
3p
Te = 3|I2 |
=
s
(X10 + X2 )2 + ( Rs2 )2
2 R2

Now, if we note that what counts is the absolute slip of the rotor, we might define a slip with
respect to base frequency:
r
r B
B
s=
=
= sB

Then, if we assume that voltage is applied proportional to frequency:

V 0 = V00

and with a little manipulation, we get:

Te =

2
|V00 |2 R
3p
sB
2 2
B (X10 + X2 )2 + ( R
sB )

This would imply that torque is, if voltage is proportional to frequency, meaning constant applied
flux, dependent only on absolute slip. The torque-speed curve is a constant, dependent only on the
difference between synchronous and actual rotor speed.
This is fine, but eventually, the notion of volts per Hz runs out because at some number of
Hz, there are no more volts to be had. This is generally taken to be the base speed for the drive.
Above that speed, voltage is held constant, and torque is given by:
2
|V 0 |2 R
3p
sB
Te =
2 2
B (X10 + X2 )2 + ( R
s )
B

The peak of this torque has a square-inverse dependence on frequency, as can be seen from Figure 5.
Induction Motor Torque

250

Nm

200

150

100

50

0
0

500

1000
1500
Speed, RPM

2000

Figure 5: Idealized Torque-Speed Curves: Zero Stator Resistance

2.2

Peak Torque Capability

Assuming we have a smart controller, we are interested in the actual capability of the machine. At
some voltage and frequency, torque is given by:
Te = 3|I2 |2

3 p |V 0 |2 Rs2
R2
=
s
((X10 + X2 )( B ))2 + (R10 +

R2 2
s )

Now, we are interested in finding the peak value of that, which is given by the value of R 2 s which
maximizes power transfer to the virtual resistance. This is given by the matching condition:
R2
=
s

R102 + ((X10 + X2 )(

2
))
B

Then maximum (breakdown) torque is given by:

Tmax =

3p
0 2
|V |

R102 + ((X10 + X2 )( B ))2

((X10 + X2 )( B ))2 + (R10 +

R102 + ((X10 + X2 )( B ))2 )2

This is plotted in Figure 6. Just as a check, this was calculated assuming R 1 = 0, and the
results are plotted in figure 7. This plot shows, as one would expect, a constant torque limit region
to zero speed.
Breakdown Torque
300

250

NewtonMeters

200

150

100

50

0
0

20

40

60
80
Drive Frequency, Hz

100

120

Figure 6: Torque-Capability Curve For An Induction Motor

Field Oriented Control

One of the more useful impacts of modern power electronics and control technology has enabled
us to turn induction machines into high performance servomotors. In this note we will develop a
5

Breakdown Torque
300

NewtonMeters

250

200

150

100

50
0

20

40

60
80
Drive Frequency, Hz

100

120

Figure 7: Idealized Torque Capability Curve: Zero Stator Resistance

picture of how this is done. Quite obviously there are many details which we will not touch here.
The objective is to emulate the performance of a DC machine, in which (as you will recall), torque
is a simple function of applied current. For a machine with one field winding, this is simply:
T = GIf Ia
This makes control of such a machine quite easy, for once the desired torque is known it is easy to
translate that torque command into a current and the motor does the rest.
Of course DC (commutator) machines are, at least in large sizes, expensive, not particularly
efficient, have relatively high maintenance requirements because of the sliding brush/commutator
interface, provide environmental problems because of sparking and carbon dust and are environmentally sensitive. The induction motor is simpler and more rugged. Until fairly recently the
induction motor has not been widely used in servo applications because it was thought to be hard
to control. As we will show, it does take a little effort and even some computation to do the
controls right, but this is becoming increasingly affordable.

3.1

Elementary Model:

We return to the elementary model of the induction motor. In ordinary variables, referred to the
stator, the machine is described by flux-current relationships (in the d-q reference frame):
"

dS
dR

"

qS
qR

"

LS
M

M
LR

#"

idS
idR

"

LS
M

M
LR

#"

iqS
iqR

Note the machine is symmetric (there is no saliency), and since we are referred to the stator,
the stator and rotor self-inductances include leakage terms:
LS

= M + LS`

LR = M + LR`
The voltage equations are:
ddS
qS + rS idS
dt
dqS
vqS =
+ dS + rS iqS
dt
ddR
0 =
s qR + rR idR
dt
dqR
0 =
+ s dR + rR iqR
dt
Note that both rotor and stator have speed voltage terms since they are both rotating with
respect to the rotating coordinate system. The speed of the rotating coordinate system is w with
respect to the stator. With respect to the rotor that speed is s = m , where m is the rotor
mechanical speed. Note that this analysis does not require that the reference frame coordinate
system speed w be constant.
Torque is given by:
3
T e = p (dS iqS qS idS )
2
vdS

3.2

Simulation Model

As a first step in developing a simulation model, see that the inversion of the flux-current relationship is (we use the d- axis since the q- axis is identical):
idS

idR =

LR
M
dS
dR
LS LR M 2
LS LR M 2
M
LS
dS
dR
LS LR M 2
LS LR M 2

Now, if we make the following definitions (the motivation for this should by now be obvious):
Xd = 0 LS
Xkd = 0 LR
Xad = 0 M
Xd0

M2
LS
LR

= 0

the currents become:

idS

idR =

0
Xad 0
dS
dR
Xd0
Xkd Xd0
Xad 0
Xd 0
dS 0
dR
0
Xk d Xd
Xd Xkd
7

The q- axis is the same.

Torque may be, with these calculations for current, written as:
3
3 0 Xad
Te = p (dS iqS qS idS ) = p
(dS qR qS dR )
2
2 Xkd Xd0
Note that the usual problems with ordinary variables hold here: the foregoing expression was
written assuming the variables are expressed as peak quantities. If RMS is used we must replace
3/2 by 3!
With these, the simulation model is quite straightforward. The state equations are:
ddS
dt
dqS
dt
ddR
dt
dqR
dt
dm
dt

= VdS + qS RS idS
= VqS dS RS iqS
= s qR RR idR
= s dR RS iqR
=

1
(Te + Tm )
J

where the rotor frequency (slip frequency) is:

s = pm
For simple simulations and constant excitaion frequency, the choice of coordinate systems is
arbitrary, so we can choose something convenient. For example, we might choose to fix the coordinate system to a synchronously rotating frame, so that stator frequency = 0 . In this case, we
could pick the stator voltage to lie on one axis or another. A common choice is V d = 0 and Vq = V .

3.3

Control Model

If we are going to turn the machine into a servomotor, we will want to be a bit more sophisticated
about our coordinate system. In general, the principle of field-oriented control is much like emulating the function of a DC (commutator) machine. We figure out where the flux is, then inject
current to interact most directly with the flux.
As a first step, note that because the two stator flux linkages are the sum of air-gap and leakage
flux,
dS

= agd + LS` idS

qS

= agq + LS` iqS

This means that we can re-write torque as:

3
T e = p (agd iqS agq idS )
2
8

Next, note that the rotor flux is, similarly, related to air-gap flux:
agd = dR LR` idR
agq = qR LR` iqR

Torque now becomes:

3
3
T e = p (dR iqS qR idS ) pLR` (idR iqS iqR idS )
2
2
Now, since the rotor currents could be written as:
dR
M

idS
LR
LR
qR
M

iqS
LR
LR

idR =
iqR =
That second term can be written as:
idR iqS iqR idS =

1
(dR iqS qR idS )
LR

So that torque is now:

Te =

3.4

3
LR`
p 1
2
LR


3 M
(dR iqS qR idS ) = p
(dR iqS qR idS )
2 LR

Field-Oriented Strategy:

What is done in field-oriented control is to establish a rotor flux in a known position (usually this
position is the d- axis of the transformation) and then put a current on the orthogonal axis (where
it will be most effective in producing torque). That is, we will attempt to set
dR = 0
qR = 0
Then torque is produced by applying quadrature-axis current:
Te =

3 M
p
0 iqS
2 LR

The process is almost that simple. There are a few details involved in figuring out where the
quadrature axis is and how hard to drive the direct axis (magnetizing) current.
Now, suppose we can succeed in putting flux on the right axis, so that qR = 0, then the two
rotor voltage equations are:
0 =
0 =

ddR
s qR + rR IdR
dt
dqR
+ s dR + rR IqR
dt
9

Now, since the rotor currents are:

idR =
iqR =

dR
M

idS
LR
LR
qR
M

iqS
LR
LR

The voltage expressions become, accounting for the fact that there is no rotor quadrature axis
flux:


ddR
dR
M
0 =
+ rR

idS
dt
LR
LR
M
0 = s dR rR
iqS
LR
Noting that the rotor time constant is
TR =

LR
rR

we find:
TR

ddR
+ dR = M idS
dt
M iqS
s =
TR dR

The first of these two expressions describes the behavior of the direct-axis flux: as one would
think, it has a simple first-order relationship with direct-axis stator current. The second expression,
which describes slip as a function of quadrature axis current and direct axis flux, actually describes
how fast to turn the rotating coordinate system to hold flux on the direct axis.
Now, a real machine application involves phase currents i a , ib and ic , and these must be derived
from the model currents idS and iqs . This is done with, of course, a mathematical operation which
uses a transformation angle . And that angle is derived from the rotor mechanical speed and
computed slip:
Z
=

(pm + s ) dt

A generally good strategy to make this sort of system work is to measure the three phase currents
and derive the direct- and quadrature-axis currents from them. A good estimate of direct-axis flux is
made by running direct-axis flux through a first-order filter. The tricky operation involves dividing
quadrature axis current by direct axis flux to get slip, but this is now easily done numerically (as
are the trigonometric operations required for the rotating coordinate system transformation). An
elmentary block diagram of a (possbly) plausible scheme for this is shown in Figure 8.
In this picture we start with commanded values of direct- and quadrature- axis currents, corresponding to flux and torque, respectively. These are translated by a rotating coordinate transformation into commanded phase currents. That transformation (simply the inverse Parks transform)
uses the angle q derived as part of the scheme. In some (cheap) implementations of this scheme
the commanded currents are used rather than the measured currents to establish the flux and slip.
We have shown the commanded currents i a , etc. as inputs to an Amplifier. This might be
implemented as a PWM current-source, for example, and a tight loop here results in a rather high
performance servo system.
10

dR

N
D

M
1 + STa

M
TR

i a*

*
d

ib*

-1

T
i q*

ia
o

Amp

o
o

*
c

Figure 8: Field Oriented Controller

11

ic

Motor
Load

ib

MIT OpenCourseWare
http://ocw.mit.edu

6.685 Electric Machines

Fall 2013

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.