Modeling and Simulation of A 3 Inducti
Modeling and Simulation of A 3 Inducti
Modeling and Simulation of A 3 Inducti
Corresponding Author:
Ahmed. J. Ali
Department of Power Technichal Engineering
Northern Technichal University, Iraq
Email: [email protected]
1. INTRODUCTION
Due to their robustness, reliability and low requirement for maintenance, three phase induction
motors have a many application in industrial processes. Different methods had been used to compute the self
and mutual inductances of electric machines. The winding function approach (WFA) had been extensively
used to calculate these parameters. The principle of this approach is depending on solving the differential
equations that explain the electromagnetic coupling between stator and rotor of the induction motor [1-3].
Calculation of coupled inductances are based on the basic laws of electrical circuit. Firstly, WFA had been
applied to compute the inductances of healthy induction motors, then the formulas of this method had been
modified to be suitable n-phase conventional induction motors. In the last decades of previous century it had
been applied for analysis of linear induction motor [2, 4, 5]. Classical methods that had been used to simulate
the electric machines have a limitation to add the effect of winding connection. Generally it assume that both
stator and rotor are consists of three phases only and ignore the effect of winding distribution. This drawback
had been removed by introducing the WFA where all types of winding can be introduced easily. Also the effects
of materials types and saturation can be taken into consideration [6-9]. Then the WFA had been modified in
order represents the actual distributions of magnetic flux which gives the correct simulation of electric
machine parameters. However, obtaining lower harmonics and simple manufacturing process is a challenge in
the design and optimization of analyst [10, 11]. Due to ability of WFA to simulate any type of winding
distribution or a variable air-gap width around rotor, it has been effectively used for modeling and simulation
of any type of abnormality of electrical machine operation. Fay model. The ability of WFA makes it more
suitable to simulate the inter turn sort circuit between two windings or with other phases in addition to simulate
the broken bars of electric machines. Finally it has been successfully for simulation the static, dynamic and mixed
rotor eccentricity [10, 12]. Multiple coupled circuit is assumed that the stator of the induction motor consists of
n-circuits while the rotor bars of m-circuits. These electrical circuits are mutually coupling and used to
determine the inductances of the motor using simple mathematical formulas. WFA give an appropriate way to
simulate all machine parameters, stator winding turns, rotor bar and air-gap eccentricities. According to
principle of this theory, the mutual and self-inductances between stator winding and rotor bars (windings)
stator-rotor circuits can be considered. The relation between these winding is time varying because it depends
of the rotor position with respect each other and with stator winding. The leakage inductance of both stator
winding and rotor bars are considered as constants and can be computed using special design formulas. In
current work two procedure of WFA have been considered. The slot effect is not included in the first type
where the coupling mmfs have a sinusoidal shape. While in second approach a real winding distribution in the
slots is assumed and mmfs have a non-sinusoidal waveform. A complete model of 3-Ф cage induction motor
which has an m-stator electrical circuits and n-rotor bars is implemented using matlab software [13]. An epic
strategy to display electrical machines dependent on WFA is proposed. In traditional WFA, it is accepted that
transition lines in air-gap are in outspread heading. The primary downside of this supposition that is
incorrectness of the reproduction results. In this technique WFA is improved by limited component strategy in
which the definite permeance work is determined by finite element method. This new methodology has a decent
exactness adjacent to of low estimation preparing time. In this approach the simulation results are presented and
show the performance of the proposed method and these results are compared to experimental results [14-16].
Int J Elec & Comp Eng, Vol. 11, No. 2, April 2021 : 1105 - 1113
Int J Elec & Comp Eng ISSN: 2088-8708 1107
where Figure 1 display the currents distribution of the rotor loops. Both loop shares two rotor bars with
the two national rings and two end ring sections with the two end-ring loops. Not the higher shape codes are
explained later and according to previous assumptions the complete model of the cage rotor induction motor
can be represented by:
where, Vs is the vector of the stator voltages, Is is the vector of the stator currents, Ir is the vector of the rotor
loops currents, Rs is the matrix of the stator windings resistances, Rr is the matrix of the rotor resistances, Ls
is the matrix of the stator windings inductances, Lr is the matrix of the rotor inductances, Lsr and Lrs are the
matrix of the mutual inductances between stator and rotor respectively φs and φr are the stator and rotor
windings linkage flux respectively [12, 18, 19].
From (1) and (3), the voltage equations of the stator per phase may be written as:
where Va, Vb and Vc are the per phase stator voltages ia, ib and ic are the per phase stator currents La, Lb and
Lc are the per phase stator self-inductances, Lab , Lbc and Lca are the per phase stator mutual inductances Lls is
the leakage inductance of the stator per phase. These equations are rewritten in suitable way to suit the basic
m.file. For example refer to phase "a" the solve one phase current ia (t) rewritten of equation in the following from:
𝑑𝑖𝑎 1 𝑑
= ((va− R s ∗ 𝑖𝑎) − (Lab 𝑖𝑏 + Lac 𝑖𝑐 + L𝑎𝑟 𝑖𝑟) ) (8)
𝑑𝑡 La 𝑑𝑡
where, Lab Lac , Lar are inductances and their value can be computed based on the position of rotor. Similarly
the above equations can be arranged for the other two phases (dib/dt and dic/dt) currents and two damper
windings currents. The complete m. file for three phase dynamic model of induction machine. The rotor
loops currents (ir1, ir2,…, irn) , (Lar1, Lar2,…, Larn) are mutual inductances of stator phase A to rotor loops,
(Lbr1, Lbr2,…, Lbrn) are mutual inductances of stator phase B to rotor loops, and (Lcr1, Lcr2,…, Lcrn) are
mutual inductances of stator phase C to rotor loops. All rotor circuits are mutually coupled with each other
and with stator phase’s windings. Based on 2 and 4 and with depend on Figure 1, the rotor loop voltage
equation k may be given as:
0 = 2(R b ) Irk − R b (Ir(k−1) + Ir(k+1) ) + 2(Lb + Le ) dIrk ⁄dt − Lb (dIr(k−1) ⁄dt + dIr(k+1) ⁄dt ) +
(Lk1 dIr1 ⁄dt + Lk2 dIr2 ⁄dt + ⋯ + Lkk dIrk ⁄dt + ⋯ + Lkn dIrn ⁄dt ) − d⁄dt (Lak Ia + Lbk Ib + Lck Ic (9)
where, Re, rb, Lb ,Le, are the resistance and inductance of each rotor bar. Lkk is the self inductance of each
rotor loop (Lk1, Lk2, …, Lk(k-1), Lk(k+1) ..., Lkn) are the mutual inductances loop k to the other rotor loops
(Lak, Lbk, and Lck) are the stator mutual inductances with rotor loop k.
𝑑𝐿𝑆𝑟
Te = [𝑖𝑠 ]𝑇 [ ] [𝑖r ] (10)
𝑑𝐿𝜃𝑟
where, θr is the angular position of rotor , Te is the developed torque produced by the motor.
Modeling and simulation of a 3-ф induction motor based on two types of WFA (Ahmed J. Ali)
1108 ISSN: 2088-8708
where,
J is moment of inertia of rotor,
wr is rotor angular speed,
B is the friction coefficient.
TL is the load torque,
Ѳr is the angular position of the rotor.
where: µ𝑜 is the permeability of air and =4π.10-7 H, 𝑔 is the width of the air gap, 𝜃 is the angular position of
the rotor, 𝑟 is the average radius of the air gap, 𝑙 is the effect length of the motor stack, 𝜗 is a certain point
along the air gap, 𝜆(𝜃, 𝜗) is called the winding function and represents the distribution of the magneto
motive force along the air gap when a current of '1' A flow in the winding. In current paper a three phase
induction motor, 36 stator slots and symmetrical concentrate winding. While the cage rotor of the motor
compromise a 28-equally spaced bars. The geometric construction of the motor has been taken into
consideration, where the slots effect on the stator winding is considered. The self any two stator phases Laa ,
Lbb , Lcc have been calculated based on WFM as given in (14). Where the motor inductances formulas can be
written as [25-27]:
π0 Lr αr
Lkk = ( ) (1 − ) αr (15)
g 2π
π0 Lr
Lk1 = Lk2 = ⋯ = Lkn = Lki ( ) (αr )2 /2 π )61(
g
where: Ns is the number of stator turns per phase, p is the number of pole pairs, αr is the angle between two
neighboring bars 𝐿𝑚 =(4/π𝑁𝑠 ) 𝐿𝑚𝑠 sin (pδ), where δ=𝛼𝑟 /2.
Int J Elec & Comp Eng, Vol. 11, No. 2, April 2021 : 1105 - 1113
Int J Elec & Comp Eng ISSN: 2088-8708 1109
4. SIMULATION RESULTS
The flux equations have been used to simulate the proposed model. The simulated steps that have
been used to solve the differential equations of the induction motor are explained in the algorithm shown in
Figure 2. The integration of flux linkage derivatives with respect to time are calculated to obtain flux linkages,
as follows:
dφ
= vsource − RI (20)
dt
Matlab m-file has been used to simulate the dynamic equations of a 3-Ф cage rotor induction motor.
The proposed mathematical model has been used to simulate the time response performance of a 4 kW
induction motor with the parameters given in Table 1.
According to the two suggested approaches, the model of motor has been simulated under operation
conditions by applying a load torque of TL=3.5 Nm on the shaft, starting time from zero, applying of
220 volt and 50 Hz. The time response of the stator current are shown in Figure 3 at Non-sinusoidal shapes
self and mutual inductance. Also the dynamic characteristics (Mutual inductance, rotor speed and developed
torque) are shown in Figures 4-6. The WFA principle is based on the coupling between stator and rotor as
shown Figure 4. Then for the second approach, the model has been simulated under the same operating
conditions mentioned for first approach. It's clear from the simulation results of stator and rotor currents are
shown in Figure 7 and Figure 8 that waveform have a smooth shape compare with Figure 3. This unifomlly and
reular current is due to sinusoidal assumption of mutual inductance that couple between the stator and rotor.
(a) (b)
Figure 3. The time response of the stator current, (a and b) stator current at non-sinusoidal shapes self and
mutual inductance using WFA
Int J Elec & Comp Eng, Vol. 11, No. 2, April 2021 : 1105 - 1113
Int J Elec & Comp Eng ISSN: 2088-8708 1111
(a) (b)
Figure 7. The simulation results of stator, (a and b) stator currents due to sinusoidal shapes self and mutual
inductance using WFA
(a) (b)
Figure 8. The simulation results of rotor, (a and b) rotor currents due to sinusoidal shapes self and mutual
inductance using WFA
5. CONCLUSION
Two types of winding function approaches have been modeled and simulated in this paper. It's
obvious from the results that the first model has results different slightly from the second approach. For the
stator, rotor currents and torque they have a non-smooth and ripple shapes compared with the results obtained
from the second model. This non-linearity is due to the actual winding distribution on the motor slots. The
stator MMF have a non-sinusoidal waveform due to slots effect which make this model more suitable to
simulate the faulty motor compared with first type.
REFERENCES
[1] J. Carrera, R. Franceschi, G. Gonzalez, “Induction Motor Fault Modeling Based on the Winding Function,” ESTEC
Conference Proceedings 6th Engineering, Science and Technology Conference, vol. 2018, 2017, pp. 758-767.
[2] N. L. Shmitz and D. W. Novotny, “Introductory Electromechanics,” NewYork: Roland, 1965.
[3] M. Ghods, Z. Nasiri-Gheidari, F. Tootoonchian and A. Oraee, "Performance Evaluation of Permanent Magnet
Vemier Generator Based on Modified Winding Function Method," Electrical Engineering (ICEE), Iranian
Conference on, Mashhad, 2018, pp. 1107-1112.
[4] Basil M. Saied, Ahmed. H. Ahmed, “Modeling and Simulation of Salient Pole Synchronous Machine with
Comparison between Actual and d-q Models,” Second Scientific Engineering Conference, Mosul, vol. 2, pp. 41-56, 2013.
[5] R. Krall, J. Krenn and A. Schmid, “Fractional slot winding versus distributed winding using winding function
method,” 2014 International Conference on Optimization of Electrical and Electronic Equipment (OPTIM), Bran,
2014, pp. 437-444.
[6] Ah-Jaco A., “Modélisation des moteurs asynchrones triphasés en régime Transitoire avec saturation et harmoniques
d’espace: Application au diagnostic,” thèse de doctorat de l’université de Lyon juillet, 1997.
[7] J. C. Moreria and T. A. Lipo, “Modeling of saturated AC machines including air gap flux harmonic components,”
in IEEE Transactions on Industry Applications, vol. 28, no. 2, pp. 343-349, 1992.
Modeling and simulation of a 3-ф induction motor based on two types of WFA (Ahmed J. Ali)
1112 ISSN: 2088-8708
[8] K. C. Deekshit Kompella, M. Venu Gopala Rao, R. Srinivasa Rao, “Baring Fault Diagnoses in 3 phase Induction
Machine Using Current Spectral Subtraction With Different Wavelet Transform Techniques,” Journal of Electrical
System, vol. 13, no. 1, pp. 143-159, 2017
[9] A. Doulatshah and P. Naderi, “Eccentricity fault detection of a salient-pole synchronous machine using modified
winding function approaches and finite element method,” 2016 7th Power Electronics and Drive Systems
Technologies Conference (PEDSTC), Tehran, 2016, pp. 29-34.
[10] H. A. Toliyat, M. M. Rahimian, and T. A Lipo, “Transient analysis of cage induction machines under internal faults
using winding function,” 3rd Int. Conf. Electrical Rotating Machines-ELROMA, 2006.
[11] R. Alipour-Sarabi, Z. Nasiri-Gheidari, F. Tootoonchian and H. Oraee, “Improved Winding Proposal for Wound
Rotor Resolver Using Genetic Algorithm and Winding Function Approach,” in IEEE Transactions on Industrial
Electronics, vol. 66, no. 2, pp. 1325-1334, 2019.
[12] H. A. Toliyat and T. A Lipo, “Transiant analysis of cage induction machines under stator, rotor bar and end ring
faults,” in IEEE Transactions on Energy Conversion, vol. 10, no. 2, pp. 241-247, 1995.
[13] S. M. Raziee, O. Misir and B. Ponick, “Winding Function Approach for Winding Analysis,” in IEEE Transactions
on Magnetics, vol. 53, no. 10, pp. 1-9, 2017.
[14] M. Nourmohammadpour, M. Khalilzadeh and K. Abbaszadeh, “Accurate modeling of switched reluctance motor
by using improved winding function method,” 2014 22nd Iranian Conference on Electrical Engineering (ICEE),
Tehran, 2014, pp. 793-797.
[15] F. Boubakar, A. Aissa, K. Abdellah and Z. Abdellim, “A three dimensional finite element method analysis and
winding function theory for torque ripple reducing in reluctance machine,” 2016 8th International Conference on
Modelling, Identification and Control (ICMIC), Algiers, 2016, pp. 418-423.
[16] F. Boubakar, K. Abdellah, A. Mehdi and A. Aissa, “Torque ripple reducing in SRM based on winding function
theory and 3D Finite Element Methods,” 2015 4th International Conference on Electrical Engineering (ICEE),
Boumerdes, 2015, pp. 1-5.
[17] S. Nandi, S. Ahmed, and H. A. Toliyat, “Detection of rotor slot and other eccentricity related harmonics in a three
phase induction moto with different rotor cages,” in IEEE Transactions on Energy Conversion, vol. 16, no. 3,
pp. 253-260, 2001.
[18] H. A. Toliyat and M. S. Arefeen, “A method for dynamic simulation of air-gap eccentricity in induction machines,”
in IEEE Transactions on Industry Applications, vol. 32, no. 4, pp. 910-918, 1996.
[19] G. M. Joksimovic, J. Penman, and N. Arthur, “Dynamic simulation of dynamic eccentricity in induction machines-
winding function approach,” in IEEE Transactions on Energy Conversion, vol. 15, no. 2, pp. 143-148, 2000.
[20] Petar Matić, Aleksandar Rakić, Slobodan N. Vukosavić, “Space Vector Representation of Induction Motor Model
in Field Weakening Regime,” Serbian Journal of Electrical Engineering, vol. 9, no. 1, pp. 53-61, 2012.
[21] Toliyat H. A., Lipo T. A., and White J. C., “Analysis of a concentrated winding induction machine for adjustable
speed drive application. I. Motor analysis,” IEEE Transaction on Energy Conversion, vol. 6, no. 4, pp. 679-684, 1991.
[22] Luo X., Liao Y., Toliyat H.A., El-Antably A. and Lipo T.A., “Multiple coupled circuit modeling of induction
machines,” IEEE Transaction on Industry apllications, vol. 31, no. 2, pp. 311-318, 1995.
[23] S. Shah, A. Rashid, M. K. L. Bhatti, "Direct Quadrate (D-Q) Modeling of 3-Phase Induction Motor Using
MatLab/Simulink,” Canadian Journal on Electrical and Electronics Engineering, vol. 3, no. 5, pp. 237-243, 2012.
[24] Ahmed. J. Ali, Ahmed. H. Ahmed and B. M. Saied, "Cogging torque Mitigation for PMSM using stator slots
design and Magnets skewing," 2019 2nd International Conference on Electrical, Communication, Computer,
Power and Control Engineering (ICECCPCE), Mosul, Iraq, 2019, pp. 240-245.
[25] D. Uma, K. Vijayarekha, “Modeling and Simulation of VSI Fed Induction Motor Drive in Matlab/Simulink,”
International Journal of Electrical and Computer Engineering (IJECE), vol. 7, no. 2, pp. 584-595, 2017.
[26] M. Zuo and H. Guo, "Integrative simulation and analysis of induction motor and control system," Proceedings of
2011 International Conference on Electronic & Mechanical Engineering and Information Technology, Harbin,
2011, pp. 2412-2414.
[27] K. E. Amin, R. M. Hamouda, A. D. Elkoshairy and M. A. El-Dessouki, “Dynamic Analysis of Electromechanical
Oscillations in Cascaded Induction Motors Driving Different Loads,” 2019 21st International Middle East Power
Systems Conference (MEPCON), Cairo, Egypt, 2019, pp. 473-477.
BIOGRAPHIES OF AUTHORS
Int J Elec & Comp Eng, Vol. 11, No. 2, April 2021 : 1105 - 1113
Int J Elec & Comp Eng ISSN: 2088-8708 1113
Modeling and simulation of a 3-ф induction motor based on two types of WFA (Ahmed J. Ali)