Next Article in Journal
Heat Transfer Mechanism Investigation of Bubble Growth on the Superhydrophilic Nano-Structured Surface Using Moving Particle Semi-Implicit Method
Next Article in Special Issue
A Quantitative Stability Evaluation Method of VSC for Optimizing Control Parameters and Extending Stability Boundary to Avoid Oscillations
Previous Article in Journal
Numerical Simulation of Ice Crystal Trajectory and Its Influencing Factors Based on Lagrangian Method
Previous Article in Special Issue
Calculation of Main Circuit Steady-State Parameters for Capacitor Commutated Converter System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Online Adaptive Set of Virtual Voltage Vectors for Stator Current Regulation of a Six-Phase Induction Machine Using Finite State Model Predictive Controllers

by
Manuel R. Arahal
1,*,†,
Manuel G. Satué
1,†,
Federico Barrero
2,† and
Cristina Martín
3,†
1
Departamento de Ingeniería de Sistemas y Automática, University of Seville, Camino de los Descubrimientos S/N, 41092 Seville, Spain
2
Departamento de Ingeniería Electrónica, University of Seville, Camino de los Descubrimientos S/N, 41092 Seville, Spain
3
Departmento de Ingeniería Eléctrica, University of Seville, Camino de los Descubrimientos S/N, 41092 Seville, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Appl. Sci. 2023, 13(7), 4113; https://doi.org/10.3390/app13074113
Submission received: 17 February 2023 / Revised: 16 March 2023 / Accepted: 17 March 2023 / Published: 23 March 2023
(This article belongs to the Special Issue Electric Power Applications II)

Abstract

:
Virtual voltage vectors (VVV) have been used for the control of multi-phase induction machines, where different sub-spaces appear related to the torque production and losses generation. In the literature, several sets of VVV have been used, aiming at reducing harmonic content while maintaining a low computational burden. This paper proposes the use of different sets of VVV to regulate the stator current of multi-phase drives using finite-state model predictive controllers. In the proposal, only one set is active at each control period. This active set is obtained through a preliminary analysis using performance maps. As a result, a method is derived for the online selection using the current operating point. The selection is based on a simple computation from variables usually measured on variable-speed drives. Results are provided for a symmetrical six-phase IM, showing that the proposal improves the closed-loop performance of the multi-phase drive with a low computational cost.

1. Introduction

Model-based predictive control (MBPC) has recently been favored by the research community as a versatile control technique for power converters. Its flexibility to consider different control criteria and constraints has spurred a large number of publications. Control strategies for the tracking of torque [1], and speed [2] have been implemented. In the realm of multi-phase drives, Finite State Model Predictive Stator Current Control (FSMPSCC) is arguably the most popular technique [3]. It easily allows treating the extra number of phases considering torque production, losses and other factors pertinent for different applications [4].
In most published reports, the design of FSMPSCC faces a trade-off between conflicting criteria [5,6,7]. This fact comes from the use of a cost function (CF) with different terms. The instantaneous (discrete-time-wise) minimization of the cost function imprints in the drive of certain global behavior. However, figures of merit are related to each other. This means that using only CF tuning, it is not possible to reduce one figure of merit without increasing the others [8]. The usual practice in MBPC for drives is to select a CF structure and then tune its parameters. The goal of the tuning phase is a compromise solution that is considered best globally [9,10].
One problem related to multi-phase drives is that of harmonic content [11]. In this context, Virtual Voltage Vectors (VVV) have been proposed as a solution. The rationale of VVV is that harmonic content can be reduced by combining Voltage Vectors (VV) within a single sampling period. The basic VVV scheme combines two Voltage Source Inverter (VSI) configurations sharing the sampling period T s . The combination is such that the average x y voltage is zero. For example, a Medium-Large VV is applied during T 1 seconds and a Large VV is applied for T 2 seconds. In 5-phase and 6-phase systems T 1 = 0.73 × T s and T 2 = 0.27 × T s as explained in [12].
Many variants of this concept have been presented in the multi-phase drives’ literature. Among these, the following are useful to put the contribution into context. In [13], a nine-phase induction machine (IM) is considered. Several strategies are proposed using two and four VVs. In [14], VVV are combined with modulated predictive current control considering space vector modulation. A bi-subspace predictive current control based on VVV is presented in [15], where two groups of VVV are considered to regulate the fundamental and harmonic components, respectively. In this way, the x y components due to non-linearities are regulated. Furthermore, different VVV sets are considered for the bi-subspace scheme. In [16], VVV of two different magnitudes are used with a lookup table to exclude some combinations. CF tuning is avoided by transforming the torque and flux amplitude into an equivalent reference flux vector. The work in [17] utilizes VVV that are synthesized using a transformation method to consider both flux and torque. In [18], the increase in the predicted current caused by all VVs is divided into a combination of a group of basic VVs. In this way, the global optimization is converted into the optimization of just the basic voltage vector. In [19], a group of four candidate VVs is generated based on a lookup table. The selection allows for just one commutation. The method is applied to an asymmetrical six-phase IM. Finally, additional auxiliary vectors of two different magnitudes are used in [20] to reduce the torque ripple.
Regarding CF tuning, the usual practice is the use of Weighting Factors (WF) to give more relevance to some terms over others. Early publications use simulations and/or experiments to guide trial-and-error procedures for WF tuning, but more recently a Pareto analysis combined with a non-linear optimization method has been proposed in [21], and similar approaches have been presented in [8,22], to tackle the complexity of WF and CF tuning in conventional (not VVV-based) FSMPSCC for six-phase IM drives.
In the previously cited literature, several figures of merit are considered. In this work, the objective is to minimize current tracking errors by keeping the switching frequency below limits. The criterion is relevant, as it leads to reduced losses without compromising dynamic performance. For a thorough assessment, the performance map analysis is carried out for a six-phase motor drive using FSMPSCC and considering all operating points, where the typical Medium-Large + Large VVV set is contemplated in conjunction with the Small + Medium-Large VVV set. The main contribution of our work is the proposal of a novel VVV structure that combines two VVV sets, where the optimal tuning for each VVV set is obtained for every operating point following the previous work in [22].
In order to manage the continuous operating space, it is discretized by taking 12 sub-divisions of the load and 12 subdivisions of speed. Other choices are possible but this one gives a good trade-off between improved tuning and increased tuning effort. The proposal is then experimentally analyzed on a six-phase IM test rig. The analysis considers all operating points and possible tunings. This contribution is rarely encountered in the realm of FSMPSCC. It allows for a thorough assessment and also helps illustrate the obtained benefits.
The manuscript is organized as follows. The next section presents the scheme of FSMPSCC for a six-phase IM using VVV and introduces the experimental setup utilized for the evaluation and comparisons. Section 3 details the performance maps for the two sets of VVV, while the proposed control structure, which combines both VVV sets, is introduced and experimentally analyzed in Section 4. The conclusions are detailed in Section 5.

2. FSMPSCC of a Six-Phase Induction Drive Using VVV

Figure 1 shows the FSMPSCC diagram of a VSI-driven 6 phase IM using a set of VVV. An outer loop, indicated as the PI controller in the diagram, is responsible for torque/speed regulation. The reference currents in the d q plane are translated into the α β subspace using the Park transformation given by the matrix D as follows
D = cos θ s sin θ s sin θ s cos θ s ,
where the flux position θ a is obtained as θ a = ω e d t . Flux and electrical torque are calculated using reference currents i s d * and i s q * . These quantities are set by the outer speed loop using typically a PI controller, although different control methods, such as MBPC, can be used [23]. The electrical frequency can be found by considering the mechanical rotor speed and the slip frequency from ω e = ω r + ω s l , where field orientation is assumed and ω s l = R r L r i s q * i s d * .
The measured phase currents i s = ( i s a , i s b , i s c , i s d , i s e , i s f ) are converted to α β using the Clarke transformation matrix T based on power-invariant decoupling and defined as
T = 1 3 1 1 / 2 1 / 2 3 / 2 3 / 2 0 0 3 / 2 3 / 2 1 / 2 1 / 2 1 .
At each discrete time k, the reference signal for the tracking of the stator current i * ( k ) uses an amplitude I s * calculated as I s * = i s d * 2 + i s q * 2 . As a result, the references for stator currents are i s α * ( t ) = I s * sin ω e t , i s β * ( t ) = I s * cos ω e t , i s x * ( t ) = 0 , i s y * ( t ) = 0 . In this way, the FSMPSCC receives as input the α β reference currents i α β * , the actual α β currents i α β and the sensed mechanical speed ω m . With this information, the FSMPSCC obtains the optimal VVV to be applied during the next sampling period V V V o p t ( k + 1 ) . This VVV is achieved by the appropriate selection of the VSI states S for their respective application times T 1 , T 2 over the duration of the sampling period k + 1 . As a result, the stator currents of the IM, i s , evolve so that the α β components move to their reference values i α β * ( k + 2 ) . The selection of the optimal VVV is made by exhaustive optimization over all possible control moves (elements of the VVV set). The procedure is repeated at each sampling period following the receding-horizon rule typical of predictive controllers.
In the present case, the multi-phase drive incorporates an asymmetrical six-phase IM. The machine has distributed windings where two independent and isolated neutral points are created. Two three-phase two-level VSIs are connected to a single DC link. The switching state of every VSI leg is associated with the can be defined using a binary variable S i j , being S i j = 0 if the lower switch is ON and the upper switch is OFF, and S i j = 1 if the opposite situation occurs. In Figure 2, voltages are represented in α β and x y sub-spaces. To obtain those voltages aligned with the different phases. A decimal representation is used to identify the stator voltage in the α β plane. Such a number is derived from the binary number corresponding to the switching state S, where S = ( S a 1 , S b 1 , S c 1 , S a 2 , S b 2 , S c 2 ) , where the six phases have been numbered as a 1 , a 2 , b 1 , b 2 , c 1 , and c 2 .
The predictive model needed by the FSMPSCC is a discrete-time model for the evolution of i s derived from the IM equations in the phase variables and from the VSI structure. The model uses the vector space decomposition approach to map the 6-phase space into sub-spaces α β (energy conversion) and x y (non-torque producing). The predictive model often is presented in state space. The state space vector i = i α , i β is obtained from measurements of phase stator currents. Control actions for the actual and next sampling periods u ( k ) , u ( k + 1 ) are the voltages corresponding to the VVV.
i ^ ( k + 2 | k ) = A i ( k ) + B 1 u ( k ) + B 2 u ( k + 1 ) + G ( k )
where matrix A and vectors B 1 , B 2 are obtained as A = ( I + T s A c ) 2 , B 1 = ( I + T s A c ) T s B c , B 2 = T s B c . Matrices A c and B c are defined as
A c = a s 2 a m 4 0 0 a r 4 a l 4 a m 4 a s 2 0 0 a l 4 a r 4 0 0 a s 3 0 0 0 0 0 0 a s 3 0 0 a s 4 a m 5 0 0 a r 5 a l 5 a m 5 a s 4 0 0 a l 5 a r 5 , B c = c 2 0 0 0 0 c 2 0 0 0 0 c 3 0 0 0 0 c 3 c 4 0 0 0 0 c 4 0 0
The coefficients used are: c 1 = L s L r L m 2 , c 2 = L r / c 1 , c 3 = 1 / L l s , c 4 = L m / c 1 , c 5 = L s c 1 , a s 2 = R s c 2 , a s 3 = R s c 3 , a s 4 = R s c 4 , a r 4 = R r c 4 , a r 5 = R r c 5 , a l 4 = L r c 4 ω r , a l 5 = L r c 5 ω r , a m 4 = M c 4 ω r and a m 5 = M c 5 ω r , ϑ = 2 π / 6 .
Vector G accounts for the dynamics resulting from rotor currents that are usually not measured [24]. Note that alternative methods can be found in the scientific literature [25].
The selection of u ( k + 1 ) in discrete time k is made by minimizing CF, which can incorporate a number of different terms to convert a multi-objective optimization problem into a one with only one objective [26]. The simplest CF penalizes the predicted control error e ^ = ( i * i ^ ) , where i * represents the reference for the state space vector i, and i ^ the prediction two steps ahead. More complex objective functions can be used; in particular, penalization of VSI commutations is achieved by computing the number of the switch changes S C produced in the VSI after configuration u ( k ) is changed to u ( k + 1 ) as
S C ( k + 1 ) = i = 1 6 | S i 2 ( k + 1 ) S i 1 ( k + 1 ) | + | S i 1 ( k + 1 ) S i 2 ( k ) | + | S i 2 ( k ) S i 1 ( k ) |
where S is the state of the VSI switches. It must be taken into account that during a sampling period, there are different VSI states corresponding to the VVV in use. These are denoted as S 1 for the first sub-period and S 2 for the second sub-period. Recall from the introduction that each sub-period uses a fraction of the sampling period. In the present case, the lengths of each sub-period are T 1 = 0.73 T s and T 2 = 0.27 T s . With these considerations, the particular cost function J ( k + 2 ) that will be used is defined as
J ( k + 2 ) = e ^ α β ( k + 2 ) 2 + λ s c S C ( k + 2 )
where . denotes vector modulus and λ s c is the WF that allows one to give more or less importance to the penalization of switching changes over the current tracking error.
The diagram in Figure 1 can be used with different sets of VVV. Several possibilities can be applied, as explained below. Figure 2 shows the voltage vectors, numbered from 0 to 63, for a six-phase VSI in the α β and x y planes. Four coronas (or dodecagons) are clearly visible. These are usually identified by the module of their voltage vectors in the α β plane. For each VV, the attribute (S, M, ML or L) is derived from its modulus, comparing it to the modulus of all the other VVs. For example, the VVs receiving the attribute Small are those whose modulus is the smallest possible. Those lie in the innermost corona in the alpha-beta sub-space. Similarly, the medium (M), medium-large (ML), and large (L) attributes for VV are defined (see Figure 2).
Similarly, VV combinations can be formed using S and ML vectors of the same angle. In this case, the sub-periods are 0.42 × T s for the S vector and 0.58 × T s for the ML vector. The set can be referred to as S + ML and again contains 12 elements. This set of VVV produces voltages that have on average less modulus than those of ML + L. This means that the DC link is not fully used and that some operating points may not be reachable despite the installed DC link voltage. Thus, the combination of S and ML voltage vectors has the potential to be used at low and medium load operating points. However, if one considers the set (S + ML) ∪ (ML + L) ∪ Z then 25 configurations are possible covering the whole operating interval.
With these definitions and observations, it is proposed to use the following sets:
  • Set 1. It contains 13 VVV obtained from the (ML + L) ∪ Z combinations.
  • Set 2. It contains 12 S + ML VVV in addition to the 13 VVV of the previous set, resulting in 25 VVV. It must be noted that set 2 contains set 1.
Table 1 summarizes the main aspects of the two VVV sets that will be used in this work, where C1 and C2 represent FSMPSCC using set 1 and set 2, respectively. Note that C1 and C2 need different computation times, resulting in different attainable sampling times. The sampling times depend on the hardware used, so the particular values of Table 1 have been obtained from the multi-phase system presented below.

Multi-Phase System under Analysis

The proposal will be implemented in a multi-phase system based on an asymmetrical six-phase IM, whose parameters are shown in Table 2. The electrical machine can be loaded independently of the speed, thanks to a coaxial DC machine that acts as a DC generator. A resistive variable load is connected to it. The DC link powers two three-phase Semikron VSI (SKS22F power converters), and the FSMPSCC technique using VVV is performed by a C program running on a TMS320F28335 DSP. Stator currents are measured using LEM LAH25-NP sensors, while mechanical speed is calculated using a GHM510296R/2500 digital encoder.
It is interesting to mention that identification plays a crucial role in this system. Model accuracy has a great impact on performance, as reported in several research articles [27,28]. Although model-free methods have been proposed [29], an IM model is used here. To generate the performance maps for many different operating points and WF values, it is necessary first to correctly identify the machine model. In this case, different techniques have been used to derive the IM model, including the AC time domain, and stand still with inverter supply tests. Additionally, results for FSMPSCC with VVV, using different WF tunings, have been generated for assessment. Figure 3 presents an example of the trajectories gathered in the experimental setup. The results correspond to i s α , i s β for a particular operating point with the C1 controller. From these trajectories, one can extract the values for the figures of merit that make up the performance maps.
After modeling and identification, a simulation environment is devised to be able to search for each operating point and WF value. This environment uses ODE solvers in MATLAB (R2017a version) for the continuous part of the dynamics, interacting with programs to emulate the controllers running at their specific sampling times. The computing time needed by the controller is carefully considered and tested against the Digital Signal Processor (DSP) used for the experiments.

3. FSMPSCC Assessment with Performance Maps

Two requirements are needed to use performance maps. First, a performance index or figure of merit must be defined, and second, a partition of the operating space must be provided. Regarding the performance index, in the case of FSMPSCC, the stator’s current tracking error and the switching rate at the VSI are usually considered. This would lead to several performance maps, one for each of the following variables: e α , e β , e x , e y , and f s w . However, the use of VVV takes care of the x y components to a satisfactory level, and the number of performance maps can be reduced to three variables: e α , e β , and f s w . Furthermore, the components e α and e β should be quite similar for normal drive operation, so they can be merged in a single quantity taking the mean. In this way, we will use e α β ( e α + e β ) / 2 , leaving the number of performance indices to two: e α β , and f s w . These quantities can be computed for each operating point by averaging over a time window ( k 1 , k 2 ) . The duration of this time window, L = k 2 k 1 + 1 , should include some electrical cycles to be meaningful. In mathematical notation, one gets the following:
e α β = 1 ( k 2 k 1 + 1 ) k = k 1 k 2 ( e α ( k ) + e β ( k ) ) / 2 2
f s w = 1 / 6 T s ( k 2 k 1 + 1 ) k = k 1 k 2 S C ( k ) .
Regarding the partition of the operating space, for the particular case of FSMPSCC, the amplitude for the stator reference I * and the electrical frequency f e define the operating point. Depending on the physical characteristics of the system, these variables take values at certain intervals: I * I a = [ 0 , A ¯ ] , f e I f = [ 0 , F ¯ ] . The set of all operating points Ω is included in the Cartesian product I a × I f .

Optimal Tuning

As previously indicated, tuning in FSMPSCC consists of selecting the WF for the CF (6). In [22], the authors propose the use of a different tuning for each operating point of the system. In this way, the compromise between conflicting criteria is tackled independently for each operating point. This is in contrast with the traditional FSMPSCC setup, where the WF is computed considering the range of operation of the system as a whole [9]. This method can be improved as it is realized that for different regions of operation, the terms present in (6) take different relative values. In the case at hand, it will be assumed that the VSI limitations impose a limit U on f s w . Then, the WF tuning should provide a f s w below the limit for all possible operating points, that is, f s w < U for all ( I * , f e ) Ω . On the other hand, the tuning should provide the minimum tracking error in α β for each operating point. This can be expressed as
e α β o = min λ s c e α β
As an example, the performance maps of Figure 4 are presented. They correspond to controller C1 with λ s c = 0.014 (left) and with an optimized λ s c (right). The first tuning ( λ s c = 0.014 ) provides a value of f s w 11 kHz for all operating points using a single value of the WF for the whole operating space; this represents the standard method of tuning the FSMPS. It is interesting to compare this case with the optimized λ s c (right), where a value of f s w 11 kHz is also achieved. The values of e α β are clearly better for the optimized case. This kind of analysis is allowed by using performance maps since all operating points are easily visualized.
In the following, optimized tuning will always be considered. The criterion is that of Equation (8). The method is the independent tuning of λ s c for each VVV set and for each operating point. Since there is only one tuning parameter, almost any optimization method can be used [30]. In our case, a first-order Newton-like method has been used. In order to manage the continuous operating space it is discretized by taking 12 sub-divisions of the load and 12 subdivisions of speed.
It must be stressed that in FSMPSCC, the controller configuration (selection of λ s c ) is not a trivial task, since the various performance indices ( e α β , f s w ) are affected by tuning in a way that depends on the particular operating point. As a result, some operating points might get a reduction in either index at the cost of an increase for other operating points. As an example, consider Figure 5 where the optimal values of λ s c are given for controller C1 for two scenarios: to attain f s w 11 kHz (left) and to attain f s w 8 kHz (right). A non-linear relationship between the optimal WF and the operating points is clearly visible. Please notice that for each operating point, there is a unique value for the WF. The existence of one or several peaks must not be interpreted as a choosing problem.
Performance maps can be used to compare different controllers even if they use different VVV sets. Using Equation (9), the performance map of Figure 6 is found. An optimal value of λ s c is used for each operating point while keeping f s w < 11 (kHz), a value used for illustrative purposes. A comparison with previous methods is provided. Firstly, comparing the proposal to the previous case of Figure 4, one gets a drastic reduction of e α β . Table 3 supplements the graphical information by including the minimum, mean, and maximum values for each case. Of course, the reduction in e α β comes at the expense of an increase in f s w . It has been pointed out in previous publications that this must be the case [8,21]. However, it must be remarked that the merit of this tuning is not the simultaneous reduction of all indices, but an optimized solution considering the objectives.
Figure 6. Performance map for the C2 controller, e α β and each operating point, with an optimal value for λ s c (a) and values of λ s c (b).
Figure 6. Performance map for the C2 controller, e α β and each operating point, with an optimal value for λ s c (a) and values of λ s c (b).
Applsci 13 04113 g006
Table 3. Summary of Results for Comparison with Previous Methods.
Table 3. Summary of Results for Comparison with Previous Methods.
e α β (mA) f sw (kHz)
Figureminmeanmaxminmeanmax
C1Figure 447.455.058.30.986.012.4
C1 (opt.)Figure 419.324.059.51.668.511
C2 (opt.)Figure 615.823.344.32.5110.011
C3 (opt.)Figure 715.822.144.32.519.711
Figure 7. Performance map for the proposed controller (C3) using the best set of VVV for each operating point.
Figure 7. Performance map for the proposed controller (C3) using the best set of VVV for each operating point.
Applsci 13 04113 g007
A similar procedure can be followed to assess the capabilities of C2. This controller uses a set of VVV composed of the S + ML and ML + L coronas, plus zero VV. Figure 6 shows the performance map for the C2 controller, e α β and each operating point, with an optimal value for λ s c .
These results are quite interesting and remarkable. First, it should be noted that most articles in the literature dealing with VVV focus on C1, and the C2 controller clearly reduces e α β at many points in the operating space. Second, one might argue that the C2 controller should always outperform C1, given the fact that it incorporates all the VVV that C1 uses plus the extra S + ML combinations. In other words, C2 has more degrees of freedom than C1. However, the reduced e α β is not verified for all operating spaces. The reason is that C2 has extra computations, compared to C1, and therefore cannot achieve the same sampling frequency, as shown in Table 1.

4. FSMPSCC with Two VVV Sets

In view of the results of the previous section, it is natural to consider the possibility of using the best VVV set for each point of operation. In this section, the idea is shown to be feasible, and the regions assigned to each controller in the operating space Ω are also stated to be easy to describe and use online.
The previously computed performance maps for C1 and C2 (Figure 4 and Figure 6) can be used to find the operating points where C2 outperforms C1. For these operating points, the optimal WF yields a lower value of e α β o for C2, that is, e α β o ( C 2 ) < e α β o ( C 1 ) . The region, from now on named Ω 2 , corresponds approximately to the operating points where 65 < I * × f e < 180 , as shown in Figure 8. This allows an easy implementation of the proposed FSMPSCC with two VVV sets.
Figure 9 shows a diagram of the proposed FSMPSCC with the two VVV sets. The additions to the diagram in Figure 1 are marked with a different color, to clarify the novelty and added complexity of our proposal. As stated previously, the outer controller is responsible for speed/torque regulation and usually operates on a larger time scale than the stator current control. This controller is responsible for issuing the reference values for the stator’s current control, I * and f e .
At this point, it is possible to check whether the actual pair ( I * , f e ) belongs to Ω 2 simply by computing I * × f e and comparing it with the limits indicated previously.
The extra operations needed to find the VVV set to be used at each sampling time are multiplication and two comparisons. These are easily implemented in any DSP. In addition, these additional computations impose a burden that is negligible compared to the rest of the operations: measurements, PI control, and optimization of the cost function. As a result, the proposed control scheme has almost the same computational load as the standard VVV method. For the particular hardware used in the experimental part (a TMS320F28335 DSP), the additional computational cost of the proposal is less thaN 0.5 μ s. Notice that the standard FSMPSCC algorithm takes on the order of 50 μ s. Furthermore, the operating point does not change every sampling period. This means that the outer controller can select which of the two constellations will be used.

Comparison with Previous Methods

The proposed controller using two VVV sets (from now on termed C3) will be compared with the help of performance maps. An optimal value for λ s c will be used for each point of operation. Figure 7 shows that C3 outperforms both C1 and C2. Furthermore, from Table 3 one can see that C3 retains the positive traits of C2 with respect to the minimum and maximum values of e α β and a reduced mean value for said e α β . The switching frequency remains below the limits required by the objectives.
The proposed C3 controller is now compared to the C1 controller using the standard VVV set (set 1). Figure 10 presents the trajectories for mechanical speed, stator currents on the d q axes, and stator phase currents for a particular operating point that falls in Ω 2 . The plot on the left corresponds to the standard C1 controller, and the plot on the right corresponds to C3. An improvement in current regulation is observed when using the proposal.
Regarding dynamic performance, Figure 11 shows the experimental results for a speed ramp test. It can be seen that the proposed scheme is able to deal with transient states. During the test, the reference speed is changed from 200 (rpm) to 400 (rpm). The tracking of the reference speed is satisfactory (see the upper plot in Figure 11), while stator phase currents perfectly evolve in the expected sinusoidal states.

5. Discussions

Modern electric machines are gaining great interest from the research community as a result of an increasing number of present-day applications, such as electric vehicles or wind power generators, among others. Different machines and control technologies are under study, being the combination of predictive controllers and multi-phase drives a particularly interesting case, where VVV has been applied to diminish the generated harmonic content in the stator current control. Although different combinations of VV can be used to define the applied VVV, those based on large and medium-large VV (set 1 in our analysis) are generally applied due to the exploitation of the DC link voltage. However, a different set of VVV can be defined, or combined with set 1 of VVV, to obtain an improvement in the stator current control of the multi-phase machine in all operating points. Our work focuses on the analytical study of different sets of VVV that can be applied to multi-phase drives, where performance maps and indices are applied and defined, and a particular six-phase IM is used as a case example. The results obtained conclude that the best stator current regulation is achieved from the definition of a novel control structure based on an adaptive online selection of the applied set of VVV. The proposed controller is finally implemented on an experimental test rig, where a digital signal processor (TMS320F28335) can execute the proposal at a negligible computational cost and experimental results are obtained for its validation. Note that the introduced method is extensible to different multi-phase systems, being the 6-phase IM just a particular one that requires dealing with a particular extra number of phases.

Author Contributions

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

Funding

This work is part of grant TED2021-129558B-C22 funded by MCIN/AEI/10.13039/ 501100011033 and by the “European Union NextGenerationEU/PRTR” and grant I+D+i/PID2021-125189OB-I00, funded by MCIN/AEI/10.13039/501100011033/ by “ERDF A way of making Europe”.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
C1Controller using the standard VVV set from the ML and L circles.
C2Controller using a VVV set with extra voltage vectors from the S and ML circles.
C3Proposed predictive controller using a VVV set adjusted for each operating point.
CFCost Function
DSPDigital Signal Processor
FSMPSCCFinite State Model Predictive Stator Current Control
IMInduction Machine
MBPCModel Based Predictive Control
PMPerformance Map
VVVoltage Vector
VVVVirtual Voltage Vector
VSIVoltage Source Inverter
WFWeighting Factors

References

  1. Wang, W.; Liu, C.; Liu, S.; Zhao, H. Model predictive torque control for dual three-phase PMSMs with simplified deadbeat solution and discrete space-vector modulation. IEEE Trans. Energy Convers. 2021, 36, 1491–1499. [Google Scholar] [CrossRef]
  2. Preindl, M.; Bolognani, S. Model predictive direct speed control with finite control set of PMSM drive systems. IEEE Trans. Power Electron. 2013, 28, 1007–1015. [Google Scholar] [CrossRef]
  3. Frikha, M.A.; Croonen, J.; Deepak, K.; Benômar, Y.; El Baghdadi, M.; Hegazy, O. Multiphase Motors and Drive Systems for Electric Vehicle Powertrains: State of the Art Analysis and Future Trends. Energies 2023, 16, 768. [Google Scholar] [CrossRef]
  4. Xue, Z.; Niu, S.; Chau, A.M.H.; Luo, Y.; Lin, H.; Li, X. Recent Advances in Multi-Phase Electric Drives Model Predictive Control in Renewable Energy Application: A State-of-the-Art Review. World Electr. Veh. J. 2023, 14, 44. [Google Scholar] [CrossRef]
  5. Duran, M.J.; Gonzalez-Prieto, I.; Gonzalez-Prieto, A.; Aciego, J.J. The Evolution of Model Predictive Control in Multiphase Electric Drives: A Growing Field of Research. IEEE Ind. Electron. Mag. 2022, 16, 29–39. [Google Scholar] [CrossRef]
  6. Elmorshedy, M.F.; Xu, W.; El-Sousy, F.F.; Islam, M.R.; Ahmed, A.A. Recent Achievements in Model Predictive Control Techniques for Industrial Motor: A Comprehensive State-of-the-Art. IEEE Access 2021, 9, 58170–58191. [Google Scholar] [CrossRef]
  7. Liu, C.; Luo, Y. Overview of advanced control strategies for electric machines. Chin. J. Electr. Eng. 2017, 3, 53–61. [Google Scholar]
  8. Arahal, M.R.; Satué, M.G.; Barrero, F.; Ortega, M.G. Adaptive Cost Function FCSMPC for 6-Phase IMs. Energies 2021, 14, 5222. [Google Scholar] [CrossRef]
  9. Lim, C.S.; Levi, E.; Jones, M.; Rahim, N.; 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. 2014, 61, 3867–3878. [Google Scholar] [CrossRef] [Green Version]
  10. Shawier, A.; Habib, A.; Mamdouh, M.; Abdel-Khalik, A.S.; Ahmed, K.H. Assessment of predictive current control of six-phase induction motor with different winding configurations. IEEE Access 2021, 9, 81125–81138. [Google Scholar] [CrossRef]
  11. Yepes, A.G.; Shawier, A.; Abdel-Azim, W.E.; Abdel-Khalik, A.S.; Ahmed, S.; Doval-Gandoy, J. General Online Current-Harmonic Generation for Increased Torque Capability With Minimum Stator Copper Loss in Fault-Tolerant Multiphase Induction Motor Drives. IEEE Trans. Transp. Electrif. 2023. [Google Scholar] [CrossRef]
  12. Xue, C.; Song, W.; Feng, X. Finite control-set model predictive current control of five-phase permanent-magnet synchronous machine based on virtual voltage vectors. IET Electr. Power Appl. 2017, 11, 836–846. [Google Scholar] [CrossRef]
  13. Garcia Entrambasaguas, P.; Gonzalez-Prieto, I.; Duran, M.; Bermudez, M.; Barrero, F. Direct torque control based on virtual voltage vector for a six-phase induction machine. Rev. Iberoam. Autom. Inform. Ind. 2018, 15, 277–285. [Google Scholar]
  14. Gonzalez, O.; Ayala, M.; Romero, C.; Delorme, L.; Rodas, J.; Gregor, R.; Gonzalez-Prieto, I.; Duran, M. Model Predictive Current Control of Six-phase Induction Motor Drives using Virtual Vectors and Space Vector Modulation. IEEE Trans. Power Electron. 2022, 37, 7617–7628. [Google Scholar] [CrossRef]
  15. Gonçalves, P.F.; Cruz, S.M.; Mendes, A.M. Bi-subspace predictive current control of six-phase PMSM drives based on virtual vectors with optimal amplitude. IET Electr. Power Appl. 2019, 13, 1672–1683. [Google Scholar] [CrossRef]
  16. Luo, Y.; Liu, C. A flux constrained predictive control for a six-phase PMSM motor with lower complexity. IEEE Trans. Ind. Electron. 2018, 66, 5081–5093. [Google Scholar] [CrossRef]
  17. Ren, Z.; Ji, J.; Tang, H.; Tao, T.; Huang, L.; Zhao, W. Model Predictive Torque Control for a Dual Three-phase PMSM Using Modified Dual Virtual Vector Modulation Method. Chin. J. Electr. Eng. 2022, 8, 91–103. [Google Scholar] [CrossRef]
  18. Zhou, H.; Li, H.; Xiang, X.; Yuan, B.; Zhou, T.; Li, W. Model Predictive Control Algorithm of Dual Three-Phase Motor Considering Global Single Vector. In Proceedings of the 2022 25th International Conference on Electrical Machines and Systems (ICEMS), Chiang Mai, Thailand, 29 November–2 December 2022; pp. 1–6. [Google Scholar]
  19. Mamdouh, M.; Abido, M.A. Simple Predictive Current Control of Asymmetrical Six-Phase Induction Motor with Improved Performance. IEEE Trans. Ind. Electron. 2022, 70, 7580–7590. [Google Scholar] [CrossRef]
  20. Yao, M.; Peng, J.; Sun, X. Model predictive flux control of six-phase permanent magnet synchronous motor with novel virtual voltage vectors. Electr. Eng. 2022, 104, 2835–2847. [Google Scholar] [CrossRef]
  21. Fretes, H.; Rodas, J.; Doval-Gandoy, J.; Gomez, V.; Gomez, N.; Novak, M.; Rodriguez, J.; Dragičević, T. Pareto Optimal Weighting Factor Design of Predictive Current Controller of a Six-Phase Induction Machine based on Particle Swarm Optimization Algorithm. IEEE J. Emerg. Sel. Top. Power Electron. 2021, 10, 207–219. [Google Scholar] [CrossRef]
  22. Arahal, M.R.; Kowal, G.; Barrero, F.; del Mar Castilla, M. Cost function optimization for multi-phase induction machines predictive control. Rev. Iberoam. Autom. Inform. Ind. 2019, 16, 48–55. [Google Scholar] [CrossRef]
  23. Li, L.; Zhou, W.; Bi, X.; Sun, X.; Shi, X. Second-Order Model-Based Predictive Control of Dual Three-Phase PMSM Based on Current Loop Operation Optimization. Actuators 2022, 11, 251. [Google Scholar] [CrossRef]
  24. Yaramasu, V.; Wu, B. Model Predictive Control of Wind Energy Conversion Systems; John Wiley & Sons: Hoboken, NJ, USA, 2016. [Google Scholar]
  25. Yao, M.; Peng, J.; Sun, X. Robust model predictive current control for six-phase PMSMS with virtual voltage vectors. J. Power Electron. 2022, 22, 1885–1896. [Google Scholar] [CrossRef]
  26. Camacho, E.F.; Bordons, C. Model Predictive Control; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  27. Arahal, M.R.; Barrero, F.; Satué, M.G.; Ramírez, D.R. Predictive Control of Multi-Phase Motor for Constant Torque Applications. Machines 2022, 10, 211. [Google Scholar] [CrossRef]
  28. Sawma, J.; Khatounian, F.; Monmasson, E.; Idkhajine, L.; Ghosn, R. Analysis of the impact of online identification on model predictive current control applied to permanent magnet synchronous motors. IET Electr. Power Appl. 2017, 11, 864–873. [Google Scholar] [CrossRef]
  29. Zhang, Y.; Jin, J.; Huang, L. Model-free predictive current control of PMSM drives based on extended state observer using ultralocal model. IEEE Trans. Ind. Electron. 2020, 68, 993–1003. [Google Scholar] [CrossRef]
  30. Arahal, M.R.; Martín, C.; Kowal, A.; Castilla, M.; Barrero, F. Cost function optimization for predictive control of a five-phase IM drive. Optim. Control. Appl. Methods 2020, 41, 84–93. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Diagram of FSMPSCC of a 6-phase IM using a single VVV set.
Figure 1. Diagram of FSMPSCC of a 6-phase IM using a single VVV set.
Applsci 13 04113 g001
Figure 2. Voltage vectors, numbered from 0 to 63, for a six-phase VSI shown in α β plane (left) and x y plane (right). The S, M, ML and L coronas have been highlighted in the α β plane.
Figure 2. Voltage vectors, numbered from 0 to 63, for a six-phase VSI shown in α β plane (left) and x y plane (right). The S, M, ML and L coronas have been highlighted in the α β plane.
Applsci 13 04113 g002
Figure 3. Example of trajectories gathered in the experimental setup for model tuning.
Figure 3. Example of trajectories gathered in the experimental setup for model tuning.
Applsci 13 04113 g003
Figure 4. Performance maps for e α β for the C1 controller, with λ s c = 0.014 (a) and with an optimized λ s c (b).
Figure 4. Performance maps for e α β for the C1 controller, with λ s c = 0.014 (a) and with an optimized λ s c (b).
Applsci 13 04113 g004
Figure 5. Optimal values of λ s c to attain f s w 11 kHz (a) and to attain f s w 8 kHz (b).
Figure 5. Optimal values of λ s c to attain f s w 11 kHz (a) and to attain f s w 8 kHz (b).
Applsci 13 04113 g005
Figure 8. The space between the lines is region Ω 2 , where C2 outperforms C1 according to the selected criteria.
Figure 8. The space between the lines is region Ω 2 , where C2 outperforms C1 according to the selected criteria.
Applsci 13 04113 g008
Figure 9. Diagram of the proposed FSMPSCC with two VVV constellations.
Figure 9. Diagram of the proposed FSMPSCC with two VVV constellations.
Applsci 13 04113 g009
Figure 10. Comparison with previous methods in terms of trajectories of control variables for several speeds. In all cases the the standard C1 controller is on the (left) and the proposed C3 on the (right).
Figure 10. Comparison with previous methods in terms of trajectories of control variables for several speeds. In all cases the the standard C1 controller is on the (left) and the proposed C3 on the (right).
Applsci 13 04113 g010
Figure 11. Trajectories of speed and currents for a ramp test (left) and a speed reversal test (right).
Figure 11. Trajectories of speed and currents for a ramp test (left) and a speed reversal test (right).
Applsci 13 04113 g011
Table 1. VVV constellations.
Table 1. VVV constellations.
ControllerVVV SetSet Number N v T s ( μ s)
C1(ML + L) ∪ Z11336
C2(S + ML) ∪ (ML + L) ∪ Z22551
Table 2. Parameters of the 6-phase IM drive.
Table 2. Parameters of the 6-phase IM drive.
V D C 200 (V) I p e a k 4.1 (A) n m 1000 (rpm)
R s 4.2 ( Ω ) R r 2 ( Ω ) L m 420 (mH)
L l s 1.5 (mH) L l r 55 (mH)P3 (pairs)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Arahal, M.R.; Satué, M.G.; Barrero, F.; Martín, C. Online Adaptive Set of Virtual Voltage Vectors for Stator Current Regulation of a Six-Phase Induction Machine Using Finite State Model Predictive Controllers. Appl. Sci. 2023, 13, 4113. https://doi.org/10.3390/app13074113

AMA Style

Arahal MR, Satué MG, Barrero F, Martín C. Online Adaptive Set of Virtual Voltage Vectors for Stator Current Regulation of a Six-Phase Induction Machine Using Finite State Model Predictive Controllers. Applied Sciences. 2023; 13(7):4113. https://doi.org/10.3390/app13074113

Chicago/Turabian Style

Arahal, Manuel R., Manuel G. Satué, Federico Barrero, and Cristina Martín. 2023. "Online Adaptive Set of Virtual Voltage Vectors for Stator Current Regulation of a Six-Phase Induction Machine Using Finite State Model Predictive Controllers" Applied Sciences 13, no. 7: 4113. https://doi.org/10.3390/app13074113

APA Style

Arahal, M. R., Satué, M. G., Barrero, F., & Martín, C. (2023). Online Adaptive Set of Virtual Voltage Vectors for Stator Current Regulation of a Six-Phase Induction Machine Using Finite State Model Predictive Controllers. Applied Sciences, 13(7), 4113. https://doi.org/10.3390/app13074113

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