Design and Simulation of A Three-Phase Induction M
Design and Simulation of A Three-Phase Induction M
Design and Simulation of A Three-Phase Induction M
[email protected]|[email protected]|[email protected]
Abstract- Three-phase Induction motors are increasingly used in a variety of applications such as fans, milling machines, transportation, etc.
In applications requiring precise speed control such as robotics, centrifugal pumps, mills and other high-performance applications, it is
essential that speed is maintained constant at a desired fixed value. Many speed control techniques exist in literature but this paper presents
the design of a simple PI controller and its use to control the speed of a three-phase induction motor. A controller design that applies the
Phase margin (PM) as the stability criterion is employed. The cross-over frequency and phase margin of the open loop gain of motor and
controller are specified in order to develop a stable control system. The three-phase induction motor model is presented and modified for the
purpose of control. A comparative analysis between the motor performance with and without the PI controller was performed. The unit step
response of the speed control loop is characterized by rise time, settling time, steady state error and peak overshoot of 0.0253s, 0.19s, 2.2e-
14% and 24.03% respectively. Simulation results show that the speed of the uncontrolled motor changed whereas that of the controlled motor
returned quickly to its initial value after the motor is subjected to load disturbances in steady state.
2 MATHEMATICAL MODEL ′
𝜆𝑑𝑠 = 𝐿𝑠 𝑖𝑑𝑠 + 𝐿𝑚 𝑖𝑑𝑟 (5)
In modelling the three-phase induction motor, the three
phase windings of the stator and rotor are usually
′
represented by equivalent orthogonal set of d-q axis 𝜆𝑞𝑠 = 𝐿𝑠 𝑖𝑞𝑠 + 𝐿𝑚 𝑖𝑞𝑟 (6)
windings which produce the same MMF in the airgap.
The schematic diagram of a line-operated three-phase 𝜆′𝑑𝑟 = 𝐿′𝑟 𝑖𝑑𝑟
′
+ 𝐿𝑚 𝑖𝑑𝑠 (7)
induction motor is illustrated in Figure 1. Its equivalent in
d – q coordinates is drawn in Figure 2. The a-b-c phase 𝜆′𝑞𝑟 = 𝐿′𝑟 𝑖𝑞𝑟
′
+ 𝐿𝑚 𝑖𝑞𝑠 (8)
winding voltages 𝑉𝑎𝑠 , 𝑉𝑏𝑠 and 𝑉𝑐𝑠 are first transformed to
their equivalent d-q winding voltages, 𝑉𝑑𝑠 and 𝑉𝑞𝑠 in order 3 𝑃
to analyze motor operation. The d-q reference frame 𝑇𝑒 = ( ) (𝜆′𝑑𝑟 𝑖𝑞𝑠 − 𝜆′𝑞𝑟 𝑖𝑑𝑠 ) (9)
2 2
model of the three-phase induction motor is provided in
equations (1)-(10) (Krause et al., 2013; Metwaly et al., 2019). 𝑝𝑤𝑚 =
𝑇𝑒 −𝑇𝐿
(10)
They are the voltage, flux linkage and torque equations of 𝐽𝑒𝑞
′ ′
quantities. 𝑖𝑑𝑠 , 𝑖𝑞𝑠 , 𝑖𝑑𝑟 and 𝑖𝑞𝑟 are the stator and rotor modified because of the alignment of the rotor flux
currents. 𝜆𝑞𝑠 , 𝜆𝑑𝑠 , 𝜆′𝑑𝑟 and 𝜆′𝑞𝑟 are the flux linkages in the magnitude with the d-axis to the following;
stator and rotor. Rs and Rr are the stator and rotor winding
resistances respectively. 𝐿𝑚 is the magnetizing 3𝑃 𝐿𝑚
𝑇𝑒 = (𝜆′𝑟 𝑖𝑞𝑠 ) (15)
2 𝐿′𝑟
inductance. 𝐿𝑠 is self-inductance of the stator windings
while 𝐿′𝑟 is the self-inductance of the rotor windings. It is important to note that the rotor voltages 𝑉𝑑𝑟 and 𝑉𝑞𝑟
Primed quantities are quantities of the rotor referred to of the squirrel cage motor are each equal to zero because
the stator. 𝜔𝑒 and 𝜔𝑠𝑙 are the supply and slip angular the rotor windings are short-circuited. This point has been
frequency respectively. 𝐽𝑖 is the inertia of the rotor and the taken into consideration in the modeling. The dynamics
connected load while 𝑇𝐿 is the load torque. 𝑇𝑒𝑚 is the of the rotor flux is developed from (3) and (7) to yield;
electromagnetic torque produced in the motor.
𝑑𝜆′𝑟
𝜏𝑟 + 𝜆′𝑟 = 𝐿𝑚 𝑖𝑑𝑠 (16)
𝑑𝑡
3 MODIFICATIONS TO THE MODEL
𝐿′𝑟
In order to effectively decouple the flux from torque in the where 𝜏𝑟 =
𝑅𝑟
motor, the d-axis is now made to align with the rotor flux
linkage space vector as in Figure 3. The Figure illustrates The term, 𝜏𝑟 , is the rotor circuit time constant. The
the relationship between the synchronously rotating d-q dependence of the rotor flux on 𝑖𝑑𝑠 is evident in Equation
reference frame and the stationary reference frame (in (16). The equation describes a first order system with a
which the real machine currents are represented) when time constant given by 𝜏𝑟 . If the rotor flux, 𝜆′𝑟 , is
the rotor flux vector is properly aligned with the d-axis. maintained constant, then according to equation (15) the
The d-axis component of this current is seen to be in phase torque, 𝑇𝑒 will depend entirely on 𝑖𝑞𝑠 , which is the
with the rotor flux. component of the stator current in the q-axis. In this
manner, the control of torque and flux may be achieved
independently. The modified model which achieves the
decoupling of the control of torque and flux in the three-
phase induction motor is therefore stated as follows;
′
𝜆𝑑𝑠 = 𝐿𝑠 𝑖𝑑𝑠 + 𝐿𝑚 𝑖𝑑𝑟 , (22)
′
𝜆𝑞𝑠 = 𝐿𝑠 𝑖𝑞𝑠 + 𝐿𝑚 𝑖𝑞𝑟 , (23)
Fig. 3: d-axis alignment with the rotor flux vector ′
𝜆′𝑟 = 𝐿′𝑟 𝑖𝑑𝑟 + 𝐿𝑚 𝑖𝑑𝑠 , (24)
This means that the q-axis component of the rotor flux 0= 𝐿′𝑟 𝑖𝑞𝑟
′
+ 𝐿𝑚 𝑖𝑞𝑠 , (25)
linkage space vector, 𝜆′𝑞𝑟 is zero (Umar, Akbar and Kazmi,
2018). The d-axis component of the rotor flux space vector, The rotor flux angle is calculated as follows;
𝜆′𝑑𝑟 is now equal to the rotor flux magnitude, 𝜆𝑟 . The rotor 𝑡
𝜃𝑒 = ∫0 𝜔𝑒 (𝜏)𝑑𝜏 (26)
flux vector is computed as follows;
Where 𝜏 is the variable of integration. 𝜃𝑒 is extremely
crucial to the realization of the control scheme as its value
𝜆′𝑟 = √𝜆′2
𝑞𝑟 + 𝜆′2
𝑑𝑟 (11)
at every time step is required in order to calculate the
command currents that will drive the motor at the desired
Mathematically, aligning the d-axis with the rotor flux
speed. 𝜃𝑒 is also required for the transformation of the a-
vector means equating 𝜆′𝑞𝑟 to zero in equations (3), (4), (8)
b-c currents to their equivalent d-q values at every time
and (9). Therefore, the following modifications result;
step in accordance to Equations (27) and (28). This action
maintains the alignment of the d-axis with the rotor flux
𝜆′𝑟 = 𝜆′𝑑𝑟 (12) vector. (Abad, 2016).
1 1
′
𝑖𝑞𝑟 =−
𝐿𝑚
𝑖 (13) 1 − 𝑖𝑎
𝐿′𝑟 𝑞𝑠
𝑖𝛼 2 2 2
[𝑖 ] = [ √3 √3
] [ 𝑖𝑏 ] (27)
𝛽 3
𝑅𝑟 𝐿𝑚 𝑖𝑞𝑠 0 − 𝑖𝑐
𝜔𝑠 − 𝜔𝑟 = (14) 2 2
𝐿′𝑟 𝜆′𝑟
𝑖𝑑𝑠 cos 𝜃𝑒 sin 𝜃𝑒 𝑖𝛼
′
With 𝑖𝑞𝑟 substituted as provided in (13), the torque [𝑖 ] = [ ][ ] (28)
𝑞𝑠 − sin 𝜃𝑒 cos 𝜃𝑒 𝑖𝛽
produced in the motor represented by Equation (9) is
The d-q quantities are constants (or DC values) in The behaviour of the motor speed when the load is
balanced sinusoidal steady state and therefore it is easy to suddenly reduced and increased within its rated value is
design a PI controller for motor speed control. examined and discussed in section IV.
converter is employed to deliver desired current to the using the inverse transformations of Equations (28) and
three phase-induction motor. The objective of the control (27). The reference currents are translated to the motor as
design is to make the motor speed constant at rated value command inputs by the converter.
regardless of any load disturbances. A simplified
equivalent block diagram of the speed control loop is
provided in Figure 6.
∗
𝝎𝒎 𝒊∗𝒒𝒔 𝒊∗𝒂𝒃𝒄 Converter 𝒊𝒂𝒃𝒄 Motor 𝝎𝒎
+_ PI 𝐾𝑠−1 model
(Figure 4)
𝒊∗𝒒𝒔 1
3𝑃 𝜆′𝑟 𝐿𝑚 𝑇𝑒𝑚 𝝎𝒎
∗
𝝎𝒎
+_ PI 𝑠𝐽𝑒𝑞
2 𝐿′𝑟
𝑒𝑞𝑛 (10)
where,
1 1 PM = 1800+ Phase at crossover frequency (35)
1 −
𝑐𝑜𝑠 𝜃
2 𝑠𝑖𝑛 𝜃 2 2
𝐾𝑠 = [ ][ ] (29)
3 − 𝑠𝑖𝑛 𝜃 𝑐𝑜𝑠 𝜃 0 √3 √3
− Equations (36) and (37) satisfy this condition for the speed
2 2
control loop.
3.1 CONTROLLER DESIGN 𝑘
𝑘𝑝 (𝑗𝜔𝑐 + 𝑖 )𝑐
The control system block diagram of Figure 6 is reduced 𝑘𝑝
| |=1 (36)
as shown in Figure 7. −𝜔𝑐2 𝐽𝑖
𝑘
∗
𝜔𝑚 ∗ 𝑘𝑝(𝑗𝜔𝑐 + 𝑖 )𝑐
𝐼𝑞𝑠 𝜔𝑚 𝑃𝑀 = ∠ (
𝑘𝑝
) + 1800 (37)
𝐺𝐶 (𝑠) 𝐺𝑃 (𝑠) −𝜔𝑐2 𝐽𝑖
+_
Bode Diagram
Gm = -Inf dB (at 0 rad/s) , Pm = 60 deg (at 50 rad/s)
4 RESULTS AND DISCUSSION
80
The induction motor was operated in steady state. The
60 load on the motor was suddenly reduced by half from
rated value of 1 pu to 0.5 pu and stayed at that value for
Magnitude (dB)
40
-135
performance is plotted in Figure 10 (b).
1.5
-180
Torque (Nm)
10
0
10
1
10
2 3
10 1
Frequency (rad/s)
-5 1820
-10
1800
Speed (rev/min)
-15
1780
-20
-25 1760
-30 1740
0
1720
0 1 2 3 4 5 6
(b) Time (s)
Phase (deg)
-45 Fig. 10 (a): Load torque and Induced torque (b): Speed
performance of the uncontrolled three-phase induction motor
System: G1
Phase Margin (deg): 120
Delay Margin (sec): 0.0419
-90
At frequency (rad/s): 50 Transients can be observed in the induced torque at the
Closed loop2 stable? Yes
10
0
10
1
10 10
3
instant of load change in Figure 10 (a). The induced torque
Frequency (rad/s)
is equal to the load torque after the transients die out at
Fig. 9: Closed loop response of the speed control loop steady state. The speed of the motor, on the other hand,
rises from its rated value of 1770 rev/min to 1785 rev/min
after the change in load from 1 pu to 0.5 pu and then from
The parameters of the selected three-phase induction 1785 rev/min to 1793 rev/min after the reduction in load
motor are provided in Table 1. from 0.5 pu to 0.25 pu. The speed then falls to 1785 rev/min
with a rise in load torque at 3.5 s of simulation time. The
speed is seen to fall back to the rated value of 1770
Table 1. Parameters of The Three-Phase Induction Motor
rev/min after the rise in load torque back to its rated
value. The speed is observed to undergo some transients
Parameters Values at the instant of the changes. Oscillations can be seen in
Rated Power 2.4kW the speed as well during these changes.
Line Voltage 460V
In an un-controlled system such as this one, speed
Line frequency 60Hz
changes after load disturbances lead to inefficient motor
Rated Speed 1770 rpm operation. In high performance applications such as
Stator resistance 1.77 𝞨 robotics and factory automation, speed changes can lead
Rotor resistance 1.34 𝞨 to loss of precision. Neither of these conditions is
desirable. Hence the need for speed control. The unit step
Stator leakage inductance 0.01392 H
response of the controlled system is plotted in Figure 11.
Rotor inductance 0.0126 H It can be observed that the step response tracks its
Mutual inductance 0.369 H reference after a single peak overshoot and a slight
undershoot. Associated system performance indices were
calculated using the control system tool box in MATLAB
and are provided in Table 2.
1.4
(i.e., rated) value of 1770 rev/mins. It maintains this value
till end of simulation time despite the changes in load
1.2 torque. Oscillations in speed have reduced significantly
as well and this reduction can be observed in Figure 12
1 (b). There is an overshoot and a slight undershoot in the
value of motor speed following the disturbances. These
0.8
results validate the prediction from the step response plot
Amplitude
of Figure 11.
0.6
5 CONCLUSION
0.4
The speed and torque responses of an un-controlled and
controlled three phase induction motor were examined
and discussed. The d-q reference frame model of the three-
0.2
phase induction motor was presented and modified to
allow for an independent control of torque and flux in the
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 motor. A PI controller was designed in the light of
Time classical control theory and analytical expressions were
Fig. 11. Unit step response of the speed control loop obtained for controller parameters 𝑘𝑖 and 𝑘𝑝 . In carrying
out the controller design, the chosen cross-over frequency
Table 2. Command-Tracking Performance of the Speed was imposed at a desired phase margin that produces a
Control Loop stable system. The crossover frequency was selected as 50
Performance Index Value rad/s so as to obtain a very fast speed response and a
Rise Time 0.0253s phase margin of 600 was chosen for a stable closed loop
Maximum Overshoot 24.033% operation. The parameters 𝑘𝑖 and 𝑘𝑝 of the PI controller
Peak Time 0.04s were determined analytically from two closed form
Settling Time 0.19s
equations derived from the imposition of the chosen
Steady State Error 2.22e-14 %
cross-over frequency at the desired phase margin that
makes for a stable closed loop. A step response of the
The controlled system was also implemented in controlled motor was obtained and characterized.
MATLAB/Simulink with controller parameters 𝑘𝑖 and 𝑘𝑝
calculated as 14.12 and 0.489 respectively. The speed A significant change in speed of the uncontrolled motor
performance of this controlled Induction motor is plotted was observed when the motor was subjected to load
in Figure 12. disturbances in contrast to that of the controlled motor
which returned quickly to its initial value after similar
1.5
Load torque disturbances, thereby validating the prediction from the
Induced torque (Uncontrolled)
step response. The results obtained demonstrate the
Torque (Nm)
1780
REFERENCES
1760
Abad, G. (2017). Power Electronics and Electric Drives for Traction
Uncontrolled speed
1740 Reference speed Applications. John Wiley & Sons, United Kingdom.
Controlled speed Alahmad, A. and Kaçar, F. (2021). Simulation of Induction Motor
1720
0 1 2 3 4 5 6
(b) Time (s) Driving by Bridge Inverter at 120°, 150°, and 180° Operation. 8th
International Conference on Electrical and Electronics Engineering
Fig. 12 (a): Torque performances (b): Speed performances of the (ICEEE), pp. 121-125. Doi: 10.1109/ICEEE52452.2021.9415930.
uncontrolled and controlled three-phase induction motor Biswal, P. and Satpathy, S. (2021). Vector Control of 3-Phase
Induction Motor. 1st Odisha International Conference on Electrical
Given same load disturbances at the same simulation Power Engineering, Communication and Computing Technology
times as before, it can be observed in Figure 12 (a) that, (ODICON). pp. 1-4. Doi: 10.1109/ODICON50556.2021.9428930
although the induced torque response of the controlled El-Zohri, E. H. and Mosbah, M. A. (2019). Comparative Evaluation
motor undergoes some transient, the oscillations are of PWM Techniques Used at Mega 328/p with PI Control for
significantly reduced. The speed of the controlled system Inverter-Fed Induction Motor. Proceedings of the International
in Figure 12(b), on the other hand, rises or falls sharply Conference on Innovative Trends in Computer Engineering (ITCE),
depending on whether the load torque is reduced or Aswan, Egypt, pp. 219-224.
increased at the instant of the disturbances but tracks the Hareesh, A. and Jayanand, B. 2021. Scalar and Vector Controlled
commanded or reference speed and returns to its initial Infinite Level Inverter (ILI) Topology Fed Open-Ended Three-
Phase Induction Motor. IEEE Access. 9: 98433-98459.