Simple Hybrid Model for Efficiency Optimization of Induction Motor Drives with Its Experimental Validation

New hybridmodel for efficiency optimization of inductionmotor drives (IMD) is presented in this paper. It combines two strategies for efficiency optimization: loss model control and search control. Search control technique is used in a steady state of drive and loss model during transient processes. As a result, power and energy losses are reduced, especially when load torque is significant less related to its rated value. Also, this hybrid method gives fast convergence to operating point of minimal power losses and shows negligible sensitivity to motor parameter changes regarding other published optimization strategies. This model is implemented in vector control induction motor drive. Simulations and experimental tests are performed. Results are presented in this paper.


Introduction
Induction motor is a widely used electrical motor and a great energy consumer.The vast majority of induction motor drives are used for heating, ventilation, and air conditioning (HVAC).These applications require only low dynamic performance, and in most cases only voltage source inverter is inserted between grid and induction motor as cheapest solution.The classical way to control these drives is constant V/f ratio, and simple methods for efficiency optimization can be applied [1].From the other side there are many applications where, like electrical vehicles, electric energy has to be consumed in the best possible way and use of induction motors.These applications require an energy optimized control strategy [2].
One interesting algorithm which can be applied in a drive controller is algorithm for efficiency optimization.
In a conventional setting, the field excitation is kept constant at rated value throughout its entire load range.If machine is underloaded, this would result in overexcitation and unnecessary copper losses.Thus in cases where a motor drive has to operate in wider load range, the minimization of losses has great significance.It is known that efficiency improvement of IMD can be implemented via motor flux level and this method has been proven to be particularly effective at light loads and in a steady state of drive.Also flux reduction at light loads gives less acoustic noise derived from both converter and machine.From the other side low flux makes motor more sensitive to load disturbances and degrades dynamic performances [3].
Drive loss model is used for optimal drive control in loss model control (LMC) [3][4][5][6][7].These algorithms are fast because the optimal control is calculated directly from the loss model.But power loss modeling and calculation of the optimal control can be very complex.Often the loss model is not accurate enough.
Search strategy methods have an important advantage compared to other strategies [8][9][10][11].It is completely insensitive to parameter changes, while effects of the parameter variations caused by temperature and saturation are very expressed in two other strategies.The online efficiency optimization control on the basis of search, where the flux is decremented in steps until the measured input power settles down to the lowest value, is very attractive.Algorithm is applicable universally to any motor.Besides all good characteristics of search strategy methods, there is an outstanding Advances in Power Electronics problem in its use.For many applications flux convergence to its optimal value is too slowly.Also, flux is never reached its optimal value then in small steps oscillates around it.
For electrical drives that work in periodic cycles, it is possible to calculate the optimal trajectory of magnetization flux, using optimal control theory, so that power losses in one working cycle are minimal [12].These methods give good results if the working conditions do not change.
Hybrid method combines good characteristics of two optimization strategies SC and LMC [3,[13][14][15].It was enhanced attention as interesting solution for efficiency optimization of controlled electrical drives.During transient process LMC is used, so fast flux changes and good dynamic performances are kept.Search control is used for efficiency optimization in a steady state of drive.Loss model of IM in d-q rotational system and procedure for parameter identification in a loss model based on Moore-Penrose pseudoinversion is given in Section 2. New hybrid model is presented in Section 3. Qualitative analyses of this method with simulation and experimental results are given in Section 4. At the end, obtained results are presented in Conclusion.

Power Loss Modeling
The process of energy conversion within motor drive converter and motor leads to the power losses in the motor windings and magnetic circuit as well as conduction and commutation losses in the inverter [6].
The losses in the motor consist of hysteresis losses and eddy current losses in a magnetic circuit (iron losses), losses in the stator and rotor windings (copper loss), and stray losses.In nominal operating conditions the iron losses are typically 2-3 times smaller than the copper losses, but at low loads, these losses are dominant.These losses consist of hysteresis and eddy current losses.Eddy current losses are proportional to the square of supply frequency, and hysteresis losses are proportional to supply frequency.Both components of iron losses are dependent of stator flux level, so next expression is suitable to represent these losses: where  eddy is eddy current and  hys is hysteresis loss coefficients.
Copper losses are appeared as a result of the passing the electric current through the stator and rotor windings.These losses are proportional to the square of current through stator and rotor windings, and they are given by The total additional losses typically do not exceed 5% when the drive works with light loads.This case is the most important for the power loss minimization algorithms, so stray losses are not considered as a separate loss component.
Losses in the drive converter consist of the losses in the rectifier and the conductive and switching losses in the inverter.The losses in the rectifier are independent of the magnetizing flux and not specifically taken into account.Only conductive losses in the inverter are dependent on the magnetizing flux, and these can be presented in the next form: Based on the previous considerations, the losses in the induction motor drive, dependent on the magnetizing flux, can be expressed as follows: Take into account expression for output power: where  out is output power and  is constant which depends on the characteristics of Park's rotating transformation (in this case it is 1.5), and the relation expression for input power can be written in the next form: and  is positive constant.
Two typical cases are differed: (1) linear dependence of magnetizing flux from the magnetizing current, (2) nonlinear dependence of magnetizing flux from the magnetizing current.
In the algorithms for loss minimization, magnetizing flux is less than or equal to the nominal value, so it is used in the linear part of magnetization characteristics.
Starting from the expressions (4) and taking into account expression (5), the power losses can be expressed as a function of   ,  em , and   as follows: Slip angular speed can be defined as follows: where   is time rotor constant.
Based on expression (8) power losses can be given as function of magnetizing current   and operating conditions (,  em ): ), and  3 = (2  +  2   ) em /  (10) can be written as follows: Parameter  3 is a function of   which is time variant especially due to temperature changes.Time rotor constant is continuously updated in the algorithm for parameter identification, so the parameter  3 too.First and second derivations of   in a function of   are Value   max is chosen to be   , and   min is determined based on characteristics of drive and requirements for electromagnetic torque reserve.In this case   min = 0.5  , where   is nominal value of   current.On the basis of expressions in (12), it can be concluded that in the steady state of drive and known operating conditions exists accurate value of   which gives minimal losses as follows: Based on the specified optimal value of  component stator current vector and command of electromagnetic torque,  component of the stator current is determined as follows:

Parameter Identification in a Loss
Model.This procedure of parameter identification is a mainly based on the procedures described in [6,7].In order to determine the parameters in the model losses, , ,  1 , and  2 , it is necessary to measure the input power and know the exact value of the output power.The exact value of the output power is needed to determine correct information of the losses in the drive and to avoid coupling between load pulsation and work of optimization controller.Also, the samples of the values that define input powers V DC ,  DC ( in = V DC  DC ) are available in the controller, and these values are used for protect and control functions (dynamic braking, the soft-start operation, PWM compensation, etc.).The proposed procedure of parameter identification is shown in Figure 1.The samples that exist in the expression ( 7) are memorized at each period.Usually the   is 100-200 s.Due to high-frequency components that do not contribute to identification of the relevant parameters  = [   1  2 ]  , input signal and the input power are averaged in the period  =   : Values   ,   ,  1 ,  2 , and   ,  = 1, . . .,  successively form the columns (:, 1), (:, 2), (:, 3), (:, 4), and (:, 5) of matrix  ×5 .
Vector  is formed from  averaged values of input power   ,  = 1, . . ., .Vector   is calculated by applying Moore-Penrose pseudoinversion (Figure 1) as follows: That is   is approximate solution of the matrix equations  = , which gives a minimum value ‖ − ‖ or minimum mean-square error.
The choice of  is of great importance for the process of parameters identification (see (15)) because the input power fluctuations should be sufficient in the period  =   .Otherwise    matrix will be singular or aspire singularity, that is, det (  ) will be close to 0. The process of determining the parameters in the loss model is repeated during the operation of electric drive.

Hybrid Model for Efficiency Optimization
Electrical drive with block for efficiency optimization is shown in Figure 2. Electric drive is supplied from the primary power network 3×380 V.This voltage is rectified, and the voltage and current in DC link are measured.The drive inverter is current regulated (CR) voltage three-phase inverter.Reference and measured  and  components of stator current are kept in the current controllers.These controllers are realized in rotational ,  coordinate system as a linear PI controller.Outputs of controller are  and  component of stator voltage.These voltages are transformed (B-transformation) in threephase reference voltages V *  , V *  , and V *  .These voltages are scaled and lead to PWM modulator, where control signals of inverter switches are generated.This drive works as speed controlled drive.Speed reference and measured speed are led into speed controller.Speed controller is realized as PI controller in incremental form, with proportional coefficient in a feedback local branch.The output of speed controller is reference value of electromagnetic torque.Reference value of stator current  component is determined in a block for efficiency optimization and  component of stator current  in Figure 1: A method for determining the parameters , ,  1 ,  2 and,  from the input signal.
vector on the basis of  stator current vector component and electromagnetic torque reference (see ( 14)).Position of rotor flux vector is determined in indirect vector control block (IVC).Hybrid model for efficiency optimization consists from 3 blocks, LMC, SC, and steady state control (SSC) block.LMC is used during transient states caused of external speed or torque demand [3].Optimal control ( * LMC ,  * LMC ) is calculated directly from loss model for a given operational conditions what obtains power loss optimization and good dynamic performances.SC is used in a steady state for a constant output power.
On the basis of speed reference and measured speed, SSC block defines its output and controls switches (Figure 2).If transient states is detected, LMC is active, and its outputs are forwarded to indirect vector control (IVC) block and current regulators.When steady state is detected in SSC block, last value of magnetizing current during transient state is used as starting point for search algorithm.

Loss Model Controller.
Optimal control calculation in LMC for a given operational conditions is described in Section 2.
Expressions for  and  component of stator current are defined by ( 8) and (10).
This method is sensitive to parameter changes due to temperature changes and magnetic circuit saturation, what consequently leads to error in a current references calculation.So, algorithm for parameter identification is always active, and parameters in the loss model are continuously updated (Figure 2).

Search Controller.
Search algorithm is used in steady state, which is detected in the SSC block.Error that exists between the current reference   that is generated in the LMC model and in the search model appears as a consequence of inverter and stray losses which are not included in the model.The applied search algorithm is simple.Since the current   is very close to the value which gives minimal losses, small step of magnetization current Δ  = 0.01  is chosen.
For two successive values of the   current, power losses are determined.Sign of Δ  is maintained if power losses are reduced.Otherwise, the sign of Δ  is opposite in the next step: When the two values of magnetization current  1 and  2 were found, so the sign of power loss is changed between these values, and new reference of   current is specified as: In this way, there are no oscillations of   current, air gap flux, and electromagnetic torque, which are characteristics of the search algorithm.    of testing the drive behavior in the dynamic mode and its robustness for a sudden load perturbations.Graphs of magnetization current (  ) and active component (  ) of stator current vector for a given operating conditions and applied hybrid method are shown in Figure 4.

Simulation and Experimental Results
Graphs of power loss for nominal flux and hybrid method are given in Figure 5.It is obvious that the use of hybrid methods gives significantly less power losses than when the flux is maintained at nominal value for a light load of drive.

Experimental
Results.The experimental tests have been performed on the setup which consists of the following: (i) induction motor (3 phases, Δ380 V/220 V, 3.7/2.12A, cos  = 0.71, 1410 r/min, 50 Hz), (ii) incremental encoder connected with the motor shaft, (iii) three-phase drive converter (DC/AC converter and DC link), (iv) PC and dSPACE1102 controller board with TMS320C31 floating point processor and peripherals, (v) interface between controller board and drive converter.Parameters of used induction motor are given in the appendix of the paper.The algorithm observed in this paper used the MATLAB-Simulink software, dSPACE real-time interface, and C language.Handling real-time applications is done in ControlDesk.
All experimental tests and simulations have been done in the same operating conditions of the drive, and some comparisons between algorithms for efficiency optimization are made through the experimental tests.Graphs of power losses for nominal flux and when the hybrid method is applied are shown in Figures 6 and 7, respectively.Graphs of magnetization current for a given operational conditions are shown in Figure 8, and speed response is shown in Figure 9.  Based on graph in Figure 6, it can be concluded that the speed drop for a step change of load is small and speed response is strictly aperiodic.

Conclusion
Algorithms for efficiency optimization of induction motor drive are briefly described.Detailed theoretical analysis of power losses in induction motor is presented.Algorithm for parameter identification in the loss model based on Moore-Penrose pseudoinversion is presented.Also, hybrid algorithm for efficiency optimization has been applied.According to the theoretical analysis, performed simulations, and experimental tests we have arrived to the following conclusions.
If load torque has a value close to nominal or higher, magnetizing flux is also nominal regardless of whether an algorithm for efficiency optimization is applied or not.For a light load, hybrid method for efficiency optimization gives significant power loss reduction (Figures 5, 6, and  7).
There is no oscillation of magnetization flux which is characteristic of the search algorithms (Figures 4 and 8).
Also, this hybrid method shows good dynamic performances and no sensitivity to parameter changes (Figure 9).Implementation of presented algorithm is simple, and it can be universally applied to any electrical motor.Changes are only related to different models of power losses.

4. 1 .
Simulation Results.Hybrid method for efficiency optimization of IMD has been studied through computer simulation in MATLAB-Simulink.Speed reference and load torque are shown in Figure 3.The steep change of load torque appears with the aim  in

Figure 2 :
Figure 2: Overall proposed block diagram of efficiency optimization controller in IMD.

𝑖
and   stator currents Current   p.u.Current   p.u.

Figure 4 :Figure 5 :
Figure 4: Graphs of  and  component of stator current for a given operating conditions.

Figure 6 :
Figure 6: Power losses for nominal magnetization flux and operational conditions shown in Figure 3.

Figure 7 :
Figure 7: Power losses when the hybrid method is applied and operational conditions shown in Figure 3.

Figure 8 : 3 Figure 9 :
Figure8: Magnetization current   and active current when the hybrid method is applied and operational conditions shown in Figure3.