1. Introduction
Nowadays, power converters application is attracting attention in academy due to the global energy demand growth and the search for new types of renewable energy. In a classical distribution system, the consumers are passive loads; however, due to technological improvement and advances in power electronics and power generation, by using renewable sources, these consumers can generate their own energy. They might operate connected on-grid or islanded, so a set of renewable sources and loads that are possible to operate islanded are called microgrids. The DC form is the most common type of power generation by renewable sources, DC systems are focused on several investigations. Thus, the investigation of the DC-DC power converters applications has a wide area of studies and functions, such as photovoltaic systems [
1], energy management, microgrids [
2,
3,
4,
5], electric vehicles, and others applications.
Meanwhile, many studies focused on stability and performance improvement, so several control strategies applied on DC-DC power converter were studied, such as: robust control [
4,
5,
6,
7,
8,
9]; sliding mode control [
10,
11]; adaptive control [
12]; model predictive control [
13]; fault-tolerant control [
14,
15]; classical pole-zero cancellation [
16]; designing the controller in frequency domain [
17]; the fractional-order PID controllers [
18,
19]; and lastly, based on the GPI Observer [
20].
In [
4,
5], robust control methodologies were applied on the cascade DC-DC buck converter, to assess the oscillations in the feeder converter caused by the power variation of the constant power load (CPL). Additionally, the control design based on Kharitonov’s theorem considered the linear parameter varying (LPV) model of the DC-DC converter taking into account an uncertainty family in the type of hyper box; the control strategy was designed in indirect form, in other words, it was designed in s-domain, in sequence the Tustin conversion was used to digital implementation of the controller.
In [
6], the sampled data output feedback control applied to the DC-DC buck converter was investigated. In sequence, a reduced-order observer and a robust output feedback controller were designed to ensure the robustness and stability of the system. In [
7], a regulation of DC-DC boost converter considering an uncertain family was investigated. The LPV model was developed when considering uncertainties, such as load variation, non-modeled dynamics, and input voltage variation. Subsequently, the proposed methodology merges the cascade controllers and nested in a reduced-order PI observers design, in order to maintain the desirable voltage regulation performance.
Two control methodologies were presented in [
8] for mismatched DC-DC buck converters: the first method used a multiple surface sliding mode control to handle mismatched load uncertainty and the other method is based on simultaneous state and disturbance and observer. In addition, the sensorless structures combined with the second-order sliding mode control is a common strategy developed in several reported works [
21,
22].
In [
18], a fractional sliding mode control (FRSMC) applied in a practical buck converter system based on PI and PID structures was investigated; besides that, the voltage setpoint variation was performed. In [
19], the effectiveness of the fractional order PID (FOPID) controller applied in a DC-DC buck converter feeding a DC motor was investigated, aiming to compare the performance and robustness with the classical control approaches.
A cascade control structure with internal sliding mode control was applied to the multiphase synchronous DC-DC converter [
10]. In the outer voltage loop, a feedforward compensation was designed. Furthermore, disturbance observers in both the inner and outer loop were used in order to compensate the difference between the model and rigged system developed.
According to Chan et al. [
11], adaptive current-mode control for DC-DC boost converter is given; therefore, the adaptive law is formed combining the existing current-mode control law with adaptive law that generates the inverse of the load resistance. In [
12], a predictive control that is based on both digital average voltage and current (DAV/DAC) for switching the DC-DC converter is proposed, aiming to ensure the stability of the system, besides enabling the system exhibit high output voltage accuracy and overcurrent protection.
A re-configurable structure of the resonant control applied on the DC-DC converter was investigated [
13]. Besides that, it was used as a soft-switching technique to acquire high performance. Thus, another resonant element was introduced into the LLC converter. In [
14], an analysis of fault-tolerant operation capabilities of an isolated bidirectional full-bridge current source DC-DC converter was investigated. This study regarded the system reconfiguration in another simplified structure.
DC microgrids have been widely used for critical applications, several studies in control use them to enhance the performance and stability of these systems investigated in recent years. In [
17], a novel decentralized output constrained control algorithm to ensure the high performance of the single bus DC microgrid is presented; in addition, the goals of this control technique are: high performance of the voltage control of the DC bus, user-defined the load sharing, and minimize the circulating current. Meanwhile, in [
23] a unified distributed control strategy applied in a stand alone DC microgrid was proposed. This control strategy integrates three types of control solution for DC microgrid, such as: power control, voltage control, and droop control. Thus, to design the controller first formulate a solution of the optimal power flow to stand-alone DC microgrid, in sequence the dynamic solving algorithm based on primal-dual decomposition method is proposed, the convergence and optimal conditions of its equilibrium point were proved. Six DC microgrid systems based on the DC microgrid benchmark are used to evaluate the effectiveness of the proposed methodology.
In [
24], a novel control structure applied in a hybrid DC/AC microgrid is investigated. This technique integrates concepts of the Lyapunov theory and input-output feedback linearization (IOFL) to design the controller and apply in a specific microgrid that take into account several energy sources, i.e. wind genearation, photovoltaic generation, besides that they use a storage device composed by lithium battery and buck-boot converter, each energy source and device has a control subsystem that is designed by using the proposed technique to ensure the stability of the system. Several simulation tests were performed in order to evaluate the effectiveness of the proposed technique.
A recent literature review presents a lack of contribution in the direct form to design digital controllers as well as perform the experimental investigations of the robust control methodologies applied in the DC-DC converters or in DC microgrids. Therefore, this article proposes a direct form to design the digital RST robust control (i.e., name given because of the R, S, and T polynomials computed) when considering the set of parametric interval uncertainties, i.e., the family of uncertainties is a hyper box representation to work with the LPV model of the DC-DC buck converter. This control methodology is applied on a DC-DC buck converter while taking a parametric variation, such as load resistance variation and voltage reference setpoint variation into account; furthermore, the proposed methodology uses the Chebyshev sphere theorem as constraints of the optimization problem to compute the robust controller gains aiming to guarantee the performance and stability of the system. The proposed digital RST robust controller is designed to provide robust performance and robust stability when the system is subject to an uncertain family. Summarizing, this study proposes the following main contributions:
a new direct form to design the digital robust control by using convex optimization based on the Chebyshev sphere circumscribed in the Kharitonov’s rectangle, with the purpose to relax the space of solutions and ensure the robust performance and robust stability; and,
simulations and experiments were performed to demonstrate the effectiveness of the proposed methodology. The results showed that the proposed approach outperforms a classical control approach as well as the robust control based on Kharitonov’s rectangle.
The remaining of this paper is organized, as follows:
Section 2 describes the mathematical model of the DC-DC buck converter as well as the operating conditions;
Section 3 discuss the digital robust control based on convex optimization constrained to a Chebyshev sphere design as well as presents a robust control methodology based on Kharitonov rectangle;
Section 4 performs the proposed methodology, shows the board system developed, besides that presents a brief description of the tests performed;
Section 5 discusses results performed by the performed tests; and finally,
Section 6 concludes this study.
2. Mathematical Model of DC-DC Buck Power Converter
Figure 1 shows a well-known DC-DC buck converter topology, when considering a passive LC filter. This system is divided into three parts: hardware, microcontroller, and computer. The hardware part comprises all power elements (green box), the microcontroller part (blue box), and the computer part corresponds to graphical user interface (GUI) in order to have a user-friendly operation (red box).
It is necessary to consider how the power converter operation mode functions in order to obtain the converter mathematical model, so, in this paper, the continuous conduction mode (CCM) is adopted, in other words, the current that passes in the inductor element is not null during any time of the switching period [
25].
Figure 2 depicts the equivalent circuit of each stage of the switching element.
Figure 3 illustrates the system’s waveform of all electrical variables, when the DC-DC buck power converter operates in CCM mode; therefore, considering that the condition above presented, the following equations can be obtained to each switch condition, i.e., switch on (Equation (
1)) and off (Equation (
2)), respectively.
When the switching element (
) is activated, the diode element (
) is disabled; therefore, the source feed directs the resistance load and charges the capacitive and inductive elements. On the other hand, when the switching element (
) is not activated, the source does not feed the load; therefore, the diode element conducts and the active components (inductive and capacitive elements) discharge to supply the load.
Figure 3 represents these effects during the switching period [
25].
Re-writing the (
1) and (
2) to easy represents in matrix form the system dynamics, as follows.
Subsequently, adopting the average state-space model (ASM) that considers (
3) and (
4), it is possible write the following non-linear system.
When considering that the input source
is fixed, and adopting the controlled duty cycle
as input, thus
,
and
in (
5) are decomposed into a fixed term and a time-varying term:
where, the operation duty cycle
is chosen by
, where
is the desired output voltage. In addition,
and
denote, respectively, the current in the inductor and the voltage in the capacitor. Thus, the linear state space model (SSM) around the operational point (OP)
and its transfer function
are presented below.
3. Digital Robust Rst Control Design by Using Convex Optimization
An essential goal of control systems theory is to guarantee the stability and desired performance, despite the uncertainties due to parametric variations and disturbances. However, classical control design methodologies are not able to handle model uncertainties, forcing the use of local models around an operational point, which represents a nominal condition.
The uncertainties can be classified into two categories: nonstructural uncertainties (nonparametric uncertainties) and structured uncertainties (parametric uncertainties). In this study, control techniques for mathematical models considering parametric uncertainties are investigated; in addition, the extension to discrete domain the main theorems of robust control aiming to design the robust controller in a direct form by using RST structure was considered [
26,
27,
28,
29,
30].
Several control design methodologies that consider a family of uncertainties uses the Kharitonov theorem to guarantee polynomial stability [
26,
27,
28]. In view of the interval uncertainties, i.e., hyperbox, the polynomial stability can be evaluated by testing four Kharitonov’s polynomials, and when considering the following interval variation of each parameter obtained by a polynomial. Subsequently, the uncertainty region can be delimited by (
10), where the superscripts + and − denote the upper and lower bound limit of each uncertainty interval, respectively.
According to Bhattacharyya et al. [
27], when assuming that the family of the polynomial has a constant degree, it is possible to ensure the Schur’s stability of these uncertainties, they must be computed for the Kharitonov polynomials, as following [
26,
27].
According to Kharitonov’s theorem, the family of polynomials (
11) is robustly stable if and only if all Kharitonov’s polynomials are Schur stable, for more details about this property see [
26,
27,
29,
30].
3.1. Robust Controller Design via Interval Pole Placement Based on Kharitonov’s Theorem
To design the controller, a region of uncertainties must be defined based on the interval parameter variation of the plant model. According to [
26,
27], the controller can be designed based on the Kharitonov’s theorem applying linear programming to a set of linear inequality constraints. It was considering an
n-th order of the strictly proper uncertain plant
(i.e., the order of numerator is lesser the denominator, in other words,
) and that presents an
-th and
-th the controller order that has the following polynomial structure
,
, respectively, and
can be computed by the sum of all coefficients of the
polynomial, in other words, it is only a static gain, due to the polynomial
presents an integrator, when considering that the uncertain plant is the 0-type (i.e., it does not have an integrator), another way to compute the
can be used the Internal Model Control (IMC) methodology to ensure that the system performs the desired behavior, such as that proposed by [
31]. Furthermore, it is essential to define the desired constraints that allow the choice of the control structure in order to consider this control structure and order of the controller.
Thus, (
12) and (
13) can be computed as a polynomial closed loop, as follows.
Adopting the hyperbox region to represent a desired closed-loop polynomial containing a required performance, this region can serve each coefficient of the desired polynomial as a range depending on uncertain parameters, as described in (
17).
When considering that closed-loop polynomial (cf. (
14) and (
15)) is inside of the region defined by the desired polynomial region (cf. (
16)), the problem can be rewritten as a set of matrix inequalities, as follows.
where,
Thus, the problem of pole placement can be rewritten as a linear optimization problem and any linear programming methodology can be used to solve it. According to [
26,
27], any cost function can be chosen; however, a proper cost function is the sum of controller parameters, aiming at reducing controller gains and control effort. Hence, the problem presented in (
17) can be rewritten as a linear optimization problem, as presented in (
18).
3.2. Robust Design Methodology Based on Chebyshev Sphere
The optimization problem presented in the previous Subsection can be rewritten based on the Chebyshev sphere to represent the set of inequalities that describes the plant uncertainties. The Chebyshev theorem describes that it is possible to determine the largest compact set
of center and maximum radius, whose Euclidian norm, which is inscribed in a polytope
, is described by the set of the linear inequalities constraints [
32]. Hence, the ball set can be observed in
Figure 4, which the ball is inscribed in convex set
. Therefore, the set of solutions is more relaxed when compared with the one that is provided by Kharitonov’s rectangle.
To rewrite this optimization problem, suppose that
is defined by a set of convex inequalities, i.e.,
. If
it can be found an
(center of Chebyshev sphere) that solve the linear optimization problem, according to the following relations:
where,
where,
is the norm of coefficients of matrix A.
Figure 5 provides a simplified flowchart of the design of the digital robust RST controller methodology for robust pole placement based on Chebyshev sphere constraints in the linear optimization problem (cf.
Figure 5, (
19) and (
20)). The main difference between the Chebyshev convex optimization and the classical form of this problem is the conservativeness that is reduced for a methodology that is based on Chebyshev convex optimization, such that all the solutions to the compact set defined is possible [
27,
32]. Furthermore, this region ensures that the system obeys the constraints set defined, resulting in a less conservative solution.
Figure 5 presents the flowchart of the proposed methodology for direct form to design the digital RST robust control by using the Chebyshev sphere as constraints to solve this problem.
The first step is to obtain the generic model of the plant, and the control structure (cf. showed in (
12) and (
13)), as shown in
Figure 5. In the second step, the uncertainty region must be described. The uncertainty region is made of intervals of each plant parameters, with each coefficient of the polynomial being a range of the values (cf. presented in (
9) and (
10)). In the third step, the desired performance region is defined (cf. presented (
16)) and the closed-loop polynomial is presented to evidence the dependence between each coefficient and controller parameters (cf. (
14) and (
15)). In the fourth step, the optimization problem is rewritten in order to maximize the radius of the Chebyshev sphere [
32]; the constraints of the optimization problem are the same as those that are presented by Bhattacharyya et al. [
26,
27]. Subsequently, the solution must be tested; if the solution is not feasible, it is necessary to relax the constraints, that is, relaxing the desired performance to augment the uncertainty region. However, if the solution is feasible, a sixth step must be used, where it is written as the controller model using the feasible solution.
4. Methodology
In this section, the procedure to apply the proposed methodology is discussed, all controller gains are presented, and the adopted topology is described. Furthermore, the experimental environment is presented, and a brief description of the developed tests is given.
4.1. Presentation of the Experimental Environment
To perform the experiments, a DC-DC Buck converter test board was developed, c.f. shown in
Figure 6.
Table 1 presents the main elements that compound the experimental system as well as shows the main values of the parameters.
When considering the data presented in
Table 1, it is possible to obtain the mathematical model of the buck converter (cf. Equation (
8)), in addition, was considered the digital conversion by using the ZOH method and use the sampling frequency to discretize. The following equation presents the math model of the system in discrete domain.
When considering the uncertainty range of the load variation, the ZOH method was used to discretize the plant model, so the parametric interval transfer function is obtained, as follows.
To design the controller, the following requirements are desired, such as settling time of 0.5 s and considering the maximum overshoot less than 15%. This condition is validated for classical pole placement; however, when the system is subject to a specific family uncertainty, this point of performance can be unfeasible. Therefore, to design the digital RST robust control for all methodologies investigated in this study, the following desired region was adopted: settling time-varying in the range 0.1–0.9 s and the maximum overshoot varying in the range 10–30%. Therefore, it is possible to determine the desired performance region (cf. presented in (
16)). The control structure adopted can be displayed in
Figure 7, which shows a generic block diagram to develop the system as well as the control signal can be computed by using the (
23) as well as the structure adopted for each control polynomial. The controller gain values for each methodology are presented in
Table 2. To compute the controller gains, an automatic computational routine is developed and executed by using the computer with the following specifications: Intel core i7-3770 processor (4 cores, 3.4–3.9 GHz of clock) and 16 Gb of RAM. Using that computer, the control design routine takes less than 6 s to compute the gains.
4.2. A Brief Description of the Experiments
Several tests were performed in order to evaluate the performance of the controllers designed; all the tests developed presented numerical results (simulation) and experimental results (i.e., obtained by the board system developed). Three types of tests were prepared, as follows:
Test 1—the first test focused on voltage setpoint variation. To perform this test, the system starts with the reference voltage of 5.0 V; at every 2 s, the pulse voltage setpoint is varied for 2 s, the first variation is the 1 V, then the pulse setpoint is increased 1 V in variation until 3 V beyond the setpoint was reached. In sequence, the same procedure was performed to decrease the pulse setpoint variation.
Test 2—the second test aimed to evaluate the performance of the digital RST controller when the system was subject to a load resistance variation. First, the DC-DC buck converter was tuned at the operational point of 5.0 V. Additionally, the load resistance value was 4.0 . Then, at each 1.0 s, the load was varied by using the following values of load resistance variation around the nominal value of load resistance . For each variation, the oscillation in the capacitor voltage was observed and an evaluation was done to determine which methodology outperforms the others approaches when the system was subject to parametric variation.
Test 3—the last test had the goal of quantifying the performance approaches by using integral indices. For this test, the integral square error (ISE) index was computed in order to assess the performance of the output signal and integral index ISCS to assess the energy of the control effort.