Current Control and Reactive Power Minimization of a Direct Matrix Converter Induction Motor Drive With Modulated Model Predictive Control
Current Control and Reactive Power Minimization of a Direct Matrix Converter Induction Motor Drive With Modulated Model Predictive Control
Abstract— This paper investigates a novel model- predictive above all in high performance drives that require more and
control strategy to control an induction motor fed by a matrix more performing regulation.
converter with multi-objective control capability. The proposed
control method combines the features of the classical Model
Rd Matrix converter
Predictive Control and the Space Vector Modulation technique VA Isa
into a Modulated Model Predictive Control. This new solution
maintains all the characteristics of Model Predictive Control Cf
(such as fast transient response ,multi-objective control using VB Isb SAa SAb SAc
only one feedback loop, easy inclusion of nonlinearities and
constraints of the system, the flexibility to include other system Lf
requirements in the controller) adding the advantages of working VC Isc SBa SBb SBc
at fixed switching frequency and improving the quality of the
controlled waveforms. Simulation results employing the control
method to control the stator current and input reactive power of Bidirectional
switch SCa SCb SCc
a direct matrix converter induction motor drive are presented.
Vs, Is Rd
As shown in Fig.2, the control strategy consists of two Vi, Ii
cascaded loops for speed and current-reactive power control. Lf
The speed of the machine is conrolled by a classical
proportional-integral (PI) controller whereras the stator Il Cf
currents in the dq coordinate frame and input reactive power
are controller by the M2PC algorithm.
In order to achieve a modulation-like behavior the M2PC uses
Fig. 3 Schematic of the input filter
the vector sequence employed in the SVM. Unlike MPC
where one voltage vector is applied for the entire sampling
The state space model in continuous time domain is given by co-ordinates and 𝑇𝑠 is the sampling time for the discretization.
the set of equations in (3) The stator currents can then be predicted in order to form the
cost function. The current prediction is also done in the dq
reference frame since a field oriented model is used.
−1 1 1 1
−
𝑅𝑑 𝐶𝑓 𝐶𝑓 𝑅𝑑 𝐶𝑓 𝐶𝑓
𝑥̇ (𝑡) = [ −1 ] 𝑥(𝑡) + [ −1 ] 𝑢(𝑡) (3) 𝑖𝑠𝑑 (𝑘 + 1) = 𝜖𝑉𝑠𝑑 (𝑘) + 𝛾 𝑖𝑠𝑑 (𝑘)
0 0
𝐿𝑓 𝐿𝑓
𝜆
where +𝑇𝑠 [𝜇 𝐼𝑠𝑞 2 + 𝜔𝑒 𝑖𝑠𝑞 (𝑘) + ] (10)
𝜏𝑟
𝑉𝑖 (𝑡) 𝑉 (𝑡)
𝑥(𝑡) = [ ] 𝑎𝑛𝑑 𝑢(𝑡) = [ 𝑠 ] 𝑖𝑠𝑞 (𝑘 + 1) = 𝜖𝑉𝑠𝑞 (𝑘) + 𝛾𝑖𝑠𝑞 (𝑘)
𝐼𝑙 (𝑡) 𝐼𝑖 (𝑡)
Substituting (8) in (7), the expression of the rotor flux is 𝑖𝑠𝑑 (0) (𝑘 + 1) = 𝛾 𝑖𝑠𝑑 (𝑘)
obtained in (9).
𝜆
𝐿𝑠 𝐿𝑟 𝐿𝑟 +𝑇𝑠 [𝜇𝐼𝑠𝑞 2 + 𝜔𝑒 𝑖𝑠𝑞 (𝑘) + ] (14)
𝜏𝑟
𝜑𝑟𝑑𝑞 = (𝐿𝑚 − ) 𝑖𝑠𝑑𝑞 + 𝜑 (9)
𝐿𝑚 𝐿𝑚 𝑠𝑑𝑞
𝑖𝑠𝑞 (0) (𝑘 + 1) = 𝛾 𝑖𝑠𝑞 (𝑘)
where 𝐿𝑚 , 𝐿𝑠 , 𝐿𝑟 are the mutual, stator and rotor inductance
respectively, 𝑅𝑠 the stator resistance, 𝜑𝑟𝑑𝑞 the rotor flux in dq −𝑇𝑠 [𝜇 𝐼𝑠𝑑 2 + 𝜔𝑒 𝑖𝑠𝑑 (𝑘) + 𝜆 𝜔𝑒 ] (15)
Similarly the input reactive power is obtained from the The application times for all the vectors are then obtained
predictions of the source current according to the discrete from (18) as though each of the vector is applied for one
model of the input filter. The instantaneous reactive power can sampling resulting in a set of linear equations. Since the
be predicted by using calculations are done for one sampling period, the sum of the
application times for all vectors should be equal to the
𝑄(𝑘 + 1) = 𝑉𝑠𝛼 (𝑘 + 1)𝐼𝑠𝛽 (𝑘 + 1) sampling interval.
−𝑉𝑠𝛽 (𝑘 + 1)𝐼𝑠𝛼 (𝑘 + 1) (16)
𝑇𝑠 = 𝜏1 + 𝜏2 + 𝜏3 + 𝜏4 + 𝜏0 (19)
where the subscripts α and β represent real and imaginary
where τ1, τ 2, τ 3, τ 4 are the times for the four active vectors and
components of the associated vector.
t0 is the time for three zero vectors. Combining (18) and (19) it
is possible to obtain the normalizing coefficient ‘K’. The
C. Cost Function Calculation application times for zero vectors can be calculated from the
Similar to the MPC, the switching states are calculated based times for active vectors as shown in (20).
on a cost function minimization. The cost function can include
different performance factors according to the controlled 𝜏0 = 𝑇𝑠 − ( 𝜏1 + 𝜏2 + 𝜏3 + 𝜏4 ) (20)
variables and to the constraints required. More than one
variable can be controlled with the same control loop, The duty cycle ‘δ’ for all the vectors is then calculated using
providing a multi-objective control approach and avoiding (21) for j=0, 1, 2, 3, 4
nested loops.
The cost function for stator current control and input reactive 𝜏(𝑗)
𝛿(𝑗) = (21)
power control are computed for all the five vectors (active and 𝑇𝑠
zero) resulting in 5 cost functions 𝐺0 , 𝐺1 , 𝐺2 , 𝐺3 and 𝐺4 with
j=0, 1, 2, 3, 4. The values of the duty cycles thus obtained can be verified by
using the relation shown in (22).
𝐺(𝑗) = |𝐼𝑠𝑑𝑟𝑒𝑓 (𝑘 + 1) − 𝐼𝑠𝑑 𝑗 (𝑘 + 1)| +
𝛿 1 + 𝛿2 + 𝛿3 + 𝛿4 + 𝛿0 = 1 (22)
|𝐼𝑠𝑞𝑟𝑒𝑓 (𝑘 + 1) − 𝐼𝑠𝑞 𝑗 (𝑘 + 1)|
+𝐴 ∗ |𝑄𝑟𝑒𝑓 (𝑘 + 1) − 𝑄(𝑘 + 1)| (17) The active and zero vectors are then applied for their
respective duty cycles within one sampling interval.
In (17) 𝐼𝑠𝑑𝑟𝑒𝑓 (𝑘 + 1) and 𝐼𝑠𝑞𝑟𝑒𝑓 (𝑘 + 1) are the predicted The obtained duty cycles are directly related to the error
values of the stator current references in dq frame. 𝐼𝑠𝑞𝑟𝑒𝑓 (𝑘 + between the reference and the control variable; thereby
1) is generated by a Proportional-Integral (PI) speed achieving modulation and control functions simultaneously.
controller. For minimization of reactive power, 𝑄𝑟𝑒𝑓 (𝑘 + The intrinsic modulation strategy in the M2PC ensures fixed
1) = 0. ‘A’ is the weighting factor for the reactive power switching frequency.
control. The selection of this variable is done by trial and III. SIMULATION RESULTS
error. The weighting factor decides the controllability of one
control variable over another when more than one control The M2PC control strategy proposed in this paper is applied to
variable exists. Once the cost functions for each of the active a 30kW induction machine fed by a direct matrix converter.
and zero vectors are calculated; the predictive control will The simulations are carried out in MATLAB-Simulink
choose the best sequence of vectors which produces the lowest environment using the schematic diagram in Fig.1 The
error. These vectors are then applied within the sampling parameters for the input filter and load are given in Table.1.
interval with different duty cycles. The duty cycles are TABLE I. SYSTEM PARAMETERS
calculated based on the cost functions values of the active and Parameter Symbol Value
zero vectors. Filter inductance Lf 0.7mH
Filter capacitance(star) Cf 24.9 µF
D. Calculation of Duty Cycles Damping resister Rd 10 Ω
The expressions for the duty cycles are derived based on the Stator resistance Rs 0.081 Ω
characteristic of MPC where the vector with the least error Rotor resistance Rr 0.0788 Ω
will be applied for the longest time. Hence the application Stator leakage inductance Lls 1.36mH
time of a vector is inversely proportional to the cost function Mutual inductance Lm 36.8mH
associated with that vector. This can be expressed by (18), Rotor leakage inductance Llr 0.916mH
Rated Speed ω 1470rpm
𝐾 Rated Voltage V 380V
𝜏(𝑗) = (18)
𝐺(𝑗) Rated Torque Te 195Nm
where ‘G(j)’ is the cost function and ‘𝜏 (j)’ is the time for which Sampling time Ts 50µs
the respective vector is applied.
Torque and stator current control are achieved by the predictive of speed at full load when controlled by M2PC are shown in
control algorithm while the speed control is implemented with Fig. 5. The simulation is carried out for different values of
a PI regulator. weighting factor, A and the performance is analysed. In Fig.5
Fig.4.a) and 4.b) shows the steady state three phase stator a), there is no input reactive power control and the stator
currents at no load for both MPC and M2PC respectively. As currents, flux, and torque are controlled using the M2PC
seen from the figure, the stator currents of the motor controlled algorithm. However the input currents related to this conditons
using MPC present higher and more spread harmonic content have a considerable amount of distortion as seen from Fig.6a).
than the one controlled by M2PC with the same simulation When the input reactive power control is added to the M2PC
parameteres and conditions. algorithm with a weighting factor of 0.0002 (VA)-1, the
50 40
performance of the machine remained the same and the input
30
current quality is considerably improved as noticeable from
20 Fig.6.b).This proves that there is good control of both control
Stator Currents (A)
-20
of the system at different value of A and it is found that as the
-30 value of A increases, the importance on the reactive power
-50
0.6 0.61 0.62 0.63 0.64 0.65 -40
0.4 0.41 0.42 0.43 0.44 0.45
control becomes dominant, which in turn affects the
Time (s)
performance of the load current controller. To demonstrate this,
Time (s)
(a) (b)
1 0.5
when A=0.0006 (VA)-1 , as seen in Fig.6.c) the input currents
THD=11.8% THD=4.8%
0.8 0.4 starts gettting more distorted and the performance of the load is
also affected. This is visible from the distortions in torque and
Amplitude (A)
0.6 0.3
stator currents in Fig. 5.c). The selection of the weighting
0.4 0.2
factor A, as mentioned before is done on a trial and error basis.
0.2 0.1 Therefore the optimum value of A which gives a satisfactory
0 0
control overt the input and output of the drive has to be
0 2000 4000 6000 8000 10000 12000 14000 16000 0 0.5 1 1.5 2 2.5 3 3.5 4
Frequency (Hz) Frequency (Hz)
x 10
4
selected.
(c) (d)
Fig.4 (a) Steady state three phase no load stator currents for M2PC, (b) Steady
state three phase no load stator currents for MPC ,(c) FFT of stator current
IV. CONCLUSIONS
controlled by M2PC , (d) FFT of stator current controlled by MPC The use of M2PC to control an induction motor fed by a direct
To investigate this further, the fast fourier transform (FFT) of matrix converter is discussed in this paper. The prospect of
both the currents shown in Fig.4.c) and 4.d) is analysed. In the using more than one control variable in the cost function of the
case of the MPC controlled drive, a total harmonic distorion control is also investigated with promising simulation results.
(THD) of 11.8% is calculated with more harmonic content in By varying the weighting factor in the cost function, the
the low frequency range. However, in the case of the M2PC relevance of the control variable can be altered according to
controlled drive the THD is about 4.8%. Fig.4.d) shows the the design requirements. A comparison of quality of the stator
harmonics in sidebands multiples of switching frequency currents controlled by MPC and M2PC indicates that M2PC
(20kHz) which is the result of a fixed switching frequency can provide better current quality with fixed switching
operation. The three phase stator currents, startor flux , speed frequency.
and electromagnetic torque of the machine during the reversal
Flux (Wb) currents (A)
0 0
Stator
0 0 0
-2 -2 -2
1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.1 1.2 1.3 1.4 1.5 1.6 1.7
2000 2000 2000
(rpm)
Speed
0 0 0
200 200
0 0 0
-200 -200 -200
-400 -400 -400
1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.1 1.2 1.3 1.4 1.5 1.6 1.7
Time (s) Time (s) Time (s)
(a) (b) (c)
Fig.5: Simulation results for M2PC of IM drive: Three Phase stator currents, Stator flux in alpha-beta , speed of the motor and electromagnetic torque during
speed reversal for (a) A=0 (VA)-1, (b) A=0.0002 (VA)-1and (c) A=0.0006 (VA)-1with a load torque of 195Nm.
100 100 100
0 0 0
0.9 0.95 1 1.05 1.1 0.9 0.95 1 1.05 1.1 0.9 0.95 1 1.05 1.1
Time (s) Time (s) Time (s)
V. REFERENCES
[9] T. Geyer, G. A. Beccuti, G. Papafotiou, and M. Morari,
[1] D. Casadei, G. Grandi, G. Serra, and A. Tani, "Space "Model Predictive Direct Torque Control of permanent
vector control of matrix converters with unity input power magnet synchronous motors," in IEEE Energy Conversion
factor and sinusoidal input/output waveforms," in Fifth Congress and Exposition (ECCE), 2010, pp. 199-206.
European Conference on Power Electronics and [10] P. Correa, M. Pacas, and J. Rodriguez, "Predictive Torque
Applications, 1993, pp. 170-175 vol.7. Control for Inverter-Fed Induction Machines," IEEE
[2] A. Alesina and M. Venturini, "Analysis and design of Transactions on Industrial Electronics, vol. 54, pp. 1073-
optimum-amplitude nine-switch direct AC-AC converters," 1079, 2007.
IEEE Transactions on Power Electronics, vol. 4, pp. 101- [11] R. Vargas, U. Ammann, B. Hudoffsky, J. Rodriguez, and P.
112, 1989. Wheeler, "Predictive Torque Control of an Induction
[3] D. Casadei, G. Serra, and A. Tani, "The use of matrix Machine Fed by a Matrix Converter With Reactive Input
converters in direct torque control of induction machines," Power Control," IEEE Transactions on Power Electronics,
IEEE Transactions on Industrial Electronics, vol. 48, pp. vol. 25, pp. 1426-1438, 2010.
1057-1064, 2001. [12] H. Jiabing and Z. Q. Zhu, "Improved Voltage-Vector
[4] P. Wheeler, J. C. Clare, M. Apap, D. Lampard, S. Sequences on Dead-Beat Predictive Direct Power Control
Pickering, K. J. Bradley, et al., "An Integrated 30kW of Reversible Three-Phase Grid-Connected Voltage-Source
Matrix Converter based Induction Motor Drive," in 36th Converters," IEEE Transactions on Power Electronics,
IEEE Power Electronics Specialists Conference, PESC vol. 28, pp. 254-267, 2013.
'05, pp. 2390-2395. [13] P. Antoniewicz and M. P. Kazmierkowski, "Virtual-Flux-
[5] T. F. Podlesak, D. C. Katsis, P. W. Wheeler, J. C. Clare, L. Based Predictive Direct Power Control of AC/DC
Empringham, and M. Bland, "A 150-kVA vector-controlled Converters With Online Inductance Estimation," IEEE
matrix converter induction motor drive," IEEE Transactions on Industrial Electronics, vol. 55, pp. 4381-
Transactions on Industry Applications, vol. 41, pp. 841- 4390, 2008.
847, 2005. [14] L. Tarisciotti, P. Zanchetta, A. Watson, J. Clare, and S.
[6] E. F. Camacho and C. B. Alba, Model predictive control: Bifaretti, "Modulated Model Predictive Control for a 7-
Springer, 2013. Level Cascaded H-Bridge back-to-back Converter," IEEE
[7] J. Rodriguez and P. Cortes, Predictive control of power Transactions on Industrial Electronics, vol. PP, pp. 1-1,
converters and electrical drives vol. 37: John Wiley & 2014.
Sons, 2012. [15] L. Tarisciotti, P. Zanchetta, A. Watson, J. Clare, M.
[8] J. Rodriguez, J. Pontt, C. A. Silva, P. Correa, P. Lezana, P. Degano, and S. Bifaretti, "Modulated model predictve
Cortes, et al., "Predictive Current Control of a Voltage control (M2PC) for a 3-phase active front-end," in 2013
Source Inverter," IEEE Transactions on Industrial IEEE Energy Conversion Congress and Exposition
Electronics, vol. 54, pp. 495-503, 2007. (ECCE),pp.1062-1069.