# A Multilevel Inverter Bridge Control Structure with Energy Storage Using Model Predictive Control for Flat Systems

Authors: Paolo Mercorelli
Journal: Journal of Engineering
 (2013)

---

## Abstract

The paper presents a novel technique to control the current of an electromagnetic linear actuator
fed by a multilevel IGBT voltage inverter with dynamic energy storage. The technique uses a “cascade model predictive control (MPC),” which consists of two MPCs. A predictive control of the trajectory position predicts the optimal current, which is considered to be the desired current for the second MPC controller in which a hysteresis control technique is also integrated. Energy is stored in a capacitor using energy recovery. The current MPC can handle a capacitor voltage higher than the source voltage to guarantee high dynamic current and disturbance compensation. The main contribution of this paper is the design of an optimal control structure that guarantees a capacitor recharge. In this context, the approach is quite new and can represent a general emerging approach allowing to reduce the complexity of the new generation of inverters and, in the meantime, to guarantee precision and acceptable switching frequency. The proposed technique shows very promising results through simulations with real actuator data in an innovative transportation technology.


---

## Body

## 1. Introduction and Motivation

The performance of a multilevel inverter is better than that of a classical inverter. The total harmonic distortion of a classical inverter is very high; that is, the total harmonic distortion for a multilevel inverter is low. This topic has become the focus of a significant amount of research. One study has simulated and implemented a multilevel inverter-fed induction motor drive [1]. The output harmonic content is reduced by using a multilevel inverter. In a symmetrical circuit, the voltage and power increase as the number of levels in the inverter increases. The switching angle for the pulse is selected to reduce the harmonic distortion. This drive system has several advantages, including reduced total harmonic distortion and higher torque. A normal neutral point potential stabilization technique using the output current polarity has been proposed [2]. The neutral point potential balancing algorithm for three-level neutral point clamped inverters using an analytically injected zero-sequence voltage has been developed [3]. Modulation schemes to eliminate common mode voltage in multilevel inverter topology have been suggested [4]. A generalized multilevel inverter topology with self-voltage balancing has also been suggested [5]. A survey of topologies, control, and applications of multilevel inverters has been published [6]. A digital modulation technique for dual three-phase alternating current (AC) machines has been presented [7]. A space vector pulse-width modulation (PWM) technique for a dual three-phase AC machine and its digital signal processor (DSP) implementation has been presented [8]. Practical medium voltage converter topologies for high power applications have been described previously [9–11]. Industrial topologies have been presented [12]. Different pulse-width modulation techniques for symmetrical cascade inverters with high and fundamental switching frequency have been shown, for instance, in [13, 14]. An improvement in terms of efficiency of the converter is presented in [15] where the authors use different DC sources to reduce the switching losses. More recently in [16], the authors proposed a novel H-bridge multilevel pulse modulation converter topology in which the structure is based on a series connection of high-voltage diode-clamped inverter and low-voltage conventional inverter. The present work implements a controlled multilevel inverter feeding an electromagnetic actuator. In particular, two MPCs in cascade structure are presented to control five-level inverter electromagnetic actuators which are becoming increasingly important in many industrial applications. Thanks to recent progress in permanent-magnet technology, very compact and high-power electromagnetic actuators are now available. Particularly in automotive systems, these actuators have often replaced conventional mechanical components due to their high efficiency, excellent dynamic behavior, and control flexibility; thus, mechatronics is one of the most notable innovative fields in the automobile industry. A recent design of such an actuator (e.g., [17]) considers optimizing with respect to dimensions and forces. Nevertheless, practical tests show that this actuator structure is not always effective for a large variety of applications because a complex control structure is needed to control such a valve [18]. Recently, for instance, a new type of electromagnetic valve drive system has been proposed and described [19]. This paper presents the control design for a novel permanent-magnet linear valve actuator for use in a variable engine control to allow short-stroke and high-dynamic motions. In a linear actuator, the electromagnetic force is totally independent of the constrained motion. Specifically, the electromagnetic force is parallel to the magnetic axis, and thus it is orthogonal to the constrained motion. These performance characteristics make this type of actuator an excellent candidate in applications requiring both high precision and high-load capacity. In order to obtain an actuator structure that could be easily controlled, various linear designs with different permanent magnet actuator topologies using NdFeB magnets were considered. In particular, the objective of this paper is to design a controller structure consisting of (i)
a flatness-based feedforward control,(ii)
a multilevel inverter configuration,(iii)
a cascade MPC structure based on optimized energetic functions to control the proposed multilevel inverter with storage energy.An important aspect of this approach is the use of geometric system properties, such as the differential flatness needed to track the desired trajectories. In this context, recently published research results have been considered (e.g., [20]) in which the author proposed an interesting method to integrate MPC and a feedforward flatness control. The second part of the paper shows the structure of the cascade MPC defined above. For trajectory-based motion cycles, there is usually a need to control the actuator current separately. The reference value for the actuator current is normally generated by the speed or position controller in an outer control loop based on the desired speed or position profiles. The current control is used to provide the force or torque required for the desired motion quickly and precisely. Generally, three main techniques are employed for current control of a voltage-source inverter (VSI):(i)
closed loop control (e.g., PID control) using pulse-width modulated terminal voltages [21],(ii)
hysteresis control techniques [22, 23],(iii)
predictive current control [24–26].The PWM technique is currently the most popular method for an inverter current control. The main advantage of PWM techniques is that the inverter switches operate at a fixed frequency, and hardware- or software-based standard modulators are available as industrial products (e.g., on-board microcontrollers). A staggered space vector modulation technique applicable to three-phase cascaded voltage-source inverter topologies has also been demonstrated with a single-phase cascaded voltage-source inverter that uses a series connection of insulated gate bipolar transistor (IGBT) H-bridge modules with isolated DC buses [27]. Among various modulation techniques for a multilevel inverter, the space vector pulse width modulation (SVPWM) is widely used. However, implementing the SVPWM for a multilevel inverter is complicated because it is difficult to determine the location of the reference vector, calculate on times, and determine and select switching states. A previous paper has proposed a general SVPWM algorithm for multilevel inverters based on standard two-level SVPWM. However, the system response is affected by the stability requirements of the feedback loop, which also depend on load parameters; moreover, appreciable phase lag may arise even in the steady state. Recently, many studies have focused on solving these problems (e.g., [28, 29]) of multilevel inverter modules with independent control of the phase angle and magnitude of the output voltage [29]. Control and new power bridge structures have been studied. Very recent works have presented a novel bridge structure [30, 31]. In both structures, capacitors are used but not recharged. In particular, an MPC strategy is proposed for a three-phase power bridge [30]. Hysteresis current control has a fast response and a good accuracy. It can be implemented with a simple hardware structure, and, in many cases, it does not require any knowledge of load parameters. However, it can sometimes cause very high switching frequencies or, if the maximum switching frequency is limited, the current waveform may vary widely and the current peaks may appreciably exceed the hysteresis band depending on the operation conditions and load parameters. Conventional predictive current control uses a simple gradient model to predict load current in vector space and determine the proper switching voltage vector based on the one-step prediction [24–26]. However, no feedback loop is applied to compensate for model uncertainties. To overcome the disadvantages of the mentioned methods, this study develops a novel approach that combines hysteresis control with a model predictive control (MPC) strategy. The proposed control system integrates flatness and two cascade MPCs. The first MPC generates an optimal desired current by minimizing position error. The second MPC considers this optimal desired current to be the reference signal to find an optimal switching law in combination with the hysteresis control. A special energy storage and charging circuit was designed to provide multilevel voltages for an effective current control. By applying the proper voltage level, which was determined by the MPC strategy, to the actuator, dynamic current changes and small current ripples are possible in spite of relatively low switching frequencies. A recent work [32] proposes “fast” algorithms based on computing the formulated problems in parallel. In fact, the optimization technique uses subsets of the state subspace offline and then optimizes them in parallel. In this paper, a cascade structure of MPCs is utilized to control a multilevel inverter. Two different optimization techniques are presented as well. The main contribution of the paper is the optimal control structure, which also guarantees that the capacitor recharges. In this context, the approach is quite new. The paper is organized in the following way.Section2 describes the model. In Section 3, the flatness property of the system is shown. In Section 4, the general MPC goal is defined and the main idea of using a cascade positional MPC followed by a current MPC is explained. In Section 5, a positional MPC around the desired trajectory is proposed. Section 6 is devoted to analyzing and optimizing the current MPC. The simulation results and some concluding remarks end the paper.


## 2. Description of the Physical Systems

The electromagnetic actuator is depicted in Figure1. The structure of this actuator is a moving coil; thus, the coils are mounted on the upper part of the stem of the valve. The moving valve was connected to the power system to feed the coil with two normal cables with fix contacts. This arrangement is possible because of the short stroke to be covered (8 mm). The valve can be modeled mathematically in the following way:
(1)∂icoil(t)∂t=-RcoilLcoilicoil(t)+Uin(t)-uq(icoil(t),s(t))Lcoil,(2)∂s(t)∂t=v(t),(3)∂v(t)∂t=Dg(icoil(t),s(t))micoil(t)+-kdv(t)-kfs(t)+F0(t)m,
where
(4)Dg(icoil(t),s(t))=Φh(icoil(t),s(t))lcoilAcoil,uq(icoil(t),s(t))=Dg(icoil(t),s(t))v(t).
It is important to note that Dg(icoil(t),s(t))>0  for all  icoil(t) and for all  s(t). Rcoil and Lcoil are the resistance and the inductance of the coil windings; Uin(t) is the input voltage; uq(t) is the induced back voltage; Φh is the magnetic flux penetrating the coil; icoil(t) is the coil current; lcoil is the coil length. s(t), v(t), and m are the position, velocity, and mass of the actuator, respectively, while kdv(t), kfs(t), and F0(t) represent the viscose friction, the total spring force, and the disturbing force acting on the valve. Equation (1) represents the electrical system of the actuator. Both (2) and (3) describe the mechanical behavior of the actuator, and (2) and (4) also represent the magnetic system. In particular, the following expression
(5)FL(t)=Dg(icoil(t),s(t))icoil(t)
describes the Lorentz force generated by the actuator. Essentially, the magnetic flux, generated by the permanent magnets, has two components in the air gap: (1) the main flux Φh(icoil(t)), which does not depend on the displacement of the mover and is responsible for the Lorentz force and the induced back voltage, and (2) the leakage flux Φσ(icoil(t),s(t)), which disperses around the coil and does not contribute to the electromagnetic force and induced back voltage:
(6)Φh(icoil(t))+Φσ(icoil(t),s(t))=ΦPM(icoil(t)).
However, due to the special actuator design, the leakage flux is almost equal to zero. Thus, it is possible to conclude
(7)Dg(icoil(t),s(t))≈Dg(icoil(t))=ΦPM(icoil(t))lcoilAcoil.Figure 1
Cross-section of the perpendicular linear actuator.


## 3. Differential Flatness of the System

Roughly speaking, a system is differentially flat if it is possible to find a set of outputs equal in number to the number of inputs such that all states and inputs are expressed in terms of those outputs and their derivatives. To be more precise, if the system has state variablesx∈ℜn and inputs u∈ℜm, then the system is flat if the outputs y∈ℜm have the following form:
(8)y=y(x,u,u˙,…,u(p)),(9)x=x(y,y˙,…,y(q));u=u(y,y˙,…,y(q+1)).Differentially flat systems are especially interesting in situations in which explicit trajectory tracking is required. Because the behavior of the flat system is given by the flat output, it is possible to plan trajectories in output space and then map them to the appropriate inputs. Concerning the flat output, for dating neither the necessary and sufficient condition nor a general method for the determination of the flat output has been provided. Typically, a flat output isguessed, and Definitions (8) and (9) are used to verify that the chosen output is flat. Once the system is proven to be differentially flat, the basic approach of two degrees of freedom controller design consists of two steps: first, separating the nonlinear controller synthesis problem into designing a feasible feedforward controlled trajectory for the nominal model of the system and second regulating that trajectory using controllers, that guarantee robust performance in the presence of uncertainties and disturbances. To explain the meaning of flatness intuitively, it is adequate to consider a system to be flat when the dynamic of at least one output is “visible” from the state and the input of the system; thus, there is at least one output that can be functionally controlled by the inputs as well by the states.Now the first step is to show that our model represented in (1), (2), and (3) is flat. If the position s of the moving part of the actuator is chosen as “guessed output,” then
(10)y(t)=s(t)⇒s(t)=y(t)s˙(t)=y˙(t)⇒v(t)=y˙(t),(11)icoil(t)Dg(icoil(t))=m(kdmv(t)+kfms(t)+∂v(t)∂t),(12)Uin(t)=Lcoil∂icoil(t)∂t+Rcoilicoil(t)+Dg(icoil(t))v(t).
To verify the flatness property it is necessary to look at (13). As evident from (7), Dg(icoil(t)) is proportional to the flux coupled with the permanent magnets which is again proportional to the flux density BPM. The operating point (BPM,HPM) of the permanent magnets is mainly determined by the magnet geometry and the demagnetization curve (Figure 2). The coil current only changes the field strength around the operating point (ΔH) algebraically. An exact analytical expression between ΔH and icoil(t) is very difficult to be derived because it also depends on the saturation level of the iron parts. However, a numerical solution exists. In other words, there is always a (locally) unique solution for icoil(t) that depends on Dgicoil(t). Thus, using (13), the coil current can always be determined uniquely by the actuator position and its derivatives. Considering additional system equations, it can be easily seen that the flatness conditions defined in (9) are satisfied. Relationships (13) and (12) define the inverse system. Based on the flatness of the system, the linearizing trajectory of the model (3) and (2) is as follows:
(13)id(t)=mDg(id(t))(kdmvd(t)+kfmsd(t)+∂vd(t)∂t),
where sd(t) is the desired trajectory and id(t) is the corresponding feed forward control current. As it was explained before, the flux does not depend on the position of the armature because the leakage flux is almost equal to zero. This means that the flux is a function of the current. Function Dg(icoil(t)) is never equal to zero, at least for our conceived structure, because of controllability. The actuator is conceived and designed in a way in order to guarantee the controllability at any point of its movement. For (1), the inverse system that determines the desired feedforward control input ud(t) is defined as follows:
(14)ud(t)=Lcoil∂id(t)∂t+Rcoilid(t)+Dg(id(t),s(t))∂sd(t)∂(t).
Based on the feed forward control presented previously, the next step is to build a feedback control law that, in the presence of external disturbances and uncertainties in the parameters, takes the system around the desired trajectory.Figure 2
Permanent-magnet demagnetization curve.


## 4. General MPC Problem Formulation

Figure3 shows the power electronic circuit feeding the permanent magnet linear actuator. The power supply is a constant voltage source US (battery). However, to realize a current change as fast as possible to ensure the required dynamic, it is desirable to have a higher voltage applied to the actuator. Therefore, a capacitor C, charged with an initial voltage Uc higher than the battery level, is used. To maintain the high voltage level of the capacitor during the valve operation, an energy storage block consisting of the capacitor C, the diode D5, and the IGBT switch V5 is used. Besides the current regulation, another task of the control strategy is to recharge the capacitor appropriately by utilizing the braking energy. On the right, the IGBT inverter bridge and the linear actuator are represented. Model predictive control has been a widely used control concept for over 15 years, especially in the process industry. Applications of MPC in the field of electrical drives are quite rare because of a high computational complexity. The central point of a model predictive controller is a process model, which is capable to predict future output signals based on future input signals and initial values. Using this process model, the future dynamic behavior of the real plant is predicted within a prediction horizon. These predicted output signals can be used to minimize an open loop performance criterion (e.g., the sum of squared control errors within the prediction horizon) and to calculate the input signals u(k) for a control horizon. Outside the control horizon, the input remains constant. The calculated input signals feed the plant until a new measurement becomes available. This procedure is repeated with a receding prediction and control horizon. The receding horizon strategy makes a closed loop control law from the originally open loop minimization. The minimization step can easily include constraints, such as input, output, or state constraints, that may be taken into account already in the controller design. The principle of the predictive control design for the current control problem considered in this paper is depicted in Figure 4 with the desired current (dashed-point line) that results from the outer position control loop and the measured current (solid line). At the present time t, all possible currents for the future two time steps t+ts⋯t+2ts are calculated. This calculation is performed based on a multimode model and assessed afterwards with the help of a defined cost function (details follow). The optimum switching configuration is also shown (bold-dashed line).Figure 3
Multilevel inverter bridge structure.Figure 4
Principle of the complete applied MPC strategy.Figure5 shows the whole proposed control structure. The presented technique can be interpreted in the following points.(i)
The positional MPC structure generates an optimal current trajectory with the desired optimal horizon through the electrical model.(ii)
The optimal current trajectory is used as a desired trajectory for the current MPC.(iii)
The current MPC makes it possible to calculate the switch control system structure.Figure 5 shows the control system structure in which it is possible to recognize the desired current id(t) which is used to linearize the system, the predicted reference current iopt(t), and current icoil(t).Figure 5
The whole control scheme.


## 5. Solving a Linear Position MPC Optimization Problem

To obtain the desired current trajectory, a position model predictive control is implemented. In this case, the sample time is relatively short, 40μs (25 kHz), to make the algorithm as fast as possible. Considering the model in (1), (2), and (3) in which F0(t)=0 and an Euler discretization with k=nts, n∈N, where ts is the sampling time, the following system is obtained:
(15)icoil(k+1)=icoil(k)+tsLcoil(-Rcoilicoil(k)+umpc(k)+ud(k)-uq(k))s(k+1)=s(k)+tsv(k)v(k+1)=v(k)+ts(Dg(icoil(k))micoil(k)v(k)+tsc+-kdv(k)-kfs(k)m),
where umpc(k) is the model predictive control input to be calculated and ud(k) and uq(k) are the discretized voltages of ud(t) and uq(t) already defined in (1), (2), and (3). If (13) is discretized, then
(16)id(k)=mDg(id(k-1))·(kdmvd(k)+kfmsd(k)+vd(k+1)-vd(k)ts).
Considering the Euler discretization of (14), which allows us to express the future value of an output as a function of the pass inputs, which is the required form of the MPC approach, the following equation is obtained:
(17)ud(k)=Rcoilid(k)+Lcoil(id(k+1)-id(k)ts)+uq(k).
If (16) and (17) are inserted into the discretized equations (1), (2), and (3), then the following linear system is obtained:
(18)icoil(k+1)=icoil(k)-id(k)+id(k+1)-tsRcoilLcoil(icoil(k)-id(k)-umpc(k)Rcoil)s(k+1)=s(k)+tsv(k)v(k+1)=v(k)+ts(kdvd(k)+kfsd(k)mv(k)+tscc-kdv(k)+kfs(k)m).
If Δumpc(k)=umpc(k)-umpc(k-1) is assumed, the system described in (18) becomes
(19)x(k+1)=Akx(k)+Bk(Δumpc(k)+umpc(k-1))+Ekd(k)y(k)=Hkx(k),
where matrix Hk=[010] is the output matrix which determines the position, and the previous notation means
(20)Ak=[(1-tsRcoilLcoil)0001ts0-tskfm(1-tskdm)],Bk=[tsLcoil00],Ek=[1(tsRcoilLcoil-1)00000000tskfmtskdm],d(k)=[id(k+1)id(k)sd(k)vd(k)].
In the considered representation, vector d(k) is the known input. In the model approach, just two samples are considered:(21)s^(k+1)=HkAkx(k)+HkBkΔumpc(k)+HkBkΔumpc(k-1)+HkEkd(k),s^(k+2)=HkAk2x(k)+HkAkBkΔumpc(k)+HkAkBkΔumpc(k-1)+HkBkΔumpc(k+1)+HkBkΔumpc(k)+HkAkEkd(k)+HkEkd(k).
If the following performance criterion is assumed,
(22)J=12∑j=1N‍(sd(k+j)-s^(k+j))T·Qp(sd(k+j)-s^(k+j))+∑j=1N‍(Δumpc(k+j))T·RpΔumpc(k+j),
where sd(k+j), j=1,2,…,N, is the position reference trajectory, N the prediction horizon, and Qp and Rp are nonnegative definite matrices, then the solution minimizing performance index (22) may be then obtained by solving
(23)∂J∂Δumpc=0.
If only two steps for the prediction horizon are considered, then
(24)Gp=[HkAkHkAk2],F1p=[HkBk0HkAkBk+HkBkHkBk],F2p=[HkBkHkAkBk+HkBk],F3p=[HkEk0HkAkEk+HkEkHkEk].
A direct computation may be obtained explicitly as
(25)Δumpc(k)=(F1pTQpF1p+Rp)-1·(F1pTQp(Ydp(k)-Gpx(k)F1pTQpccccc-F2pumpc(k-1)-F3pd(k)F1pTQp(Ydp(k)-Gpx(k))),
where Ydp(k) is the desired output column vector. Once the predicted optimized voltage is obtained, then it is possible to obtain the predicted optimized current as the reference current for the current MPC based on the model of the system.Remark 1.
It should be noted that, because of the linearization around a trajectory, the model depends on the value of the trajectory, which is formalized in vectord(k).


## 6. Inverter and Motor Electrical Model in Current MPC Structure

Given a desired current reference, the goal is to control the real actuator current following this signal by a switched voltage. Direct control of the switching states has some advantages over the conventional pulse-width modulation method: the reachable dynamic is higher, and no pulse-pattern generation is needed. However, a proper switching strategy is necessary, and chattering phenomena can occur due to measurement noise. To provide a good reference, the following control problem is stated: given the system depicted in Figure3 and the following linear system for the actuator:(26)∂icoil(t)∂t=-RcoilLcoilicoil+𝒰in-uq(k)Lcoil,
where 𝒰in∈[u1,u2,u3,u4,u5]=ℐ is the available set of input voltages applied to the actuator and uq(k) is the induced back voltage of the linear actuator. Find a sequence 𝒰* of p elements u∈ℐ such that
(27)J=min𝒰in12∑j=1Np‍(i^coil(k+j)-iopt(k+j))Qimin𝒰in12∑j=1Np‍·(i^coil(k+j)-iopt(k+j))+∑k=1Np‍12(ΔUc(k+j))Ri·(ΔUc(k+j)),
where i^coil(k) is the predicted current, iopt(k) is an optimized current which is assumed to be a reference current (this current is calculated by the positional MPC structure), and Np is the prediction horizon (equal to the control horizon in our case). Qi is the selection matrix of the weights of the current error, and Ri is a selection matrix of the weights of the input voltage.It is worthwhile to remark that the matricesQi provide adaptive weights depending on the length of the prediction horizon. The formulation in (26) and (27) is a well-known system with switching command variables. In our case, if a short prediction horizon (2–4 samples) is considered, the induced back voltage uq could, in this interval, be treated as constant. Observing Figure 3, five working phases are identified. These working phases are summarized as follows.(i)
Working phase one: IGBT 2 and IGBT 3 on; IGBT 1 and IGBT 4 off; IGBT 5 off (Figure6).(ii)
Working phase two: IGBT 2 and IGBT 3 off; IGBT 1 and IGBT 4 on; IGBT 5 off (Figure6).(iii)
Working phase three: IGBT 1 and IGBT 4 on; IGBT 2 and IGBT 3 off; IGBT 5 on (Figure7).(iv)
Working phase four: IGBT 1 and IGBT 4 off; IGBT 2 and IGBT 3 on; IGBT 5 on (Figure7).(v)
Working phase five: IGBT 2 and IGBT 3 off; IGBT 1 and IGBT 4 off; IGBT 5 on (Figure8).MPC Representation. In general, dependence on the switching configuration the discrete-time dynamic model of the considered system can be formulated in one of the following two ways:
(28)icoil(k+1)=A1icoil(k)+B1(Us-uq(k)),
where the pair (A1,B1) represents the first order “R-L” system depicted in Figure 6 or by defining
(29)x(k)=[icoil(k)Uc(k)].
The dynamic of phases three and four is
(30)x(k+1)=A2x(k)+B2(-uq(k)),y(k)=C2x(k)=icoil(k),
where the term (A2,B2,C2) represents the second order “R-L-C” system depicted in Figure 8.(i)
During the commutation from phases “one and two” to phases “three and four,” the predicted current with a prediction horizon equal to 2 is(31)icoil(k+2)=C2A2[A1icoil(k)+B1(Us-uq(k))Uc(k)]+C2B2(-uq(k)),where Uc(k) is the measured capacitor voltage, icoil(k) is the measured coil current, and uq(k) is induced voltage (back voltage), which can be estimated by a state observer.(ii)
During commutation from phases “three and four” to phases “one and two,” the predicted current with a prediction horizon equal to 2 samples is(32)icoil(k+2)=A1C2A2[icoil(k)Uc(k)]+C2B2(-uq(k))+C1B1(Us-uq(k)).(iii)
Phase five is the phase in which the capacitor is charged, and it could be represented as(33)icoil(k+2)=-C2A2[A1icoil(k)+B1(Us-uq(k))Uc(k)]-C2B2(-uq(k)).The voltage of the capacitor, Uc(k), is measurable, and the voltage uq(k) is assumed to be constant in the prediction period.Figure 6
Working phases one and two (fed byUs).Figure 7
Working phases three and four (fed byUC).Figure 8
Phase five: capacitor charging phase.


### 6.1. Optimization and Current Control Technique

Generally, MPCs suffer from high computational complexity in online applications. In the presented case, there are theoretically 32 possible configurations of the IBGT switching states, but only6 of them are used in practice; most of them cannot be used because they generate short circuits, and others are practically redundant. If a global optimum is wanted in a “one-step” prediction horizon, then there are at most 5 possible voltage combinations. In general, for N-step prediction, there are 5N possible combinations. This problem is known as an NP-complete problem (the solution time grows exponentially with the problem size). Thus, in a general solution of the optimization problem, all 5N possible states must be taken into account, which means that the calculation time is expected to be relatively high. However, for short prediction horizons, using this method to determine the controller is still reasonable. In Figure 10, this situation is depicted. Three logical inputs ul1(t), ul2(t), ul3(t)∈{0,1} which identify the configurations are defined. Besides the general optimization solution described previously, a different optimization procedure, which is actually the well-known branch and bound method, was used. In this case 0-1 combinations through a binary tree are explored. The feasible region is partitioned in subdomains systematically, and valid upper and lower bounds are generated at different levels of the binary tree. The main advantage of the branch and bound method is that it can be interrupted at any intermediate step to obtain a suboptimal solution, although, in this case, the tracking performance deteriorates. To build a procedure for finding a local optimum, the following definition is introduced.Definition 2.
Given a voltage and inverter bridge configuration at timek that corresponds to the minimum of the function defined in (27) and its logical code, a “near permutation” at time k+1 is defined as “the combination code,” which identifies the two nearest possible voltages to the minimum at time k.If the “near permutation” at timek is considered, then the following procedure is proposed.Step 1.
Letul⋆ be an initial optimum of logic inputs that is found by a total search.Step 2.
Considerul⋆  and its “near permutations” in an Np prediction horizon branch phase.Step 3.
Check the cost function defined in (27) around ul⋆ to bound the branch.Step 4.
Find the minimum in the branch.Step 5.
Brancharound the new minimum candidate and check the cost function (27) to bound the branch.Step 6.
Go toStep4 until the specified time out.The procedure described previously is suitable for current reference values that do not contain steps or very fast changes. In these cases, it is possible to assume that the optimum in (27) changes slowly and that the computational advantages are known. At the first step control horizon, all the switching possibilities are considered to manage abrupt changes that may arise. The following step horizons only consider the “adjacent switch possibility.” In general, independently of the optimal search technique, Figure 11 shows the structure of the proposed controller. As shown previously [33], this type of approach could also be useful in nonlinear control systems. In fact, as proposed elsewhere [34], it is possible to describe a nonlinear system with this technique, which normally works around known working points, with a set of corresponding linear systems. In our presented case the set has 5N models. As already discussed, current control is particularly important for positioning a mechatronic system. It is also possible to build an MPC outer loop, as depicted in Figure 5. In other words, the desired current is calculated to minimize the quadratic error between the desired position and the predicted position for the current control loop.


## 6.1. Optimization and Current Control Technique

Generally, MPCs suffer from high computational complexity in online applications. In the presented case, there are theoretically 32 possible configurations of the IBGT switching states, but only6 of them are used in practice; most of them cannot be used because they generate short circuits, and others are practically redundant. If a global optimum is wanted in a “one-step” prediction horizon, then there are at most 5 possible voltage combinations. In general, for N-step prediction, there are 5N possible combinations. This problem is known as an NP-complete problem (the solution time grows exponentially with the problem size). Thus, in a general solution of the optimization problem, all 5N possible states must be taken into account, which means that the calculation time is expected to be relatively high. However, for short prediction horizons, using this method to determine the controller is still reasonable. In Figure 10, this situation is depicted. Three logical inputs ul1(t), ul2(t), ul3(t)∈{0,1} which identify the configurations are defined. Besides the general optimization solution described previously, a different optimization procedure, which is actually the well-known branch and bound method, was used. In this case 0-1 combinations through a binary tree are explored. The feasible region is partitioned in subdomains systematically, and valid upper and lower bounds are generated at different levels of the binary tree. The main advantage of the branch and bound method is that it can be interrupted at any intermediate step to obtain a suboptimal solution, although, in this case, the tracking performance deteriorates. To build a procedure for finding a local optimum, the following definition is introduced.Definition 2.
Given a voltage and inverter bridge configuration at timek that corresponds to the minimum of the function defined in (27) and its logical code, a “near permutation” at time k+1 is defined as “the combination code,” which identifies the two nearest possible voltages to the minimum at time k.If the “near permutation” at timek is considered, then the following procedure is proposed.Step 1.
Letul⋆ be an initial optimum of logic inputs that is found by a total search.Step 2.
Considerul⋆  and its “near permutations” in an Np prediction horizon branch phase.Step 3.
Check the cost function defined in (27) around ul⋆ to bound the branch.Step 4.
Find the minimum in the branch.Step 5.
Brancharound the new minimum candidate and check the cost function (27) to bound the branch.Step 6.
Go toStep4 until the specified time out.The procedure described previously is suitable for current reference values that do not contain steps or very fast changes. In these cases, it is possible to assume that the optimum in (27) changes slowly and that the computational advantages are known. At the first step control horizon, all the switching possibilities are considered to manage abrupt changes that may arise. The following step horizons only consider the “adjacent switch possibility.” In general, independently of the optimal search technique, Figure 11 shows the structure of the proposed controller. As shown previously [33], this type of approach could also be useful in nonlinear control systems. In fact, as proposed elsewhere [34], it is possible to describe a nonlinear system with this technique, which normally works around known working points, with a set of corresponding linear systems. In our presented case the set has 5N models. As already discussed, current control is particularly important for positioning a mechatronic system. It is also possible to build an MPC outer loop, as depicted in Figure 5. In other words, the desired current is calculated to minimize the quadratic error between the desired position and the predicted position for the current control loop.


## 7. Robust Stability Considerations and Simulation Results

In the presented case the performances are tested through numerical simulations. The numerical simulations are performed considering an extension of a sufficient condition stated in [35, Lemma 3.1]. If the systems defined previously are characterized by the following dynamic matrices, Ak, A1, and A2 are stable, and if for all matrices Qp and Qi as in (22) and (27) there exist matrices Pp, P1i, and P2i such that
(34)AkTPpAk-Pp=-Qp,A1TP1iA1-P1i=-Qi,A2TP2iA2-P2i=-Qi,
then one has the perturbed systems:
(35)x(k+1)=Akx(k)+mpAkx(k),x(k+1)=A1x(k)+m1A1x(k),x(k+1)=A2x(k)+m2A2x(k),
where mp, m1, and m2 are matrices with appropriate dimensions and are open loop stable if(36)∥mpAk∥≤min{λmin(Qp)2∥AkTPp∥+∥Pp∥,1},∥m1A1∥≤min{λmin(Qi)2∥A1TP1∥+∥P1∥,1},∥m1A2∥≤min{λmin(Qi)2∥A2TP2∥+∥P2∥,1}.
To stabilize the open loop structure robustly, Q(·)=I. As described in [35, Theorem 3.1] states that it is possible to stabilize the positional MPC in a closed loop. In fact, the theorem mentioned previously states a sufficient condition:
(37)∥(F1pTQpF1p+Rp)-1(F1pTQp)∥<min{12∥ATPp+Pp∥,1}.In our case, the control structure is a cascade of controllers. No stability conditions are given in the literature. The heuristic procedure to obtain the stability of the whole system, which was adopted in this paper, consists of finding the gain of the positional MPC according to relation (37) and testing the stability of the whole loop control structure. If the loop is not stable, then the value of matrix Qi is reduced, and the method is repeated. This trial-and-error technique causes the system to achieve stability. It is always possible to consider a number of models as in [34] such that the uncertainties of the system respect the sufficient condition stated in [35, Lemma 3.1]. To justify the presented approach, the following comparison results are reported that track the current control. These results provided some essential information regarding the presented research direction. To compare the performance of a possible control strategy, the second MPC was considered. The results shown in Figures 13, 15, and 16 are obtained through the control structure presented in Figures 11, 12, and 14, using a current sinusoidal test signal. In Table 1, the results are summarized. From these results, a majority of the approximations of the MPC are superior if a minimum number of models are considered. This set of models was used to determine the best value of the uq(k) voltage at the moving working point as proposed in [34]. A sinusoidal test function was chosen based on the fundamental harmonic of the current in a real working system. Moreover, to complete the comparison, a step test, the results of which are summarized in Figure 17, was done. Also, in this case, the superiority of the MPC seems to be evident, see Table 2. With both the hysteresis controller and the PI controller, it is not possible to recharge the capacitor, at least not with our proposed bridge structure. After determining the superiority of the MPC for our task [36], an MPC that follows a given current profile was presented, and in [37], a similar control strategy was shown but with a bridge with the capacitor recharging structure and a new control law. In [37], no trajectory tracking was analyzed. Nevertheless, PI controllers and hysteresis controller are very often used in these kinds of applications. Even though we must renounce to recharge the capacitor, PI and hysteresis controllers are very often used because of their simplicity to be realized. Based on the real pressure acting in a valve actuator, motion with an exponential disturbance with an initial value equal to 400 N is simulated, even though in the control law F0(t)=0. Special attention was paid to the applicability of the control system. In fact, the sample time was 40 μs. In this time, just a limited number of arithmetical operations can be performed, and for this reason, the simulation was performed with N=2 (prediction horizon).Table 1
Comparison of the results in the sinusoidal test current.
H
MPC
Hysteresis regulator
PI regulator
Energyerror(∫e2dt)
0.23% (0.39%)
0.35%
1.25%
Calculationcomplexity
30 (5) models
—
—Table 2
Comparison of the results in the step test current.
H
MPC
Hysteresis regulator
PI regulator
Energyerror(∫e2dt)
24% (42%)
40%
45%
Calculationcomplexity
30 (5) models
—
—


### 7.1. Capacitor Recharging

Along with the tracking task, another interesting aspect is the following: a self-recharging voltage capacitor control. A capacitor can be charged in working phase 5, where the motor inductance either allows the charging current or the mechanical braking energy can be fed back. Thus, appropriately partitioning the working phases can make both tracking and recharging possible. The balance between these two tasks is achieved by the cost function defined in (27). Depending on the motion cycle dynamic, the weights Qj and Rj were set to maintain the level of the voltage capacitor within an acceptable tolerance. Thus, in particular ΔUc(k)=(Uc(k)-Ucd*), where Uc(k) is the measured voltage of the capacitor and Ucd* is its desired level, which should be kept within a tolerance band. The capacitor voltage is used to enable high dynamic changes in the actuator current, which in turn cause large variations in the capacitor voltage. Thus, the system could not be operated continuously if the capacitor control was not included in the optimization procedure. Furthermore, efficiency was improved by properly utilizing the mechanical braking energy. Simulation results were achieved using real data, and they show promising results. In this phase, some different control strategies were tested using real data from the already-built valve. The results encourage proceeding in the prototyping phase, in which the rest of the structure, which basically consists of a DSP, must have a power bridge according to the simulated data. Some interesting aspects are worth identifying. The lower part of Figure 18 shows how the self-charging phase can start during the following phase. In this case, an initial error occurs as long as the capacitor does not achieve the full voltage. In the presented simulated case, the initial voltage value in the capacitor is 48 V. It is possible to start from 0 V; in this case the charging phase lasts 6 cycles more, if the same weighting matrices are considered in the cost function. To achieve a higher charging velocity it is necessary to increase the matrix Rj in the index (27). Finally, it is interesting to note that through MPC, it is possible to obtain a self-balanced switch of the IBGT. In fact Figure 9 shows two possible freewheeling circuits. It is acceptable to balance the switching phase on these two circuits to obtain a balance frequency. A compromise between good tracking and stability of the supply voltage is achieved. In particular, good capacitor voltage stability provides the necessary condition for good follow-up results; it is always necessary to identify the trade-off between capacitor voltage stability and followup.Figure 9
Working phase six (freewheeling).Figure 10
Three of the combinations in two-step prediction.Figure 11
Scheme of the multimodel predictive control.(a) Control scheme using PI controller in the loop. (b) PWM generation.
(a)
(b)Figure 13
Current tracking obtained with the control scheme of Figure12.Figure 14
Control scheme using hysteresis controller in the loop.Figure 15
Current tracking obtained with the control scheme of Figure14.Figure 16
Current tracking obtained using an MPC (2-step prediction) with the control scheme of Figure11.Figure 17
Step test: current tracking obtained using MPC (2-step prediction) and PI controller.From the top: three cycles for position, current tracking, and capacity charging using MPC.
(a)
(b)
(c)


## 7.1. Capacitor Recharging

Along with the tracking task, another interesting aspect is the following: a self-recharging voltage capacitor control. A capacitor can be charged in working phase 5, where the motor inductance either allows the charging current or the mechanical braking energy can be fed back. Thus, appropriately partitioning the working phases can make both tracking and recharging possible. The balance between these two tasks is achieved by the cost function defined in (27). Depending on the motion cycle dynamic, the weights Qj and Rj were set to maintain the level of the voltage capacitor within an acceptable tolerance. Thus, in particular ΔUc(k)=(Uc(k)-Ucd*), where Uc(k) is the measured voltage of the capacitor and Ucd* is its desired level, which should be kept within a tolerance band. The capacitor voltage is used to enable high dynamic changes in the actuator current, which in turn cause large variations in the capacitor voltage. Thus, the system could not be operated continuously if the capacitor control was not included in the optimization procedure. Furthermore, efficiency was improved by properly utilizing the mechanical braking energy. Simulation results were achieved using real data, and they show promising results. In this phase, some different control strategies were tested using real data from the already-built valve. The results encourage proceeding in the prototyping phase, in which the rest of the structure, which basically consists of a DSP, must have a power bridge according to the simulated data. Some interesting aspects are worth identifying. The lower part of Figure 18 shows how the self-charging phase can start during the following phase. In this case, an initial error occurs as long as the capacitor does not achieve the full voltage. In the presented simulated case, the initial voltage value in the capacitor is 48 V. It is possible to start from 0 V; in this case the charging phase lasts 6 cycles more, if the same weighting matrices are considered in the cost function. To achieve a higher charging velocity it is necessary to increase the matrix Rj in the index (27). Finally, it is interesting to note that through MPC, it is possible to obtain a self-balanced switch of the IBGT. In fact Figure 9 shows two possible freewheeling circuits. It is acceptable to balance the switching phase on these two circuits to obtain a balance frequency. A compromise between good tracking and stability of the supply voltage is achieved. In particular, good capacitor voltage stability provides the necessary condition for good follow-up results; it is always necessary to identify the trade-off between capacitor voltage stability and followup.Figure 9
Working phase six (freewheeling).Figure 10
Three of the combinations in two-step prediction.Figure 11
Scheme of the multimodel predictive control.(a) Control scheme using PI controller in the loop. (b) PWM generation.
(a)
(b)Figure 13
Current tracking obtained with the control scheme of Figure12.Figure 14
Control scheme using hysteresis controller in the loop.Figure 15
Current tracking obtained with the control scheme of Figure14.Figure 16
Current tracking obtained using an MPC (2-step prediction) with the control scheme of Figure11.Figure 17
Step test: current tracking obtained using MPC (2-step prediction) and PI controller.From the top: three cycles for position, current tracking, and capacity charging using MPC.
(a)
(b)
(c)


## 8. Conclusions and Outlook

The paper implements a multilevel inverter control system by integrating flatness and two cascaded MPCs. The first MPC generates an optimal desired current by minimizing the position error. The second MPC considers this optimal desired current as a reference signal to find an optimal switching law for the proposed multilevel inverter. The current MPC considers energy storage optimization using a capacitor in the proposed inverters. Finally, simulations using real data are shown.



---

*Source: 750190-2013-03-13.xml*