Next Article in Journal
Grid-Scale Ternary-Pumped Thermal Electricity Storage for Flexible Operation of Nuclear Power Generation under High Penetration of Renewable Energy Sources
Next Article in Special Issue
Proportional Usage of Low-Level Actions in Model Predictive Control for Six-Phase Electric Drives
Previous Article in Journal
Bottled Biogas—An Opportunity for Clean Cooking in Ghana and Uganda
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Algorithm for Implementation of Optimal Vector Combinations in Model Predictive Current Control of Six-Phase Induction Machines

1
Laboratory of Power and Control Systems (LSPyC), Facultad de Ingeniería, Universidad Nacional de Asunción, Luque 2060, Paraguay
2
Facultad Politécnica, Universidad Nacional de Asunción, San Lorenzo 2111, Paraguay
*
Author to whom correspondence should be addressed.
Energies 2021, 14(13), 3857; https://doi.org/10.3390/en14133857
Submission received: 23 May 2021 / Revised: 10 June 2021 / Accepted: 23 June 2021 / Published: 26 June 2021
(This article belongs to the Special Issue Control, Operation and Protection of Multiphase Machines and Drives)

Abstract

:
The development of new control techniques for multiphase induction machines (IMs) has become a point of great interest to exploit the advantages of these machines compared to three-phase topology, for example, the reduced phase currents and lower harmonic contents. One of the most analyzed techniques is the model-based predictive current control (MPC) with a finite control set. This technique presents high xy currents because of the application of one switching state throughout the whole sampling period. Nevertheless, it is one of the most used due to its excellent dynamic response. To overcome the aforementioned drawbacks, new techniques called virtual vectors have been developed, but although there are several articles with experimental results, the algorithm for implementing the technique has not been appropriately described. This document provides a clear and detailed explanation for algorithm implementation of virtual vectors through two proposed variants VV4 and VV11, in a six-phase machine drive. The first entails lower computational cost and the second lower loss in the xy plane. According to performance indicators such as the total harmonic distortion and the mean square error for both case studies, experimental tests were evaluated to determine the implementation’s behaviour.

1. Introduction

For decades, three-phase IM was used for most variable speed applications, and control techniques are widely developed. Still, in recent years, the strengths of the multiphase IM has encouraged their use and thus the development of new sophisticated control techniques to regulate their speed, flux, torque or current. Multiphase IM stands out from three-phase IM due to the following characteristics: an inherent ability for post-fault operation without additional hardware, lower stator current per phase for the same voltage and lower torque ripple [1,2,3,4]. On the other hand, implementing more complex control strategies and the need for careful regulation of secondary currents that appear in multiphase IM also entails a higher computational cost [5,6]. However, this fact can be addressed through more powerful digital controllers such as Digital Signal Processor (DSP).
In this sense, concerning model-based control techniques, there are numerous options with their advantages and disadvantages. Among them is the MPC with finite set control, is one of the most widely used [7,8,9,10,11]. This technique employs a cost function to select an optimal vector from the 64 available vectors, for a six-phase IM, and apply it during the sampling period. The cost function consists of a configurable weighting factor to prioritize current tracking or loss reduction in the xy plane. This method has the disadvantage of high tracking error and considerable losses in the xy plane since it uses a single vector [12,13,14]. This fact produces high harmonic content in the stator currents compared to other approaches that integrate modulation strategies in MPC to face this weakness [15].
An alternative to mitigate the xy currents was proposed in [12], called virtual voltage vectors (VV). This method minimizes current losses by using a zero vector (ZV), a large vector (LV), and the corresponding medium-large collinear vector (MLV) at each sampling time ( T s ) since in the xy plane, these vectors are 180 degrees out of phase and tend to be eliminated. Thus, to obtain a zero voltage on average in the xy plane, it is necessary to apply the vectors in a certain proportion of time, which means the LV is applied 0.73 of T s and the MLV 0.27 of T s [16]. This technique uses less computational resources because it analyzes only 12 vectors; however, combining the vectors decreases the scope, covering only 93 % of the LV module in the α β plane.
Another MPC called M2PC based on space vector modulation (SVM), it uses a cost function to determine a pair of optimal LV at each T s . The selected optimal LV and the ZV a combination is performed to represent the target voltage more accurately, cover 96 % of the entire α β plane. This method allows to cover an area, therefore minimising the tracking error and naturally reducing the losses because the LV are mapped as the smallest in the xy plane [17].
Recently, other MPC techniques have been developed, namely as NM2PC, which uses up to 4 active vectors at each T s and does not use the ZV. Employing a cost function, it selects up to 4 vectors, two LV and two medium vectors (MV) to reduce the losses in the xy plane. It is essential to mention that in M2PC, the cost function implements a weighting factor between current tracking and xy loss reduction, and to some extent, it is expected that the algorithm tends to use the zero vector to reduce losses, even if the current tracking decreases in the α β plane, so the NM2PC presents an advantage by not using the ZV, but on the other hand, the tracking ability of the controller for small currents deteriorates [18]. Further, in [19], it was shown a direct relationship between the stability of modulated MPC and the values of sampling frequency and rotor speed for multiphase machines.
Although several works, as described above, have reported experimental implementations of MPC techniques for multiphase IM drives, few details have explained how to implement these techniques in DSP. In that sense, the main goal of this paper is to describe the implementation of optimal vectors combinations into MPC technique focused on two study situations called virtual vectors with four application times (VV4) and virtual vectors with 11 application times (VV11). In both cases, two active vectors with different duty cycles are used to effectively reduce the xy currents and improve the current tracking capability using an asymmetrical six-phase IM. Experimental results are presented to contrast both study situations, comparing in terms of current tracking considering the mean square error (MSE) and total harmonic distortion (THD) of the stator currents, to analyze the performance of the system. The proposed method is tested at different operating points, considering transient and steady-state conditions.
The sections that compose the rest of this document are the following: the mathematical model of the six-phase IM is shown in Section 2, as well as the voltage source inverter (VSI) used. Section 3 presents the description of the classic MPC, which will be used for the proposed MPC design. Section 4 describes the optimal vectors combinations implementation into MPC using two VV case scenarios. Section 5 shows the experimental results of both proposed methods. Conclusions and special remarks are presented in the last section.

2. Six-Phase IM Mathematical Model

The system is composed of an asymmetrical six-phase IM and a six-phase VSI based on isolated gate bipolar transistors (IGBT) connected as shown in the electrical diagram in Figure 1. The stator phase voltages of the six-phase IM depend on the switching states and the dc-link voltage source ( V d c ). In this regard, complementary signals are applied to each leg of the six-phase VSI, and thus it is possible to obtain 2 6 = 64 different switching states. These states are represented by a six-dimensional vector [ S ] = [ S a , S d , S b , S e , S c , S f ] . Then, the stator phase voltages can be calculated from the ideal six-phase VSI model as shown below:
v a s v d s v b s v e s v c s v f s = V d c 3 2 0 1 0 1 0 0 2 0 1 0 1 1 0 2 0 1 0 0 1 0 2 0 1 1 0 1 0 2 0 0 1 0 1 0 2 S a S d S b S e S c S f
It is possible to simplify the system model employing the vector space decomposition (VSD) technique. Making transformations through the amplitude invariant Clarke conversion matrix (C), in this way is easy to represent the six-dimensional space defined by the six phases ( a , d , b , e , c , f ) in three two-dimensional orthogonal planes in the stationary reference frame, α β , xy, and z 1 z 2 . The α β plane is related to the torque and flux production. The xy plane is linked to the copper losses, while the components projected in the z 1 z 2 plane are not examined due to the adopted topology (isolated neutral points). The transformation matrix is shown below:
C = 1 3 1 3 2 1 2 3 2 1 2 0 0 1 2 3 2 1 2 3 2 1 1 3 2 1 2 3 2 1 2 0 0 1 2 3 2 1 2 3 2 1 1 0 1 0 1 0 0 1 0 1 0 1
Hence, using the transformation matrix defined in (2) the stator phase voltages are converted into orthogonal planes through the following equation:
v α s v β s v x s v y s v z 1 s v z 2 s = C v a s v d s v b s v e s v c s v f s
It must be noted that the 64 switching states can be projected simultaneously in both planes as shown in Figure 2, where 49 of them are different. Vectors can be classified into (LV), (MLV), (MV), and (SV), which are marked as red, green, blue, and light blue circles, respectively.
Furthermore, the mathematical model of the system by using the state-space representation is expressed as follows:
X . ( t ) = A ( t ) X ( t ) + B ( t ) U ( t ) + H n p ( t )
defining X ( t ) as the state vector, formed by stator and rotor currents as shown in (5), U ( t ) as the input vector constituted by the stator voltages represented by (6), while A ( t ) and B ( t ) are the matrices that contain the electrical parameters of the six-phase IM expressed in (7) and (8), respectively. The Gaussian process noises ( n p ( t ) ), and its corresponding noise weighting matrix (H), defined as the identity matrix, have been contemplated in the model.
X ( t ) = x 1 , , x 6 T = [ i α s , i β s , i x s , i y s , i α r , i β r ] T
U ( t ) = [ u 1 , , u 4 ] T = [ u α s , u β s , u x s , u y s ] T
A ( t ) = a 11 a 12 0 0 a 15 a 16 a 21 a 22 0 0 a 25 a 26 0 0 a 33 0 0 0 0 0 0 a 44 0 0 a 51 a 52 0 0 a 55 a 56 a 61 a 62 0 0 a 65 a 66
B ( t ) = b 11 0 0 0 0 b 22 0 0 0 0 b 33 0 0 0 0 b 44 b 51 0 0 0 0 b 62 0 0
where the coefficients are defined as:
a 11 = a 22 = R s L r L s L r L m 2 a 12 = a 21 = L m 2 ω r L s L r L m 2 a 15 = a 26 = R r L m L s L r L m 2 a 16 = a 25 = L m L r ω r L s L r L m 2 a 33 = a 44 = R s L l s a 51 = a 62 = R s L m L s L r L m 2 a 52 = a 61 = L s L m ω r L s L r L m 2 a 55 = a 66 = R r L s L s L r L m 2 a 56 = a 65 = L r L s ω r L s L r L m 2 b 11 = b 22 = L r L s L r L m 2 b 33 = b 44 = 1 L l s b 51 = b 62 = L m L s L r L m 2
L s = L l s + L m , L r = L l r + L m , R s , and R r being the electrical parameters of the six-phase IM.
The output vector, including zero-mean Gaussian measurement noises ( n m ( t ) ) and uncorrelated process, in the model, is calculated as follows:
Y ( t ) = 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 X ( t ) + n m ( t )
Finally, the torque ( T e ) and the load torque ( T L ) of the six-phase IM are obtained using the following equations:
T e = 3 P ψ α s i β s ψ β s i α s
J i ω . m + B i ω m = T e T L
ω m = ω r P
J i , B i , ω m , ω r , ψ α s , ψ β s and P, being the inertia coefficient, the friction coefficient, the rotor mechanical speed, the rotor electrical speed, the stator fluxes and the number of pole pairs, respectively.

3. Classic MPC

To implement the classic MPC in a DSP, the mathematical model of the system must be in a discrete version. For this purpose, a discretization technique derived from the forward-Euler equation is used to obtain the predictive model. Note that the predicted variables only depending on past values of the variables and not on present values. As a result, the predictive model is described as:
X ^ ( k + 1 | k ) = X ( k ) + f X ( k ) , U ( k ) , ω r ( k ) , T s
being k the actual sample.

3.1. Reduced Order Observer

Rotor currents can be estimated through the reduced-order observer. There are two well-known techniques, the Luenberger Observer (LO) and the Kalman Filter (KF) [20,21], the latter being the more optimal as it takes into account the noise input to the sensors, thus optimizing the gain of the observer. Therefore, the KF is designed and implemented in this document, increasing the accuracy of the predictions. The equations of the mathematical model of system state space can be defined as:
X ^ ( k + 1 | k ) = A ( k ) X ( k ) + B ( k ) U ( k ) + H n p ( k )
Y ( k + 1 | k ) = 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 X ( k + 1 ) + n m ( k + 1 )
where k represents the actual sample, n p ( k ) the process noise, and H the noise weight matrix.
The discretized matrices A ( k ) and B ( k ) are represented in (16) and (17) respectively. A ( k ) depends on the present amount of the rotor electrical speed and it must be actualized at every sampling period.
A ( k ) = a 11 ( k ) a 12 ( k ) 0 0 a 15 ( k ) a 16 ( k ) a 21 ( k ) a 22 ( k ) 0 0 a 25 ( k ) a 26 ( k ) 0 0 a 33 ( k ) 0 0 0 0 0 0 a 44 ( k ) 0 0 a 51 ( k ) a 52 ( k ) 0 0 a 55 ( k ) a 56 ( k ) a 61 ( k ) a 62 ( k ) 0 0 a 65 ( k ) a 66 ( k )
B ( k ) = b 11 ( k ) 0 0 0 0 b 22 ( k ) 0 0 0 0 b 33 ( k ) 0 0 0 0 b 44 ( k ) b 51 ( k ) 0 0 0 0 b 62 ( k ) 0 0
The discretized matrix coefficients were included in the Appendix A.

3.2. Cost Function

The objective of the cost function is to obtain a numerical value that allows a decision to be made by prioritizing one or more variables. We can optimize, for instance, the switching losses, torque ripple, harmonic content, or tracking error of stator currents. The latter is usually the priority, so the current tracking error calculations are made in the α β and xy plane. The current tracking error between the reference ( i s ( k + 2 ) * ) and their predicted values ( i ^ s ( k + 2 ) ) in the α β and xy planes is calculated as follows:
J = ( e α s ) 2 + ( e β s ) 2 + λ x y [ ( e x s ) 2 + ( e y s ) 2 ]
where
e α s = i α s ( k + 2 ) * i ^ α s ( k + 2 ) e β s = i β s ( k + 2 ) * i ^ β s ( k + 2 ) e x s = i x s ( k + 2 ) * i ^ x s ( k + 2 ) e y s = i y s ( k + 2 ) * i ^ y s ( k + 2 )
Note that a second-step ahead prediction is implemented to consider the delay compensation. This fact is produced by the comparable time between the computation of the control signal and the sampling period [22]. The tuning parameter ( λ x y ) expressed in (18) is usually considered in the control of the multiphase machines to provide greater weight to the α β currents over the xy currents.

4. Proposed Algorithm

This section aims to describe the algorithm implementation in a DSP based on Optimal Vector Combinations characteristics into MPC described in [12,13,16].

4.1. Case 1: VV4

The first MPC strategy to be discussed in this paper is called VV4. The described implementation algorithm is based on the virtual vector technique [12]. It uses two vectors per sampling time, a long vector and the nearest collinear medium-large vector, to obtain the target voltage. The Figure 3 exemplifies the selection of two vectors, LV and MLV, from the I-sector of the α β plane. The advantage is taken of the fact that the LV and collinear MLV are 180 degrees out of phase in the xy plane, and although they have different modulus as shown in Figure 2, they tend to subtract each other, thus decreasing the losses in this plane.
The VV4 technique adopts strategies that take into account the ratio of the modulus of LV and MLV to calculate their application times to obtain a zero average value in the xy plane. Suitably, it is required to apply 73 % of the LV time and 27 % of the MLV time as described in [12]. The VV4 technique proposes a simple algorithm, based mainly on an interrupt and four conditional sequences, which divide the time of total application in four intervals, for which it is possible to apply 75 % of the time LV and 25 % MLV. This virtual vector is applied as follows:
V V 4 = d l L V + d m l M L V
where d l = 0.75 T s v and d m l = 0.25 T s v .
On the other hand, in VV4, there is a decrease in the maximum final voltage that can be reached since it conditions the use of the LV and MLV vectors, obtaining a resulting vector of lower modulus.
The process for the implementation of VV4 is shown in Figure 4 and detailed below:
  • Step 1: The algorithm starts with the variables initialization and calculation of the electrical parameters of IM.
  • Step 2: The basic control technique, in this case the classical MPC, acquires measurements of the stator currents every sampling time ( T s ). To implement the VV4 technique, the original sampling time ( T s ) is increased by 4, hereafter referred to as ( T s v ), consequently the sampling frequency is reduced by 4. On the other hand, an interrupt is used to control current readings, calculations and transforms, vector selection and application. The input period to the routines contained in this interrupt is every T s v /4, and they are called interrupt time ( T i ), so the interrupt is entered four times per sampling time.
  • Step 3: In the first interval ( T i ) of the sampling time ( i = 1 ), the measurements are considered to proceed to perform the transformations and calculations.
  • Step 4: The 12 external vectors (LV) are evaluated using a cost function and the optimal vector is obtained. The MLV vector corresponding to the optimal LV is identified, and both will be applied during the sampling time.
  • Step 5: Only one vector is applied at each interrupt time, making use of conditional statements and a counter to identify the four interrupts in a sampling period. To reduce the losses in the xy plane, the total application time of the active vectors is distributed between the LV and the MLV in a proportion of 75 % and 25 % respectively, so that the LV is applied three times and the MLV once. The final pattern of applied vectors is shown in Figure 5.

4.2. Case 2: VV11

The second case study is the technique called VV11, based on the MPC strategy with virtual vectors. Unlike VV4, in VV11 the LVand MLV are applied 73 % and 27 % of the time, so that the average loss in the xy plane tends to zero. In order to achieve the latter, the sampling period is divided into 11 equal intervals, using an interrupt, conditional sequences, and a counter (i), so that T i = T s v / 11 .
Similar to VV4, measurements are taken in the first interrupt of each sampling period ( i = 1 ), calculations, transformations are performed and optimal vectors are selected. Only one of them its applied in each T i , in total, eight times the LV and three times the MLV are applied, resulting in approximately 72.73 % of T s v for LV and 27.27 % of T s v for MLV. The final applied vector pattern for VV11 is shown in Figure 6. This extra process implies higher computational costs compared to VV4, but it is closer to the appropriate proportion of application of the vectors LV and MLV.

5. Experimental Results

The case studies are implemented to analyze their performance by comparing them with each other through the experimental results collected in the six-phase IM bench.

5.1. Six-Phase IM Test Bench Description

The control system is based on Six-Phase VSI, driven by a DSP (TMS320F28335 from Texas Instruments), with a Technosoft MCK28335 development environment. The IM parameters are shown in Table 1 determined by the tests described in [23,24].
The periodic measurements of the rotor currents and speed, necessary to perform the calculations and implement the control techniques, are made by a series of sensors. A 10,000 ppr encoder is used to obtain the mechanical speed of the motor.
A matching board is provided to amplify the signals coming from 4 LA 55-P sensors used to measure the stator currents. The digitalization is performed by means of a 16-bit A/D converter of the DSP. A dc-link is set to 325 V through a three-phase rectifier and the power grid with a step-down transformer. A mechanical brake pad was installed to simulate mechanical loading of approximately 120 Nm. A block diagram of the test bench is shown in Figure 7.

5.2. Figures of Merit

First, we tested experimentally VV4 and VV11 to analyze the performance under steady-state and transient conditions. The MSE between the reference and stator current measured in the α β and xy planes is obtained to compare current and loss tracking. On the other hand, the THD in the α β plane is calculated to quantify the harmonic distortion obtained with both techniques. In this context, the MSE is defined as:
MSE ( i s Φ ) = 1 N k = 1 N ( i s Φ [ k ] i s Φ * [ k ] ) 2
where N is the number of samples, i s Φ * the stator current reference, i s Φ the measured stator currents and Φ   { α , β , x , y } . At the same time, THD is defined as:
THD ( i s ) = 1 i s 1 2 j = 2 N ( i s j ) 2
where i s 1 is the fundamental stator currents and i s j is the harmonic stator currents.

5.3. Steady-State Study

For the purpose of analyzing the performance of the xy current, references are set to zero ( i x s * = i y s * = 0 ), and the d axis current ( i d s * = 1.5 ) has been considered for all cases.
For the VV4 and V11 technique, the sampling frequency is 2.5 kHz. The tests are performed at mechanical speeds of 100 rpm to 600 rpm. The results obtained from VV4 and VV11 at the different speeds are shown in Table 2 and Table 3. At the same time, Figure 8 and Figure 9 show the stator currents tracking in α β and xy planes for classic MPC, VV4 and VV11 respectively. The operations were tested at a sampling frequency of 2.5 kHz and a rotor speed of 200 rpm for the techniques.
The results show good current tracking in the α β plane for both techniques, with a slight advantage for VV11, showing lower MSE values for all speeds. We can observe higher values for VV4 in terms of MSE in the xy plane, corroborating the decrease in the resultant vector in the xy plane also for all speeds. In terms of THD, it is observed that the VV11 technique leads to higher harmonic distortion.

5.4. Transient Study

A step modification in the rotor speed is considered for a transient operation from 200 to 200 rpm (reversal condition) as shown in Figure 10. Additionally, Figure 11 presents a dynamic test (q current tracking), which reveals the transient operation of VV4 (a) and VV11 (b) for a step-change in the q current ( i q s * ). The dynamic performance is obtained through a reversal condition of the rotor speed ( ω m * ). The reaching time is approximately 0.5 ms of both techniques, and the overshoot, for VV4 and VV11, are 10% and 12.5 %, respectively. Therefore, it can be noted that both control strategies show good speed tracking in reversal conditions; however, in terms of i q s current tracking, VV4 is slightly better than VV11.

6. Conclusions

The paper describes the implementation of two virtual vector modulation techniques applied to a 15 kW six-phase asymmetrical IM. Among the advantages of both techniques based on the same control strategy are: (a) reduced amount of switching, since only the change of state of the transistor pair of one VSI branch is required to switch from the MLV to the collinear LV; (b) combining LV and their corresponding MLV reduces the losses in the xy plane. It can be seen in this work that as the combination of LV and MLV is closer to a ratio of 73 % and 27 % at each sampling time, the xy currents and the MSE in the α β plane are reduced, although these improvements are negligible compared to the 75 % and 25 % ratio used in VV4. On the other hand, the VV11, presents higher harmonic distortion and higher computational cost. These observations lead us to conclude that VV11 does not present a higher performance over VV4, so the best alternative in order to achieve good figures of merit and lower computational cost is VV4. It can be noted that due to the simple nature of the implemented controllers, they could be performed in any industrial application of a wide range of power and speed, being the computational cost the main reason to prefer VV4 over VV11.

Author Contributions

Conceptualization, C.R., L.D., O.G. and M.A.; methodology, C.R., L.D., O.G. and M.A. software, C.R., L.D., O.G. and M.A.; validation, C.R., L.D., O.G. and M.A.; formal analysis, C.R., L.D., O.G., M.A., J.R. and R.G.; investigation, C.R., L.D., O.G. and M.A.; resources, C.R., L.D., O.G., M.A., J.R. and R.G.; data curation, C.R., L.D., O.G. and M.A.; writing, original draft preparation, C.R., L.D., O.G. and M.A.; writing, review and editing, C.R., L.D., O.G., M.A., J.R. and R.G.; visualization, C.R., L.D., O.G., M.A., J.R. and R.G.; project administration, C.R., L.D., O.G., M.A., J.R. and R.G.; funding acquisition, R.G. and J.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been funded through the Consejo Nacional de Ciencia y Tecnología (CONACYT)-Paraguay, Grant Number POSG16-05.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DSPDigital Signal Processor
IGBTIsolated Gate Bipolar Transistors
IMInduction Machine
VSIVoltage Source Inverter
PIProportional-Integral
FOCField Oriented Control
MPCModel Predictive Control
PCCPredictive Current Control
PFSCCSPredictive-Fixed Switching Current Control Strategy
LVLarge Vector
MLVMedium Large Vector
MVMedium Vector
SVSmall Vector
SVMSpace Vector Modulation
VSDVector Space Decomposition
VVVirtual Vector
VV4Virtual Vectors with 4 application times
VV11Virtual Vectors with 11 application times
ZVZero Vector
LOLuenberger Observer
KFKalman Filter
MSEMean Squared Error
THDTotal Harmonic Distortion

Appendix A

a 11 ( k ) = a 22 ( k ) = 1 R s L r T s L s L r L m 2 a 12 ( k ) = a 21 ( k ) = L m 2 ω r ( k ) T s L s L r L m 2 a 15 ( k ) = a 26 ( k ) = R r L m T s L s L r L m 2 a 16 ( k ) = a 25 ( k ) = L m L r ω r ( k ) T s L s L r L m 2 a 33 ( k ) = a 44 ( k ) = 1 R s T s L l s a 51 ( k ) = a 62 ( k ) = R s L m T s L s L r L m 2 a 52 ( k ) = a 61 ( k ) = L s L m ω r ( k ) T s L s L r L m 2 a 55 ( k ) = a 66 ( k ) = 1 R r L s T s L s L r L m 2 a 56 ( k ) = a 65 ( k ) = L r L s ω r ( k ) T s L s L r L m 2 b 11 ( k ) = b 22 ( k ) = L r T s L s L r L m 2 b 33 ( k ) = b 44 ( k ) = T s L l s b 51 ( k ) = b 62 ( k ) = L m T s L s L r L m 2

References

  1. Levi, E.; Bojoi, R.; Profumo, F.; Toliyat, H.; Williamson, S. Multiphase induction motor drives—A technology status review. IET Electr. Power Appl. 2007, 1, 489–516. [Google Scholar] [CrossRef] [Green Version]
  2. Barrero, F.; Duran, M.J. Recent advances in the design, modeling, and control of multiphase machines—Part I. IEEE Trans. Ind. Electron. 2015, 63, 449–458. [Google Scholar] [CrossRef]
  3. Duran, M.J.; Barrero, F. Recent advances in the design, modeling, and control of multiphase machines—Part II. IEEE Trans. Ind. Electron. 2015, 63, 459–468. [Google Scholar] [CrossRef]
  4. Bermudez, M.; Gonzalez, I.; Barrero, F.; Guzman, H.; Kestelyn, X.; Duran, M.J. An experimental assessment of open-phase fault-tolerant virtual-vector-based direct torque control in five-phase induction motor drives. IEEE Trans. Power Electron. 2018, 33, 2774–2784. [Google Scholar] [CrossRef] [Green Version]
  5. Kali, Y.; Ayala, M.; Rodas, J.; Saad, M.; Doval-Gandoy, J.; Gregor, R.; Benjelloun, K. Current control of a six-phase induction machine drive based on discrete-time sliding mode with time delay estimation. Energies 2019, 12, 170. [Google Scholar] [CrossRef] [Green Version]
  6. Kali, Y.; Ayala, M.; Rodas, J.; Saad, M.; Doval-Gandoy, J.; Gregor, R.; Benjelloun, K. Time Delay Estimation based Discrete-Time Super-Twisting Current Control for a Six-Phase Induction Motor. IEEE Trans. Power Electron. 2020. [Google Scholar] [CrossRef]
  7. Lim, C.S.; Levi, E.; Jones, M.; Rahim, N.A.; Hew, W.P. A comparative study of synchronous current control schemes based on FCS-MPC and PI-PWM for a two-motor three-phase drive. IEEE Trans. Ind. Electron. 2013, 61, 3867–3878. [Google Scholar] [CrossRef] [Green Version]
  8. Lim, C.S.; Levi, E.; Jones, M.; Rahim, N.A.; Hew, W.P. FCS-MPC-based current control of a five-phase induction motor and its comparison with PI-PWM control. IEEE Tran. Ind. Electron. 2013, 61, 149–163. [Google Scholar] [CrossRef]
  9. Levi, E.; Barrero, F.; Duran, M.J. Multiphase machines and drives-Revisited. IEEE Trans. Ind. Electron. 2015, 63, 429–432. [Google Scholar] [CrossRef] [Green Version]
  10. Gonçalves, P.; Cruz, S.; Mendes, A. Finite control set model predictive control of six-phase asymmetrical machines—An overview. Energies 2019, 12, 4693. [Google Scholar] [CrossRef] [Green Version]
  11. Pérez-Estévez, D.; Doval-Gandoy, J. A Model Predictive Current Controller With Improved Robustness Against Measurement Noise and Plant Model Variations. IEEE Open J. Ind. Appl. 2021, 2, 131–142. [Google Scholar] [CrossRef]
  12. Gonzalez-Prieto, I.; Duran, M.J.; Aciego, J.J.; Martin, C.; Barrero, F. Model predictive control of six-phase induction motor drives using virtual voltage vectors. IEEE Trans. Ind. Electron. 2017, 65, 27–37. [Google Scholar] [CrossRef]
  13. Aciego, J.J.; Prieto, I.G.; Duran, M.J. Model predictive control of six-phase induction motor drives using two virtual voltage vectors. IEEE J. Emerg. Sel. Top. Power Electron. 2018, 7, 321–330. [Google Scholar] [CrossRef]
  14. Aciego, J.J.; Gonzalez-Prieto, I.; Duran, M.J.; Bermudez, M.; Salas-Biedma, P. Model predictive control based on dynamic voltage vectors for six-phase induction machines. IEEE J. Emerg. Sel. Top. Power Electron. 2020. [Google Scholar] [CrossRef]
  15. González, O.; Ayala, M.; Romero, C.; Rodas, J.; Gregor, R.; Delorme, L.; González-Prieto, I.; Durán, M.J.; Rivera, M. Comparative Assessment of Model Predictive Current Control Strategies applied to Six-Phase Induction Machines. In Proceedings of the 2020 IEEE International Conference on Industrial Technology (ICIT), Buenos Aires, Argentina, 26–28 February 2020; pp. 1037–1043. [Google Scholar]
  16. Durán, M.J.; Gonzalez-Prieto, I.; Gonzalez-Prieto, A. Large virtual voltage vectors for direct controllers in six-phase electric drives. Int. J. Electr. Power Energy Syst. 2021, 125, 106425. [Google Scholar] [CrossRef]
  17. Gonzalez, O.; Ayala, M.; Doval-Gandoy, J.; Rodas, J.; Gregor, R.; Rivera, M. Predictive-Fixed Switching Current Control Strategy Applied to Six-Phase Induction Machine. Energies 2019, 12, 2294. [Google Scholar] [CrossRef] [Green Version]
  18. Ayala, M.; Doval-Gandoy, J.; Rodas, J.; Gonzalez, O.; Gregor, R.; Rivera, M. A Novel Modulated Model Predictive Control Applied to Six-Phase Induction Motor Drives. IEEE Trans. Ind. Electron. 2020. [Google Scholar] [CrossRef]
  19. Ayala, M.; Doval-Gandoy, J.; Gonzalez, O.; Rodas, J.; Gregor, R.; Rivera, M. Experimental Stability Study of Modulated Model Predictive Current Controllers Applied to Six-Phase Induction Motor Drives. IEEE Trans. Power Electron 2021, 1. [Google Scholar] [CrossRef]
  20. Martin, C.; Arahal, M.R.; Barrero, F.; Durán, M.J. Five-phase induction motor rotor current observer for finite control set model predictive control of stator current. IEEE Trans. Ind. Electron. 2016, 63, 4527–4538. [Google Scholar] [CrossRef]
  21. Martín, C.; Arahal, M.R.; Barrero, F.; Durán, M.J. Multiphase rotor current observers for current predictive control: A five-phase case study. Control Eng. Pract. 2016, 49, 101–111. [Google Scholar] [CrossRef]
  22. Rodriguez, J.; Cortes, P. Predictive Control of Power Converters and Electrical Drives; John Wiley & Sons: Hoboken, NJ, USA, 2012; Volume 40. [Google Scholar]
  23. Yepes, A.G.; Riveros, J.A.; Doval-Gandoy, J.; Barrero, F.; López, O.; Bogado, B.; Jones, M.; Levi, E. Parameter identification of multiphase induction machines with distributed windings Part 1: Sinusoidal excitation methods. IEEE Trans. Energy Conv. 2012, 27, 1056–1066. [Google Scholar] [CrossRef]
  24. Riveros, J.A.; Yepes, A.G.; Barrero, F.; Doval-Gandoy, J.; Bogado, B.; Lopez, O.; Jones, M.; Levi, E. Parameter identification of multiphase induction machines with distributed windings Part 2: Time-domain techniques. IEEE Trans. Energy Conv. 2012, 27, 1067–1077. [Google Scholar] [CrossRef]
Figure 1. Six-phase IM fed by two level six-phase VSI.
Figure 1. Six-phase IM fed by two level six-phase VSI.
Energies 14 03857 g001
Figure 2. Voltage vectors in the α β and xy planes for a six-phase IM drive.
Figure 2. Voltage vectors in the α β and xy planes for a six-phase IM drive.
Energies 14 03857 g002
Figure 3. Space sectors in the α β plane for MPC with VV.
Figure 3. Space sectors in the α β plane for MPC with VV.
Energies 14 03857 g003
Figure 4. Algorithm diagram for VV4 technique.
Figure 4. Algorithm diagram for VV4 technique.
Energies 14 03857 g004
Figure 5. Timing diagram for VV4 technique.
Figure 5. Timing diagram for VV4 technique.
Energies 14 03857 g005
Figure 6. Timing diagram for VV11 technique.
Figure 6. Timing diagram for VV11 technique.
Energies 14 03857 g006
Figure 7. Test bench composed of the 2L-6PVSI, DSP unit control, six-phase IM and the mechanical load.
Figure 7. Test bench composed of the 2L-6PVSI, DSP unit control, six-phase IM and the mechanical load.
Energies 14 03857 g007
Figure 8. Performance of stator currents in α β plane for a rotor speed of 200 (rpm) at 2.5 (kHz) of sampling frequency: (a) classic MPC; (b) VV4; (c) VV11.
Figure 8. Performance of stator currents in α β plane for a rotor speed of 200 (rpm) at 2.5 (kHz) of sampling frequency: (a) classic MPC; (b) VV4; (c) VV11.
Energies 14 03857 g008
Figure 9. Performance of stator currents in xy plane for a rotor speed of 200 (rpm) at 2.5 (kHz) of sampling frequency: (a) classic MPC; (b) VV4; (c) VV11.
Figure 9. Performance of stator currents in xy plane for a rotor speed of 200 (rpm) at 2.5 (kHz) of sampling frequency: (a) classic MPC; (b) VV4; (c) VV11.
Energies 14 03857 g009
Figure 10. Reversal condition from 200 (rpm) to 200 (rpm) at 2.5 (kHz) of sampling frequency: (a) VV4; (b) VV11.
Figure 10. Reversal condition from 200 (rpm) to 200 (rpm) at 2.5 (kHz) of sampling frequency: (a) VV4; (b) VV11.
Energies 14 03857 g010
Figure 11. Transient response in q-axis of stator current for a rotor speed change from 200 (rpm) to 200 (rpm) at 2.5 (kHz) of sampling frequency: (a) VV4; (b) VV11.
Figure 11. Transient response in q-axis of stator current for a rotor speed change from 200 (rpm) to 200 (rpm) at 2.5 (kHz) of sampling frequency: (a) VV4; (b) VV11.
Energies 14 03857 g011
Table 1. Electrical and mechanical parameters of the six-phase IM.
Table 1. Electrical and mechanical parameters of the six-phase IM.
ParameterValueParameterValue
R r 0.63  Ω L s 206.2 mH
R s 0.62  Ω P3
L l s 6.4 mH P w 15 kW
L l r 3.5 mH J i 0.27 kg·m 2
L m 66.6 mH B i 0.012 kg·m 2 /s
L r 203.3 mH ω r n o m 1000 rpm
Table 2. Performance indicators for the drive operating from 100 [rpm] to 600 [rpm] for VV4.
Table 2. Performance indicators for the drive operating from 100 [rpm] to 600 [rpm] for VV4.
VV4       at        f s  = 2.5 (kHz)
ω m * MSE α MSE β MSE x MSE y THD α
1001.951.652.973.2914.22
2002.151.943.103.2615.43
3002.522.443.053.2813.64
4002.872.833.103.3516.13
5003.273.153.143.3718.04
6003.563.673.163.3816.67
Table 3. Performance indicators for the drive operating from 100 [rpm] to 600 [rpm] for VV11.
Table 3. Performance indicators for the drive operating from 100 [rpm] to 600 [rpm] for VV11.
VV11       at        f s  = 2.5 (kHz)
ω m * MSE α MSE β MSE x MSE y THD α
1001.581.262.552.8419.59
2001.531.272.502.7722.22
3001.651.382.472.7923.43
4001.631.392.562.7324.73
5002.381.862.502.7724.56
6002.602.172.572.8629.01
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Romero, C.; Delorme, L.; Gonzalez, O.; Ayala, M.; Rodas, J.; Gregor, R. Algorithm for Implementation of Optimal Vector Combinations in Model Predictive Current Control of Six-Phase Induction Machines. Energies 2021, 14, 3857. https://doi.org/10.3390/en14133857

AMA Style

Romero C, Delorme L, Gonzalez O, Ayala M, Rodas J, Gregor R. Algorithm for Implementation of Optimal Vector Combinations in Model Predictive Current Control of Six-Phase Induction Machines. Energies. 2021; 14(13):3857. https://doi.org/10.3390/en14133857

Chicago/Turabian Style

Romero, Carlos, Larizza Delorme, Osvaldo Gonzalez, Magno Ayala, Jorge Rodas, and Raul Gregor. 2021. "Algorithm for Implementation of Optimal Vector Combinations in Model Predictive Current Control of Six-Phase Induction Machines" Energies 14, no. 13: 3857. https://doi.org/10.3390/en14133857

APA Style

Romero, C., Delorme, L., Gonzalez, O., Ayala, M., Rodas, J., & Gregor, R. (2021). Algorithm for Implementation of Optimal Vector Combinations in Model Predictive Current Control of Six-Phase Induction Machines. Energies, 14(13), 3857. https://doi.org/10.3390/en14133857

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