Next Article in Journal
Multisensor RFS Filters for Unknown and Changing Detection Probability
Previous Article in Journal
Limitations of Harmonics Control in Power Converters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reduction of Induction Motor Energy Consumption via Variable Velocity and Flux References

by
Jorge Resa
1,*,†,
Domingo Cortes
1,†,
Juan Francisco Marquez-Rubio
1,† and
David Navarro
2,‡
1
Instituto Politecnico Nacional, ESIME Culhuacan, Mexico City 04440, Mexico
2
Escuela de Ingeniería y ciencias, Tecnologico de Monterrey, CCM, Mexico City 14380, Mexico
*
Author to whom correspondence should be addressed.
Avenida Santa Ana 951, Coyoacan 04440, CDMX, Mexico.
Calle del Puente 222, Tlalpan 14380, CDMX, Mexico.
Electronics 2019, 8(7), 740; https://doi.org/10.3390/electronics8070740
Submission received: 31 May 2019 / Revised: 19 June 2019 / Accepted: 26 June 2019 / Published: 29 June 2019
(This article belongs to the Section Systems & Control Engineering)

Abstract

:
Induction motors (IM) have been a fundamental part of industrial applications for over a century and the number of their applications continues to expand. A significant amount of the world’s total energy expense is consumed by this kind of motor. Hence, it is very important to increase the energy efficiency of these machines. Due to its good performance, field-oriented control (FOC) is the most common strategy to control IM. FOC requires references for stator current and rotor magnetic fluxes. For velocity regulation, a velocity reference is used instead of a stator current reference. However, at motor start-up or when a change of torque is required, it would be convenient for these references to be variable in order to reduce energy consumption. In this work, it is shown that this is indeed the case, and a technique to find optimal time-variable references for stator currents and magnetic rotor fluxes to reduce energy consumption is proposed. It is shown that, depending on the mechanical load, an energy reduction of 20–45% can be achieved.

1. Introduction

The reduction of energy consumption is a very important issue nowadays. Great effort is being dedicated to reducing the energy used by electric and electronic devices, particularly in two general cases: (a) in commonly used devices and (b) in devices fed by batteries. The first case is important because any improvement in energy expenditure in these devices has a multiplicative effect. The importance of the second case is obvious, as it is always desirable to maximize battery charge duration.
Induction motors are ubiquitous in industry, to the extent that more than 30% of the world’s total energy consumption comes from this kind of motor [1,2]. In addition, their use in electric transportation is increasing. Hence, any improvement in the energy expenditure of these machines is important.
Energy consumption optimization in three-phase induction motors has been a widely discussed topic since the first electronic converters were used to drive this type of motor [3,4]. With the introduction of the Field-Oriented Control (FOC), attention on this topic became deeper. This control strategy seeks to directly modifying the dynamics of the current that drives the motor. The FOC acts over the motor flux (d-axis) and current vector (q-axis). The significance of these vectors is due to the resultant torque from their product, and the energy consumption is proportional to their magnitude. Hence, it is natural that many works on optimization have relied on the FOC as a base [5].
Many previously published works have focused on obtaining optimal values for gains of controllers of the FOC. Some works have shown that a change in the gain modifies the resulting performance of the motor [6,7,8]. Other works have focused on the selection of an appropriate set-point for the optimal magnetic flux, given a certain load. Once the set-point is found, it is kept constant, even if the motor starts or stops numerous times [9,10,11]. In [12], the approach of minimizing the power losses in the stator to find the correct references for the controllers was employed. To find the set-point value, some studies have used numerical methods [13]. Others found the value using more heuristic techniques such as genetic algorithms or swarm algorithms [14].
A significant amount of energy is consumed by the motor when torque or velocity changes, particularly at the start-up of the motor or when an increase in torque is required. In this dynamic situation, to reduce energy consumption, it would be convenient for fluxes and current to be variable. However, most works keep the fluxes constant. Taking this observation into account, in this paper, a procedure to find references for flux and current vectors that minimize the motor energy consumption is proposed. The proposed algorithm seeks two time-variable functions that act as references for the flux and current vectors.
The proposed optimization procedure uses the whole dynamic system equations, including the motor dynamical model, a static or dynamical model for the mechanical load, and dynamic controller equations. The cost function to be minimized is constructed from the supplied currents, which are directly related to the total energy consumption.
It is important to point out that, for the optimization procedure, it is a controller strategy is assumed to be given. Thus, the controller parameters and final values of the set points are assumed to be known. There are some proposals to find optimal values for these parameters [15]. The procedure proposed here does not substitute, but rather complements, these proposals.
A control strategy that has proven to be very robust in simulations and in experiments [16] is employed in this paper. This control strategy includes proportional-integral (PI) controllers for the flux vector and the current vector. It also includes a third PI controller for the speed. Note that, although a control strategy is necessary to obtain the results, a different control strategy could be used in the optimization procedure.
As a result of the optimization procedure, time-variable vectors of flux and current are obtained as references. Using such references, the actual variables have peak values that are notably lower than those achieved without optimization. More significantly, average power at motor start-up is 20–45% lower.
The proposed procedure to find variable references requires a great number of calculations which must be performed before the motor operations. However, there are many applications where velocity changes are previously known. These applications include battery-fed motors, high peak starting current motors, and motors that continuously start and stop.
The rest of the paper is organized as follows: in Section 2, the overall setting for control of the induction motor is described. In Section 3, the optimization procedure is explained. Simulation results are presented in Section 4. The performance of the system before and after optimization is also analyzed in Section 4. In the last section, the benefits and limitations of the proposed procedure are discussed.

2. System Model

To obtain optimized references, it is necessary to build a complete model of the system. Such model can be derived by dividing the system into three parts: the induction motor model, the controller model including the magnetic flux estimator model, and the mechanical load model. In this work, a constant mechanical load is considered, but it would not be difficult to include a dynamic model for the load. In the following subsections, the motor model and the controller model are established.

2.1. The Model of the Induction Motor

The induction motor is an electro-mechanical device that has been studied extensively [4,5,17,18]. Basically, it consists of three stator coils, respectively fed by three (usually) sinusoidal voltages with 120 ° of phase between them. The current through the coils generated by the stator voltages results in a rotating electromagnetic field. This electromagnetic field induces currents in the rotor formed by short-circuited bars. Currents through bars on the rotor cause a second field to be generated. The interaction between the stator field and rotor field produces torque.
The mathematical model of the induction motor can be derived from relations between stator current equations, the magnetic rotor flux equations, and a final equation that relates the electric torque to the mechanical torque [17,18,19]. Some works have developed models from relations between the three phase stator voltages V a , b , c , stator currents i a , b , c , and magnetic fluxes ψ a , b , c .
It is well known that stator voltages and currents are not independent, and transformations can be made to describe a system using two instead of three phases. This is known as the Clarke transformation [20]. The two voltages (and currents) obtained after the Clarke transformation are separated by 90 ° and are known as the d q -frame.
Considering all the previous observations given in [17], the following model was obtained:
i ˙ d s s = L m 2 R r + L r 2 R s σ L r 2 L s i d s s + L m R r σ L r 2 L s ψ d r s + L m ω r σ L r L s ψ q r s + V d r s σ L s ,
i ˙ q s s = L m 2 R r + L r 2 R s σ L r 2 L s i q s s + L m R r σ L r 2 L s ψ q r s L m ω r σ L r L s ψ d r s + V q r s σ L s ,
ψ ˙ d r s = L m R r L r i d s s R r L r ψ d r s ψ q r s ω r ,
ψ ˙ q r s = L m R r L r i q s s R r L r ψ q r s + ψ d r s ω r ,
ω ˙ r = P 3 P L m ( i q s s ψ d r s i d s s ψ q r s ) ( 2 J ) 4 L r T L 2 J .
The symbols of model (1) are physical motor parameters and variables. They are defined in Table 1, except for σ , L s , and L r , which are introduced to simplify the notation and are given by
σ = 1 L m 2 L s L r ,
L s = L l s + L m ,
L r = L l r + L m .

2.2. Field-Oriented Control of the Induction Motor

The voltages, currents, and magnetic fluxes of induction motors are time-variable. However, if a rotating coordinated frame synchronized with the rotor magnetic flux is defined, for these coordinates, electric and magnetic variables would be constants. Having electric and magnetic variables as constants, the magnitude of these variables could be controlled with the usual PI controllers to achieve maximum torque. These are the key ideas of FOC. The overall scheme of FOC is depicted in Figure 1. In what follows, the main components of this schema are described.

2.2.1. Estimation of the Magnetic Flux Angle

To construct the rotating coordinate frame, it is necessary to know the direction of the magnetic flux of the rotor. The angle between the magnetic flux and the motor stator describes this direction. This angle is denoted by t h e t a (see Figure 1). However, this angle can not be measured; hence, it must be estimated.
Several flux angle estimators have been proposed. The Voltage Model uses supplied voltages and stator currents. The Current Model uses stator currents and the rotor angular velocity. There are also several sensor-less techniques that do not need the angular velocity of the rotor, such as observers or slip estimators. Each estimation method has its advantages and disadvantages. Sensor-less methods do not need the velocity or position sensor, but they require a higher number of calculations.
The estimation model used in this work is known as the Voltage Model [17]. In the voltage model, the stator magnetic flux d q components are obtained from
ψ ^ ˙ d s s = V d s s R s i d s s ,
ψ ^ ˙ q s s = V q s s R s i q s s .
Then, the rotor magnetic flux d q components are calculated using
ψ ^ d r s = L r ψ ^ d s s L l s i d s s L m L l r i d s s ,
ψ ^ q r s = L r ψ ^ q s s L l s i q s s L m L l r i q s s .
Having the d q components of rotor flux estimation, the magnitude of the rotor flux is calculated using
ψ ^ r s = ψ ^ d r s 2 + ψ ^ q r s 2 .
Finally, sin ( θ ) and cos ( θ ) are obtained from rotor magnetic flux estimation components and its magnitude:
cos ( θ ) = ψ ^ d r s ψ ^ r s ,
sin ( θ ) = ψ ^ q r s ψ ^ r s .
Expressions (6) are used to construct the rotating coordinate transformation and its inverse.

2.2.2. The Park Transformation

In order to use Proportional-Integral controllers, a rotating frame needs to be used. In this rotating frame, all of the electric and magnetic variables can be seen as stationary vectors. This frame can be aligned to any vector in particular. Generally, this rotation frame is aligned with a magnetic flux vector, which can be the rotor flux or the stator flux.
In summary, the motor model has a stationary coordinate frame (denoted with super index s), and the controller has a rotating coordinate frame. In order to be used in the controller, the variables from the motor model have to be transformed with the Park transformation given by
i q s i d s = cos ( θ ) sin ( θ ) sin ( θ ) cos ( θ ) i q s s i d s s .

2.2.3. The PI Controllers

To control the motor velocity via FOC, it is necessary to have three control loops. One control loop acts over the magnetic flux, a second control loop acts over the rotor current, and the third is a velocity control loop. The output of the velocity control loop can act as a reference for the rotor current control loop.
The control structure C ψ (see Figure 1) acts over the rotor flux, and it is described by
ψ e r = ψ r e f i d s ,
ψ ˙ i e = ψ e r ,
V d s = k a ψ ψ e r + k b ψ ψ i e ,
where ψ r e f is the reference of the magnetic flux.
Controller C ω acts over the angular velocity while controller structure C i acts over the rotor current, which produces torque. The output signal of C ω is the reference for C i .
The structure of controller C ω is given by
ω e r = ω r e f ω r ,
ω ˙ i e = ω e r ,
i s p = k a ω ω e r + k b ω ω i e ,
where ω r e f is the reference for the rotor angular velocity.
Controller C i uses i r e f as the reference, hence
i e r = i r e f i q s ,
i ˙ i e = i e r ,
V q s = k a i i e r + k b i i i e .
Finally, controller output voltages V d s and V q s must be converted from the rotating frame to the fixed frame using the inverse of the Park transformation given by
V q s s V d s s = cos ( θ ) sin ( θ ) sin θ ) cos ( θ ) V q s V d s .
The gains of the three PI controllers depend on motor parameters such as the inductance and resistance as well as other performance parameters such as the damping factor and bandwidth.
The gains for controller C i are given by
k a i = L s r ,
k b i = R r + R s L s r ,
where
L s r = L l s + L m 1 L m 2 L l r + L m L l s + L m .
The gains for controller for controller C ψ are
k a ψ = L s r ,
k b ψ = R s L s r ,
where
k p r = 3 P i ψ c L m 2 4 L l r + L m .
Finally, the gains for controller C ω are given by
k a ω = 1 k t k p r ,
k b ω = 1 δ 2 ,
where the parameters δ and k t are proposed by the user based on the damping ratio and velocity sensor time constant.

3. The Proposed Procedure to Find Variable References

This section describes how the references for velocity ω r e f and flux ψ r e f are found. A variational approach is used. Optimization based on variational tools has been studied in many references (see, for example, the excellent books [21,22,23])
The main result of variational approach used in this work is summarized below, and then the specific procedure is detailed.

3.1. A Variational Tool

Consider the system
x ˙ = f ( x , u ) ,
where x R n , u R m , y f : R n + m R n . To find an input u * that causes the system (17) to follow a trajectory x * that minimizes the functional
J ( u ) = h x ( t f ) + t o t f F ( x , u ) d t ,
it is necessary that, for all t t 0 , t f , expressions
x ˙ * = p H x * , u * , p * ,
p ˙ * = x H x * , u * , p * ,
0 = u H x * , u * , p * ,
and boundary conditions
x h x * ( t f ) , t f p * ( t f ) T δ x f = 0
hold. Function p ( t ) is a vector function in R n . The function H is called the Hamiltonian and is defined by
H x , u , p = F ( x , u ) + p T f ( x , u ) .
There are several methods to find u * , x * and p * that satisfy optimal conditions (19). The steepest descent is an iterative method that is easy to understand. The method is summarized in Figure 2. First, an initial discretized u ( t ) for t [ t 0 , t f ] is proposed. Using u, the estate x * is found by forward simulation, and the co-estate is found by backward simulation. At each iteration, small improvements are made by changing u according to
u ( k + 1 ) = u ( k ) τ u H ( k ) ,
where τ is the step size and can be constant or variable. If τ is constant, it must be smaller than the smallest valley. Changes in u given by Equation (22) seek to achieve
u H ( x , u , p ) = 0 .
Note that, at each iteration, conditions (19a) and (19b) are satisfied. Then, if condition (19c) is fulfilled, it means that all optimal conditions (19) are accomplished. That is, when Equation (23) is accomplished, conditions (19) are fulfilled.
The algorithm stops when the left-hand-side of Equation (23) is close enough to zero.

3.2. Finding Variable References

The optimization procedure seeks to find references for the velocity and flux; hence, it is convenient to define
u = u ω u ψ T = ω r e f ψ r e f T .
The main objective of this work is the reduction of energy consumption. Consequently, it is convenient to propose a function cost that depends on motor stator currents because these currents are directly related to the motor energy consumption.
Thus, the cost function is proposed as
J ( u ) = g ω ω r e f ω r ( t f ) 2 + t o t f i d s s + i q s s 2 d t .
Note that a big value of g ω reinforces the fact that the final value of the velocity has been reached.
The overall system model is composed of the motor model and the controller model. Hence, the complete system dynamical equation is formed by the five differential equations of the motor given by Equation (1), the two flux estimator differential equations given by Equation (3), and the three differential equations of the PI controllers given by Equations (8)–(10). It is convenient to define the motor state using
x m = i d s s i q s s ψ d r s ψ q r s ω r T .
The estimator state is defined by
x e = ψ ^ d r s ψ ^ q r s T
and the controllers state is defined by
x c = ω i e i i e ψ i e T .
Hence, the complete system state can be defined by
x = x m x e x c T .
According to Section 3.1, from Equations (25)–(28), the Hamiltonian can be expressed by
H ( x , u , p ) = i d s s 2 + i q s s 2 + p x , u T x ˙ .
The co-state dynamical equations are obtained using
p ˙ = H x .
It is important to point out that, in addition to dynamical equation, the system description includes algebraic equations such as the Park and Inverse Park transformations (see Figure 1). Although, theoretically, these algebraic equations could be considered as additional restrictions, it is not convenient to do so. Including algebraic equations as restrictions means that the numerical algorithm has to solve them at each iteration. This significantly increases the time required for algorithm execution, and means that it is prone to errors and could result in numerical instability.
Instead of including algebraic relations as restrictions, they were solved and included in the system dynamical equations. Such dynamical equations were used in Equations (30) and to find the co-state equations from Equation (31). Note that the inclusion of algebraic equations in dynamical equations gives rise to very long expressions that are not written here because they do not add clarity. However, they are not difficult to obtain using symbolic computation software, like Wolfram Mathematica (WM) (version 10.0.0.0, Wolfram Research Inc., Champaign, IL, USA).
According to Equation (20), the boundary values of co-states are zero, except for
p ω ( t f ) = ω r e f ω r ( t f ) .
After determining all system and co-state dynamical equations and their boundary conditions, the steepest descent was applied to find the optimal functions for u. The initial time was zero and the final time t f was proposed from a simulation without optimization. The norm employed was
t 0 t f H u ω 2 + H u ψ 2 d t .

4. Simulation Results

Using WM, the ten differential equations of the complete system were obtained in explicit form. The developed code used the motor model, the flux estimator, and the controller equations. The ability of WM to manipulate symbolic expressions allows different dynamic systems with different parameters to be obtained with almost no effort. To actually obtain the references, optimization and simulation procedures were developed in Java. Methods to evaluate the cost functions and norm functions were also programmed. A translator between expressions produced by WM and expressions with the syntax of Java was also developed. The optimization algorithm was programmed in Java because calculations take less time to perform on this platform than in WM, and has more tools for debugging. Hence, it is easier to find and solve potential problems. Once the optimized references were obtained in Java, these functions were fed back to the motor model in WM. With this strategy, we ensured that well performing optimization and simulation procedures were developed.
The motor parameters used in the simulations are shown in Table 2. Such parameters were obtained from the 3/4 HP induction motor. The optimization parameters were τ = 0.01 and m a x ( H / u ) = 5 ( 10 6 ) .
Several simulations were performed to observe the behaviors of the supply currents and the angular velocity. The final angular velocity was restricted to keep voltage supply below 220VAC. In all simulation conditions, the current consumption and peak value reduction were observed. The percentage of energy consumption reduction was between 20 and 45% depending on the angular velocity and mechanical load.
In Figure 3 and Figure 4, two of these simulations are shown. As can be observed from these simulations, the optimization procedure changed the references. As a consequence of this change, peak values were notably reduced. Particularly, it can be observed from the graphs that optimization of references u ω ( t ) and u ψ ( t ) reduced the mechanic stress by reducing the peak angular velocity value and smoothing the acceleration. Optimization also reduced the peak values of i d s s and i q s s , which are the optimization target variables. If the motor demands less current, the required voltage will be less. Reductions in the current and voltage peak values mean less peak power is needed.
The level of reduction achieved is summarized in Table 3 and Table 4. To obtain the values reported in these tables, the current consumption was calculated using
t 0 t f i d s s 2 + i q s s 2 d t
and the energy consumed was calculated by
3 2 t 0 t f i d s s V d s s + i q s s V q s s d t .
The average apparent power delivered by the power source was calculated from
A P A v = 3 2 V d 2 + V q 2 i d 2 + i q 2 ,
where V d , V q , i d and i q are the root mean square (RMS) values [24].
Note that, instead of using the commonly used acceleration ramp, the proposed algorithm seeks the optimization of current consumption by applying large initial values at the set points and then rapidly converging them to the final value, which is equal to the original set-point value. The large initial set-point values help the system to have a faster response, but this faster response does not imply big peak values.

5. Conclusions

It has been shown in this work that, if a Field-Oriented Controlled induction motor has variable references for angular velocity and rotor magnetic flux, it can reduce the energy consumed by 20–45% and decrease the current consumption by 10–15%. A procedure to find these references was proposed. Variable current and flux references also reduce peak current values. This effect can reduce the stress over the electric installation; furthermore, acceleration peaks are lower than in the initial simulation reducing the mechanic stress of the motor and the load.
Due to its massive number of calculations, the proposed procedure to find variable references cannot be done during motor operation. If there is a change in parameters or mechanical load, the calculations must be redone. A plausible strategy is to perform the calculations for several common cases and use interpolation if another case arises during the operation. Nevertheless, there are significant motor applications where changes in velocity are known and the motor starts and stops continuously.

Author Contributions

Author contribution: conceptualization, J.R. and D.C.; methodology, J.R. and D.C.; software, J.R. and D.N.; validation, J.R. and D.N.; formal analysis, J.R. and D.C.; investigation, J.R. and D.C.; resources, D.C.; data curation, F.M.-R.; writing–original draft preparation, J.R.; writing–review and editing, F.M.-R. and D.N.; visualization, J.R.; supervision, D.C. and F.M.-R.; project administration, D.C.; funding acquisition, D.C.

Funding

This research was funded by Instituto Politecnico Nacional under grant number SIP-20181542.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Waide, P.; Brunner, C.U. Energy-Efficiency Policy Opportunities for Electric Motor-Driven Systems; IEA Energy Papers 2011/7; OECD Publishing: Paris, France, 2011. [Google Scholar]
  2. Ferreira, F.J.T.E.; de Almeida, A.T. Overview on energy saving opportunities in electric motor driven systems—Part 1: System efficiency improvement. In Proceedings of the 2016 IEEE/IAS 52nd Industrial and Commercial Power Systems Technical Conference (I CPS), Detroit, MI, USA, 1–5 May 2016; pp. 1–8. [Google Scholar] [CrossRef]
  3. Stumper, J.; Kennel, R. Real-time dynamic efficiency optimization for induction machines. In Proceedings of the 2013 American Control Conference, Washington, DC, USA, 17–19 June 2013; pp. 6589–6594. [Google Scholar] [CrossRef]
  4. Bazzi, A.M.; Krein, P.T. Review of Methods for Real-Time Loss Minimization in Induction Machines. IEEE Trans. Ind. Appl. 2010, 46, 2319–2328. [Google Scholar] [CrossRef]
  5. Ghozzi, S.; Jelassi, K.; Roboam, X. Energy optimization of induction motor drives. In Proceedings of the 2004 IEEE International Conference on Industrial Technology, 2004, Hammamet, Tunisia, 8–10 December 2004; Volume 2, pp. 602–610. [Google Scholar] [CrossRef]
  6. Thangaraj, R.; Chelliah, T.R.; Pant, M.; Abraham, A.; Grosan, C. Optimal gain tuning of PI speed controller in induction motor drives using particle swarm optimization. Log. J. IGPL 2011, 19, 343–356. [Google Scholar] [CrossRef]
  7. Mitsukura, Y.; Yamamoto, T.; Kaneda, M. A design of self-tuning PID controllers using a genetic algorithm. In Proceedings of the 1999 American Control Conference (Cat. No. 99CH36251), San Diego, CA, USA, 2–4 June 1999; Volume 2, pp. 1361–1365. [Google Scholar] [CrossRef]
  8. Mohamed, H.A.F.; Lau, E.L.; Yang, S.S.; Moghavvemi, M. Fuzzy-SMC-PI Flux and Speed Control for Induction Motors. In Proceedings of the 2008 IEEE Conference on Robotics, Automation and Mechatronics, Chengdu, China, 21–24 September 2008; pp. 325–330. [Google Scholar] [CrossRef]
  9. Taheri, A.; Rahmati, A.; Kaboli, S. Efficiency Improvement in DTC of Six-Phase Induction Machine by Adaptive Gradient Descent of Flux. IEEE Trans. Power Electron. 2012, 27, 1552–1562. [Google Scholar] [CrossRef]
  10. Hung, N.T.; Thien, N.C.; Nguyen, T.P.; Le, V.S.; Tuan, D.A. Optimization of Electric Energy in Three-Phase Induction Motor by Balancing of Torque and Flux Dependent Losses; AETA 2013: Recent Advances in Electrical Engineering and Related Sciences; Zelinka, I., Duy, V.H., Cha, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2014; pp. 497–507. [Google Scholar]
  11. Melo, P.; de Castro, R.; Araujo, R.E. Evaluation of an Energy Loss-Minimization Algorithm for EVs Based on Induction Motor. In Induction Motors; IntechOpen: Rijeka, Croatia, 2012; Chapter 17. [Google Scholar] [CrossRef] [Green Version]
  12. Thanh Hung, N.; Chi Thien, N.; Phuong Nguyen, T.; Le, V.S.; Tuan, D.A. Optimization of Electric Energy in Three-Phase Induction Motor by Balancing of Torque and Flux Dependent Losses. Lect. Notes Electr. Eng. 2014, 282, 497–507. [Google Scholar] [CrossRef]
  13. Borisevich, A. Numerical Method for Power Losses Minimization of Vector-Controlled Induction Motor. Int. J. Power Electron. Drive Syst. 2015, 6, 486–497. [Google Scholar] [CrossRef]
  14. Ta, C.M.; Hori, Y. Convergence improvement of efficiency-optimization control of induction motor drives. IEEE Trans. Ind. Appl. 2001, 37, 1746–1753. [Google Scholar] [CrossRef]
  15. Ismail, M. Efficiency optimisation with PI gain adaptation of field-oriented control applied on five phase induction motor using AI technique. Int. J. Model. 2013, 20, 344–360. [Google Scholar] [CrossRef]
  16. Wilson, D. Teaching Your PI Controller to Behave. 2015. Available online: http://www.webcitation.org/78mq5YtqU (accessed on 31 May 2019).
  17. Bose, B. Modern Power Electronics and AC Drives; Prentice Hall: Upper Saddle River, NJ, USA, 2001. [Google Scholar]
  18. Krishnan, R. Electric Motor Drives: Modeling Analysis In addition, Control; Prentice-Hall Of India Pvt. Limited: New Delhi, India, 2008. [Google Scholar]
  19. De Doncker, R.; Pulle, D.; Veltman, A. Advanced Electrical Drives: Analysis, Modeling, Control; Power Systems; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  20. Duesterhoeft, W.C.; Schulz, M.W.; Clarke, E. Determination of Instantaneous Currents and Voltages by Means of Alpha, Beta, and Zero Components. Trans. Am. Inst. Electr. Eng. 1951, 70, 1248–1255. [Google Scholar] [CrossRef]
  21. Kirk, D. Optimal Control Theory: An Introduction; Dover Books on Electrical Engineering; Dover Publications: Englewood Cliffs, NJ, USA, 2004. [Google Scholar]
  22. Cassel, K. Variational Methods with Applications in Science and Engineering; Cambridge University Press: New York, NY, USA, 2013. [Google Scholar]
  23. Liberzon, D. Calculus of Variations and Optimal Control Theory: A Concise Introduction; Princeton University Press: Princeton, NJ, USA, 2012. [Google Scholar]
  24. Reginatto, R.; Ramos, R.A. On electrical power evaluation in dq coordinates under sinusoidal unbalanced conditions. IET Gener. Transm. Distrib. 2014, 8, 976–982. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Velocity Field-Oriented Control.
Figure 1. Velocity Field-Oriented Control.
Electronics 08 00740 g001
Figure 2. Flow diagram of the steepest descent method.
Figure 2. Flow diagram of the steepest descent method.
Electronics 08 00740 g002
Figure 3. Simulation with load. Final velocity: 30 rads/s. Final time: 0.5 s. (a) motor power; (b) motor angular velocity; (c) motor supply current d-axis; (d) motor supply current q-axis; (e) angular velocity reference; (f) magnetic flux reference.
Figure 3. Simulation with load. Final velocity: 30 rads/s. Final time: 0.5 s. (a) motor power; (b) motor angular velocity; (c) motor supply current d-axis; (d) motor supply current q-axis; (e) angular velocity reference; (f) magnetic flux reference.
Electronics 08 00740 g003
Figure 4. Simulation with load. Final velocity 50 rads/s. Final time: 0.5 s. (a) motor power; (b) motor angular velocity; (c) motor supply current d-axis; (d) motor supply current q-axis; (e) angular velocity reference; (f) magnetic flux reference.
Figure 4. Simulation with load. Final velocity 50 rads/s. Final time: 0.5 s. (a) motor power; (b) motor angular velocity; (c) motor supply current d-axis; (d) motor supply current q-axis; (e) angular velocity reference; (f) magnetic flux reference.
Electronics 08 00740 g004
Table 1. Induction motor parameters.
Table 1. Induction motor parameters.
ParameterDescription
i d s s Supply current d-axis
i q s s Supply current q-axis
ψ d r s Rotor magnetic flux d-axis
ψ q r s Rotor magnetic flux q-axis
V d r s Supply voltage d-axis
V q r s Supply voltage q-axis
ω r Motor angular velocity
T L Mechanic load torque
R r Rotor resistance
R s Stator resistance
L l s Stator inductance
L l r Rotor inductance
L m Mutual inductance
JRotor inertia
PNumber of poles
Table 2. Motor parameters used in simulations.
Table 2. Motor parameters used in simulations.
ParameterValue
P o w e r 559.27 W
P2
R r 21.34 Ω
R s 4.19 Ω
L m 1.37 H
L l s 0.05 H
L l r 0.05 H
T L 1.0 Nm
J5.89 × 10 4 kg·m 2
Table 3. Simulation with load. Final velocity: 30 rads/s. Final time: 0.5 s.
Table 3. Simulation with load. Final velocity: 30 rads/s. Final time: 0.5 s.
ParameterOriginalOptimized%
Peak Current3.62 A1.56 A 56.9 %
Peak Power352.35 W70.18 W 80.08 %
Energy34.83 J18.50 J 46.88 %
Apparent Average Power103.73 VA51.59 VA 50.26 %
Current Consumption1.32 A 2 · s1.13 A 2 · s 14.4 %
Table 4. Simulation with load. Final velocity: 50 rads/s. Final time: 0.5 s.
Table 4. Simulation with load. Final velocity: 50 rads/s. Final time: 0.5 s.
ParameterOriginalOptimized%
Peak Current3.54 A2.62 A 25.9 %
Peak Power378.01 W199.73 W 47.16 %
Energy42.52 J32.58 J 23.38 %
Apparent Average Power143.76 VA110.01 VA 23.47 %
Current Consumption1.34 A 2 · s1.24 A 2 · s 7.46 %

Share and Cite

MDPI and ACS Style

Resa, J.; Cortes, D.; Marquez-Rubio, J.F.; Navarro, D. Reduction of Induction Motor Energy Consumption via Variable Velocity and Flux References. Electronics 2019, 8, 740. https://doi.org/10.3390/electronics8070740

AMA Style

Resa J, Cortes D, Marquez-Rubio JF, Navarro D. Reduction of Induction Motor Energy Consumption via Variable Velocity and Flux References. Electronics. 2019; 8(7):740. https://doi.org/10.3390/electronics8070740

Chicago/Turabian Style

Resa, Jorge, Domingo Cortes, Juan Francisco Marquez-Rubio, and David Navarro. 2019. "Reduction of Induction Motor Energy Consumption via Variable Velocity and Flux References" Electronics 8, no. 7: 740. https://doi.org/10.3390/electronics8070740

APA Style

Resa, J., Cortes, D., Marquez-Rubio, J. F., & Navarro, D. (2019). Reduction of Induction Motor Energy Consumption via Variable Velocity and Flux References. Electronics, 8(7), 740. https://doi.org/10.3390/electronics8070740

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop