Chaos in Power Electronics: An Overview
Chaos in Power Electronics: An Overview
Chaos in Power Electronics: An Overview
Mario di Bernardo
and
Chi K. Tse
Facolt
a di Ingegneria
Universit
a del Sannio in Benevento
82100 Benevento, Italy
(on leave from Dept. Eng. Maths,
University of Bristol, U.K.)
dibernardo@unisannio.it
Abstract
Power electronics is rich in nonlinear dynamics. Its operation is
characterized by cyclic switching of circuit topologies, which gives
rise to a variety of nonlinear behavior. This chapter provides an
overview of the chaotic dynamics and bifurcation scenarios observed
in power electronics circuits, emphasizing the salient features of the
circuit operation and the modelling strategies. This chapter covers
the modelling approaches, analysis methods, and a classication of
the common types of bifurcations observed in power electronics.
1 Financial support for C.K. Tse to undertake this work has been provided by the Hong
Kong Research Grants Council under Grant PolyU-5131/99E.
317
318
16.1
Introduction
+
D A
319
+
C
Vo
Vin
(a)
C
Vo
+
(b)
L
Vin
H
H
H
H
C
+
Vo
(c)
FIGURE 16.1
Simple dc/dc converters. (a) Buck converter; (b) buck-boost converter; (c)
boost converter.
16.2
Regardless of its specic function, a power electronics circuit operates by toggling its topology among a set of linear or nonlinear circuit topologies, under
the control of a feedback system. As such, they can be regarded as piecewiseswitched circuits.
For example, in simple dc/dc converters, such as the ones shown in Fig. 16.1,
an inductor is switched between the input and the output through an appropriate switching element (labelled as S in the gure). The way in which the
inductor is switched determines the output voltage level and transient behavior. Usually, a semiconductor switch and a diode are used to implement the
said switching, and through the use of a feedback control circuit, the relative
durations of the various switching intervals can be continuously adjusted. Such
feedback action eectively controls the dynamics and steady-state behavior of
the circuit. Thus, both the circuit topology and the control method determines
the dynamical behavior of a power electronics circuit.
320
16.2.1
Many power electronics converters are constructed on the basis of the three
simple converters shown in Fig. 16.1. In a typical period-1 operation, the
switch S and the diode D are turned on and o in a cyclic and complementary
fashion, under the command of a pulse-width modulator. When the switch
S is closed (the diode D is open), the inductor current ramps up. When the
switch S is open (the diode D closed), the inductor current ramps down. The
duty ratio, dened as the fraction of a repetition period during which S is
closed, is continuously controlled by a feedback circuit that aims to maintain
the output voltage at a xed level even under input and load variations. Two
typical feedback arrangements are shown in the next subsection.
The behavior of a dc/dc converter is greatly inuenced by its operating mode.
Typically, we can distinguish two dierent modes of operation, namely, continuous conduction mode (CCM) and discontinuous conduction mode (DCM).
In CCM, the inductor current is maintained non-zero throughout the entire
repetition cycle. This happens when the inductance L is relatively large, or
the load current demand is relatively high. Moreover, in DCM, the inductor
current is zero for an interval of time within a cycle. This happens when the
inductance L is relatively small, or the load current demand is low, causing the
inductor current to fall to zero as the inductor is being discharged. During the
interval of zero inductor current, both switch S and diode D are open.
For both CCM and DCM, the output voltage has a xed relationship with
the input voltage, as determined by the duty ratio. A steady-state relationship
can be easily found from a volt-time balance consideration of the inductor.
Table 16.1 shows the steady-state output voltage expressions for stable period1 operation [5]. Here, we denote by s the steady-state duty ratio, and by T
the repetition period.
As we will see in this chapter, although the period-1 stable operation is the
preferred operation for most industrial applications, it represents only one particular operating regime. Because of the existence of many possible operating
regimes, it would be of practical importance to have a thorough understanding
of what determines the behavior of the circuit so as to guarantee a desired
operation or to avoid an undesirable one.
16.2.2
Most dc/dc converters are designed to deliver a regulated output voltage. The
control of dc/dc converters usually takes on two approaches, namely voltage
feedback control and current-programmed control, also known as voltage-mode
and current-mode control, respectively [6]. In voltage-mode control, the output
voltage is compared with a reference to generate a control signal which drives
S
Vin
321
6
+
D A
+
C
Vo
Zf
H
HH
+
COMP
+
PWM signal HH
Vramp
R1
Vref
R2
(a)
L
Vin
L
i
R
H
H
C
Vo
Q
S
clk
HH
+
COMP H
Zf
H
HH
+
R1
Vref
R2
(b)
FIGURE 16.2
Typical control approaches for dc/dc converters: (a) Voltage-mode control; (b)
current-mode control.
322
Converter type
CCM
DCM
Buck
Vo = s Vin
Vo = 2Vin 1 + 1 +
Buck-boost
Vo =
Boost
1
Vo =
Vin
1 s
s
Vin
1 s
RT
2L
8L
RT s2
1
Vo = Vin s
Vin
Vo =
2
1+
22 L
1+ s
RT
TABLE 16.1
Steady-state output voltage expressions for stable period-1 operation [5]. s
denotes the steady-state duty ratio and T denotes the switching period.
the pulse-width modulator via some typical feedback compensation conguration. For current-mode control, an inner current loop is used in addition to the
voltage feedback loop, the aim of which is to force the peak inductor current
to follow a reference signal which is derived from the output voltage feedback
loop. The result of current-mode control is a faster response. This kind of
control is mainly applied to boost and buck-boost converters which suer from
an undesirable non-minimum phase response [5, 6]. The simplied schematics
are shown in Fig. 16.2.
16.3
Conventional Treatments
As mentioned previously, power electronics circuits are essentially piecewiseswitched circuits [7]. The number of possible circuit topologies is usually xed,
and the switching is done in a cyclic manner (but not necessarily periodically
because of the feedback action). This results in a nonlinear time-varying operating mode, which naturally demands the use of nonlinear methods for analysis
and design. Indeed, researchers and engineers who work in this eld are always dealing with nonlinear problems and have attempted to explore methods
not normally used in other circuit design areas, e.g., state-space averaging [4],
phase-plane trajectory analysis [8], Lyapunov based control [9], Volterra series
approximation [10], etc. However, in order to expedite the design of power
electronics systems, adequate simplifying models are imperative. In the process of deriving models, accuracy is often traded o for simplicity for many
good practical reasons. Since closed-loop stability and transient responses are
basic design concerns in practical power electronics systems, models that can
permit the direct application of conventional frequency-domain approaches will
present obvious advantages. Thus, much research in modelling power electron-
323
ics circuits has been directed towards the derivation of a linear model that is
appropriate to a frequency-domain analysis; the limited validity being the price
to pay. (The fact that most engineers are trained to use linear methods is also
a strong motivation for developing linearized models.) For example, the averaging approach [4], one of the most widely adopted modelling approaches for
switching converters, initially yields simple nonlinear models that contain no
time-varying parameters and hence can be used more conveniently for analysis
and design. Essentially, an averaged model discards the switching details and
focuses only on the envelope of the dynamical motion. This is well suited to
characterize power electronics circuits in the low-frequency domain.
In practice, moreover, such so-called averaged models are often linearized
to yield linear time-invariant models that can be directly studied in a standard
Laplace transform domain or frequency domain, facilitating design of control
loops and evaluation of transient responses in ways that are familiar to practitioners.
16.4
Power electronics engineers frequently encounter such phenomena as subharmonic oscillations, jumps, quasi-periodic operations, sudden broadening of
power spectra, bifurcations and chaos, despite not knowing what causes them
[11, 12]. Most power supply engineers would have experienced bifurcation phenomena and chaos in switching regulators when some parameters (e.g., input
voltage and feedback gain) are varied, but usually do not examine the phenomena in detail. The usual reaction of the engineers is to avoid these phenomena
by adjusting component values and parameters, often through some trial-anderror procedure. Thus, the phenomena remain somewhat mysterious and rarely
examined in a formal manner.
What has been said so far may be regarded as a kind of stereotypical development in application-driven disciplines. Indeed, if linear models can be
protably used in design (and as long as the restricted validity of the models does not adversely compromise the design integrity), there seems to be no
immediate needs for investigating such nonlinear phenomena as chaos and bifurcation. However, as the eld of power electronics gains maturity and as the
demand for better functionality, reliability and performance of power electronics systems increases, in-depth analysis into nonlinear behavior and phenomena
becomes justiable and even mandatory. On the one hand, the study of nonlinear phenomena oers the opportunity of rationalizing the commonly observed
behavior. Thus, knowing how and when chaos occurs, for instance, will certainly help avoid it, if avoiding it is what the engineers want. On the other
hand, many previously unused nonlinear operating regimes may be protably
324
exploited for useful engineering applications, provided that such operations are
thoroughly understood. For these reasons, the study of bifurcations and chaos
in power electronics has recently attracted much attention from both the power
electronics and the circuits and systems communities.
In the next section, we present a chronological survey of the recent ndings
in the identication, analysis and modelling of nonlinear phenomena in power
electronics circuits and systems.
16.5
The occurrence of bifurcations and chaos in power electronics was rst reported
in the literature in the late eighties by Hamill et al. [13, 14]. Experimental
observations regarding boundedness, chattering and chaos were also made by
Krein and Bass [15] back in 1990. Although these early reports did not contain
any rigorous analysis, they seriously pointed out the importance of studying
the complex behavior of power electronics and its likely benets for practical
design. Since then, much interest has been aroused in the power electronics
and circuits research communities in pursuing formal studies of the complex
phenomena commonly observed in power electronics.
In 1990, Hamill et al. [16] presented a paper at the IEEE Power Electronics Specialists Conference, reporting an attempt to study chaos in a simple
buck converter which became a subject of intensive research in the following decade (and still is!). Using an implicit iterative map, the occurrence
of period-doublings, subharmonics and chaos in a simple buck converter was
demonstrated by numerical analysis, PSPICE simulation and laboratory measurement. The derivation of a closed-form iterative map for the boost converter
under a current-mode control scheme was presented later by the same group of
researchers [17]. This closed-form iterative map allowed the analysis and classication of bifurcations and structural instabilities of this simple converter.
Since then, a number of authors have contributed to the identication of bifurcation patterns and strange attractors in a wider class of circuits and devices of
relevance in power electronics. Some key publications are summarized below.
See also [1820] for alternative reviews.
The occurrence of period-doubling cascades for a simple dc/dc converter
operating in DCM was reported in 1994 by Tse [21, 22]. By modelling the
dc/dc converter operating in DCM as a rst-order iterative map, the onset of
period-doubling bifurcations can be located analytically. The idea is based on
evaluating the Jacobian of the iterative map about the xed point corresponding to the solution undergoing the period-doubling, determining the condition
for which a period-doubling bifurcation occurs (i.e., Jacobian equals 1). Simulations and laboratory measurements have conrmed the ndings. Formal
325
1.2
1.1
0.9
0.9
0.7
0.8
0.7
0.6
1
0.8
0.6
0.5
0.5
0.4
0.4
0.3
0.3
0.2
0.4
0.2
0.5
0.6
0.7
0.8
0.9
0.1
0.3
0.4
0.5
0.6
0.7
Iref
(a)
0.8
0.9
1.1
1.2
Iref
(b)
FIGURE 16.3
Bifurcation diagrams from a current-mode controlled boost converter with
L = 1.5 mH, R = 40 and T = 100 s. For (a), T /CR = 0.125, and for (b),
T /CR = 0.625.
326
7.8
10.5
10
7.6
9.5
7.4
7.2
8.5
7
8
6.8
7.5
6.6
6.4
6.5
6.2
0.18
0.2
0.22
0.24
0.26
0.28
6
0.2
0.3
0.25
0.3
0.35
0.4
0.45
0.5
0.7
0.8
(a)
(b)
12
16
15
11
14
13
12
10
9
11
10
9
8
7
6
6
0.2
0.25
0.3
0.35
0.4
0.45
0.5
0.55
0.6
5
0.2
0.3
0.4
0.5
(c)
0.6
(d)
FIGURE 16.4
Trajectories from a current-mode controlled boost converter. (a) Stable
period-1 operation; (b) stable period-2 operation; (c) stable period-4 operation;
and (d) chaotic operation.
Modelling Strategies
327
When external clocks are absent and the system is free-running, for example, under a hysteretic control scheme, the system is autonomous and does
not have a xed switching period. Such free-running converters were indeed
extremely common in the old days when xed-period integrated-circuit controllers were not available. For this type of autonomous converters, chaos cannot occur if the system order is below three. A representative example is the
free-running Cuk
converter which has been shown by Tse et al. [37] to exhibit
Hopf bifurcation and chaos.
Power electronics circuits other than dc/dc converters have also been examined in recent years. Dobson et al. [38] reported switching time bifurcation
of diode and thyristor circuits. Such bifurcation manifests as jumps in the
switching times. Bifurcation phenomena from induction motor drives were reported separately by Kuroe [39] and Nagy et al. [40]. Finally, some attempts
have been made to study higher order parallel-connected systems of converters
which are becoming popular design choice for high current applications [41].
16.6
Modelling Strategies
Cuk
converter, Hopf bifurcation that gives birth to a limit cycle consisting of
many switching periods (low-frequency behavior) can be clearly observed from
a suitable averaged model [37].
An eective approach for modelling power electronics circuits with a high degree of exactness is to use appropriate discrete-time maps obtained by uniform
or non-uniform sampling of the system states. Essentially the aim is to derive
an iterative function that expresses the state variables at one sampling instant
in terms of those at an earlier sampling instant. In what follows, the basic
concepts of discrete-time maps and how they can be used to explore nonlinear
phenomena in power electronics are explained.
328
16.6.1
Discrete-time maps
Modelling Strategies
329
(stroboscopic instants). Note that the system states are sampled at each
stroboscopic instant irrespectively of whether the system conguration
switches or not at that instant.
S-switching map (synchronous switching): obtained by sampling the system states at those stroboscopic time instants when the system commutes
from phase 2 to phase 1.
A-switching map (asynchronous switching): obtained by sampling the
system states at time instants within each ramp cycle when the system
commutes from phase 1 to phase 2.
In order to derive the relevant maps, we introduce the following two simplifying assumptions. These assumptions can be later removed.
1. No more than one commutation takes place during each period of the
modulating signal.
2. The commutation from phase 2 to phase 1 can only take place at time
instants which are multiples of the ramp cycle T .
Note that the second assumption is always true when the converter is under
current-mode control, but is not so for voltage-mode control.
16.6.2
The stroboscopic map is the most widely used type of discrete-time maps for
modelling dc/dc converters. It can be obtained by sampling the system dynamics every T seconds, say, at the beginning of each ramp cycle. To avoid
confusing with other types of maps, we use sux k as the counting index for
this map, as shown in Fig. 16.5. Applying a simple iterative procedure to the
solutions of the state equations for phases 1 and 2, the stroboscopic map can
obtained as
xk+1 = N2 (1 k )N1 (k )xk + [N2 (1 k )M1 (k ) + M2 (1 k )] E, (16.6.2)
where xk denotes x(kT ), k is the duty ratio in the kth period, and the solution
to the state equation in each phase is given by x() = Ni () + Mi ()E with
Ni () = eAi T
(16.6.3)
Ai T
Mi () = A1
I)Bi .
i (e
(16.6.4)
Moreover, under certain drastic control conditions, the converter may stay in
phase 1 or phase 2 for the entire period. In such cases (normally called skipped
cycles), one must assume k = 1 or k = 0, as appropriate.
330
vc (t) 6
A-switching
?
e u
e
tm+2
tm+1
tm
?
e
?
e
k+1 = 1 m+1 T-
k T -
tk
6
tn
6
tk+1
tk+2
tk+3
6stroboscopic
tn+1
tn+2
S-switching
FIGURE 16.5
Two iterations of the stroboscopic, S-switching and A-switching maps under
voltage-mode control.
(16.6.5)
Modelling Strategies
331
The A-switching map is obtained when the state vector is asynchronously sampled at switching times internal to the modulating period, as shown in Fig. 16.5.
Specically we can readily write down the A-switching map by noting that the
converter enters phase 2 after an A-switching, i.e.,
xm+1 = N1 (m+1 + 1 )N2 (1 m + 2 )xm
(16.6.6)
(16.6.7)
gT xm+1 = + m+1 T.
(16.6.8)
(16.6.9)
where T = gT /(T ) and a = /(T ). Note that (16.6.9) has a closed form,
although it is an implicit map.
332
16.7
System formulation
As discussed earlier, because of their switching nature, power electronics systems are often modelled by sets of ordinary dierential equations (ODEs) or
maps which are intrinsically piecewise-smooth (PWS). Specically, whenever
the circuit under investigation switches to a dierent conguration (for example, from ON to OFF, or vice versa), the vector eld of the corresponding
model changes from one functional form to another. Thus, mathematically,
power electronics circuits can be modelled by dynamical systems of the form:
...
333
Bifurcation possibilities
Power electronics systems can exhibit standard bifurcations such as perioddoubling or saddle-node. Such bifurcations are indeed frequently observed both
analytically and experimentally, as surveyed earlier in Sec. 16.5. Nevertheless,
some of the most common dynamical transitions observed in power electronics circuits, such as the sudden jump to chaos mentioned earlier, cannot be
explained in terms of standard bifurcations [4345]. In fact, power electronics
systems, being switched dynamical systems, are known to exhibit an interesting
class of bifurcations which cannot be observed in their smooth counterparts.
Specically, for switched dynamical systems, a dramatic change of the system
behavior is usually observed when a part of the system trajectory hits tangentially one of the boundaries between dierent regions in phase space. When
this occurs, the system is said to undergo a grazing bifurcation which is also
known as C-bifurcation in the Russian literature [4649].
For example, in the case of the buck converter, such an event corresponds
to the feedback signal grazing the tip of the ramp signal at a stroboscopic
point. As studied by Banerjee et al. [31] and shown further by Di Bernardo et
al. [50], it is possible to analyze the system behavior associated with grazing
by considering an appropriate piecewise linear normal-form map of the form:
334
L0
D
S
M
M0
xn
M **
D+
FIGURE 16.6
Derivation of normal-form map near a grazing bifurcation.
where A1 , A2 , B are the linear system matrices valid for describing the local behavior near the point of grazing, and is a system parameter. In fact, when the
system undergoes a grazing bifurcation, a xed point of such a normal-form map
crosses transversally some boundary in the phase space (see Fig. 16.6). This
is also called border-collision bifurcation whose occurrence in one-dimensional
and two-dimensional maps was rst reported in the western literature by Nusse
et al. [33, 34].
It is worth noting that such normal-form maps are not always piecewise
linear. In fact, it can be shown [51] that the form of the normal-form map
at a grazing depends on the discontinuity of the system vector eld, and is
piecewise linear only if the discontinuity boundary between the ON and OFF
zones is itself discontinuous such as in the cases of many power electronics
circuits. Thus, knowing the form of the normal-form map associated with a
border collision becomes important when the aim is to predict the dynamical
behavior of the system following such a bifurcation. In the next subsection,
we describe a classication method to identify the scenario following a border
collision.
16.7.3
335
Classification
We begin with a brief description of the derivation of the normal-form map associated with a grazing or border-collision bifurcation. Referring to Fig. 16.6,
suppose that as the value of some parameter increases, a periodic orbit of system (16.7.10), say L0 , becomes tangent (grazing) to the switching hyperplane,
0 , when =
. The grazing limit cycle L0 is associated with a xed point, say
M0 , on the Poincare section D. As the parameter is varied, such a xed point
moves from M to M , i.e., from a xed point associated with an orbit which
does not cross the boundary to one corresponding to a solution which crosses
the boundary. Linearizing the system ow about each of these xed points, we
can then obtain a piecewise linear map of the form (16.7.11). This approximate
map is valid if the switching hyperplane 0 is itself non-smooth [51].
An eective method for classifying and predicting the dynamical scenarios
following a border-collision bifurcation is given in Di Bernardo et al. [52]. The
available results (up to now) for the n-dimensional case are summarized as
follows. Let 1+ and 2+ be the numbers of eigenvalues, respectively, of A1
and A2 in (16.7.11), which are greater than 1. Likewise, let 1 and 2 be
the numbers of eigenvalues that are less than 1. Specically, a periodic orbit
undergoing a border collision will
smoothly change into one containing an additional section on the other
side of the switching hyperplane, if
1+ + 2+ is even;
(16.7.12)
(16.7.13)
undergo a period-doubling, if
1 + 2 is odd.
(16.7.14)
336
on the discontinuity boundary. The presence of sliding can give rise to the
formation of so-called sliding orbits, i.e., periodic solutions characterized by
sections of sliding motion (or chattering). These solutions can play an important role in organizing the dynamics of a given power electronics circuit [36].
Research is still on-going in identifying a novel class of bifurcations, called sliding bifurcations, which involve interactions between the system trajectories and
discontinuity sets where sliding motion is possible [56].
16.8
References
337
References
[1] B. K. Bose (Editor), Modern Power Electronics: Evolution, Technology, and
Applications, IEEE Press, New York, 1992.
[2] Proceedings of Annual IEEE Applied Power Electronics Conf. Exp. (APEC),
IEEE, New York, since 1986.
[3] Records of Annual IEEE Power Electronics Specialists Conf. (PESC), IEEE,
New York, since 1970.
338
References
References
339
340
References
[49] M. I. Feigin, On the structure of C-bifurcation boundaries of piecewise continuous systems, PMM, Vol. 2, pp. 820829, 1978.
[50] M. di Bernardo, C. J. Budd and A. R. Champneys, Corner collision implies
border collision, Physica D, Vol. 154, pp. 171194, 2001.
[51] M. di Bernardo, C. J. Budd and A. R. Champneys, Grazing and border-collision
in piecewise-smooth systems: a unied analytical framework, Physical Review
Letters, Vol. 86, pp. 25532556, 2001.
[52] M. di Bernardo, M. I. Feigin, S. J. Hogan and M. E. Homer, Local analysis of
C-bifurcations in N -dimensional piecewise smooth dynamical systems, Chaos,
Solitons and Fractals, Vol. 10, No. 11, pp. 18811908, 1999.
[53] M. I. Feigin, The increasingly complex structure of the bifurcation tree of a
piecewise-smooth system, J. Appl. Math. Mech., Vol. 59, pp. 853863, 1995.
[54] S. Banerjee and C. Grebogi, Border collision bifurcations in two-dimensional
piecewise smooth maps, Physical Review E, Vol. 59, pp. 40524061, 1999.
[55] L. E. Nusse, E. Ott and J. A. Yorke, Border collisions bifurcations an explanation for observed bifurcation phenomena, Physical Review E, Vol. 49,
pp. 10731076, 1994.
[56] M. di Bernardo, K.H. Johansson and F. Vasca, Self-oscillations in relay feedback
systems: symmetry and bifurcations, Int. J. Bifur. Chaos, Vol. 11, April 2001.
[57] T. Kousaka, T. Ueta and H. Kawakami, Bifurcation in switched nonlinear dynamical systems, IEEE Trans. Circ. Syst. II, Vol. 46, No. 7, pp. 878885, July
1999.
[58] M. di Bernardo, The complex behaviour of switching devices, IEEE CAS
Newsletter, Vol. 10, No. 4, pp. 113, December 1999.
[59] F. Ueno, I. Oota and I. Harada, A low-noise control circuit using Chuas circuit
for a switching regulator, Proc. European Conf. Circuit Theory and Design,
pp. 11491152, 1995.
[60] J. H. B. Deane and D. C. Hamill, Improvement of power supply EMC by chaos,
Electronics Letters, Vol. 32, No. 12, p. 1045, June 1996.
[61] P. J. Aston, J. H. B. Deane and D. C. Hamill, Targeting in systems with
discontinuities, with applications to power electronics, IEEE Trans. Circ. Syst.
I, Vol. 44, No. 10, pp. 10341039, October 1997.
[62] C. Batlle, E. Fossas and G. Olivar, Stabilization of periodic orbits of the buck
converter by time-delayed feedback, Int. J. Circuit Theory Appl., Vol. 27,
pp. 617631, 1999.