Advanced Electromechanical Systems Course Manual PDF
Advanced Electromechanical Systems Course Manual PDF
Advanced Electromechanical Systems Course Manual PDF
ELE847
Advanced Electromechanical Systems
Course Notes
2009 Edition
ELE847 Advanced Electromechanical Systems
Table of Contents
Course Description
A course on modeling and simulation of electromechanical systems. The main topics include: modeling of dc
motors, dc motor dynamic performance, reference frame theory, modeling of induction and synchronous machines,
small signal (linearized) analysis, solid state converters, advanced motor speed control schemes, and simulation
techniques. The modeling and simulation techniques developed in this course provide a useful tool for the analysis
and design of electric machines, power electronics circuits and dc/ac motor drives.
Prerequisite
All required third year courses.
Course Organization
This course consists of three hours of lecture and one hour of laboratory per week.
Course Material
Text: "Analysis of Electric Machinery and Drive Systems" by P.C. Krause, O.Wasynczuk and S.D.Sudhoff,
published by Wiley-IEEE Press, 2002. ISBN 0-471-14326-X
Instructor
Bin Wu, Ph.D., P.Eng., Professor
Room ENG328, 245 Church Street, Toronto
Department of Electrical and Computer Engineering
Ryerson University
(416) 979-5000 ext: 6484
In order to achieve a passing grade, the student must achieve an average of at least 50% in both theoretical and
laboratory components.
Course Material
1. Course Outline
4. Lecture Slides
A1. Objectives
- Build a Simulink model for a separately exited DC motor; and
- Study dynamic performance of the motor.
File: Dcm2.mdl
Ia ia
Va Load
Armature Current [A] Scope To Workspace
Step Voltage [V]
Te Te
Motor
Torque [N.m] Scope1 To Workspace1
0 TL
Load
Const Torque [N.m] n n
Wm -K-
Motor To Workspace2
Gain Speed [rpm] Scope2
Vf
if
Field
Field
Step1 Voltage [V]
Current Scope3
DC Motor2
(Masked subsystem) t
Clock To Workspace3
2) Run the Simulink model and plot the waveforms of va, vd1, vd2, vdc and idc assuming that the phase voltage
of the three phase ac supply is 110V (60Hz).
3) Calculate the average value of the rectifier output voltage (vdc) based on the simulated waveforms.
2) Assume that the phase voltage of the ac supply is 110V (60Hz). Run the Simulink model and plot the
waveforms of of va, vd1, vd2, vdc and idc with the delay angle of 30 and 90 degrees respectively.
3) Derive an expression which can be used to calculate the average value of the rectifier output voltage (vdc).
Verify this expression using simulated waveforms.
Dcdrv1.mdl
Va Vdc
Armature Current
Vb Ia
3 phase Va
Vdc Ia
Vc
Electromagnetic Torque
Te
AC Supply (N.m)
Te
SCR Rectifier 29.2 TL
t_d
Speed (rad/s)
Load (N.m) Wm -K-
Firing rpm
acos fcn included n
Circuit Gain
Dcdrv2.mdl
Va
Vb Ia
3 phase Mux
Vdc Va
Vc Ia
Mux
Te Torque
AC Supply (N.m)
Te
SCR Rectifier 29.2 TL
t_d
Load (N.m) Wm -K-
Firing rpm
acos fcn included n
Circuit Gain
Sum2
Current
32.4
Ref (A)
Va
Vb Ia
3 phase Mux
Vdc Va
Vc Ia
Mux
Te
AC Supply
Te
SCR Rectifier 10 TL
t_d rpm
Load (N.m) Wm -K-
Firing
acos fcn included n
Circuit Gain
240 Vf
V_alpha if
Vf (V)
Current If
PI
PI DC Motor2
(with limiters)
(Separately Excited)
Sum1
Ia
Ia (Ref)
Sum2
Wm (rad/s)
Speed Ref
63.9
(rad/s)
2) Run the model and plot the waveforms of armature current and rotor speed;
3) Based on simulation results, calculate the speed overshoot (%).
4) Determine the speed PI regulator parameters such that the speed overshoot is approximately 5% and the
speed settling time as short as possible.
C2.3 Print the Simulink model including all subsystems.
Post-lab Report
Refer to Lab 1 for general instruction on post-lab report.
Objectives
- To build a Simulink model for three phase induction motors; and
- To investigate induction motor dynamic performance.
AC Supply
A three-phase power supply with phase voltage of 127V (rms) and frequency of 60Hz.
File: Lab3.mdl
Mux Ia
Va
Vqs iqs
Vb 3-phase 2-phase
Ia
3 phase To To
2-phase 3-phase Te
Vc Vds ids
To Workspace1
AC Supply 3-2 Transform
2-3 Transform
0 Tl Te
Load Torque
Te
0 W Wrm 60/(2*pi)
Stator Frame
n
Gain
IM_dq_Arbi
n
t
Induction motor model To Workspace2
Clock To Workspace3 in the arbitrary frame
Build the induction motor d-q model using the equivalent circuit given in Fig. 4.5-1 (Page 151). The zero-axis
equivalent circuit can be neglected since this is a 3-phase balance system. The torque-speed relationship is
described by Eq. 4.3-8 and the electromagnetic torque generated by the motor can be calculated according to
Eq. 4.6-4.
t st
Te1
Te,start
t
Te 2
Post-lab Report
Refer to Lab 1 for general instruction on post-lab report.
Objectives
- To build Simulink models for three phase voltage source inverters; and
- To investigate PWM inverter performance.
Lab4a.mdl
Van
G1
Van Ia
G3
Gating
Vbn Ib
VSI RL Load Mux
G5
Load
Current
Square Wave Vcn Ic
Gating Generator Mux1
Vdc
A note on simulation parameters. You can either use fixed-step or variable-step differential equation solver. If you
use a fixed-step differential equation solver with a large time step, you may not be able to obtain accurate results or
the results may even be wrong. This is mainly due to the switching operation of the inverter and small time
constants that the drive system may contain. If you choose Runge Kutta (ode4) method, you may try to use a step
size of 10 s or smaller.
A2.2 Plot the harmonic spectrum of the waveforms of van and ia. Frequency range for the plot: 0 to 2kHz.
Note: This is a common task for electrical engineers working in the area of power electronics and motor drives.
Van
20
G1
Freq Sine Wave G1 Van Ia
Sine
0.8
G3
G3
Md
Vbn Ib
3-phase VSI RL Load Mux
SW Generator Carrier G5
G5 Load
Current
PWM Vcn Ic
Carrier Wave Mux
Generator
250
Vdc
Other Blocks
The specifications for the other blocks are given in Part A.
B2.2 Plot the harmonic spectrum of the waveforms of van and ia in Part B2-4). Frequency range for the plot: dc to
2kHz. Compare the harmonic spectrum with that in Part A2.2, and then make conclusions.
Post-lab Report
Refer to Lab 1 for general instruction on post-lab report.
Objectives
To investigate characteristics of two induction motor speed control systems.
Part A Induction Motor Speed Control Using a Six-step Voltage Source Inverter
A.1 Model Building
Build the Simulink model according to the block diagram shown in Fig. 1. You may use some of the models you
built in the previous lab sessions. The parameters of the induction motor speed control system are as follows.
Induction motor Nameplate data: 3, 3hp, 220V, 8.4A (rated) and 1710rpm. Use the motor parameters given
in Table 4.10-1, P165, textbook.
DC link resistor This resistor represents the power loss of the SCR rectifier and dc link bus. The equivalent
resistance is 0.5.
Low pass filter This is a second order lower pass filter with dc gain k = 1 and quality factor Q = 1. The
corner frequency of the filter should be the same as the reference frequency of the drive. The
LP filter is used to extract the fundamental component from the six-step inverter output
voltage Van.
Power Supply 60Hz, 127V per phase.
Load torque 12.5N.m (rated torque)
Other constants Kv = 0.95/60 and Vcomp = 0.
Proja.mdl Mux
Vdc
Vdc Van
Vqs iqs
Va 3-phase 2-phase Ia
G1 To To
Vb 2-phase 3-phase
3 phase DC Link Vds ids
Vc VSI
Filter G3 Vbn
Vdo
2-3 Transform
AC Supply 12.5 Tl Te
0.05
LP Filter
Vcomp
Kv
-K- (2nd Order)
Van1
60 Ref [Hz]
2) Start the drive system until a steady state operation is reached. Complete the following table.
V / f Ratio (Volts,rms/Hz)
3) To compensate the voltage drop on stator winding resistance at low frequencies, let Kv = 0.9/60 and
Vcomp = 0.05. Run the drive system and complete the following table.
V / f Ratio (Volts,rms/Hz)
4) Start the drive system at 60Hz under the operating conditions given in 2) until a steady state operation is
reached. Plot the transient waveforms of motor speed n, stator current Ia, motor torque Te, diode rectifier
output voltage Vdo and inverter input voltage Vdc.
Report
The formal report should include the following items:
1) Cover page (including project title, your name, student ID, date)
2) Abstract (a paragraph of about 200 words)
3) Theory (two full pages, 1.5 line space)
4) All required waveforms, tables and calculations
5) Comment on the size of the dc link capacitors used in both systems
6) Compare the simulation results obtained in Part A.2-2 and A.2-3 by answering the following questions:
- Is the V/f ratio constant? You may draw V versus f curves for comparison.
- Is the slip speed constant? Why?
7) Compare the simulation results obtained in Part A.2-2 and B.2
8) Comments on harmonic issues of the two systems
9) Conclusions (300 400 words)
10) Appendix: Simulink model in Part B including block diagrams of all subsystems.
1.1 Consider the dynamic equivalent circuit of a shunt dc motor given in Fig. 2.4-2 (page 79, textbook). Drive
a block diagram for this motor. It is assumed that the armature voltage and load torque are input variables
while the armature current, rotor speed and electromagnetic torque are output variables. Show these
variables on the diagram.
1.2 Repeat Problem 1.1 for a series dc motor using the equivalent circuit given in Fig. 2.4-6 (page 83,
textbook). Answer: Discussed in the lecture class.
1.3 Formulate the following transfer functions for a shunt dc motor under the assumption that the field
current If is constant:
I (s)
(a) a assuming that armature voltage is zero. This transfer function can be used to study
T L (s)
the dynamic response of armature current due to changes in load torque.
Kv
I (s) J a ra
Answer: a =
S + (1/ a + B m /J )S + (B m /J + 1/ m )
T L (s) 2 1
a
J ra
where m = 2
Kv
I a (s)
(b) assuming that the load torque is zero. This transfer function can be used to study the dynamic
V a (s)
response of armature current due to changes in armature voltage.
J S + Bm
Answer: a =
I (s) ra J a
S + (1 / a + B m /J )S + (B m /J + 1/ m )
V a (s) 2 1
a
T e (s)
(c) assuming that the load torque is zero.
V a (s)
J S + Bm
Kv
T (s) ra J a
Answer: e =
S + (1 / a + B m /J )S + (B m /J + 1/ m )
V a (s) 2 1
a
1.4 Derive an expression for motor speed r(s) in terms of armature voltage Va(s) and load torque TL(s) for a
separately excited motor with a constant field current. Refer to Page 97 of textbook for answers.
1.5 Sketch to scale the dc side waveforms ( v d 1 , v d 2 and v dc = v d 1 v d 2 ) of a single phase full-wave
thyristor rectifier with a delay angle of 30 and 90 degrees respectively.
1.6 Repeat Problem 1.5 for a three phase full-wave thyristor rectifier.
1.7 Using standard Simulink blocks, derive a Simulink model for a single phase full-wave thyristor rectifier.
1.9 Derive an expression which can be used to calculate the average dc output voltage of a three phase full-
wave thyristor rectifier.
2.2 Assume that i and i are variables in the stationary reference frame and iq and id are variables in the arbitrary
frame which rotates in space at an arbitrary speed of as shown in Fig. 2.1 below. Verify the following
equations which can be used to transform the variables in the stationary frame to the arbitrary frame.
iq = i cos i sin ; and id = i sin + i cos (Two-phase to two-phase transformation).
q-axis
axis
(Stationary Frame)
d-axis
axis
(Stationary Frame)
2.4 Derive a coefficient matrix which can be used to transform abc variables in the stationary reference frame to
qdo variables also in the stationary reference frame, assuming that the q-axis is coincident with the a-axis.
2.5 Derive an equation which can be used to transform abc variables in the stationary reference frame to qdo
variables in the arbitrary reference frame which rotates in space at a speed of .
2.6 Derive an equation which can be used to transform qdo variables in the arbitrary reference frame to abc
variables in the stationary reference frame.
2.7 Derive arbitrary-frame (q-d) equivalent circuits for a three-phase balanced capacitor bank.
(Hint: Refer to Pages 119-120, textbook).
2.8 Derive arbitrary-frame (q-d) equivalent circuits for a three-phase RL circuit. It is assumed that
1) the RL circuit is three-phase balanced;
2) the resistors and inductors are connected in series; and
3) no mutual inductances exist between any two phases.
(Hint: Refer to Pages 120-122, textbook).
3.1 A simplified version of induction motor dq model in the arbitrary reference frame is shown in Fig. 3.1, where
LL is the total leakage inductance of the stator and rotor windings. It is also called the equivalent circuit of
the induction motor.
a) Express dq-axis flux linkages in terms of motor inductances and currents in a matrix form.
b) Express dq-axis currents in terms of motor inductances and flux linkages. Based on the derived equations,
draw a block diagram using standard Simulink blocks.
rs ds LL (-r )dr rr
q-axis
rs qs LL (-r )qr rr
d-axis
Fig. 3.1 Induction motor dq model in the arbitrary reference frame, where the stator
and rotor winding leakage inductances are lumped together (the equivalent circuit).
3.2 Derive equations for the calculation of the dq voltages specified in Fig. 3.1. Assume that this is a wound rotor
induction motor where the rotor winding is open. Following the same procedure discussed in the lecture class,
derive a Simulink model for the induction motor.
rs sds Lls Llr ( s-r )dr rr
q-axis
d-axis
Fig. 3.2 Induction motor dq model in the synchronous reference frame.
3.3 The induction dq model in the synchronous frame is shown in Fig. 3.2, where s is the speed of the
synchronously rotating frame. Assuming that the rotor winding is shorted, repeat the questions given in 3.1
and 3.2.
3.5 Repeat Problem 3.4 except the induction motor model is in the rotor reference frame.
4.1 A three-phase induction motor has the following nameplate data: 10hp, 60Hz, 220V and 1150rpm. The
maximum torque (Tmax) of the motor is 155N.m, which occurs at the motor speed of 970rpm. The starting
torque of the motor is 70N.m. Assuming that the air gap flux of the motor is kept the constant, sketch to scale
the toque versus speed curves at the stator frequency of 60Hz, 40Hz and 20Hz.
4.2 Assuming that the input dc voltage of a three-phase IGBT-based voltage source inverter is 200V and the duty
cycle of the IGBTs is 50% (i.e., the IGBT conduction angle is 180 degrees).
a) Determine the amplitude of the fundamental component, 5th, 7th and 11th harmonics of the inverter output
voltage (line-to-line).
b) Sketch to scale the line-to-neutral voltage waveform (Phase a) and determine the amplitude of the
fundamental component, 5th, 7th and 11th harmonics.
4.3 Draw the inverter output voltage waveform (line-to-line) assuming that the IGBT devices in Question 4.2 have
a conduction angle of 120 degrees per cycle. Assume that the inverter is loaded with a three-phase balanced
resistor. Hint: the IGBT gating signals for the upper and lower IGBTs in the same inverter leg are no longer
complementary.
4.4. Derive a Simulink model for a three phase voltage source inverter, assuming the conduction angle of the
switching devices is 180 degrees.
4.5 Draw a block diagram (not Simulink Model) for a three-phase induction motor drive using Volts per Hertz
control scheme. To ensure a constant flux operation, voltage feedback should be used. The voltage drop on the
stator winding resistance at low operating frequencies should also be compensated. The drive system is
implemented with a sine pulse width modulator.
4.6 Fig. 4.1 shows the block diagram of an induction motor field oriented control scheme. Derive equations and
then Simulink block diagrams for the following subsystems:
Is* & T* Resolver, Current Reference Generator and Flux & Torque Estimator.
f ias
ibs
Te Flux & ics
Torque
vas
r Estimator
vbs
Tachometer
IM
5.1 Use the synchronous machine model discussed in the lecture class or given in Figure 5.5-1 (Page 202,
textbook). Note: Superscripts associated with machine variables may be ignored.
(a) Express flux linkages (qs, ds, kq1, kq2, kd, and fd) in terms of machine currents
(qs, ds, kq1, kq2, kd, and fd).
(b) Derive expressions for the following voltages: Vqs, Vds, Vkq1, Vkq2, Vkd, and Vfd.
5.2 For small size synchronous machines, the damper windings may be omitted to reduce manufacturing cost.
Assuming that for a synchronous machine, the kq2 and kd windings are not equipped, repeat questions given
in 5.1.
5.3 State the reasons why the rotor reference frame is often used for synchronous machine analysis.
5.4 Briefly explain the main functions of damper windings. What are the main differences between the
synchronous machine and induction machine.
D1 D3 D5
va
ia
vb
ib
vdc RL
vc
ic
D4 D6 D2
vd 2
(Free acceleration)
(Free acceleration)
(Free acceleration)
(Free acceleration)
ias ( pu ) nr ( rpm )
t (sec )
T1 D1 T3 D3 T5
D5
G1 G3 G5
ia
a
ib
Vdc b n
c ic
T4 D4 va T6
D6
T2
D2
G4 G6 G2
GND
G1 t
G2
G3 t
G4
G5
G6 t
2 3
/3 /3
va Vdc
t
vb Vdc
t
vc Vdc
t
vab Vdc
t
vbc Vdc
t
vca
Vdc
t
I II III IV V VI
Fig. 4-1 Three phase voltage source inverter with square wave operation.
T1 D1 T3 D3 T5
D5
G1 G3 G5
ia
a
ib
Vdc b n
c ic
T4 D4 va T6
D6
T2
D2
G4 G6 G2
GND
va
(G1 ) Vdc
t
vb
(G3 ) Vdc
t
vab1
vab
Vdc
t
2
vab = va vb
r (rpm)
r*
r
T (N m)
e
ias ( A)
(
r
a
d
)
(Wb),
r f
r
f
Fig. 4-3 Simulated waveforms of field oriented control for induction motor drives.
ikq1
rkq1
Vkq1
Llkq1
iqs rs rds Lls pkq1
imq
Llkq2
vqs pqs Lmq
rkq2
ikq2
pkq2
q-axis
Vkq2
ikd
rkd
Vkd
Llkd
ids rs rqs Lls pkd
imd
Llfd
vds pds Lmd
rfd
ifd
pfd
d-axis
Vfd
Fig. 5-2 dq-axis model of synchronous generator in the rotor reference frame.