Next Article in Journal
Analysis of Primary Energy Factors from Photovoltaic Systems for a Nearly Zero Energy Building (NZEB): A Case Study in Lithuania
Previous Article in Journal
A Review of Battery Technology in CubeSats and Small Satellite Solutions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Different Features of Control Systems for Single-Phase Voltage Source Inverters

by
Zbigniew Rymarski
* and
Krzysztof Bernacki
Department of Electronics, Electrical Engineering and Microelectronics, Faculty of Automatic Control, Electronics and Computer Science, Silesian University of Technology, Akademicka 16, 44-100 Gliwice, Poland
*
Author to whom correspondence should be addressed.
Energies 2020, 13(16), 4100; https://doi.org/10.3390/en13164100
Submission received: 23 July 2020 / Accepted: 4 August 2020 / Published: 7 August 2020
(This article belongs to the Section L: Energy Sources)

Abstract

:
The aim of this paper is to determine which type of control loop is better for each particular type of standard voltage source inverter (VSI) load (static, dynamic, nonlinear rectifier resistive-capacitive RC load). A comparison of three different types of controllers for single-phase VSIs is presented. The first two are of the single input single output (SISO)/proportional-integral-derivative controller (PID) and coefficient diagram method (CDM) types, and the third is of the multi-input single output (MISO)/passivity-based control (PBC) type. The selections of the gains for SISO and MISO controllers are presented, including the problem in the choice of PBC controller gains caused by the imperfection of the pulse width modulation (PWM) modulator. For a standard nonlinear rectifier RC load, the new control quality factor (CQF) is defined to distinguish the properties of the controllers. The results show the superiority of the MISO-PBC controller for the RC load; however, for a linear dynamic load, the MISO-PBC performed worse than the SISO-PID or CDM. Therefore, the choice between either a simple and inexpensive SISO with one measurement trace or an expensive MISO with three measurement traces depends on the purpose of the VSI. The initial simulations and measurements of the experimental models are presented and discussed.

Graphical Abstract

1. Introduction

There are many different types of control systems that are used in voltage source inverters (VSIs) for uninterruptible power supply (UPS) systems. A serious problem is selecting the best control for dedicated applications. The solution can be a single input single output (SISO) system, in which only the output voltage is measured and controlled, while the output current is treated as an independent disturbance [1,2,3]. However, the output voltage has an influence on the output current, which creates an additional feedback loop. It is possible to calculate the shift of the poles of the characteristic equation of a closed-loop system when we consider this additional feedback loop and determine the final stability of an actual system [4,5]. The distortion of the input DC voltage of the inverter that results from operating different types of DC/DC converters (e.g., impedance networks that cooperate with photovoltaic (PV) modules) should also be considered [6,7]. The main advantage of an SISO solution is its simplicity and low cost because only one trace measuring output voltage with galvanic isolation is required. There have been some improvements to SISO using double control loops—a fast inner loop, e.g., a proportional-integral-derivative controller PID [8,9] or the coefficient diagram method (CDM) [10,11,12], and an outer loop that is only used to damp harmonic distortions, e.g., a repetitive controller (RPC) [13,14,15,16,17,18]. A repetitive controller is a discrete-time harmonic generator [14] that is plugged into the outer feedback loop, which works perfectly in the steady state with a standard nonlinear RC load [19] while minimizing the output voltage static error (it damps the fundamental harmonic). However, it has one feature that can be a disadvantage for dynamic load changes: it remembers a disturbance from the previous fundamental cycle (the basic RPC has a register in which all of the previous fundamental cycle’s values are stored) and even after the output disturbance vanishes, it still tries to damp it by distorting the output voltage. Therefore, a solution with an RPC cannot be compared with the fast multi-input instantaneous controllers. The inner feedback loop cooperates with the RPC best in the outer loop if its magnitude Bode plot is almost flat up to the Nyquist frequency.
A much more sophisticated system is a multi-input single output (MISO) system (Figure 1), in which the output current, the output filter inductor current, and the output voltage are measured and all of them are used as the input variables of the output voltage controller. In this case, the output current is also treated as a disturbance because the state-space matrix of the VSI does not depend on the particular load. However, in the case of MISO, we measure the output current (treated as a space variable or a measured disturbance [20]), which renders the problem of the missed feedback loop (output voltage—output current) unimportant.
The assignment of a linear or nonlinear plant model is the basic problem in control loop design. The inverter plant has some nonlinearity. The dead time Tdt of switching transistors in the H-bridge is one of the reasons for distortions of the output voltage [4] (a little step decrease in output voltage occurs when the inductor current changes direction) and can be taken into account in the nonlinear model of the inverter [21]. However, for the standard Tdt ≤ 0.5 µs (for Metal-Oxide Semiconductor Field-Effect Transistors MOSFET, it can be tens of ns) and for the switching frequency fs = 20 kHz (switching cycle Ts = 50 µs), the decrease in the fundamental harmonic of the output voltage is approximately only Tdt/Ts ≤ 1%. The problem of the influence of the nonlinear characteristic of the inductor in the inverter output filter on the design of the adaptive control loop (with calculation of the nonlinear inductance characteristic) is presented in [22]. However, the change in the coil’s inductance and its equivalent serial resistance as a function of the amplitude and the frequency of its magnetizing current depends on the core material. For contemporary alloy-powder materials (e.g., Sendust (MS)/Super-MSS™ [23]), the exemplary change of inductance is about 5%, and the equivalent serial resistance is low in the operating point [24]. Therefore, it has no serious impact on the VSI model. The solution of the state-space equations of the pulse width modulation (PWM)-controlled VSI results in the nonlinear (exponential) dependency of the state-space variables on the switching on (duty) time [25]. In [26], the nonlinear control function is approximated by the Fourier series. The nonlinear modeling of the inverter using Hammerstein’s approach (by means of the black-box identification method) is presented in [27]. All the approaches that take into account the nonlinearity of the inverter slightly increase the accuracy of its modeling; however, it has been shown that a linear approximation of VSI results in quite a satisfactory compatibility of the linear theory and the experimental model measurements and enables the easily understandable design of the control loop. Modeling the inverter as linear was practiced for many years (e.g., [1]) and resulted in the control design keeping low distortions of the output voltage for standard loads. We do not consider the thermal dependencies of the inductance, assuming that we work in standard ambient temperature conditions and that self-heating results in a stabilized temperature of the coil.
Estimating the quality of the VSI output voltage is also a challenge. The most common method for estimating static load distortions of the output voltage is to use the total harmonic distortion (THD) coefficient, which is simply the ratio of the root mean square (RMS) value of a set of the higher than the first harmonics of the output voltage to the RMS value of the fundamental harmonic. The problem is that this is an averaged value and the spectrum of harmonics cannot be discussed using it. The distortions of the VSI output voltage after dynamic load changes [19] are defined by means of the over- or undershoots of the output voltage and the settling time. There are some approaches that use wavelet analysis to present more complex distortions of the output voltage [28]. In systems with a discrete-time pulse width modulation (PWM) control, in which the delay of some (at least one) switching cycles is unavoidable, the output voltage over- or undershoots are mainly the result of an insufficient capacitance in the VSI output filter (the feedback loop is important but to a lesser degree) because the system works in at least one switching period after the step decrease/increase in the load without control. For a sufficiently high output filter capacitance value, all of the properly designed SISO and MISO controllers reduced the distortions of the output voltage for the dynamic load to a similar degree. However, the settling time was different. Therefore, it seems that the most significant comparison of controller features can be done using a nonlinear RC load. A comparison of the distorted input inverter DC voltage in the case of the cooperation of a VSI with an impedance network for the discontinuous current mode (DCM) of a network was previously presented in [7]. Because the THD or the weighted THD (WTHD) [29] coefficients only show the averaged values, we will now analyze the spectrum of the harmonics of a single-phase VSI with selected control systems and attempt to determine which spectrum is the most appropriate for controlling a VSI. The newly defined control quality factor (CQF) will be used in order to permit a clear estimation of the SISO and MISO controller properties (Figure 1). Section 2 presents an analysis of the frequency range, in which the features of the controllers will be compared when considering the magnitude of the VSI output filter. Section 3 presents the approach to the discrete-time modeling of the VSI. Section 4 presents the exemplary SISO controllers of the VSI, and the MISO passivity-based control (PBC) PBC-IPBC2 is described in Section 5. Special attention was paid to selecting the PBC gains that were used in the simulations. Section 6 is devoted to a comparison of the SISO and MISO controllers for standard loads using an experimental VSI model. The main aim of the paper is to determine which type of control loop is more suitable for each particular type of standard load (static, dynamic, nonlinear rectifier RC load).

2. Analysis of the Frequency Range of the Controller Features

The transfer function of a simple continuous SISO model (Figure 1) of VSI (Figure 2) is calculated in Equation (1). Let us assume that in the steady state VOUT = VREF and kV = 1.
V O U T ( s ) = G ( s ) exp ( s T s ) F L C ( s ) 1 + k v F l f ( s ) G ( s ) exp ( s T s ) F L C ( s ) V R E F ( s ) + M F L C ( s ) 1 + k v F l f ( s ) G ( s ) exp ( s T s ) F L C ( s ) Δ V D C ( s ) + + Z O U T ( s ) 1 + k v F l f ( s ) G ( s ) exp ( s T s ) F L C ( s ) I O U T ( s )
In the simplified continuous model (Figure 2), for Ts << Tm (switching period Ts, fundamental period Tm), we can omit the delay, exp(−sTs), of the PWM modulator (Figure 1). The transfer function of the changes ΔVDC in the VDC voltage supplying the VSI can decrease the result of the distortion of the DC voltage from the impedance network [6,7]. The output impedance of the VSI without a feedback loop ZOUT is calculated in Equation (2), where RLF is an equivalent serial parasitic resistance of the VSI bridge and the filter choke. In further calculations, we will assign kV = 1.
Z O U T ( j ω ) = R L F + j ω L F ( 1 ω 2 L F C F ) + j ω C F R L F
We neglect the equivalent serial resistance (ESR) RCF of the capacitor CF. There is a problem with estimating the damping coefficient ζF when a load current is modeled as an independent current source. Therefore, the damping coefficient ζF will be calculated in Equation (3) for a case in which the nominal load resistance RLOAD is present in the transfer function of the VSI. For RLF << RLOAD, the output filter transfer function is calculated in Equation (3),
F L C ( j ω ) 1 ( 1 ω 2 / ω F 0 2 ) + j 2 ξ F ω 2 / ω F 0 2 ,   where   ω F 0 = 1 L F C F ,   ξ F = 1 2 ( R L F C F L F + 1 R L O A D L F C F )
For the low-frequency range below ωF0 (the resonant frequency of the output filter LFCF), the output impedance of the VSIs without feedback loop ZOUT, the transfer functions FLC, and Flf (Figure 2) are calculated in Equation (4),
Z O U T ( j ω ) | ω < < ω F 0 R L F + j ω L F ,   F L C ( j ω ) | ω < < ω F 0 1 ,   F l f ( j ω ) | ω < < ω F 0 1
The transfer function of the output current, which is treated as an independent disturbance for the low-frequency range, lower than ωF0 is calculated in Equation (5),
V O U T ( j ω ) I O U T ( j ω ) | ω < < ω F 0 ( R L F + j ω L F ) G ( j ω )
The transfer function of DC voltage change ΔVDC for the low-frequency range ω << ωF0 and kV = 1 is calculated in Equation (6),
V O U T ( j ω ) Δ V D C ( j ω ) | ω < < ω F 0 M G ( j ω )
The specific case shown in Equation (7) is for ω = ωF0, and the values of LF = 1 mH, CF = 50 μF, and RLF < 1 were used.
Z O U T ( j ω ) | ω = ω F 0 1 R L F L F C F , F L C ( j ω ) | ω = ω F 0 1 j 2 ξ F ,   V O U T ( j ω ) I O U T ( j ω ) | ω = ω F 0 1 R L F L F C F 1 j G ( j ω F 0 ) / ( 2 ξ F ) 1
The transfer function of ΔVDC for the low-frequency range ω = ωF0 and kV = 1 is calculated in Equation (8),
V O U T ( j ω ) Δ V D C ( j ω ) | ω = ω F 0 M j 2 ξ F + G ( j ω F 0 )
For the high-frequency range over ωF0, the output impedance ZOUT of the VSI without the feedback loop and the transfer function of the output current (which is treated as an independent disturbance) for the high-frequency range over ωF0 is calculated in Equation (9). We assume that filter FLC() for ω >> ωF0 damps the signals of whatever is the transfer function Flf() (|Flf()| ≤ 1).
Z O U T ( j ω ) | ω > > ω F 0 1 j ω C F ,   F L C ( j ω ) | ω > > ω F 0 ω F 0 2 ω 2 , V O U T ( j ω ) I O U T ( j ω ) | ω > > ω F 0 = 1 j ω C F .  
The transfer function of ΔVDC for the high-frequency range ω >> ωF0 and kV = 1 is calculated in Equation (10),
V O U T ( j ω ) Δ V D C ( j ω ) | ω > > ω F 0 M ω F 0 2 ω 2
Equations (9) and (10) mean that for ω >> ωF0, the output impedance of the VSI with the feedback loop and the transfer function of ΔVDC with the feedback loop are not dependent on it.
It can be seen from the presented Equations (1)–(10) that the gain of the controller is important in damping the disturbances for the resonant frequency ωF0 and below this frequency of the output filter. For the higher frequencies, it is not important because of the high suppression of the output filter. This discussion, which is based on Figure 2, led to the further investigation of controller properties being carried out in the 2ωF0 frequency range.

3. Discrete-Time Modeling of the Voltage Source Inverter

Today, only the digital control of VSI is important. Every PWM discrete-time modulator implements one switching period delay Ts (the data are written for the modulator in one switching cycle, and the pulse width is set in the next cycle). For VSIs that are used in UPS systems, the requirement of low output voltage harmonics has made the sinusoidal PWM the most popular. In the case of a single-phase VSI with a four-transistor, three-level H-bridge, a double-edge PWM is usually a sufficient solution. The design of the discrete-time controller requires that a discrete-time or continuous model of a VSI be used. The first approach is based on the continuous to discrete-time transformation (the discretized plant) using the zero order hold (ZOH) model or the Tustin model [4,30]. The second approach is based on solving the state-space equations (the discrete plant) and linearizing the control function [1,21,30]. The third approach is based on the continuous model of a VSI, adjusting the quasi-continuous-time counterpart of the discrete-time controller, and finally, using a discrete-time controller (the discretized controller) [8,9]. The first and third solutions are based on the averaged models. The second solution considers the PWM features that distinguish the single and double-edge PWM. However, for a relatively high switching frequency, all of them gave similar results of control. The state variables x = [vOUT iLF iOUT]T, the control (input) variable u = vFIN, the output variable y = vOUT (Figure 1), and the state A and control B matrixes were calculated in Equation (11) when the ESR of the capacitor was neglected. The state-space equations x ˙ = A x + B u were solved for the time period of the switching on TONk (VFIN = VDC) and switching off TOFFk (VFIN = 0). However, we obtained nonlinear equations for the control, so the discrete-time control matrix GD was approximated exp(ATONk/2) ≈ I + ATONk/2 to linearize them [1,16,25]. Finally, the discrete-time state-space equations were obtained, as shown in Equations (11)–(15) [4,21,30].
A = [ 0 1 C F 1 C F 1 L F R L F L F 0 0 0 0 ] , B = [ 0 1 L F 0 ]
x k + 1 = A D x k + G D T O N k
A D = e A T s = Φ ( T s ) = L 1 [ ( s I A ) 1 ] | t = T s
G D = e A T s / 2 B V D C = Φ ( T s / 2 ) B V D C
A D = Φ ( T s ) = [ ϕ 11 ϕ 12 ϕ 13 ϕ 21 ϕ 22 ϕ 23 ϕ 31 ϕ 32 ϕ 33 ] , G D = [ g 11 g 21 g 31 ]
Where
ξ F = 1 2 R L F C F / L F ϕ 11 = [ cos ( ω F 0 T s ) + ξ F sin ( ω F 0 T s ) ] exp ( ξ F ω F 0 T s ) , ϕ 12 = 1 ω F 0 C F sin ( ω F 0 T s ) exp ( ξ F ω F 0 T s ) ,   φ 13 = φ 12 + R L F ( φ 11 1 ) , ϕ 21 = C F L F φ 12 ,   ϕ 22 = [ cos ( ω F 0 T s ) ξ F sin ( ω F 0 T s ) ] exp ( ξ F ω F 0 T s ) ,   ϕ 23 = 1 φ 11 , ϕ 31 = 0 ,   ϕ 32 = 0 ,   ϕ 33 = 1 g 11 = V D C ω F 0 sin ( ω F 0 T s 2 ) exp ( ξ F ω F 0 T s 2 ) , g 21 = V D C L F [ cos ( ω F 0 T s 2 ) ξ F sin ( ω F 0 T s 2 ) ] exp ( ξ F ω F 0 T s 2 ) ,   g 31 = 0
The damping ζF calculated from the space Equation (14) was not dependent on the load ZLOAD (Figure 1) because if IOUT was treated as an independent current source (Figure 2), the load was absent in the space equations. This is the difference in the Equations (16) and (3) describing ζF value.
The output voltage discrete-time transfer function for the open loop is a control function of the duty cycle, which is controlled by vCTRL and a function (−ZOUT, where ZOUT is an output impedance) of the output current Equation (17). The discrete-time control transfer function of a VSI Equation (24) considers an additional delay of one Ts of the discrete-time modulator H(z).
V O U T ( z ) = z 1 T s V D C G i n v V C T R L ( z ) Z O U T ( z ) I O ( z )
For the open-loop system, Equations (18) and (19) are used:
K V S I = V O U T ( z ) V C T R L ( z ) = z 1 T s V D C G i n v = T s V D C g 11 z 2 + T s V D C ( φ 12 g 21 φ 22 g 11 ) z 3 1 ( φ 11 + φ 22 ) z 1 + ( φ 11 φ 22 φ 12 φ 21 ) z 2
Z O U T ( z ) = φ 13 z 1 + ( φ 12 φ 23 φ 13 φ 22 ) z 2 1 ( φ 11 + φ 22 ) z 1 + ( φ 11 φ 22 φ 12 φ 21 ) z 2
Similar results of the discrete-time model for a VSI can be achieved using the c2d MATLAB function with the simplest ZOH discretization [30] or the discretization in the analytic way [8,9].

4. SISO Control of a VSI

The first presented example of an SISO control is a modern discretized PID controller (presented in detail in [8,9]), which uses the simplest SISO model from Figure 2 but includes an RLOAD. The quasi-continuous transfer function [8,9] of the plant Kp(s) (including RLOAD), considering a further ZOH discretization and one Ts delay of the modulator, is shown in Equation (20),
K p ( s ) = exp ( s T s 2 ) exp ( s T s ) F L C ( s ) ( 1 s T s 2 ) 1 s T s 2 1 + s T s 2 ω F 0 2 s 2 + 2 ξ F ω F 0 s + ( 1 + R L F R L O A D ) ω F 0 2
The discrete-time PID controller is as in Equation (21),
C ( z 1 ) = b 0 + b 1 z 1 + b 2 z 2 1 z 1
The difference control law for kV = 1 is as in Equation (22),
v C T R L ( k ) = v C T R L ( k 1 ) + b 0 ( v O U T ( k ) v R E F ( k ) ) + b 1 ( v O U T ( k 1 ) v R E F ( k 1 ) ) + b 2 ( v O U T ( k 2 ) v R E F ( k 2 ) )
For the experimental model: Ts = 1/fs = 1/2,5600 s, LF = 1 mH, CF = 50 μF, RLF = 1 Ω, RLOAD = 50 Ω; the PID controller transfer function from [8] will be b0 = 18.014, b1 = 33.495, b2 = 16.094. These coefficients were analyzed based on the root locus [8]. In Reference [9], these coefficients were slightly adjusted: b0 = 18.381, b1 = 34.179, b2 = 16.423.
A more sophisticated SISO system is presented in Figure 3 [10,11,31,32].
The concept of a Manabe CDM [10,11,12,31] design of the controller (T, S, R) is such that the coefficients of the closed-loop characteristic polynomial have assumed values, in the simplest case, of Standard Manabe Form. For kv = 1 (vREF is properly scaled), the output voltage of a closed-loop system is as in Equation (23),
v O U T ( z 1 ) = T N R D + S N v R E F ( z 1 ) Z O U T R D R D + S N I O U T ( z 1 )
The characteristic equation of a closed-loop system is as in Equation (24),
P ( z 1 ) = R ( z 1 ) D ( z 1 ) + S ( z 1 ) N ( z 1 ) = i = 0 n p z i z i
From Equation (18), Equation (25) can be derived:
K V S I = N ( z 1 ) D ( z 1 ) = a 2 z 2 + a 3 z 3 1 + b 1 z 1 + b 2 z 2
For a system with a disturbance (in our case IOUT), the degrees of R and S are equal to or higher than n−1, where n is the degree of D. In our case, we will assume Equation (26):
S ( z 1 ) = i = 0 2 s i z i ,   R ( z 1 ) = i = 0 2 r i z i ,   r 0 = 1
The Diophantine Equation (27) should be solved for r0 = p0 = 1.
( 1 + r 1 z 1 + r 2 z 2 ) ( 1 + b 1 z 1 + b 2 z 2 ) + ( s 0 + s 1 z 1 + s 2 z 2 ) ( a 2 z 2 + a 3 z 3 ) = i = 0 5 p z i z 1
To solve Equation (27), Equation (28) should be solved to obtain the ri and si coefficients,
[ 1 0 0 0 0 b 1 1 a 2 0 0 b 2 b 1 a 3 a 2 0 0 b 2 0 a 3 a 2 0 0 0 0 a 3 ] [ r 1 r 2 s 0 s 1 s 2 ] = [ p z 1 b 1 p z 2 b 2 p z 3 p z 4 p z 5 ]
where Equation (29),
a 2 = T s V D C g 11 ,   a 3 = T s V D C ( φ 12 g 21 φ 22 g 11 ) ,   b 1 = ( φ 11 + φ 22 ) ,   b 2 = φ 11 φ 22 φ 12 φ 21
The coefficients pi of Standard Manabe Form in a continuous system for the 5th degree of P(s) are the following:
p0(s0) = 1, p1(s1) = p0τ, p2(s2) = 0.4 p0τ2, p3(s3) = 0.08 p0τ3, p4(s4) = 0.008 p0τ4, p5(s5) = 0.0004 p0τ5,
where τ is the time constant of a closed-loop system. Satisfactory results of the control of the experimental model were achieved for τ = 4 Ts and for τ = 5 Ts.
Let us define the following discrete-time transfer function using the c2d function with the discretization cycle Ts = 1/2.5600 s and the default ZOH method Equation (30) [4]:
K ( z ) = c 2 d ( 1 i = 0 5 p i ( s ) s i , T s ) = i = 0 5 w i ( z ) z i i = 0 5 p z i ( z 1 ) z i
For τ = 4 Ts:
pz0(z0) = 1, pz1(z−1) = −1.327, pz2(z−2) = 0.6811, pz3(z−3) = −0.1826, pz4(z−4) = 0.0381, pz5(z-5) = −0.006738.
This last calculation enables vOUT = vREF to be kept in the steady state Equation (31).
t 0 = P ( z = 1 ) N ( z = 1 ) = V D C T s 1 + p z 1 + p z 2 + p z 3 + p z 4 + p z 5 φ 12 g 21 + ( 1 φ 22 ) g 11
For the experimental model, the solutions of Equations (28) and (31) were:
r0 = 1, r1 = 0.5898, r2 = 0.4218, s0 = 29.5050, s1 = −24.2037, s2 = −0.4607, t0/VDC = 0.1713.
However, in many cases t0 is adjusted individually. The difference control law for kv = 1 is as in Equation (32),
v C T R L ( k ) = r 1 v C T R L ( k 1 ) r 2 v C T R L ( k 2 ) + t 0 v R E F s 0 v O U T ( k ) s 1 v O U T ( k 1 ) s 2 v O U T ( k 2 )

5. MISO Control of VSI

MISO controllers with multiple controller inputs, i.e., the output voltage, filter inductor current, and output current, and a single output, i.e., voltage (Figure 1), effectively reduce the distortions of the output voltage for the different types of loads that are defined in [19]. One of the MISO controls that is described often lately is the passivity-based control (PBC) presented by Ortega [33] in 1989. PBC seems to be one of the best solutions for power conversion systems such as VSIs [20,34,35,36]. A VSI is presented as an energy transformation multiport device [37]. If the stored energy is less than the supplied energy, the system is passive. The “injection” of the appropriate damping [34] is basic for the control of the inductor and output currents; however, in the so-called improved PBC [34], there is direct feedback from the output voltage. The state variables that define the energy of the VSI are defined as in Equations (33) and (34). The output current, which does not exist in the function of the system energy, is treated as an independent disturbance [20].
x = [ L F i L F C F v O U T ] T = P [ i L F v O U T ] T
P = [ L F 0 0 C F ] ,   P 1 = [ 1 / L F 0 0 1 / C F ]
The total energy that is stored in the system is described by the Hamiltonian function H(x) in Equations (35) and (36),
H ( x ) = 1 2 ( L F i L F 2 + C F v O U T 2 ) = 1 2 x T P 1 x
[ i L F v O U T ] T = P 1 x = H ( x ) / ( x )
The error vector e is defined as Equations (37)–(39),
e = [ L F ( i L F i L F r e f ) C F ( v O U T v O U T r e f ) ] = P [ i L F i L F r e f v O U T v O U T r e f ]
H ( e ) = 1 2 ( L F ( i L F i L F r e f ) 2 + C F ( v O U T v O U T r e f ) 2 ) = 1 2 e T P 1 e
[ i L F i L F r e f v O U T v O U T r e f ] = P 1 e = H ( e ) / e
vOUTref is the reference, sinusoidal output voltage waveform;
iLFref is the calculated reference current of the inductor.
The equilibrium of a closed-loop system is asymptotically stable [38] and is achieved if H(e) has a minimum in x = xref in Equation (40),
H ( e ) x | x = x r e f = 0 , 2 H ( e ) x 2 | x = x r e f > 0
The system is passive if the time derivative H(e) is negative in Equation (41),
d H ( e ) d t < 0
Two equations are used—the first for a closed-loop system Equation (42) [20] and the second for an open-loop system, as in Equation (43),
e ˙ = [ J ( R + R a ) ] P 1 e
x ˙ = [ J R ] P 1 x + [ V D C 0 ] m + [ 0 1 ] i O U T
where the interconnection matrix, J, and the damping matrix, R, can be defined as in Equation (44),
J = [ 0 1 1 0 ] ,   R = [ R L F 0 0 0 ]
Ra (the PBC controller) is the matrix Equation (45) of the injected damping, Ri is the gain of the current error, and Kv is the conductive gain of the voltage error.
R a = [ R i 0 0 K v ]
Subtracting Equation (42) from Equation (43) results in the control law Equation (46),
e ˙ x ˙ = [ J R ] P 1 ( e x ) R a P 1 e [ V D C 0 ] m [ 0 1 ] i O U T
For vCTRL = mVDC, the difference control law of a single-phase VSI with a PBC is shown in Equations (47) and (48),
v C T R L ( k ) = R i i L F ( k ) + ( R i + R L F ) i L F r e f ( k ) + L F i L F r e f ( k ) i L F r e f ( k 1 ) T c + v O U T r e f ( k )
i L F r e f ( k ) = K v [ v O U T r e f ( k ) v O U T ( k ) ] + C F v O U T r e f ( k ) v O U T r e f ( k 1 ) T c + i O U T ( k )
The inductor current is an integral of a function of the output voltage, which is important for the steady-state error in the control law. The adjustment of the PBC gains, Ri and Kv, is a problem. The matrix R + Ra should be positively defined to fill this requirement Equation (41) [20]. In practice, this means that Ri + RLF and Kv should be positive. Calculating the upper restrictions of their values can be difficult [34]. The higher their values are, the higher the convergence of the error tracking is. The characteristic polynomial of a closed-loop system with a PBC is shown in Equation (49),
det { [ J ( R + R a ) ] P 1 λ 1 } = [ ( R L F R i ) 1 / L F λ 1 / C F 1 / L F ( K v ) 1 / C F λ ]
The roots λ1,2 Equation (50) of the characteristic polynomial Equation (49) of a closed-loop system Equation (42) will always have the real part negative (Figure 4a,b) for the positively defined R + Ra matrix.
λ 1 , 2 = { [ ( R L F + R i ) C F + L F K V ] ± [ ( R L F + R i ) C F + L F K v ] 2 4 L F C F [ 1 + ( R L F + R i ) K v ] } 2 L F C F
However, in the simulations and in the experimental model, overly high values of Ri and Kv led to oscillations in the control voltage and output voltage (Figure 5). An overly high value of modulation index M [39,40] for the nonlinear rectifier RC load led to the saturation of the PWM modulator (Figure 6). Both of these cases produced the same result. The reason can be the disability of fast current changing in the filter inductor [40]. Therefore, the best adjustment of the PBC controller gains Ri and Kv is to determine the lowest THD (for the nonlinear rectifier RC load).
Figure 6 presents two problems. The first is saturation in the PWM modulator due to overly high values of the PBC gains (Ri and KV). The second case can be presented comparing the root loci from Figure 4a,b and Figure 6. For Ri = 10, oscillations can occur for KV = 0–1, for Ri = 20, and for KV = 0.5–1.5. The second problem was not so important. The value of M = 0.3 was too low in practice. It was experimentally checked in the VSI model and finally set at M = 0.8, Ri = 5, and KV = 0.5. The equations from the literature did not solve the problem of restricting the upper Ri and KV values. From Figure 5, it is obvious that an MISO control with a low M value is able to efficiently control the output current. The presented PBC control was further named IPBC2 because compared to the basic PBC theory [34], it used the output voltage as in [34] and its derivative in the final control law. This derivative was absent in the improved PBC (IPBC) that was presented in [34].
The idea of the calculation of the upper limit of Ri gain was presented in [34]. The derivative of the control voltage vCTRL in one switching cycle should be lower than the maximum carrier slope in the PWM modulator. For the double-edge modulation and recalculating it to the DC voltage level, it is equal to VDC/(Ts/2) [34]. However, it is shown in Figure 5 and Figure 6 that the mutual dependency of the two PBC controller gains is very important. In one sampling period, we can assume d(vOUTref)/dt ≈ 0. Therefore, from Equation (48) we can calculate:
i L F r e f ( k T s ) = ( 1 R L O A D K v ) v O U T ( k T s ) + v O U T r e f ( k T s )
d i L F R E F ( k T s ) d t ( 1 R L O A D K v ) d v O U T ( k T s ) d t
From Equations (47), (51) and (52), we can calculate Equation (53):
d v C T R L ( k T s ) d t L F ( 1 R L O A D K v ) d 2 v O U T ( k T s ) d t 2 + ( R i + R L F e ) ( 1 R L O A D K v ) d v O U T ( k T s ) d t R i d i L F ( k T s ) d t
In one switching cycle, for RLOAD = ∞, we can assume Equation (54):
d i L F ( k T s ) d t | max , min ± V D C L F ,   d v O U T ( k T s ) d t | max i L F C F ,   d 2 v O U T ( k T s ) d t 2 | max d d t ( i L F C F ) | max ± V D C L F C F
| d v C T R L ( k T s ) d t | max K v [ L F + ( R i + R L F e ) T s ] V D C L F C F + R i V D C L F
The final restriction of PBC gains, Ri and Kv, from Equation (55) is Equation (56),
K v [ L F + ( R i + R L F e ) T s ] 1 C F + R i < 2 L F T s
The restrictions of the operating area, shown in Figure 7, can be one of the reasons for the THD increase in the higher values of PBC gains presented in Figure 5 and Figure 6. However, Equation (56) does not consider the modulation index M value because the maximum carrier slope does not depend on M (the reference voltage depends on M). Equation (56) is calculated for RLOAD = ∞, while the restrictions of the Kv value for the existing load resistance will be slightly lower.

6. Comparison of Experimental Voltage Source Inverters with PID, CDM, and PBC Control Systems

The comparison of the SISO-PID (the version from [8,9]), SISO-CDM, and MISO-PBC-IPBC2 systems was based on measurements of the VSI output voltage for standard loads. Figure 8a–d presents the output voltage for the standard [19] nonlinear rectifier RC load for R = 100 Ω and C = 430 μF. The vertical axis of Figure 8a–d is scaled in units of the 13th bit bipolar Analog To Digital Converter ADC. The actual voltage amplitude was about 60 V. The output filter (LF = 1 mH, CF = 50 μF) resonant frequency was about 712 Hz. The initial analysis (Section 2) showed that the feedback loop is only important for the neighborhood of this frequency (about 14th harmonic) and below. Therefore, a sufficient frequency range for comparing the controllers was from the 2nd up to the 30th harmonic. Figure 8a–d presents the output voltage waveforms without feedback and with the PID, CDM, and PBC-IPBC2 controllers. The THD coefficient was calculated for each case. The PID and CDM controllers had similar THD, while for the PBC, it was a bit lower. However, the output voltage waveforms were different. The harmonics spectra are presented in Figure 9a–d. It can be noticed that the harmonics spectrum is different for the different controllers. Some of the harmonics are increased (e.g., the 2nd harmonic for the CDM controller increased because the voltage waveform was not ideally symmetrical). The importance of a decrease or increase in a specific harmonic was dependent on its relative amplitude. Therefore, the influence of the feedback loop on the particular harmonic can be defined as in Equation (57),
x n = h N F B n 20 log 10 h N F B n h C T R L n
where hNFBn is the n-th relative harmonic amplitude of the output voltage without feedback and hCTRLn is the n-th relative harmonic using the selected controller. The sum of the xn values is a measure of the feedback quality. The better the quality, the higher the value of the coefficient. Therefore, the final definition of the control quality factor (CQF) is as in Equation (58),
C Q F = 1 n h N F B n 20 log 10 h N F B n h C T R L n
where m is the harmonic order above the filter resonant frequency, in our case m = 30.
The higher the CQF, the better the control. It can be seen from Figure 8a–d that for the nonlinear rectifier RC load, the best controller was the PBC (CQF = 1.51), the second best was the CDM (CQF = 1.41), and the third best was the PID (CQF = 1.2). The CQF measure enabled us to distinguish the quality of PID and CDM controls that had almost the same THD. The nonlinear load caused harmonic disturbances, and the CQF was the proper measure of their damping.
The quality of the control should also be determined for a dynamic load change. Figure 10a–d shows the over-and undershoots after a dynamic load decrease and increase (45 Ω to 500 Ω and 500 Ω to 45 Ω) and the settling time. The lowest over- and undershoots (about 3% and −2%) were for the CDM control, the shortest settling time (about 1 ms in both cases) was for the PID control, and the worst parameters were for the PBC-IPBC2 (about 6% and −5.1%, 3.5 and 2 ms). It would be difficult to improve the dynamic properties of the PBC because increasing the gains Ri and KV would lead to higher oscillations and a longer settling time, while decreasing them would lead to higher over- and undershoots. In the PBC, we directly control the filter inductor current as in Equations (47) and (48) but the inductance disables fast current changes, which can lead to the saturation of the modulator. The capacitance in the output filter has a significant influence on the over- and undershoots.

7. Materials and Methods

The presented simulations were performed using MATLAB-Simulink for the VSI model with LF = 1 mH, RLF = 1 Ω, CF = 50 µF and the switching frequency fs = 25,600 Hz (Figure 5, Figure 6 and Figure 7). The modulation index in simulations was M = 0.3. The first PWM scheme [41,42] of 3-level modulation was used. The VSI experimental model had the same nominal parameter values as the MATLAB model except that the modulation index was M = 0.8. The filter coil core was made of Super-MSS™ material [23,42], and the output capacitor was of the Metallized Polypropylene Capacitor MKP type. The transistor switches of the H-bridge were International Rectifier HEXFETs - IRFP 360, which have RDS(on) = 0.20 Ω (@10 A). The PID [8,9], CDM, and IPBC2 control laws were programmed in an STM32F405VG microcontroller by means of Keil µVision C++. The static load resistance varied from 50 to 2000 Ω, and the dynamic load resistance was a step load from 50||500 to 500 Ω and from 500 Ω to 50||500. The nonlinear rectifier RC load was for R = 100 Ω and C = 430 µF. The amplitude of the experimental VSI output voltage was 60 V, and the values of the voltages in Figure 8 and Figure 10 are presented in units of the ADC converter (the −4093 to 4093 range). The waveforms (two fundamental periods) were transmitted from the inverter model to the PC with the visualization software (Figure 11), then computed and plotted using MATLAB (Figure 8 and Figure 10).

8. Results

The results of the measurements of the VSI experimental model are presented in Table 1. With the inclusion of the new CQF factor, an increase in the resolution of the estimation of the three different controls’ results should be enabled. We can use simple and cheap SISO-PID or SISO-CDM controllers (without the filter coil current and the VSI output current sensors) in the case of the static and dynamic linear loads of the VSI. However, in the case of the nonlinear rectifier RC load, a better solution is MISO-IPBC2. The restrictions of the Kv and Ri gains in the MISO-IPBC2 (Figure 5, Figure 6 and Figure 7) complicate its further improvement. The requirement of the restriction of the product of the modulation index M and the filter inductor value LF (the higher the LF, the lower the M) is a problem that exists for all the controller types and was described, e.g., in [39,40].

9. Conclusions

A comparison of SISO-PID, SISO-CDM, and MISO-PBC-IPBC2 systems showed that for the nonlinear rectifier RC load defined in the IEC 62040-3 standard [19], the best results (the lowest distortions of the output voltage) were for the PBC-IPBC2 controller that had the gains adjusted after the initial simulations (Figure 5, Figure 6 and Figure 7). The measurement of the currents enabled them to be better shaped (Figure 5). A new control quality factor was defined and enabled a better resolution in the description of the control results. An analysis of the root locus (Figure 4a,b) showed that the voltage source inverter with the PBC-IPBC2 control should always be stable for any positive gains Ri and KV; however, the simulations (Figure 5, Figure 6 and Figure 7) and measurements of the experimental model showed that oscillations in the output voltage existed for the higher values of gains. It is obvious that the MISO-PBC-IPBC2 system with the output current and filter inductor current measurements better damped the important harmonics (Figure 9a–d). The output voltage parameters (the over- and undershoots, the settling time) for the linear dynamic load were worse for the PBC-IPBC2, and we were not able to improve them by changing the values of the controller gains. An increase in the gains decreased the over- and undershoots but increased the settling time because of oscillations. The reason for the worse dynamic response of the PBC could be the direct control the filter inductor current, during which the inductance disables fast current changes, which can lead to the saturation of the modulator [40]. The static error, which was calculated as the decrease in output voltage amplitude after an increase in the load (500 to 45 Ω), was low and was similar for all of the controllers (about 2% vs. 4% without feedback), which is understandable because the static error was mainly reduced by the integration action of the output voltage. All of the controllers had such an action (the PBC-IPBC2 controller due to the control of the inductor current, which was dependent on the integration of the function of the VSI output voltage). All of the controllers cause small oscillations (Figure 10a–d) in the VSI output voltage (in the experimental model) close to the zero crossing, which could have been the result of the wind up effect of the controller (the voltage error was integrated, and the three-level PWM for the voltage crossing zero was not able to control the output voltage efficiently with very short pulses). To summarize, in the case of a nonlinear load, the more expensive MISO-PBC-IPBC2 controller had much lower distortions of the VSI output voltage, while in the case of dynamic loads, the less expensive SISO-PID and SISO-CDM controllers were sufficient and the static error was the same for all the controllers. Therefore, the final selection of either a simple and inexpensive SISO with one measurement trace or an expensive MISO with three measurement traces depends on the proposed function of the VSI.

Author Contributions

Conceptualization, Z.R. and K.B.; methodology, Z.R. and K.B.; software, Z.R.; validation, Z.R., and K.B.; formal analysis, Z.R. and K.B.; investigation, Z.R. and K.B.; resources, Z.R. and K.B.; writing—original draft preparation, Z.R.; writing—review and editing, Z.R. and K.B.; visualization, Z.R.; supervision, Z.R.; project administration, Z.R. and K.B.; funding acquisition, Z.R. and K.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by the Polish Ministry of Science and Higher Education funding for statutory activities (BK-284/RAU11/2020). Authors were supported by the Polish National Centre for Research and Development, grant no. TANGO3/427467/NCBR/2019.

Acknowledgments

The calculations were performed using IT infrastructure that was funded by the GeCONiI project (POIG.02.03.01-24-099/13).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kawamura, A.; Chuarayapratip, R.; Haneyoshi, T. Deadbeat control of PWM inverter with modified pulse patterns for uninterruptible power supply. IEEE Trans. Ind. Electron. 1988, 35, 295–300. [Google Scholar] [CrossRef]
  2. Kusko, A.; Galler, D.; Medora, N. Output impedance of PWM UPS inverter-feedback vs. filters. In Conference Record of the 1990 IEEE Industry Applications Society Annual Meeting, 7–12 October; Institute of Electrical and Electronics Engineers: New York, NY, USA, 1990; Volume 2. [Google Scholar] [CrossRef]
  3. Deng, H.; Srinivasan, D.; Oruganti, R. Modeling and control of single-phase UPS inverters: A survey. In Proceedings of the International Conference on Power Electronics and Drives Systems, PEDS 2005, Kuala Lumpur, Malaysia, 28 November–1 December 2005; Volume 2, pp. 848–853. [Google Scholar] [CrossRef]
  4. Rymarski, Z. Design Method of Single-Phase Inverters for UPS Systems. Int. J. Electron. 2009, 96, 521–535. [Google Scholar] [CrossRef]
  5. Rymarski, Z. The influence of the output filter parameters tolerance and the load on the PID/CDM control of the single phase VSI. Przegląd Elektrotechniczny 2011, 87, 114–117. [Google Scholar]
  6. Rymarski, Z.; Bernacki, K. Drawbacks of impedance networks. Int. J. Circuit Theory Appl. 2018, 46, 612–628. [Google Scholar] [CrossRef]
  7. Rymarski, Z.; Bernacki, K.; Dyga, Ł. Decreasing the single phase inverter output voltage distortions caused by impedance networks. IEEE Trans. Ind. Appl. 2019, 55, 7586–7594. [Google Scholar] [CrossRef]
  8. Blachuta, M.; Rymarski, Z.; Bieda, R.; Bernacki, K.; Grygiel, R. Design, Modeling and Simulation of PID Control for DC/AC Inverters. In Proceedings of the 24th International Conference on Methods and Models in Automation and Robotics, Międzyzdroje, Poland, 26–29 August 2019; pp. 428–433. [Google Scholar] [CrossRef]
  9. Blachuta, M.; Rymarski, Z.; Bieda, R.; Bernacki, K.; Grygiel, R. In Proceedings of the Continuous-time approach to discrete-time PID Control for UPS inverters—A case study. In Proceedings of the Asian Conference on Intelligent Information and Database Systems, ACIIDS, Phuket, Thailand, 23–26 March 2020; pp. 1–12. [Google Scholar]
  10. Manabe, S. Coefficient diagram method. IFAC Proc. Vol. 1998, 31, 211–222. [Google Scholar] [CrossRef]
  11. Manabe, S. Importance of coefficient diagram in polynomial method. In Proceedings of the 42nd IEEE Conference on Decision and Control, Maui, HI, USA, 9–12 December 2003; pp. 3489–3494. [Google Scholar] [CrossRef]
  12. Coelho, J.P.; Pinho, T.M.; Boaventura-Cunha, J. Controller system design using the coefficient diagram method. Arab. J. Sci. Eng. 2016, 41, 3663–3681. [Google Scholar] [CrossRef]
  13. Tomizuka, M.; Tsao, T.C.; Chew, K.K. Analysis and synthesis of discrete-time repetitive controllers. Asme J. Dyn. Syst. Meas. Control 1989, 111, 353–358. [Google Scholar] [CrossRef]
  14. Luo, F.L.; Ye, H.; Rashid, M. Digital Power Electronics and Applications; Elsevier Academic Press: Cambridge, MA, USA, 2006. [Google Scholar] [CrossRef]
  15. Gangling, Z.; Guan, M.; Wang, Y. Application of repetitive control for aeronautical static inverter. In Proceedings of the 2007 2nd IEEE Conference on Industrial Electronics and Applications, ICIEA, Harbin, China, 23–25 May 2007; pp. 121–125. [Google Scholar] [CrossRef]
  16. Rech, C.; Pinheiro, H.; Gründling, H.A.; Hey, H.L.; Pinheiro, J.R. Comparison of digital control techniques with repetitive integral action for low cost PWM inverters. IEEE Trans. Power Electron. 2003, 18, 401–410. [Google Scholar] [CrossRef]
  17. Lidozzi, A.; Ji, C.; Solero, L.; Crescimbini, F.; Zanchetta, P. Load-adaptive zero-phase-shift direct repetitive control for stand-alone four-leg VSI. IEEE Trans. Ind. Appl. 2016, 52, 4899–4908. [Google Scholar] [CrossRef]
  18. Liu, Z.; Zhang, B.; Zhou, K.; Yang, Y.; Wang, J. Virtual variable sampling repetitive control of single-phase DC/AC PWM converters. IEEE J. Emerg. Sel. Top. Power Electron. 2019, 7, 1837–1845. [Google Scholar] [CrossRef] [Green Version]
  19. IEC 62040-3:2011. Uninterruptible Power Systems (UPS)-Part 3: Method of Specifying the Performance and Test Requirements. Available online: https://webstore.iec.ch/publication/6344 (accessed on 21 July 2020).
  20. Serra, F.M.; De Angelo, C.H.; Forchetti, D.G. IDA-PBC control of a DC–AC converter for sinusoidal three-phase voltage generation. Int. J. Electron. 2017, 104, 93–110. [Google Scholar] [CrossRef]
  21. Xie, R.; Hao, X.; Yang, X.; Chen, W.; Huang, L.; Wang, C. An exact discrete-time model considering dead-time nonlinearity for an H-bridge grid-connected inverter. In Proceedings of the 2014 International Power Electronics Conference (IPEC-Hiroshima 2014—ECCE ASIA), Hiroshima, Japan, 8–12 May 2014; pp. 2950–2953. [Google Scholar] [CrossRef]
  22. Yuan, W.; Wang, Y.; Liu, D.; Deng, F.; Chen, Z. Robust droop control of ac microgrid against nonlinear characteristic of inductor. In Proceedings of the IEEE 10th International Symposium on Power Electronics for Distributed Generation Systems (PEDG), Xi’an, China, 3–6 June 2019; pp. 642–647. [Google Scholar] [CrossRef]
  23. 2019 Alloy Powder Products Catalog, Micrometals Powder Core Solutions. Available online: https://www.micrometals.com/design-and-applications/literature/ (accessed on 21 July 2020).
  24. Bernacki, K.; Rymarski, Z.; Dyga, Ł. Selecting the coil core powder material for the output filter of a voltage source inverter. Electron. Lett. 2017, 53, 1068–1069. [Google Scholar] [CrossRef]
  25. Rymarski, Z. The discrete model of the power stage of the voltage source inverter for UPS. Int. J. Electron. 2011, 98, 1291–1304. [Google Scholar] [CrossRef]
  26. Cheng, X.; Chen, Y.; Chen, X.; Zhang, B.; Qiu, D. An extended analytical approach for obtaining the steady-state periodic solutions of SPWM single-phase inverters. In Proceedings of the IEEE Energy Conversion Congress and Exposition (ECCE), Cincinnati, OH, USA, 24 August 2017; pp. 1311–1316. [Google Scholar] [CrossRef]
  27. Cheng, Y.; Zha, X.; Liu, Y. Nonlinear modeling of inverter using the Hammerstein’s approach. In Proceedings of the INTELEC 2009—31st International Telecommunications Energy Conference, Incheon, Korea, 18–22 October 2009; pp. 1–4. [Google Scholar] [CrossRef]
  28. Dyga, Ł.; Rymarski, Z.; Bernacki, K. The wavelet-aided methods for evaluating the output signal that is designated for uninterruptible power supply systems. Przegląd Elektrotechniczny 2020. accepted paper. [Google Scholar]
  29. van der Broeck, H.W. Analytical calculation of the harmonic effects of single phase multilevel PWM inverters. In Proceedings of the 29th Annual Conference of the IEEE Industrial Electronics Society IECON, 2–6 November 2003; pp. 243–248. [Google Scholar] [CrossRef]
  30. Rymarski, Z.; Bernacki, K. Different approaches to modelling single-phase voltage source inverters for uninterruptible power supply systems. IET Power Electron. 2016, 9, 1513–1520. [Google Scholar] [CrossRef]
  31. Hamamci, S.E.; Kaya, I.; Koksal, M. Improving performance for a class of processes using coefficient diagram method. In Proceedings of the 9th Mediterranean Conference on Control and Automation, MED’01, Dubrovnik, Croatia, 27–29 June 2001; pp. 1–6. [Google Scholar]
  32. Astrom, K.J.; Wittenmark, B. Computer-Controlled Systems: Theory and Design, 3rd ed.; Dover Publications Inc.: Upper Saddle River, NJ, USA; Mineola, NY, USA, 2011. [Google Scholar]
  33. Ortega, R.M.; Spong, M.W. Adaptive motion control of rigid robots: A tutorial. Automatica 1989, 25, 877–888. [Google Scholar] [CrossRef]
  34. Komurcugil, H. Improved passivity-based control method and its robustness analysis for single-phase uninterruptible power supply inverters. IET Power Electron. 2015, 8, 1558–1570. [Google Scholar] [CrossRef]
  35. Gui, Y.; Wei, B.; Li, M.; Guerrero, J.M.; Vasquez, J.C. Passivity-based coordinated control for islanded AC microgrid. Appl. Energy 2018, 229, 551–561. [Google Scholar] [CrossRef] [Green Version]
  36. Meshram, R.V.; Bhagwat, M.; Khade, S.; Wagh, S.R.; Aleksandar, M.; Stankovic, A.M.; Singh, N.M. Port-controlled phasor hamiltonian modeling and IDA-PBC control of solid-state transformer. IEEE Trans. Control Syst. Technol. 2019, 27, 161–174. [Google Scholar] [CrossRef]
  37. Hill, D.; Zhao, J.; Gregg, R.; Ortega, R. 20 Years of Passivity–Based Control (PBC): Theory and Applications; CDC Workshop: Shanghai, China, 15 December 2009. [Google Scholar]
  38. Wang, Z.; Goldsmith, P. Modified energy-balancing-based control for the tracking problem. IET Control Theory Appl. 2008, 2, 310–312. [Google Scholar] [CrossRef]
  39. Rymarski, Z. The analysis of output voltage distortion minimization in the 3-phase VSI for the nonlinear rectifier RoCo. Przegląd Elektrotechniczny 2009, 85, 127–132. [Google Scholar]
  40. Rymarski, Z.; Bernacki, K.; Dyga, Ł.; Davari, P. Passivity-based control design methodology for UPS systems. Energies 2019, 12, 4301. [Google Scholar] [CrossRef] [Green Version]
  41. Bernacki, K.; Rymarski, Z. Electromagnetic compatibility of voltage source inverters for uninterruptible power supply system depending on the pulse-width modulation scheme. IET Power Electron. 2015, 8, 1026–1034. [Google Scholar] [CrossRef]
  42. Rymarski, Z.; Bernacki, K.; Dyga, Ł. Measuring the power conversion losses in voltage source inverters. AEU-Int. J. Electron. Commun. 2020, p. 153359, in press. Available online: https://doi.org/10.1016/j.aeue.2020.153359 (accessed on 17 July 2020).
Figure 1. Single-phase voltage source inverter with single input single output (SISO) or multi-input single output (MISO) control.
Figure 1. Single-phase voltage source inverter with single input single output (SISO) or multi-input single output (MISO) control.
Energies 13 04100 g001
Figure 2. Simplified continuous SISO model of a voltage source inverter (VSI) with switching frequency fs.
Figure 2. Simplified continuous SISO model of a voltage source inverter (VSI) with switching frequency fs.
Energies 13 04100 g002
Figure 3. An advanced discretized SISO system.
Figure 3. An advanced discretized SISO system.
Energies 13 04100 g003
Figure 4. The root locus of the characteristic polynomial of a closed-loop system with passivity-based control (PBC) for LF = 1 mH, CF = 50 μF: (a) RLF + Ri = 10 Ω, (b) RLF + Ri = 20 Ω.
Figure 4. The root locus of the characteristic polynomial of a closed-loop system with passivity-based control (PBC) for LF = 1 mH, CF = 50 μF: (a) RLF + Ri = 10 Ω, (b) RLF + Ri = 20 Ω.
Energies 13 04100 g004
Figure 5. The simulation of the output voltage and current waveforms for rectifier R = 100 Ω and C = 430 μF for M = 0.3 and different values of Ri and KV.
Figure 5. The simulation of the output voltage and current waveforms for rectifier R = 100 Ω and C = 430 μF for M = 0.3 and different values of Ri and KV.
Energies 13 04100 g005
Figure 6. The adjustment of the Ri and Kv gains of the PBC (simulation) for M = 0.3.
Figure 6. The adjustment of the Ri and Kv gains of the PBC (simulation) for M = 0.3.
Energies 13 04100 g006
Figure 7. The restrictions (worst case) of the maximum Kv gain due to the Ri gain in the function of LF, for the constant CF to keep the derivative of the control voltage below the maximum carrier slope in the PWM modulator for LF = 1–5 mH, CF = 10 and 50 µF, RLF = 1 Ω and RLOAD = ∞, according to Equation (56).
Figure 7. The restrictions (worst case) of the maximum Kv gain due to the Ri gain in the function of LF, for the constant CF to keep the derivative of the control voltage below the maximum carrier slope in the PWM modulator for LF = 1–5 mH, CF = 10 and 50 µF, RLF = 1 Ω and RLOAD = ∞, according to Equation (56).
Energies 13 04100 g007
Figure 8. Measured VSI output voltage for the different types of control systems for the standard RC load, (a) without the feedback loop; (b) the SISO-PID controller; (c) the SISO-CDM controller; (d) the MISO-PBC-IPBC2 controller
Figure 8. Measured VSI output voltage for the different types of control systems for the standard RC load, (a) without the feedback loop; (b) the SISO-PID controller; (c) the SISO-CDM controller; (d) the MISO-PBC-IPBC2 controller
Energies 13 04100 g008
Figure 9. Measured harmonics spectrum of the different types of control systems, (a) without the feedback loop; (b) the SISO-PID controller; (c) the SISO-CDM controller; (d) the MISO-PBC-IPBC2 controller.
Figure 9. Measured harmonics spectrum of the different types of control systems, (a) without the feedback loop; (b) the SISO-PID controller; (c) the SISO-CDM controller; (d) the MISO-PBC-IPBC2 controller.
Energies 13 04100 g009
Figure 10. The over- and undershoots after a decrease and increase in the dynamic load (45 Ω to 500 Ω and 500 Ω to 45 Ω) and the settling time, (a) without the feedback loop; (b) the SISO-PID controller; (c) the SISO-CDM controller; (d) the MISO-PBC-IPBC2 controller.
Figure 10. The over- and undershoots after a decrease and increase in the dynamic load (45 Ω to 500 Ω and 500 Ω to 45 Ω) and the settling time, (a) without the feedback loop; (b) the SISO-PID controller; (c) the SISO-CDM controller; (d) the MISO-PBC-IPBC2 controller.
Energies 13 04100 g010
Figure 11. The experimental VSI model and PC with dedicated visualization software.
Figure 11. The experimental VSI model and PC with dedicated visualization software.
Energies 13 04100 g011
Table 1. Parameters of the VSI output voltage with standard loads for the different controllers.
Table 1. Parameters of the VSI output voltage with standard loads for the different controllers.
Type of ControlNo FeedbackSISO-PIDSISO-CDMMISO-PBC-IPBC2
Rectifier RC load (R = 100 Ω, C = 430 μF)THD = 4.714%THD = 1.78% CQF = 1.2THD = 1.88% CQF = 1.41THD = 1.33% CQF = 1.51
Load decrease (45 Ω to 500 Ω)ΔVOUT = 10% Settling: 4 msΔVOUT = 3.4% Settling: 1 msΔVOUT = 3.0% Settling: 1.5 msΔVOUT = 6% Settling: 3.5 ms
Load increase (500 Ω to 45 Ω)ΔVOUT = −9% Settling: 3.5 msΔVOUT = −3.2% Settling: 1 msΔVOUT = −2.1% Settling: 1.2 msΔVOUT = −5.1% Settling: 2 ms
Static error of the amplitude (500 Ω to 45 Ω)4%2%2%2%

Share and Cite

MDPI and ACS Style

Rymarski, Z.; Bernacki, K. Different Features of Control Systems for Single-Phase Voltage Source Inverters. Energies 2020, 13, 4100. https://doi.org/10.3390/en13164100

AMA Style

Rymarski Z, Bernacki K. Different Features of Control Systems for Single-Phase Voltage Source Inverters. Energies. 2020; 13(16):4100. https://doi.org/10.3390/en13164100

Chicago/Turabian Style

Rymarski, Zbigniew, and Krzysztof Bernacki. 2020. "Different Features of Control Systems for Single-Phase Voltage Source Inverters" Energies 13, no. 16: 4100. https://doi.org/10.3390/en13164100

APA Style

Rymarski, Z., & Bernacki, K. (2020). Different Features of Control Systems for Single-Phase Voltage Source Inverters. Energies, 13(16), 4100. https://doi.org/10.3390/en13164100

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