Mathematical Modeling and Analysis of Different Vector Controlled CSI Fed 3-Phase Induction Motor Drive

The main objective of this paper is to build a simple mathematical competent model that describes the circuits and interconnections of a 3-phase squirrel cage induction motor used for industrial applications. This paper presents the detailed analysis of theoretical concepts used in mathematical modeling, simulation’ and hardware implementation. The objective of this work is to compare the dynamic performances of the vector control methods for CSI fed IM drives. Based on the results, dynamic performances of the proposed drives are individually analysed using the sensitivity tests. The tests that are chosen for the comparison are step changes in the reference speed and torque of the motor drive. Here the IM is mathematically modeled in different reference frames for input output linearization (IOL) control, field oriented control (FOC), and direct torque control method (DTC) which are designed using hardware equivalent mathematical equations. The most important contributions in this paper are mathematical simulation structure of IM model in rotor flux frame using current and speed that were developed and implemented in MATLAB-Simulink. The operation and performance of the different vector control methods are verified by simulation using MATLAB/SIMULINK and experimental results.


Introduction
The alternating current (AC) motor, especially asynchronous three-phase induction motor (IM), has been the motor of choice in industrial settings for about the past half century as power electronics can be used to control its output behavior.Before that, the direct current (DC) motor, especially the separately excited one was widely used because of its easy speed and torque controllability.The IM is a rugged structure motor because it is brushless and has very fewer internal parts that need maintenance or replacement.This makes it cheaper in comparison to other motors, such as the DC motor.Thus, IM and its drive system have been gaining market share in industry and even in alternative applications such as hybrid electric vehicles.For variable speed electric motor applications in low to moderate power, pulse width modulation (PWM) with voltage source inverter (VSI) is usually used.However, the switched voltages produce high voltage slopes over the stator windings, which stress the insulations and cause bearing current problems.A possible solution is to use PWM with current source inverter (CSI).Phillips [1], Lipo and Cornell [2], Palaniappan et al. [3], Kaimoto et al. [4], Krishnan et al. [5], Kazmierkowski and Koepcke [6], and Hombu et al. [7] tried to develop the CSI and make it suitable for ac motor (especially IM) variable speed drives; since both stator current and voltage waveforms are close to the sinusoidal waveform, the above problems are reduced.Further, CSI can be used for high power applications.In the recent years, the high performance CSI fed variable speed drives, especially IM drives have been dominated by various control methods.The important methods are vector based / control, input output linearization (IOL) control, indirect field oriented control (IFOC), and direct torque control (DTC) as each method has its own unique features./ is by far the simplest drive since it requires no parameter knowledge and is essentially an open loop drive.IOL is one type of nonlinear state feedback control which is completely input-output decoupled at all times even in transients.IFOC offers similar dynamic torque control performance to that of DC motors, giving fast, near step changes in machine torque 2 Journal of Applied Mathematics and DTC provides simpler control architecture with a similar dynamic performance as that of IFOC.
Similar studies of vector controlled inverter fed drives are discussed mostly with the voltage source inverters and very few papers are available with the study of vector controlled current source inverters; Veerachary [8], Salo and Tuusa [9], and Babaei and Heydari [10].This paper presents a comparative analysis of different mathematical modeling of vector controlled CSI fed induction motor drives and in particularly the dynamic responses of the drive motor speed against change in speed as well as the torque references in a single simulation cycle and verified experimentally and to ascertain different mathematical dynamic IM modeling for CSI fed drives as Holmes et al. [11] and to apply control methods on them.The MATLAB s-functions are described to design IM modeling and drives.Finally, the above mentioned drive strategies are compared as fairly as possible by comparing their dynamic responses when there is step change in the references (usually stator currents, motor speed, and load torque).

Mathematical Modeling of 3-Phase Induction Motor Drive
In this section, Mathematical model of IM is shown in Figure 1 and described.In general, IM is used to transform electrical energy to mechanical energy.It consists of electric circuitry, electromagnetic circuitry, and electromechanic circuitry Krishnan [12].The main objective of the motor modeling is to build a simple but competent model that describes these circuits and their interconnections.The main path in the model is in relation with the voltage or current and the stator phases in input and magnetic flux inside the motor and electromagnetic torque in the output.As mentioned earlier, the squirrel cage IM is prepared for industry applications.The IM can either be supplied by a voltage source or a current source.These sources are assumed to be ideal; this means that they can supply any desired voltage or current without losses.According to Krause and Thomas [13], the basic equations of the IM are the stator voltage equation (1), the rotor voltage equation ( 2), the equations of the flux linkages of the stator, and the rotor winding systems, (3) and ( 4), and the torque equation (5).Consider the following: The equations are given in vector format Bose [14].Every vector variable has two components: the first component is parallel to the reference frame axis or "" axis and the second component is perpendicular to the reference axis or "" axis.
The symbols are V-voltage, -current,   -flux linkage, -resistance, and   -leakage inductance; subscript "" indicates the stator winding system, "", the rotor winding system; superscript "" indicates the stator reference frame (stator coordinates), "", the rotor reference frame, and "", any arbitrary reference frame   is the electromagnetic torque; (/2) is a rotation matrix over /2 radians; the dot ( ̇) indicates a time derivative.The stator voltage equation is restricted to the stator reference frame because of the differential of the stator flux vector, whereas the rotor voltage equation is restricted to the rotor reference frame because of the differentiation.The flux equations and the torque equation can be defined in any reference frame, denoted by superscript "." It is evident from (1) to ( 5) that the motor voltage, flux, and torque can be controlled easily by controlling the current of the motor.So, the CSI is the right choice to control the motor directly because the inputs of the current fed IM model are the stator currents.
Before these equations are completed to obtain the IM model, the reference frame of the model in relation to the problem being investigated and the type of computer (analog or digital) are chosen.In order to evaluate the usefulness of each reference frame, the above equations are simulated using Matlab software to predict the transient performance of an IM.
The mathematical modeling of 3-phase induction motor drive which can be run at any of the three rotating reference frame models, namely, stator rotating reference frame, rotor rotating reference frame, and synchronous rotating reference frame in Matlab/Simulink is shown below in Figure 2.

Transient Performance of Different Mathematical Reference Frame Models
Figure 3 shows the transient performance of "" phase and -axis stator currents in each reference frame when there is a step change in the reference torque from 0 to full load at 1.4 seconds.From Figure 3(a), it is evident that in stator or stationary reference frame stator -axis current behaves exactly the same way as do the stator "" phase current of the motor.If the stationary reference frame is used, then the stator -axis current is identical to the stator phase current.This would be useful when interest is specifically confined to stator variables, for example, variable speed stator-fed IM drives.Similarly, if the rotor reference frame is used then the rotor -axis current will behave in exactly the same manner as the rotor "" phase current.This would be useful when interest is confined to rotor variables only, as, for example, variable speed rotor-fed IM drives.From Figure 3(b), it is evident that in rotor rotating reference, frame stator -axis current does not behave as the stator phase current.From Figure 3(c), it is apparent that in synchronously rotating reference frame, stator -axis current behaves like DC quantity.If the synchronously rotating reference frame is used, then the stator and rotor current behave like DC quantity.Thus, the controller design is simple and is widely used in AC drives.Figure 4 shows the transient performance of -axis stator current and -rotor flux in each reference frame when there is a step change in the reference torque from 0 to full load at 1.4 seconds.From Figures 4(a) and 4(b), it is apparent that there is a coupling effect between stator -axis current and rotor -axis flux.So, independent control is difficult.Figure 4(c) shows that the rotor -axis flux and stator axis current (control variables) have linear relationship with torque (speed).Thus, it is useful in the rotor flux oriented (RFO) control methods.
Figure 5 shows the transient performance of -axis stator current and -stator flux in each reference frame when there is a step change in the reference torque from 0 to full load at 1.4 seconds.Figures 5(a  The basic IM equations from Krause et al. [15] are simulated for three different reference frames using MATLAB 7.6 software and the results obtained are compared with each reference frame.From the comparison, it is observed that if the stator rotating reference frame is used, then stator -axis current is identical to that of the stator "" phase current as shown in Figure 3(a), which is used in stator fed IM drives.Similarly, in rotor reference frame, rotor -axis current is identical to the rotor "" phase current and useful in rotor fed IM drives.If the synchronous reference frame is used, then the stator -axis current behaves like DC quantity as shown in Figure 3(c).From Figures 4(c) and 5(c), it is evident that synchronous reference frame is useful in RFO and SFO drive control methods.In this paper, synchronous reference frame is preferred to model the IM in /, IOL control methods.RFO control is used in IFOC method and SFO is used in DTC method.

IOL Control for CSI Fed IM Drives
In this section, a type of IOL control method for the CSI fed IM drive is described.The controller is designed based on pole placement technique of Wonham [16] to a - axis state space linearized model of the drive which includes a reduced order rotor current observer.A control law has been developed to improve the dynamic response of the drive by using the concepts from Veerachary [8] which includes feed forward as well as feedback controls.

Mathematical Modeling of IOL Control Method.
The equivalent circuit for the CSI, capacitor, and IM in a - axis rotating reference frame is shown in Figure 6.The corresponding state equation (differential equations) with motor flux linkages as variables is derived from Wu et al. [17].
Veerachary [8] presented the linearized CSI-IM drive state space model in synchronously rotating reference frame as follows: where where, With usual state feedback control, there can be the tracking error in steady state.To eliminate this error, Smith and Davidson [18] introduced an integral control which includes feedback as well as feed forward control.In this work, an optimal control law is derived by using the concepts of Smith and Davidson [18] and Veerachary [8].

Design of Optimal Control Law.
A new multivariable optimal controller incorporating state feedback as well as feed forward control for fast regulation and stability of a CSI fed IM drive is presented in this section.The design of the state feedback controller is based on the industrial regulator theory together with pole placement technique applied on a linearized - axes state space model of the drive and includes a reduced order observer to estimate the inaccessible states like - axes rotor currents.A feed forward control in terms of reference and disturbance inputs has been added to the feedback controller-observer to obtain faster dynamic response either from (10) or (11) as follows (Figure 7): or

IFOC for CSI Fed IM Drives
In this section, a high performance current fed indirect RFO control method is described for IMs.The IM is mathematically modeled in the rotor flux reference frame that is presented.The rotor flux orientation is also obtained.In the low speed region, the rotor flux components can be synthesized more easily with the help of speed and current signals.Design of indirect vector controller Blaschke [19,20] and rotor flux estimator is discussed.Implementation of IFOC to CSI fed IM drive and transient response analysis of other simulation results are also discussed.

Mathematical Modeling of IFO Control Method.
The rotor equations of the IM containing flux linkages as variables are given by Krause et al. [15] as follows: where then (12) becomes The resultant rotor flux linkage,   , also known as the rotor flux linkages phasor, is assumed to be on the direct axis to reduce the number of variables in the equations by one.Moreover, it relates to reality that rotor flux linkages are a single variable.Hence, aligning d axis with rotor flux phasor yields Substituting ( 15) into ( 14) results in the new rotor equations as follows: The rotor currents in terms of the stator currents are derived from (18) as The and -axes currents are relabeled as torque (  ) and flux producing (  ) components of the stator-current phasor, respectively.  denotes rotor time constant.Similarly, by the same substitution of the rotor currents from ( 19) into torque expression, the electromagnetic torque is derived as where the torque constant   is defined as Note that the torque is proportional to the product of rotor flux linkage and the stator -axis current.This resembles the air gap torque expression of the DC motor which is proportional to the product of the field flux linkages and the armature current.If the rotor flux linkage is maintained as constant, the torque is proportional to the torque producing component of the stator current.This is in relation to the separately excited DC motor with armature current control, where the torque is proportional to the armature current when the field current is constant.Further, here too, the time constant is measured in the order of a few milliseconds.The rotor flux linkages and air gap torque given in ( 21) and ( 22), respectively, complete the transformation of the IM into an equivalent separately-excited dc motor from a control point of view.

Design of IFO Controller.
The IFO controller was designed using the concepts of Vas [21] and Salo and Tuusa [9].
From that the stator-current phasor is the phasor sum of the and -axis stator currents in any frames and it is given as and the  axes (2) to (3) abc phase current relationship is obtained from And it can be written as where where   ,   , and   are the three phase stator currents.Note that the elements in the  matrix are cosinusoidal functions of electrical angle,   .The electrical field angle in this case is that of the rotor flux-linkages phasor and is obtained as the sum of the rotor and slip angles as follows: and the slip angle is obtained by integrating the slip speed and is given as IFO controller, developed from these derivations, accepts the torque and flux requests and generates the torque and flux producing components of the stator-current phasor and the slip-angle commands.The command values are denoted with asterisks throughout this thesis.From ( 20), (21), and (23), the command values of   ,   , and  sl are obtained as follows: The command slip angle,  * sl , is generated by integrating  * sl .The torque angle command is obtained as the arctangent of  *  and  *  .The field angle is obtained by summing the command slip angle and rotor angle.With the torque and flux producing components of the stator current commands and rotor field angle, the  axes current commands ( phase current commands) are obtained as follows.The relevant steps involved in the realization of the IFO controller are as follows: where By using ( 35) to (37), the stator and -axes and abc current commands are derived as where The implementation of IFOC on a CSI fed IM is shown in Figures 8 and 9.The torque command  *  is generated as a function of the speed error signal generally processed through a PI controller.The flux command  *  can be given directly or as a function of speed.The rotor position   can be measured with an encoder and converted into necessary digital information for feedback.There has been a substantial amount of research in the development of rotor flux observers for field orientation that are compensated for variations in parameters by their feedback corrections.Digital implementation of integrators for the estimation of rotor flux of an IM from the stator voltages and stator currents poses problems associated with the offset in the sensor amplifiers.Traditional low-pass filters can replace the integrator.So, in this work, rotor flux estimated using low pass filter is described.The instantaneous flux linkage can be computed using the measured d-axis stator current using (40), which is referred to as the current model as follows: The theoretical concepts of IFOC method for CSI fed IM drive discussed in Sections 5.1 and 5.2 are verified through simulation using the Simulink toolbox of MATLAB.The detailed simulations results are presented in Section 7.1.

DTC for CSI Fed IM Drives
Variable speed drives are used in all industries to control precisely the speed of electric motors driving loads ranging from pumps and fans to complex drives on paper machines, rolling mill cranes, and similar drives.The most modern method for these drives is direct torque and stator flux vector control method usually called DTC.The VSI fed version has been realized in an industrial way by ABB, by using the theoretical background proposed by Takahashi and Ohmori [22] and Depenbrock [23].This solution is based both on field oriented control (FOC) as well as on the direct selfcontrol theory.The idea is that motor flux and torque are used as primary control variables which is contrary to the way in which traditional AC drives control input frequency and voltage, but is in principle similar to what is done with a DC drive, where it is much more straightforward to achieve.
Casadei et al. [24] has described that by controlling motor torque directly, DTC provides dynamic speed accuracy equivalent to closed loop AC and DC systems and torque response times that are 10 times faster.It is also claimed that the DTC does not generate noise like that produced by conventional PWM AC drives.This work describes a new modular approach to implement DTC method on a CSI fed IM drive, the simulation implementation procedure.The dynamic response of the drive is analyzed and presented.

Mathematical Modeling of DTC Control Method.
The DTC method is developed from Babaei and Heydari [10].From Krause et al. [15], in the stationary reference frame, the stator voltage and flux equation can be written as where   is the stator resistance and   is the stator flux vector in the stationary reference frame.|  | and ∠  are the amplitude and position of stator flux vector, respectively.In DTC of CSI fed IM, it is desirable to determine the stator current reference so that the torque and stator flux follow their reference values, Babaei and Heydari [10].The proposed DTC system is based on the stator flux oriented (SFO) reference frame.In this rotating reference frame, it is written as where  and  are the real and imaginary axes in the SFO reference frame.Therefore, the torque equation is rewritten as where,   is the imaginary part of the stator current vector in the SFO reference frame.Therefore, if | *  | and  *  are stator flux magnitude and torque references, respectively, the imaginary part of the stator current reference vector which leads to torque and stator flux magnitude references is as follows: Because of the existence of a capacitor in the control system and its corresponding losses, when (47) is used to control the torque for all rotor speeds, the control system will not be accurate enough.In order to have a more accurate torque control, the block diagram shown in Figure 10 is used to calculate  *  .This block diagram is a combination of (47) and a simple PI controller.This controller removes the steady-state error and leads to an accurate control.The coefficients of this controller must be small enough so that it remains stable, and the ripples in  *  are limited, because a ripple of  *  results in a ripple of the generated torque.If only a PI is used,  *  will have a large value of ripple, and the torque response time will be very large.Therefore, in order to have a faster transient response of torque, it is necessary to use (47).In DTC method with VSI, the real part of voltage vector (V  ) in the stator flux reference frame leads to direct control of stator flux magnitude.
Therefore, using the real part of the stator current vector in the stator flux reference frame, it is possible to control the flux magnitude indirectly.In order to obtain the real part of the stator current reference, a simple PI controller is used.The input of this PI controller is the error that resulted from the comparison between the reference stator flux magnitude and the estimated flux value.In order to control the motor current in all cases, especially in the starting mode, the output of this PI controller must be limited between proper values.The output signal of the PI controller is the real part of reference stator current  *  .Therefore, the stator current reference vector that leads to control the IM stator flux and torque is as follows: The obtained reference current vector is in the stator flux reference frame.In order to generate it by SVM of the CSI, it must be transferred to the stationary reference frame.To do that, the stator flux position is used to transfer this vector from the stator flux reference frame to the stationary reference frame as follows: (50) This reference current vector must be generated in the CSI first and then applied to the motor after removing its harmonics with filtering capacitor.This current has a nearly sinusoidal waveform that leads to a nearly sinusoidal voltage in motor terminals after feeding the motor and passing through motor impedances.The speed and stator flux control of an IM with DTC method is shown in Figure 11.Separate PI controllers are used to control the speed and flux.The input to the speed controller is the error caused by comparing the reference speed and the rotor speed.The output of this controller is the imaginary part of stator current vector  *  , which is limited to a proper value.

Results and Discussion
7.1.Simulation Results.The simulation models of IOL, IFOC, and DTC control methods for CSI fed IM drive systems are developed and simulated using MATLAB simulation software.From the simulation results, dynamic performances of the developed drives have been individually analyzed mainly using the sensitivity tests.Here, the results of all the drives are compared with one another.The tests that are chosen for the comparison are step changes in the reference speed and torque of the motor drive.
The first test is step change in reference speed which is useful when trying to change from one speed to another.The next test is step change in reference torque that would be useful in identifying the controller with a high torque response.Change in speed reference is the first test that is done to compare the motor drives.
The results of this test can be seen in the Figures 12, 13, and 14.The speed command is varied from 0 to its rated value.The speed response of IOL drive is faster and has the settling time of 0.5 second range, while DTC and IFOC drives perform similarly but slower than IOL (settling time of less than 1 second).
The second test, as mentioned earlier, is the torque step test.The results of this test can be seen in Figures 12 to 14.The torque reference is varied from 0 to its rated value.The speed response of IFOC drive is faster and has the settling time of 0.2 seconds, while IOL and DTC drives perform similarly (settling time of less than 0.5 second).
During this test, the motor torque response of DTC drive is almost instantaneous within the 0.2 second range, while IFOC drive closely follows behind it.IOL and vector based / have slower torque response than all other Vector Control techniques.The above experimental result of IOL method of vector control technique in Figure 16 shows the fast settling time with feed forward control compared to feedback technique.Similarly the experimental result shown in Figure 17 depicts a greater torque response within 0.2 seconds compared to other vector control techniques.With this experimental verification, it is clear that IOL is adaptive for speed response and DTC is more suitable for torque response controlled drives.
From Tables 2 and 3, it seems that the dynamic speed and torque performance of all four drives using change in speed reference.It is evident from the table that IOL drive has the best speed response, with DTC and IFOC occupying positions next to it.It is also evident that DTC drive has the best instantaneous torque response, with IFOC and IOL next to it.

Conclusion
This paper presents a new approach in mathematical modeling of a 3-phase induction motor drive.Mathematics concepts are used in developing the simulation and implementation of IOL, IFOC, and DTC methods for CSI fed IM drives, In this paper, a rotor flux based reference frame control method for PWM CSI fed IM with simplified simulation circuits has been developed in which the motor speed and torque follow the references closely in all the three vector control methods.However a CSI fed induction motor drive has some drawbacks such as slower dynamic response and bulky, overall size due to the requirement of large smoothing Step increase in speed command Step decrease in speed command Step increase in torque command Step decrease in torque command Reference torque Step decrease in speed command Step increase in speed command Step decrease in torque command Reference torque Step increase in torque command inductor at its dc source end, and also a CSI drive system with on-line control strategy is more complex compared to a VSI drive due to the CSI gating requirements.But depending on some particular applications, its advantages can at times outweigh its disadvantages, hence rendering it as a suitable choice for certain applications.From the results, it is concluded that DTC drive is the top most performer for wide range of torque control applications (load varying  applications) with respect to change in speed as well as torque references.IOL drive is the good performer for speed control applications with respect to change in speed reference, and IFOC drive is the good performer for speed control applications with respect to change in torque reference.Hence for industrial drive applications with the requirement of torque Lm ) /(Lsl * Lrl + Lsl * Lm + Lrl * Lm) Lm /(Lsl * Lrl + Lsl * Lm + Lrl * Lm) Lm /(Lsl * Lrl + Lsl * Lm + Lrl * Lm) (Lsl + Lm)/(Lsl * Lrl + Lsl * Lm + Lrl * Lm)

Figure 1 :
Figure 1: Induction motor modeling in arbitrary reference frame.

Figure 3 :Figure 4 :Figure 5 :
Figure 3: Simulated results of "" phase and -axis stator currents in different rotating reference frames.

Figure 5 (
Figure5(c), it is evident that the stator -axis flux and stator axis current (control variables) have linear relationship with torque (speed).It is used in the stator flux oriented (SFO) control methods.The basic IM equations from Krause et al.[15] are simulated for three different reference frames using MATLAB 7.6 software and the results obtained are compared with each reference frame.From the comparison, it is observed that if the stator rotating reference frame is used, then stator -axis current is identical to that of the stator "" phase current as shown in Figure3(a), which is used in stator fed IM drives.Similarly, in rotor reference frame, rotor -axis current is identical to the rotor "" phase current and useful in rotor fed IM drives.If the synchronous reference frame is used, then the stator -axis current behaves like DC quantity as shown in Figure3(c).From Figures4(c) and 5(c), it is evident that synchronous reference frame is useful in RFO and SFO drive control methods.In this paper, synchronous reference frame is preferred to model the IM in /, IOL control methods.RFO control is used in IFOC method and SFO is used in DTC method.

Figure 7 :
Figure 7: Combined feed forward and integral feedback controller.

Figure 11 :
Figure 11: Block diagram schematic of the proposed DTC.

Figure 12 :
Figure 12: Speed, torque response of IOL control method.

Figure 15 :
Figure 15: Schematic diagram of the hardware setup of CSI fed IM.

Figure 16 :
Figure 16: Experimental result torque-speed response of IOL control method.

Table 1 :
Prototype motor parameters.1 Hp, 3 phase, star connected, 4 pole, 415 V, 1.8 A, 50 Hz The induction motor ratings are given in Table1.The inverter consists of CM75TL-12NF.The gate driving signal is developed by using Texas instruments DSP processor TMS320F2812 DSP.The software used to develop programs for TMS320F2812 DSP is by code composer studio and Matlab/Simulink.The processor contains the program that is downloaded to the computer.The clock speed of the DSP is 150 MHz, and it is capable of 32-bit operations.The on board available flash memory is 2.048 Mb.It was created specifically for motor control operation, and therefore Park's and Clark's transformations are conveniently built in.Another convenient feature is that it has sixteen 12-bit ADC pins that allow for a high degree of precision while taking many possible measurements.