Ashutosh 2
Ashutosh 2
Ashutosh 2
1
PG Scholar, S R Engineering College, Warangal, Telangana, India. Pin: 506371.
2
Department of Electrical and Electronics Engineering, S R University, Warangal, Telangana,
India. Pin: 506371.
3
Department of Electrical and Electronics Engineering, Sumathi Reddy Institute of Technology,
Warangal, Telangana, India. Pin: 506371
Email: Saikumar.vamshi@hotmail.com
Abstract: in this paper, the mathematical model of three-phase and five-phase induction motors
has been developed using PARK’S Transformation and implemented in MATLAB SIMULINK.
This Mathematical modelling plays a key role in analyzing the system and its algorithms. Both
the models have been tested under loaded and No load condition to find a comparative analysis.
The main focus of this paper is to determine the increased capability of motor in switching from
three phase to five phase.
1. INTRODUCTION
With Rapid spike in industrialization, there has been a huge demand for heavy and robust equipment.
Mostly in industries, the rating of induction motor required is very high. As the structure of induction
motors have been rugged and most suitable for rough environments, several research has been done on
them leading to buildout of multi-phase induction motors such as five-phase and six-phase induction
motors.
And one of the main reason for developing multi-phase induction motors is, notably in varying speed
electrical drives the interlinking circuit between supply side and motor is inverter.
As such increasing of the number of limbs in inverter is not restricted technically, and thus any number
of phases can be acquired from an inverter. So this has paved path for improvement of multi-phase
induction motors and their drive systems.
The various applications of five phase induction motors is rolling mills, ship propulsion, Electric
vehicles, Hybrid electric vehicles [1], Wind energy systems, Traction.
The main reason for switching from three-phase induction motor to five-phase induction motor is
because of the following advantages.
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution
of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Published under licence by IOP Publishing Ltd 1
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
Mathematical model is the very fundamental and vital step to get an efficient and reliable electrical
machine. More the accurate mathematical model, more the efficiency and reliability. Despite the fact
that mathematical modeling of induction motor has very old history and there is still lot of work needed
to be done, to get a better replica of induction motor. Different models has been proposed for different
applications such as for prediction of behaviour of system, analysis of the system, fault analysis, flux
level control. A reliable motor depends on its accurate design and perfect parameter estimation.
The mathematical machine modeling in state space form is significant for transient analysis. Despite the
fact that rotating frame model is typically opted, but the model of stationary frame has been using now
extensively. The parameters in the model can be selected as currents, fluxes, or a concoction of both.
Given below the derived state space equations for motor in rotating frame with flux linkages as the
primary variables. The five phase voltages are assumed as
2𝜋
𝑉𝑏𝑠 = 𝑉𝑚 sin (𝜔𝑒 𝑡 − 5
) (2)
4𝜋
𝑉𝑐𝑠 = 𝑉𝑚 sin (𝜔𝑒 𝑡 − ) (3)
5
6𝜋
𝑉𝑑𝑠 = 𝑉𝑚 sin (𝜔𝑒 𝑡 − ) (4)
5
8𝜋
𝑉𝑒𝑠 = 𝑉𝑚 sin (𝜔𝑒 𝑡 − 5
) (5)
Where, 𝑉𝑚 is the peak value of terminal voltage 𝜔𝑒 is the supply frequency in radians per second
These stationary reference frame variables is converted into two phase stationary reference frame
variables using decoupling transformation matrix [4]. The decoupling transformation matrix is given
below.
𝑠
𝑉𝑑𝑠 𝑉𝑎𝑠
1 𝑐𝑜𝑠𝛼 𝑐𝑜𝑠2𝛼 𝑐𝑜𝑠3𝛼 𝑐𝑜𝑠4𝛼
𝑉𝑞𝑠𝑠 0 − 𝑠𝑖𝑛𝛼 − 𝑠𝑖𝑛2𝛼 − 𝑠𝑖𝑛3𝛼 − 𝑠𝑖𝑛4𝛼 𝑉 𝑏𝑠
2
𝑉𝑥𝑠𝑠 = 1 𝑐𝑜𝑠3𝛼 𝑐𝑜𝑠6𝛼 𝑐𝑜𝑠9𝛼 𝑐𝑜𝑠12𝛼 𝑉𝑐𝑠 (6)
5
𝑠 0 − 𝑠𝑖𝑛3𝛼 − 𝑠𝑖𝑛6𝛼 − 𝑠𝑖𝑛9𝛼 − 𝑠𝑖𝑛12𝛼 𝑉
𝑉𝑦𝑠 𝑑𝑠
𝑠 [ 0.5 0.5 0.5 0.5 0.5 ] [ 𝑉𝑒𝑠 ]
[𝑉0𝑠 ]
2𝜋
Where 𝛼 = 5
𝑠
𝑉𝑑𝑠 and 𝑉𝑞𝑠𝑠 are the 𝑑 𝑠 and 𝑞 𝑠 axis stator fundamental voltages. 𝑉𝑥𝑠 ,𝑉𝑦𝑠 ,𝑉0𝑠 are the extra components of
stator considered so as it will be useful during inverse transformation from dq to abc.
2
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
These two phase stationary reference frame variables are transformed into synchronously rotating
reference frame, ignoring the zero sequence components as shown below [5].
𝑠
𝑉𝑞𝑠 = 𝑉𝑞𝑠𝑠 𝑐𝑜𝑠𝜃𝑒 − 𝑉𝑑𝑠 𝑠𝑖𝑛𝜃𝑒 (7)
𝑠
𝑉𝑑𝑠 = 𝑉𝑞𝑠𝑠 𝑠𝑖𝑛𝜃𝑒 + 𝑉𝑑𝑠 𝑐𝑜𝑠𝜃𝑒 (8)
Figure 1. conversion of stationary reference frame variables to synchronously rotating frame variables
Where, 𝑉𝑑𝑠 and 𝑉𝑞𝑠 are the 𝑑 and 𝑞 axis synchronously rotating voltages, which rotate at synchronous
speed 𝜔𝑒 with respect to 𝑑 𝑠 and 𝑞 𝑠 axis and the angle 𝜃𝑒 = 𝜔𝑒 𝑡.
The block diagram for this conversion is shown in Figure 1.
𝑑
𝑉𝑞𝑠 = 𝑅𝑠 𝑖𝑞𝑠 + 𝑑𝑡 𝜓𝑞𝑠 + 𝜔𝑒 𝜓𝑑𝑠 (9)
𝑑
𝑉𝑑𝑠 = 𝑅𝑠 𝑖𝑑𝑠 + 𝑑𝑡 𝜓𝑑𝑠 − 𝜔𝑒 𝜓𝑞𝑠 (10)
The last terms in above equations is defined as speed emf due to rotation of the axes, that is,
when 𝜔𝑒 = 0, the equations reverse back to stationary form.
𝑑
𝑉𝑞𝑟 = 𝑅𝑟 𝑖𝑞𝑟 + 𝑑𝑡 𝜓𝑞𝑟 + (𝜔𝑒 − 𝜔𝑟 )𝜓𝑑𝑟 (11)
𝑑
𝑉𝑑𝑟 = 𝑅𝑟 𝑖𝑑𝑟 + 𝑑𝑡 𝜓𝑑𝑟 − (𝜔𝑒 − 𝜔𝑟 )𝜓𝑞𝑟 (12)
3
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
Where, 𝜔𝑏 is the base frequency of the machine. Substituting Equations (19) to (22) in equations (9) to
(12), We get,
1 𝑑𝐹𝑞𝑠 𝜔𝑒
𝑉𝑞𝑠 = 𝑅𝑠 𝑖𝑞𝑠 + + 𝐹 (23)
𝜔𝑏 𝑑𝑡 𝜔𝑏 𝑑𝑠
1 𝑑𝐹𝑑𝑠 𝜔
𝑉𝑑𝑠 = 𝑅𝑠 𝑖𝑑𝑠 + 𝜔 − 𝜔𝑒 𝐹𝑞𝑠 (24)
𝑏 𝑑𝑡 𝑏
4
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
𝐹𝑑𝑠 −𝐹𝑑𝑚
𝑖𝑑𝑠 = 𝑋𝑙𝑠
(35)
𝐹𝑑𝑟 −𝐹𝑑𝑚
𝑖𝑑𝑟 = (36)
𝑋𝑙𝑟
𝑋∗ 𝑋∗
𝐹𝑞𝑚 = 𝑋𝑚 𝐹𝑞𝑠 + 𝑋𝑚 𝐹𝑞𝑟 (38)
𝑙𝑠 𝑙𝑟
Similarly,
𝑋∗ 𝑋∗
𝐹𝑑𝑚 = 𝑋𝑚 𝐹𝑑𝑠 + 𝑋𝑚 𝐹𝑑𝑟 (39)
𝑙𝑠 𝑙𝑟
5
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
∗ 1
Where, 𝑋𝑚 = 1 1 1
(40)
( + + )
𝑋𝑚 𝑋𝑙𝑠 𝑋𝑙𝑟
𝑅 1 𝑑𝐹𝑞𝑠 𝜔
𝑉𝑞𝑠 = 𝑋 𝑠 (𝐹𝑞𝑠 − 𝐹𝑞𝑚 ) + 𝜔 + 𝜔𝑒 𝐹𝑑𝑠 (41)
𝑙𝑠 𝑏 𝑑𝑡 𝑏
𝑅 1 𝑑𝐹𝑑𝑠 𝜔
𝑉𝑑𝑠 = 𝑋 𝑠 (𝐹𝑑𝑠 − 𝐹𝑑𝑚 ) + 𝜔 − 𝜔𝑒 𝐹𝑞𝑠 (42)
𝑙𝑠 𝑏 𝑑𝑡 𝑏
𝑑𝐹𝑞𝑠 𝜔𝑒 𝑅𝑠
= 𝜔𝑏 [𝑉𝑞𝑠 − 𝐹 − (𝐹 − 𝐹𝑞𝑚 )] (45)
𝑑𝑡 𝜔𝑏 𝑑𝑠 𝑋𝑙𝑠 𝑞𝑠
𝑑𝐹𝑑𝑠 𝜔 𝑅
𝑑𝑡
= 𝜔𝑏 [𝑉𝑑𝑠 + 𝜔𝑒 𝐹𝑞𝑠 − 𝑋 𝑠 (𝐹𝑑𝑠 − 𝐹𝑑𝑚 )] (46)
𝑏 𝑙𝑠
6
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
7
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
Development of torque by interaction of airgap flux and rotor mmf is expressed in generalized form as,
𝑛 𝑃
𝑇𝑒 = ( ) ̅̅̅̅
𝜓𝑚 × 𝐼̅𝑟 (49)
2 2
5 𝑃
𝑇𝑒 = 2 ( 2 ) (𝜓𝑑𝑠 𝑖𝑞𝑠 − 𝜓𝑞𝑠 𝑖𝑑𝑠 ) (50)
5 𝑃 1
𝑇𝑒 = 2 ( 2 ) (𝜔 ) (𝐹𝑑𝑠 𝑖𝑞𝑠 − 𝐹𝑞𝑠 𝑖𝑑𝑠 ) (51)
𝑏
2 𝑑𝜔𝑟
𝑇𝑒 = 𝑇𝐿 + 𝑃 𝐽 𝑑𝑡
(52)
Using equations (33) to (36), those synchronously rotating reference frame currents are converted into
stationary frame two phase variables
8
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
These stationary two phase variables are converted into five phase variables using the below conversion
technique.
𝑖𝑎 1 0 1 0 1 𝑖𝑑
𝑖𝑏 cos 𝛼 sin 𝛼 cos 2𝛼 sin 2𝛼 1 𝑖𝑞
2
𝑖𝑐 = √ × cos 2𝛼 sin 2𝛼 cos 4𝛼 sin4 𝛼 1 𝑖𝑥 (55)
5
𝑖𝑑 cos 3𝛼 sin 3𝛼 cos 6𝛼 sin 6𝛼 1 𝑖𝑦
[ 𝑖𝑒 ] [cos 4𝛼 sin 4𝛼 cos 8𝛼 sin 8𝛼 1] [ 𝑖0 ]
3. SIMULATION MODEL
The simulation model of five phase induction motor designed is shown in Fig.11
9
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
2𝜋
𝑉𝑏𝑠 = 𝑉𝑚 sin (𝜔𝑒 𝑡 − 3
) (57)
4𝜋
𝑉𝑐𝑠 = 𝑉𝑚 sin (𝜔𝑒 𝑡 − 3
) (58)
Where, 𝑉𝑚 is the peak value of terminal voltage 𝜔𝑒 is the supply frequency in radians per second.
These stationary reference frame variables is converted into two phase stationary reference frame
variables using decoupling transformation matrix. The decoupling transformation matrix is given below.
𝑠
𝑉𝑑𝑠 1 𝑐𝑜𝑠𝛼 𝑐𝑜𝑠2𝛼 𝑉𝑎𝑠
2
[ 𝑉𝑞𝑠𝑠 ] = 3
[ 0 − 𝑠𝑖𝑛𝛼 − 𝑠𝑖𝑛2𝛼 ] [𝑉𝑏𝑠 ] (59)
𝑠
𝑉𝑛𝑠 0.5 0.5 0.5 𝑉𝑐𝑠
2𝜋
Where 𝛼 = 3
𝑠
𝑉𝑑𝑠 and 𝑉𝑞𝑠𝑠 are the 𝑑 𝑠 and 𝑞 𝑠 axis stator fundamental voltages. 𝑉𝑛𝑠 is the zero sequence component of
stator.
These two phase stationary reference frame variables are transformed into synchronously rotating
reference frame, ignoring the zero sequence components as shown below.
𝑠
𝑉𝑞𝑠 = 𝑉𝑞𝑠𝑠 𝑐𝑜𝑠𝜃𝑒 − 𝑉𝑑𝑠 𝑠𝑖𝑛𝜃𝑒 (60)
𝑠
𝑉𝑑𝑠 = 𝑉𝑞𝑠𝑠 𝑠𝑖𝑛𝜃𝑒 + 𝑉𝑑𝑠 𝑐𝑜𝑠𝜃𝑒 (61)
Where, 𝑉𝑑𝑠 and 𝑉𝑞𝑠 are the 𝑑 and 𝑞 axis synchronously rotating voltages, which rotate at synchronous
speed 𝜔𝑒 with respect to 𝑑 𝑠 and 𝑞 𝑠 axis and the angle 𝜃𝑒 = 𝜔𝑒 𝑡.
The process will remain same as above until the torque equation, Where the final torque equation
becomes [6-22],
3 𝑃
𝑇𝑒 = 2 ( 2 ) (𝜓𝑑𝑠 𝑖𝑞𝑠 − 𝜓𝑞𝑠 𝑖𝑑𝑠 ) (62)
3 𝑃 1
𝑇𝑒 = 2 ( 2 ) (𝜔 ) (𝐹𝑑𝑠 𝑖𝑞𝑠 − 𝐹𝑞𝑠 𝑖𝑑𝑠 ) (63)
𝑏
10
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
5. PARAMETERS
Lr
Tr = (65)
𝑅𝑟
Base speed,
wb = 2 ∗ π ∗ fb (66)
X ls = wb ∗ Lls (67)
X lr = wb ∗ Llr (68)
X m = wb ∗ Lm (69)
∗ 1
Xm = 1 1 1 (70)
+ +
Xlr Xls Xm
The above created models are tested under various conditions with machine parameters as shown in
table 1. Full load torque and speed of three phase induction motor thus simulated is shown in figure.13
and figure.14
The stator currents of three phase induction motor is shown in figure 15.
Similarly for five phase induction motor as shown in figure 16, 17, and 18.
11
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
12
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
The torque and speed of both the machines are compared using simulation data inspector, shown in
figures 19 and 20.
13
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
Figure 19. Torque comparison of three phase and five phase Induction Motor
Figure 20. Speed comparison of three phase and five phase Induction Motor
14
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
7. CONCLUSION
As it is evident from the simulation results that, under fully loaded condition five phase induction motor
took around 1 second to reach the steady state value. Whereas the three phase induction motor took around
1.5 seconds to reach the steady state value. Under constant changing loads, this parameter plays a very
important role in analyzing the system behaviour as well as it improves the reliability of the system. Also
under same machine parameters, the maximum torque of machine increased by 67 times approximately
with change in three phase to five phase. And we can say that per phase power distribution is improved.
The above model can be implemented for simulating multi-phase induction motors for higher number of
phases.
8. REFERENCES
[1] M. G. Simoes and P. Vieira, "A high-torque low-speed multiphase brushless machine-a perspective
application for electric vehicles," in IEEE Transactions on Industrial Electronics, vol. 49, no. 5, pp. 1154-
1164
[2] H. Guzmán, M. J. Durán, F. Barrero and S. Toral, "Fault-tolerant current predictive control of five-
phase induction motor drives with an open phase," IECON 2011 - 37th Annual Conference of the IEEE
Industrial Electronics Society, Melbourne, VIC, 2011, pp. 3680-3685.
[3] S. Williamson and S. Smith, “Pulsating torque and losses in multiphase induction machines,” IEEE
Trans. Ind. Appl., vol. 39, no. 4, pp. 986–993
[5] Lee, R.J., Pillay, P., and Harley, R.G., “D, Q Reference frame for the simulation of induction motors”,
Electric Power Systems Research Journal, vol. 8, pp. 15-16
[6] K. L. Shi, T. F. Chan and Y. K. Wong, “Modelling of the three-phase induction motor using
SIMULINK,” 1997 IEEE International Electric Machines and Drives Conference Record, Milwaukee,
WI, USA
[7] Chinthamalla, R., Karampuri, R., Jain, S., Sanjeevikumar, P., & Blaabjerg, F. (2018). Dual Solar
Photovoltaic Fed Three-Phase Open-End Winding Induction Motor Drive for Water Pumping System
Application. Electric Power Components and Systems, 46(16-17), 1896-1911.
[8] Kumar, B. J., & Chandramouli, A. (2017, November). Modeling and simulation of nine-level cascaded
H-bridge inverter based shunt active power filter for single-phase distribution system. In 2017
International Conference on Inventive Computing and Informatics (ICICI) (pp. 675-680). IEEE.
[9] Arulmurugan, R., & Chandramouli, A. (2019). Modeling of PV Powered Seven-Level Inverter for
Power Quality Improvement. In Smart Intelligent Computing and Applications (pp. 113-121). Springer,
Singapore.
[10] Karampuri, R., Jain, S., & Somasekhar, V. T. (2018). Common-Mode Current Elimination PWM
Strategy along with Current Ripple Reduction for Open-Winding Five-Phase Induction Motor Drive.
IEEE Transactions on Power Electronics, 34(7), 6659-6668.
15
ICRAEM 2020 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 981 (2020) 042059 doi:10.1088/1757-899X/981/4/042059
[11] Chinthamalla, R., Karampuri, R., Jain, S., Sanjeevikumar, P., & Blaabjerg, F. (2018). Dual Solar
Photovoltaic Fed Three-Phase Open-End Winding Induction Motor Drive for Water Pumping System
Application. Electric Power Components and Systems, 46(16-17), 1896-1911.
[12] Vedik B, Ritesh K, Deshmukh R and Shiva C K 2020 Renewable energy based load frequency
stabilization of interconnected power systems using quasi-oppositional dragonfly algorithm J Control
Autom Electr Syst. https://doi.org/10.1007/s40313-020-00643-3
[13] Vedik B, Naveen P and Shiva C K 2020 A novel disruption based symbiotic organisms search to
solve economic dispatch Evol. Intel. https://doi.org/10.1007/s12065-020-00506-5
[14] Kumar R, Sahu B, Shiva C K and Rajender B 2020 A control topology for frequency regulation
capability in a grid integrated PV system Archives of Electrical Engineering 69(2) 389–401
[15] Basetti V, Chandel A K and Subramanyam K B 2018 Power system static state estimation using
JADE-adaptive differential evolution technique Soft Computing 22(21) 7157-76
https://doi.org/10.1007/s00500-017-2715-3
[16] Shiva C K, Vedik B and Kumar R 2019 Integration of distributed power sources to hydro-hydro
power system subjected to load frequency stabilization International Journal of Engineering and
Advanced Technology 8(2) 128-32
[17] Pasha S N, Harshavardhan A, Ramesh D and Md S Shabana 2019 Variation analysis of artificial
intelligence machine learning and advantages of deep architectures International Journal of Advanced
Science and Technology 28(17) 488-495
[18] Swathi N, Padmaja Ch and Navya Jyothi G 2020 Audio assistive for blind people to identify the
cloth patterns and colors Journal of Critical Reviews 7(17) 154-158 10.31838/jcr.07.17.23
[19] Prakash T C, Mamatha M and Samala S 2020 An IoT based under weather monitoring system
Journal of Critical Reviews 7(17) 148-153 10.31838/jcr.07.17.22
[20] Kumar C N and Satyanarayana N 2015 Hybrid loss recovery technique for multipath load balancing
in MANETs 2nd International Conference on Electronics and Communication Systems ICECS 2015
1294-1301 10.1109/ECS.2015.7124793
[21] Sekhar V M Rao, KVG Rao N S and Chand M G 2016 Comparing the capacity NCC and fidelity of
various quantization intervals on DWT Advances in Intelligent Systems and Computing 413 45-55
10.1007/978-981-10-0419-3_6
[22] Mudi J, Shiva C K, Vedik B and Mukherjee V 2020 Frequency stabilization of solar thermal-
photovoltaic hybrid renewable power generation using energy storage devices Iran. J. Sci. technol. trans.
electr. eng. https://doi.org/10.1007/s40998-020-00374-w
16