Dynamic Model Analysis of Three Phase Induction Motor Using Matlab Simulink
Dynamic Model Analysis of Three Phase Induction Motor Using Matlab Simulink
ISSN 2229-5518
Abstract This paper reproduce the dynamic model of three phase induction motor derived from detailed mathematical equation in
MATLAB/Simulink and study the transient performance as well as steady state response of machine including the constant moment of
inertia and load torque for the external system. Dynamics generalized approach [1] of digital computer simulation of three phase induction
machine is universally applicable method to propose for simulating a dynamic behavior of a large, non linear, time varying physical system.
The model can account for unbalanced on both sides of the air gap provided the structure of winding remains unaltered. The simulation in
the three-phase Induction motor is with variations in moment of inertia [2]. Transient performance of any electrical machine is greatly
affected by sudden changes in its supply system operating speed, shaft load including any variation of moment of inertia due to gear
arrangement application. Rotor reference frame is used for the simulation study of this machine. MALTLAB/ Simulink based model is
adopted to compare a transient performance of induction machine including main flux saturation with and without the moment of inertia.
Simulated results have been compared and verified with experimental results on a test machine setup. During the development of dynamic
models of induction motors, most of the researchers neglected the effect of magnetic saturation and assumed inductances to be constant
[2]. The computer simulation for these various modes of operation is conveniently obtained from the equations which describe the
symmetrical induction machine in an arbitrary reference frame.
Index Terms d-q axis, Dynamic model, Leakage inductance, Magnetic saturation, Reference frame, Transient analysis,
IJSER
1 INTRODUCTION
IJSER
as
dFds R
=
V rs iabcs + pabcs = b [vds + e Fqs + s ( Fmd + Fds )] (7)
(1) b
abcs
dt X ls
=
V rr iabcr + pabcr
( r )
abcr
For a magnetically linear system, the flux linkages may be ex- dFdr R
= b [vdr + e Fqr + r ( Fmd Fdr )] (8)
pressed as dt b X lr
abcs Ls Lsr iabcs dFqr (e r )
= ( L )T (2) R
Lr iabcr =b [vqr Fdr + r ( Fmq Fqr )] (9)
abcr sr dt b X lr
The voltage equations expressed in terms of machine variables For a Squirrel Cage Induction Machine as in the case of this
referred to the stator windings may be written as paper, V qr and V dr are set to zero. The modeling equations of a
Vqs Rs + Ls p r Ls Lm p r Lm iqs Squirrel cage Induction motor in state space become:
Vds = r Ls Rs + Ls p r Lm Lm p ids (3) dFqs e Rs X *ml X ml *
0 Lm p 0 Rr + Lr p 0 iqr = b [vqs Fds + ( Fqr + ( 1) Fqs )]
Rr + Lr p idr
dt b X ls X lr X ls
0 0 Lm p 0
(10)
The torque and rotor speed are related by
dFds R X X * *
-200
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
time (in Second)
IJSER
Fig.5.1 Stator Current in q- axis versus time curves (iqs vs. time)
200
150
mag (in Ampere)
100
50
-50
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
time (in Second)
Fig.5.2 Stator Current in d- axis versus time curves (ids vs. time)
Stator Current in d-q-axis Iqds
180
160
140
120
100
80
mechanical sub-model to yield the rotor angular speed. In electrical time (in Second)
Table-5.2
Fig 5.6 Speed - Torque Characteristics
Absolute Peak value and Average Value of rotor current in
d-q axis:
Type Rotor Current Rotor Current Rotor Current
in q-axis (I qr ) in d-axis (I dr ) in q-d axis 6 RESULT & DISCUSSION
(amp) (amp) (I qdr ) (amp) To analyze the performance of induction motor with
Average 0.0984 2.913 1.983 Matlab/Simulink, the model is simulated when pure sinusoi-
dal three phase voltage is applied. In this paper, the dynamic
Value
responses of speed and stator currents for the starting process
with constant torque applied and a constant command flux are
IJSER
Absolute 173.33 215 160.4594
shown in simulation results. The peak values and the average
Peak Value values of speed and torque are shown in above. It can be said
that the rotor angular velocity is one fourth () of the rated
angular speed.
Speed-Time Characteristic In this paper, the comparison between simulink imple-
mented model and Matlab/Simulink model in explained. The
350
300
simulink implemented model means which is implemented by
250 us based on detail mathematical equation, as discussed in
above. On the other hand the Matlab/Simulink demo model
Speed (in rad/sec)
200
150
means the model which is already exists in the
MATLAB/Simulink software. For this comparison of three
phases IM, same parameters for both models should be taken,
100
50
which are already mentioned in table-6.1. The fundamental
0 supply voltage is applied to the both models for this compari-
-50
son.
Table-6.1
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Time (in Sec)
400
0.754
Mutual inductance ()
300
200
26.13
100
Load torque (N-m) 11.9
0 Moment of inertia (Kg-m2) 0.089
-100
-200
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Matlab simulink model of three phase asynchronous machine
Time (in Sec) is taken from the MATLAB/Simulink software. After getting
the model, the parameters are putted, which are already used
Fig5.5 Speed Characteristics to the implemented model. These detail natures of comparison
IJSER 2016
http://www.ijser.org
International Journal of Scientific & Engineering Research, Volume 7, Issue 3, March-2016 576
ISSN 2229-5518
are shown in below figure. Here firm line denotes the Matlab
implemented simulink model, whereas, dotted line denotes
the Matlab / Simulink demo model.
150
Current (in Ampere)
100
50
700
Te(demo)
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Time (in Second) 600
Te(main)
500
Fig.6.1 Stator Current in q- axis versus time curves (iqs vs. time)
200
200
Ids(demo)
150
IJSER
Iqs(demo) 100
100 0
Current (in Ampere)
50 -100
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Time (in Second)
0
-100
-150
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Time (in Second)
Fig.6.2 Stator Current in d- axis versus time curves (ids vs. time)
Stator Current in q-d-axis Characteristics
150
Iqds(demo)
Iqds(main)
100
Current (in Ampere)
50
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Time (in Second) CONCLUSION
In this paper, a simulation approach was proposed to analyze
Fig.6.3 Stator Current in dq- axis versus time curves (iqds vs. the transient behavior as well as steady state behavior of three
time) phase induction motor. Electrical and mechanical parts of in-
duction motor were implemented by detail mathematical
equations which are explained. The dynamic model is pre-
sented in this paper has shown very accurate simulation re-
sults. The main disadvantage of this implementing IM model
is the value of the rotor angular velocity is times of actual
rated rotor angular velocity.
IJSER 2016
http://www.ijser.org
International Journal of Scientific & Engineering Research, Volume 7, Issue 3, March-2016 577
ISSN 2229-5518
FUTURE WORK [9] P.C.Krause,Analysis of Electrical Machinery 1rd ed. Vol.2, 4, 5,
McGraw-Hill: New York, 1986.
Future work will be on the simulation of induction motor [10] Burak Ozpineci, Leon M. Tolbert,Simulink implementation of In-
model by taking the effect of saturation of flux in supply side duction Machine Model-A Modular Approach IEEE, 0-7803-7817-
and variable load torque, with and without MOI and friction 2/03/2003.
co-efficient. Also the work will be on operation of the Induc- [11] M.L.De Aguiar, M.M.Cadd,The concept of complex transfer func-
tion motor model to simulate for generator responses at a var- tions applied to the modeling of induction motors, Power engineer-
iable load torque. Apart from this, we will consider the effect ing Society Winter Meeting, 2000, pp.387-391.
of higher order odd harmonics on induction motor. The total [12] A. Dumitrescu, D. Fodor, T. Jokinen, M. Rosu, S. Bucuren-
MATLAB/Simulink model will be implemented by the State cio,Modeling and simulation of electric drive systems using
space equation of the IM. Matlab/Simulink environments, International Conference on Elec-
tric Mechines and Drives (IEMD), 1999, pp.451-453.
ACKNOWLEDGEMENT [13] Scott. Wade, M.W. dunnigor, B.W. Williams, Modeling and simula-
It is a pleasant task to express our gratitude to all those who tion of induction machine vector control with rotor resistance identi-
have accompanied and helped us in this work. fication,IEEE Transactions on power electronics, Vol.12, no.
First and foremost, we take this opportunity to express deep 3,May1997,pp.495-506.
sense of gratitude to our Guide (late) Dr. Prof. A. K. Mukha- [14] H. Le-Huy, Modeling and Simulation of electric drives using Matlab
padhaya, Department of Electrical Engineering, Narula Insti- /Simulink and power system Blockset, The 27th Annual conference
tute of Technology, Kolkata- India, for his invaluable sugges- of the IEEE Industrial Electronics Society (IECON01), Denver/ Colo-
tions and encouragement throughout the project which helped rado, pp.1603-1611.
[15] Paul C. Krause, O. Wasynczuk and S. D. Sudhoff. 2004. Analysis of
us a lot to improve this project work. We are also indebted to
Electric Machinery and Drive Systems. IEEE Press Series on Power
Principal and HOD-EE of Narula Institute of Technology, Kol-
Engineering, John Wiley and Sons Inc. Publication.
kata- India their kind help and co-operation during this work.
[16] K.L. Shi, T.F. Chan, Y.K Wong and S.L Ho department of electrical
Our sincere thanks to all our family members for their moral
IJSER
engineering, Hong Kong Polytechnic University, Modeling and
support and encouragement, without which, the work would
Simulation of the Three phase Induction motor using Simulink.
not have been possible.
[17] R. Krishnan. 2007. Electric Motor Drives. Pearson Prantice Hall.
Last but not least, we extend our thanks and appreciation to
[18] S.I. Moon and A. Keyhani. 1994. Estimation of Induction Machine
our friends, colleagues, batch-mates and everyone who have
Parameters from Standstill Time-Domain Data. IEEE Transactions on
helped us directly or indirectly to get this work done.
Industry Applications. 30(6): 1609-1615.
REFERENCES
[1] Sayeed Nurul Ghani, The Digital computer simulation of three
phase induction machine dynamics a generalized approach. 1990.
Improved models for the simulation of deep bar induction motor.
IEEE Trans. on Energy Conversion. 5(2): 393-400. June.
[2] K.S Sandhu and Vivek Pahwa, simulation study of Three Phase
induction motor with variations in moment of inertia, ARPN
Journal of Engineering and Applied Sciences, vol. 4, No. 6,
August 2009, pp-72-77.
[3] Paul C Krause, O Wasynczuk and S.D. Sudhoff. 2004,Analysis of
electric Machinery and Drive systems.IEEE press series on Power
Engineering John Wiley and Sons Inc Publication.
[4] J.R. Willis, G.J. Brock and J.S. Edmonds. 1989. Derivation of Induc-
tion Motor Models from standstill Frequency response tests.IEEE
transactions on energy conversion. Vol.4, Dec. pp.608-613.
[5] M. Ikeda, and T. Hiyama, Simulation studies of the transients of
squirrel cage induction motors, IEEE trans. Energy Conv. Sec.,
Vol.22, No.2, pp.233-239, June 2007.
[6] C.F. Landy, the prediction of transient torques produced in induc-
tion motors due to rapid reconnection of the supply, SAIEE. Trans,
Vol.112, No.7, 1965, pp.1348-1354.
[7] A.K. chottopadhyay,Digital computer Simulation of an adjustable
speed induction motor drive with a cycloconverter type thyristor
commutator in the rotor ,IEEE Trans Industrial electronic,vol.2,pp-
86-92,Feb 1976
[8] O.I. Okoro,Matlab Simulation of induction machine with saturable
leakage and magnetizing inductances.PJSTS, Vol.5, No.1, 2003, pp-
398-400
IJSER 2016
http://www.ijser.org