Direct Torque Control with Space Vector Modulation (DTC-SVM) of ...

POLITECHNIKA
WARSZAWSKA
WARSAW UNIVERSITY OF TECHNOLOGY
Faculty of Electrical Engineering
ROZPRAWA DOKTORSKA
Ph.D. Thesis
Dariusz Świerczyński, M. Sc.
Direct Torque Control with
Space Vector Modulation (DTC-SVM) of Inverter-Fed
Permanent Magnet Synchronous Motor Drive
WARSZAWA
2005

WARSAW UNIVERSITY OF TECHNOLOGY
Faculty of Electrical Engineering
Institute of Control and Industrial Electronics
Ph.D. Thesis
M. Sc. Dariusz Świerczyński
Direct Torque Control with Space Vector
Modulation (DTC-SVM) of Inverter-Fed
Permanent Magnet Synchronous Motor Drive
Thesis supervisor
Prof. Dr Sc. Marian P. Kaźmierkowski
Warsaw, Poland - 2005

Contents
Table of Contents
Chapter 1 1
INTRODUCTION
Chapter 2 8
MODELING AND CONTROL MODES OF PM SYNCHRONOUS DRIVES
2.1 Mathematical model of PM synchronous motor 8
2.1.1 Voltage and flux-current equations 9
2.1.2 Instantaneous power and electromagnetic torque 17
2.1.3 Mechanical motion equation 22
2.2 Static characteristic under different control modes 25
2.3 Summary 33
Chapter 3 34
VOLTAGE SOURCE PWM INVERTER FOR PMSM SUPPLY
3.1 Introduction 34
3.2 Voltage source inverter (VSI) 35
3.3 Space vector based pulse width modulation (PWM) methods 46
3.4 Summary 52
Chapter 4 53
CONTROL METHODS OF PM SYNCHRONOUS MOTOR
4.1 Introduction 53
4.2 Field oriented control (FOC) 54
4.3 Direct torque control (DTC) 57
4.4 Summary 64
Chapter 5 65
DIRECT TORQUE CONTROL WITH SPACE VECTOR MODULATION
(DTC-SVM)
5.1 Introduction 65
5.2 Cascade structure of DTC–SVM scheme 66
5.2.1 Digital flux control loop 68
5.2.2 Digital torque control loop 82
5.3 Parallel structure of DTC–SVM scheme 91
5.3.1 Digital flux control loop 92
5.3.2 Digital torque control loop 102
5.4 Speed control loop for DTC–SVM structure control 113
5.5 Summary 122
Chapter 6 121
DIRECT TORQUE CONTROL WITH SPACE VECTOR MODULATION (DTC-
SVM) OF PMSM DRIVE WITHOUT MOTION SENSOR
6.1 Introduction 121
6.2 Initial rotor position estimation method 123
6.3 Stator flux estimation methods 127
6.3.1 Overview 127
6.3.2 Current model based flux estimator 127
6.3.3 Voltage model based flux estimator with ideal integrator 128

Contents
6.3.4 Voltage model based flux estimator with low pas filter 129
6.3.5 Improved voltage model based flux estimator 130
6.4 Electromagnetic torque estimation 132
6.5 Rotor speed estimation methods 132
6.5.1 Overview 132
6.5.2 Back electromotive force (BEMF) technique 133
6.5.3 Stator flux based technique 133
6.6 Summary 136
Chapter 7 137
DSP IMPLEMENTATION OF DTC-SVM CONTROL
7.1 Description of the laboratory test-stand 137
7.2 Steady state behaviour 140
7.3 Dynamic behaviour 143
7.3.1 Flux and torque control loop 143
7.3.2 Speed control loop 151
Chapter 8 161
SUMMARY AND CLOSING CONCLUSIONS
Appendices 163
Picture of rotor and stator of PMSM machine
Basic transformation
Model of PM synchronous motor- SABER
Parameters of PMSM machine
Parameters of voltage source inverter
PI speed controller
PWM technique - overmodulation
List of Symbols 170
References 172

Introduction
Chapter 1 INTRODUCTION
Recently, an increased interest in application of permanent magnet synchronous motors
(PMSM) in speed controlled drives has been observed. This is stimulated mainly by:
• development of modern high switching frequency semiconductor power devices (as
for example IGBT modules of 5-th generation),
• new rare earth magnetic materials as samarium-cobalt (Sm-Co) or neodymium-ironboron
(Nd-Fe-B),
• specialized digital signal processor (DSP) for AC drive applications with integrated
PWM function, A/D converters as well as processing of encoder signals (e.g
ADMC401, TMS320FL24XX, TMS320FL28XX).
Synchronous motors with an electrically excited rotor winding have a conventional threephase
stator winding (called armature) and an electrically excited field winding on the rotor,
which carries a DC current. The armature winding is similar to the stator of induction motor.
The electrically excited field winding can be replaced by permanent magnet (PM) [1]. The use
of permanent magnets has many advantages including the elimination of brushes, slip rings,
and rotor copper losses in the field winding. It leads to higher efficiency. Additionally since
the copper and iron losses are concentrated in the stator, cooling of machines through the
stator is more effective. The lack of field winding and higher efficiency results in reduction of
the machine frame size and higher power/weight ratio.
Figure. 1.1. General classification of AC synchronous motors.
1

Introduction
Generally, the permanent magnet AC machines can be classified into two types (Fig.1.1):
trapezoidal type called “brushless DC machine” (BLDCM) and sinusoidal type called
permanent magnet synchronous machine (PMSM). The BLDC machines operate with
trapezoidal back electromagnetic force (EMF) and require rectangular stator phase current.
The PMSM’s generate sinusoidal EMF and operate with sinusoidal stator phase current.
The PMSM can be further divided into two main groups in respect how the magnet bars have
mounted in the rotor [6,7]. In the first group magnets are mounted in the rotor (Fig. 1.2 c-d)
and this type is called interior permanent magnet synchronous motors (IPMSM). The second
group is represented by surface permanent magnet synchronous motors (SPMSM). In the
SPMSM magnet bars are mounted on the rotor surface (Fig. 1.2 a-b).
q
q
SPMSM
S N
S N d
a )
b)
S
N
S
N
d
q
q
IPMSM
S
N
c )
d)
S
N
d
d
Fig. 1.2. The cross section of the PMSM rotor shaft and the magnet bars placements:
a),b),c) axial field direction, d) radial field direction.
The magnets can be placed in many ways on the rotor (Fig. 1.2). In radial field fashion the
magnet bars are along the radius of the machine and this arrangement provides the highest air
gap flux density, but it has the drawback of lower structural integrity and mechanical
robustness. Machines with this arrangement of magnets are not preferred for high-speed
applications (higher than 3000 rpm). In axial field manner the magnets are placed parallel to
the rotor shaft. This arrangement of magnets is much more robust mechanically as compared
2

Introduction
to surface-mounted machine. It makes possible to use IPMSM for higher-speed applications
(contrary to SPMSM’s).
Regardless of the fashion of mounting the PM, the basic principle of motor control is the same
and the differences are only in particularities. An important consequence of the method of
mounting the rotor magnets is the difference in direct and quadrature axes inductance values.
The direct axis reluctance is greater than the quadrature axis reluctance, because the effective
air gap of the direct axis is multiple times that of the actual air gap seen by the quadrature
axis. As consequence of such an unequal reluctance, the quadrature inductance is higher than
direct inductance L q
> L d
. It produces reluctance torque in addition to the mutual torque.
Reluctance torque is produced due to the magnet saliency in the quadrature and the direct axis
magnetic paths. Mutual torque is produced due to the interaction of the magnet field and the
stator current. In case where the magnets bars are mounted on the rotor surface the quadrature
inductance is equal direct inductance L q
= L d
, because of the same flux paths in d and q axis.
As result the reluctance torque disappears.
Among the main advantage of PM machines are [12]:
• high air gap flux density,
• higher power/weight ratio,
• large torque/inertia ratio,
• small torque ripples,
• high speed operation,
• high torque capability (quick acceleration and deceleration),
• high efficiency and high cosφ (low expense for the power supply),
• compact design.
Thanks to this advantages the PMSM’s are usually used in high performance servo drives, in
special applications as computer peripheral equipment, robotics, ect. However, recently the
PMSM are also used as adjustable–speed drives in variety of application such as fans, pumps,
compressors, blowers. Another area is automotive application as an alternative drive in hybrid
mode with classical engine. The power of offered synchronous motors is in the range several
kW to MW.
3

Introduction
The main requirements for high performance PWM inverter-fed PMSM drive can be
formulated as follows:
• operation with and without mechanical motion sensor,
• fast flux and torque response,
• available maximum output torque in wide range of speed operation region,
• constant switching frequency,
• uni-polar voltage PWM,
• low flux and torque ripples,
• robustness to parameters variation,
• four quadrant operation.
To meet the above requirements, different control methods can be used [3,4,10].
Variable
Frequency
Control
Scalar based
controllers
Vector based
controllers
V/Hz=const
with stabilization
loop
Field
Oriented
(FOC)
Direct Torque
Control
(DTC)
PM (rotor)
Flux Oriented
(RFOC)
Stator Flux
Oriented
(SFOC)
Direct Torque
Control with Space
Vector Modulation
(DTC-SVM)
Circular flux
trajectory
(Takahashi)
Figure 1.3 Classification of PMSM control methods.
The general classification of the variable frequency control for PMSM is presented in Fig. 1.3.
The PMSM control methods can be divided into scalar and vector control. According to [3],
in scalar control, which based on a relation valid for steady states, only the magnitude and
frequency (angular speed) of voltage, currents, and flux linkage space vectors are controlled.
Thus, the control system does not act on space vector position during transient. Therefore, this
control is dedicated for application, where high dynamics is not demanded. Contrary, in
4

Introduction
vector control, which is based on relation valid for dynamics states, not just magnitude and
frequency (angular speed), but also instantaneous position of voltage, current and flux space
vectors are controlled. Thus, the control system adjust the position of the space vectors and
guarantee their correct orientation for both steady states and transients.
The scalar constant V/Hz control for PMSM without damper winding (squire cage) is not
simple as for induction motor. It requires additional stabilization control loop, which can be
provide by feedback from: rotor velocity perturbation, active power or DC-link current
perturbation [9].
The most popular vector control method developed in 70s, known as field oriented control
(FOC) [31] gives the permanent magnet synchronous motor high performance. In this method
the motor equation are transformed in a coordinate system that rotates in synchronism with
permanent magnet flux. It allows separately and indirectly control flux and torque quantities
by using current control loop with PI controllers like in well known DC machine control [3].
In search of a simpler and more robust high performance control system in 80s new vector
control called direct torque control (DTC) was developed [50]. It was innovative studies at
this time and completely different approach which depart from the idea of coordinate
transformation and the analogy with DC motor control. It allows direct control flux and torque
quantities without inner current control loops. Using bang-bang hysteresis controllers for flux
and torque control loops made this control concept very fast and not complicated. However,
the main disadvantage of DTC is fast sampling time required and variable switching
frequency, because of hysteresis based control loops. In order to eliminate above
disadvantages and kept basic control rules of classical DTC, at the beginning of 90’s a new
developed control technique called direct torque control with space vector modulator (DTC-
SVM) has been introduced [54,55]. However, from the formal consideration this method can
also be viewed as stator flux oriented control (SFOC). This control employed instead of
hysteresis controller as for classical DTC, the PI controllers and space vector modulator
(SVM). It allows to achieve fixed switching frequency, what considerably reduce switching
losses as well as torque and current ripples. Also requirement of very fast sampling time is
eliminated [113,115,117]. Therefore, this new method is subject of this thesis. In spite of
many control strategies there is no one which may be considered as standard solution.
5

Introduction
Therefore, the following thesis can be formulated:
“In the view of commercial manufacturing process the most convenient control scheme
for voltage source inverter-fed permanent magnet synchronous motor (PMSM) drives is
direct torque control with space vector modulator (DTC-SVM)”.
To prove the above thesis, the author used methodology based on an analyze and simulation
as well as experimental verification on the laboratory setup with 3kW PMSM motor.
Moreover, the presented control algorithm DTC-SVM has been introduced and used in serial
commercial product of Polish manufacture TWERD, Toruń.
In the author’s opinion the following results of the thesis are his original achievements:
• development of a simulation algorithm in SABER package for the investigation
of PWM inverter-fed PMSM control,
• elaboration and experimental verification of digital flux and torque controller
design based on the Z-transform approach for series (cascade) and parallel
structure of DTC-SVM schemes,
• implementation and verification of series (cascade) and parallel DTC-SVM
schemes on experimental laboratory setup with 3kW PM synchronous motor
drive controlled by floating point DS1103 board.
• bringing into production and testing of developed DTC-SVM algorithm in Polish
industry.
The thesis consists of eight chapters. Chapter 1 is an introduction. In Chapter 2 mathematical
model of PM synchronous motor and his basic control modes are presented. Chapter 3 is
devoted to voltage source inverter, his nonlinear characteristics and different PWM
techniques. Chapter 4 gives brief review of PM synchronous motor control method such as
FOC and classical DTC. In Chapter 5 two kind of DTC-SVM control schemes are presented.
Also, the analysis and synthesis of digital flux, torque and speed controllers based on Z
transform approach are given. Chapter 6 is devoted to initial rotor detection methods, stator
flux vector and rotor speed estimation algorithms. In Chapter 7 experimental results are
6

Introduction
presented and studied. Chapter 8 includes the finally conclusions. Description of the SABER
based control algorithm, basic coordinate transformations and parameters of used PM
synchronous machine as well as inverter are given in Appendices.
7

Modeling and control modes of PM synchronous motor drives
Chapter 2
MODELING AND CONTROL MODES OF PM SYNCHRONOUS DRIVES
2.1 Mathematical model of PM synchronous motor
Development of the machine model through the understanding of physics of the
machine is the key requirement for any type of electrical machine control. Since in this
project a Surface type Permanent Magnet Synchronous Motor (SPMSM) is used for the
investigation [9,13,14,15,16]. The development of those models is under bellow
assumptions as [3]:
• three-phase motor is symmetrical,
• only a fundamental harmonic of the magneto motive force (MMF) is taking in to
account,
• the spatially distributed stator and rotor winding are replaced by a concentrated
coil,
• an anisotropy effects, magnetic saturation, iron loses and eddy currents are not
taking into considerations,
• the coil resistances and reactances are taking to be constant,
• in many cases, especially when is considered steady state, the currents and
voltages are assumed to be sinusoidal,
• thermal effect for permanent magnets is omitted.
The synchronous motor model will be presented in space vector notation. Space vector
form of the machine equations has many advantages such as compact notation, easy
algebraic manipulation, and very simple graphical interpretation. Specially, this notation
is very useful when analyzing the vector control based technique of the AC machines.
The space vector representation of AC machine equations has been discussed in detail
in number of text books ([3,4,12]).
The instantaneous value of a three-phase system KA, KB,
K
C
(such as currents, voltages
and flux linkages) can be replaced by one resultant vector called the space vector,
2
K = ⎡ 1 ⋅ K + a⋅ K + a K
3 ⎣
2
A B C
⎤
⎦
(2.1)
8

Modeling and control modes of PM synchronous motor drives
2π
2π
4π
j 1 3
where:1,
3 2
− j j 1 3
a= e =− + j , a = e
3
= e
3
=− − j - complex vectors, 2/3 –
2 2
2 2
normalization factor (guarantee that for balanced sinusoidal waveforms the magnitude
of the space vector is equal to the amplitude of that phase waveforms).
The elements of this space vector satisfy the condition:
KA + KB + KC
= 0
(2.2)
and it means that we have three-phase system without neutral wire.
2.1.1 Voltage and current equations
For idealized motor (Fig. 2.1), the following equations of the instantaneous stator phase
voltages can be written [3]:
B
b
Z sB
I sB
a
U sB
S
N
N
S
U sA
γ m
Z sA
A
I sA
I sC
Z sC
U sC
N
S
C
c
Figure 2.1. Layout and symbols for three-phase PMSM electric motor windings.
dΨ
dt
sA
sA
=
sA sA
+ (2.3a)
U I R
U I R
dΨ
dt
sB
sB
=
sB sB
+ (2.3b)
U I R
dΨ
dt
sC
sC
=
sC sC
+ (2.3c)
9

Modeling and control modes of PM synchronous motor drives
where
U ,
sA, U
sB
U
sC
are the instantaneous stator voltage values, I
sA
I
sB
, I
sC
, are
instantaneous values of the current,
R = R = R = R is the resistance of the stator
s
sA
sB
sC
windings, and ΨsA,
Ψ
sB
and Ψ
sC
are magnetic flux linkages stator windings A, B and
C , respectively.
Using the space vector theory to voltage equations we can written in vector form
where:
U
d Ψ
dt
sABC
sABC
= Rs
I
sABC
+ (2.4)
2 2
U
sABC
= (1 UsA + aUsB + a UsC
) ,
2 2
2 2
I
sABC
= (1 IsA + aIsB + a IsC
) ,
sABC (1
sA
a
sB
a
sC )
3
3
3
stator voltage, current and flux space vectors, respectively.
Ψ = Ψ + Ψ + Ψ are the
The stator winding flux consist of rotor flux and stator flux linkages:
where,
Ψ
sABC
=Ψ
ABC ( s) +Ψ
ABC( r)
(2.5)
⎡ LsA MsAB MsAC⎤⎡IsA⎤
⎢
M L M
⎥⎢
I
⎥
Ψ
ABC( s)
= ⎢ sBA sB sBC ⎥⎢ sB ⎥
⎢⎣ M
sCA
MsCB L ⎥⎢
sC ⎦⎣I
⎥
sC ⎦
(2.6)
⎡
⎤
⎢ cosθ
⎥
r
⎢
⎥
2π
Ψ
ABC ( r)
=Ψ ⎢
PM
cos( θr
− ) ⎥
⎢ 3 ⎥
⎢
2π
⎥
cos( θr
+ )
⎢⎣
3 ⎥⎦
(2.7)
and, θr
is electrical rotor position. Mechanical rotor position is defined as:
θ = p γ
(2.8)
r b m
where: pb
- number of pole pairs, γ
m
- mechanical position.
In equation (2.6) LsA
is the self-inductance of phase A winding, M
sAB
and M
sAC
are the
mutual inductances between A and B phase, A and C phase, respectively. For self and
mutual inductances of B and C phase the same notations used. In (2.7),
Ψ
PM
is the
10

Modeling and control modes of PM synchronous motor drives
amplitude of the flux linkages established by the permanent magnet on the rotor. The
inductances are described below.
Due to the rotor saliency in IPMSM the air gap is not uniform and, therefore, the self
and mutual inductances of stator windings are a function of the rotor position.
The derivation of these rotor position dependent inductances is available in details in
[5]. The results are summarized here as follows:
The stator winding self-inductances are
L = L + L − L cos 2θ
(2.9a)
sA ls A B r
2π
4π
LsB = Lls + LA −LB cos2( θr − ) = Lls + LA −LB cos(2 θr
− )
3 3
(2.9b)
2π
4π
LsC = Lls + LA − LB cos 2( θr + ) = Lls + LA − LB cos(2 θr
+ )
3 3
(2.9c)
where, Lls
is stator-winding leakage inductance and LA,
LB
are given by
L
A
⎛ms
⎞
= ⎜
2 ⎟
⎝ ⎠
2
πµ rlε
0 1
(2.10a)
L
B
1 ⎛ms
⎞
=
2 ⎜
2 ⎟
⎝ ⎠
2
πµ rlε
0 2
(2.10b)
where,
m
s
is number of turns of each phase winding, r is radius, which is from center
of machine to the inside circumference of the stator, and l is the axial length of the air
gap of the machine, µ
0
is permeability of the air, ε
1
and ε
2
are defined as s:
1 1 1
ε
1
= ( + )
(2.11a)
2 g g
min
max
1 1 1
ε
2
= ( − )
(2.11b)
2 g g
min
max
where,
gmin
is minimum air gap length and
max
g
is maximum air gap length.
The mutual inductances between stator phase are:
1 π 1 2π
MsAB = MsBA =− LA −LB cos 2( θr − ) =− LA −LB cos(2 θr
− ) (2.12a)
2 3 2 3
11

Modeling and control modes of PM synchronous motor drives
1 π 1 2π
MsAC = MsCA =− LA − LB cos 2( θr + ) =− LA − LB cos(2 θr
+ ) (2.12b)
2 3 2 3
1 1
MsBC = MsCB =− LA − LB cos 2( θr + π) =− LA − LB cos(2θr
+ 2 π)
2 2
1
=− LA −LBcos 2θ
r
2
(2.12c)
Using the space vector theory, the flux linkage
Ψ
sABC
space vector can be written as:
3 3
∗ j2θr
jθr
Ψ
sABC
= ( Lls + LA)
IsABC − LB IsABC
e +Ψ
PM
e
(2.13)
2 2
where,
2 2
(1
sA sB sC )
IsABC
= I + aI + a I ,
3
space vector and conjugate stator current space vector.
2 2 (1
sA sB sC )
I ∗ sABC
= I + a I + aI are the stator current
3
Taking into account that:
Ld = Lls + Lmd
(2.14a)
Lq = Lls + Lmq
(2.14b)
3 3
where, Lmd = ( LA + LB
) , Lmq = ( LA − LB
) are d and q magnetizing inductances and
2
2
are defined as [5].
Finally, equations (2.13) comes as:
where, L
d
, L
q
are d and q inductances.
Ld + Lq Lq −Ld ∗ j2θr
jθ
r
Ψ
sABC
= ( ) I
sABC
− ( ) IsABC
e +Ψ
PM
e
(2.15)
2 2
Space vector form of machine equations (2.4, 2.15) becomes more compact, but the
rotor position dependent parameters still exist in that form of expressions for the stator
flux linkage space vector. Therefore, the space vector model is still not simple to use for
the analysis. A simplification can be made if the space vector model is referred to a
suitably selected rotating frame.
12

Modeling and control modes of PM synchronous motor drives
Figure 2.2 shows axes of reference for the three-stator phase ABC. , , It also shows a
rotating set of x,
y axes, where the angleθ K
is position of x -axis in respect to the stator
A phase axis. Variables along the AB , and C axes can be referred to the x − and
y − axes by the expression:
⎡K
A ⎤
⎡K
x ⎤ 2 ⎡ cosθK cos( θK − 2 π / 3) cos( θK
+ 2 π / 3) ⎤⎢
K
⎥
⎢ B
K
⎥ =
y 3
⎢
sinθK sin( θK 2 π / 3) sin( θK
2 π / 3)
⎥⎢ ⎥
⎣ ⎦ ⎣− − − − + ⎦
⎢⎣K
⎥
C ⎦
(2.16)
y
K B
K ABC
K y
K x
Ω K
x
θ K
K A
K C
Figure 2.2. Stator fixed three phase axes (A,B,C) and general rotating reference frame ( x,
y ).
Finally, the space vector in general rotating frame can be written as:
j K
K = K (cosΘ + jsin Θ ) = K e θ
(2.17)
ABCs K K
K K
In this case the voltage equation (2.4) using (2.17) can written as:
jθK jθ d
K jθK
U
sKe Rs IsKe (
sKe
)
dt
= + Ψ (2.18)
Using chain rule, equation. (2.17) and divided by term
j K
e θ
can be written as:
d Ψ
U
sK
R I j
dt
= + + Ω Ψ (2.19)
sK
s sK K sK
where
U
sK
rotating frame.
, I sK , Ψ sK is the stator voltage, current and flux space vector in general
Making similar arrangement like for the voltage equation the flux linkage vector in
general reference frame can be expressed as:
13

Modeling and control modes of PM synchronous motor drives
Ld + Lq Lq −Ld j2( θr θK) j( θr θK)
Ψ
sK
= ( ) IsK − ( ) I ∗ sK
e − +Ψ
PM
e
− (2.20)
2 2
Stator fixed system ( α,
β )
Taking the angular speed of the reference frame to be Ω
K
= 0 and θ
K
= 0 , the set of
synchronous machine vector equations (2.19) and (2.20) my be written as:
U
d Ψ
sαβ
sαβ
= Rs
Isαβ
+ (2.21)
dt
Ld + Lq Lq −Ld ∗ j2θr
jθ
r
Ψ
sαβ = ( ) I
sαβ − ( ) Isαβ
e +Ψ
PM
e
(2.22)
2 2
Substituting to above equations the following expressions for complex vectors
U U jU
= s sα
+ , αβ sβ
s sα
sβ
I αβ = I + jI , Ψ αβ =Ψ
sα
+ jΨ sβ
and splitting into real and
s
imaginary parts one can obtain the scalar form of the machine equations in stationary
α,
β reference frame:
U
U
dΨ
dt
sα
sα
= RsIsα
+ (2.23a)
dΨ
sβ
sβ
= RsIsβ
+ (2.23b)
dt
Ld + Lq Lq −Ld Lq −Ld
( cos2 θ ) Is
( )sin2θ Is
cosθ
2 2 2
Ψ = − − +Ψ (2.24a)
sα r α r β PM r
Lq − Ld Ld + Lq Lq −Ld
Ψ
sβ =− ( )sin 2 θrIsα + [( ) + ( )cos2 θr] Isβ
+Ψ
PM
sinθ
(2.24b)
r
2 2 2
Note, that in the flux-current equations (2.24a and b) still we can observe that value of
inductances depends on rotor position θ
r
.
Stator flux fixed system ( x,
y )
In order to take advantage of the set of equations (2.19) and (2.20) in rotating coordinate
system, one assumes that the coordinate system rotates with the stator flux linkage
angular speed Ω
K
=Ω
Ψ
and θK
= θ Ψ
. As a
sx s
Ψ =Ψ , δ = −( θ − θ )
Ψ
r
Ψ
d Ψ
U R I j
dt
sxy
sxy
=
s sxy
+ + ΩΨsΨ (2.25)
sxy
Ld + Lq Lq −Ld ∗ − j2δΨ
− jδΨ
Ψ
sxy
= ( ) Isxy − ( ) Isxy
e +Ψ
PM
e
(2.26)
2 2
14

Modeling and control modes of PM synchronous motor drives
Substituting to above equations the following expressions for complex vectors
U U jU
sxy
=
sx
+
sy, Isxy
Isx jIsy
= + , Ψ =Ψ
s
and splitting into real and imaginary parts
sxy
one can obtain the scalar form of the machine equations in stationary x,
y reference
frame:
U
dΨ
dt
s
sx
= Rs Isx
+ (2.27a)
U = R I +ΩΨ
Ψ (2.27b)
sy s sy s
s
1 1
Ψ
s
= [( Ld + Lq) −( Lq − Ld)cos 2 δΨ ] Isx
+ ( Lq − Ld)sin 2δΨIsy
+Ψ
PM
cosδΨ
2 2
(2.28a)
1 1
0 = ( Lq − Ld)sin2 δΨIsx
+ [( Ld + Lq) + ( Lq −Ld)cos2 δΨ] Isy
−Ψ
PM
sinδΨ
2 2
(2.28b)
The current-flux equations can be expressed also in simplest form as:
2 2
s
( Ld cos δΨ Lqsin δΨ) Isx
( Lq Ld)sinδΨcosnδΨIsy
PM
cosδΨ
Ψ = + + − +Ψ (2.29a)
2 2
q d
δΨ δΨ sx d
δΨ q
δΨ sy PM
δΨ
0 ( L L )sin cos I ( L sin L cos ) I sin
= − + + −Ψ (2.29b)
Rotor flux fixed system ( dq) ,
In order to take advantage of the set of equations (2.19) and (2.20) in rotating coordinate
system, one assumes that the coordinate system rotates with the rotor flux angular speed
Ω
K
= pbΩ m
and θK = pbγm = θr
d Ψ
U
sdq
R I jp
dt
sdq
s sdq
b m
= + + Ω Ψ (2.30)
sdq
Ld + Lq Lq −Ld
∗
Ψ
sdq
= ( ) Isdq − ( ) Isdq
+Ψ
PM
(2.31)
2 2
Substituting the following expressions for complex vectors U = Usd + jUsq,
I I jI
sdq
=
sd
+
sq, sdq sd
j
sq
Ψ =Ψ + Ψ to (2.30) and (2.31), and splitting for real and
imaginary parts the scalar form of the machine equations in rotational fixed reference
frame can be obtained:
dΨ
sd
Usd = Rs Isd + − pbΩmΨ sq
(2.32a)
dt
sdq
15

Modeling and control modes of PM synchronous motor drives
dΨsq
Usq = Rs Isq + + pbΩmΨ sd
(2.32b)
dt
where,
Ψ
sd
= LI
d sd
+Ψ (2.33a)
PM
Ψ
sq
= LqI
(2.33b)
sq
It should be noted that when transforming the flux linkage vector
Ψ s
to the dq ,
reference frame the rotor position θ r
dependent terms disappear it can be seen from
equation (2.31). This is the main advantage of rotor-oriented representation.
Substituting the relationship of (2.33a-b) into (2.32a-b), and also considering
dΨ PM = 0 , the most common scalar form of the machine voltage equations in the rotor
dt
reference frame can be obtained as:
dIsd
Usd = RsIsd + Ld − pbΩ mLqIsq
(2.34a)
dt
dIsq
Usq = Rs Isq + Lq + pbΩmΨ PM
+ pbΩ mLd Isd
(2.34b)
dt
Based on the above voltage-current equations it is possible to draw the equivalent
electrical circuit separately for d and q axes (Fig. 2.3).
R s
p Ω
L I
b m q sq
R s
p Ω
L I
b m d sd
U
sd
I sd
L d
U
sq
I sq
p Ω Ψ
b m PM
L q
Figure 2.3. Equivalent circuit model of PMSM in the rotor reference frame. (a) Rotor d-axis
equivalent circuit, (b) Rotor q-axis equivalent circuit.
16

Modeling and control modes of PM synchronous motor drives
2.1.2 Instantaneous power and electromagnetic torque
The three-phase star-connection system without neutral wire is shown in Fig. 2.4. This
is classical configuration for AC motor windings connections.
A
I sA
U sAC
U sAB
B
U sA
Z sC
Z sA
U sAB
Z sB
C
U sBC
I sC
U sC U sB
I sB
Figure 2.4. Three-phase star connection system without neutral wire.
For this configuration the expression for instantaneous active power supplied to load
can be expressed as:
P= UsAIsA+ UsBIsB+ UsCIsC
(2.35)
Introducing space vector definition, after some arrangement and taking into account the
relation: I + I + I = 0, the equation (2.35) can be written as:
sA sB sC
3 ∗
P= Re[ U
sABC
IsABC
]
(2.36)
2
For dq , frame, the equation (2.35) for the active power can be expressed as:
3
P= ( UsdIsd + UsqIsq
)
(2.37)
2
Substituting voltage equation (2.4) into (2.36), and adopting Ω
K
= pbΩmone obtains
3 ∗ d Ψ
[Re(
sABC ∗ ∗
P= RsIsABC IsABC + IsABC − jpbΩmΨ sABCIsABC
)] (2.38)
2
dt
Note that
sABC
sABC
2
s
I I ∗ = I and:
3 2 d Ψ
[ Re(
sABC ∗
∗
P= Rs Is + IsABC ) + Re( −jpbΩmΨ sABC
IsABC
)] (2.39)
2
dt
d ΨsABC
Hence, neglecting the losses in stator resistance R
s
and assuming that = 0 , the
dt
electromagnetic power is expressed:
17

Modeling and control modes of PM synchronous motor drives
3
∗
Pe = pbΩmIm( Ψ
sABC
IsABC
)
(2.40)
2
In dq , frame the active power can be written:
3
Pe = pbΩm( ΨsdIsq −Ψ
sqIsd)
(2.41)
2
For the presented system (Fig. 2.4) the expression for instantaneous reactive power
supplied to the three-phase load system without neutral wire can be calculated as:
1
Q= ( IsAUsBC + IsBUsCA + IsCUsAB
)
(2.42)
3
Introducing the space vector definition into equation (2.42), after some arrangement,
and taking into account the relation: I + I + I = 0, one obtains:
sA sB sC
3 ∗
Q= Im[ U
sABC
IsABC
]
(2.43)
2
In dq , frame the reactive power is expressed as:
3
Q= ( UsqIsd − UsdIsq
)
(2.44)
2
Substituting voltage equation (2.4) into (2.43), adopting Ω
K
= pbΩ m
and made similar
arrangements like for active power calculation, the final expression for reactive power
is:
3
∗
Q= pbΩmRe( Ψ
sABC
I
sABC
)
(2.45)
2
In dq , frame the expression (2.45) for the reactive power becomes:
3
Q= pbΩm ( Ψ
sd
Isd +Ψ
sqIsq
)
(2.46)
2
The important quantity of the drive is the power factor cosφ , which can be calculated
as:
Q
cosφ = (2.47)
S
where S is module of apparent power vector S = P+ jQ:
2 2
S = P + Q
(2.48)
18

Modeling and control modes of PM synchronous motor drives
The instantaneous electromagnetic torque developed by an electric motor can be defined
as:
P
e
M
e
= (2.49)
Ω
m
where, P
e
is the electromagnetic power and
Ωm
is the mechanical angular rotor speed.
Finally, taking into account equation (2.49) the expression for electromagnetic torque
can be obtained as:
and in dqframe: ,
3
∗
Me = pbIm( Ψ
sABC
IsABC
) , (2.50)
2
3
M
e
= pb ( Ψsd Isq −Ψ
sqIsd
)
(2.51)
2
Substituting Ψ , Ψ from (2.33a-b), the torque expression of equations (2.47)
becomes:
sd
sq
3
M
e
= pb ( ΨPM Isq −( Lq − Ld ) Isd Isq
)
(2.52)
2
It can be seen from (2.52), that developed torque consist of two parts, one produced by
the permanent magnet flux called synchronous torque ( M
reluctance torque ( M
er
es
) and the second called
), which is produced by the difference of the inductance in rotor
d- and q-axes. Expressions for those two torque components are:
3
M
es
= pbΨ PM
Isq
(2.53a)
2
3
M
er
=− pb ( Lq − Ld ) Isd Isq
(2.53b)
2
It should be mentioned that for SPMSM ( L d
= L ) the reluctance torque does not exist
q
due to the same inductance paths in rotor d- and q-axes.
The torque expression (2.52) can also be written in polar form using the current vector
amplitude
vector (Fig. 2.5.).
I
s
and the torque angle δ I
, i.e. angle between rotor d-axis and current
19

Modeling and control modes of PM synchronous motor drives
q−
axis
I sq
I s
Ω s
δ I
d
− axis
I sd
Ψ PM
Figure 2.5. Stator current vector in rotor reference frame.
For two current components using trigonometrical rules we can write:
I = I cosδ
(2.54a)
sd
s
I
I = I sinδ
(2.54b)
sq
s
I
Substituting I , I into equation (2.52), the torque expression can be obtain as:
sd
sq
M 3 1
2
e = b[ PM sin I ( q d) sin2 I]
2 p Ψ I s δ −
2
L − L I s δ
(2.55)
M
es
For given current amplitude the synchronous and reluctances torque varies according to
the sine of torque angle δ
I
. The variation of
M
er
M
es
and M
er
and resultant torque M
e
with
torque angle are illustrated in Fig. 2.6. The IPMSM parameters used for this calculation
are given in the Appendices.
e [ ] M Nm
M [ ] er
Nm M
es
[ Nm ]
[deg] δ I
Figure 2.6. Variation of synchronous torque M es
, reluctance torque M
er
and resultant
torque M as a function of torque angle (for rated current amplitude).
e
20

Modeling and control modes of PM synchronous motor drives
Referring to Fig. 2.7, stator flux components in rotor reference frame can be written as:
where:
Ψ = Ψ cos = LI +Ψ (2.56a)
sd s
δ Ψ d sd PM
Ψ = Ψ sin = LI
(2.56b)
sq s
δ Ψ q sq
Ψ
s
is stator flux linkage amplitude, Ψ
PM
is rotor permanent magnet and δ Ψ
is
torque angle (angle between stator flux linkage vector and rotor permanent magnets flux
vector).
q−
axis
I sq
Ψ sd
Ψ PM
I s
Ω s
Ψ s
Ψ
sq
= LqI
sq
LI
d
sd
δ Ψ
I sd
d
− axis
Ψ
sd
= LI
d sd
+ΨPM
Figure 2.7. Rotor permanent magnet flux vector and stator flux linkage vector in rotor reference
frame.
From (2.56a) and (2.56b) the I
sd
and Isq
can be obtained as:
I
I
sd
sd
Ψs
cosδ Ψ
−ΨPM
= (2.57a)
L
d
Ψs
sinδ = Ψ
(2.57b)
L
q
Substituting current components (2.57a), (2.57b) into equation (2.51), one can obtain
another useful torque expressions:
M
e
3 Ψ sin
2 s
ΨPM
δΨ Ψs
( Lq
−Ld)sin2δ
Ψ
= pb[ −
]
2 Ld 2LdL
q
Mes
M
er
(2.58)
where:
Ψ
s
stator flux linkage amplitude, and Ψ
PM
rotor flux, δ Ψ
is torque angle, M
es
-
synchronous torque,
M
er
- reluctance torque.
21

Modeling and control modes of PM synchronous motor drives
For the PM synchronous motor the amplitude of stator flux
Ψ
s
is established by
permanent magnet. Operation with stator flux amplitude belong the nominal value of
rotor flux amplitude
Ψ
PM
increases the amplitude of stator phase current. Please note
that maximum amplitude of the stator current vector is calculated as:
higher value may damage the PM (complete demagnetization).
I
s
Ψ
≤
L
PM
d
, and
From the Fig. 2.8 it can be observed that rated torque is achieved for torque angle
0< < 25
electrical degree.
δ Ψ
e[ ] M Nm
M
er
[ Nm ]
M
es
[ Nm ]
δ Ψ
[deg]
Figure 2.8. Variation of synchronous torque M es
, reluctance torque M
er
and resultant
torque M as a function of torque angle (for constant stator flux equal value of PM).
e
2.1.3 Mechanical motion equation
The equation of rotor motion dynamics describes the mechanical equilibrium of a drive
system. Taking the moment of inertia to be constant ( J = const.
) and neglecting friction
and elastic torque we can write:
where,
M
e
= Ml + Md
(2.59)
Ml
is the external torque on the motor shaft, and
M
dΩ
dt
where: J is total moment of inertia,
M
d
the dynamic torque
m
d
= J
(2.60)
Ω
m
angular speed of the rotor.
22

Modeling and control modes of PM synchronous motor drives
In general, for a drive system,
J = J + J
(2.61)
m
where: J
m
- motor inertia, J
l
load moments of inertia.
l
From equation (2.55) and (2.56) one can write:
and, with (2.50),
dΩ m 1 = ( M
e − M
l )
(2.62)
dt J
dΩ m 1 3
∗
= ( pb
Im( Ψ
s
Is
) − Ml
)
(2.63)
dt J 2
Finally, the full mathematical model of PM synchronous machine which is used in
simulation studies [Appendices] is described in dq , reference frame as:
dΨ
sd
Usd = Rs Isd + − pbΩmΨ sq
(2.64a)
dt
dΨ
sq
Usq = Rs Isq + + pbΩmΨ sd
(2.64b)
dt
Ψ
sd
= LsdIsd +Ψ
PM
(2.65a)
Ψ
sq
= LsqIsq
(2.65b)
dΩ m 1 = ( M
e − M
l )
(2.66)
dt J
3 ∗ 3
M
e
= pbIm( Ψ
s
Is) = pb( ΨsdIsq −Ψ
sqIsd)
(2.67)
2 2
Based on above equations we can create the block scheme of the PMSM machine (Fig.
2.9), where the input signals are the voltage components in dq , reference frame
U , U and the output signal is the mechanical speed of the rotor Ω
m
. As the external
sd
sq
load torque
M
l
is disturbance.
23

Modeling and control modes of PM synchronous motor drives
M l
U sd
+
−
pbΩmΨsd
∫
Ψ sd
Ψ PM
−
1
L d
I sd
Ψ
sq
I
sd
p b
R s
−
3
p
2 b
M e
−
1
J
∫
Ω m
R s
U sq
−
−
∫
Ψ sq
1
L q
I sq
Ψ
sd
I
sq
pbΩmΨsq
p b
Figure 2.9. Block scheme of PM synchronous machine in rotating dq , frame.
Based on equations (2.64-2.67) we can also draw the vector diagram of PM
synchronous motor (Fig. 2.10). From this vector representation it can see the positions
of the vectors (currents, voltages and fluxes). Especially, power angle φ (angle between
voltage and current vectors) and torque angle defined in two manners: as an angle
between current and rotor flux vectors -δ I
, or as angle between stator flux and rotor
flux vectors -δ Ψ
.
q−axis
β
U s
RI
s
s
I s
Ψ sq
I sq
LI
d sd
ΩΨ
s
s
φ
Ψ s
LqI
sq
I sd
δ I
δ Ψ
θ r
Ψ sd
θ Ψs
Ψ PM
d −axis
rotor
α ( A)
stator
Figure 2.10. Vector diagram of PM synchronous motor in rotor reference frame dq. ,
24

Modeling and control modes of PM synchronous motor drives
2.2 Static characteristic under different control modes
In this section, basic steady state properties of the PMSM under different control mode
strategies will be study [6,9]. The key control strategies for the PMSM can be listed as
follows:
• Constant torque angle control (CTAC).
• Maximum torque per ampere control (MTPAC)
• Unity power factor control (UPFC)
• Constant stator flux control (CSFC)
Constant torque angle (CTA) control
This control strategy for PMSM keeps the torque angle δ
I
(angle between stator current
vector and rotor permanent magnet flux) at constant value 90 .
q−
axis
Isq
= I
s
δ = 90
I
d
− axis
Figure 2.11. Current vector and permanent magnet flux vector for constant torque angle
operation (CTAC)
Hence, this control can be achieved by controlling the d-axis current components to
zero leaving the current vector on the rotor q-axis (see Fig. 2.11). Therefore, this
strategy is also referred to as I
sd
= 0 control. The amplitude of rotor flux vector is
constant and also the torque angle is constant. So, the torque depends only on the value
of stator current amplitude. Therefore, this control strategy is not recommended for
IPMSM with high saliency ratio. However, for SPMSM, this strategy is commonly
used.
Ψ PM
The torque equation in this mode of operation becomes:
3 3
M
e
= pbΨ PMIsq = pbΨ PM
Is
(2.68)
2 2
25

Modeling and control modes of PM synchronous motor drives
The steady state voltage components based on the equations (2.34a) and (2.34b) are:
Usd =−pbΩ mLq Iqs =−pbΩ mLq Is
(2.69a)
U = R I + p Ω Ψ = R I + p Ω Ψ (2.69b)
sq s qs b m PM s s b m PM
The amplitude of stator voltage vector can be calculated as:
s
2 2
sd sq
U = U + U
(2.70)
The stator flux vector amplitude can be calculated from equations (2.65a-b) as:
s
2 2
sd sq
Ψ = Ψ +Ψ (2.71)
The active and reactive power and also the power factor can be obtained from equations
(2.41),(2.46), (2.47).
Maximum torque per ampere (MTPA) control
The main idea of this control is develop the torque using minimum value of stator
current amplitude. In this case the I sd
components is not equal zero, and may cancel the
reluctance torque produced by high saliency ratio. Therefore, this control strategy is
recommended for IPMSM.
q−
axis
I s
I sq
δ >= 90
I
d
− axis
I sd
Figure 2.12. Current vector I s and permanent magnet flux vector Ψ
PM
for maximum torque
per ampere operation (MTPAC).
In order to obtain the maximum torque per ampere we should solve the derivative of
torque equations (2.55) in respect to torque angle. Solving for torque angle α and taking
into account that only negative sign should be considered for the solution, we can
calculate torque angle as:
Ψ PM
−1 −1 1 1 2
δ
I
= cos [ − + ( ) ]
4( L −L ) I 2 4( L −L ) I
d q s
d q
s
(2.72)
26

Modeling and control modes of PM synchronous motor drives
From Fig. 2.8, it can be seen that
M is maximum when torque angle is 90 < < 180
e
The relevant torque equation in this mode of operation becomes from (2.55).
δ
I
.
The steady state voltage equations can be written using the current vector amplitude
I
s
and the torque angle δ
I
as:
U = R I cosδ
+ p Ω L I sinδ
(2.73a)
sd s s I b m q s I
U = R I sinδ
− p Ω L I cosδ
+ p Ω Ψ (2.73b)
sq s s I b m d s I b m PM
The amplitude of stator voltage vector can be calculated from equation (2.70) and
amplitude of stator flux vector from (2.71). The active and reactive power and also the
power factor can be obtained from equations (2.41),(2.46), (2.47).
Unity power factor (UPF) control
Under this control strategy there is no phase different between the current vector and the
voltage vector. Hence, power factor angle φ (see Fig. 2.13) becomes zero. Since only
active power is supplied to the machine under unity power factor operation, the VA
rating requirement of the inverter can be reduced.
q − axis
U s
φ = 0
I s
δ I
d − axis
Figure 2.13. Current vector and permanent magnet flux vector under unity power factor
operation (UPFC).
Ψ PM
In this case when φ = 0 we have the relationship:
U
U
sq
sd
Isq
= = tanδ
I
(2.74)
I
sq
Substituting the voltage equations (2.69a-b) into (2.71) and made some simplifying, we
can obtain:
I L −L −Ψ + L I = (2.75)
2
s
(
d q)cos δI PM
cosδI q s
0
27

Modeling and control modes of PM synchronous motor drives
Solving for the torque angle δ
I
:
2 2
PM PM
4 Is
( Ld Lq)
Lq
Ψ − Ψ − −
−1
δ
I
= cos [ ]
2( L − L ) I
d q s
(2.76)
only positive sign should be take into consideration.
After obtaining δ I
the amplitude of stator voltage vector can be calculated from
equation (2.70) and amplitude of stator flux vector from (2.71). The active and reactive
power and also the power factor can be obtained from equations (2.41),(2.46), (2.47).
Constant stator flux (CSF) control
As it can be see from the torque expression (2.58) for a given stator flux amplitude
the electromagnetic torque
amplitude
M
e
is a function of torque angle δ Ψ
. The stator flux linkage
Ψ
s
is kept constant of the permanent magnet flux amplitude Ψ
PM
.
I s
q − axis
Ψs
Ψ
s
δ I
δ Ψ
d − axis
Figure 2.14. Flux vector and permanent magnet flux vector under constant stator flux operation
(CSFC).
The amplitude of the stator flux linkage vector is
Ψ PM
2 2 2 2
Ψ
s
= Ψ
sd
+Ψ
sq
= ( LI
q sq
) + ( LI
d sd
+Ψ
PM
)
(2.77)
Equating
Ψ =Ψ (2.78)
s
PM
can be obtain the relationship for rotor frame currents as:
2 2
q sq d sd d PM sd
( LI ) + ( LI ) + 2LΨ I = 0
(2.79)
This condition is true if I
sd
< 0 , because expression
always positive values.
2 2
q sq
LI
d sd
( LI ) + ( ) and L , Ψ are
d
PM
28

Modeling and control modes of PM synchronous motor drives
2 2 2 2 2 2
d q s
δI d PM s
δI q s
( L − L ) I cos + 2L Ψ I cos + L I = 0 (2.80)
Solving for the torque angle δ
I
2 2 2 2 2 2
d PM d PM s
(
d q
)
q
−L Ψ ± L Ψ − I L −L L
−1
δ
I
= cos [ ]
( L L ) I
2 2
d
−
q s
(2.81)
For givenδ I
, the amplitude of stator voltage vector we can be calculated from equation
(2.70) and amplitude of stator flux vector from (2.71). The active and reactive power
and also the power factor can be obtained from equations (2.41),(2.46),(2.47),
respectively for defined speed.
Comparison study
In order to compare the control strategies and to cancel dependence of machine power,
per unit values defined as shown in Table 2.1 below have been introduced [3,9].
The value of current vector:
I
sN
= Is
Is
I
= 2I
(2.82)
b
srms( rated )
The value of voltage vector: U
s
U
s
U
sN
= = (2.83)
Ub ΩΨ
b PM
where: Ω b
= 2 π f and b
fb
is rated frequency
of the PM motor.
The value of flux vector is: Ψ
s
Ψ
s
Ψ
sN
= =
(2.84)
Ψb
ΨPM
The value of torque is:
Me
Me
M
eN
= =
(2.85)
M 3
b pbΨ
PMIb
2
The value of apparent power vector
S S
S = N
S
= 3
(2.86)
b UI
b b
2
The value of active power
P
PN
= (2.87)
Sb
The value of reactive power
Q
QN
= (2.88)
S
b
Table 2.1. Per unit values definition.
In order to compare the steady state performance characteristic of the above discussed
control strategies, for each of the control strategy some important quantities of the
machine have been plotted as a function of the torque. The PMSM parameters, which
are used for the calculations are given in Appendices.
29

Modeling and control modes of PM synchronous motor drives
The current requirement versus torque is illustrated in Fig. 2.15 for the different control
strategies. It can be seen, that up to 1 pu torque, the requirement for current is lowest for
CSF control. Highest than 1 pu torque the low current needs MTPA control requirement
lowest current for a given torque.
3
I
sN
[ pu]
2.5
2
1.5
UPF
CSF
CTA
MTPA
1
0.5
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
M [ ] eN
pu
Figure 2.15. Stator current amplitude under different control strategies versus electromagnetic
torque.
The voltage requirement versus torque for the different control strategies is illustrated in
Fig. 2.16. It can be seen, that CSF requires the highest value of stator voltage.
U
sN
[ pu]
2.5
2
1.5
1
CSF
CTA
MTPA
0.5
UPF
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
M [ ] eN
pu
Figure 2.16. Stator voltage amplitude versus electromagnetic torque under different control
strategies (at 1 pu rotor speed).
3
PN
2
[ pu]
1
UPF
CSF
CTA
MTPA
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
M [ ] eN
pu
Figure 2.17. Active power versus electromagnetic torque under different control strategies.
30

Modeling and control modes of PM synchronous motor drives
The active power requirement as a function of torque is illustrated in Fig. 2.17 for the
different control strategies. It can be seen, that all control strategies require
approximately the same value of active power for a given torque. CSF control needs
less active power in the region up to 1.3 pu torque.
The reactive power requirement as a function of torque is illustrated in Fig. 2.18. It can
be seen, that CTA control requires the highest value of active power for a given torque
and the CSF control lowest.
4
QN
[ pu]
3
2
CTA
MTPA
1
UPF
CSF
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
M [ ] eN
pu
Figure 2.18. Reactive power versus electromagnetic torque under different control strategies.
1.1
1
0.9
UPF
CSF
cosφ
0.8
MTPA
0.7
0.6
CTA
0.5
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
M [ ] eN
pu
Figure 2.19. Power factor as a function of electromagnetic torque under different control
strategies.
The power factor as a function of torque is illustrated in Fig. 2.19. It can be seen, that as
it could be expected, UPF control requires constant power factor for a given torque.
CSF control is very close to the unity power factor up to 1 pu torque.
The above analysis can be summarized as shown in Table. 2.2.
31

Modeling and control modes of PM synchronous motor drives
Table. 2.2. Summary of voltage, power, power factor requirements under control modes.
Requirement
Control
method
Voltage
Current
Power
factor
CTA middle low
low
MTP
A
low
low
middle
UPF low high 1
CSF high lowes up to 1.1
pu torque
Close to 1 up to
1 pu torque
From this comparison study it can be concluded that CSF control appears to be superior
in terms of steady state performance characteristics compared to other methods under
discussion.
32

Modeling and control modes of PM synchronous motor drives
2.3 Summary
‣ There are different forms to express the PMSM equations, but the rotor
reference frame equations are the most widely used. The simplification in rotor
dq , reference frame equations results from the disappearance of position
dependent inductances.
‣ The electromagnetic torque of the IPMSM is not only produced by the
permanent magnet flux, but also by the reluctance difference in rotor d- and q-
axes.
‣ Electromagnetic torque as cross vector product of the stator flux linkage and
current space vectors or rotor and stator flux linkages is independent of
coordinate system selected. Therefore, can be expressed in stationary ( α,
β ) or
rotated ( dq) , coordinates.
‣ For further control strategies consideration it is convenient to express the
electromagnetic torque of PMSM machine by:
• vector product of stator current and rotor flux vectors. The rotor flux
vector in PMSM machine is constant, because of the PM. Therefore, to
increase and decrease the torque, the current amplitude and the torque
angle δ
I
should be changed (see Fig. 2.20a),
• vector product of stator flux vector and rotor flux vector. Generally, the
value of the stator flux amplitude is kept constant at value of rotor flux
produced by permanent magnets. So, in this case to change the torque we
should adjust the torque angle δ Ψ
(see Fig. 2.20b).
q
β
a )
q β b)
I s
Ψ PM
d
Ψ s
Ψ
PM
d
δ I
γ m
α
δ Ψ
γ m
α
Fig. 2.20 Torque production: a) current control, b) flux control
‣ Taking into account discussion regarding static characteristic under different
control strategies it can be said that –depart from special requirements- the most
suited for general application PMSM drives is constant stator flux (CSF)
operation.
33

Voltage source PWM inverter for PMSM supply
Chapter 3 VOLTAGE SOURCE PWM INVERTER FOR PMSM SUPPLY
3.1 Introduction
The block scheme of an adjustable speed drive (ASD) commonly used in industrial
applications to supply three-phase AC motor is presented in Fig. 3.1.
Three-phase
grid
Rectifier
DC link
filter
Inverter
AC motor
Choke
Figure 3.1 Basic scheme of adjustable speed AC motor system.
An ASD is supplied from three or single phase grid. It consists of a diode rectifier, DC
link filter and an inverter. The rectifier converts supply AC voltage into DC voltage.
The DC voltage is filtered by a capacitor in the DC link. The inverter converts the DC
to an variable voltage, variable frequency AC for motor speed or (torque/current)
control.
The rectifier section of an ASD, called the front end, is responsible for generating
current harmonics into the power supply system. Therefore, to reduce the total harmonic
distortion (THD) of phase current it is necessary to add additional choke inductances.
There are generally to way how insert choke inductances (see Fig. 3.2).
a)
b)
L F
D 1
D3
D5
D 1
D3
D5
L F
L F
C F
U DC
C F
U DC
L F
D 2
D4
D6
D 2
D4
D6
L F
Fig. 3.2 Three phase diode rectifier with smoothing choke: a) at the input b) at the DC link side.
By adding a choke inductance at the input of rectifier gives the significant harmonic
reduction. Some drive manufactures are starting to include this choke inductance in the
DC link of the drive, providing the same harmonic current reduction benefit.
Regarding to power electronics standards IEEE Std 519 it is recommended that
production of harmonics should be less than 5%. It is a trend to replace the diode
34

Voltage source PWM inverter for PMSM supply
rectifier by fully controllable active rectifier [8] (see Fig. 3.3), which guaranties
following futures:
• power flow from AC/DC or DC/AC side (there is no need of break resistor ),
• significant reduction of phase current THD,
• unity power factor (phase voltage is in phase with current),
• reduction of DC link capacitor,
• controllable DC link voltage.
Three-phase
grid
Active
Rectifier
DC link
filter
Inverter
AC motor
Choke
Fig. 3.3 Modern AC/DC/AC converter topology of adjustable speed drives.
In the next part of this thesis the author will be focus on voltage source inverter.
3.2 Voltage source inverter (VSI)
The made constant DC voltage by rectifier is delivered the to the input of inverter (Fig.
3.4), which thanks to controlled transistor switches, converts this voltage to three-phase
AC voltage signal with wide range variable voltage amplitude and frequency [3].
Voltage Source Inverter
C F
T 1
D 7
T 3
D 9
T 5
D 11
U DC
O
C F
T 2
D 8
T 4
D 10
T 6
D 12
A B C
U
sAN
U
sBN
U
sCN
N
Three-phase motor windings
Figure 3.4 Basic scheme of voltage source inverter circuit.
35

Voltage source PWM inverter for PMSM supply
The one leg of inverter consists of two transistor switches. A simple transistor switch
consist of feedback diode connected in anti-parallel with transistor. Feedback diode
conducts current when the load current direction is opposite to the voltage direction.
Assuming that the power devices are ideal: when they are conducting the voltage across
them is zero and they present an open circuit in their blocking mode. Therefore, each
inverter leg can be represented as an ideal switch. Its gives possibility to connect each
of the three motor phase coils to a positive or negative voltage of the dc link ( U ).
Thus the equivalent scheme for three-phase inverter and possible eight combinations of
the switches in the inverter are shown in Fig. 3.5.
DC
U
7
= 111
U
0
= 000
U DC
1 1
S
A
SB
1
S C
U DC
S
A
0 0
SB
0
S C
U
1
= 100
A B C
U
2
= 110
A B C
U
3
= 010
U DC
1
SA
0
SB
0
S C
U DC
1 1
S
A
SB
0
S C
U DC
0
S
A
1
SB
0
S C
U
4
= 011
A B C
U
5
= 001
A B C
U
6
= 101
A B C
U DC
0
SA
1
SB
1
S C
U DC
0
S
A
0
SB
1
S C
U DC
1
S
A
0
SB
1
S C
A B C
A B C
A B C
Figure 3.5 Possible switches state in VSI.
The six positions of switches ( U1 − U6) produce an output phase voltage equal ± 1/3 or
± 2/3 of the DC voltage. The last two ( U0,
U
7) give zero output voltage. The output
phase voltages produced by inverter are shown in Fig. 3.6a. and adequate line to line
voltage calculated in bellow formula also are presented in Fig. 3.6b.
UsAB = UsAN − UsBN
(3.1a)
UsBC = UsBN − UsCN
(3.1b)
UsCA = UsCN − UsAN
(3.1c)
36

Voltage source PWM inverter for PMSM supply
a )
b)
U sAN
U sAB
2
U
3 DC
U DC
2
− U
3 DC
2π
Ωt
−U DC
2π
Ωt
2
U
3 DC
U sBN
U DC
U sBC
2
− U
3 DC
2π
Ωt
−U DC
2π
Ωt
2
U
3 DC
U sCN
U DC
U sCA
2
− U
3 DC
U1
U
2
U3
U
4
U5
U6
2π
Ωt
−U DC
a )
b)
U1
U
2
U3
U
4
U5
U6
Figure. 3.6 Three voltage waveforms generated by the inverter: a) phase voltages, b) line to line
voltages.
2π
Ωt
Form the Fourier analysis for phase voltage produced by inverter (Fig. 3.7) the
maximum amplitude of fundamental phase voltage for a given DC link voltage is given
by:
U
_ amp
2
= UDC
(3.2)
π
U out
2
U
3 DC
1
U
3 DC
1
− U
3 DC
2
− U
3 DC
2
U DC
π
2π
ωt
Figure 3.7 Inverter phase voltage generated during six step operation (solid line), corresponding
fundamental component of output voltage (dashed line) and harmonic spectrum of phase
voltage.
37

Voltage source PWM inverter for PMSM supply
The three-phase output voltage of the inverter can be described by space vector
definition as:
U
k
⎧2
⎪ U
3
⎪
= ⎨
⎪⎪⎪ 0
⎩
DC
e
π
j( k−1) 3
for k =1,2...,6.
for k=0,7.
, (3.3)
where k denotes numbers vector.
Vectors from 1-6 are called active vectors, whereas vectors 0,7 are called zero vectors
or non active vectors. The voltage space vector
U
k
in complex plane forms a regular
hexagon and divides in into six equal sectors (one sector takes 60 electrical degree) Fig.
3.8.
Im
U 4
(011)
U 3 (010)
sector3
sector4
sector2
sector5
U (110) 2
sector1
U (000) U 1
(100)
0 U (111) 7 Re
sector6
2
U
3 DC
U (001)
U (101)
5 6
Figure 3.8 Representation of the inverter states in the complex space.
In practice the real voltage source inverter has non-linear characteristic due to [19]:
• the dead-time,
• a voltage drop across the power switches,
• pulsation of the DC link voltage.
Dead time effect [17,20,27,30]
Semiconductors power switches of voltage source inverter operate not ideally. They do
not turn-on or turn-off instantaneously. Therefore it is necessary to include a protection
time to avoid a short circuit in the DC link, when two switching devices are in the same
leg (see Fig. 3.9). This time T d
is included in the control signals and it is called “dead
38

Voltage source PWM inverter for PMSM supply
time”. It guarantees safe operation of the inverter. The typical value is from 1µ s - 5µ s .
When the lowest value is for small power IGBT and is growing in respect to increasing
of IGBT power. More details about real IGBT module you can find in Appendices.
The effect of dead time can be examined from one phase of PWM inverter. The basic
configuration is shown in Fig. 3.9. Consist of upper and lower power devices T
1
and T
2
,
and reverse recovery diodes D
1
and D
2
, connected between the positive and negative
rails of power supply. The gating signals S
A
and S
Ai
come from control block. Output
voltage terminal U
0
is connected to motor phase.
S A
T 1
D 1
U DC
I sA
S Ai
T 2
D 2
LOAD
U 0
S
Dead time
Td
Figure 3.9 Circuit diagram of one inverter leg.
Fig. 3.10 shows the ideal control signals and real control signals with inserted dead time
T
d
. As can be observed the time duration of real drive signal for upper transistor is
shorted than ideal drive signal and for lower transistor is longer than ideal.
Ideal drive signals
Real drive signals
S A
S A
T d
S Ai
S Ai
T d
Figure 3.10 Gate signals control of one inverter leg.
39

Voltage source PWM inverter for PMSM supply
As a consequence when the phase current I sA
is positive, the output voltage is reduced,
and when the current I
sA
is negative the output voltage is increased (see Fig. 3.11).
U0
I
sA
> 0
ideal voltage
real voltage
U I < 0
0
sA
U DC
2
U DC
2
−U DC
2
deacresing
−U DC
2
increasing
Figure 3.11 Dead time effect on the inverter output voltage: (fat line real voltage, doted line
ideal voltage).
Voltage drop across power devices
In real voltage source inverter power switches do not conduct ideally. When they are
conducting the voltage across them is not zero and equal the voltage drop on the
conducted transistor V T
. Also in blocking mode the power switches have voltage drop
on the conducted diode V D
. More details about real IGBT module you can find in
Appendices.
The voltage drop across the power devices is dependent on the direction of the phase
current. It has influence on the output voltage, especially at low speed operation of
motor and high load current [17,20,27,30]. Fig. 3.12 shows the voltage drop influence
on the output voltage. Also shows that the output voltage is asymmetric (with offset)
and the voltage drop decreases the output voltage when the phase current is positive and
increases the output voltage when the phase current is negative.
U 0
I > 0
sA
U 0
I < 0
sA
U DC
2
V T
U DC
2
V D
V D
V T
−U DC
2
−U DC
2
Figure 3.12 Output voltage in voltage source inverter due to voltage drop across the power
devices a) for I
sA
> 0 , b) I
sA
< 0 (fat line real voltage, doted line ideal voltage).
The influence of dead time effect and voltage drop across power devices on the output
voltage from inverter is illustrated in block diagram (Fig. 3.13). The ideal reference
*
voltage components in stationary reference frame ( U α
,
_ ideal
*
U β _ ideal
) are equal real
40

Voltage source PWM inverter for PMSM supply
*
( U α
,
_ real
*
U β _ real
) and delivered to pulse width modulation (PWM) modulator block
with real non-linear inverter. As a result the output voltages ( U α _ out
, U β _ out
) are distorted
(green signals) and as consequence the phase currents ( I α _ out
, I β _ out
) in the load (red
signals) are also distorted.
*
U *
α _ideal
U α _ real
*
*
U β
U
_ideal
β _ real
PWM
Modulator
S
A
S B
S C
Inverter
U α _out
U β _out
I α _ out
I β _out
Load
Figure 3.13 Block diagram illustrated the dead time effect and voltage drop across power
devices in three phase motor supplied from non-ideal voltage source inverter.
Pulsation of the DC link voltage [19]
In practice it should be take into account that the real input dc-link voltage required for
supply VSI is not ideal. It has ripples and fluctuation, because of not ideal filtering and
disadvantages of diode rectifier. Therefore, the quality of dc-link voltage has impact on
the output voltage from inverter. If dc-link voltage will change we can observed the
changing at the output of inverter. In order to overcome this problem:
• in PWM modulator we can not assume a constant dc-link voltage and we should
measured this voltage in order to calculate the modulation index (see subchapter
3.3),
41

Voltage source PWM inverter for PMSM supply
• instate of diode rectifier will be use the active rectifier, which provided
controllable DC-link voltage.
• or used bigger capacitor in the DC-link side in order to increase possibility of
filtering.
Let us summarize, adding influence of non-linear VSI causes by:
• serious distortion in the inverter output voltage,
• distorted machine currents,
• torque pulsation,
Additionally, also causes motor instability due to the interaction between motor and the
PWM inverter, or the choice of the PWM strategy [25].
Based on simulated and experimental observation one can say that the dead time effect
is:
• more visible in low speed operation of the motor,
• may become significantly in drives where high switching frequency is required
for good dynamics performances.
In some applications such as sensor-less vector control, the inverter output voltages are
needed to calculate the rotor or stator flux vectors. Unfortunately, it is very difficult to
measure the output voltage and requires additional hardware. The most desirable
method to obtain the output voltage feedback signal is to use the reference voltages
instead. However, the relation between the output and reference voltage is nonlinear due
to the dead-time effect and voltage drop across power devices. Thus, unless the properly
dead-time and voltage drop compensation will be applied, the reference voltage can not
be used instead of the inverter output voltage. Several compensation method were
proposed to overcome this problem. One of them will be present bellow.
Compensation based on modification of reference voltage waveform [17]
The compensation process of dead time effect and voltage drop across power devices on
the inverter output voltage from is illustrated in Fig. 3.14.
42

Voltage source PWM inverter for PMSM supply
Compensation
block
IsA
IsB
Compensation of
inverter
I sC
U DC
T comp
U α _comp
U β _comp
*
U α _ideal
*
U β _ideal
*
U α _ real
*
U α _ real
*
U β _ real
PWM
Modulator
S
A
S B
S C
Inverter
*
U β _ real
U α _out
U β _out
I α _out
I β _out
Load
Figure 3.14 Block diagram illustrating the dead time and voltage drop across power devices
compensation method in three phase motor supplied from non-ideal voltage source inverter.
In order to compensate the inverter non-linearity to the ideal reference voltage
*
components in stationary reference frame ( U α
,
(
_ comp
_ ideal
*
U β _ ideal
) a compensation signal
U α
, U β _ comp
) is added. As a consequence the real reference voltage components
*
( U α
,
_ real
*
U β _ real
) are pre-distorted. Further those signals are delivered to PWM
modulator with non-ideal inverter. As a result the output voltages are not distorted in
Fig. 3.14 and thus phase currents in the load (red signals in Fig. 3.14) are almost
sinusoidal.
To calculate an average compensation voltages ( U α _ comp
, U β _ comp
), the parameters of
IGBT modules as:
• dead time T
d
,
• turn on T
ON
and turn off T
OFF
of IGBT transistors,
• and also on a voltage drop on diode V
D
and transistor
should be know.
V
T
,
43

Voltage source PWM inverter for PMSM supply
The total compensation time to compensate the non-linearity of inverter can be
calculated as:
U
T = T + T − T + T
(3.4)
dp
comp d ON OFF s
U
DC
Where Udp = VD −ton ( VD − VT )/ Ts
and t on
is conducting time of IGBT devices in one
sampling time.
The compensation voltage vector can be obtained as:
Tcomp
Ucomp = 2 UDC s ign( Is ) = 2Uth s ign( Is
)
(3.5)
Ts
where
and
2
2
sign( Is) = ( sign( IsA) + asign( IsB) + a sign( IsC))
,
3
⎧ sign( IsA) = 1 if IsA
> 0
sign( IsA)
= ⎨
⎩sign( IsA) = 0 if IsA
< 0
(3.6)
The sign function for remain phase currents are calculated similarly.
Solving equations (3.5) for real and imagine part in stationary frame, one obtains:
1
U U sign I sign I sign I
3
α _ comp
= 2
th
(2 (
sA) −0.5 (
sB
) − 0.5 (
sC
)) (3.7a)
1
U U sign I sign I
3
β _ comp
= 2
th
( (
sB
) − (
sC
))
(3.7b)
The waveform of compensation voltages in stationary frame are shown in Fig. 3.15.
44

Voltage source PWM inverter for PMSM supply
4
U
3 th
8
U
3 th
2
U
3 th
Figure 3.15 Voltage compensation components in stationary reference frame.
From the top reference α and β components.
In order to illustrate the effectiveness of the proposed compensation a simulation study
has been performed. Fig. 3.16a shows the phase current in α,
β frame and their
hodograph without compensation and Fig. 3.16b with proposed compensation method.
a)
a)
b)
b)
Figure 3.16 Nonlinearity effect of voltage source inverter on phase current of AC machine:
a) without compensation, b) with compensation.
45

Voltage source PWM inverter for PMSM supply
3.3 Space vector based pulse width modulation (PWM) methods
In voltage source inverter the transistors are controlled in a on-off fashion. In order to
obtain a suitable duty cycle for each switches the technique pulse with modulation is
used. The modulation methods [18,21,22,23,24,26,28,29,30] have the influence on:
wide range of linear operation, low content of higher harmonics in voltage and current,
low frequency harmonics, minimal number of switching to decrease switching losses in
the power components.
The most important factor in PWM mode is modulation index defined as the ratio of the
reference voltage amplitude value to the maximum voltage amplitude value during sixstep
operation (see Fig 3.17) and is given by:
M
U
ref
= (3.8)
2
U
DC
π
where the
U
DC
is the DC link voltage (for three phase six diodes rectifier is 560 V ).
The modulation index varies between 0-1 and can be divided into two regions: the
linear ( 0< M ≤ 0.907) and the nonlinear ( 0.907 < M ≤ 1) as is shown in Fig. 3.17.
Im
U3(010)
U 2(110)
End of overmodulation
region (six-step mode)
2
π
Umax
= U DC
⇒ M = 1
sector = 2
sector = 3
2
U
3 dc
U ref
End of linear region
U
max
= U DC
⇒ M = 0.907
3
U 4
(011)
sector = 4
U 0
(000)
U 7
(111)
sector = 5
θ ref
T1
U
1
Ts
sector = 1
sector = 6
Re
U (100) 1
nonlinear
(overmodulation) region
linear region
U 5 (001)
U 6 (101)
Figure 3.17 Space vector diagram of the available switching vectors.
46

Voltage source PWM inverter for PMSM supply
Linear range of operation ( 0> M

Voltage source PWM inverter for PMSM supply
Next from
U
ref
, α
ref
it is necessary to calculate the time interval for particular vectors.
U 2
U ref
120°- α
UU
0, 7
t 0
t 2
α ref
t 1
Using the low of sine it is possible to write:
U 1
120°
Figure 3.18 One sector in voltage plane.
U
ref U1 U2
= =
sin120° sin(60 °−α
) sinα
ref
ref
(3.12)
From this relations calculated value of vectors
sin(60 °−αref
) 2
U1
= Uref = Uref sin(60 °−αref
)
sin120°
3
(3.13a)
sin
ref 2
U2
= U α = U
sin120°
3
sinα
ref ref ref
(3.13b)
and respectively the normalized times are given:
3 U
1
ref
1
= = sin(60 °− αref
)
(3.14a)
2
U
U
DC
DC
t
3
U
3 U
2
ref
2
= = sinα
ref
(3.14b)
2
U
U
DC
DC
t
3
U
Putting Eq. (3.13a-b) in to Eq. (3.14a-b) the normalized value can be presented as:
2 3
t1
= Msin(60 °− α ref
)
(3.15a)
π
2 3
t2
= Msinα
ref
(3.15b)
π
or in other form:
48

Voltage source PWM inverter for PMSM supply
2 3
t2
= Msinα
ref
(3.16a)
π
3 1
t1 = M cosαref
− t2
(3.16b)
π
2
After t 1
and t 2
calculation, the remaining normalized time is reserved for zero vectors
U
0
, U7
with condition t 1
+ t 2
≤ 1.Therefore, the normalized total time for zero vectors
becomes:
t = t + t = 1 − ( t + t )
(3.17)
07 0 7 1 2
The equations (3.15a-b) for time interval of active vectors and equation (3.17) for total
time interval of zero vectors are identical for all variants of space vector modulation
(SVM) techniques.
The absence of neutral wire in star connected load provides a degree of freedom in
selecting the partitioning (zero sequence signals -ZSS) time of the two zero vectors. It is
equivalent to the freedom of injected signals in to phase signals. Therefore, it gives
different equations of t 0
and t 7
for each PWM method, but normalized duration time of
must fulfill condition in Eq. 3.17. As a consequence is only in different placement of
zero vectors U
0
, U
7
. Therefore, we can introduce the portioning factor of zero vectors,
which is defined as:
t t
k = =
t + t t
7 7
0 7 07
(3.18)
Please note that, the zero sequence signals does not change the inverter output line-toline
voltage.
From knowledge of the neutral voltage U
N 0
(see Fig. 3.4) and information what kind of
zero sequence signal (ZSS) will be injected in each phase of motor it is possible to
calculate normalized duration time of zero vectors t 0
and t 7
. In bellow Table 3.1 are
summarized different three-phase modulation techniques and remarks. Also graphical
interpretation are shown in Fig. 3.19.
49

Voltage source PWM inverter for PMSM supply
Zero Sequence Signal Time interval of zero vectors Remark
SPWM
ZSS =0
for even sectors
4
t0
= (1 − ( M cos θ )) / 2
π
for odd sectors
4
t0 = t07
−(1 −( M cos θ )) / 2
π
and t7 = t07 − t0
U
DC
/2 π
M = = = 0.785
2U
DC 4
π
complicated calculation
for time interval of
vectors
THIPWM
ZSS =sinusoidal signal
with triple harmonic
for even sectors
4 1
t0
= (1 −( M cos θ) − cos 3 θ)) / 2
π 6
for odd sectors
4 1
t0 = t07
−(1 −( M cos θ) − cos 3 θ)) / 2
π 6
and t7 = t07 − t0
U
DC
/ 3 π
M = = = 0.907
2U
DC 2 3
π
Increase the linearity of
inverter grater than 15%
of SPWM
SVPWM
ZSS=triangle signal
with triple harmonic
0.906
for all sectors
M =
t = t t t /2
time interval
7 07 0 07 t = t /2
0 07 Simple calculation for
k=0.5
Table 3.1. Variants of three-phase space vector modulation techniques.
a) b)
c)
the portioning
factor of zero vectors
Fig. 3.19 PWM techniques with various zero signal sequence shape: a) SPWM, b) THIPWM,
c) SVPWM. The upper part of figure: phase voltage U
AN
(green), pole voltage U
AO
(blue),
voltage between neutral points U
NO
(black). The lower part of figure shows the portioning
factor of zero vectors for all PWM techniques.
50

Voltage source PWM inverter for PMSM supply
Except SPWM technique all PWM method guarantee that the ZSS extends the range of
modulation index from 0.78 to 0.906, i.e. 15% greater than that obtained with standard
version of SPWM.
The duty time cycles in sector 1 for each phase can be written:
d = t + t + kt
(3.19a)
a
b
1 2 07
d = t + kt
(3.19b)
dc
2 07
= kt
(3.19c)
07
for all sectors the duty time calculations for each phase can be calculated:
sec tor1 sec tor 2 sector 3 sec tor 4 sector5 sec tor 6
⎡
⎤
⎡d 1 1 1 0 0 0 0 0 0 1 1 1
a ⎤ k k k k k k t
⎢
⎥ ⎡ 1 ⎤
⎢
d
⎥
= ⎢0 1 k 1 1 k 1 1 k 1 0 k 0 0 k 0 0 k ⎥ ⎢
t
⎥
b
2
⎢ ⎥ ⎢
⎥ ⎢ ⎥
⎢⎣d ⎥ 0 0 k 0 0 k 0 1 k 1 1 k 1 1 k 1 0 k ⎢t
⎥
c ⎦ ⎣ 07 ⎦
⎢
⎣
T
⎥
⎦
(3.20)
Depending on the location of the space vector, the basic vectors must be chosen in order
to get the minimum number of changes in the switches of the converter. The switching
sequence for each sector and suitable pulse pattern for first sector are shown in Fig.
3.20.
Sector
Three-phase Modulation
T S
/2
T S
t
0 t1
2 7
t
t
1 U0, U1, U2, U7, U2, U1,
U
0
2 U0, U3, U2, U7, U2, U3,
U
0
3 U0, U3, U4, U7, U4, U3,
U
0
4 U0, U5, U4, U7, U4, U5,
U
0
5 U0, U5, U6, U7, U6, U5,
U
0
6 U0, U1, U6, U7, U6, U1,
U
0
0 1 1 1 1 1 1 0
0 0 1 1 1 1 0 0
0 0 0 1 1 0 0 0
U0
U1
U2
U7
U
2
U1
U0
d A
d B
d C
Fig. 3.20 Switching sequence for three-phase PWM techniques (on the left ) and pulse pattern
of three-phase vector modulator in sector 1 (on the right).
51

Voltage source PWM inverter for PMSM supply
3.4 Summary
The general conclusions from this chapters can be summarized as follows:
‣ to supply the voltage source inverter can be used diode rectifier or active
rectifier with IGBT transistors,
‣ supplied PMSM machine from VSI do not required mechanical comutator. It is
thanks to electronic commutation. This overcome the problem with brushes and
periodical service,
‣ The voltage source inverter is non-linear power amplifier in respect to:
o dead time effect,
o voltage drop across the power devices,
o DC link voltage pulsation.
‣ Without appropriate dead-time and voltage drop compensation, the sensorless
operation in low speed range is not possible.
‣ The quality of DC link voltage has influence on proper operation of AC drive,
‣ Sinusoidal modulation technique (SPWM) guarantees the inverter output voltage
U
DC
amplitude from 0 until V. This correspond to changes modulation index
2
from 0-0.785.
‣ Modulation techniques with zero sequence signals not equal zero guarantees the
U
DC
inverter output voltage amplitude from 0 until V (it is 15% grater than
3
SPWM). This correspond to changes modulation index from 0-0.907.
‣ In this work the PWM modulator with zero sequence signal of triple harmonics
(SVPWM) will be used. Mainly because of simple calculation of zero vector
duration and placement.
52

Control methods of PM Synchronous motor
Chapter 4 CONTROL METHODS OF PM SYNCHRONOUS MOTOR
4.1 Introduction
The basic block scheme of adjustable speed drive with control block for PMSM is
presented in Fig. 4.2. It consists of two parts: power (fat line) and control part
employed microprocessor (thin line). The first one previously have been explained in
chapter 3. The second one will be described bellow.
Figure 4.2. The basic block scheme of PMSM drive supplied voltage source inverter.
The main task of control block is follow demand reference speed by motor and provide
proper operation in static (insight of the limits) and dynamic states without any
instability. This is ensured through suitable generated gate signals for the IGBT
transistor inside of the inverter. Therefore, to make good decision how to control power
transistors in the inverter, the following feedback signals are measured and used:
• DC link voltage,
• motor phase currents,
• speed or position of the rotor.
This significantly improve dynamic behavior of the system (good performance of the
torque and speed response, very fast dynamics response with fully controllable torque in
wide speed range).
The scalar control for PMSM without damper winding (squire cage) is not simple as for
induction motor [39,40]. It requires additional stabilization loop, which can be provide
by feedback loop from: rotor velocity perturbation, active power or DC-link current
perturbation [9].
The vector control method will be described bellow.
53

Control methods of PM Synchronous motor
4.2 Field oriented control (FOC)
During many years a DC motor has been mostly used. Because of simple control
method, which based on fact that flux and torque can be controlled separately using
current control loop with PI controllers. However the weak point of this drive was DC
motor, which could not worked in aggressive or volatile environment and required
cyclical maintenance. This disadvantages has been eliminated, when instead of DC
machine a three phase PMSM motor were used.
In searching new control method for induction machine in 1971 was developed vector
control method known as field oriented control (FOC) [31,38,49]. This method allows
control the flux and the torque in the AC machine in similar way as for DC motor. It
was achieved by transform current vector in stationary reference frame ( α,
β ) into new
coordinate system ( dq) , with respect to rotor (magnet) flux vector. So the flux
produced by permanent magnet is frozen to the direct axis of the rotor (see Fig. 4.4).
q−axis
β −axis
I s
I sq
δ I
γ I
I sd
γ m
Ω s
Ψ PM
d−axis
α −axis
Figure 4.4 Vector diagram illustrated the principle of FOC.
Further, stator current vector can be split into two current components: flux current I sd
and torque producing current I sq
. In analogy to separate commutator motor, the flux
current components corresponds to excitation current and torque-producing current
corresponds to the armature current. Therefore, the goal of the control system is to
reference the
I
sd_
ref
,
sq_
ref
I stator current components in respect to requirement of
references torque and flux. The flux and torque producing stator current references are
obtained on the output of the reference current generation block (see Fig. 4.5).
54

Control methods of PM Synchronous motor
Reference Current
Calculation
Reference Current
Calculation
reference
torque
M e _ ref
FG1
FG2
I sd _ ref
reference
d-axis current
I sq _ ref
reference
q-axis current
reference
torque
M e _ ref
FG3
FG3
δ I
I s
γ I
Figure 4.5. Reference current generator block for FOC technique
a) in cartesian form, b) in polar form.
γ m
Function generation FG1 gives the relationship between the torque and the direct axis
stator current component
I
sd_
ref
, and function generator FG2 gives the relationship
between the torque and the quadrature axis stator current
I
sq_
ref
. His graphical
illustration in Fig. 4.6 are presented.
0.5
0
CTA
I [ . ]
sdN
pu
0.5
1
CSF
MTPA
1.5
2
UPF
2.5
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
M
eN
[ pu . ]
2.5
I [ . ]
sqN
pu
2
1.5
1
CSF
CTA
MTPA
0.5
UPF
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
M
eN
[ pu . ]
Figure 4.6. Generated current components I
sd
and I sq
dependent on required electromagnetic
torque under difference control strategies. CTA-constant torque angle control, MTPA-maximum
torque per ampere control, CSF-constant stator flux control, UPF – unity power factor control.
55

Control methods of PM Synchronous motor
For many years the CTA control ( I
_
= 0 ) method has been a popular technique for a
sd
ref
long time because of simple control. This method was dedicated for surface permanent
magnet synchronous motor (SPMSM), where the magnetic saliency does not exist
[126]. So the maximum torque per ampere is obtained, when stator current vector is
shifted in respect to rotor flux vector 90 degree. However, in IPMSM, maximum torque
per ampere is obtained with torque angle more than 90 degree. This is because of
existence of reluctance torque component due to magnetic saliency (see subchapter
2.2.2). Therefore, the I
sd_
ref
should be negative value [44].
The main question is how or in which manner produce the reference currents in d-q
frame. Its leads to many realization of current control structure. Among them generally
we can distinguish two structures of current control loop. One of them is hysteresis
based control (Fig.4.7a) [3,52] and the second one is PI based current controllers
(Fig.4.7b).
Hysteresis based current control has following disadvantage such as [3]:
• measurement of three phase currents are required,
• three independent hysteresis current controllers are required,
• variable switching frequency is achieved,
• fast sampling time is required.
All this listed above disadvantage can be eliminate, when the PI current control are
used. This structure are mostly used in industrial application (Fig. 4.7b).
a)
U DC
M e _ ref
Fig. 4.5
I d _ ref
I q _ ref
d,q/ABC
I A _ ref
I B _ ref
I C _ ref
−
−
S A
S B
S C
Inverter
current
feedback
−
I A current
I
sensors
B
I C
rotor
position
sensor
PMSM
56

Control methods of PM Synchronous motor
b)
U DC
M e _ ref
Fig. 4.5
I sd _ ref
I sq _ ref
-
e Isd
-
e Isq
PI
PI
U sd
U sq
Reference
VoltageVector
Calculation
U s _ ref
ϕ Us _ ref
Space
Vector
Modulator
S A
S B
S C
Inverter
γ s
I sd
I sq
dq/ABC
current
sensors
I s
γ m
rotor
position
sensor
PMSM
Figure 4.7 Vector control structure for PM synchronous motor with: a) hysteresis current
control, b) synchronous PI current control
4.3 Direct torque control (DTC)
The name direct torque control is deliver by the fact that, on the basis of the errors
between the reference and the estimated values of torque and flux, it is possible to
directly control the inverter states without inner current control loop as for FOC
[32,34,35,50] and [57-66].
The basic idea of this control rely on stator voltage vector equation of AC motor.
U
dΨ
s
s
= RsIs
+ (4.1)
dt
Making the assumption that ohmic voltage drop on the stator resistor can be neglected
the equation for stator flux vector takes the form:
Ψ =∫ ( U ) dt
(4.2)
s
s
It can be said that the stator voltage vector has directly influence on control stator flux
vector. Using a three phase voltage source inverter to supply the AC motor, there are six
non-zero vectors and two zero voltage vectors. The active vectors change the amplitude
and position of stator flux vector, while the zero vectors stop the stator flux vector as
shown in Fig. 4.8.
57

Control methods of PM Synchronous motor
y
moment with active
forward vector
stops with zero
vector
Ψ s
moment with active
backward vector
γ s
δ Ψ
Ψ =Ψ
r
PM
x
rotates
continuously
α
stator
Figure. 4.7 Stator flux vector Ψ
s
movement relative to rotor flux vector Ψ
r
=ΨPM
under the
influence of active and zero inverter voltage vectors.
Therefore, it is possible control the torque angle δ Ψ
across control stator flux vector
position in respect to rotor flux vector produced by permanent magnet Ψ
r
=Ψ
PM
Ψ
s
, what
further allows to have impact on control the electromagnetic torque in accordance with
following formula:
M
e
3 Ψ sin
2 s
ΨPM
δΨ Ψs
( Lq
−Ld)sin2δ
Ψ
= pb[ − ]
(4.3)
2 L
2L L
d d q
Generally the DTC technique operate at constant stator flux amplitude
Ψ
s
, what
correspond to CSF operation, because of simple reference stator flux amplitude equal
nominal value of permanent magnet. For DTC technique can be also apply all control
strategies discussed in Chapter 2.
Using the block generator for reference stator flux amplitude and electromagnetic
torque as is shown in Fig 4.8. it is possible to draw the relationship between required
Reference
Flux
Calculation
reference
flux
Ψ s _ ref
M
e_
ref
reference
torque
reference
torque
M e _ ref
Figure 4.8. Reference flux and torque generator block for DTC technique.
reference flux and torque. His graphical illustration in Fig. 4.9 is presented.
58

Control methods of PM Synchronous motor
Ψ
sN
[ pu]
2
1.5
1
CTA
MTPA
CSF
0.5
UPF
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
M
eN
[ pu]
Figure 4.9. Generated stator flux amplitude Ψ s _ ref dependent on required electromagnetic
torque under difference control strategies. CTA-constant torque angle control, MTPA-maximum
torque per ampere control, CSF-constant stator flux control, UPF – unity power factor control.
The basic structure of direct flux and torque control voltage-sourced PWM inverter-fed
permanent magnet synchronous motor is shown in Fig.4.10.[3,36,37,41,45]
Flux
hysteresis
U DC
Ψ s _ ref
M e _ ref
−
−
H Ψ s
H m
d Ψ s
d M e
Switching
Table
Torque
sector γ
histeresis s ( N )
S A
S B
S C
Inverter
Ψ s
M e
Flux and
Torque
Estimation
I s
γ m
PMSM
Figure. 4.10. Block diagram of switching table based direct torque control ST-DTC.
The command stator flux amplitude
Ψ
s _ ref
and electromagnetic torque
e_
ref
M values
are compared with the actual
Ψ s
and M e
values, in hysteresis flux and torque
controllers, respectively. The flux and the torque controllers are a two-level
comparators.
The digitized outputs signals of the flux controllers are defined as:
d = 1 (increase flux) for Ψ > Ψ
_
+ (4.4a)
Ψ s
s s ref
H Ψ
d = 0 (decrease flux) for Ψ < Ψ
_
− (4.4b)
Ψ s
s s ref
H Ψ
59

Control methods of PM Synchronous motor
and those of the torque controller as
where
d = 1(increase torque) for M > M
_
+ H
(4.5a)
M e
e e ref m
d = 0 (decrease torque) for M < M
_
− H
(4.5b)
M e
e e ref m
H
m
and H Ψ
are hysteresis bands for torque and flux, respectively.
The digitized variables d Ψ
, d and the stator flux position sector γ M
s( N )
information
s
e
create a digital word, which select appropriate voltage vector from the switching table.
Next, from the selection table the proper voltage vectors are selected and the gate pulses
SA, SB,
S
C
to control the power switches in the inverter are generated.
The circular stator flux vector trajectory can be divided into six symmetrical sectors
(according to the non zero voltage vectors), which are defined as (see Fig. 4.11):
sector 1 − 30°≤ γ s
< 30°
V 3 (010)
V 2 (110)
sector 2 30°< γ s
≤ 90°
sector 3 90°< γ s
≤ 150°
sector 4 150°< γ s
≤ 210°
sector 5 210°< γ s
≤ 270°
sector 6 270°< γ s

Control methods of PM Synchronous motor
β
Ψ
s
↓
M
e
↑
Ψ
s
↑
M
e
↑
Ψ
s
↓
Ψs
↑
Ψ s
γ
s γ
s
α
Ψ
s
↓
M
e
↓
Ψ
s
↑
M
e
↓
Figure. 4.12 Voltage vector effects in sector 1 on stator flux and torque.
The presented rule for first sector can be extended for other sectors, what further help
construct the switching Tables 4.1 and 4.2 as below [56]
Flux Torque sector 1 sector 2 sector 3 sector 4 sector 5 sector 6
d ψs =1
d me. =1 V
2
V
3
V
4
V
5
V
6
V
1
d me =0 V
6
V
1
V
2
V
3
V
4
V
5
d ψs =0
d me. =1 V
3
V
4
V
5
V
6
V
1
V
2
d me =0 V
5
V
6
V
1
V
2
V
3
V
4
Table 4.1. Switching table for DTC with active vectors.
Flux Torque sector 1 sector 2 sector 3 sector 4 sector 5 sector 6
d ψs =1
d me. =1 V
2
V
3
V
4
V
5
V
6
V
1
d me =0 V
7
V
0
V
7
V
0
V
7
V
0
d ψs =0
d me. =1 V
3
V
4
V
5
V
6
V
1
V
2
d me =0 V
0
V
7
V
0
V
7
V
0
V
7
Table 4.2. Switching table for DTC with zero and active vectors.
Tables 4.1 and 4.2 represent the eight and six voltage-vectors switching tables.
61

Control methods of PM Synchronous motor
The DTC has lesser parameter dependence and fast torque when compare with the
torque control via PWM current control.
Among the well-know advantages of the DTC scheme are the following:
• Simple control,
• Excellent torque dynamics,
• Absence of coordinate transformations,
• Absence of separate voltage modulation block,
• Absence of voltage decoupling circuits,
• There are no current control loops, hence, the current is not regulated directly,
• Stator flux vector and torque estimation is required.
Among the well-know disadvantages of the DTC scheme are the following:
• variable switching frequency (difficulties of LC input EMI filter design),
• high sampling time is required (fast microprocessor and A/D converter ),
• inverter switching frequency depending on: flux and torque hysteresis bands,
machine parameters, sampling frequency,
• violence of polarity consistency rules (huge voltage stress for IGBT transistor),
• current and torque distortion caused by sector changes,
• start and low speed operation problems,
• high sampling frequency needed for digital implementation of hysteresis
comparators,
• high noisy level,
• high current and torque ripple.
Many modifications of the basic switching table based direct torque control (ST-DTC)
scheme at improving starting, overload condition, very low speed operation, torque
ripple reduction, variable switching frequency functioning, and noise level attenuation
have been proposed during last decade.
In the last five years, many researcher have been carried out to try solve the above
mentioned problems of ST-DTC scheme. Therefore, the following solutions have been
developed in order to eliminated mentioned before problems:
62

Control methods of PM Synchronous motor
• Use of improved switching table,
• Use of comparators with and without hysteresis, at two or three levels,
• Use of multi-level inverter,
• Introduction of fuzzy or neuro-fuzzy techniques,
• Use of sophisticated flux estimators to improve the low speed behavior,
• Implementation of DTC schemes with constant switching frequency operation
In multi-level inverter there will be more voltage vectors available to control the flux
and torque. Therefore, a smoother torque can be expected. However, more power
switches are needed to achieved a lower ripple, which will increase the system cost and
complexity.
All this contributions allow the DTC performance to be improved, but at the same time
they lead to more complex schemes. As expected, conventional DTC is growing in
field-oriented control area and the so-called improved DTC with space vector
modulation (SVM). Let us call it DTC-SVM. This control concept will be deeply
discussed in the next Chapter.
63

Control methods of PM Synchronous motor
4.4 Summary
‣ In FOC drive flux linkage and electromagnetic torque are controlled indirectly
and independently by PI controllers with space vector modulator (SVM). In this
control concept current control loop is required,
‣ In DTC drive, flux linkage and electromagnetic torque are controlled directly
and independently by hysteresis controllers and selection of optimum inverter
switching modes. In this control concept flux and torque control loop is required,
‣ In DTC all switch changes of the inverter are based on the electromagnetic state
of the motor.
‣ The DTC technique is different from the traditional methods of controlling
torque, where the current controllers in the rotor reference frame are used. It is
completely different control concept (approach) from FOC. The new control
technology was characterized by simplicity, good performance and robustness,
because of bang-bang hysteresis control. Using DTC it is possible to obtain a
good dynamic control of the torque without current controllers and any
mechanical transducers on the machine shaft. Moreover, in this control structure
the PWM modulator is not required. Its is occupied by variable switching
frequency.
‣ The flux weakening control becomes easier because stator flux linkage can be
controlled directly in the DTC system of PMSM.
64

Direct Torque Control with Space Vector Modulation
Chapter 5
DIRECT TORQUE CONTROL WITH SPACE VECTOR MODULATION
(DTC-SVM)
5.1 Introduction
The DTC-SVM greatly improves torque and flux performance by:
• Achieved fixed switching frequency,
• Reduced torque and flux ripples.
The main idea of DTC-SVM is based on analyze of the torque equation
M
Assuming that the Ld = Lq = Ls
e
3 Ψ sin
2 s
ΨPM
δΨ Ψs
( Lq
−Ld)sin2δ
Ψ
= pb[ − ]
(5.1)
2 L
2L L
d d q
M
e
3 Ψs
ΨPM
sinδ = pb[ Ψ
]
(5.2)
2 L
From equation (5.2) we can see that for constant stator flux amplitude
by permanent magnet
Ψ
PM
s
Ψ
s
and flux produced
, the electromagnetic torque can be changed by control of the
torque angle δ Ψ
. This is the angle between the stator and rotor flux linkage, when the stator
resistance is neglected. The torque angle, in turn, can be changed by changing position of the
stator flux vector θ Ψ s
in respect to PM vector using the actual voltage vector supplied by
PWM inverter.
In the steady state, δ Ψ
is constant and corresponds to a load torque, whereas stator and rotor
flux rotate at synchronous speed. In transient operation, δ Ψ
varies and the stator and rotor
flux rotate at different speeds (Fig. 5.1).
65

Direct Torque Control with Space Vector Modulation
q − axis
β
δ Ψ
Ψ s _ ref
∆δ Ψ
θ Ψ
s
Ψ s
Ψ PM
d − axis
θ r
Figure 5.1. Space vector diagram illustrating torque control conditions.
α
5.2 Cascade structure of DTC–SVM scheme
The structure of proposed control scheme is shown in the Fig. 5.2. [11,33,42,48,51,53,54]
U DC
Ψ s _ ref
U
s α _ ref
S A
S B
M e _ ref
e M
∆δ Ψ
U
s β _ ref
S C
θ Ψs Ψs
I s
M e
I s
γ m
Figure 5.2. Cascade structure of DTC-SVM scheme.
The error between reference and measured torque can be expressed as:
Ψ Ψ sin( δ +∆δ ) Ψ Ψ sinδ
eM = M − M = p
− (5.3)
3 s_
ref PM Ψ Ψ
_
[
s PM Ψ
e ref e b
]
2
Ls
Ls
From equation (5.3) we can see that the relation between torque error and increment of load
angel
∆δ Ψ
is nonlinear. Therefore, we used PI controller which generates the load angel
66

Direct Torque Control with Space Vector Modulation
increment required to minimize the instantaneous error between reference
M and actual
e_
ref
M
e
torque.
In control scheme of Fig. 5.2 the torque error signal e M
is delivered to the PI controller,
which determines the increment of torque angle
∆ δ Ψ
. Based on this signal and reference
amplitude of stator flux
Ψ , the reference voltage vector in stator coordinates α,
β is
s _ ref
calculated. The calculation block of reference voltage vector also uses information about the
actual stator flux vector (amplitude
Ψ
s
and position θ Ψ s
) as well as measured current vector
I
s
. The reference stator voltage vector is delivered to space vector modulator (SVM), which
generates the switching signals
S , S , S for power transistors of inverter.
A
B
C
The calculation block of reference voltage vector is shown in Fig. 5.3.
Ψ s _ ref
RsI s α
Ψ
s α _ ref
∆Ψ sα
U α
_
s
_ ref
∆δ Ψ
−
Ψ
s β _ ref
∆Ψ sβ
−
U
s β ref
θ Ψs
Ψ sα
Ψ sβ
Rs
Is
β
Figure 5.3. Calculation block of reference voltage vector.
Based on
∆ signal, reference of stator flux amplitude Ψ
s _ ref
and measured stator flux
δ Ψ
vector position
θ Ψ s
(Fig. 5.3.), the reference flux components
Ψ Ψ in stator
s α _ ref
,
sβ
_ ref
coordinate system are calculated as:
Ψ = Ψ cos( θ +∆δ
)
sα
_ ref s_
ref Ψs
Ψ = Ψ sin( θ +∆δ
)
sβ
_ ref s _ ref Ψs
Ψ
Ψ
(5.4)
Pleas note that for constant flux operation region the reference value of stator flux amplitude
Ψ is equal flux amplitude of permanent magnet Ψ
PM
.
s _ ref
67

Direct Torque Control with Space Vector Modulation
The references of stator flux components (see Fig. 5.3) are compared with estimated value:
Ψ
Ψ
= Ψ
cos θ ,
sα
s Ψs
= Ψ
sin θ ,
sβ
s Ψs
The command voltage can be calculated from flux errors in
(5.5)
α, β coordinate system as
follows:
U
U
∆Ψ
= + R I
sα
sα
_ ref s sα
Ts
∆Ψ
= + R I
sβ
sβ
_ ref s sβ
Ts
(5.6)
Where: T
s
is sampling time,
∆Ψ = Ψ − Ψ , ∆Ψ =Ψ
_
−Ψ .
sα sα_
ref sα
sβ sβ ref sβ
The presented bellow design methodology for flux and torque control loops based on the
approach presented in literature [11,43].
5.2.1 Digital flux control loop
The flux control loop is based on the voltage equations of PMSM machine in stator
coordinates.
U
U
dΨ
dt
sα
sα
= RsIsα
+ (5.7a)
dΨ
sβ
sβ
= RsIsβ
+ (5.7b)
dt
Using Laplace transformation the above equations can be written as:
U = sΨ + R I
(5.8a)
sα sα s sα
U = sΨ + R I
(5.8b)
sβ sβ s sβ
It corresponds to flux model of PMSM machine in α,
β system presented in Fig. 5.4.
68

Direct Torque Control with Space Vector Modulation
R I α
U
s α
U
s β
s
−
s
1
s
1
s
−
R I β
s
s
Ψ sα
Ψ sβ
Figure 5.4. Flux model of PMSM in stator coordinates.
In order to control the flux components in α,
β frame the bellow control structure can be
applied.
Ψ
s α _ ref
Ψ
s β _ ref
Flux control
part
−
∆Ψ sα
−
∆Ψ sβ
P block
P=
1/ Ts
P block
P=
1/ Ts
∆Ψ
s α
T s
∆Ψ
s β
T s
RI
s s α
U
s α _ ref
U
s β _ ref
RI
s s β
RI
s s α
−
1
s
1
s
−
RI β
s
s
Flux PMSM
model
Ψ
sα
Ψ sβ
Ψ sα
Ψ sβ
Figure 5.5. Flux control loop with two P controller in α,
β reference frame.
Pleas note that regarding to Fig. 5.1 the following rules are keeping:
Ψ
sα _ ref
= Ψ
s _ ref
cos( θΨs
+∆ δ )
(5.9a)
Ψ
sβ _ ref
= Ψ
s _ ref
sin( θΨs
+∆ δ )
(5.9b)
and
sα
s_
refcosθ Ψ s
Ψ = Ψ (5.10a)
sβ
s_
refsinθ Ψ s
Ψ = Ψ (5.10b)
In order to find the formula for tuning the P controllers in the flux loop, the following
assumptions should be made:
69

Direct Torque Control with Space Vector Modulation
• increment of torque angle ∆δ Ψ
coming from torque control loop (see Fig.5.2.) is equal
zero. It means that the torque is not produced,
• stator flux vector position θ Ψ s
and rotor flux vector position θ r
are equal zero. It
corresponds to situation, where the those two flux vectors lie along the α axis.
In this special case the reference stator flux amplitude Ψ
_
=Ψ
α _
can be controlled
s ref s ref
trough the reference stator voltage component U
α _
= U
_
, when the voltage drop on the
s ref s ref
stator resistances in α,
β axes are neglected (see Fig. 5.5). Therefore, the simplified flux
control loop can be shown in Fig. 5.6.
Ψ =Ψ α
s _ ref s _ ref
s sα
Ψ =Ψ
−
∆Ψ sα
controller
P
U
sα _ ref
= Us _ ref
U β
=
s _ ref
0
PMSM
Ψ sα
Ψ sβ
Continuous s-domain
Figure 5.6. Simplified flux control loop in α,
β coordinates.
Simplified flux control loop in s domain is shown in Fig. 5.7, where CΨ ( s)
is a transfer
function of the P controller given by:
CΨ() s = KpΨ
(5.11)
The transfer function between stator flux amplitude Ψ
s
=Ψ sα and stator voltage amplitude
U
s
can be expressed as:
G
Ψ
Ψ
s 1
() s = = (5.12)
U s
s
70

Direct Torque Control with Space Vector Modulation
Control Plant
P controller
Ψ s_ref
C
Ψ
() s
U s
G
Ψ
() s
Ψ s
Figure 5.7. Block diagram of flux controller in s domain.
Hence the transfer function of the closed stator flux amplitude control loop is obtained as:
G
Ψ _ closed
Ψ
s
() s
_ ref CΨ() s GΨ()
s
() s = =
Ψ () s 1 + C () s G () s
s
Ψ
Ψ
(5.12)
Substituting transfer function for CΨ ( s)
and GΨ ( s)
one becomes:
⎛1⎞ ⎛1⎞
KpΨ⎜ ⎟ KpΨ⎜ ⎟
s s K
p
GΨ _ closed
() s =
⎝ ⎠
=
⎝ ⎠
=
⎛1
⎞ s+ KpΨ
s+
K
1+ K
pΨ
⎜ ⎟
⎝s
⎠ s
Ψ
pΨ
(5.13)
Discrete design
The transfer function for P controller in discrete system is expressed as:
CΨ( z) = KpΨ
(5.14)
Ψ
s_ ref( z)
C
Ψ
( z)
U
= sα U
s
Dz ( )
z −1
GΨ
( z)
}
ZOH
1
s
Ψ
s
( z)
Figure 5.8. Block diagram of flux controller in discrete domain.
Where, CΨ ( z)
is discrete transfer function for P controller, Dz ( )
delay for voltage generation from PWM .
1
z −
is for one sampling time
71

Direct Torque Control with Space Vector Modulation
The GΨ ( z)
is discrete transfer function for voltage-flux relationship with zero hold order
(ZOH) can be calculated as:
−1
GΨ
() s z−1 1
GΨ
( z) = (1 − z ) Z[ ] = Z[ ]
(5.15)
2
s z s
Using table of Z transformation [2]. Finally, it gives
G
Ψ
( z −1)
zTs
AΨ
d
2
z ( z −1) ( z−1)
( z)
= =
(5.16)
Where A = Ψ d
T and s
T
s
is sampling time of the discrete system.
Hence, the closed loop transfer function between Ψ
s
( z)
and Ψ ( z)
is obtained as:
_ ref
s
G
Ψ _ closed
K
pΨ
Ψ
s
( z) _ ref CΨ( z) GΨ( z) D( z)
( z)
= =
Ψ ( z) 1 + C ( z) G ( z) D( z)
A
Ψd
( −1)
z z
K A
= =
K A z z K A
1+
z z
pΨ
Ψd
2
pΨ Ψd − +
pΨ Ψd
( −1)
s
Ψ
Ψ
(5.16)
The flux step response depended on poles placement of closed flux control loop. The pole
placement can be selected by setting the
K
p Ψ
.
Assuming, that CΨd = KpΨAΨd
the GΨ _ closed
( z)
expressed by equations (5.16) will take the
following form:
G
CΨd
( z) =
z − z+
C
Ψ _ closed
2
Ψd
(5.17)
The nomogram of Fig. 5.9 shows the relationship between overshoot M
p[%]
, rise time t r
and
settling time t s
in respect to
C Ψ d
.
Please not that t r
is time calculate from 10% to 90% of output signals and t s
is the time it
takes the system transient to decay +-1%.
72

Direct Torque Control with Space Vector Modulation
Figure 5.9. The relationship between overshoot, rise time and settling time versus to
flux amplitude control loop.
C Ψ d
for stator
Now from few values of C Ψ d
=[0.4046, 0.2688, 0.1720] we can choose C Ψ d
=0.2688, which
guaranties overshoot about 0% and settling time about 10 times of sampling time.
1
2
3
1
Figure 5.10. Step flux response for different to
C Ψ d
=0.2688, black line (3)
d
C Ψ d
: red line(1) C Ψ d
=0.4046, blue line (2)
C Ψ
=0.1720.
It corresponds to the transfer function of closed stator flux control loop as:
G
C
0.2688
( z)
= =
z − z+ C z −z+ 0.2688
Ψd
Ψ _ closed
2 2
Ψd
(5.18)
Using digitalized motor parameters A = Ψ d
T and chosen s
C Ψ d
value we can calculate the
parameters of P digital flux controller as:
73

Direct Torque Control with Space Vector Modulation
K
pΨ
C
A
Ψd
Ψd
= = (5.19)
Ψd
C
T
s
For example let us assume that sampling time in digital flux control loop is equal T = 200µ
s.
The gain of P controller is:
s
K
pΨ
0.2688 0.2688
= = = 1344
(5.20)
A 0.0002
Ψd
In digital control when the sampling time changes the parameters of digitalized plant control
A Ψ d
will also change. Therefore, to keep closed loop transfer function as close as possible to
G
C
0.2688
( z)
= =
, the gain of P flux controller should also be
z − z+ C z −z+ 0.2688
Ψd
Ψ _ closed
2 2
Ψd
changed (see Table 5.1.).
Keeping constant transfer function GΨ _ closed
( z)
the flux step response for different sampling
times
T
s
= 50µ s , 100µ s , 200µ s , 400µ s , which correspond to switching frequency f s
=
20kHz, 10kHz, 5kHz, 2.5kHz are presented in Fig. 5.11.
Figure 5.11. Flux tracking performance for different sampling times T
s
= 50µ s (blue line -1), 100µ
s
(green line -2), 200µ s (red line -3), 400µ s (light blue line -4).
74

Direct Torque Control with Space Vector Modulation
From Fig. 5.11 it can observed that overshoot is around 0% and the settling time took 10
times of microprocessor sampling time. So, it is possible control the flux amplitude in 10
samples
The settings of P flux controller for different sampling time T
s
= 50 µ s,100 µ s,
200 µ s,400µ
s
are summarized in Table 5.1.
75

Direct Torque Control with Space Vector Modulation
76

Direct Torque Control with Space Vector Modulation
The behavior of the flux control loop was tested using SABER simulation package. The
model created in SABER takes into account the whole control system, which include
real models of inverter and permanent magnet synchronous motor.
The flux step response is shown in Fig. 5.12., when parameters of P flux controller
designed for sampling time T s
= 200µ s were used for control plant for different
sampling times T s
= 50µ s , 100µ s , 200µ s , 400µ s , which correspond to switching
frequency f
s
= 20kHz, 10kHz, 5kHz, 2.5kHz.
Figure 5.12. Flux tracking performance for different sampling time T
s
= 50µ s , 100µ s ,
200µ s , 400µ s ( f
s
= 20kHz, 10kHz, 5kHz,2.5kHz) using gain of P controller designed
for T
s
= 200µ s ( f
s
= 5kHz).
After modification of P flux controller gain according to Table 5.1 it is possible to
achieve better results as shown in Fig. 5.13, what confirms proper flux tracking
performance in steady and dynamics state.
77

Direct Torque Control with Space Vector Modulation
Figure 5.13. Flux tracking performance for different sampling time T s
= 50µ s , 100µ s ,
200µ s , 400µ s ( f
s
= 20kHz, 10kHz, 5kHz,2.5kHz) using designed gain of P controller
calculated individually (see table 5.1.)
The simulation results in SABER package for Ts=200us is presented in Fig. 5.14.
Figure 5.14. Simulated (SABER) flux tracking performance for step change from 70% - 100%
of nominal flux.
As we can observed from Fig. 5.14 that overshoot is around M
p
= 0% and settling time
took about 10 sampling time of microprocessor, what proved the design procedure of P
digital flux controller gain.
78

Direct Torque Control with Space Vector Modulation
5.2.2 Digital torque control loop
The considered torque control loop is shown in Fig. 5.15.
Ψ s _ ref
M e _ ref
−
Torque
controller
PI
∆δ Ψ
Reference
flux generator
In stator frame
Ψ
s α _ ref
Ψ
s β _ ref
−
∆Ψ sα
−
∆Ψ sβ
Flux
controllers
In stator frame
U
s α _ ref
U
s β _ ref
M e
θ Ψs
Ψ sα
Ψ sβ
Figure 5.15. Torque control loop with PI controller.
Based on the equation (5.9a-b and 5.10a-b) the stator flux errors in α,
β coordinates can
be calculated as:
∆Ψ
α
= Ψ cos( θ + ∆δ ) − Ψ cos θ = Ψ [cos( θ + ∆δ ) − cos θ ]
s s_ ref Ψs Ψ s_ ref Ψs s_
ref Ψs Ψ Ψs
(5.21a)
∆Ψ
β
= Ψ sin( θ +∆δ ) − Ψ sin θ = Ψ [sin( θ +∆δ ) − sin θ ]
s s_ ref Ψs Ψ s_ ref Ψs s_
ref Ψs Ψ Ψs
(5.21b)
Assuming that for small changes of ∆δ Ψ
the cos ∆δ Ψ
≅ 1and sin ∆δ Ψ
≅∆ δ
Ψ
, the equations
(5.21a) and (5.21b) are given by:
∆Ψ
α _
= − Ψ
_
∆ δ sinθ
(5.22a)
s ref s ref Ψ Ψs
∆Ψ
β _
= Ψ
_
∆ δ cosθ
(5.22b)
s ref s ref Ψ Ψs
In order to design the PI torque controller the following assumption are made:
• stator flux vector position θ Ψ s
and rotor flux vector position θ
r
are equal zero. It
correspond to situation, when those two flux vectors lie along theα axis,
• the reference stator flux amplitude is equal value of permanent magnet flux
Ψ
s _ ref
=Ψ
PM
,
• stator resistance is neglected.
Therefore, the error stator fluxes in α,
β coordinates are calculated as:
∆Ψ = , (5.23a)
sα
_ ref
0
79

Direct Torque Control with Space Vector Modulation
∆Ψ = Ψ ∆ . (5.23b)
sβ
_ ref PM
δ Ψ
Stator voltage components using equations (2.23a-b) can be expressed as:
U α
= (5.24a)
s _ ref
0
U
= ∆Ψ K
(5.24b)
sβ _ ref sβ _ ref pΨ
And further because of U
β _
= U
_
s ref s ref
U
=Ψ ∆ δ K
(5.25)
s PM pΨ
So, the transfer function between stator voltage amplitude
angle ∆ δ Ψ
can be written as:
U
s
G
δ
() s = = K
ΨΨ
∆δ
M p PM
Ψ
U
s
and increment of torque
(5.26)
Where
K
p Ψ
is the gain of stator flux P controller.
For example for sampling time T = 200µ
s, calculated
s
0.2688
K = pΨ 1344
T
= (see
s
Table5.1.) and nominal value of Ψ
PM
= 0.264Wb
the calculated GMδ () s = 354,82V / rad .
The obtained transfer function between electromagnetic torque
amplitude
U
s
is (see equation 5.72):
M
e
and stator voltage
G
M
M () s A s
() s = =
U s s + B s+ C
(5.27)
e
M
2
s
()
M M
Where
A
M
2
3pbΨ
PM
RsΨ
3 Ψ
PM
s
ΨPM pb
= and BM
= CM
=
2L
Ψ L
2JL
s
Using the motor parameters (see Appendices), one obtains:
A = 198 and B = 115.3 C = 9065
M
M
M
s
s
s
Continuous s-domain
The torque control loop is shown in Fig. 5.16, where C ( s ) is a transfer function of the
PI controller given by [105]:
M
80

Direct Torque Control with Space Vector Modulation
C
M
⎛ K ⎞
iM
KpM
s+
⎜
K ⎟
pM
() s =
⎝ ⎠
(5.28)
s
where
K
iM
=
K
T
pM
iM
M e _ ref
() ∆δ U = sβ Us
CM s Ψ
G ()
GM
() s
Mδ
s
M e
Figure 5.16. Block diagram of torque control loop represented in s-domain.
Hence, the transfer function between M
_
() s and M ( s ) is obtained as:
e
ref
e
G
M _ closed
() s
M () s C () s G () s G () s
e_
ref M M Mδ
= = (5.29)
M
e() s 1 + CM() z GM() s GMδ
() s
Substituting transfer function for CM
( s)
and GM
( s ) equation (5.29) becomes:
KiM
K Ψ p PMAMKpM( s +
Ψ
)
K
pM
GM
_ closed()
s = =
2
s + ( B + K A ) s+ C + K A
M pM M M iM M
(5.30)
Discrete design
The transfer function (equations 5.28) for PI controller in discrete system using
z −1
backward difference method for discretization process ( s = ) [2] is expressed as:
Tz
s
K
pM
( z − )
KpM
+ KiM
CM( z) = ( KpM + KiM)
( z −1)
(5.31)
where:
K
iM
K
pM
= Ts
,
s
TiM
T - sampling time, C ( z ) is the discrete transfer function of
M
torque PI controller, Dz ( )
1
z −
is one sampling time delay for voltage generation from
PWM, and
81

Direct Torque Control with Space Vector Modulation
D( z)
} G ( ) M
z
M
e_ ref( z)
C ( ) M
z
∆δ Ψ
G
Mδ
( z)
U s
z −1
AM
s
ZOH
2
s + B s+
C
M
M
M ( ) e
z
Figure 5.17. Block diagram of torque control loop in discrete domain.
G ( z)
= K Ψ (5.32)
Mδ
pΨ
PM
is discrete transfer function for relation between stator voltage amplitude
increment of torque angle ∆ δ Ψ
(see Fig. 5.17)
U s
and
The discrete transfer function G ( z ) for voltage-torque relationship with zero order
hold (ZOH) can be calculated as:
M
G () s z −1
A
G z = − z Z = Z
−1
M
M
M
( ) (1 ) [ ] [ ]
2
s z s + BMs+
CM
(5.33)
Finally, the discrete transfer function of controlled plant G ( z ) can be written as:
M
⎛ A ⎞
Md
GM
( z) = ( z−1) ⎜ 2
⎟
⎝ z − BMd
z+
CMd
⎠
(5.34)
BM
2
A
− T
M
s
B
2
M
Where: AM _ d
= e sin( T )
2
s
CM
− ,
B
4
M
CM
−
4
BM
2
− Ts
B
2
M
BM _ d
= 2e cos( Ts CM
− ), CM
_
4
d
BM
Ts
= e − , and
s
T is sampling time.
Hence, the transfer function of closed torque control loop is obtained in the following
form:
G
M _ closed
( z)
M ( z) C ( z) G ( z) D( z) G ( z)
M ( z) 1 + C ( zG ) ( zDzG ) ( ) ( z)
e_
ref M M Mδ
= = =
e M M Mδ
82

Direct Torque Control with Space Vector Modulation
K
K Ψ A ( K + K )( z− )( z−1)
pM
pΨ
PM M _ d pM iM
KpM
+ KiM
3 2
M _ d M _ d pM iM M _ d M _ d pM
= =
[ z − B z + [ A ( K + K ) + C ] z− A K ]( z−1)
K
K Ψ A ( K + K )( z−
)
pM
pΨ
PM M _ d pM iM
KpM
+ KiM
3 2
−
M _ d
+ [
M _ d( pM
+
iM) +
M _ d]
−
M _ d pM
z B z A K K C z A K
(5.35)
Selecting
K Ψ
, K Ψ
will influence poles placement of closed torque control loop and as
p
i
a consequence also torque step responses can be selected.
The transfer function of closed torque control loop is more complicated than flux
control loop (see design of P-flux controller – section 5.2.1). One possibility is use to
the SISO tools from Matlab package to tune parameters of PI torque controller [106].
a) b)
Figure 5.18. a) Torque step response for sampling time T = 200µ
s, b) with denoted rise time,
overshoot and settling time.
As can be observed in (Fig. 5.18) torque response is characterized by overshoot about
40%, rise time 4 samples and settling time 17 samples.
s
To eliminate high overshoot it is recommended to insert at the input prefilter (see
Fig.5.19 ) with transfer function:
z−
b z−
0.6878
PM
( z)
= K = K
K
pM z − 0.855
( z − )
K + K
1
where K =
=0.466 is gain of the prefilter.
z − 0.6878
lim
z→1
z − 0.855
pM
iM
(5.36)
83

Direct Torque Control with Space Vector Modulation
Dz ( )
G ( ) M
z
M
e_ ref( z)
P ( ) M
z
C ( ) M
z
∆δ Ψ
G
Mδ
( z)
U s
z −1
AM
s
ZOH
2
s + B s+
C
M
M
M ( ) e
z
Figure 5.19. Block diagram of torque control loop with prefilter (discrete domain).
Finally, the step response of closed torque control loop with prefilter at the input is
presented bellow:
a) b)
Figure 5.20. a) Torque step response for sampling time T = 200µ
s, b) with denoted rise time,
overshoot and settling time.
s
As it can be observed, the response is now characterized by overshoot about 2%, rise
time 5 samples and the settling time 15 samples.
In digital control system when the sampling time is changed the parameters of
digitalized control plant AMd, BMd,
C
Md
will also change. Therefore, the parameters of PI
torque controller will change also (see Table 2).
Simulation results for digital torque control loop in SABER package for 5KHz with and
without prefilter are shown in Fig. 5.21. Also, the torque step response for different
level of reference torque are presented in Fig. 5.22.
The settings for PI torque controller for different sampling time T s
= 50µ s , 100µ s ,
200µ s , 400µ s are summarized in Table 5.2.
84

Direct Torque Control with Space Vector Modulation
a)
b)
Figure 5.21. Torque step response: a) without prefilter, b) with prefilter.
Figure 5.22. Torque step response with prefilter (from 0 to 25%, 50%, 75% and 100% of
nominal torque).
85

Direct Torque Control with Space Vector Modulation
86

Direct Torque Control with Space Vector Modulation
The torque step response is shown in Fig. 5.23, when parameters of PI torque controller
designed for sampling time T s
= 200µ s were used for control plant for different sampling
times
T
s
= 50µ s , 100µ s , 200µ s , 400µ s , which correspond to switching frequency f s
=
20kHz, 10kHz, 5kHz, 2.5kHz.
Figure 5.23. Torque tracking performance for different sampling time T s
= 50µ s , 100µ s , 200µ s ,
400µ s ( f
s
= 20kHz, 10kHz, 5kHz,2.5kHz) using parameters of PI controller designed for
T
s
= 200µ s ( f
s
= 5kHz). (Please note that for 2.5kHz the system was unstable).
After modification of PI torque controller parameters according to Table 2 it is possible to
achieve better results as shown in Fig. 5.23, what confirms and proper torque tracking
performance in steady and dynamics state.
Figure 5.24. Torque tracking performance for different sampling time T
s
= 50µ s , 100µ s , 200µ s ,
400µ s ( f
s
= 20kHz, 10kHz, 5kHz,2.5kHz) using designed parameters of PI controller calculated
individually (see table 5.2)
87

Direct Torque Control with Space Vector Modulation
5.3 Parallel structure of DTC–SVM scheme
Block scheme of the control structure is shown in Fig. 5.25. Two PI controllers are used for
regulation torque and flux magnitude loops [11,55].
U DC
Ψ s _ ref
e Ψ
ss
U sx _ ref
s
_ ref
S A
S B
M e _ ref
e M
U sy _ ref
U α
U
s β _ ref
S C
θ Ψs
Ψ
I s
M e
γ
m
Figure 5.25. Parallel structure of DTC-SVM scheme.
In this control scheme the reference stator flux magnitude
Ψ and reference
s _ ref
electromagnetic torque
M are compared with estimated values, respectively. The flux and
e_
ref
torque errors eΨ , e are delivered to PI controllers, which generate command value the stator
s
M
voltage components in stator flux coordinates
U , U
_
sx_
ref
sy
ref
. This voltage signals are
transformed to stationary coordinates using the stator flux position angle θ Ψ s
. The reference
stator voltage vector ( U
s α _ ref
, U
s β _
ref
) is delivered to space vector modulator (SVM), which
generates the switching signals
S , S , S to control power transistors of the inverter.
A
B
C
The presented control strategy is based on simplified stator voltage equations described in
stator flux oriented x-y coordinates (equations 2.27a-b):
U
d
Ψ
s
sx
= Rs Isx
+ (5.37)
dt
U = R I +Ω Ψ = R I + E = k M + E
(5.38)
sy s sy Ψs s s sy sy s e sy
88

Direct Torque Control with Space Vector Modulation
where:
k
s
=
2
3
R
p
= Ψ and Esy = ΩΨs Ψ
s
.
s
b
Ψ , 2
M
e
pb s
Isy
s
3
The above equations show that the
U
sx
component has influence only on the change of stator
flux magnitude
Ψ
s
, and the component
Ω
U
sy
– if the term
Ψs s
Ψ is decoupled – can be
used for torque adjustment. Therefore, the flux and torque quantities can be controlled as
shown in Fig. 5.26.
Note, that this DTC-SVM scheme formally corresponds to the stator flux oriented voltage
source inverter-fed drive induction motor. The block diagram of the DTC-SVM scheme with
two PI controllers is shown in Fig. 5.26 The dashed line represents the PMSM part [124,125].
Ψ s _ ref
Ψ s
_
R I
s
sx
U sx
R I
s
sx
∫
Ψ s
R I
s
sy
Ω Ψs
⊗
Ω
Ψ
Ψs s
M e _ ref
_
E sy
U sy
1 Isy
R s
3
2
p
b
Ψ
s _ ref
M e
M e
5.3.1 Digital flux control loop
Figure 5.26. Block diagram of the scheme presented in Fig. 5.27
Putting the stator x-axis current expression from equation 2.28a under the assumption Ld
= L
q
into equation (5.37) one can obtaines
U
I
sx
Ψs
−ΨPM
cos
= (5.39)
L
s
δ Ψ
Ψs −ΨPM cosδΨ
d Ψs R d Ψ
s s RsΨ
PM
= R ( ) + = Ψ + − cosδΨ
(5.40)
L dt L dt L
sx s s
s s s
Using Laplace transformation to equation 5.40 can be written as:
U
sx
Rs RsΨ
PM
=Ψ
s
( s+ ) − cosδ Ψ
. (5.41)
L L
s
s
89

Direct Torque Control with Space Vector Modulation
Assuming small changes of δ Ψ
, the cosδ Ψ
≅ 1, and equations (5.41) reduces to:
Rs
Usx
=Ψ
s
( s+ ) + WΨ
(5.42)
L
s
RsΨ
PM
Where WΨ
=
L
s
The transfer function between stator flux amplitude
Ψ
s
and x-axis of stator voltage is:
G
Ψ
Ψ
s Ls
1 1
() s = = = =
U R
sx
+ WΨ
sLs + Rs
s
s +
s+
A
L
s
Ψ
(5.43)
Where W
Ψ
Rs
Rs
≅ ΨPM
cosδΨ
≅ Ψ
L
L
s
Rs
For motor parameters (see Appendices ): A = Ψ
115.333
L
= .
s
PM
Continuous s-domain
The flux control loop is shown in Fig. 5.27, where CΨ ( s)
is a transfer function of the PI
controller given by [105]:
C
Ψ
s
⎛ K ⎞
iΨ
KpΨ
s+
⎜
K ⎟
pΨ
() s =
⎝ ⎠
(5.45)
s
where
K
iΨ
=
K
T
pΨ
iΨ
Ψ s_ref
C
Ψ
() s
U sx
W Ψ
()
Rs
Ψ
L
s
PM
cos
δ Ψ
+
G
Ψ
s
Ψ s
Figure 5.27. Block diagram of flux control loop in s-domain.
Hence, the transfer function of the closed stator flux amplitude control loop is obtained as:
90

Direct Torque Control with Space Vector Modulation
G
Ψ _ closed
Ψ
s
() s
_ ref CΨ() s GΨ()
s
() s = =
Ψ () s 1 + C () z G () s
s
Ψ
Ψ
(5.46)
Substituting transfer function for CΨ ( s)
and GΨ ( s)
one becomes
G
_ closed
() s
⎛ K ⎞
iΨ
KpΨ
s+
⎜ K ⎟
⎝ pΨ
⎠ ⎛ 1 ⎞
⎛ K ⎞
iΨ
⎜ ⎟ KpΨ
s+
s s A
⎜ K ⎟
⎝ + ⎠ pΨ
=
⎝ ⎠
⎛ K ⎞
s + ( A + Kp
) s+
K
iΨ
KpΨ
s+
⎜ K ⎟
pΨ
⎛ 1 ⎞
1+ ⎝ ⎠
⎜ ⎟
s ⎝s+
AΨ
⎠
Ψ
Ψ
=
2
Ψ Ψ iΨ
(5.47)
Discrete design
z −1
Using backward difference method for discretization process ( s = ) [2] the transfer
Tz
function of equation (5.45) for flux PI controller in discrete system is expressed as:
s
K
pΨ
( KpΨ
+ KiΨ)( z−
)
Tz
K
s
pΨ
+ KiΨ
CΨ( z) = KpΨ(1 + ) =
T ( z−1) ( z−1)
iΨ
(5.48)
Where:
K
s
K
pΨ
iΨ
= Ts; s
Ti
Ψ
Ψ
_
( z)
ref
T - sampling time.
C
Ψ
( z)
W( z)
U sx
Dz ( )
z −1
GΨ
( z)
}
ZOH
1
s+
A Ψ
Ψ
s
( z)
Figure 5.28. Block diagram of flux control loop in discrete domain.
Where: CΨ ( z)
discrete transfer function of PI controller, Dz ( )
1
z − - one sampling time delay
for voltage generation from PWM, and W( z)
- disturbance voltage due to cross coupling
between x-y axis (see Fig. 5.28).
The GΨ ( z)
is discrete transfer function of voltage-flux relationship with zero order hold
(ZOH) block can be calculated as:
91

Direct Torque Control with Space Vector Modulation
−1 GΨ() s z −1
AΨ
GΨ
( z) = (1 − z ) Z[ ] = Z[ ] =
s z A s( s+
A )
( z − )
1 1 AΨ
Z[ ]
z A s( s+
A )
Ψ
Ψ
Ψ
Ψ
(5.49)
Using table of Z transformation [2] one can calculate:
G
Ψ
( z − )
− AΨTs
1 1 z(1 − e ) AΨ
d
( z)
=
− A T
z A
Ψ s
=
z 1( z−e
) z − B
Ψ
( − )
Ψd
(5.50)
Where:
A
Ψd
− AΨT (1 − e s
)
− A Ts
= , BΨ d
e Ψ
A
= and T
s
is sampling time.
Ψ
Hence, the transfer function of closed stator flux control loop can be expressed in the
following form:
G
Ψ _ closed
Ψ
s
( z) _ ref CΨ( z) GΨ( z) D( z)
( z)
= =
Ψ ( z) 1 + C ( z) G ( z) D( z)
s
K
pΨ
( KpΨ + KiΨ) AΨd( z−
)
KpΨ
+ KiΨ
=
K
pΨ
zz ( −1)( z− BΨd) + ( KpΨ + KiΨ) AΨd( z−
)
K + K
Ψ
Ψ
pΨ
iΨ
(5.51)
Now selecting
K Ψ
, K Ψ
is possible to obtain poles placement, which define the dynamic of
p
i
closed torque control loop.
Assuming, that
B
Ψd
=
K
K
pΨ
pΨ
+ K
iΨ
KpΨ − BΨdKpΨ KpΨ
⇒ KiΨ
= = (1 − BΨd)
B B
Ψd
Ψd
and the transfer function of closed stator flux control loop will take the following form:
G
( KpΨ + KiΨ)
AΨd
( z)
=
z − z+ ( K + K ) A
Ψ _ closed
2
pΨ iΨ Ψd
(5.52)
Putting into above equation
K
K
pΨ
pΨ
+ KiΨ
= one obtains:
BΨd
G
K
A
CΨd
K z − z+
C
pΨ
Ψd
BΨd
Ψ _ closed
( z)
= =
2
2
pΨ
z − z+
AΨ
d
BΨd
Ψd
(5.53)
92

Direct Torque Control with Space Vector Modulation
K
pΨ
where CΨd
= AΨd
BΨd
The bellow diagrams shows the relationship between overshoot
time t s
as function
C Ψ d
value.
M
p
, rise time t r
and settling
Please note that t r
is time calculated from 10% to 90% of output signals and t s
is the time in
witch the system transient decay to +-1%.
Figure 5.29. The relationship between overshoot, rise time and settling time versus
amplitude control loop.
From a few values of
C Ψ d
=[0.4046, 0.2688, 0.1720] we can selected C Ψ d
guaranties overshoot 0% and settling time about 10 samples.
C Ψ d
for stator flux
=0.2688, which
1
2
3
Figure 5.30. Flux step response for different values of
C Ψ d
=0.2688, black line (3) C Ψ d
=0.1720.
C Ψ d
: red line (1) C Ψ d
=0.4046, blue line (2)
93

Direct Torque Control with Space Vector Modulation
It corresponds to the transfer function of closed stator flux control loop:
G
C
0.2688
( z)
= =
z − z+ C z −z+ 0.2688
Ψd
Ψ _ closed
2 2
Ψd
(5.54)
Using digitalized motor parameters A
can calculate the parameters of digital PI flux controller as:
Ψd
− AΨT (1 − e s
)
− A Ts
= , BΨ d
e Ψ
A
= and chosen C Ψ d
value, we
Ψ
K
pΨ
C B
A
Ψd
Ψd
= (5.55a)
Ψd
T
iΨ
K
T
pΨ
s
= (5.55b)
K
iΨ
K
iΨ
K
pΨ
= (1 − BΨd)
(5.55c)
B
Ψd
For example with sampling time T = 200µ
s, parameters of PI controller are:
s
K
pΨ
0.2688BΨd
0.2688*0.9772
= = = 1328.64
(5.56a)
A 0.0001977
Ψd
T
iΨ
KpΨTs
1328.64*200µ
s
= = = 8572µ
s
(5. 56b)
K 30.999
iΨ
K
iΨ
K
pΨ
1328.64
= (1 − BΨd) = *(1-0.9772) = 30.999 (5. 56c)
B
0.9772
Ψd
For different sampling time the closed transfer function GΨ _ closed
( z)
of digital flux control
loop should be kept to:
G
C
0.2688
( z)
= =
z − z+ C z −z+ 0.2688
Ψd
Ψ _ closed
2 2
Ψd
(5.57)
In order to find the original function of Z transfer function GΨ _ closed
( z)
using the Z properties
as (sum transformations) [2]:
n
z
z
∑ f kT Z F z Z G z
k = 0
z−1 z−1
−1
z a
= Z [ ]
2
z−1( z − z+
a)
−1 −1
(
s) = [ ( )] = [
Ψ _ closed( )]
(5.58)
94

Direct Torque Control with Space Vector Modulation
As an example the calculated response for 4 samples are given:
n
∑
k = 0
f kT az az a a az
−2 −3 −4
(
s
) = + 2 − (1 + ) + 2 + .......
which gives:
f(0) = f( k) = 0
f(1 T ) = f( k+ 1) = 0
s
f(2 T ) = f( k+ 2) = a=
0.2688
s
f(3 T ) = f( k+ 3) = a( b+ 1) = 2a=
0.5376
s
f T f k a b c a b a a a
2
(4
s
) = ( + 4) =− ( + ) + ( + 1) =− (1 + ) + 2 = 0.7342
Keeping constant
C Ψ d
in equation (5.57) the flux step response for different sampling time
T
s
= 50 µ s,
100 µ s,
200 µ s,
400µ s , which correspond to switching frequency f s
= 20kHz,
10kHz, 5kHz, 2.5kHz are presented.
1
2
3
4
Figure 5.31.Flux step response for different sampling time T
s
= 50 µ s,
100 µ s,
200 µ s,
400µ
s
(switching frequency f
s
20kHz (1 -blue line), 10kHz (2 -green line), 5kHz (3 -red line),2.5kHz
(4 -light blue line).
We may observe from Fig. 5.31 that overshoot is 0% and the settling time is 10 samples.
Selected parameters of PI flux controller for sampling time T s
= 50 µ s,
100 µ s,
200 µ s,
400µ s are summarized in Table 5.3
95

Direct Torque Control with Space Vector Modulation
96

Direct Torque Control with Space Vector Modulation
The behavior of the flux control loop was tested using SABER simulation package. The
model created in SABER takes into account the whole control system, which include real
models of inverter and permanent magnet synchronous motor.
The flux step response is shown in Fig. 5.32., when parameters of PI flux controller designed
for sampling time T s
= 200µ s were used for control plant for different sampling times
T
s
= 50µ s , 100µ s , 200µ s , 400µ s , which correspond to switching frequency f
s
= 20kHz,
10kHz, 5kHz, 2.5kHz.
Figure 5.32. Flux tracking performance for different sampling time T
s
= 50µ s , 100µ s , 200µ s ,
400µ s ( f
s
= 20kHz, 10kHz, 5kHz,2.5kHz) using parameters of PI controller designed for
T
s
= 200µ s ( f
s
= 5kHz).
After modification of PI flux controller parameters according to Table 5.31 it is possible to
achieve better results as shown in Fig. 5.33, what confirms proper flux tracking performance
in steady and dynamics state.
97

Direct Torque Control with Space Vector Modulation
Figure 5.33. Flux tracking performance for different sampling time T
s
= 50µ s , 100µ s , 200µ s ,
400µ s ( f
s
= 20kHz, 10kHz, 5kHz,2.5kHz) using parameters of PI controller calculated individually
(see table 5.3.)
The simulation results in SABER package for Ts=200us is presented in Fig. 5.34.
Figure 5.34. Simulated (SABER) flux tracking performance for step change from 70% - 100% of
nominal flux.
As we can observed from Fig. 5.34 that overshoot is around M
p
= 0% and settling time took
about 10 sampling time of microprocessor, what proved the design procedure of PI digital
flux controller parameters.
98

Direct Torque Control with Space Vector Modulation
5.3.2 Digital torque control loop
The PMSM equations (2.27a,b-2.28a,b) in stator flux coordinates under the assumption
L
d
= L can be written as:
q
U = R I +ΩΨ
Ψ (5.59)
sy s sy s s
0 LI sin
= −Ψ (5.60)
s sy PM
δ Ψ
3
M
e
= pb Ψ I
s sy
(5.61)
2
dΩ m 1 = ( M
e − M
l )
(5.62)
dt J
The load angle can be expressed (Fig. 5.1):
δ = θ − p γ , (5.63)
Ψ
Ψs b m
Where: δ Ψ
is torque angle, θ Ψ s
is stator flux vector position, and γ
m
is rotor position in stator
α,
β coordinates,
p
b
is number of pole pars.
After differentiation equation (5.63) can be written as:
dδ
dθ
dγ
dt dt dt
Ψ Ψs
m
= − pb
(5.64)
δ Ψ
d
dt
δ Ψ
d
=ΩΨ
−p
s bΩm
⇒Ω Ψ s
= + pbΩ m
(5.65)
dt
Putting equations (5.64) and (5.65) into voltage equation (5.59) one obtains:
dδ Usy = RsIsy +Ω ( )
s s
RsI Ψ
Ψ
Ψ =
sy
+ Ψ
s
+ pbΩ m
(5.66)
dt
From equation 0= LI −Ψ sin with assumption that for small angle δ = sinδ
, the
s sy PM
torque angle can be expressed as:
δ Ψ
Ψ
Ψ
L I
s sy
δ Ψ
= (5.67)
Ψ
PM
So, the voltage equation (5.59) becomes:
L dI
s sy
Usy = Rs Isy +ΩΨ
Ψ ( )
s s
= Rs Isy + Ψ
s
+ pbΩm
Ψ dt
PM
(5.68)
99

Direct Torque Control with Space Vector Modulation
After differentiating of the above equation one obtains:
2
sy sy Ls
sy
= Rs + Ψ
s
+ pb
Ψ
PM
dU dI d I dΩm
( )
dt dt dt dt
Take into account that from equation (5.61) the y-axis current is equal
I
sy
(5.69)
2M
e
=
3p
Ψ ,
b
s
dΩ m 1 = ( M
e − M
l ) and under assumption that the motor is no loaded equation (5.69) takes
dt J
form:
dU dM L d M p
( )
dt p dt p dt J
2
sy 2 e s 2
e b
= Rs + Ψ
s
+ M
e
3
b
Ψs ΨPM 3
b
Ψs
(5.70)
Using Laplace transformation and after some arrangements the equation (5.70) can be written:
2Ls
2 2R
Ψs
p
s
b
sU
sy
= M
e
( s + s + )
3p Ψ 3p Ψ J
b PM b s
(5.71)
Hence, the transfer function between electromagnetic torque
be obtained as:
G
M
e
M
2
sy()
M M
M
e
and y-axis voltage
U
sy
can
M () s A s
() s = =
U s s + B s+ C
(5.72)
Where:
A
M
2
3pbΨ
PM
RsΨ
3 Ψ
PM
s
ΨPM pb
= and BM
= CM
=
2L
Ψ L
2JL
s
s
s
s
Using the motor parameters (see Appendices) we may calculates:
A = 198 , B = 115.3 and C = 9065
M
M
M
Continuous s-domain
The torque control loop of the block scheme DTC-SVM from Fig. 5.25 is shown in Fig. 5.35,
where CM
( s ) is a transfer function of the PI controller given by equation 5.28:
100

Direct Torque Control with Space Vector Modulation
M e _ ref
U sy
C () M
s GM
() s
M e
Figure 5.35. Block diagram of torque control loop in s-domain.
The transfer function of torque control loop is obtained as:
G
M _ closed
() s
M () s C () s G () s
e_
ref M M
= = (5.73)
M
e() s 1 + CM() z GM()
s
Substituting in equation (5.73) transfer function for CM
( s)
-Eq.5.28 and GM
( s)
- Eq.5.72 we
may calculate:
KiM
AMKpM( s+
)
K
pM
GM
_ closed()
s = =
2
s + ( B + K A ) s+ C + K A
M pM M M iM M
(5.74)
Discrete design
z −1
Using backward difference method for discretization process ( s = ) the transfer function
Tz
for discrete PI controller is expressed as:
s
K
pM
( z − )
Tz
KpM
+ K
s
iM
CM( z) = KpM(1 + ) = ( KpM + KiM)
T ( z−1) ( z−1)
iM
(5.75)
Where:
K
iM
K
pM
= Ts
- integration gain;
s
TiM
T - sampling time
101

Direct Torque Control with Space Vector Modulation
PI controller
Time delay
G ( ) M
z
Dz ( )} Plant
M
e_ ref( z)
C ( ) M
z
U sy
z −1
AM
s
ZOH
2
s + B s+
C
M
M
M ( ) e
z
Figure 5.36. Block diagram of torque control loop in discrete domain.
Where: CM
( z)
- discrete transfer function for PI controller, Dz ( )
delay for voltage generation from PWM block (see Fig. 5.36).
1
z − - one sampling time
The discrete transfer function G ( z ) for voltage-torque relationship with zero order hold
(ZOH) can be calculated as:
M
G () s z −1
A
G z z Z Z
−1
M
M
M
( ) = (1 − ) [ ] = [ ]
2
s z s + BMs+
CM
⎡
⎤
⎡ ⎤ ⎢ ⎥
⎢ ⎥
z−1 AM
z−1
⎢ A
⎥
M
= Z ⎢
⎥ = Z
2
⎢
2
2
⎥ =
z ⎢⎛ B 2
M ⎞ B ⎥ z
M ⎢ B ⎛
M 2 B ⎞ ⎥
⎢⎜s+ ⎟ + CM
− ⎥
M
2 4 ⎢( s+ ) + CM
−
⎥
⎣⎝ ⎠ ⎦
⎢ 2 ⎜ 4 ⎟
⎣ ⎝ ⎠ ⎥⎦
.
⎡
⎤
2
⎢
B ⎥
M
CM
−
z −1 A ⎢
M
Z
4
⎥
= ⎢
2
2
⎥
z B 2
M
⎢
C B ⎛
M 2 B ⎞ ⎥
M
M
−
4
⎢( s+ ) + CM
−
⎥
⎢ 2 ⎜ 4 ⎟
⎣ ⎝ ⎠ ⎥⎦
(5.76)
Assuming that
have:
B
a = M
and
2
2
BM
b= CM
− , and using table of Z transformation [2] we
4
−aTs
⎡ b ⎤ ze sin( bTs
)
Z ⎢ 2 2 2 aTs
( s a) b
⎥ =
−
⎣ + + ⎦ z − 2 e (cos( bTs
)) z+
e
−2aTs
(5.77)
102

Direct Torque Control with Space Vector Modulation
Finally, the discrete transfer function of controlled plant G ( z ) can be written as:
M
Where:
C
M _ d
⎛ A ⎞
Md
GM
( z) = ( z−1) ⎜ 2
⎟
⎝ z − BMd
z+
CMd
⎠
2
A
BM
A = e sin( T C − ) ,
4
BM
− T
M
s
2
M _ d
2
s M
BM
CM
−
4
BM
Ts
= e − and T
s
is sampling time.
(5.78)
BM
2
− Ts
B
2
M
BM _ d
= 2e cos( Ts CM
− )
4
Hence, the transfer function of closed torque control loop is obtained as:
G
M _ closed
( z)
M ( z) C ( z) G ( z) D( z)
M ( z) 1 + C ( z) G ( z) D( z)
e_
ref M M
= = =
e M M
K
A ( K + K )( z−
)
pM
M _ d pM iM
KpM
+ KiM
3 2
−
M _ d
+ [
M _ d( pM
+
iM) +
M _ d]
−
M _ d pM
=
z B z A K K C z A K
(5.79)
Selecting
K Ψ
, K Ψ
will influence poles placement of closed torque control loop and as a
p
i
consequence also torque step responses can be selected.
The transfer function of closed torque control loop is more complicated than flux control loop
(see design of PI-flux controller – section 5.3.1). One possibility is use to the SISO tools from
Matlab package to tune parameters of PI torque controller [106].
a) b)
Figure 5.37. a) Torque step response for sampling time T = 200µ
s, b) with denoted rise time,
overshoot and settling time.
As can be observed the response is characterized by overshoot about 40%, rise time 4 samples
and settling time 17 samples.
s
103

Direct Torque Control with Space Vector Modulation
The torque control loop should be as fast as possible even with some overshoot. This improve
response to disturbance (for example from flux control loop –see Fig. 5.38).
M
e_ ref( z)
C ( ) M
z
U sy
Dz ( )
z −1
G ( ) M
z
}
AM
s
ZOH
2
s + B s+
C
M
M
M ( ) e
z
Figure 5.38. Block diagram of torque controller in discrete domain with disturbance.
a )
b)
Figure 5.39. Disturbance rejection in torque control loop: a) short voltage impulse, b) voltage
step.
To improve reference tracking performance (without any overshoot) it is recommended to
insert a input prefilter (see Fig. 5.40 ) described by transfer function:
z−b z−b z−0.6663
PM
( z)
= K = K = K
(5.80)
K
pM ( z−a) z−0.855
( z − )
K + K
pM
iM
Where:
K
1
=
=0.43413 is gain of the prefilter.
z − 0.6663
lim
z→1
z − 0.855
104

Direct Torque Control with Space Vector Modulation
M
e_ ref( z)
P ( ) M
z
C ( ) M
z
U sy
Dz ( )
z −1
G ( ) M
z
}
AM
s
ZOH
2
s + B s+
C
M
M
M ( ) e
z
Figure 5.40. Block diagram of torque control loop with prefilter in discrete domain.
Finally, the reference tracking performance of closed torque control loop with prefilter is
presented in Fig. 5.41.
Figure 5.41. a) Reference tracking performance of the torque control loop for sampling time
Ts
= 200µ
s, b) with denoted rise time, overshoot and settling time.
In digital control when the sampling time is changing the parameters of digitalized plant
control AMd, BMd,
C
Md
will be also changes. It is normally that the parameters of PI torque
control will be changes also (see Table 5.4.).
Simulation for PI flux calculated parameters with and without prefilter are shown in Fig.
5.42a-b. Also torque step response for different level of reference torque are presented in Fig.
5.43.
105

Direct Torque Control with Space Vector Modulation
a)
b)
Figure 5.42. Torque step response: a) without prefilter, b) with prefilter.
Figure 5.43. Torque step response with prefilter (from 0 to 25%, 50%, 75% and 100% of nominal
torque).
106

Direct Torque Control with Space Vector Modulation
107

Direct Torque Control with Space Vector Modulation
The behavior of the torque control loop was tested using SABER simulation package.
The model created in SABER takes into account the whole control system, which
include real models of inverter and permanent magnet synchronous motor.
The torque step response is shown in Fig. 5.44, when parameters of PI torque controller
designed for sampling time T s
= 200µ s were used for control plant for different
sampling times T s
= 50µ s , 100µ s , 200µ s , 400µ s , which correspond to switching
frequency f
s
= 20kHz, 10kHz, 5kHz, 2.5kHz.
Figure 5.44. Torque tracking performance for different sampling time T s
= 50µ s , 100µ s ,
200µ s , 400µ s ( f
s
= 20kHz, 10kHz, 5kHz,2.5kHz) using parameters of PI controller
designed for T
s
= 200µ s ( f
s
= 5kHz). Pleas note that for 2.5kHz the system was unstable.
After modification of PI flux controller parameters according to Table 5.4 it is possible
to achieve better results as shown in Fig. 5.45, what confirms proper torque tracking
performance in steady and dynamics state.
108

Direct Torque Control with Space Vector Modulation
Figure 5.45. Torque tracking performance for different sampling time T s
= 50µ s , 100µ s ,
200µ s , 400µ s ( f
s
= 20kHz, 10kHz, 5kHz,2.5kHz) using parameters of PI controller
calculated individually (see Table 5.4.)
109

Direct Torque Control with Space Vector Modulation
5.4 Speed control loop for DTC-SVM structure
The structure of speed control loop for cascade (a) and parallel (b) DTC-SVM scheme
is shown in the Fig. 5.46.
a)
U DC
Ψ s _ ref
U
s α _ ref
S A
S B
Ω m _ ref
M e _ ref
e M
∆δ Ψ
U
s β _ ref
S C
θ Ψs Ψs
I s
M e
I s
d
dt
γ m
U DC
b)
Ψ s _ ref
e Ψ
ss
U sx _ ref
s
_ ref
S A
S B
Ω m _ ref
M e _ ref
e M
U sy _ ref
U α
U
s β _ ref
S C
θ Ψs
Ψ
I s
M e
d
dt
Figure 5.46. Speed control loop for: a) cascade structure of DTC-SVM scheme, b) parallel
structure of DTC-SVM scheme.
γ
m
The mechanical speed equation (2.66) for the PMSM is:
dΩ
M M J dt
m
e
−
L
= (5.81)
Taking Laplace transformation to equation (5.81) one obtains:
M () s − M () s = JsΩ () s
(5.82)
e L m
The transfer function between mechanical rotor speed
M
e
can be expressed as:
Ω
m
and electromagnetic torque
110

Direct Torque Control with Space Vector Modulation
Ωm
1 1
GΩ
() s = = = (5.83)
M Js Js
e
Continuous s-domain
The block diagram of speed control loop is shown in Fig. 5.47, where CΩ ( s)
is a
transfer function of the PI speed controller given by:
KiΩ
KpΩ( s+
)
1
K
pΩ
CΩ() s = KpΩ(1 + ) = (5.84)
T s s
iΩ
and DΩ ( s)
is approximated transfer function of closed torque control loop for cascade
or parallel DTC-SVM structure.
M L
Ω m _ ref
M e ref
_
CΩ () s
D () s
G () s
Ω
M e
Ω
Ω m
Figure 5.47. Block diagram of speed control loop in s-domain.
Discrete design
The transfer function for PI controller in discrete system using backward difference
method for discretization process is expressed as:
K
pΩ
( z − )
KpΩ
+ KiΩ
CΩ( z) = ( KpΩ + KiΩ)
(5.85)
( z −1)
K
pΩ
Where: KiΩ
= Ts- integration and K
p Ω
proportional gain of speed controller, T s
-
T
sampling time.
iΩ
111

Direct Torque Control with Space Vector Modulation
Ω m _ ref
( z)
C
Ω
( z)
Me_ ref()
z
D ( z)
Ω
M ( ) e
z
M ( L
z )
GΩ
( z)
}
ZOH
1
Js
Ω ( z m
)
Figure 5.48. Block diagram of speed control loop in discrete domain.
The closed torque control loop transfer function (see for cascade DTC-SVM -Table 5.2
or parallel DTC-SVM - Table 5.4) is:
0.466*0.37267 a
DΩ ( z) = GM
_ closed( z)
≅ = (5.86)
2 2
z - 1.289z + 0.4633) z +bz + c
The GΩ ( z)
is discrete transfer function for torque-speed relationship with zero order
hold (ZOH). The G ( z)
can be calculated as:
Ω
G () s
1
G z z Z z Z
s
( Js)
s
−1 Ω
−1
Ω( ) = (1 − ) [ ] = (1 − ) [ ]
−1
1 1 1 Ts
= (1 − z ) Z [ ] = J s
2 J ( z −1)
(5.87)
Finally, it can be expressed as:
AΩ
d
GΩ ( z)
=
( z −1)
(5.88)
Ts
Where: A = Ω d
J
, and T
s
is sampling time.
Using the sampling time T = 200µ
s and motor parameters J = 0.0173 the GΩ ( z)
can
be calculated as:
s
G ( z)
Ω
=
0.01156
( z −1)
Hence, the transfer function of closed speed control loop can be written:
G
Ω _ closed
Ωm( z) CΩ( z) DΩ( z) GΩ( z)
( z)
= =
Ω ( z) 1 + C ( z) D ( z) G ( z)
m_
ref
Ω Ω Ω
(5.89)
(5.90)
And finally
112

Direct Torque Control with Space Vector Modulation
G
Ω _ closed
( z)
=
K
pΩ
( K + K )( z−
) aA
K + K
pΩ iΩ Ωd
pΩ
iΩ
2
pΩ
( z−1)( z − bz+ c)( z− 1) + ( KpΩ + KiΩ)( z−
) aAΩ
d
KpΩ
+ KiΩ
K
(5.91)
In many practical cases the digital filter is used in speed measurement loop (see Fig.
5.49).
Ω m _ ref
( z)
PI controller GΩ
( z)
C ( z)
Ω
Me_ ref()
z
D ( z)
Ω
M ( ) e
z
M ( ) L
z
ZOH
Control Plant
}
1
Js
Ω ( z m
)
Torque control loop
F ( z)
Ω
Digital Filter
Figure 5.49. Block scheme of speed control with digital filter in speed measurement loop
(discrete domain).
The transfer function FΩ ( s)
of first order low pass filter in s domain is expressed as:
FΩ () s =
1
s
+ 1
ω
i
(5.92)
2 ω
Where ωc
= 2π
fc
and f
c
is cut off frequency and tan(
cTs
ω
i
= ) [2]. In practice f
c
T 2
is selected in the range 20-250Hz
2( z −1)
Using the Tutsins’s approximation method s = for discretization process, the
Ts
( z+
1)
discrete transfer function of first order low pass filter can be expressed as:
Tsωi
2 −Tsωi
Where: a1
= , b1
=
2 + T ω 2 + T ω
s
i
Tsωi
( z + 1)
2 + T ω a ( z+
1)
FΩ
( z)
= =
z
2 T ω
s
s
s i
1
(5.93)
2 −Tsωi
−
z−
b1
+
i
s
i
113

Direct Torque Control with Space Vector Modulation
The discrete transfer function of digital filter for T = 200µ
s and f = 25Hz
can be
calculated as:
0.015466 (z+1)
FΩ ( z)
= (5.94)
(z-0.9691)
Hence, the transfer function speed control loop with digital filter:
s
c
G
Ω _ closed
Ωm( z) CΩ( z) DΩ( z) GΩ( z)
( z)
= =
Ω ( z) 1 + C ( z) D ( z) G ( z) F( c)
m_
ref
Ω Ω Ω
(5.95)
And finally
G
=
Ω _ closed
( z)
=
K
( K + K )( z−b)( z−
) aA
pΩ
pΩ iΩ 1
Ωd
KpΩ
+ KiΩ
K
z bz c z z z b K K z aA a z
2
pΩ
( − + )( −1)( −1)( −
1) + (
pΩ +
iΩ)( − )
Ωd
1( + 1)
KpΩ
+ KiΩ
(5.96)
Selecting
K Ω
, K Ω
will influence poles placement of closed speed control loop and as
p
i
a consequence also speed step responses can be selected.
In order to select the best value of PI speed controller it is recommended to use the
SISO tools from Matlab package to tune the parameter of PI speed controller.
The speed response with digital filter simulated in SIMULINK is shown in Fig. 5.50
and simulated in SABER in Fig. 5.51 is presented.
114

Direct Torque Control with Space Vector Modulation
Figure 5.50. Simulated (SIMULINK) speed response with digital filter
Figure 5.51. Simulated (SABER) speed response with digital filter in feedback. From the top
reference torque, measured speed.
However, the speed respond is characterized by large overshoot. Therefore, the prefilter
will be applied in order to reduce overshoot (see Fig. 5.52). The discrete transfer
function of prefilter Pz ( ) can be expressed as”
KK _ s 0.005
Pz ( ) = =
z−bb_ s z−0.995
(5.97)
where
KK _ s = lim( z − 0.995)=0.005 is gain of the prefilter.
z→1
115

Direct Torque Control with Space Vector Modulation
Ω m _ ref
( z)
P ( z)
Ω
C ( z)
Ω
M
_
() z
e ref
D ( z)
Ω
M ( ) e
z
M ( L
z )
GΩ
( z)
}
ZOH
1
Js
Ω ( z m
)
F ( z)
Ω
Figure 5.52. Speed response with digital filter in feedback FΩ ( z)
and prefilter PΩ ( z)
at the
input.
The speed response with and without prefilter are shown in Fig. 5.53.
Without prefilter
With prefilter
Figure 5.53. Speed response: blue signal without prefilter and green signal with prefilter at the
input.
Design parameters of PI speed controller for sampling time T
s
= 50 µ s,
100 µ s,
200 µ s,
400µ s are summarized in Table 5.5.
116

Direct Torque Control with Space Vector Modulation
117

Direct Torque Control with Space Vector Modulation
Simulated results for speed tracking performance for different reference speed level in
Fig. 5.54 are shown.
Figure 5.54. Simulated speed tracking performance with prefilter at the input for 10%, 20%,
50%, 100% of nominal speed. From the top actual speed, reference torque.
Investigation for influence of load torque in Fig. 5.55 is presented.
Figure 5.55. Simulated disturbance rejection performance of speed control loop for step load
change 50% of nominal torque. From the top electromagnetic torque, measured speed
Simulation results for speed control loop in Saber package for sampling time
Ts
= 200µ
sand PI speed parameters controller: K Ω
= 1.1940 , T Ω
= 0.0398 (see Table
5.5) in Fig. 5.56 is shown.
p
i
118

Direct Torque Control with Space Vector Modulation
Figure 5.56. Simulated speed tracking performance to step of speed from 0 to 1000rpm.
The presented simulation results confirm well proper operation and design methodology
for digital speed control loop.
119

Direct Torque Control with Space Vector Modulation
5.5 Summary
This Chapter presents design of discrete control loops for two DTC-SVM schemes:
series (cascade) and parallel structures of flux and torque controllers, Fig. 5.2 and
Fig. 5.25, respectively. The cascade structure operates with P-flux controller and PItorque
controller whereas in parallel structure two PI controllers are used. In the first
step of design calculation of discrete Z- transfer function from continuous s- domain
transfer function using zero order hold (ZOH) method of discretization has been
performed. The continuous PI controller transfer function has been discretized using
backward difference approach. Secondly, a SISO tool from Matlab package for
digital controller parameter calculation has been applied. The results of design were
verified by Simulink (using simplified discrete transfer function) and Saber (using
full motor and inverter model) simulation. Also, the influence of sampling time
selection on controller parameters have been discussed. Finally, also the speed
control loop was synthesized using similar methodology.
120

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
Chapter 6
DIRECT TORQUE CONTROL WITH SPACE VECTOR MODULATION (DTC-
SVM) OF PMSM DRIVE WITHOUT MOTION SENSOR
6.1 Introduction
Many motion control applications, such as material handling, packaging and hydraulic
or pneumatic cylinder replacement, require the use of a position transducer for speed or
position feedback, such as an encoder or resolver. In addition, permanent magnet
synchronous motors require position feedback to perform commutation. Some of
systems utilize velocity transducer as well. These sensor add cost, weight, and reduce
the reliability of the system. Also, a special mechanical arrangement needs to be made
for mounting the position sensors. An extra signal wires are required from the sensor to
the controller. Additionally, some type of position sensors are temperature sensitive and
their accuracy degrades, when the system temperature exceed the limits. Therefore, the
research in the area of sensorless speed control of PMSM is beneficial because of the
elimination of the feedback wiring, reduced cost, and improved reliability.
Sensorless speed DTC-SVM control block scheme is presented in Fig. 6.1.
Figure 6.1. Block scheme of DTC-SVM for PMSM drive without motion sensor.
As we can see the operation of speed controlled PMSM drive without mechanical
motion sensor is based only on measurement of following signals, which are available
in every PWM inverter-fed drive system as:
• DC link voltage,
• motor phase currents,
121

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
Based on this signals, the state variable of the drive can be indirectly calculated or
estimated what further allow to achieve the estimated (actual) rotor speed of PMSM.
Two motion sensorless control schemes of DTC-SVM for PMSM drive are presented in
Fig. 6.2
a)
U DC
Ω m _ ref
Ψ s _ ref
M e _ ref
e Me
∆δ
U
s α _ ref
U
s β _ ref
S
A
S B
S C
θ Ψ s
Ψ
s
I
s
U s
M e
I s
Ω m _ est
b)
U DC
Ω m _ ref
Ψ s _ ref
M e _ ref
e Ψ s
e M e
U sx _ ref
U sy _ ref
U
s α _ ref
U
s β _ ref
S A
S B
S C
Ψ
s
M
e
θ Ψs
U s
I
s
Ω m _ est
Figure 6.2. The DTC-SVM block schemes of PMSM without motion sensor: a) cascade
structure and b) parallel structure.
In motion sensorless PMSM drives, as shown in Fig.6.2, the position or speed
transducer (see Fig. 5.52) is replaced by a speed estimation block, which generates the
speed feedback signal Ω into the control systems and stator flux model.
m_
est
122

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
The problem associated with speed sensorless operation of PMSM drive sourced from
VSI is listed below:
• initial rotor flux detection at start up of PMSM controlled drive,
• stator flux estimation without measured speed or position signal from sensor,
• rotor speed estimation based on state variable of the PMSM, especially in low
speed operation region.
Therefore, in this Section the initial rotor position detection method of permanent
magnets, as well as stator flux and rotor speed estimation techniques will be discussed.
6.2 Initial rotor detection method
In a PMSM drive the detection of initial flux position is an important task. The initial
position of the rotor must be detected correctly in order to initialize the flux estimation
procedure. In case of wrong detection the control algorithm has incorrect information
and the rotor shaft can be rotated through few second in positive or negative direction.
This situation is not acceptable in any drive system. Therefore, for the stable starting of
PMSM drive without the temporary reversal rotation, the initial rotor position
estimation is proposed.
The simplest method to achieved the initial rotor flux position is based on the following
rule. For short time the stator winding is supplied by the DC voltage. It impress the DC
current, which generates the magnetic field. The permanent magnet of PMSM sets
accordance with this field line. This position of PM flux is used to set initial values for
the stator flux estimation algorithm.
This method is very simple and not complicate. However, has disadvantage that during
this process the rotor can be moved in unknown direction depending on:
• position of PM before initial detection procedure,
• direction of DC voltage supply into the motor phase.
In order to make the initial rotor flux position correct without any movement (at
standstill) the following algorithms can be used in the literature [78,90,93,95,99,101]
123

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
To estimate the initial rotor position before starting DTC-SVM control, two kind of the
rectangular pulsewise voltages are applied from the inverter to the motor:
• one is the short pulsewise voltage,
• and another is the longer one.
Short pulsewise voltage test
This test based on general principle that the three-phase winding inductances of PMSM
are a function of the mechanical rotor position. Therefore, from the line current
responses in stator oriented coordinates α,
β under the short pulse wise voltage (see
Fig. 6.3) the position of PM can be estimated.
Figure 6.3. Voltage pulse wise during short time voltage test.
During the short time (100µ s ) the vector V 1
=(100) and opposite V 4
=(011) is
generated by voltage source inverter. The achieved current responses in α,
β system for
two type of PMSM during this test in respect to mechanical rotor position are presented
in Fig. 6.4.
a) b)
Figure 6.4. Current components in stator oriented coordinates α,
β under supplied voltage
vector to the motor for very short time : a) for Ld = Lq
and b) for Lq > Ld.
124

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
The stator current trajectory can be divided in to eight sectors (see Fig. 6.5a).
7π
8
o
157 .5
Sector 4
− +
5π
8
o
112 .5
Sector 3
− −
3π
8
+ −
o
67 .5
Sector 2
I s
o
22.5
π
8
+ +
Sector 5
o
202 .5
9π
8
Sector 6
+ −
o
247 .5
11π
8
− −
Sector 7
− +
o
292 .5
13π
8
+ +
Sector 8
Sector 1
15π
8 π
−
8
o
337 .5
Figure 6.5. Stator current trajectory.
Based on the measured response of phase currents in α,
β coordinates, the
( ) sign I s
+ and sign( I s
)
− is calculated from following formulas:
sign( I )
+ = I + I − I
(6.1)
s sα
sβ
s
sign( I )
− = I −I − I
(6.2)
s sα
sβ
s
The possible combinations of sign( I s
)
+ and sign( I s
)
− are shown in Fig. 6.6:
Figure 6.6. Possible combination of sign( I s
)
+
sign( I s
) + + − −
−
sign( I s
) + − − +
+ and sign( I s
)
− under short pulse supply.
Let us assuming, for example, the case where the sign( I s
)
+ and sign( I s
)
− have positive
π π
sign. The position γ m
exist in the domain of − ~ or 7 π 9 π
~ and two estimated
8 8 8 8
position can be obtained.
The mathematical analysis of I , I
sα
sβ waveforms leads to following equations:
I = I +∆ I cos 2γ
(6.3)
sα
s s m
125

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
I
=∆ I sin 2γ
(6.4)
sβ
s m
where I
s
is DC component in Is
α
and
∆ Is
amplitude of fluctuated component.
The current components in α,
β can be modeled as an average value of I
s
plus some
offset value
∆
Is
, as a function of the mechanical position γ m
. Fig. 6.4 shows the
current components given in Eq. (6.1-2), which are function of the phase angle 2γ m
.
Solving those Eqs. (6.1 and 6.2) in respect to mechanical rotor position, two domains of
mechanical rotor position can be obtained as:
Isβ
γ
m1
=
(6.5)
2( I − I )
sα
s
= Isβ
γ
m2
2( I − I )
+ π (6.6)
sα
s
The estimated rotor positions for other combination of sign( I s
)
summarized in Table 6.7.
+ and sign( I s
)
− are
15π
π
−
8 8
7π
9π
−
8 8
π 3π
−
8 8
9π
11π
−
8 8
3π
5π
−
8 8
11π
13π
−
8 8
5π
7π
−
8 8
13π
15π
−
8 8
I sα − I s + I sβ
+ +
+ −
−
−
−
I sα
− I s − I sβ
γ m
+
I sβ
2(
I sα
− I s )
I sβ
+ π
2(
I sα
− I s )
I sα
+ I s π
− +
2I
sβ
4
I sα
+ I s 5π
− +
2I
sβ
4
I sβ
π
+
2( I s α − I s ) 2
I sβ
3π
+
2( I s α − I s ) 2
I sα
+ I s 3π
− +
2I
sβ
4
I sα
+ I s 7π
− +
2I
sβ
4
Table 6.7. Mechanical rotor position calculations.
Long pulsewise voltage test
This test help us to choose the proper estimated value of mechanical rotor position from
two values calculated during the short pulse wise voltage test.
126

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
This test is based on the saturation effect of magnetic circuit. If the long pulse will send
in direction of north pole of PM, the current response will be slower, and, if current
response will be faster that the previous, it means that it is south pole.
6.3 Stator flux estimation methods
Flux estimation is an important task in implementation of high-performance DTC-SVM
motor drives. Vector control method of PMSM drive needs knowledge about actual
value of the stator flux magnitude and position as well electromagnetic torque. Also, the
flux estimation is needed to calculate the actual rotor speed for sensorless operation.
6.3.1 Overview
Many different technique has been developed for PMSM flux estimation [107].
Generally, they may be divided into two groups: open loop estimators and closed loop
estimators/observers. Most of these method are based on so called “current model” or
“voltage model” [110,113]. In fact closed loop estimators/observers are based on the
current or voltage model with an error correction loop, which drives error between two
flux models to zero in steady state. However, an observer has its own dynamics, is
sensitive to parameter changes, and has to be carefully designed for individual drives.
Therefore, for commercially manufactured drives is to complicated and impractical.
This is the reason why in this Chapter only open loop flux estimators will be
considered.
6.3.2 Current model based flux estimator
The block scheme of the current based flux model is presented in Fig. 6.7. It requires:
• knowledge of PMSM machine inductance L , L ,
• speed or position signal,
• PMSM phase currents.
d
q
This kind of flux estimator served in experimental test as a master (standard) to run the
DTC-SVM scheme with speed sensor.
127

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
Ψ PM
I sA
ABC
Is
α
αβ
Isd
L d
Ψ sd
dq
Ψ sα
Ψ s
I sB
αβ
Is
β
dq
I sq
L q
Ψ sq
αβ
Ψ sβ
θ Ψs
γ m
Figure 6.7. Current model based stator flux estimator.
6.3.3 Voltage model based flux estimator with ideal integrator
The stator flux linkage can be obtained by using terminal voltages and currents. It is the
integral of terminal voltages minus the resistance voltage drop:
dΨ sα
= ( Us α − RsIs
α
)
(6.7)
dt
dΨ sβ
= ( Usβ
− RsIsβ
)
(6.8)
dt
However, at low speed (frequencies) some problems arise, when this technique is
applied, since the stator voltage becomes very small and the resistive voltage drops
become dominant, requiring very accurate knowledge of the stator resistance R s
and
very accurate integration. The stator resistance can vary due to temperature changes.
This effect can also be taken into consideration by using the thermal model of the
machine. Drifts and offsets can greatly influence the precision of integration. The
overall accuracy of the estimated flux linkage vector will also depend on the accuracy
of the monitored voltages and currents.
The most know classical voltage model obtains the flux components in stator
coordinates ( α,
β ) by integrating the motor back electromotive force E , E
sα
sβ (see Fig.
6.8). The method is sensitive for only one motor parameter, stator resistance R s
.
However, the application of pure integrator is difficult because of dc drift and initial
value problems.
128

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
I sA
ABC
Is
α
I sB
αβ
Is
β
R s
U sA
U sB
ABC
αβ
U
s α
U
s β
R s
Es
α
Es
β
∫
∫
Ψ sα
Ψ sβ
Ψ s
θ Ψs
Figure 6.8. Voltage model based estimator with ideal integrator.
There are proposed many improvements of the classical voltage model. Some of them
are presented bellow.
6.3.4 Voltage model based flux estimator with low pas filter
A common way to improve the stator flux voltage based model is to use a first-order
low-pass filter (LP) instead of the pure integrator. The equations (6.7 and 6.8) are
transferred to the form:
dΨ sα
= ( Us α − RsIs α)
+ Fc Ψ
sα
(6.9)
dt
dΨ sβ
= ( Usβ − RsIsβ)
+ Fc Ψ
sβ
(6.10)
dt
The block diagram of the estimator is presented in Fig. 6.9. Discrete time
implementation of the integrator becomes:
zΨ ( z) =Ψ ( z) + ( U − R I ) T
(6.11)
sα sα sα s sα
s
zΨ ( z) =Ψ ( z) + ( U − R I ) T
(6.12)
sβ sβ sβ s sβ
s
A LP filter does not give high accuracy at frequencies lower than cutoff frequency
ω = 2π F . There will be errors both in the magnitude and in the phase angle. As
c
c
results, the proposed voltage estimator with LP filter can be used successfully only in a
limited speed range above cutoff frequency
129

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
I sA
ABC
Is
α
I sB
αβ
Is
β
R s
F c
U sA
U sB
ABC
αβ
U
s α
U
s β
-
-
R s
Es
α
Es
β
-
-
∫
∫
Ψ sα
Ψ sβ
Cartesian
To
Polar
Ψ s
θ Ψs
Stator flux estimator (improved voltage model)
F c
Figure 6.9. Voltage model based estimator with low-pass filter.
Discrete time implementation of the LP filter becomes:
zΨ ( z) =Ψ ( z) + ( U − R I ) T + FΨ ( z)
(6.13)
sα sα sα s sα s c sα
zΨ ( z) =Ψ ( z) + ( U −R I ) T −FΨ ( z)
(6.14)
sβ sβ sβ s sβ s c sα
6.3.5 Improved voltage model based flux estimator
Many other methods were developed in order to eliminate dc-offset and initial values
problems [107]. In general, the output Y of these new integrators (Fig. 6.10) is
expressed as:
ωc
1
Y = X + Y
s+ ω s+
ω
c
c
lim
(6.15)
Where X is the input and Y is output of the integrator respectively. The Y lim
is a
compensation signal used as a feedback and ω
c
is cutoff frequency.
X
1
s + ω c
Y
Compensation
signal
ωc
s +ω
c
Y lim
− lim
lim
Figure 6.10. Improved integration method with saturation block.
130

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
The first part of the equation represents a LP filter. The second part realizes a feedback,
which is used to compensate the error in the output. The block diagram of new
integration algorithm with saturation block is shown in see Fig. 6.11.
I sA
ABC
Is
α
I sB
αβ
Is
β
R s
ωc
s +ω
c
Y lim
− lim
lim
U sA
U sB
ABC
αβ
U
s α
U
s β
R s
Es
α
Es
β
1
s + ω c
1
s + ω c
Ψ α _ comp
Ψ s α
Ψ sβ
Ψ s
θ Ψs
Ψ
s β _ comp
ωc
s +ω
c
Y lim
− lim
lim
Figure 6.11. Full block diagram of voltage model based estimator with saturation block on the
α,
β components.
The main task of saturation block is to stop the integration when the output signal
Ψ
sα
or
Ψ
sβ
exceed the reference value of stator flux amplitude. Please note that if the
compensation signal is set to zero, the improved integrator represents a first-order LPfilter.
If the compensation signal
integrator operates as a pure integrator.
Ψ or Ψ
β _
is not zero the improved
sα
_ comp
s
comp
Discrete time implementation of the improved integrator becomes:
ω
zΨ ( z) =Ψ ( z) + ( U − R I ) T + ( Ψ ( z) −Ψ ( z))
(6.16)
c
sα sα sα s sα s sα _lim
sα
Ts
ω
zΨ ( z) =Ψ ( z) + ( U − R I ) T + ( Ψ ( z) −Ψ ( z))
(6.17)
c
sβ sβ sβ s sβ s sβ _lim
sβ
Ts
The output of saturation block can be described as:
⎧Ψ
sαβ
( z) if (Ψ
sαβ(z))=lim
sαβ
131

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
Where lim is the limited value. Please note that lim should be set at reference stator
flux amplitude
lim =Ψ
s _ ref
equal Ψ
PM
.
From Eq. 6.18 it can be observed that when one of the stator flux linkage components
Ψ
sα
or Ψ
sβ
exceeds the limit, it causes in distortion of the output waveform.
6.4 Electromagnetic torque estimation
The PMSM motor output torque is calculated based on the equations (2.51), (2.52),
(2.55), (2.58) presented in Chapter 2, which for stator oriented coordinate system can be
written as follows:
3
M
e
= pb( Ψs αIsβ −Ψ
sβIsα
)
(6.19)
2
It can be seen that calculated torque is dependent on the current measurement accuracy
and stator flux estimation method. In practice current measurement is performed with
high accuracy ( ≤ 1% with 150kHz frequency bandwidth) using, for example, LEM
sensors.
6.5 Rotor speed estimation methods
6.5.1 Overview
High performance operation of motion sensorless PMSM drives depends mainly on
accurate knowledge of rotor PM flux magnitude, position and speed. The rotor position
estimation methods can be classified into two major groups:
• motor model based,
• rotor saliency based techniques.
The rotor saliency based approach is suitable only for the Interior PMSM (see Fig. 1.2 c
and d). Motor model based approach detect the back EMF vector, which includes
information about position and speed, using either open loop models/estimators
[81,85,86] or closed loop estimators/observers [70, 73,74,96,97,100]. Also adaptive
observers [72,92,98,83] and Extended Kalman Filters (EKF) [67,73] have been
proposed for motor position and speed estimation. These methods, however, are
computationally intensive and require careful design and proper initialization.
132

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
Therefore, for commercially manufactured drives are impractical and further a simple
open loop based techniques will be considered.
6.5.2 Back electromotive force (BEMF) technique
This technique uses the back electromotive force to estimate the rotor speed [70]. The
velocity signal could be integrated to generate a position estimate. However, this signal
is sensitive to parameter variations and tends to drift and have offset problem. Another
problem with using BEMF to estimate position is that at zero speed the BEMF goes to
zero and at low speed the signal to noise ratio can not be ignored.
6.5.3 Stator flux based technique
Generally, the calculation of rotor speed is based on the simple relationship:
θ = θ − δ , (6.20)
r Ψs Ψs
where θ
r
is electrical position, θ Ψ s
is stator flux position and δ Ψ s
is torque angle.
After differentiation equation (6.20) and taking into account that θ r
= p b
γ m
the
mechanical speed of PMSM rotor can be expressed as:
⎛dθ
⎜
⎝ dt
Ψs
Ω
m
= −
dδΨ
⎞
⎟/
pb
, (6.21)
dt ⎠
dθ s
where Ω Ψ
Ψs
= is angular speed of stator flux vector and δ Ψ
is torque angle.
dt
As we can observe form equation (6.21) in order to calculate the mechanical rotor speed
it is necessary to calculate separately two components. One of them is angular speed of
stator flux vector Ω
Ψs
and the second one is change of the load angle d δ Ψ (see Fig.
dt
6.12).
133

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
Is
α
Is
β
U
s α
Ψ
sα
Ψ sβ
Ω =
Ψs
d
dt
θ Ψ s
dδ Ψ
dt
−
1
p b
Ω m
U
s β
Figure 6.12. Block diagram of stator flux vector based rotor speed estimator.
The synchronous speed
Ω
Ψs
is calculated based on the stator flux estimator:
θ Ψ s
Ψ
sβ
= arctg( )
Ψ
sα
(6.22)
and the calculation of
ΩΨs
can be done as:
The estimation of the synchronous speed
dθ θ
s Ψs( k) −θ
Ψ
Ψs( k−1)
Ω
Ψs
= = (6.23)
dt T
s
Ω
Ψs
based on the derivative of the position of
stator flux space vector can be modified taking in to account equation (6.22), which
finally gives [12]:
Ψ
Ω =
Ψs
sα
dΨ
sβ
−Ψ
dt
Ψ +Ψ
sβ
2 2
sα
sβ
dΨ
dt
sα
(6.24)
Digital implementation of equation (6.24) can be written as:
Ψ Ψ −Ψ Ψ
Ω =
Ψ T
sα( k−1) sβ( k) sβ( k−1) sα( k)
Ψs( k) 2
s
s
(6.25)
Also from equation (2.27b) in stator flux coordinate system the synchronous speed
can be obtained as:
Ω
Ψs
U − R I E
Ω
Ψs
= =
Ψ Ψ
sy s sy sy
s
s
(6.26)
134

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
Calculation of d δ Ψ
dt
Based on the flux-current (Eq. 2.29b) and torque (Eq. 2.58) equations in stator flux
coordinates under consideration that for surface PMSM machine L d
= L q
and making
assumption that for small changes of torque angleδ Ψ
, the sinδ
be written as:
Ψ
= δ , the equations can
Ψ
0= LI −Ψ sin
(6.27)
s sy PM
δ Ψ
3
M
e
= pbΨ PMIsy
(6.28)
2
the torque angle δ Ψ
can be calculated as:
LI
s sy 2ML
e s
δ Ψ
= =
Ψ 3p
Ψ Ψ
PM b s PM
(6.29)
Further, the d dt
δ Ψ is calculated as:
( ) ( 1)
δΨ
Ψ k Ψ k−
d
dt
δ −δ
= (6.30)
T
s
Figure 6.13. Simulated oscillograms of rotor speed estimation according to block scheme from
Fig. 6.12 (in Saber). From the top: synchronous speed Ω
Ψs
, the d signal, the measured and
dt
estimated rotor speed, the measured and estimated electromagnetic torque.
δ Ψ
135

Sensorless Speed Direct Torque Control with Space Vector Modulation (DTC-SVM)
The speed estimation problem is still open, especially at low and zero speed operations.
The accuracy of presented method depends on accuracy of applied stator flux estimation
and differentiation algorithm. It allows, however, for robust start, closed loop operation
above 10% of nominal speed, and braking the drive to zero speed.
6.5 Summary
The main problems associated with PMSM sensorless speed operation are presented
in this Chapter. For robust start of PMSM without the temporary rotor reversal a
special initialization algorithm has been used. This algorithm performs two test:
short and the longer voltage generated by the PWM inverter. The used speed
estimation algorithm is based on stator flux vector and torque angle estimation and
does not operate accurately around zero speed region. However, it allows robust
start and closed speed operation in the speed range above 10% of nominal speed.
The effectiveness of DTC-SVM with and without motion sensor has been proved by
simulation and experimental results (see Chapter 7)
136

DSP implementation of DTC-SVM control
Chapter 7 DSP IMPLEMENTATION OF DTC-SVM CONTROL
7.1 Description of the laboratory test-stand
The basic structure of laboratory setup is presented in Fig. 7.1 and the photo of
laboratory setup is shown in Fig. 7.2. The motor setup consists of 3kW permanent
magnet synchronous motor and DC motor, which is used as a load. The PMSM machine
is supplied by PWM inverter, which is controlled by digital signal processor (DSP)
based on DS1103 board. The voltage inverter is supplied from three-phase diodes
rectifier. The DSP interface is used in order to separate the high power from the low
power circuit (computer part). Please note that the DS1103 is inserted inside the PC
computer.
Figure 7.1. Block scheme of laboratory setup.
Figure 7.2. Laboratory setup. 1-voltage inverter, 2-control for DC motor, 3- PMSM machine, 4
– DSP interface.
137

DSP implementation of DTC-SVM control
The detailed power circuit of the laboratory setup is shown in Fig. 7.3.
Rectifier
Inverter
three-phase
supply
network
U A
U B
U C
D 1
D3
D5
CF
T 1
D 1
T 3
D 3
T 5
D 5
current
sensors
A
B
C
PMSM
D 2
D 4
D 6
T 2
D 2
T 4
D 4
T 6
D 6
speed or position
sensor
SA
SA
SB
SB
SC
SC
Reference speed
DTC-SVM
DS1103
microprocessor
Figure 7.3. Power circuit of the laboratory setup.
In presented system the actual two currents and DC link voltage are measured by LEM
sensors and processed by A/D converter. The rotor position and speed are obtained with
an encoder of 2500 pulse per revolution. All internal data of DSP can be sent through a
D/A converter and displayed in the scope. All data of the PMSM and inverter are given
in the Appendices.
The control algorithm for PMSM machine was written in C language and was
implemented inside the processor. Also, a simple dead-time compensation method and
voltage drop on the semiconductor elements are implemented.
The phase voltage of the motor are reconstructed inside the processor using the
measured DC-link voltage and duty cycles of PWM for each phases. Motor and PI
controller model are given in Appendices.
Various tests have been carried out in order to investigate the drive performance and to
characterize the steady-state and transient behavior.
138

DSP implementation of DTC-SVM control
Control Desk experiment software run on the PC computer provides all function for
controlling, monitoring and automation of real-time experiments and makes the
development of controllers more effective. A Control Desk experiment layout for
control the PMSM machine using DTC-SVM control method is presented in Fig. 7.4.
Figure 7.4. Performed Control Desk experimental layout for control of PMSM drive.
139

DSP implementation of DTC-SVM control
7.2 Steady state behaviour
The experimental steady state no load operation at 25Hz stator frequency for
conventional ST-DTC (Fig. 4.10) and DTC-SVM (Fig. 5.46b) control is presented in
Fig. 7.5. The sampling time has been set at T = 200µ
s for DTC-SVM and T = 25µ
s for
hysteresis based ST-DTC method, respectively.
a)
s
s
b)
Figure 7.5. No load experimental steady state oscillograms at stator frequency 25Hz.
(a) ST-DTC for T = 25µ
s (b) DTC-SVM for T = 200µ
s.
s
From the top: line to line voltage, phase current, amplitude of stator flux, motor torque.
s
As we can observed from Fig. 7.5a the motor phase current characterized by high
current ripples. This is mainly because the inductances of the PMSM is smaller than an
equivalent power induction motor IM. In order to reduce the current ripples the
140

DSP implementation of DTC-SVM control
sampling time of microcontroller should be decrease. However there is hardware
limitation. The loaded program to microprocessor can not run faster.
Using the space vector modulation (SVM) based DTC much better results can be
obtained. Note, that in spite of lower switching frequency DTC-SVM guarantees less
current and torque ripple. This is mainly because contrary to hysteresis operation with
SVM operation, the inverter output voltage is unipolar (compare output voltage
waveform in Fig. 7.5a with Fig. 7.5b). This also reduces semiconductor device voltage
stress and instantaneous current reversal in DC link.
The presented experimental results (Fig. 7.6-7.9) are measured in the system with
measured speed taken to the feedback. These investigations have been performed to
show the behaviour of the DTC-SVM system without influence of the speed estimation.
In Fig. 7.6 and Fig. 7.7 steady state operation for different values of the mechanical
speed and load torque are shown.
Figure 7.6. Experimental steady state operation of PMSM controlled via DTC-SVM with the
encoder speed signal taken to the feedback ( Ω = 300rpm
, M = 0 ).
m
From the top: 1) amplitude of stator flux (0.25Wb/div), 2) electromagnetic torque (10Nm/div),
3) line to line voltage (1000V/div), 4) stator phase current (10A/div).
l
141

DSP implementation of DTC-SVM control
Figure 7.7. Experimental steady state operation of PMSM controlled via DTC-SVM with the
encoder speed signal taken to the feedback ( Ω = 300rpm
, M = 10Nm
-50% of nominal
torque). From the top: 1) amplitude of stator flux (0.25Wb/div), 2) electromagnetic torque
(10Nm/div), 3) line to line voltage (1000V/div), 4) stator phase current (10A/div).
m
l
Figure 7.8. Experimental steady state operation of PMSM controlled via DTC-SVM with the
encoder speed signal taken to the feedback ( Ω = 1500rpm
, M = 0 ). From the top: 1)
amplitude of stator flux (0.25Wb/div), 2) electromagnetic torque (10Nm/div), 3) line to line
voltage (1000V/div), 4) stator phase current (10A/div).
m
l
142

DSP implementation of DTC-SVM control
Figure 7.9. Experimental steady state operation of PMSM controlled via DTC-SVM with the
encoder speed signal taken to the feedback ( Ω = 1500rpm
, M = 10Nm
-50% of nominal
torque). From the top: 1) amplitude of stator flux (0.25Wb/div), 2) electromagnetic torque
(10Nm/div), 3) line to line voltage (1000V/div), 4) stator phase current (10A/div).
m
l
7.3 Dynamic behaviour
The experimental results of flux and torque control loop obtained in dynamic states for
PMSM machine controlled via two different DTC-SVM schemes are presented.
7.3.1 Flux and torque control loop
Cascade DTC–SVM control scheme (Fig. 5.46a)
In order to show behaviour of the system the dynamic testes for the flux and torque
controllers has been carried out for sampling time, T = 200µ
s. It corresponds to
switching frequency f = 5kHz
. Please note that the flux digital controller parameters
s
were selected according to Table 5.1 and the torque digital controller parameters were
selected from Table 5.2. (see Chapter 5.2).
In Fig. 7.10 stator flux tracking performance is presented. This result is comparable
with simulation results presented in Fig. 5.15.
s
143

DSP implementation of DTC-SVM control
Figure 7.10. Experimental flux tracking performance of PM synchronous motor for zero speed
at sampling time Ts
= 200µ
s. Reference flux from70% to 100% of nominal value . From the
top: 1-reference stator flux amplitude(0.05Wb/div), 2- stator flux amplitude (0.05Wb/div), 3-
electromagnetic torque (2Nm/div), 4- motor phase current (10A/div).
In Fig. 7.11 torque tracking performance is presented. The achieved result is
comparable with simulation results presented in Fig. 5.21b.
Figure 7.11. Experimental torque tracking performance of PM synchronous motor for zero
speed at sampling time Ts
= 200µ
s. Reference torque from 0 to nominal value. From the top:1-
reference torque (10Nm/div), 2 - estimated torque (10Nm/div), 3 - stator flux amplitude
(0.1Wb/div), 4- motor phase current (10A/div)
144

DSP implementation of DTC-SVM control
Figure 7.12. Experimental torque tracking performance of PM synchronous motor for zero
speed (zoom) at sampling time Ts
= 200µ
s. Reference torque from 0 to nominal value. From
the top:1- reference torque (10Nm/div), 2 - estimated torque (10Nm/div), 3 - stator flux
amplitude (0.1Wb/div), 4- motor phase current (10A/div).
All experimental results presented in Fig. 7.10-7.12 confirm very well proper and stable
operation of flux and torque control loops for cascade DTC-SVM structure.
145

DSP implementation of DTC-SVM control
Influence of sampling time for torque control loop in cascade DTC-SVM
The influence of sampling time on experimental torque tracking performance is
illustrated in Fig. 7.13-7.15. The dynamic test has been carried out for the same
condition ( Ω = 0rpm
) as for simulation shown in Fig. 5.24. The controller parameters
m
has been set according to Table 5.2. In all oscilograms we may see proper operation of
the torque control loop for different sampling time used in practice.
Figure 7.13. Experimental torque tracking performance of PM synchronous motor for zero
speed at sampling time Ts
= 100µ
s( fs
= 10kHz
). Reference torque from 0 to nominal value.
From the top:1- reference torque (4Nm/div), 2 - estimated torque (4Nm/div)
Figure 7.14. Experimental torque tracking performance of PM synchronous motor for zero
speed at sampling time Ts
= 200µ
s( fs
= 5kHz
). Reference torque from 0 to nominal value.
From the top:1- reference torque (4Nm/div), 2 - estimated torque (4Nm/div)
146

DSP implementation of DTC-SVM control
Figure 7.15. Experimental torque tracking performance of PM synchronous motor for zero
speed at sampling time T = 400µ
s( f = 2.5kHz
). Reference torque from 0 to nominal value.
s
From the top:1- reference torque (4Nm/div), 2 - estimated torque (4Nm/div).
s
Parallel structure of DTC–SVM scheme (Fig. 5.46b)
Dynamic testes for the flux and torque controllers were carried out for sampling time
Ts
= 200µ
s, which corresponds to switching frequency f = 5kHz
. Please not that the
digital flux controller parameters were selected according to Table 5.3 and the digital
torque controller parameters were selected from Table 5.4. (see Chapter 5.3).
In Fig. 7.16 stator flux tracking performance is presented. This result is comparable
with simulation results presented in Fig. 5.33 for flux and Fig. 5.45 for torque loop,
respectively.
s
147

DSP implementation of DTC-SVM control
Figure 7.16. Experimental flux tracking performance of PM synchronous motor for zero speed
at sampling time Ts
= 200µ
s. Reference flux from70% to 100% of nominal value . From the
top: 1-reference stator flux amplitude(0.05Wb/div), 2- stator flux amplitude (0.05Wb/div), 3-
electromagnetic torque (2Nm/div), 4- motor phase current (10A/div).
It can be observed that achieved result is comparable with simulation results presented
in Fig. 5.34.
Figure 7.17. Experimental torque tracking performance of PM synchronous motor for zero
speed at sampling time Ts
= 200µ
s. Reference torque from 0 to nominal value. From the top:1-
reference torque (10Nm/div), 2 - estimated torque (10Nm/div), 3 - stator flux amplitude
(0.1Wb/div), 4- motor phase current (10A/div)
148

DSP implementation of DTC-SVM control
Figure 7.18. Experimental torque tracking performance of PM synchronous motor for zero
speed (zoom) at sampling time Ts
= 200µ
s. Reference torque from 0 to nominal value. From
the top:1- reference torque (10Nm/div), 2 - estimated torque (10Nm/div), 3 - stator flux
amplitude (0.1Wb/div), 4- motor phase current (10A/div).
The achieved result is comparable with simulation results presented in Fig. 5.42b.
Experimental results presented in Fig. 7.16-7.18 confirm very well the effectiveness of
controller design and proper operation of flux and torque control loops for DTC-SVM
structure.
149

DSP implementation of DTC-SVM control
Influence of sampling time for torque control loop in parallel DTC-SVM
The influence of sampling time on experimental torque tracking performance is
illustrated in Fig. 7.19-7.21. The dynamic test has been carried out for the same
condition ( Ω = 0rpm
) as for simulation shown in Fig. 5.45. The controller parameters
m
has been set according to Table 5.4. In all oscilograms we may see proper operation of
the torque control loop for different sampling time used in practice.
Figure 7.19. Experimental torque tracking performance of PM synchronous motor for zero
speed at sampling time Ts
= 100µ
s. Reference torque from 0 to nominal value. From the top:1-
reference torque (4Nm/div), 2 - estimated torque (4Nm/div)
Figure 7.20. Experimental torque tracking performance of PM synchronous motor for zero
speed at sampling time Ts
= 200µ
s. Reference torque from 0 to nominal value. From the top:1-
reference torque (4Nm/div), 2 - estimated torque (4Nm/div).
150

DSP implementation of DTC-SVM control
Figure 7.21. Experimental torque tracking performance of PM synchronous motor for
zero speed at sampling time Ts
= 400µ
s. Reference torque from 0 to nominal value.
From the top:1- reference torque (4Nm/div), 2 - estimated torque (4Nm/div).
In this Chapter the results of experimental verification 3kW PMSM drive with two
DTC-SVM schemes has been presented. As shown the drive performance confirms
applied design methodology. The performance of both cascade and parallel DTC-SVM
control structure are similar. However, parallel structure has been selected for industrial
manufacturing because of :
• less noisy control algorithm (differentiation required in cascade structure –
equation (5.6) is eliminated),
• stator flux control in closed loop,
• the same structure can be used for IM and PMSM control (universal control for
AC motors).
7.3.2 Speed control loop
0peration with speed sensor
Dynamic testes for the speed control loop were measured for sampling timeT
= 200µ
s.
Please note that the digital speed controller parameters were selected according to Table
5.5 (see Chapter 5.4). In Fig. 7.22-7.26 rotor speed tracking performance are presented.
s
151

DSP implementation of DTC-SVM control
Figure 7.22. Experimental start up and breaking to zero speed of PMSM motor controlled via
DTC-SVM with the encoder speed signal taken to the feedback ( Ω
m
= 0rpm → 300rpm
). From
the top: 1- reference speed (180rpm/div), 2- measured speed (180rpm/div), 3 - electromagnetic
torque (20Nm/div), 4- motor phase current (20A/div).
Figure 7.23. Experimental speed reversal for PMSM motor controlled via DTC-SVM with the
encoder speed signal taken to the feedback ( Ω
m
=−300rpm → 300rpm
). From the top: 1-
reference speed (360rpm/div), 2- measured speed (360rpm/div), 3 - electromagnetic torque
(20Nm/div), 4- motor phase current (20A/div).
152

DSP implementation of DTC-SVM control
Figure 7.24. Experimental start up and breaking to zero speed of PMSM motor controlled via
DTC-SVM with the encoder speed signal taken to the feedback ( Ω
m
= 0rpm → 1500rpm
).From
the top: 1- reference speed (900rpm/div), 2- measured speed (900rpm/div), 3 - electromagnetic
torque (20Nm/div), 4- motor phase current (20A/div).
Figure 7.25. Experimental speed reversal for PMSM motor controlled via DTC-SVM with the
encoder speed signal taken to the feedback ( Ω
m
=−1500rpm → 1500rpm
). From the top: 1-
reference speed (1800rpm/div), 2- measured speed (1800rpm/div), 3 - electromagnetic torque
(20Nm/div), 4- motor phase current (20A/div).
153

DSP implementation of DTC-SVM control
Figure 7.26. Experimental speed reversal for PMSM motor controlled via DTC-SVM with the
encoder speed signal taken to the feedback ( Ω
m
=−1200rpm → 1200rpm
). From the top: 1-
stator flux component Ψ
sα
(0.25Wb/div), 2- measured speed (360rpm/div), 3 - electromagnetic
torque (20Nm/div), 4- motor phase current (20A/div).
Figure 7.27. Experimental response to load torque step change of PMSM motor controlled via
DTC-SVM with the encoder speed signal taken to the feedback at Ω = 0rpm
.
( M = 0Nm → 10Nm
) From the top: 1- reference speed (180rpm/div), 2- measured speed
l
(180rpm/div), 3 - electromagnetic torque (5Nm/div), 4- motor phase current (10A/div).
m
154

DSP implementation of DTC-SVM control
Figure 7.28. Experimental response to load torque step change of PMSM motor controlled via
DTC-SVM with the encoder speed signal taken to the feedback at Ω = 300rpm
.
( M = 0Nm → 10Nm
).From the top: 1- reference speed (180rpm/div), 2- measured speed
l
(180rpm/div), 3 - electromagnetic torque (5Nm/div), 4- motor phase current (10A/div).
m
Figure 7.29. Experimental response to load torque step change of PMSM motor controlled via
DTC-SVM with the encoder speed signal taken to the feedback at Ω = 1500rpm
.
( M = 0Nm → 10Nm
). From the top: 1- reference speed (900rpm/div), 2- measured speed
l
(900rpm/div), 3 - electromagnetic torque (5Nm/div), 4- motor phase current (10A/div).
m
155

DSP implementation of DTC-SVM control
Experimental results presented in Figures 7.22-7.29 well confirm the effectiveness of
developed controller synthesis methodology and proper operation of speed control loop
for parallel DTC-SVM structure.
Sensorless speed operation
Dynamic tests for the speed control loop without motion sensor were measured for
sampling time T = 200µ
s. Please note that the digital speed controller parameters were
s
selected exactly like for operation with speed sensor according to Table 5.5 (see
Chapter 5.4).
The results of speed estimator dynamic test are presented in Fig. 7.30. In this test speed
controller operates with the encoder signal in feedback and speed estimator works in
open loop fashion.
Figure 7.30. Experimental dynamic test of the speed estimation. Speed reversal Ω =± 300rpm
From the top: 1- reference speed (180rpm/div), 2- measured speed (180rpm/div), 3 - estimated
speed (180rpm/div), 4- error of estimated speed (5%/div).
m
The typical dynamic performance tests of sensorless DTC-SVM drive has been
illustrated in Fig. 7.31-7.36. Start up and breaking to zero speed for different speed level
are shown in Fig. 7.31 and 7.33. Also, the speed reversal for low (Fig. 7.32) and
nominal (Fig. 7.34) speed are presented. The half load torque step change tests are
shown in Fig. 7.36 and 7.37.
156

DSP implementation of DTC-SVM control
Figure 7.31. Experimental start up and breaking to zero speed of PMSM motor controlled via
DTC-SVM with the estimated speed signal taken to the feedback ( Ω
m
= 0rpm → 300rpm
).
From the top: 1- reference speed (180rpm/div), 2- measured speed (180rpm/div), 3 -
electromagnetic torque (20Nm/div), 4- motor phase current (20A/div).
Figure 7.32. Experimental speed reversal for PMSM motor controlled via DTC-SVM with the
estimated speed signal taken to the feedback ( Ω
m
=−300rpm → 300rpm
). From the top: 1-
reference speed (360rpm/div), 2- measured speed (360rpm/div), 3 - electromagnetic torque
(20Nm/div), 4- motor phase current (20A/div).
157

DSP implementation of DTC-SVM control
Figure 7.33. Experimental speed step response for PMSM motor controlled via DTC-SVM with
the estimated speed signal taken to the feedback ( Ω
m
= 0rpm → 1500rpm
).From the top: 1-
reference speed (900rpm/div), 2- measured speed (900rpm/div), 3 - electromagnetic torque
(20Nm/div), 4- motor phase current (20A/div).
Figure 7.34. Experimental speed reversal for PMSM motor controlled via DTC-SVM with the
estimated speed signal taken to the feedback ( Ω
m
=−1500rpm → 1500rpm
). From the top: 1-
reference speed (1800rpm/div), 2- measured speed (1800rpm/div), 3 - electromagnetic torque
(20Nm/div), 4- motor phase current (20A/div).
158

DSP implementation of DTC-SVM control
Figure 7.35. Experimental speed reversal for PMSM motor controlled via DTC-SVM with the
estimated speed signal taken to the feedback ( Ω
m
=−1200rpm → 1200rpm
). From the top: 1-
stator flux component Ψ
sα
(0.25Wb/div), 2- measured speed (360rpm/div), 3 - electromagnetic
torque (20Nm/div), 4- motor phase current (20A/div).
Figure 7.36. Experimental response to load torque step change of PMSM motor controlled via
DTC-SVM with the estimated speed signal taken to the feedback at Ω = 300rpm
.
( M = 0Nm → 10Nm
). From the top: 1- reference speed (180rpm/div), 2- measured speed
l
(180rpm/div), 3 - electromagnetic torque (10Nm/div), 4- motor phase current (10A/div).
m
159

DSP implementation of DTC-SVM control
Figure 7.37. Experimental response to load torque step change of PMSM motor controlled via
DTC-SVM with the estimated speed signal taken to the feedback at Ω = 1500rpm
.
( M = 0Nm → 10Nm
). From the top: 1- reference speed (900rpm/div), 2- measured speed
l
(900rpm/div), 3 - electromagnetic torque (10Nm/div), 4- motor phase current (10A/div).
m
Experimental results presented in Figures 7.31-7.37 confirm very well the effectiveness
of speed estimation algorithm of speed control loop for parallel DTC-SVM structure.
160

Summary and closing conclusion
Chapter 8
SUMMARY AND CLOSING CONCLUSIONS
This thesis studied basic problems related to selection, investigation and implementation
of the PWM inverter-fed permanent magnet synchronous motor (PMSM) drives suitable
for serial manufacturing. The selected method should provide: robust start and operation
in wide speed control range including zero speed, with and without mechanical motion
sensor; guarantee good and repeatable parameters of PMSM drive for wide power range
(1-100kW). The control and protection algorithm should be implemented in simple and
cheap microprocessor.
To solve so formulated general task several related problems had to be solved. At first,
the space vector based mathematical description and static characteristic of PMSM
under different control modes were studied (Chapter 2). Secondly, the three phase
voltage source inverter model including nonlinearities (dead time, semiconductor
voltage drop, DC link pulsation) and pulse with modulation PWM techniques were
presented (Chapter 3). Next, based on the study of most important high performance
control methods as field oriented control (FOC), and direct torque control (DTC), the
method called direct torque control with space vector modulator (DTC-SVM) has been
selected for further consideration (Chapter 4). This methods combines main advantage
of FOC (space vector modulator and fixed switching frequency) and DTC (simple
structure, rotor parameter independent), as well as eliminates disadvantages like:
coordinate transformation, the need of internal current control loops, high sampling
time, high torque and current ripple at steady state operation, etc.
Consequently the most important contribution of this thesis is included in the Chapter 5,
where the two basic variants of DTC-SVM schemes: series (cascade) and parallel
structures of flux and torque controllers are presented (Fig. 5.2 and Fig. 5.25). Also, a
systematic methodology of digital controller design for these both DTC-SVM variants
have been given. This methodology has been verified by Simulink (using simplified
discrete transfer function) and Saber (using full motor and inverter model) simulation
studies. The influence of sampling time selection on controller design has also been
discussed.
The main problems associated with PMSM sensorless speed operation are presented in
the Chapter 6. It should be noted that PMSM differ from IM drives mainly in:
• PMSM parameters strongly depend on construction ,
161

Summary and closing conclusion
• position of PM flux has to be known prior to start up to achieve smooth
operation.
Therefore, for robust starting of PMSM without the temporary rotor reversal a simple
initialization algorithm has been used. This algorithm performs two tests: short and the
longer voltage pulses generated by the PWM inverter. The used speed estimation
algorithm is based on stator flux vector and torque angle estimation and does not
operate accurately in zero speed region. However, it allows robust start and closed
speed operation in the speed range above 10% of nominal speed. For application where
high performance operation around zero speed are required, the DTC-SVM drive with
motion sensor (encoder) is recommended. The effectiveness of DTC-SVM scheme with
and without motion sensor has been proved by simulation and experimental results
(Chapter 7).
Simulation study and experimental results have shown that from two variants of DTC-
SVM schemes the parallel structure is more flexible in torque and flux controller
design. Also, because of lack of the differentiation in the main control path (compare
Fig. 5.2 and Fig. 5.25), it is less sensitive to noise which inherently associates signal
processing in power electronic converters. Therefore, the parallel structure has been
selected for industrial manufacturing and implemented on digital signal processor
(DSP).
Thanks to direct torque and flux control structure the described control is suitable to
almost all – industrial applications including electrical vehicles (for example hybrid
cars).
Finally, it should be stressed that the developed system was brought into serial
production. Presented algorithm DTC-SVM has been used in new generation of inverter
drives produced by Polish company Power Electronic Manufacture – “TWERD”,
Toruń.
162

Appendices
APPENDICES
A1 Rotor and stator of PMSM machine
A1.1. View of rotor (on the left side) and stator armature (on the right side) of PMSM.
A2 Basic transformation
y
β
K B
K
θ K
K β
K y
K
x
K α
θ K
α
x
K A
K C
Fig. A2.1. Space vector representation in stationary α,
β coordinates and synchronous
rotating x,
y coordinates.
A, BC , ⇒ xy ,
2
Kx = [ KA cos θ
K
+ KB cos(2 π / 3 − θK ) + KC cos(2 π / 3 + θK
)]
3
2
Ky = [ KA cos( π /2 + θK ) + KB cos(2 π /3 − ( π /2 + θK )) + KC cos(2 π /3 + π /2 + θK
)]
3
163

Appendices
2
Kx = [ UA cos θK + KB cos( θK − 2 π / 3) + KC cos( θK
+ 2 π / 3)]
3
2
Ky = [ −KA sin θK −KB sin( θK − 2 π / 3) − KC sin( θK
+ 2 π / 3)]
3
⎡K
A ⎤
⎡K
x ⎤ ⎡ cosθK cos( θK − 2 π / 3) cos( θK
+ 2 π / 3) ⎤⎢
K
⎥
⎢ B
K
⎥ = ⎢
y sinθK sin( θK 2 π / 3) sin( θK
2 π / 3)
⎥⎢ ⎥
⎣ ⎦ ⎣− − − − + ⎦
⎢⎣
K ⎥
C ⎦
ABC , , ⇒ α,
β θ = 0
2 1 1
K = ( KA B C)
3 − α
K K
2 − 2
2 3 3 1
K = ( KB C) (
B C)
3 2 − β
K K K
2 = 3
−
⎡ 1 1 ⎤
1 − − ⎡ K
A ⎤
⎡Kα
⎤ 2 ⎢ 2 2 ⎥
⎢
K
⎥
⎢ B
K
⎥ = ⎢
⎥
β 3 3 3 ⎢ ⎥
⎣ ⎦ ⎢
⎥
⎢0
− ⎢⎣K
⎥
C ⎦
⎣ 2 2
⎥
⎦
x, y⇒
A, B,
C
K
KA = Kxcosθ K
+ Kycos( π / 2 + θK) = KxcosθK − KysinθK
KB = Kxcos(2 π / 3 − θK) + Kycos(2 π / 3 − ( π / 2 + θK)) = Kxcos( θK −2 π / 3) −Kysin( θK
−2 π / 3)
K = K cos(2 π / 3 + θ ) + K cos(2 π / 3 + ( π / 2 + θ )) = K cos( θ + 2 π / 3) − K sin( θ + 2 π / 3)
C x K y K x K y K
⎡K
A⎤ ⎡ cosθK −sinθK
⎤
⎢ K
x
K
⎥ ⎢
B
cos( θK 2 π / 3) sin( θK
2 π / 3)
⎥ ⎡ ⎤
⎢ ⎥
=
⎢
− − −
⎥⎢ K
⎥
y
⎢K
⎥ ⎢
C
cos( θK + 2 π / 3) − sin( θK
+ 2 π / 3) ⎥ ⎣ ⎦
⎣ ⎦ ⎣ ⎦
α, β ⇒ ABC , , θ = 0
K
KA
=
K α
1 3
KB
=− Kα
+ K
2 2
1 3
KC
=− Kα
− K
2 2
β
β
164

Appendices
⎡ ⎤
⎢ 1 0 ⎥
⎡K
A ⎤ ⎢ ⎥
⎢ 1 3 K
K
⎥ ⎢ ⎥ ⎡ α ⎤
⎢ B ⎥
=
⎢
−
2 2 ⎥⎢ K
⎥
β
⎢K
⎣ ⎦
⎣ ⎥
C ⎦ ⎢ ⎥
⎢ 1 3⎥
⎢
− −
⎣ 2 2 ⎥⎦
xy , ⇒ α,
β
Kα = Kxcosθ K
+ Kycos( π / 2 + θK) = KxcosθK − KysinθK
Kβ = Kxcos( π / 2 − θK) + KycosθK = KxsinθK + KycosθK
⎡Kα
⎤ ⎡cosθK
−sinθ
K
K⎤⎡ x ⎤
⎢
K
⎥ = ⎢
β sinθK
cosθ
⎥⎢ K
⎥
⎣ ⎦ ⎣
K ⎦⎣ y ⎦
α, β ⇒ x,
y
Kx = Kα cosθ K
+ Kβ cos( π / 2 − θK) = Kα cosθK + Kβ
sinθK
Ky = Kα cos( π / 2 + θK) + Kβ cosθK =− Kα sinθK + Kβ
cosθK
⎡Kx ⎤ ⎡ cosθK sinθ
K
K ⎤⎡ α ⎤
⎢
K
⎥ = ⎢
y sinθK cosθ
⎥⎢ K
⎥
⎣ ⎦ ⎣−
K⎦⎣ β ⎦
A3 Model of PM synchronous motor
# This template models the permanent magnet synchronous motor(pmsm)
# t1,t2 and t3 are motor input terminals
# rotor speed (in rad/s) is the output
# rs-Stator windings' resistence per phase(in Ohms)
# ld-d_axis inductance(in H)
# lq-q_axis inductance(in H)
# pm-Rotor magnet flux(Wb)
# j-Moment of inertia(in kgm2)
# d- Damping constant (Nm/rad/s)
# tl-motor load (Nm)
# p-Number of pole pairs
# power-The total input power (W)
# Assumptions:No core losses,no saturation,thermal effects
#(rs,ld,lq and pm values are constants)
element template pmsm_dtc t1 t2 t3 t0 speed out_me out_psi out_thetam out_ualf out_ubet out_psia
out_psib out_ialf out_ibet out_theta
out_tl=rs,ld,lq,pm,d,tl,j,p,init_theta_m,init_omega_m,omega_m_const
electrical t1,t2,t3,t0 # motor input terminals and stator neutral point
output nu speed,
out_me,out_psi,out_thetam,out_ualf,out_ubet,out_psia,out_psib,out_ialf,out_ibet,out_theta,out_tl
number rs=0.692,ld=6m,lq=6m,pm=0.26379,d=0.002044,tl=0.0,j=0.003,p=3.0,
init_theta_m=0.0,init_omega_m=0.0,omega_m_const=0.0
{
165

Appendices

Appendices
i(t1->t0)+=it1
it1: it1=id*cos(p*theta_m)-iq*sin(p*theta_m)
i(t2->t0)+=it2
it2: it2=id*cos(p*theta_m-y)-iq*sin(p*theta_m-y)
i(t3->t0)+=it3
it3: it3=id*cos(p*theta_m+y)-iq*sin(p*theta_m+y)
omega_m: (te-tl-d*omega_m)/j=d_by_dt(omega_m)
#omega_m: omega_m=omega_m_const*math_pi/30.0
theta_m: omega_m=d_by_dt(theta_m)
speed: speed=omega_m
out_me: out_me=te
out_psi: out_psi=psi
out_thetam:out_thetam=theta_m
out_ualf: out_ualf=va
out_ubet: out_ubet=vb
out_psia: out_psia=psia
out_psib: out_psib=psib
out_ialf: out_ialf=ialf
out_ibet: out_ibet=ibet
out_theta: out_theta=theta
out_tl: out_tl=tl
}
}
A4 Motor parameters
Surface type motor
Power P 3kW
Number of pole pairs p 3
Phase current I(rms) 6.9A
Phase voltage U(rms) 70V
Magnetic flux-linkage Ψ 0.264 Wb
Rotor speed
Ω
PM
m
3000rpm
Nominal torque Me 20Nm
Moment of the inertia J 0.0174kgm 2
Stator winding resistance Rs 0.692 Ω
Stator d-axis inductance Ld 6mH
Stator d-axis inductance Lq 6mH
Interior type motor
Power P 2,2kW
Number of pole pairs p 3
Phase current I(rms) 4.1A
Rated voltage U(rms) 380V
Magnetic flux-linkage Ψ 0.4832 Wb
Rotor speed
Ω
PM
m
1750rpm
Nominal torque Me 12Nm
Moment of the inertia J 0.010074kgm 2
Stator winding resistance Rs 3.3 Ω
Stator d-axis inductance Ld 41.59mH
Stator d-axis inductance Lq 57.06mH
167

Appendices
A5 Voltage Source Inverter parameters
Detailed date of IGBT transistors (module TOSHIBA M675Q2YS50):
U = 1200V
, I = 75A
CE
C
UCEsat
= 2.8 − 3.6V
, forward diode voltage 2.4 − 3.5V
Turn on time tON
= 0.2µ
s, Turn off time tOFF
= 0.6µ
s
Delay of IGBT drivers tONd
= 0.5µ
s TOFFd
= 1µ
s
TON = tON + tONd
= 0.7µ
s total turn on time of IGBT
TOFF = tOFF + tOFFd
= 1.6µ
s total turn off time of IGBT
Dead time T = 2.5µ
s
A6 PI speed controller
d
The commonly used in industrial application speed controller is a Proportional-Integral
PI controller thanks to possibility to reduce the speed error between the reference (
and actual rotor speed (
X
ref
)
X
m
) to zero (see Fig. A6.1). The output signal of controller is a
reference torque, which has upper and lower limitation for this value equal the nominal
torque or more than 130% of nominal torque. The output of the speed controller acts as
a current reference command for the current controllers. This current command is
limited to a nominal current of the motor.
The speed controller demands produce proper electromagnetic torque.
X_
ref
Reference signal
−
error signal
Kp
Y L
Controller output
a)
Feedback signal
X_
m
Kp
T
i
∫
X _ ref
Reference signal
b)
−
Feedback signal
error signal
X _ m
Kp
Kp
T
i
∫
1
Y NL
lim_max
lim_ min
−
Y E
Y L
Controller output
T i
A.6.1. General structure of Proportional- Integral controller without antiwindup (a) and with
antwindup (b).
168

Appendices
A7 PWM technique – six step mode
Six-stepped-voltage waveforms are rich in harmonics. These time harmonics produce
respective stator current harmonics, which in turn interact with fundamental air gap
flux, generating harmonics torque pulsations. The torque pulsations are undesirable:
they generate audible noise, speed pulsations, and losses. In case of supplied motor by
using only active vectors (six step mode) we can observed non sinusoidal current, which
generates torque ripples with frequency of six time fundamental frequency of supplied phase
voltages.
Fig. A7.1. Experimental operation in six step mode. From the left side stator voltages in α , β
coordinates, From right side voltage trajectory.
Fig. A7.2. Experimental operation in six step mode. From the left side stator currents in α , β
coordinates, From right side stator current trajectory.
169

Appendices
Fig. A7.3. Experimental operation in six step mode. From the top: α stator voltage,
electromagnetic torque in machine, phase current.
170

List of symbol
List of symbols
Symbol (general)
X - instantaneous value
X
N
- normalized value
X - vector
X ∗ - conjugate vector
X - amplitude of vector
Re( X ) – real part of X
Im( X ) – imaginary part of X
Symbol (special)
α,
β - stator fixed system
dq , - rotor reference system
x,
y - general reference system
Ls
- stator inductance
Zs
- stator impedance
M
s
- mutual inductance
Is
- phase current value
U
s
- phase voltage value
Ψ
s
- phase flux value
P -active power
Q - reactive power
S - apparent power
Pe
- electro-magnetic power
Ωm
- mechanical rotor speed
Ωs
- synchronous speed
cosφ - power factor
δ
I
, δ Ψ
- torque angle
φ - power angle
R
s
- stator resistance
L
d
, L
d
- direct and quadrature inductances
θr
- electrical rotor position
γ
m
- mechanical rotor position
pb
- number of pole pairs
Ψ
PM
- rotor flux of permanent magnets
M
e
- electromagnetic torque
M
es
- synchronous torque
M
er
- reluctance torque
M
l
- load torque (external load torque)
M
d
- dynamic torque
Jm
- motor moment of inertia
Jl
- load moment of inertia
J - moment of inertia of total system (sum of J
m
and J
l
)
171

List of symbol
Subscripts
A , B , C - denote arbitrary phase quantities in a system of natural coordinate ABC. , ,
d , q - arbitrary direct and quadrature components in a system of rotor coordinate dq. ,
α , β - arbitrary alpha and beta components in a system of stator coordinate α,
β .
x , y - denote arbitrary components in a system of general coordinate x,
y .
.. r - denotes value of rotor
.. s - denotes value of stator
.. _max - maximum value
.. _min – minimum value
.. _ ref - reference value
.. _ est - estimated value
.. _ amp -amplitude value
.. _ rms - root mean square value
.. _ LL - line to line value
* - reference value
^ - estimated value
Abbreviations
RSM – reluctance synchronous motor
BLDCM – blushless DC motor
PMSM – permanent magnet synchronous motor
IPMSM - interior permanent magnet synchronous motor
SPMSM - surface permanent magnet synchronous motor
EMF – electro-magnetic force
VSI - voltage source inverter
SVM – space vector modulator
PWM – pulse width modulation
PWM-VSI – voltage source inverter with PWM
DTC - direct torque control
DTC-SVM - direct torque control with space vector modulator
RFOC - rotor field oriented control
SFOC - stator field oriented control
CTAC - constant torque angle control
MTPAC - maximum torque per ampere control
UPFC - unity power factor control
CSFC - constant stator flux control
172

References
Books and PhD-Thesis
1. I. Boldea, S.A. Nasar, “Electrical drives”, CRC Press, 1999.
2. G. F. Franklin, J. D. Powell, “Feedback Control of Dynamic Systems”, Addison-
Wesley Publishing Company, 1994.
3. M. P. Kaźmierkowski, H. Tunia, “Automatic Control of Converter-Fed drives”,
Elsevier, 1994.
4. M. P. Kaźmierkowski, R. Krishnan, F. Blaabjerg, “Control in Power Electronics”,
Academic Press, 2002.
5. P.C. Krause, O. Wasynczuk, S.D. Sudhoff, “Analysis of Electric Machinery”, IEEE
Press, 1995.
6. R. Krishnan, “Electric Motor Drives – Modeling, Analysis, and Control”, Prentice Hall,
2001.
7. T. Kaczmarek, K. Zawirski, “Układy napędowe z silnikiem synchronicznym”,
Politechnika Poznańska, 2000.
8. M. Malinowski, "Sensorless Control Strategies for Three-Phase PWM Rectifiers", PhD
Thesis, Warsaw University of Technology, 2001.
9. Ch. Perera, “Sensorless Control of Permanent Magnet Synchronous Motor Drives”,
PhD thesis, Institute of Energy Technology, 2002.
10. K. Rajashekara, A. Kawamura, K. Matsuse, “Sensorless Control of AC Motor Drives”,
IEEE Press,1996, USA.
11. M. Zelechowski, “Space Vector Modulated – Direct Torque Controlled (DTC-SVM)
Inverter-Fed Induction Motor Drive”, Phd Thesis , Warsaw University of Technology,
2005.
12. P. Vas, “Sensorless Vector and Direct Torque Control”, Oxford University Press, 1998.
Model of PMSM
13. C. Bowen; Z. Jihua; R. Zhang; “Modeling and simulation of permanent magnet
synchronous motor drives”, Electrical Machines and Systems, ICEMS 2001, Vol. 2, 18-
20 Aug. 2001, pp. 905 – 908.
14. P. Pillay, R. Krishnan, “Modeling of permanent magnet motor drives”, Transactions on
Industrial Electronics, Vol. 35, Nov. 1988, pp. 537 – 541.
15. N. Urasaki; T. Senjyu; K Uezato, “An Accurate Modeling for Permanent Magnet
Synchronous Motor Drives”, Applied Power Electronics Conference and Exposition,
APEC 2000, Vol. 1, 6-10 Feb. 2000, pp:387 – 392.
16. A.H. Wijenayake, P.B. Schmidt, “Modeling and analysis of permanent magnet
synchronous motor by taking saturation and core loss into account”, International
Conference on Power Electronics and Drive Systems, Vol. 2, 26-29 May 1997, pp. 530
– 534.
Inverter and PWM Modulation
17. L. Ben-Brahim, "The analysis and compensation of dead-time effects in three phase
PWM inverters", Industrial Electronics Society, 1998. IECON '98. Proceedings of the
24th Annual Conference of the IEEE, Vol. 2d, 31 Aug.-4 Sept. 1998, pp.792-797.
18. V. Blasko, "Analysis of a hybrid PWM based on modified space-vector and trianglecomparison
methods", IEEE Transactions on Industry Applications, Vol. 33, Issue: 3,
May-June 1997, pp.756-764.
19. F.Blaabjerg, J.K. Pedersen, P. Thoegersen, “Improved modulation techniques for
PWM-VSI drives”, IEEE Transactions on Industrial Electronics, Feb. 1997, Vol. 44, pp.
87 – 95.
20. J.W. Choi; S.K. Sul, "Inverter output voltage synthesis using novel dead time
compensation", IEEE Transactions on Power Electronics, Vol. 11, Issue: 2, March
1996, pp.221-227.
173

References
21. D.W. Chung, J. Kim, S.K. Sul, "Unified voltage modulation technique for real-time
three-phase power conversion", IEEE Transactions on Industry Applications, Vol. 34,
Issue: 2, March-April 1998, pp.374-380.
22. A. Diaz, E.G. Strangas, "A novel wide range pulse width overmodulation method [for
voltage source inverters]", Applied Power Electronics Conference and Exposition,
APEC 2000.
23. A.M. Hava, R.J. Kerkman, T.A. Lipo, "A high performance generalized discontinuous
PWM algorithm", Applied Power Electronics Conference and Exposition, APEC '97
Conference Proceedings 1997, Twelfth Annual, Vol. 2, 23-27 Feb. 1997, pp.886-894.
24. D.G. Holmes, "The significance of zero space vector placement for carrier-based PWM
schemes", IEEE Transactions on Industry Applications, Vol. 32, Issue: 5, Sept.-Oct.
1996, pp.1122-1129.
25. J. Holtz, "Pulsewidth modulation for electronic power conversion", Proceedings of the
IEEE, Vol. 82, Issue: 8, Aug. 1994, pp.1194-1214.
26. J. Holtz, W. Lotzkat, A.M. Khambadkone, "On continuous control of PWM inverters in
the overmodulation range including the six-step mode", IEEE Transactions on Power
Electronics, Vol. 8, Issue: 4, Oct. 1993, pp.546-553.
27. J.L. Lin, "A new approach of dead-time compensation for PWM voltage inverters",
IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications,
Vol. 49, Issue: 4, April 2002, pp.476-483.
28. D.C. Lee, G.M. Lee, "A novel overmodulation technique for space-vector PWM
inverters", IEEE Transactions on Power Electronics, Vol. 13, Issue: 6, Nov. 1998,
pp.1144-1151.
29. M. Malinowski "Adaptive modulator for three-phase PWM rectifier/inverter", in proc.
EPE-PEMC Conf., Kosice, 2000, pp.1.35-1.41.
30. N. Urasaki, T. Senjyu, K. Uezato, T. Funabashi, “On-line dead-time compensation
method for permanent magnet synchronous motor drive”, IEEE International
Conference on Industrial Technology, CIT '02, 11-14 Dec. 2002, Vol. 1, pp. 268 – 273.
FOC, DTC and DTC-SVM
31. F. Blaschke, "The principle of fields-orientation as applied to the Transvector closedloop
control system for rotating-field machines", in Siemens Reviev 34, 1972, pp.217-
220.
32. G.S. Buja, M.P. Kazmierkowski, "Direct Torque Control of PWM Inverter-Fed AC
Motors-A Survey", IEEE Transactions on Industrial Electronics, Vol. 51, Aug. 2004,
pp.744-757.
33. D. Casadei, G. Serra, A. Tani, L. Zarri, F. Profumo, “Performance analysis of a speedsensorless
induction motor drive based on a constant-switching-frequency DTC
scheme”; IEEE Transactions on Industry Applications, March-April 2003, Volume 39,
pp. 476 – 484.
34. D. Casadei, F. Profumo, G. Serra, A. Tani, “FOC and DTC: two viable schemes for
induction motors torque control”; IEEE Transactions on Power Electronics, Sept. 2002,
Vol. 17, pp. 779 – 787.
35. T. Chun; W.Y. Hu, “Research on the direct torque control in electromagnetic
synchronous motor drive”, Power Electronics and Motion Control Conference, PEMC
2000, 15-18 Aug. 2000. Vol. 3, pp. 1262 – 1265.
36. S.K. Chung, H.S. Kim, Ch. G. Kim, M.-J. Youn, “A New Instantaneous Torque Control
of PM Synchronous Motor for High-Performance Direct-Drive Applications” IEEE
Transaction on Power Electronics, May 1998, Vol. 12, No. 3.
37. Ch. French, P. Acarnley, “Direct Torque Control of Permanent Magnet Drivers”, IEEE
Transaction on Industrial Application, Vol. 32, No. 5, 1996, pp. 1080-1088.
38. L.-H. Hoang, “Comparison of field-oriented control and direct torque control for
induction motor drives”, Industry Applications Conference, 1999. Thirty-Fourth IAS
Annual Meeting, 3-7 Oct. 1999, Vol. 2, pp.1245 – 1252.
174

References
39. J.-I. Itoh, N. Nomura, H. Ohsawa, “A comparison between V/f control and positionsensorless
vector control for the permanent magnet synchronous motor”, Proceedings of
the Power Conversion Conference, PCC Osaka 2-5 April 2002, Vol. 3, pp. 1310 – 1315.
40. T.M. Jahns, “Motion control with permanent-magnet AC machines”, Proceedings of the
IEEE, Aug. 1994, Vol. 82, pp. 1241 – 1252.
41. C. Lascu, I. Boldea, F. Blaabjerg, “Variable-structure direct torque control - a class of
fast and robust controllers for induction machine drives”, IEEE Transactions on
Industrial Electronics, Aug. 2004,Vol. 51, pp. 785 – 792.
42. T. Lixin; Z. Limin; M.F. Rahman, Y. Hu, “A novel direct torque control for interior
permanent-magnet synchronous machine drive with low ripple in torque and flux-a
speed-sensorless approach”, IEEE Transactions on Industry Applications, Nov.-Dec.
2003,Vol. 39, pp. 1748 – 1756.
43. J. Lee; Ch.-G. Kim; M.-J. Youn, “A dead-beat type digital controller for the direct
torque control of an induction motor”, IEEE Transactions on Power Electronics, Sept.
2002, Vol. 17, pp. 739 – 746.
44. H. Rasmussen, P. Vadstrup, H. Børsting, "Adaptive sensorless field oriented control of
PM motors including zero speed", IEEE International Symposium on Industrial
Electronics, Ajaccio, France, May 2004.
45. M. F. Rahman, L. Zhong, K. W. Lim, “A direct torque controlled interior permanent
magnet synchronous motor drive incorporating field weakening”, IEEE-IAS’97, 1997 ,
p. 67 - 74.
46. M.F. Rahman, L. Zhong, K.W. Lim, “An Investigation of Direct and Indirect Torque
Conttrollers for PM Synchronous Motor Drivers”, Conference PEDS’97, May,
Singapure, pp. 519-523.
47. M. F. Rahman, L. Zhong, “Comparison of torque response of the interior permanent
magnet motor under PWM current and direct torque control”, IEEE, 1999 , p. 1464 -
1470.
48. J. Rodriguez, J. Pontt, C. Silva, S. Kouro, H. Miranda, “A novel direct torque control
scheme for induction machines with space vector modulation”; Power Electronics
Specialists Conference, PESC’04, 20-25 June 2004, Vol. 2, pp. 1392 – 1397.
49. D. Telford, M.W. Dunnigan, B.W. Williams, “A comparison of vector control and
direct torque control of an induction machine”, 2000 IEEE 31st Annual Power
Electronics Specialists Conference, 2000. PESC 00, 18-23 June 2000, Vol. 1, pp. 421 –
426.
50. L. Takahashi, T. Noguchi, “A new quick response and high efficiency strategy of
induction motor”, IAS, 1985 , p. 495-502.
51. L. Tang; L. Zhong; F. Rahman, “Modeling and experimental approach of a novel direct
torque control scheme for interior permanent magnet synchronous machine drive”,
IECON 02, 5-8 Nov. 2002, Vol. 1, pp. 235 – 240.
52. M.N. Uddin, T.S. Radwan, G.H. George, M.A. Rahman, “Performance of current
controllers for VSI-fed IPMSM drive”, IEEE Transactions on Industry Applications,
Volume 36, Issue 6, Nov.-Dec. 2000 Page(s):1531 – 1538.
53. L. Xu, M. Fu “A Novel Sensorless Control Technique for Permanent Magnet
Synchronous Motors (PMSM) using Digital Signal Processor (DSP)”, NEACON’97,
Dayton, Ohio, July 14-17, 1997.
54. L. Xu, M. Fu “A Sensorless Direct Torque Control Technique for Permanent Magnet
Synchronous Motors”, IEEE Industrial Applications Conferences,1999, Vol. 1.
55. Y. Xue, X. Xu, T.G. Habetler, D.M. Divan, “A low cost stator flux oriented voltage
source variable speed drive”, Industry Applications Society Annual Meeting, 7-12 Oct.
1990, Vol.1, pp. 410 – 415.
56. H. Yuwen,T. Cun, G. Yikang,Y. Zhiqing, L.X. Tang, M.F. Rahman “In-depth Reserch
on Direct Torque Control of Permanent Magnet Synchronous Motor”, IEEE Conference
IECON 2002 ,November, Sevilla, Spain.
175

References
57. L. Zhong, M.F. Rahman, W.Y. Hu, Lim, K.W., Rahman, M.A., “A direct torque
controller for permanent magnet synchronous motor drives”, IEEE Transactions on
Energy Conversion, Sept. 1999, Vol. 14, pp. 637 – 642.
58. L. Zhong, M.F. Rahman, W.Y.Hu, K.W. Lim, “Analysis of Direct Torque Control in
Permanent Magnet Synchronous Motor Drives”, IEEE Transaction on Power
Electronics, Vol. 12, No. 13, May, 1997.
59. L. Zhong, M.F. Rahman, „Voltage Switching Tables for Controlled Interior Permanent
Magnet Motor“, IEEE Conference IECON’99, Vol. 3, pp. 1445-1451.
60. L.Zhong, M. F. Rahman, “A direct torque controller for permanent magnet synchronous
motor drives without a speed sensor”, IEEE Trans. On Energy Conversion, Vol.14,
No.3 ,September , 1999.
61. M.R Zolghardi, D.Diallo, D.Roye, “Direct Torque Control System for Synchronous
Machine”, Europen Conference on Power Electronics and Application,
EPE’97,Trondheim, Norway,Vol. 3, pp. 694-699.
62. M.R. Zolghardi, E. M. Olasagasti, D. Roye, “Steady State Torque Correction of Direct
Torque Controlled PM Synchronous Machine”, IEEE International Conference on
Electric Machines Drivers, 1997.
63. M.R Zolghardi, J.Guiraud, J.Davoine, D. Roye, “A DSP Direct Torque Controller for
Permanent Magnet Synchronous Motor Drivers”, IEEE Power Electronics Conference,
PESC’98, Vol.2, pp. 2055-2061.
64. M.R Zolghardi, D. Roye, “Sensorless Direct Torque Control of Synchronous Motor
Drive” Proceedings of International Conference on Electrical Machines ICEM’ 98, Vol.
2, 1998, Istanbul, Turkey.
65. M.R Zolghardi, D. Diallo, D. Roye, “Start up Strategies for Direct Torque Controlled
Synchronous Machine ”, Europen Conference on Power Electronics and Application,
EPE’97,Trondheim, Norway,Vol. 3, pp. 689-693.
66. M.R Zolghardi, C. Pelissu, D. Roye, “Start up of a global Direct Torque Control
System”, IEEE Power Electronic Conference PESC’96, Vol. 1, pp. 370-374.
Sensorless Control and Initial Position Estimation
67. S. Bolognani, R. Oboe, M. Zigliotto, “Sensorless full-digital PMSM drive with EKF
estimation of speed and rotor position”, IEEE Transactions on Industrial Electronics,
Vol. 46, Issue 1, Feb. 1999, pp. 184 – 191.
68. S. Brock, J. Deskur, K. Zawirski, , “Robust speed and position control of PMSM”,
Proceedings of the IEEE International Symposium on Industrial Electronics, 1999, ISIE
'99. 12-16 July 1999, Vol. 2, pp. 667 – 672.
69. M.J. Corley, R.D. Lorenz, “Rotor position and velocity estimation for a salient-pole
permanent magnet synchronous machine at standstill and high speeds”, IEEE
Transactions on Industry Applications, July-Aug. 1998, Vol. 34, pp. 784 – 789.
70. Z. Chen, M. Tomita, S. Ichikawa, S. Doki, S. Okuma, “Sensorless control of interior
permanent magnet synchronous motor by estimation of an extended electromotive
force”; Industry Applications Conference, 2000, Vol. 3, 8-12 Oct. 2000 pp. 1814 –
1819.
71. Z. Chen, M. Tomita, S. Doki, S. Okuma, “An extended electromotive force model for
sensorless control of interior permanent-magnet synchronous motors”, IEEE
Transactions on Industrial Electronics, Vol. 50, Issue 2, April 2003, pp. 288 – 295.
72. Z. Chen, M. Tomita, S. Doki, S. Okuma, “New adaptive sliding observers for
position- and velocity-sensorless controls of brushless DC motors”, IEEE Transactions
on Industrial Electronics, Vol. 47, Issue 3, June 2000, pp. 582 – 591.
73. R. Dhaouadi, N. Mohan, L. Norum, “Design and implementation of an extended
Kalman filter for the state estimation of a permanent magnet synchronous motor”, IEEE
Transactions on Power Electronics, Vol. 6, Issue 3, July 1991, pp. 491 – 497.
176

References
74. N. Ertugrul, P. Acarnley, “A new algorithm for sensorless operation of permanent
magnet motors”; IEEE Transactions on Industry Applications, Jan.-Feb. 1994, Vol. 30,
pp. 126 – 133.
75. C. French, P. Acarnley, I. Al-Bahadly, “Sensorless position control of permanent
magnet drives”; Industry Applications Conference, IAS '95, 8-12 Oct. 1995, Vol. 1, pp.
61 – 68.
76. M. Fu, L. Xu, “A novel sensorless control technique for permanent magnet synchronous
motor (PMSM) using digital signal processor (DSP)”, Proceedings of the IEEE 1997
National Aerospace and Electronics Conference, 14-17 July 1997, NAECON 1997,Vol.
1, pp. 403 – 408.
77. J. Hu, B. Wu, “New integration algorithms for estimating motor flux over a wide speed
range”, IEEE Transactions on Power Electronics, , Sept. 1998, Vol. 13, pp. 969 – 977.
78. Y. Jeong, R.D. Lorenz, T.M. Jahns, S. Sul, “Initial rotor position estimation of an
interior permanent magnet synchronous machine using carrier-frequency injection
methods”, IEEE International Electric Machines and Drives Conference, IEMDC'03, 1-
4 June 2003, Vol. 2,pp. 1218 – 1223.
79. L.A. Jones, J.H. Lang, “A state observer for the permanent-magnet synchronous
motor”, IEEE Transactions on Industrial Electronics, Vol. 36, Issue 3, Aug. 1989 pp.
374 – 382.
80. J.P. Johnson, M. Ehsani, Y. Guzelgunler, “Review of sensorless methods for
brushless DC”, Industry Applications Conference, 1999. Thirty-Fourth IAS Annual
Meeting. Conference Record of the 1999 IEEE, 3-7 Oct. 1999, Vol. 1, pp. 143 – 150.
81. J. S. Kim; S. K. Sul, “New stand-still position detection strategy for PMSM drive
without rotational transducers”, Applied Power Electronics Conference and Exposition,
APEC '94, 13-17 Feb. 1994, Vol.1, pp. 363 – 369.
82. J.-S. Kim, S.-K. Sul, “New approach for high-performance PMSM drives without
rotational position sensors”, IEEE Transactions on Power Electronics, Sept. 1997, Vol.
12, pp. 904 – 911.
83. T.-S. Low, T.-H. Lee, K.-T. Chang, “A nonlinear speed observer for permanentmagnet
synchronous motors”, IEEE Transactions on Industrial Electronics, Vol. 40,
Issue 3, June 1993, pp. 307 – 316.
84. M. Linke, R. Kennel, J. Holtz, “Sensorless position control of permanent magnet
synchronous machines without limitation at zero speed”, IECON 02, 5-8 Nov. 2002,
Vol. 1, pp. 674 – 679.
85. N. Matsui, “Sensorless PM brushless DC motor drives”, Transactions on Industrial
Electronics, IEEE , Vol. 43, April 1996, pp. 300 – 308.
86. N. Matsui, M. Shigyo, “Brushless DC motor control without position and speed
sensors”; IEEE Transactions on Industry Applications, Jan.-Feb. 1992, Vol. 28 pp. 120
– 127.
87. N. Matsui, “Sensorless PM brushless DC motor drives”, IEEE Transactions on
Industrial Electronics, April 1996, Vol. 43, pp. 300 – 308.
88. S. Nakashima, Y. Inagaki, I. Miki, “Sensorless initial rotor position estimation of
surface permanent-magnet synchronous motor”, IEEE Transactions on Industry
Applications, Nov.-Dec. 2000, Vol. 36, pp. 1598 – 1603.
89. T. Noguchi, K. Yamada, S. Kondo, I. Takahashi, “Rotor position estimation method of
sensorless PM motor at rest with no sensitivity to armature resistance”, Industrial
Electronics, Control, and Instrumentation, Proceedings of the 1996 IEEE IECON 1996,
5-10 Aug. 1996,Vol. 2, pp. 1171 – 1176.
90. S. Ostlund, M. Brokemper, “Initial rotor position detections for an integrated PM
synchronous motor drive”, Industry Applications Conference, IAS '95, 8-12 Oct. 1995,
Vol. 1, pp. 741 – 747.
91. G. Pesse, T. Paga, “A permanent magnet synchronous motor flux control scheme
without position sensor”, Conference proceedings on the EPE’97, Trondheim 1997,
Vol. 4, pp. 553 – 556.
177

References
92. R.B. Sepe, J.H. Lang, “Real-time observer-based (adaptive) control of a permanentmagnet
synchronous motor without mechanical sensors”, IEEE Transactions on
Industry Applications, Vol. 28, Issue 6, Nov.-Dec. 1992, pp. 1345 – 1352.
93. P.B. Schmidt, M.L. Gasperi, G. Ray, A.H. Wijenayake, “Initial rotor angle detection
of a nonsalient pole permanent magnet synchronous machine”, Industry Applications
Conference, IAS '97., 5-9 Oct. 1997, Vol. 1, pp. 459 – 463.
94. J.X. Shen, Z.Q. Zhu, D. Howe, “Improved speed estimation in sensorless PM
brushless AC drives”, Transactions on Industry Applications, July-Aug. 2002, Vol. 38,
pp. 1072 – 1080.
95. Yen-Shin Lai; Fu-San Shyu; Shian Shau Tseng, “New initial position detection
technique for three-phase brushless DC motor without position and current sensors”,
Transactions on Industry Applications, , March-April 2003, Vol. 39, pp. 485 – 491.
96. J. Solsona, M.I. Valla, C. Muravchik, “A nonlinear reduced order observer for
permanent magnet synchronous motors”, IEEE Transactions on Industrial Electronics,
Aug. 1996, Vol. 43, pp. 492 - 497
97. J.A. Solsona, M.I. Valla, “Disturbance and nonlinear Luenberger observers for
estimating mechanical variables in permanent magnet synchronous motors under
mechanical parameters uncertainties”, IEEE Transactions on Industrial Electronics, Vol.
50, Issue 4, Aug. 2003 pp. 717 – 725.
98. J. Solsona, M.I. Valla, C. Muravchik, “Nonlinear control of a permanent magnet
synchronous motor with disturbance torque estimation”, International Conference on
Industrial Electronics, Control and Instrumentation, 1997. IECON 97, Vol. 1, 9-14 Nov.
1997, pp. 120 – 125.
99. T. Takeshita, N. Matsui, “Sensorless control and initial position estimation of salientpole
brushless DC motor”, Advanced Motion Control, 1996. AMC '96, 18-21 March
1996, Vol. 1, pp. 18 – 23.
100. M. Tomita, T. Senjyu, S. Doki, S. Okuma, “New sensorless control for brushless DC
motors using disturbance observers and adaptive velocity estimations”, IEEE
Transactions on Industrial Electronics, Vol. 45, Issue 2, April 1998, pp. 274 – 282.
101. M. Tursini, R. Petrella, F. Parasiliti, “Initial rotor position estimation method for PM
motors”; IEEE Transactions on Industry Applications, Nov.-Dec. 2003, Vol. 39, pp.
1630 – 1640.
102. R. Wu, G.R. Slemon, “A permanent magnet motor drive without a shaft sensor”,
Industry Applications Society Annual Meeting, 1990., Conference Record of the 1990
IEEE, 7-12 Oct. 1990, Vol.1, pp. 553 – 558.
103. D. Yousfi, M. Azizi, A. Saad, “Sensorless position and speed detection for permanent
magnet synchronous motor”; Power Electronics and Motion Control Conference,
PIEMC 2000, 15-18 Aug. 2000,Vol. 3, pp. 1224 – 1229.
104. D. Yousfi, M. Azizi, A. Saad, “Robust position and speed estimation algorithm for
permanent magnet synchronous drives”, Industry Applications Conference, 8-12 Oct.
2000, Vol.3, pp. 1541 – 1546.
Others
105. H.-B. Shin, “New antiwindup PI controller for variable-speed motor drives”
IEEE Transactions on Industrial Electronics, June 1998, Vol. 45, pp. 445 – 450.
106. MathWorks, Inc, "Matlab® The Language of Technical Computing", Release 12,
2000.
Papers written during work on this thesis
107. D. L. Sobczuk, D. Świerczyński "DSP implementation of vector controlled inverter
fed induction motor with neural network speed estimator" SENE '99 tom. 2 str. 601-
606.
178

References
108. D. Świerczyński, Marian P. Kaźmierkowski, Paweł Z. Grabowski "Nowe algorytmy
estymacji strumienia stojana w silnikach indukcyjnych klatkowych zasilanych z
falowników PWM" SENE '99 tom. 2 str. 649-654.
109. M. Angielczyk, M. Malinowski, D. Świerczyński, „Stanowisko laboratoryjne oparte
o kartę DSP do sterowania prostownikiem PWM oraz silnikiem PMSM zasilanym z
falownika PWM”, Kamień Śląski 2000, Międzynarodowe XII Sympozjum
Mikromaszyny i Serwonapędy, str. 480-487.
110. D. Świerczyński, M. P. Kaźmierkowski.:„Bezpośrednie sterowanie momentu silnika
synchronicznego o magnesach trwałych zasilanym z falownika PWM”, SENE 2001
,tom II, str. 629
111. D. Świerczyński, M. P. Kaźmierkowski."Direct torque control of permanent magnet
synchronous motor (PMSM) using space vector modulation (DTC-SVM)", SME'2002,
pp. 245
112. D. Świerczyński, M. P. Kaźmierkowski, Frede Blaabjerg "DSP Based Direct Torque
Control of Permanent Magnet Synchronous Motor (PMSM) Using Space Vector
Modulation (DTC-SVM)", IEEE-ISIE 2002, L’Aquila, Italy, pp. 723-727.
113. D. Świerczyński, M. Żelechowski "Studia symulacyjne nad universalną strukturą
bezpośredniego sterowania momentem i strumieniem dla silników synchronicznych o
magnesach trwałych oraz silników asynchronicznych ", Modelowania i Symulacja MiS
2002, Zakopane.
114. D. Świerczyński, M. P. Kaźmierkowski, Frede Blaabjerg "Direct Torque Control
with Space Vector Modulation (DTC-SVM) for Permanent Magnet Synchronous
Motor (PMSM)", EPE-PEMC 2002, Dubrovnik & Cavtat, Croatia.
115. D. Świerczyński, M. Żelechowski "Universalną strukturą bezpośredniego sterowania
momentem i strumieniem dla silników synchronicznych o magnesach trwałych oraz
silników asynchronicznych ", Mikromaszyny i Serwonapędy MiS 2002, Krasiczyn.
116. D. Świerczyński, M. P. Kaźmierkowski " Direct Torque Control of Permanent
Magnet Synchronous Motor (PMSM) Using Space Vector Modulation (DTC-SVM) –
Simulation and Experimental Results", IECON 2002, Sevilla, Spain.
117. D.Świerczyński, M. Żelechowski „Bezpośrednie sterowanie momentem i
strumieniem dla silników synchronicznych i asynchronicznych”, II Krajowa
Konferencja MiS-2 Modelowanie i Symulacja, Kościelisko, 24-28 czerwca 2002, pp.
187-194.
118. D. Świerczyński, M. Żelechowski ”Universal structure of direct torque control for
AC motor drives”, III Summer Seminar on Nordick Network for Multi Disciplinary
Optimised Electric Drives, June 21-23, 2003, Zegrze, Poland
119. M. Jasiński D. Świerczyński, M. P. Kazmierkowski, M. Malinowski ”Direct Control
of an AC-DC-AC converter ”, III Summer Seminar on Nordick Network for Multi
Disciplinary Optimised Electric Drives, June 21-23, 2003, Zegrze, Poland
120. M. Żelechowski, D. Świerczyński, M. P. Kazmierkowski, J. Załęski „Uniwersalne
napędy falownikowe ze sterowaniem mikroporcesorowym nowej generacji”,
Elektro.info Nr 6 (17) 2003, str. 26-28
121. D. Świerczyński, M. Żelechowski, „Universal Structure of Direct Torque Control for
AC Motor Drives”, V Power Electronics Seminar, 7-8 July, 2003, Berlin, Niemcy
122. M. Jasinski, D. Swierczynski, M. P. Kazmierkowski "Novel Sensorless Direct Power
and Torque Control of Space Vector Modulated AC/DC/AC Converter", In Proc. of the
ISIE 2004 Conf., 4-7 May, Ajaccio, France, pp.1147- 1152.
123. D. Świerczyński, M. Żelechowski ”Universal structure of direct torque control for
AC motor drives”, Przegląd Elektrotechniczny, May, 2004, Poland, pp. 489-492.
124. M. P. Kaźmierkowski , M. Żelechowski, D. Świerczyński „Simple DTC-SVM
Control for Induction and PM Synchronous Motor”, In Proc. of the ICEM 2004 Conf.,
5-8 September, Krakow, Poland, on CD.
125. M. P. Kaźmierkowski , M. Żelechowski, D. Świerczyński „DTC-SVM an Efficient
Method for Control Both Induction and PM Synchronous Motor”, In Proc. of the EPE-
PEMC 2004 Conf., 2-4 September, Riga, Latvia, on CD.
179

References
126. B. W. Mehdawi, D. Świerczyński, M. P. Kaźmierkowski „Design and Investigation
of Vector Control Scheme PWM Inverter-fed Nonsalient Permanent Magnet
Synchronous Motor”, In Proc. of the PPEE 2005 Conf., 2-5 April, Poland, Wisla, on
CD.
180