VECTOR CONTROL DRIVE OF PERMANENT

MAGNET SYNCHRONOUS MOTOR USING

MATLAB/SIMULINK
A project report submitted
in partial fulfillment of the requirement for the award of the Degree

BACHELOR OF TECHNOLOGY
IN
ELECTRICAL AND ELECTRONICS ENGINEERING
BY
Ms. R. RADHIKA (07BE1A0213)
Ms. B. USHA RANI (07BE1A0220)

Under the esteemed guidance of

Mr. H. SUDHEER
Associate Professor,
KITE
GHATKESAR

Department of Electrical and Electronics Engineering

Krishna Murthy Institute of Technologies and Engineering
(Affiliated to J.N.T.U.H, approved by AICTE)
Edulabad (v), Ghatkesar (M), RR dist- 501301
2010-2011

1
 KRISHNA MURTHY COLLEGE OF TECHNOLOGY AND ENGINEERING

DEPARTMENT OF ELECTRICAL AND

ELECTRONICS ENGINEERING

CERTIFICATE
This is to certify that the project entitled “VECTOR CONTROL
DRIVE OF PERMANENT MAGNET SYNCHRONOUS MOTOR
USING MATLAB/SIMULINK” that is being submitted by
Ms. R. Radhika (07BE1A0213) Ms. B. Usha Rani (07BE1A0220) in
partial fulfillment of the requirement for the award of Degree of Bachelor
of Technology in “Electrical and Electronics Engineering” to
JAWAHARLAL NEHRU TECHNOLOGICAL UNIVERSITY is a
record of bonafide work carried out by them under the guidance and
supervision. The results embodied in this project have not been submitted
to any other university or institute for the award of any degree or
diploma.

Internal Guide Head Of The Department

Mr. H. SUDHEER Mr. B. SRINIVAS

Associate Professor, K.I.T.E

EXTERNAL EXAMINER

i
 ACKNOWLEDGEMENT
It gives immense pleasure to acknowledge the perennial inspiration of

principal DR. S. F. KODAD of KRISHNA MURTHY INSTITUTE OF

TECHNOLOGY AND ENGINEERING and Prof. B.SRINIVAS, Head of

department of Electrical and Electronics for the kind co-operation and

encouragement in bringing out this task. This rightfully belonged to him to facilitate

the completion of this academic task.

I am grateful to H.SUDHEER, Associate professor, Department of Electrical

and Electronics Engineering, for taking me as his student and providing me the

opportunity to work in the area of Power electronics. I am eternally indebted to him

for the facilities of learning that he has provided to me. I am also thankful to him for

his valuable guidance during the course of my research work. I have learnt a lot of

lessons in work ethics, professional behavior and meticulous approach to problem

solving from him which will inspire me for the rest of my life.

One of the reasons for the time I spent in the power electronics lab as a

student being enjoyable is the people. The technical and non-technical discussions I

had with my friends here will be cherished by me for a very long time.

Lastly, I thank my parents for having given me the strength and moral

support during the course of my stay here. I am indebted to them eternally for all

that they have done for me.

R.RADHIKA
B.USHARANI

ii
 ABSTRACT

Permanent magnet synchronous motor (PMSM) has been widely used in high

performance drive applications for its advantages such as compactness, high

efficiency, reliability and suitability to environment. Due to its high power density

and smaller size, PMSM has evolved as the preferred solution for speed and position

control drives on machine tools and robots. A Permanent Magnet Synchronous Motor

(PMSM) is a motor that uses permanent magnets to produce the air gap magnetic field

rather than using electromagnets. These motors have significant advantages, attracting

the interest of researchers and industry for use in many applications. Permanent

magnet synchronous motors are widely used in low and mid power applications such

as computer peripheral equipments, robotics, adjustable speed drives and electric

vehicles.

In order to overcome the inherent coupling effect and the sluggish response of

scalar control the vector control is employed. By using the vector control, the

performance of the AC machine can be made similar to that of a DC machine. In this

work to achieve high performance the vector control of the Permanent magnet

synchronous motors drive is employed. The simulation of PMSM is developed using

SIMULINK. The effectiveness of the proposed control method is verified by

simulation based on MATLAB. The simulation results along with the case study is

presented and explained in detail.

i
CONTENTS
ACKNOWLEDGEMENT
ABSTRACT i
LIST OF FIGURES v-
vi
SYMBOLS vii
1 INTRODUCTION 1-3
1.1 Introduction 2
1.2 Objective 2
1.3 Research methodology 3
1.4 Scope of project 3

2 LITERATURE SURVEY 4-10

2.1 Introduction 5
2.2 Literature survey 5
2.3 Literature conclusion 10
2.4 Problem definition 10

3 THEORITICAL ANALYSES 11-18

3.1 Introduction 12
3.2 Basic constructional details of PMSM 12
3.3 Principle of operation of PMSM 15
3.4 Scalar control 17
3.5 Vector control 17

4 DESCRIPTION OF THE DRIVE SYSTEM 19-28

4.1Permanent Magnet Synchronous Motor Drive System 20
4.2 Permanent Magnet Synchronous Motor 21
4.2.1 Permanent Magnet Materials 21
4.2.2 Classification of Permanent Magnet Motors 22
4.2.2.1 Direction of field flux 22
4.2.2.2 Flux density distribution 23
4.2.2.3 Permanent magnet radial field motors 23
4.3 Position Sensor 25
4.3.1 Optical Encoders 25
4.3.1.1 Incremental encoders 26
4.3.1.2 Absolute encoders 26
4.3.2 Position resolver 27

5 MODELING OF PM DRIVE SYSTEM 29-39

5.1 Introduction 30
5.2 Detailed Modeling of PMSM 30

ii
 5.2.1 Parks Transformation and Dynamic d q Modeling 32
5.2.2 Equivalent Circuit of PMSM 32
5.3 PM Motor Control  32
5.3.1 Field Oriented Control of PM Motors 33
5.3.1.1 Constant torque operation 35
5.3.1.2 Flux-weakening 35
5.4 Speed Control of PM Motor 37
5.4.1 Implementation of the Speed Control Loop 37

6 SPEED CHECK IMPLEMENTATION OF VECTOR CONTROLLED PMSM

DRIVE IN MATLAB/SIMULIMK 40-49
6.1 Introduction 41
6.2 Advantages of SIMULINK 41
6.3 Simulation Tools   42
6.4 SIMULINK Simulation of PMSM Drive 42
6.4.1 Vabc to Vdqo block and Idqo to Iabc block 43
6.4.2 d-axis circuit and q-axis circuit 44
6.4.3 Load Torque Block 45
6.4.4 Speed Block 45
6.4.5 Vector Control Block 46
6.4.6 PM Motor Drive System in SIMULINK 47
6.4.7 Vector control of PMSM drive 48

7 SIMULATION RESULTS 50-68

7.1 Case 1 51-54
7.1.1 Reference torque and actual torque 51
7.1.2 Actual speed and Reference speed 51
7.1.3 Iabc currents and Vabc voltages 52
7.1.4 Vafa and Vbeta 53
7.1.5 Error Signal between Reference Speed and Actual Speed 54
7.1.6 Vd, Vq voltage and Id-ref, Iq-ref 54
7.2 Case 2 56-59
7.2.1 Actual speed and Reference speed 56
7.2.2 Reference torque and actual torque 56
7.2.3 Id, Iq , Id-ref and Iq-ref 57
7.2.4 Iabc currents and Vabc voltages 58
7.2.5 Error Signal between Reference Speed and Actual Speed 59
7.2.6 Vd voltage and Vq voltage 59
7.3case 3 60-63
7.3.1 Actual speed and Reference speed 60
7.3.2 Reference torque and actual torque 61
7.3.3 Iabc currents and Vabc voltages 61
7.3.4 Error Signal between Reference Speed and Actual Speed 62
7.3.5 Vd voltage and Vq voltage 62

iii
 7.3.6 Id, Iq , Id-ref and Iq-ref 63

7.4 Case 4 64-67

7.4.1 Actual speed and Reference speed 64
7.4.2 Reference torque and actual torque 64
7.4.3 Error Signal between Reference Speed and Actual Speed 65
7.4.4 Iq, Id, Id_ref and Iq_ref 66
7.4.5 Vd voltages and Vq voltages 67
7.5 Discussion of results 68

8 CONCLUSION AND FUTURE SCOPE 69

8.1 Conclusion 70
8.2 Future scope 70

Appendices 71-76

Bibliography 77

iv
LIST OF FIGURES

NAME OF THE FIGURE PAGE NO.

Figure 3.1 Three basic configurations of PMSMs 13

Figure 3.2 Torque Establishment 16
Figure 3.3 (a) Rotor (b) Stator (c) Phase Currents 16
Figure 3.4 (a) Separately Excited Dc Motor (b) Vector Controlled Ac Motor 18
Figure 4.1 Drive Systems Schematic 20
Figure 4.2 Permanent Magnets Synchronous Motor 21
Figure 4.3 Flux Density Vs Magnetizing Field of Permanent Magnetic
Materials 22
Figure 4.4 Surface Permanent Magnet Motor 24
Figure 4.5 Interior Permanent Magnet Motor 25
Figure 4.6 Optical Encoder 25
Figure 4.7 Quadrature Encoder Channels 26
Figure 4.8 Absolute Encoder 27
Figure 4.9 Resolver 27
Figure 4.10 Excitation and Output Signal of the Resolver 28
Figure 5.1 Motor Axis 30
Figure 5.2 PMM Electric Circuits without Damper Windings 32
Figure 5.3 Self Controls Synchronous Motor 33
Figure 5.4 Steady State Torque versus Speed 33
Figure 5.5 Block Diagram 38
Figure 5.6 Proportional integral Controller 38
Figure 5.7 Block Diagram of Speed Loop 39
Figure 6.1 Vabc to Vdqo block 43
Figure 6.2 Idqo to Iabc block 43
Figure 6.3 d-axix circuit and q-axis circuit 44
Figure 6.4 Load Torque Block 45
Figure 6.5 Speed Block 45
Figure 6.6 Vector Control Block 46
Figure 6.8 PM Motor Drive System in SIMULINK 48
Figure 6.9 Vector control of PMSM drive 48
Figure 7.1 Reference torque and actual torque 51
Figure 7.2 (a) Actual speed and (b) Reference speed 52
Figure 7.3 (a) I abc currents and (b) Vabc voltages 53
Figure 7.4 (a)Vafa and (b)Vbeta 53
Figure 7.5 Error Signal between Reference Speed and Actual Speed 54
Figure 7.6 Vd voltage and Vq voltage 55
Figure 7.7 (a) Id-ref and (b) Iq-ref 55
Figure 7.8 (a) Actual speed and (b) Reference speed 56
Figure 7.9 (a) Reference torque and (b) actual torque 57
Figure 7.10 (a) Id current and (b) Iq current 57
Figure 7.11 (a) Id-ref and (b) Iq-ref 58
Figure 7.12 (a)Iabc currents and (b)Vabc voltages 58
Figure 7.13 Error Signal between Reference Speed and Actual Speed 59
Figure 7.14 Vd voltage and Vq voltage 60

v
Figure 7.15 (a) Actual speed and (b) Reference speed 60
Figure 7.16 (a) Reference torque and (b)actual torque 61
Figure 7.17 (a) Iabc currents and (b)Vabc voltages 61
Figure 7.18 Error Signal between Reference Speed and Actual Speed 62
Figure 7.19 (a) Vd voltage and (b) Vq voltage 62
Figure 7.20 (a) Id current and (b) Iq current 63
Figure 7.21 (a) Id-ref and (b) Iq-ref 63
Figure 7.22 (a) Actual speed and (b) Reference speed 64
Figure 7.23 (a) Reference torque and (b) actual torque 65
Figure 7.24 Error Signal between Reference Speed and Actual Speed 65
Figure 7.25 (a) Id current and (b) Iq current 66
Figure 7.26 (a) Id-ref and (b) Iq-ref 66
Figure 7.27 (a) Vd Voltage and (b) Vq Voltage 67

vi
Symbols

iabc abc phase current vector

idrr,iqrr Fictitious rotor currents in d and q axis, A
ids,iqs d and q axis stator currents, A
iα,iβ Two phase instantaneous currents, A
iαm,iβm Currents in the ne rotor frames, A
J Total moment of inertia, Kg-m2
Lq, Ld Quadrature and direct axis self inductances, H
Lqn, Ldn Normalized quadrature and direct axis self inductances, p.u
Lqq, Ldd Self inductance of the stator q and d axis windings, H
Lαα Self inductance of α the rotor axis windings, H
Lββ Self inductance of β the rotor axis windings, H
Rd, Rq Stator d and q axis winding resistance, Ω
P Differential operator
Tabc Transformation from abc to qd0 axes
Vabc abc voltage vector

Subscripts
d Direct-axis
q Quadrature-axis

Acronyms
PM Permanent Magnet
PMDC Permanent Magnet Direct Current
PMSM Permanent Magnet Synchronous Machine
PWM Pulse Width Modulation
MMF Magneto-Motive Force
EMF Electro-Motive Force
BLDC BrushLess Direct Current

vii
 CHAPTER I

INTRODUCTION

1
 1 INTRODUCTION

1.1 INTRODUCTION

With the advent of switching power transistor and silicon controlled rectifier
devices in later part of 1950’s, and the replacement of the mechanical commutator
with an electronic commutator in the form of an inverter was achieved. These two
developments have contributed to the development of the PM synchronous and
brushless DC machines. A Permanent Magnet Synchronous Motor (PMSM) is a
motor that uses permanent magnets to produce the air gap magnetic field rather than
using electromagnets. These motors have significant advantages, attracting the
interest of researchers and industry for use in many applications.

Permanent magnet synchronous motors are widely used in low and mid power
applications such as computer peripheral equipments, robotics, adjustable speed
drives and electric vehicles.

The growth in the market of PMSM motor drives has demanded the need of
simulation tool capable of handling motor drive simulations. Simulations have helped
the process of developing new systems including motor drives, by reducing cost and
time. Simulation tools have the capabilities of performing dynamic simulations of
motor drives in a visual environment so as to facilitate the development of new
systems.

In this work, the simulation of PMSM is developed using SIMULINK. The

vector control is one of the high performance control strategies for ac machine. The
aim of the project is to study the implementation of the vector control in Permanent
Magnet Synchronous Motor (PMSM).

1.2 OBJECTIVES

The objectives of the project are

i) To stimulate the vector control of permanent magnet synchronous motor.

ii) To analyze the simulation results.

2
1.3 RESEARCH METHODOLOGY

The research work is undertaken in the following stages:

i) Studied the application of MATLAB/SIMULINK.

ii) Studied the theoretical basis of the vector control for permanent magnet
synchronous motor drives.
iii) Simulation of vector control of permanent magnet synchronous motor is
performed using SIMULINK.
iv) Analyzed the simulation results.

1.4 SCOPE OF PROJECT

The scope of work for this project

i) PMSM with saliency is considered.

ii) Simulation is performed using MATLAB/SIMULINK.
iii) The performance of vector control of PMSM is discussed based on the
simulation results.

3
 CHAPTER-II

LITERATURE SURVEY AND PROBLEM

DEFINITION

4
 2 LITERATURE SURVEY AND PROBLEM

DEFINITION

2.1 INTRODUCTION

A literature survey forms the basis on which a project can be built or

developed. It forms the core to which ideas can be added and developed into a
comprehensive system, which will be able to cover the deficiencies of some of the
existing systems.

This chapter deals with the data and information accumulated after referring to
many books, articles and technical papers written by well-known authors and the
problem definition of the project.

2.2 LITERATURE SURVEY

[1] T. Sebastian, G. Slemon, and M. Rahman, "Modeling of permanent magnet

synchronous motors," Magnetics, IEEE Transactions on, vol. 22, pp. 1069-1071,

1986.

[2] T. M. Jahns, G. B.Kliman, and T. W. Neumann, "Interior Permanent-Magnet

Synchronous Motors for Adjustable-Speed Drives," Industrial Applications, IEEE

Transactions on, vol. IA-22, pp. 738-746, 1986.

[3] P. Pillay and R. Krishnan, "Modeling of permanent magnet motor drives,"

Industrial Electronics, IEEE Transactions on, vol. 35, pp. 537-541, 1988.

[4] P. Pillay and R. Krishnan, "Modeling, simulation, and analysis of permanent-

magnet motor drives. I. The permanent-magnet synchronous motor drive," Industry

Applications, IEEE Transactions on, vol. 25, pp. 265-273, 1989.

[5] B. K. Bose, Modern power electronics and AC drives: Prentice Hall, 2002

[6] A. H. Wijenayake and P. B. Schmidt, "Modeling and analysis of permanent

magnet synchronous motor by taking saturation and core loss into account," 1997.
5
[7] K. Jang-Mok and S. Seung-Ki, "Speed control of interior permanent magnet

synchronous motor drive for the flux weakening operation," Industry Applications,

IEEE Transactions on, vol. 33, pp. 43-48, 1997.

[8] Weera Kaewjind and Mongkol Konghirun “Vector Control Drive of Permanent

Magnet Synchronous Motor Using Resolver Sensor” ECTI transactions on electrical

eng., electronics, and communications vol.5, no.1 february 2007.

PM motor drives have been a topic of interest for the last twenty years.
Different authors have carried out modeling and simulation of such drives. Some of
them have been discussed in detail.

[1] T. Sebastian, G. Slemon, and M. Rahman, "Modeling of permanent magnet

synchronous motors"

In 1986 Sebastian, T., Slemon, G. R. and Rahman, M. A. [1] reviewed

permanent magnet synchronous motor advancements and presented equivalent
electric circuit models for such motors and compared computed parameters with
measured parameters. Experimental results on laboratory motors were also given.

[2] T. M. Jahns, G. B.Kliman, and T. W. Neumann, "Interior Permanent-

Magnet Synchronous Motors for Adjustable-Speed Drives,"

In 1986 Jahns, T.M., Kliman, G.B. and Neumann, T.W. [2] discussed that
interior permanent magnet (IPM) synchronous motors possessed special features for
adjustable speed operation which distinguished them from other classes of ac
machines. They were robust high power density machines capable of operating at
high motor and inverter efficiencies over wide speed ranges, including considerable
range of constant power operation. The magnet cost was minimized by the low
magnet weight requirements of the IPM design. The impact of the buried magnet
configuration on the motor’s electromagnetic characteristics was discussed. The rotor
magnetic saliency preferentially increased the quadrature-axis inductance and
introduced a reluctance torque term into the IPM motor’s torque equation. The
electrical excitation requirements for the IPM synchronous motor were also discussed.

6
The control of the sinusoidal phase currents in magnitude and phase angle with
respect to the rotor orientation provided a means for achieving smooth responsive
torque control. A basic feed forward algorithm for executing this type of current
vector torque control was discussed, including the implications of current regulator
saturation at high speeds. The key results were illustrated using a combination of
simulation and prototype IPM drive measurements.

[3] “Modeling, Simulation, And Analysis of Permanent-Magnet Synchronous

Motor Drive” by P. Pillay and R. krishnan

In 1988 Pillay and Krishnan, R.[3] presented PM motor drives and classified
them into two types such as permanent magnet synchronous motor drives (PMSM)
and brushless DC motor (BDCM) drives. The PMSM has a sinusoidal back EMF and
requires sinusoidal stator currents to produce constant torque while the BDCM has a
trapezoidal back EMF and requires rectangular stator currents to produce constant
torque. The PMSM is very similar to the wound rotor synchronous machine except
that the PMSM that is used for servo applications tends not to have any damper
windings and excitation is provided by a permanent magnet instead of a field winding.
Hence the d, q model of the PMSM can be derived from the well known model of the
synchronous machine with the equations of the damper windings and field current
dynamics removed. Equations of the PMSM are derived in rotor reference frame and
the equivalent circuit is presented without dampers. The damper windings are not
considered because the motor is designed to operate in a drive system with field-
oriented control. Because of the non-sinusoidal variation of the mutual inductances
between the stator and rotor in the BDCM, it is also shown in this paper that no
particular advantage exists in transforming the abc equations of the BCDM to the d, q
frame. As an extension of his previous work, Pillay, P. and Krishnan, R. in 1989
presented the permanent magnet synchronous motor (PMSM) which was one of
several types of permanent magnet ac motor drives available in the drives industry.
The motor had a sinusoidal flux distribution. The application of vector control as well
as complete modeling, simulation, and analysis of the drive system were given. State
space models of the motor and speed controller and real time models of the inverter
switches and vector controller were included. The machine model was derived for the
PMSM from the wound rotor synchronous motor. All the equations were derived in

7
rotor reference frame and the equivalent circuit was presented without dampers. The
damper windings were not considered because the motor was designed to operate in a
drive system with field-oriented control. Performance differences due to the use of
pulse width modulation (PWM) and hysteresis current controllers were examined.
Particular attention was paid to the motor torque pulsations and speed response and
experimental verification of the drive performance were given.

[4] P. Pillay and R. Krishnan, "Modeling, Simulation, And Analysis Of

Permanent-Magnet Motor Drives".

As an extension of his previous work, Pillay, P. and Krishnan, R. in 1989 [4]

presented the permanent magnet synchronous motor (PMSM) which was one of
several types of permanent magnet ac motor drives available in the drives industry.
The motor had a sinusoidal flux distribution. The application of vector control as well
as complete modeling, simulation, and analysis of the drive system were given. State
space models of the motor and speed controller and real time models of the inverter
switches and vector controller were included. The machine model was derived for the
PMSM from the wound rotor synchronous motor. All the equations were derived in
rotor reference frame and the equivalent circuit was presented without dampers. The
damper windings were not considered because the motor was designed to operate in a
drive system with field-oriented control. Performance differences due to the use of
pulse width modulation (PWM) and hysteresis current controllers were examined.
Particular attention was paid to the motor torque pulsations and speed response and
experimental verification of the drive performance were given.

[5] “Modern Power Electronics And Ac Drives” by B. K. Bose

Bose, B. K., in 2001 [5], presented different types of synchronous motors and
compared them to induction motors. The modeling of PM motor was derived from the
model of salient pole synchronous motor. All the equations were derived in
synchronously rotating reference frame and was presented in the matrix form. The
equivalent circuit was presented with damper windings and the permanent magnet
was represented as a constant current source. Some discussions on vector control
using voltage fed inverter were given.

8
[6] “Modeling And Analysis Of Pmsm” by Wijenayake, A.H. and Schmidt, P.B.

The paper in 1997 by Wijenayake, A.H. and Schmidt, P.B. [6], described the
development of a two-axis circuit model for permanent magnet synchronous motor
(PMSM) by taking machine magnetic parameter variations and core loss into account.
The circuit model was applied to both surface mounted magnet and interior permanent
magnet rotor configurations. A method for on-line parameter identification scheme
based on no-load parameters and saturation level, to improve the model, was
discussed in detail. Test schemes to measure the equivalent circuit parameters, and to
calculate saturation constants which govern the parameter variations were also
presented.

[7] K. Jang-Mok and S. Seung-Ki, "Speed control of interior permanent magnet

synchronous motor drive for the flux weakening operation,"

In 1997 Jang-Mok, K. and Seung-Ki, S. [7], proposed a novel flux-weakening

scheme for an Interior Permanent Magnet Synchronous Motor (IPMSM). It was
implemented based on the output of the synchronous PI current regulator reference
voltage to PWM inverter. The on-set of flux weakening and the level of the flux were
adjusted inherently by the outer voltage regulation loop to prevent the saturation of
the current regulator. Attractive features of this flux weakening scheme included no
dependency on the machine parameters, the guarantee of current regulation at any
operating condition, and smooth and fast transition into and out of the flux weakening
mode. Experimental results at various operating conditions including the case of
detuned parameters were presented to verify the feasibility
of the proposed control scheme.

[8] Weera Kaewjind and Mongkol Konghirun “Vector Control Drive of

Permanent Magnet Synchronous Motor Using Resolver Sensor”

The rotor position is necessary to achieve the vector control drive system of
Permanent Magnet Synchronous Motor (PMSM). In this paper, the resolver sensor
detecting the rotor position of PMSM is fo- cused. The outstanding features of this
sensor are its robust structure and noise insensitivity. The resolver algorithm is
proposed and implemented in the vector control drive system of PMSM. The pro

9
posed scheme has been verified by both simulation and experiment using
MATLAB/Simulink and the TMS320F2812 based digital controller, respectively . The
proposed resolver algorithm has been verified in the current controlled drive system
of PMSM. Both simulation and experimental results are presented. According to these
results, the re-solver algorithm can force the angle error to zero. Thus, the computed
angle can eventually match with the actual rotor angle. Then, the correct rotor speed
computation is guaranteed. In the future works, this algorithm will be extensively
tested in the speed controlled drive system of PMSM.

2.3 LITERATURE REVIEW CONCLUSION

PM motor drives have been a topic of interest for the last twenty years.
Different authors have carried out modeling and simulation of such drives. This thesis
gives a brief note of the special features possessed by the interior PMSM drives
classification, modeling, comparison of induction and synchronous motor drives and
vector control technique of PMSM using resolver sensor.

2.4 PROBLEM DEFINITION

The main objective is to implement the vector control technique in permanent

magnet synchronous motor and to observe the performance of the drive under
different conditions of speed and torque. in order to overcome inherent coupling
effect in scalar control and to increase the performance of the permanent magnet
synchronous motor drive, we are making use of the vector control technique. In order
to save time and money the vector control implementation is carried out in
MATLAB/SIMULINK.

10
 CHAPTER III
THEORITICAL ANALYSIS

11
 3 THEORETICAL ANALYSIS OF PERMANENT
MAGNET SYNCHRONOUS MOTOR

3.1 INTRODUCTION

Permanent magnet synchronous motors are increasing applied in several areas

such as traction, automobiles, robotics and aerospace technology. The power density
of permanent magnet synchronous motor is higher than one of induction motor with
the same ratings due to the no stator power dedicated to the magnetic field production.
Nowadays, permanent magnet synchronous motor is designed not only to be more
powerful but also with lower mass and lower moment of inertia.

3.2 BASIC CONSTRUCTION DETAILS OF PERMANENT

MAGNET SYNCHRONOUS MOTOR

A PMSM consists of a magnetic rotor and wound stator construction. Its

wound stators can rapidly dissipate heat to the motor housing and environment. In
contrast, a brush motor traps the heat under a non-conductive air gap, resulting in
greater efficiency and power density for the PMSM design and providing high torque-
to-inertia ratios. A PMSM motor generates magnetic flux using permanent magnets in
the rotors, which are driven by the stators applying a synchronous rotational field. On
the other hand, the flux that is applied by the stators (the armature-reaction flux)
generates torque most effectively when it is perpendicular to flux generated by the
rotors. To maintain near-perpendicularity between stator flux and rotor flux, two
control methods with position-speed feedback loop are popularly used for controlling
a PMSM: Field-Oriented Control and Brushless DC Control.
In principle, the construction of a permanent magnet synchronous machine
does not differ from that of the BLDC, although distributed windings are more often
used. However, while the excitation current waveform was rectangular with a BLDC,
sinusoidal excitation is used with PMSMs, which eliminates the torque ripple caused
by the commutation. PMSMs are typically fed by voltage source inverters, which
cause time-dependent harmonics on the air gap flux. Permanent magnet synchronous
machines can be realized with either embedded or surface magnets on the rotor, and

12
the location of the magnets can have a significant effect on the motor’s mechanical
and electrical characteristics, especially on the inductances of the machine. As the
relative permeability of the modern rare-earth magnets, such as the NdFeB is only
slightly above unity, the effective air gap becomes long with a surface magnet
construction. This makes the direct-axis inductance very low, which has a substantial
effect on the machine’s overloading capability, and also on the field weakening
characteristics. As the pull-out torque is inversely proportional to the d-axis
inductance, the pull-out torque becomes very high. Typically, the per-unit values of
the d-axis synchronous inductances of the SMPMSM servos vary between 0.2−0.35
p.u., and consequently the pull-out torque is in the range of 4−6 p.u., which makes
them well suitable in motion control applications. The drawback of a low Ld –value is
the very short field weakening range, as the armature reaction with a surface magnet
construction is very weak. This means that a high demagnetizing stator current
component would be required to decrease the air gap flux, and consequently, there
would be very little current left on the q-axis to produce the torque. Direct-axis
inductance of a machine having embedded magnets becomes high, as the rotor
magnets per pole form a parallel connection for the flux, while with a surface magnet
construction they are connected in series. With equivalent magnets, the rotor
reluctance of the surface-magnet construction is therefore double compared to an
embedded-magnet construction, and the inductance is inversely proportional to the
reluctance. With embedded-magnets, the direct-axis inductance is further increased
because of the higher rotor leakage flux. Three basic configurations of PMSMs are
shown in Fig. 3.2.

Figure 3.1 Three basic configurations of PMSMs

13
The most common PM rotor constructions are
a) Non-salient surface magnet rotor, due to high d axis reluctance, Ld is low and
consequently the pull-out torque high.
b) Salient pole surface magnet rotor with inset magnets, which is basically the same
as a), but this type produces also some reluctance torque.
c) Embedded magnets in the rotor, which has a high Ld value, and consequently a
poor overloading capability, but a lot better field weakening characteristics than with
the surface magnet constructions.
Typically the construction of the PMSM servomotor is somewhere between (a)
and (b), and the q-axis inductance is larger. Industrial PMSMs often represent the type
(c).
In addition to the good overloading capability, another reason that makes the
surface magnet construction favorable in servo applications is the lower inertia. With
multi-pole machines, the rotor and the stator yokes can be made very thin, and all the
additional iron can be removed from the rotor to provide a lower inertia. These large
holes also improve the heat transfer from the rotor, as the high frequency flux
pulsations generate heat on the magnets and on the rotor iron. As the servomotors
must typically rotate very fast, gluing does usually not suffice in attaching the
magnets on the surface of the rotor, and some non-magnetic material, such as a
stainless steel cylinder or a fibre-glass band must be used to support the magnets. The
problem in using steel is that it is a highly conductive material, and the air gap
harmonics strongly generate losses and consequently heat in it. Therefore, a fibre-
glass band or a plastic cylinder is more often used for the magnet retaining.
Unfortunately, electrical insulators are also thermal insulators, which mean that their
thermal conductivity for the heat generated in the rotor iron and in the magnets is
poor. The temperature rise of the magnets decreases their remanence flux density, and
consequently the torque production.
The rotor in Fig. 3.1 (b) with inset surface magnets has better mechanical
characteristics, but on the other hand, it has higher leakages between two adjacent
magnets. In addition to the higher leakage, the torque production decreases more as
the motor must operate at higher pole angle due to increased q-axis inductance
compared to a non-salient rotor. Typically, the construction of commercial
servomotors is somewhere between (a) and (b) in Fig. 3.1, that is, the magnets are
slightly embedded in the rotor. This improves the mechanical strength of the rotor and
14
introduces a reluctance difference-based term in the torque. According to
measurements made at LUT for eight different commercial servomotors in the power
range of 3−5 kW, the values for the q-axis inductances were 10−20 % higher than the
values in the d-direction.
With buried magnets and flux concentration, a sinusoidal air gap flux density
distribution is possible with simple rectangular magnets. A sinusoidal air gap flux
distribution significantly decreases the cogging torque especially with low-speed
multi-pole machines that have a low number of slots per poles per phase number q.
Also, it is possible to increase the air gap flux density beyond the remanence flux
density of the magnets with a flux concentration arrangement, and the machine can
produce more torque at a given volume. This is especially desirable in low speed
applications, such as in wind generators and in propulsion motors (ABB Azipod®)
where the space is limited. As the direct-axis inductance is typically high with a
buried magnet construction, the overloading capability will be poor, which makes this
motor type incompetent in motion control applications. Typically, the embedded v-
shape magnet machine can have Ld approx. 0.7 p.u, which means only 1.4 p.u.
overloading capability according to the load-angle equation of a synchronous machine
with the assumption that EPM = us = 1 p.u. and Ld = Lq. If there is a reluctance
difference in the machine, the maximum torque can be somewhat larger. It must,
however, be borne in mind that despite the embedded magnets, it is of course possible
to increase the physical air gap large enough, and thereby to decrease the direct axis
inductance of the machine remarkably from the value given above. However, the
consumption of the magnet material is increased remarkably in such a case.

3.3 PRINCIPLE OF OPERATION OF PMSM

The PMSM rotate because of the magnetic attraction between the rotor and the
stator poles. When the rotor poles are facing stator poles of the opposite polarity, a
strong magnetic attraction is set up between them. The mutual attraction locks the
rotor and the stator poles together and the rotor is literally yanked into step with the
revolving stator magnetic field. At no-load conditions, rotor poles are directly
opposite to the stator poles and their axes coincide. At load conditions the rotor poles
lag behind the stator poles, but the rotor continues to turn at synchronous speed, the

15
mechanical angle “θ” between the poles increases progressively as we increase the
load.

Figure 3.2 Torque Establishment

• Permanent magnet rotor (a)

• Three-phase Y-connected stator (b)
• Sinusoidal phase currents (c)
• Each phase is 120º displaced from the others
• Phase currents must sum to 0

Figure 3.3 (a) Rotor (b) Stator (c) Phase Currents

16
3.4 SCALAR CONTROL

Scalar control is based on relationships valid in steady state. It is simple but

due to the inherent coupling effect (i.e., torque and flux are proportional to the voltage
or current and frequency) gives sluggish response and the system can be easily prone
to instability. In this only magnitude and frequency of voltage, current, etc are
controlled. Scalar control is used where several motors are driven in parallel by the
same inverter. In order to overcome these problems we are going for vector control.

3.5 VECTOR CONTROL

The vector or field-oriented control technique brought on a renaissance in

modern high-performance control of ac drives. This control method has found wide
acceptance in applications such as paper mills, textile mills, steel rolling mills,
machine tools, servos, and robotics. With vector or decoupling control, the dynamics
of ac drives is similar to that of dc drives, and with current control, the conventional
stability limit of ac machine does not arise. This is indeed a remarkable
accomplishment. The direct or feedback method, which was developed by Blaschke,
depends on unit vector generation from the machine terminal voltages. As usual,
harmonic noise becomes a problem in feedback signal processing, and the method is
difficult to use near zero speed because of the dominance of stator drop. Of course,
for servo-type applications, the unit vectors can be computed from stator currents and
speed signals. In the indirect or feed forward method, which was developed by Hasse,
the above problems do not exist, but the controller is highly dependent on machine
parameters. This method has gained popularity in industrial applications. At present,
significant R&D efforts have been focused on parameter identification techniques.
The so-called slip gain tuning in order to have decoupling between the rotor flux and
torque component of current has been attempted by reactive power balancing,
injecting a pseudo-random binary sequence, Kalman filter estimation, and MRAC
balancing of reactive power, torque, and voltages. While the on-line controller tuning
with initial parameters is not difficult, tracking of controller parameters with machine
parameters during system operation is always a challenge. Recently, a hybrid or
universal vector control method has been suggested, where the indirect vector control
operates in the lower speed range and is switched to parameter-independent direct
vector control in the higher speed range. It should be mentioned here that, the vector

17
control can be applied to both induction and synchronous machines and, in fact can be
applied to the general AC system for independent active and reactive power control.
By using the Vector control, the performance of the ac machine can be made similar
to that of a separately excited DC motor by the orientation of the stator MMF or
current vector in relation to the rotor flux to achieve a desired objective.

Figure 3.4 (a) Separately Excited Dc Motor (b) Vector Controlled Ac Motor

In vector control we are making the ia and if which are responsible for producing
the fluxes of Ψa and Ψf orthogonal to each other.
In case of DC machine, the construction is such that i a and if are perpendicular
to each other. Ia is responsible for torque and if is responsible for the production of
flux. When torque is controlled by ia the Ψf is not affected. When flux is controlled by
if, there will be no effect on Ψa. Therefore a DC motor has fast transient response.
But due to the inherent coupling effect the ac machine cannot give such fast
response. In order to exhibit the DC machine characteristics, the machine is controlled
in the synchronously rotating reference frame (d e-qe), where the sinusoidal machine
variables appear as DC quantities in the steady state.
The performance of vector control is comparable to DC machine. It produces
less ripples but the system is more complex and less robust when compared to DTC
control.

18
 CHAPTER IV

DESCRIPTION OF THE DRIVE SYSTEM

19
 4 DESCRIPTION OF THE DRIVE SYSTEM

This chapter deals with the description of the different components of the drive
system such as permanent magnet motors, position sensors, inverters and current
controllers of the drive system. A review of permanent magnet materials and
classification of permanent magnet motors is also given.

4.1 PERMANENT MAGNET SYNCHRONOUS MOTOR DRIVE

SYSTEM

The motor drive consists of four main components, the PM motor, inverter,
control unit and the position sensor. The components are connected as shown in
Figure 4.1.

Figure 4.1 Drive System Schematic

PM Motor: It is a motor that uses permanent magnets to produce the air gap magnetic
field. The most commonly used magnetic materials are rare earth magnets such as
NdFeB, SmCo, Strontium Ferrite or Barium ferrite etc.

Inverter: The stator windings of the motor are fed by an inverter that generates a
variable frequency variable voltage.

Position sensor: Instead of controlling the inverter frequency independently, the

frequency and phase of the output wave are controlled using a position sensor. These
are mainly used for determining the position of the rotor. The most commonly used
position sensors are encoders and resolvers. Depending on the application and

20
performance desired by the motor a position sensor with the required accuracy can be
selected.
Control unit: The control input and the rotor position signal is given to the controller
and depending upon both the signals it will generate the output which is given to the
inverter.
4.2 PERMANENT MAGNET SYNCHRONOUS MOTOR

A permanent magnet synchronous motor (PMSM) is a motor that uses

permanent magnets to produce the air gap magnetic field rather than using
electromagnets. These motors have significant advantages, attracting the interest of
researchers and industry for use in many applications.

Figure 4.2 Permanent Magnet Synchronous Motor

The use of permanent magnets (PMs) in construction of electrical machines
brings the benefits of no electrical energy is absorbed by the field excitation system
and thus there are no excitation losses which means substantial increase in the
efficiency, higher torque and/or output power per volume than when using
electromagnetic excitation, better dynamic performance than motors with
electromagnetic excitation (higher magnetic flux density in the air gap),simplification
of construction and maintenance, reduction of prices for some types of machines.

4.2.1 Permanent Magnet Materials

The properties of the permanent magnet material will affect directly the
performance of the motor and proper knowledge is required for the selection of the

21
materials and for understanding PM motors.
The earliest manufactured magnet materials were hardened steel. Magnets
made from steel were easily magnetized. However, they could hold very low energy
and it was easy to demagnetize. In recent years other magnet materials such as
Aluminum Nickel and Cobalt alloys (ALNICO), Strontium Ferrite or Barium Ferrite
(Ferrite), Samarium Cobalt (First generation rare earth magnet) (SmCo) and
Neodymium Iron-Boron (Second generation rare earth magnet) (NdFeB) have been
developed and used for making permanent magnets.
The rare earth magnets are categorized into two classes: Samarium Cobalt
(SmCo) magnets and Neodymium Iron Boride (NdFeB) magnets. SmCo magnets
have higher flux density levels but they are very expensive. NdFeB magnets are the
most common rare earth magnets used in motors these days. A flux density versus
magnetizing field for these magnets is illustrated in Figure 4.3.

Figure 4.3 Flux Density Vs Magnetizing Field of Permanent Magnetic Materials

4.2.2 CLASSIFICATION OF PERMANENT MAGNET MOTORS

The permanent magnet synchronous motors are broadly classified according to

the direction of field flux density distribution and permanent magnet radial field
motors as follows.
4.2.2.1 Direction of field flux
PM motors are broadly classified by the direction of the field flux. The
first field flux classification is radial field motor meaning that the flux is along the

22
radius of the motor. The second is axial field motor meaning that the flux is
perpendicular to the radius of the motor. Radial field flux is most commonly used in
motors and axial field flux have become a topic of interest for study and used in a few
applications.

4.2.2.2 Flux density distribution

PM motors are classified on the basis of the flux density distribution and the
shape of current excitation. They are PMSM and PM brushless motors (BLDC). The
PMSM has a sinusoidal-shaped back EMF and is designed to develop sinusoidal back
EMF waveforms.
They have the following:
1. Sinusoidal distribution of magnet flux in the air gap
2. Sinusoidal current waveforms
3. Sinusoidal distribution of stator conductors.
BLDC has a trapezoidal-shaped back EMF and is designed to develop trapezoidal
back EMF waveforms. They have the following:
1. Rectangular distribution of magnet flux in the air gap
2. Rectangular current waveform
3. Concentrated stator windings.

4.2.2.3 Permanent magnet radial field motors

In PM motors, the magnets can be placed in two different ways on the rotor.
Depending on the placement they are called either as surface permanent magnet motor
or interior permanent magnet motor.

23
 Surface mounted PM motors have a surface mounted permanent magnet rotor.
Each of the PM is mounted on the surface of the rotor, making it easy to build, and
specially skewed poles are easily magnetized on this surface mounted type to
minimize cogging torque. This configuration is used for low speed applications
because of the limitation that the magnets will fly apart during high-speed operations.
These motors are considered to have small saliency, thus having practically equal
inductances in both axes .The permeability of the permanent magnet is almost that of
the air, thus the magnetic material becoming an extension of the air gap.
For a surface permanent magnet motor Ld = Lq .
The rotor has an iron core that may be solid or may be made of punched
laminations for simplicity in manufacturing. Thin permanent magnets are mounted on
the surface of this core using adhesives. Alternating magnets of the opposite
magnetization direction produce radially directed flux density across the air gap. This
flux density then reacts with currents in windings placed in slots on the inner surface
of the stator to produce torque.

Figure 4.4 shows the placement of the magnet.

Figure 4.4 Surface Permanent Magnet Motor

Interior PM Motors have interior mounted permanent magnet rotor as

shown in Figure 4.4. Each permanent magnet is mounted inside the rotor. It is not as
common as the surface mounted type but it is a good candidate for high-speed
operation. There is inductance variation for this type of rotor because the permanent
magnet part is equivalent to air in the magnetic circuit calculation. These motors are
considered to have saliency with q axis inductance greater than the d axis inductance
(Lq > Ld).
24
 Figure 4.5 Interior Permanent Magnet Motor

4.3 Position Sensor

Operation of permanent magnet synchronous motors requires position sensors

in the rotor shaft when operated without damper winding. The need of knowing the
rotor position requires the development of devices for position measurement. There
are four main devices for the measurement of position, the potentiometer, linear
variable differential transformer, optical encoder and resolvers. The ones most
commonly used for motors are encoders and revolvers. Depending on the application
and performance desired by the motor a position sensor with the required accuracy
can be selected.

4.3.1 Optical Encoders

The most popular type of encoder is the optical encoder as shown in Figure
4.6, which consists of a rotating disk, a light source, and a photo detector (light
sensor). The disk, is mounted on the rotating shaft, has coded patterns of opaque and
transparent sectors. As the disk rotates, these patterns interrupt the light emitted onto
the photo detector, generating a digital pulse or output signal.

Figure 4.6 Optical Encoder

25
 Optical encoders offer the advantages of digital interface. There are two types
of optical encoder’s Incremental encoder and absolute encoder.

4.3.1.1 Incremental encoders

Incremental encoders have good precision and are simple to implement but
they suffer from lack of information when the motor is at rest position and in order for
precise position the motor most be stop at the starting point.
The most common type of incremental encoder uses two output channels (A
and B) to sense position. Using two code tracks with sectors positioned 90° degrees
out of phase, the two output channels of the quadrature encoder indicate both position
and direction of rotation as shown in Figure 4.7. If A leads B, for example, the disk is
rotating in a clockwise direction. If B leads A, then the disk is rotating in a counter-
clockwise direction. By monitoring both, the number of pulses and the relative phase
of signals A and B, it's possible to track position and direction of rotation. Some
quadrature encoders also include a third output channel, called a zero or index or
reference signal, which supplies a single pulse per revolution. This single pulse is
used for precise determination of a reference position. The precision of the encoder is
fix by its code disk but it can be increased by detecting the Up and Down transitions
on both the A and B channels.

Figure 4.7 Quadrature Encoder Channels

4.3.1.2 Absolute encoders
The absolute encoder, as shown in Figure 4.8 captures the exact position of the
rotor with a precision directly related to the number of bits of the encoder. It can
rotate indefinitely and even if the motor stops the position can be measured or
obtained.
It provides a “whole word” output with a unique code pattern representing
each position. This code is derived from independent tracks on the encoder disc (one
for each “bit” of resolution) corresponding to individual photo detectors. The output
from these detectors is HI (light) or LO (dark) depending on the code disc pattern for

26
that particular position.

Figure 4.8 Absolute Encoder

Absolute encoders are used in applications where a device is inactive for long
periods of time or moves at a slow rate, such as flood gate control, telescopes, cranes,
valves, etc. They are also recommended in systems that must retain position
information through a power outage.

4.4 POSITION REVOLVER

Position revolver as shown in Figure 4.9, also called rotary transformers works
on the transformer principle. The primary winding is placed on the rotor and
depending upon the rotor shaft angle the induced voltage at the two secondary
windings of the transformer shifted by 90° would be different. The position can be
calculated using the two voltages.
The resolver is basically a rotary transformer with one rotating reference
winding (Vref) and two stator windings. The reference winding is fixed on the rotor,
and therefore, it rotates jointly with the shaft passing the output windings, as is
depicted in Figure 4.9. Two stator windings are placed in quadrature (shifted by 90°)
with one another and generate the sine and cosine voltages (Vsin , Vcos) respectively.
Both windings will be further referred to as output windings. In consequence of the
excitement applied on the reference winding Vref and along with the angular
movement of the motor shaft θ, the respective voltages are generated by resolver
output windings Vsin , Vcos.

Figure 4.9 Resolver

27
 The frequency of the generated voltages is identical to the reference voltage and
their amplitudes vary according to the sine and cosine of the shaft angle θ.
Considering that one of the output windings is aligned with the reference winding,
then it is generated full voltage on that output winding and zero voltage on the other
output winding and vice versa. The rotor angle θ can be extracted from these voltages
shown in Figure 4.10.

Figure 4.10 Excitation and Output Signal of the Resolver

The shaft angle can be determined by an Inverse Tangent [5] function of the
quotient of the sampled resolver output voltages Vsin , Vcos . This determination can be
expressed, in terms of resolver output voltages, as follows:

U sin
θ= α tan ( )
U cos

28
 CHAPTER V
MODELING OF PM DRIVE SYSTEM

29
 5 MODELING OF PM DRIVE SYSTEM

5.1 INTRODUCTION

This chapter deals with the detailed modeling of a permanent magnet

synchronous motor. Field oriented control of the motor in constant torque and flux-
weakening regions are discussed. Closed loop control of the motor is developed using
a PI controller in the speed loop. Design of the speed controller is discussed.

5.2 DETAILED MODELING OF PMSM

Detailed modeling of PM motor drive system is required for proper simulation

of the system. The d-q model has been developed on rotor reference frame as shown
in Figure 5.1. At any time t, the rotating rotor d-axis makes and angle θ r with the fixed
stator phase axis and rotating stator MMF makes an angle α with the rotor d-axis.
Stator MMF rotates at the same speed as that of the rotor.

Figure 5.1 Motor Axis

The model of PMSM without damper winding has been developed on rotor
reference frame using the following assumptions:
1) Saturation is neglected.
2) The induced EMF is sinusoidal.
3) Eddy currents and hysteresis losses are negligible.
4) There are no field current dynamics.

30
 5) Variation in rotor temperature with time is considered to be negligible
The stator equations in the rotor reference frames using flux linkages are taken
to derive the model of the PMSM.
The stator flux-linkage equations are given by:
Vqs = Rq iqs+ pλqs+ ωrλds 5.1
Vds = Rd ids − ωr λqs+ pλds 5.2
The d and q axes Flux Linkages in rotor reference frame are given by
λqs = Lsiqs+ Lmiqr 5.3
λds = Lsids + Lmidr 5.4
The PM excitation can be modeled as a constant current source i fr. The rotor flux is
along the d-axis, so the d-axis rotor current is ifr. The q-axis current in the motor is
zero, because there is no flux along this axis in the rotor, by assumption
λqs = Lqiqs 5.5
λds = Ldsids+Lmifr 5.6
Where Lm= mutual inductance between the stator winding and rotor magnets
By substituting equations 5.5 and 5.6 into 5.1 and 5.2
Vqs = Rq iqs + ωr ( Ld ids + Lmifr ) + p Lq iqs 5.7
Vds = Rs ids − ωr Lq iqs + p ( Ld ids + Lmifr ) 5.8
Let Lmifr = λf
Arranging equations 5.7 and 5.8 in matrix form

Vq R s + ρL q ωr Ld iq ωr λ f
( )(
Vd
=
−ωr Lq Rs + ρL d )( ) ( )
id
+
ρλf

5.9
The developed torque motor is being given by

T e= 32 ( P2 ) ( λ d i q − λ q i d )
5.10
The mechanical Torque equation is
dω m
T e=T L + Bωm +J
dt
5.11
Solving for the rotor mechanical speed form equation 3.9
T e −T L −Bωm
ω m=∫ ( J ) dt
31
 5.12

And
ω m =ωr ( 2P )
5.13
In the above equations ωr is the rotor electrical speed where as ω m is the rotor
mechanical speed.
5.2.1 Parks Transformation and Dynamic d q Modeling
The dynamic d q modeling is used for the study of motor during transient and
steady state. It is done by converting the three phase voltages an currents to dqo
variables by using Parks transformation [5].
Converting the phase voltages variables Vabc to Vdqo variables in rotor reference
frame the following equations are obtained
cos θr cos( θ r −120) cos ( θr +120)
Vq

[]
Vo [
2 sin θ r
V d =3
1
2
sin ( θr −120 )
1
2
sin (θ r + 120)
1
2
][ ]
Va
Vb
Vc

5.14
Convert Vdqo to Vabc

Va cos θr sin θ r 1 Vq

[ ][
V b = cos( θ r −120)
Vc cos( θ r +120 )
sin ( θr −120 )
sin ( θr + 120) ][ ]
1 Vd
1 Vo
5.15

5.2.2 Equivalent Circuit of Permanent Magnet Synchronous Motor

Equivalent circuits of the motors are used for study and simulation of motors.
From the d-q modeling of the motor using the stator voltage equations the equivalent
circuit of the motor can be derived. Assuming rotor d axis flux from the permanent
magnets is represented by a constant current source as described in the following
equation λ f = Ldmif , Figure 5.2 is obtained.

32
 Figure 5.2 Permanent Magnet Motor Electric Circuits without Damper Windings

5.3 PM MOTOR CONTROL

Control of PM motors is performed using field oriented control for the

operation of synchronous motor as a DC motor. The stator windings of the motor are
fed by an inverter that generates a variable frequency variable voltage. Instead of
controlling the inverter frequency independently, the frequency and phase of the
output wave are controlled using a position sensor as shown in Figure 5.3.

Figure 5.3 Self Control Synchronous Motor

Field oriented control was invented in the beginning of 1970s and it
demonstrates that an induction motor or synchronous motor could be controlled like a
separately excited DC motor by the orientation of the stator MMF or current vector in
relation to the rotor flux to achieve a desired objective. In order for the motor to
behave like DC motor, the control needs knowledge of the position of the
instantaneous rotor flux or rotor position of permanent magnet motor. This needs a
resolver or an absolute optical encoder. Knowing the position, the three phase currents
can be calculated. Its calculation using the current matrix depends on the control
desired. Some control options are constant torque and flux weakening. These options
are based in the physical limitation of the motor and the inverter. The limit is
established by the rated speed of the motor, at which speed the constant torque
operation finishes and the flux weakening starts as shown in Figure 5.4.

Figure 5.4 Steady State Torque versus Speed

33
5.3.1 Field Orientated Control Of PM Motors
The Field Orientated Control (FOC) consists of controlling the stator currents
represented by a vector. This control is based on projections which transform a three
phase time and speed dependent system into a two co-ordinate (d and q co-ordinates)
time invariant system. These projections lead to a structure similar to that of a DC
machine control. Field orientated controlled machines need two constants as input
references: the torque component (aligned with the q co-ordinate) and the flux
components (aligned with d co-ordinate). As Field Orientated Control is simply based
on projections the control structure handles instantaneous electrical quantities. This
makes the control accurate in every working operation (steady state and transient) and
independent of the limited bandwidth mathematical model. The FOC thus solves the
classic scheme problems, in the following ways:
 the ease of reaching constant reference (torque component and flux component
of the stator current)
 the ease of applying direct torque control because in the (d,q) reference frame
the expression of the torque is:
m iR Sq µy
By maintaining the amplitude of the rotor flux (y R) at a fixed value we have a
linear relationship between torque and torque component (iSq). We can then control
the torque by controlling the torque component of stator current vector.
The PMSM control is equivalent to that of the DC motor by a decoupling
control known as field oriented control or vector control. The vector control separates
the torque component of current and flux channels in the motor through its stator
excitation.
The vector control of the PM synchronous motor is derived from its dynamic
model. Considering the currents as inputs, the three currents are:
i a =I m sin( ω r t +α ) 5.16
2π
i b =I m sin ( ω r t +α − 3
)
5.17
2π
i c=I m sin ( ωr t + α + 3
)
5.18

34
Writing equations 5.16 to 5.18 in the matrix form:

ia cos( ωr t +α )

()( 2π

)
i b = cos( ωr t +α − 3 ) ( I m )
ic cos( ωr t +α + 23π )

5.19
Where α is the angle between the rotor field and stator current phasor, ω r is the
electrical rotor speed.
The previous currents obtained are the stator currents that must be transformed
to the rotor reference frame with the rotor speed ωr , using Park’s transformation. The
q and d axis currents are constants in the rotor reference frames since α is a constant
for a given load torque. As these constants, they are similar to the armature and field
currents in the separately excited DC machine. The q axis current is distinctly
equivalent to the armature current of the DC machine; the d axis current is field
current, but not in its entirety. It is only a partial field current; the other part is
contributed by the equivalent current source representing the permanent magnet field.
For this reason the q axis current is called the torque producing component of the
stator current and the d axis current is called the flux producing component of the
stator current.
Substituting equation 5.19 and 5.14 is obtained id and iq in terms of Im as follows:

iq sin α
( ) (
id
=Im
cos α ) 5.20
Using equations 5.1, 5.2, 5.10 and 5.20 the electromagnetic torque equation is
obtained as given below.
2
T e= 3
2
P
2 [ 1
2
( Ld − Lq ) I m
sin 2 α + λ f I m sin α ] 5.21

5.3.1.1 Constant torque operation

Constant torque control strategy is derived from field oriented control, where
the maximum possible torque is desired at all times like the DC motor. This is
performed by making the torque producing current iq equal to the supply current Im.
That results in selecting the angle α to be 90 º degrees according to equation 5.20. By
making the id current equal to zero the torque equation can be rewritten as:

35
 3
T e=( 2 )( P2 ) λ f i q 5.22
Assuming that:
3 P
k t = ( 2 )( 2 ) λf
5.23
The torque is give by

T e=k t −i q 5.24
Like the DC motor, the torque is dependent of the motor current.

5.3.1.2 Flux-weakening

Flux weakening is the process of reducing the flux in the d axis direction of
the motor which results in an increased speed range. The motor drive is operated with
rated flux linkages up to a speed where the ratio between the induced EMF and stator
frequency (V/f) is maintained constant. After the base frequency, the V/f ratio is
reduced due to the limit of the inverter DC voltage source which is fixed. The
weakening of the field flux is required for operation above the base frequency . This
reduces the V/f ratio. This operation results in a reduction of the torque proportional to a
change in the frequency and the motor operates in the constant power region.
The rotor flux of PMSM is generated by permanent magnet which cannot be
directly reduced as induction motor. The principle of flux-weakening control of
PMSM is to increase negative direct axis current and use armature reaction to reduce
air gap flux, which equivalently reduces flux and achieves the purpose of flux-
weakening control.
This method changes torque by altering the angle between the stator MMF and
the rotor d axis. In the flux weakening region where ω r > ωrated angle α is controlled by
proper control of id and iq for the same value of stator current. Since i q is reduced the
output torque is also reduced. The angle α can be obtained as:

iq
α =Tan−1
( )
id
5.25
The current Im is related to id and iq by:

36
 2 2
I m= i √ d
+i q
5.26
Flux-weakening control realization
The realization process of equivalent flux-weakening control is as follows,
1) Measuring rotor position and speed ωr from a sensor which is set in motor rotation
axis.
2) The motor at the flux weakening region with a speed loop, Te* is obtained from the
PI controller.
3) Calculate Iq* using equation 5.22
¿
e
T

( id=
3
( )( )
2
P
λ
2 f
)
5.27
5) Calculate Id* using equation:

¿
d
λ d −λ f
i =
Ld
5.28
6) Calculate α using equation 5.25.
7) Using α and rotor position the controller will generate the reference currents as per
equation 5.19.
8) Then the current controller makes uses of the reference signals to control the
inverter for the desired output currents.
9) The load torque is adjust to the maximum available torque for the reference speed

ω rated
T L=T e ( rated )
ωr
5.29
5.4 SPEED CONTROL OF PM MOTOR

Many applications, such as robotics and factory automation, require precise

control of speed and position. Speed Control Systems allow one to easily set and
adjust the speed of a motor. The control system consists of a speed feedback system, a

37
motor, an inverter, a controller and a speed setting device. A properly designed
feedback controller makes the system insensible to disturbance and changes of the
parameters.
The purpose of a motor speed controller is to take a signal representing the
demanded speed, and to drive a motor at that speed. Closed Loop speed control
systems have fast response, but become expensive due to the need of feedback
components such as speed sensors.

5.4.1 Implementation of the Speed Control Loop

For a PM motor drive system with a full speed range the system will consist of
a motor, an inverter, a controller (constant torque and flux weakening operation,
generation of reference currents and PI controller) as shown in Figure 5.5

Figure 5.5 Block Diagram

The operation of the controller must be according to the speed range. For
operation up to rated speed it will operate in constant torque region and for speeds
above rated speed it will operate in flux-weakening region. In this region the d-axis
flux and the developed torque are reduced.
Speed controller calculates the difference between the reference speed and the
actual speed producing an error, which is fed to the PI controller. PI controllers are
used widely for motion control systems. They consist of a proportional gain that
produces an output proportional to the input error and an integration to make the
steady state error zero for a step change in the input. Block diagram of the PI
controller is shown in Figure 5.6.

38
 Figure 5.6 PI Controller
Speed control of motors mainly consist of two loops the inner loop for current
and the outer loop for speed. The order of the loops is due to their response, how fast
they can be changed. This requires a current loop at least 10 times faster than the
speed loop. Since the PMSM is operated using field oriented control, it can be
modeled like a DC motor. The design begins with the innermost current loop by
drawing the block diagram. But in PMSM drive system the motor has current
controllers which make the current loop. The current control is performed by the
comparison of the reference currents with the actual motor currents.
The design of the speed loop assumes that the current loop is at least 10 times
faster than speed loop, allowing reducing the system block diagram by considering
the current loop to be of unity gain as shown in Figure 5.7.

Figure 5.7 Block Diagram of Speed Loop

The open loop transfer function of the motor is given by:

5.30

kT = PM flux = λf
The crossover frequency has been selected an order smaller than the current
loop. To satisfy dynamic response without oscillations the phase margin ( φPM )
should be greater than 45º, preferably close to 60º. Knowing the motor parameters and
39
phase margin, the ki and kp gains can be obtained for the motor controller using
equations 5.31 and 5.32.
Phase Margin = φOL + 180°
ki kT α kp
|
Js 2 ( 1+ s
ki )
|
s = jω
=1

5.31
ki kT α kp
∠
( Js 2 ( 1+s
ki ))
=1800 + φ PM

5.32
The gains for the speed controller was obtained using the motor parameters and by
selecting a crossover frequency.
The selected values are
fc (crossover frequency)= 100Hz
J = 0.000179
kT = λf = 0.272
Using equation 3.29 and 3.30 and motor parameters the values of ki and kp are
obtained as 129.9014 and 0.3581 respectively.

CHAPTER VI
SPEED CHECK IMPLEMENTATION OF VECTOR
CONTROLLED PERMANENT MAGNET

40
SYNCHORONOUS MOTOR DRIVE IN
MATLAB/SIMULINK

41
 6 SPEED CHECK IMPLEMENTATION OF
VECTORVCONTROLLED PMSM DRIVE IN
MATLAB/SIMULINK

6.1 INTRODUCTION

This chapter describes the advantages of SIMULINK and different tools

available for electrical and electronic systems simulation and then justification is
given for selecting SIMULINK for the PMSM system. Block by block an explanation
is given for SIMULINK simulation of the drive system.

6.2 ADVANTAGES OF SIMULINK

There are two major advantages to perform simulation rather than actually
building the design and testing it. The biggest of these advantages is that it is
economical. Designing, building, testing, redesigning, rebuilding, retesting for
anything can be an expensive project. Simulations take the building/rebuilding phase
out of the loop by using the model already created in the design phase. Most of the
time, the simulation testing is cheaper and faster than performing the multiple tests of
the design each time. Considering the typical university budget cheaper is usually a
very good thing. In the case of an electric thruster the test must be run inside of a
vacuum tank. Vacuum tanks are very expensive to buy, run, and maintain. One of the
main tests of an electric thruster is the lifetime test, which means that the thruster is
running pretty much constantly inside of the vacuum tank for 10,000+ hours. This is
pouring money down a drain compared to the price of the simulation.

The second major advantage of a simulation is the level of detail that you can
get from a simulation. A simulation can give you results that are not experimentally
measurable with our current level of technology. Results such as surface interactions
on an atomic level, flow at the exit of a micro electric thruster, or molecular flow
inside of a star are not measurable by any current devices. A simulation can give these
results when problems such as it's too small to measure, the probe is too big and is
skewing the results, and any instrument would turn to a gas at those temperatures

42
come into the conversation. You can set the simulation to run for as many time steps
you desire and at any level of detail you desire the only restrictions are your
imagination, your programming skills, and your CPU.

6.3 SIMULATION TOOLS  

Study of electric motor drives needs the proper selection of a simulation tool.
Their complex models need computing tools capable of performing dynamic
simulations. Today with the growth in computational power there is a wide selection
of software titles available for electrical simulations such as ACSL, ESL, EASY5, and
PSCSP are for general systems and SPICE2, EMTP, and ATOSEC5 for simulating
electrical and electronic circuits. IESE and SABER are examples of general-purpose
electrical network simulation programs that have provisions for handling user-defined
modules. SIMULINK® is a toolbox extension of the MATLAB program. It is a
program for simulating dynamic systems.
SIMULINK has the advantages of being capable of complex dynamic system
simulations, graphical environment with visual real time programming and broad
selection of tool boxes. The simulation environment of SIMULINK has a high
flexibility and expandability which allows the possibility of development of a set of
functions for a detailed analysis of the electrical drive .Its graphical interface allows
selection of functional blocks, their placement on a worksheet, selection of their
functional parameters interactively, and description of signal flow by connecting their
data lines using a mouse device. System blocks are constructed of lower level blocks
grouped into a single maskable block. SIMULINK simulates analogue systems and
discrete digital systems.

6.4 SIMULINK SIMULATION OF PMSM DRIVE  

The PM motor drive simulation was built in several steps like abc phase
transformation to dqo variables, calculation torque and speed, and control circuit. The
abc phase transformation to dqo variables is built using Parks transformation and for
the dqo to abc the reverse transformation is used. For simulation purpose the voltages
are the inputs and the current are output. Parks transformation used for converting V abc
to Vdqo is shown in Figure 6.1 and the reverse transformation for converting I dqo to Iabc
is shown in Figure 6.2.
43
6.4.1 Vabc to Vdqo block and Idqo to Iabc block
The below Figures 6.1 and 6.2 is implemented in the MATLAB/SIMULINK.
These are used to convert V abc to Vdqo and Idqo to Iabc i.e. it used to convert rotating
frame to stationary reference frame parameters or vice versa.

Figure 6.1 Vabc to Vdqo block

cos θr cos( θ r −120) cos ( θr +120)

Vq

[]
Vo
1
[
2 sin θ r
V d =3

2
sin ( θr −120 )
1
2
sin (θ r + 120)
1
2
][ ]
Va
Vb
Vc
6.1
With the use of equation 6.1 the abc parameters are converted into dqo parameters.

Figure 6.2 Idqo to Iabc block

44
 ia cos( ωr t +α )

()(
ic
2π
i b = cos( ωr t +α − 3 ) ( I m )
2π
cos( ωr t +α + 3 ) ) 6.2
With the use of above equation the dqo parameters are converted into abc parameters.

6.4.2 d-axis circuit and q-axis circuit

The d and q axis motor circuits built using SIMULINK elements are shown in
Figure 6.3
This block is developed using equation6.3 for d-axis circuit and q-axis circuit

iq sin α
( ) (
id
=Im
cos α ) 6.3

Figure 6.3 d-axis circuit and q-axis circuit

q-axis stator current, iq 2.015 A

d-axis stator current, id 2.745e.008 A

Rotor electrical speed, ω 52.36 radians

Stator resistance, Rs 6.852 Ω

45
 d- axis stator voltage, Vds -1.277 V

q-axis stator voltage, Vqs 11.18 V

Lndd 0.02439 H

Lndq 0.0827 H

Parameters Of The d and q Circuits

6.4.3 Load Torque Block
Figure 4.4 shows the load torque block in SIMULINK. This block is
developed using equation 4.5 for load torque developed

Figure 6.4 Load Torque Block

The developed torque motor is being given by

T e= 32 ( P2 ) ( λ d i q − λ q i d )
6.5
6.4.4 Speed Block
The speed of the motor is obtained using Figure 6.4 and equation 6.6. The
developed speed block is shown in Figure 6.5.

Figure 6.5 Speed Block

46
The mechanical Torque equation is
dω m
T e=T L + Bωm +J
dt 6.6
Solving for the rotor mechanical speed form equation 6.6
T e −T L −Bωm
ω m=∫ ( J ) dt 6.7
And
ω m =ωr ( 2P ) 6.8
In the above equations ωr is the rotor electrical speed where as ωm is the rotor
mechanical speed.
6.4.5 Vector Control Block
The vector control requires a block for the calculation of the reference current
using α angle, the position of the rotor and the magnitude of the I m. The block is
shown in Figure 6.6. It is built using equation 3.17.

ia cos( ωr t +α )

()( 2π

)
i b = cos( ωr t +α − 3 ) ( I m )
ic cos( ωr t +α + 23π )

6.9
Where α is the angle between the rotor field and stator current phasor, ω r is the
electrical rotor speed.

47
 Figure 6.6 Vector Control Block
The stator currents that must be transformed to the rotor reference frame with
the rotor speed ωr, using Park’s transformation. The q and d axis currents are constants
in the rotor reference frames since α is a constant for a given load torque. As these
constants, they are similar to the armature and field currents in the separately excited
DC machine. The q axis current is distinctly equivalent to the armature current of the
DC machine; the d axis current is field current, but not in its entirety. It is only a
partial field current; the other part is contributed by the equivalent current source
representing the permanent magnet field. For this reason the q axis current is called
the torque producing component of the stator current and the d axis current is called
the flux producing component of the stator current.
The PI controller increases the order and type of the system and also improves
the steady state response of the system. There are 3 PI controllers in the vector
controller. The proportionality constant values are taken as Kp1=2.285e.007; Ki1=2.015;
Kp2=-8.235e.008; Ki2=-1.277; Kp3=-1.746e.006; Ki3=11.18.
When the input signal is within the range specified by the Lower limit and
Upper limit parameters, the input signal passes through unchanged. When the input
signal is outside these bounds, the signal is clipped to the upper or lower bound. The
upper limit value is taken as Iphmax (2.5A) and the lower limit value is taken as - I phmax (-
2.5A). When the Lower limit and Upper limit parameters are set to the same value,
the block outputs that value.
6.4.6 PM Motor Drive System in SIMULINK
Equivalent circuits of the motors are used for study and simulation of motors.
From the d-q modeling of the motor using the stator voltage equations the equivalent
circuit of the motor can be derived. Assuming rotor d axis flux from the permanent
magnets is represented by a constant current source as described in the following
equation

48
 Figure 6.8 PM Motor Drive System in SIMULINK

6.4.7 Vector Control Of PMSM Drive

Using all the drive system blocks the complete system block has been developed as
shorn in Figure6.9

FIGURE 6.9 Vector Control of PMSM Drive

49
 The operation of the controller must be according to the speed range. For
operation up to rated speed it will operate in constant torque region and for speeds
above rated speed it will operate in flux-weakening region. In this region the d-axis
flux and the developed torque are reduced.
Speed controller calculates the difference between the reference speed and the
actual speed producing an error, which is fed to the PI controller. PI controllers are
used widely for motion control systems. They consist of a proportional gain that
produces an output proportional to the input error and an integration to make the
steady state error zero for a step change in the input.
Speed control of motors mainly consist of two loops the inner loop for current
and the outer loop for speed. The order of the loops is due to their response, how fast
they can be changed. This requires a current loop at least 10 times faster than the
speed loop. Since the PMSM is operated using field oriented control, it can be
modeled like a DC motor. The design begins with the innermost current loop by
drawing the block diagram. But in PMSM drive system the motor has current
controllers which make the current loop. The current control is performed by the
comparison of the reference currents with the actual motor currents.
The design of the speed loop assumes that the current loop is at least 10 times
faster than speed loop, allowing reducing the system block diagram by considering
the current loop to be of unity gain
The operation of the controller must be according to the speed range. For
operation up to rated speed it will operate in constant torque region and for speeds
above rated speed it will operate in flux-weakening region. In this region the d-axis
flux and the developed torque are reduced.
Speed controller calculates the difference between the reference speed and the
actual speed producing an error, which is fed to the PI controller. PI controllers are
used widely for motion control systems. They consist of a proportional gain that
produces an output proportional to the input error and an integration to make the
steady state error zero for a step change in the input.

50
 CHAPTER VII

SIMULATION RESULTS

51
 7. SIMULATION RESULTS

7.1 CASE 1

In this the speed and the torque are kept constant. Hence below are the
waveforms which are obtained from simulation in MATLAB.
7.1.1 Reference Torque and Actual Torque

Figure 7.1 Reference torque and actual torque

The above Figure 7.1 waveforms are the simulation result of torque variation
which is implemented in MATLAB/SIMULINK. It give information about the
variation of torque with respective to the reference torque. In this the actual torque
follows the reference torque. Hence under steady state there is no change in actual
torque i.e. torque is not varied at any point of time.

52
7.1.2 Actual Speed and Reference Speed

Figure 7.2 (a) Actual speed and (b) Reference speed

The above Figure 7.2 waveforms are the simulation result of speed variation
which is implemented in MATLAB/SIMULINK. It give information about the
variation of speed with respective to the reference speed. In this the actual speed
follows the reference speed. Hence under steady state there is no change in actual
speed i.e. speed or load is not varied at any point of time.

7.1.3 Iabc currents and Vabc voltages

The below Figure 7.3 waveforms are the simulation result of current and
voltage variation which is implemented in MATLAB/SIMULINK. These waveforms
are the three phase current and voltages in rotating frame of reference under steady
state.

53
 Figure 7.3 Iabc currents and Vabc voltage

7.1.4 Vafa and Vbeta

The below Figure 7.4 waveforms are the simulation result of V afa and Vbeta
variation which is implemented in MATLAB/SIMULINK. Hear three phase voltages
are converted into two phase waveforms i.e. Vabc is converted into Vafa and Vbeta
(which is known as clark’s transformation or stationary reference frame)

Figure 7.4 (a) Vafa and (b) Vbeta

54
7.15 Error Signal between Reference Speed and Actual Speed

Figure 7.5 Error Signal between Reference Speed and Actual Speed

The above Figure 7.5 waveform is the simulation result of error variation
which is implemented in MATLAB/SIMULINK. This is the waveform of the error
signal between reference speed and actual speed of the permanent magnet
synchronous motor drive. In this steady state the error is nearly maintained to zero as
shown in Figure 7.6 above.

55
7.1.6 Vd, Vq Voltage and Id-ref, Iq-ref

The below Figure 7.6 and 7.7 waveforms are the simulation result of V d
voltage, Vq voltage variation and Id-ref, Iq-ref which is implemented in
MATLAB/SIMULINK. Hear two phase Vafa and Vbeta voltages are converted into two
phase Vd and Vq waveforms i.e. Vafa and Vbeta is converted into Vd and Vq as shown
in Figure 5.7 (which is known as park’s transformation or rotating frame of reference)
Under steady state at speed 100rpm. In this Id-ref, Iq-ref are the reference signals and Id-ref
is maintained at zero. As shown in Figure 7.7

Figure 7.6 (a)Vd voltage and (b) Vq voltage

56
 Figure 7.7 (a) Id-ref and (b) Iq-ref

7.2 CASE 2

In this the speed is varied and the torque is kept constant. Hence below are the
waveforms which are obtained from simulation in MATLAB.
7.2.1 Actual Speed and Reference Speed

Figure 7.8 (a) Actual speed and (b) Reference speed

57
 The above Figure 7.8 waveforms are the simulation result of speed variation
which is implemented in MATLAB/SIMULINK. It give information about the
variation of speed with respective to the reference speed. In this the actual speed
follows the reference speed. Hence under speed dynamic state there is change in
actual speed i.e. speed or load is varied at any point of time. In this Figure 7.9 the
speed is changed to 500rpm at some point of time.

7.2.2 Reference Torque and Actual Torque

The below Figure 7.9 waveforms are the simulation result of torque variation
which is implemented in MATLAB/SIMULINK. It give information about the
variation of torque with respective to the reference torque. In this the actual torque
follows the reference torque. Hence under speed dynamics there is no change in actual
torque i.e. torque is not varied at any point of time.

Figure 7.9 (a) Reference torque and (b) Actual torque

7.2.3 Iq, Id, Id-ref and Iq-ref

The below Figure 7.10 and 7.11 waveforms are the simulation result of I q, Id,
Id-ref and Iq-ref currents which is implemented in MATLAB/SIMULINK. In this I d-ref, Iq-
ref is the reference signals are maintained at zero. Thus the actual signals which are

58
generated in the control analysis are as shown in the Figure 7.11 are following the
reference signals. This is under speed dynamics.

Figure 7.10 (a) Iq current and (b) Id current

Figure 7.11 (a) Id-ref and (b) Iq-ref

7.2.4 Iabc Currents and Vabc Voltages

The below Figure 7.12 waveforms are the simulation result of current and
voltage variation which is implemented in MATLAB/SIMULINK. These waveforms
are the three phase current and voltages in rotating frame of reference under speed

59
dynamics state. The variation in current and voltages in the Figures is due to the
change in speed.

Figure 7.12 (a) Iabc currents and (b) Vabc voltages

7.2.5 Error Signal between Reference Speed and Actual Speed

The below Figure 7.13 waveform is the simulation result of error variation,
which is implemented in MATLAB/SIMULINK. This is the waveform of the error
signal between reference speed and actual speed of the permanent magnet
synchronous motor drive. In this speed dynamics state the error is nearly maintained
to zero as shown in Figure 7.13 below.

Figure 7.13 Error Signal between Reference Speed and Actual Speed

60
7.2.6 Vd Voltage and Vq Voltage

The below Figure 7.14 waveforms are the simulation result of Vd voltage and
Vq voltage variation which is implemented in MATLAB/SIMULINK. Here two
phases Vafa and Vbeta voltages are converted into two phase V d and Vq waveforms i.e.
Vafa and Vbeta is converted into Vd and Vq (which is known as park’s transformation
or rotating frame of reference) Under speed dynamics at speed 1000rpm to 500rpm.

Figure 7.14 (a)Vd voltage and (b) Vq voltage

7.3 CASE 3

In this the torque is varied and the speed is kept constant. Hence below are the
waveforms which are obtained from simulation in MATLAB.
7.3.1 Actual Speed and Reference Speed

The below Figure 7.15 waveforms are the simulation result of speed variation
which is implemented in MATLAB/SIMULINK. It give information about the
variation of speed with respective to the reference speed. In this the actual speed
follows the reference speed. Hence under torque dynamics there is no change in actual
speed i.e. speed or load is not varied at any point of time.

61
 Figure 7.15 (a) Actual speed and (b) Reference speed

7.3.2 Reference Torque and Actual Torque

The below Figure 7.16 waveforms are the simulation result of torque variation
which is implemented in MATLAB/SIMULINK. It give information about the
variation of torque with respective to the reference torque. In this the actual torque
follows the reference torque. Hence under torque dynamics state there is change in
actual torque i.e. torque is varied at any point of time.

Figure 7.16 (a) Reference torque and (b)actual torque

7.3.3 Iabc Current and Vabc Voltage

The below Figure 7.17 waveforms are the simulation result of current and
voltage variation which is implemented in MATLAB/SIMULINK. These waveforms

62
are the three phase current and voltages in rotating frame of reference under torque
dynamic state.

Figure 7.17 (a) Iabc currents and (b)Vabc voltages

7.3.4 Error Signal between Reference Speed and Actual Speed

The below Figure 7.18 waveform is the simulation result of error variation
which is implemented in MATLAB/SIMULINK. This is the waveform of the error
signal between reference speed and actual speed of the permanent magnet
synchronous motor drive. In this torque dynamics state the error is nearly maintained
to zero as shown in Figure 7.18 below.

Figure 7.18 Error Signal between Reference Speed and Actual Speed

7.3.5 Vd and Vq

The below Figure 7.19 waveforms are the simulation result of V d voltage and
Vq voltage variation which is implemented in MATLAB/SIMULINK. Here two
phases Vafa and Vbeta voltages are converted into two phase Vd and Vq waveforms i.e.

63
Vafa and Vbeta is converted into Vd and Vq (which is known as park’s transformation or
rotating frame of reference) Under torque dynamics.

Figure 7.19 (a)Vd and (b) Vq

7.3.6 Iq, Id, Id-ref and Iq-ref

The below Figure 7.20 and 7.21 waveforms are the simulation result of I q, Id,
Id-ref and Iq-ref currents which is implemented in MATLAB/SIMULINK. In this I d-ref, Iq-
ref is the reference signals in this I d-ref are maintained at zero. As shown in Figure 7.21.
Thus the actual signals which are generated in the control analysis are as shown in the
Figure 7.20 are following the reference signals. This is under torque dynamics.

Figure 7.20 (a) Iq current and (b) Id current

64
 Figure 7.21 (a) Id-ref and (b) Iq-ref

7.4 CASE 4

In this both speed and the torque are varied. Hence below are the waveforms
which are obtained from simulation in MATLAB
7.4.1 Actual Speed And Reference Speed

Figure 7.22 (a) Actual speed and (b) Reference speed

The above Figure 7.22 waveforms are the simulation result of speed variation
which is implemented in MATLAB/SIMULINK. It give information about the

65
variation of speed with respective to the reference speed. In this the actual speed
follows the reference speed. Hence speed or load is varied at any point of time. In this
Figure 5.35 the speed is changed to 500rpm at some point of time.

7.4.2 Reference Torque and Actual Torque

The below Figure 7.23 waveforms are the simulation result of torque variation
which is implemented in MATLAB/SIMULINK. It give information about the
variation of torque with respective to the reference torque. In this the actual torque
follows the reference torque. Torque is varied at any point of time.

Figure 7.23 (a) Reference torque and (b) Actual torque

7.4.3 Error Signal between Reference Speed and Actual Speed

The below Figure 7.24 waveform is the simulation result of error variation
which is implemented in MATLAB/SIMULINK. This is the waveform of the error
signal between reference speed and actual speed of the permanent magnet
synchronous motor drive. In this speed-torque dynamics state the error is nearly
maintained to zero as shown in Figure 7.24 below.

66
 Figure 7.24 Error Signal between Reference Speed and Actual Speed

7.4.4 Iq Current, Id Current and Id_ref, Iq_ref

Figure 7.25 (a) Iq current and (b) Id current

67
 Figure 7.26 (a) Id-ref current and (b) Iq-ref current

The above Figures 7.25 and 7.26, waveforms are the simulation result of I q, Id,
Id-ref and Iq-ref currents which is implemented in MATLAB/SIMULINK. In this I d-ref, Iq-
ref is the reference signals in this I d-ref are maintained at zero. As shown in Figure 7.26.
Thus the actual signals which are generated in the control analysis are as shown in the
Figure 7.25 are following the reference signals. This is under speed-torque dynamics.

7.4.5 Vd Voltage and Vq Voltage

68
 Figure 7.27 (a) Vd Voltage and (b) Vq Voltage

The above Figure 7.27 waveforms are the simulation result of V d voltage and
Vq voltage variation which is implemented in MATLAB/SIMULINK. Here two
phases Vafa and Vbeta voltages are converted into two phase Vd and Vq waveforms i.e.
Vafa and Vbeta is converted into Vd and Vq (which is known as park’s transformation or
rotating frame of reference) Under speed-torque dynamics.

7.5 DISCUSSION OF RESULTS

These simulation results show the drive operation is satisfactory. In order to

steady the speed torque characteristics of PMSM we are considering 4 cases.
In case 1, the speed and torque are kept constant. The speed is assumed to be
500rpm and the torque is considered to be 0.5 N-m. The actual and reference speed
along with the error and actual and reference torque waveforms are as shown in the
Figures. From the waveforms we can observe that the actual torque and speed follows
the reference torque and speed. Hence under steady state there is no change in actual
torque and speed i.e. speed and torque are not varied at any point of time and the error
between the reference speed and actual speed is almost maintained to zero.

In case 2, the speed is varied and the torque is kept constant. The speed is
varied upto 800 rpm in steps and the torque is taken to be 0.5 N-m. We can observe
the distortions in the waveform of the actual speed at 0.5 since torque and speed are
inversely proportional.

In case 3, the speed is kept constant and the torque is varied. The speed is
assumed to be 500 rpm and the torque is varied up to -0.2 N-m. From the waveforms
we can observe that when torque is high the speed is low and when torque is low, the

69
speed is high since speed and torque are inversely proportional to each other and the
error between the reference speed and actual speed is almost maintained to zero.

In case 4, both speed and torque are varied. The torque is varied upto -0.2 and
speed is varied upto 800 rpm and in this case we can observe that the error between
the actual and reference speed is almost maintained to zero.

CHAPTER VIII
CONCLUSIONS AND FUTURE SCOPE

70
 8 CONCLUSIONS AND FUTURE SCOPE

8.1 CONCLUSION

Thus this thesis explained the mathematical equations related to the application
of the vector control of PMSM. The SIMULINK block of vector control in PMSM
is presented. The simulation results examined the implementation of the vector
control in PMSM.
A detailed SIMULINK model for vector control of PMSM has been
developed. SIMULINK has been chosen from several simulation tools because of its
flexibility. Mathematical models can be easily incorporated in the simulation and the
present numerous tools boxes and support guides simplified the simulation of large
system compared to other software. SIMULINK is capable of showing real time
results with reduced simulation time and debugging.

71
 It is shown in the simulation results of vector control in PMSM that the
vector control implementation in PMSM is having less level of speed and toque
ripples and at the same time maintaining the good torque response.

8.2 FUTURE SCOPE

Future works needed to be carried out to improve the performance of the

motor. A vector-controlled PMSM drive with a continually on-line learning hybrid
neural-network model-following speed controller is one of the advancement in vector
control of PMSM in which the speed controller produces rapid, robust performance
and accurate response to the reference model regardless of load disturbances or
PMSM parameter variations. The other is vector control of permanent magnet
synchronous motor with surface magnet using artificial neural networks; the use of
ANN makes the drive system robust, accurate and insensitive to parameter variations.
The hardware implementation of the above drive system can also be done.

72
 APPENDICES

APPENDIX 1

MATLAB is a high-performance language for technical computing. It

integrates computation, visualization, and programming in an easy to use environment

where problems and solutions are expressed in familiar mathematical notations.

MATLAB is an interactive system whose basic data element is an array that

does not require dimensioning. This allows you to solve many technical problems,

especially those with matrix and vector formations. , in a fraction of the time it would

take to write a program in a scalar non interactive language such as C or Fortran. The

name MATLAB stands for matrix laboratory. MATLAB was originally written to

provide easy access to matrix software developed by the LINPACK and EISPACK

projects. Today a MATLAB engines incorporates the LAPACK and BLAS libraries,

73
embedding the state of the art in software for matrix computation. MATLAB has

evolved over a period of years with input from many users. In university

environments, it is the standard instructional tool for introductory and advanced

courses in mathematics, engineering and science. In industry, MATLAB is the tool of

choice for high productivity research, development, and analysis. MATLAB features

a of add on application specific solutions called toolboxes. Very important to most

users of MATLAB, toolboxes allow you to learn and apply specialized technology.

Toolboxes are comprehensive collections of MATLAB functions (M-files) that

extended the MATLAB environment to solve particular classes of problems. Areas

in which toolboxes are available include signal processing, control systems, neural

networks, fuzzy logic, wave lets, simulation, and many others.

APPENDIX 2
MATLAB CODE

1.
2.
3. clear
4.
5. Vph = 60/sqrt(3)/2; % [V], peak value, phase
voltage
6. Vs = Vph/sqrt(2); % RMS phase voltage
7.
8. nrat = 1000; % [rpm], rate speed
9.
10. Npp = 4/2; % Number of pole pairs
= 1/2 number of poles
11.
12. Omegae_rat = 2*pi*nrat/60*Npp;
i. % Rated electrical angular
frequency

74
13.
14.
15. Lndmpm = Vph/Omegae_rat;% [Web.turns], rotor
peak PM flux linkage
16.
17. J = 1*1e-4*1; % Motor inertia [J.m^2]
18. % _____________________________________________________
19.
20. Vs = 11.25; % Specified terminal
voltage, RMS value
21. Vph = sqrt(2)*Vs; % the corresponding
peak phase voltage
22.
23.
24. fs = 5e3; % switching frequency
25.
26. % Initializing the Simulink model
27. Omegae_ini = 0; % Initial motor shaft
speed, electrical value, [rad]
28. Lndd_ini = Lndmpm; % This means at t=0,
theta=0 and N-pole aligned with d-axis
29.
30. Omegae = 2*pi*5;
31.
32. % ------------------------------------------------------
33. % -------------------------------------------------------
34.
35. Rs = 3.4; % [Om]
36. Lls = 1.1e-3; % [H]
37. Lmd = 11e-3; % [H]
38. Lmq = Lmd;
39.
40. % _______________________________________________________
41. % Lmq = 0.6*Lmd;
42. % ______________________________________________________
43.
44. Ld = Lmd+Lls;
45. Lq = Lmq+Lls;
46.
47. % -------------------------------------------------------
48. % -------------------------------------------------------
49.
50. Iphmax = 2.5;
51. Omeg_vec = linspace(Omegae_rat/200,
Omegae_rat,5);

75
52.
53. Omeg_const = Omeg_vec(5); % Parameter needed by
'PMSMmotor_setp3_discrete.mdl'
54.
55. % -------------------------------------------------------
56. % -------------------------------------------------------
57.
58. switch 4
59. case 1
60. % A set of parameters that do not work! (Oscillation in
S.S.)
61. Kp1 = 0.8;
62. Ki1 = 818;
63. Kp2 = 8.16;
64. Ki2 = 7645;
65. Kp3 = 26.05;
66. Ki3 = 8014;
67. case 2
68. % A set of parameters that work well in continuous states
but resulting S.S. oscillation in
69. % discrete control system
70. Kp1 = 3.5;
71. Ki1 = 10e-3; %10e-3;
72. Kp2 = 3.5;
73. Ki2 = 3.5e3/10;
74. Kp3 = 20/2; %6.8;
75. Ki3 = 4.8e3/10; %3.8e3;
76. case 3
77. % Used for testing PI parameters one by one in
'PMSMmotor_setp3_discrete.mdl', which can
78. % both work in continuous and discrete systems
79. Kp1 = 0.025;
80. Ki1 = 0.8;
81. Kp2 = 3.0;
82. Ki2 = 5.0e3;
83. Kp3 = 5.5;
84. Ki3 = 3.0e3;
85. case 4
86. % These seem to be the best parameters!
87. Kp1 = 0.025;
88. Ki1 = 0.8;
89. Kp2 = 3.0;
90. Ki2 = 5.0e3;
91. Kp3 = 5.5;
92. Ki3 = 3.0e3;

76
93. case 5
94. % Parameters can let the continuous model work when there
is none of the 'saturation' units
95. Kp1 = 0.1;
96. Ki1 = 1.0;
97. Kp2 = 0.1;
98. Ki2 = 1.0;
99. Kp3 = 0.1;
100. Ki3 = 1.0;
101. case 6
102. Kp1 = 0.025;
103. Ki1 = 0.8;
104. Kp2 = 3.0;
105. Ki2 = 0;
106. Kp3 = 5.5;
107. Ki3 = 0;
108.
109. case 6
110.
111. end
112.
113. % _______________________________________________________
114. % Parameters used for PI anti wind-up function
115.
116. T_pi_delay = 1/fs;
117. SaturationVal_1 = 50;
118. SaturationVal_2 = 2;
119. SaturationVal_3 = 2;
120. % ______________________________________________________
121.
122. % -------------------------------------------------------
123. % -------------------------------------------------------
124.
125. if exist('t_sim')
126. switch 1
127. case 1
128. plot(t_sim,n_ref ,'Color',colorvec(mod(6-
1,6)+1,:),'Marker','none'); hold on; grid on
129.
130. colidx = 3;
plot(t_sim,n_meas,'Color',colorvec(mod(colidx-
1,6)+1,:),'Linestyle',':'); hold on; grid on
131. case 2
132. end
133. end

77
 134.
135.
136.
137.

BIBLIOGRAPHY

[1] T. Sebastian, G. Slemon, and M. Rahman, "Modeling of permanent magnet

synchronous motors," Magnetics, IEEE Transactions on, vol. 22, pp. 1069-1071,

1986.

[2] T. M. Jahns, G. B.Kliman, and T. W. Neumann, "Interior Permanent-Magnet

Synchronous Motors for Adjustable-Speed Drives," Industrial Applications, IEEE

Transactions on, vol. IA-22, pp. 738-746, 1986.

78
[3] P. Pillay and R. Krishnan, "Modeling of permanent magnet motor drives,"

Industrial Electronics, IEEE Transactions on, vol. 35, pp. 537-541, 1988.

[4] P. Pillay and R. Krishnan, "Modeling, simulation, and analysis of permanent-

magnet motor drives. I. The permanent-magnet synchronous motor drive," Industry

Applications, IEEE Transactions on, vol. 25, pp. 265-273, 1989.

[5] B. K. Bose, Modern power electronics and AC drives: Prentice Hall, 2002

[6] A. H. Wijenayake and P. B. Schmidt, "Modeling and analysis of permanent

magnet synchronous motor by taking saturation and core loss into account," 1997.

[7] K. Jang-Mok and S. Seung-Ki, "Speed control of interior permanent magnet

synchronous motor drive for the flux weakening operation," Industry Applications,

IEEE Transactions on, vol. 33, pp. 43-48, 1997.

[8] Weera Kaewjind and Mongkol Konghirun “Vector Control Drive of Permanent

Magnet Synchronous Motor Using Resolver Sensor” ECTI transactions on electrical

eng., electronics, and communications vol.5, no.1 February 2007.

79