Design and Implementation of Coupled Inductor Cuk Converter Operating in Continuous Conduction Mode
Design and Implementation of Coupled Inductor Cuk Converter Operating in Continuous Conduction Mode
.......................................................................................................... 136
4.3.4 Selection of
.......................................................................................................... 142
4.3.5 Design of Coupled-Inductor ...................................................................................... 147
5. CONTROLLER DESIGN .................................................................................................... 175
5.1 INTRODUCTION ............................................................................................................ 175
5.2 CONTROL BLOCK DIAGRAM ...................................................................................... 178
5.3 CONTROLLER DESIGN ................................................................................................. 179
5.3.1 Controller Design of the Current Loop ....................................................................... 180
5.3.2 Controller Design of the Voltage Loop ...................................................................... 187
5.4 DOMAIN CONVERSIONS .............................................................................................. 192
5.5 APPLICATION SPECIFIC POINTS ................................................................................. 194
6. EXPERIMENTAL RESULTS .............................................................................................. 197
6.1 INTRODUCTION ............................................................................................................ 197
6.2 OPEN-LOOP RESULTS .................................................................................................. 199
6.2.1 Voltage Waveforms .................................................................................................. 199
6.2.2 Current Waveforms ................................................................................................... 206
x
6.2.3 Efficiency of the Converter at Different Operating Conditions ................................... 221
6.3 CLOSED-LOOP RESULTS.............................................................................................. 222
6.3.1 Constant Input Current Mode .................................................................................... 225
6.3.2 Constant Input Power Mode ...................................................................................... 236
7. SUMMARY AND CONCLUSIONS ..................................................................................... 246
REFERENCES ......................................................................................................................... 250
xi
LIST OF TABLES
TABLES
Table 4-1 Voltage and current stresses on the switch .......................................... 132
Table 4-2 Power dissipation on the switch at end operating points ...................... 134
Table 4-3 Voltage and current stresses on the diode............................................ 134
Table 4-4 Power dissipation on the diode at end operating points ....................... 136
Table 4-5 Voltage ripple consideration of
...................................................... 139
Table 4-6 Power dissipation on
...................................................... 145
Table 4-8 Power dissipation on
and
and
and
and
and
and
waveforms in IMC,
waveform in IMC,
waveform in IMC,
waveform in IMC,
waveform in IMC,
, full load.......................... 60
Figure 2-38
and
waveforms in IMC,
in Matlab, ........................... 83
xiv
Figure 3-5: Response of
in Matlab, ........................... 85
Figure 3-7: Response of
in Matlab,
.............................. 86
Figure 3-9: Response of
....................................................................................................... 87
Figure 3-10: Step response of
in Matlab,
............................. 88
Figure 3-11: Response of
....................................................................................................... 88
Figure 3-12: Step response of
in Matlab,
.......................... 89
Figure 3-13: Response of
in Simplorer, ideal
elements,
..................................................................................... 90
Figure 3-14: Step response of
in Matlab,
........................... 91
Figure 3-15: Response of
..................................................................................................... 91
Figure 3-16 Circuit schematic of CIC with parasitic elements and directions of
currents included .............................................................................................. 92
Figure 3-17 Equivalent circuit schematic of CIC with parasitic elements
(directions of currents included) in Mode 1 ...................................................... 93
Figure 3-18: Equivalent circuit schematic of CIC with parasitic elements in Mode
2 ....................................................................................................................... 96
Figure 3-19: Step response of
in Matlab,
.......................... 113
xv
Figure 3-24: Response of
in Simplorer, parasitic
elements,
..................................................................................... 114
Figure 3-25: Step response of
in Matlab,
........................... 115
Figure 3-26: Response of
in Simplorer, parasitic
elements,
..................................................................................... 115
Figure 3-27: Step response of
in Matlab,
........................ 116
Figure 3-28: Response of
in Simplorer, parasitic
elements,
................................................................................... 117
Figure 3-29: Step response of
in Matlab,
......................... 118
Figure 3-30: Response of
in Simplorer, parasitic
elements,
................................................................................... 118
Figure 4-1 Direct connection to vehicle battery,
........................... 124
Figure 4-2 Direct connection to vehicle battery,
............ 125
Figure 4-3 Connection to vehicle battery via cigarette lighter adapter,
and
......................................................................................... 126
Figure 4-4 Connection to vehicle battery via cigarette lighter adapter,
and
......................................................................................... 127
Figure 4-5 Input current waveform of a basic uk converter with 40% peak-to-peak
ripple .............................................................................................................. 128
Figure 4-6 Simple capacitor model ..................................................................... 136
Figure 4-7
........................................................................................................ 141
Figure 4-8 Current consideration of
, a)
waveform, b)
waveform, c)
waveform ....................................................................................................... 142
Figure 4-9
........................................................................................................ 146
Figure 4-11
and
and
and
................................................................. 182
Figure 5-5 Bode plot of
......................................................................... 183
Figure 5-6 Bode plot of
..................................................................... 184
Figure 5-7 Bode plot of
................................................ 185
Figure 5-8 Bode plot of open loop transfer function,
.............................. 186
Figure 5-9 Bode plot of
.......................................................................... 187
Figure 5-10 Control block diagram of the voltage loop ....................................... 188
Figure 5-11 Pole-zero map of
.............................................................. 189
xvii
Figure 5-12 Pole-zero map of
................................................................. 191
Figure 5-13 Bode plot of
........................................................................ 191
Figure 5-14 Bode plot of
........................................................................ 192
Figure 5-15 Control block diagram with the limiters .......................................... 194
Figure 6-1 Top view of the implemented circuit ................................................. 198
Figure 6-2 Bottom view of the implemented circuit ............................................ 198
Figure 6-3
waveform at
waveform at
waveform at
waveform at
waveform at
waveform at
waveform at
waveform at
waveform at
and full
load ................................................................................................................ 205
Figure 6-13
waveform at
and full
load ................................................................................................................ 206
Figure 6-15
and full load ................................................................................................... 208
Figure 6-17
and full load ................................................................................................... 209
Figure 6-19
and 50% load.................................................................................................. 211
Figure 6-21
and 50% load.................................................................................................. 212
Figure 6-23
and 10% load.................................................................................................. 213
Figure 6-25
and full load ................................................................................................... 215
Figure 6-28
and full load ................................................................................................... 216
Figure 6-30
and 50% load.................................................................................................. 218
Figure 6-32
and 50% load.................................................................................................. 219
Figure 6-34
and 10% load.................................................................................................. 221
Figure 6-37 Noise-sensitive (left) and noise-immune (right) measurement
techniques ...................................................................................................... 223
Figure 6-38 Results of noise-sensitive (left) and noise-immune (right)
measurements ................................................................................................. 224
Figure 6-39 Output voltage waveforms in ac coupling mode at 1msec/div (left) and
2.5msec/div time scale.................................................................................... 225
Figure 6-40 Timing diagram of the control loop in CICM .................................. 226
xix
Figure 6-41
(green) and
.............................................................................................. 228
Figure 6-42 Pole-zero map of
............................................................... 229
Figure 6-43
and
229
Figure 6-44
and
....................................................................................................................... 230
Figure 6-45
and
....................................................................................................................... 231
Figure 6-46
and
.............................................................................................. 232
Figure 6-47
and
.............................................................................................. 232
Figure 6-48
and
........................................................................................... 233
Figure 6-49
and
........................................................................................... 234
Figure 6-50
and
.............................................................................................. 235
Figure 6-51
and
.............................................................................................. 235
Figure 6-52
and
......................................................................................... 236
Figure 6-53
and
......................................................................................... 236
Figure 6-54 Timing diagram of the control loop for CIPM ................................. 237
Figure 6-55
and
238
Figure 6-56
and
239
Figure 6-57
and
....................................................................................................................... 239
Figure 6-58
and
....................................................................................................................... 240
xx
Figure 6-59
and
........................................................................................... 240
Figure 6-60
and
........................................................................................... 241
Figure 6-61
and
........................................................................................... 242
Figure 6-62
and
........................................................................................... 242
Figure 6-63
and
........................................................................................... 243
Figure 6-64
and
........................................................................................... 243
Figure 6-65
and
......................................................................................... 244
Figure 6-66
and
......................................................................................... 244
xxi
LIST OF ABBREVIATIONS
ADC Analog to Digital Conversion
CCM Continuous Conduction Mode
CIC Coupled-Inductor uk Converter
CICM Constant Input Current Mode
CIPM Constant Input Power Mode
DCM Discontinuous Conduction Mode
ESR Equivalent Series Resistance
IMC Integrated Magnetic uk Converter
LVR Linear Voltage Regulator
MMF Magneto Motive Force
PI Proportional Integral
PWM Pulse Width Modulation
SMPS Switched-Mode Power Supply
SSAM State-Space Averaging Method
xxii
NOMENCLATURE
) / second (
winding (Turns)
winding (Weber)
and
windings
xxiii
) /
secondary (
) winding
The summation of the core reluctance and the air gap if it exists
and
Self inductance of
winding
and
windings
Leakage inductance of
winding
side
Self inductance of
winding
and
windings (H)
and
windings
Leakage inductance of
winding
Measured inductance of
winding while
winding is open-
circuited
Measured inductance of
winding while
winding is short-
circuited
Measured inductance of
winding while
winding is open-
circuited
Measured inductance of
winding while
winding is short-
circuited
or
or
or
or
@20
and
and
is shown in Figure 2-2.
[16, 27]
Figure 2-2 Reluctance model of a coupled-inductor
[16, 27]
Transformers have almost the same reluctance model with the only difference in the
respective polarities of MMFs
and
.
As mentioned earlier, both electrical and magnetic circuits will be analyzed
here. Voltage mesh equations will be utilized on both circuits. In magnetic circuit,
node equations will also be needed.
The analysis of the converter will be made in continuous conduction mode,
CCM in which there exist two modes: Mode 1 and Mode 2. The first applies to
analysis within a time interval in which the switch, Q is in conduction, briefly
referring as state on and the diode, D not conducting, namely state off. Mode 2
applies to analysis into the complementary time interval of Mode 1, i.e. the time
interval in which Q is off and D is on. As expected, equations governing the
9
behaviour of the converter electrical circuit in both modes differ from each other.
Hence, electrical equations are grouped into two: Mode 1 and Mode 2. Unlike
electrical equations, however, magnetic equations are same for both modes.
2.2.1.1 Electrical Mesh Equations
The following meshes are active in the equivalent circuit for CIC operating
steadily in Mode 1. The superscripts elc or mag used in loop names are for the
purpose of differentiation of the same names used in electrical and magnetic loops,
respectively.
Figure 2-3 Active meshes in Mode 1 in the equivalent electrical circuit for CIC
Recalling that the instantaneous current,
(2-1)
Since a coupled-inductor is a kind of transformer rather than an inductor, it is better
to use the simplest relationship between flux and voltage given by Faradays law of
induction as;
(2-2)
10
Throughout the analyses, voltages of the windings are replaced by the help of this
law.
In Mode 1, (2-3)-(2-5) are the electrical mesh equations for the equivalent
circuit:
(2-3)
(2-4)
(2-5)
In Mode 2, active meshes in the equivalent circuit for CIC are shown in
Figure 2-4.
Figure 2-4 Active meshes in Mode 2 in the equivalent electrical circuit for CIC
Electrical mesh equations for this mode of operation are given as (2-6)-(2-8):
(2-6)
(2-7)
(2-8)
11
2.2.1.2 Magnetic Mesh and Node Equations
The magnetic circuit for a coupled-inductor is shown in Figure 2-5. Mesh and
node equations are given in (2-9)-(2-13) for such a circuit.
Figure 2-5 Meshes and nodes in the magnetic equivalent circuit for a coupled-
inductor
(2-9)
(2-10)
(2-11)
(2-12)
(2-13)
and
(2-14)
(2-15)
(2-16)
12
(2-17)
(2-18)
Using (2-14), (2-16) and (2-17):
(2-19)
Using (2-15), (2-18) and (2-17):
(2-20)
2.2.1.3 Combination of Electrical and Magnetic Equations
Derivatives of the fluxes are obtained in terms of the derivatives of the
currents. Now, these results can be reflected to electrical mesh equations. As
mentioned earlier, magnetic equations are independent of the switch state. Hence (2-
19) and (2-20) are used for both electrical mesh equation sets.
In Mode 1, (2-21)-(2-23) are the magnetic mesh equations for the equivalent
circuit:
Using (2-3) and (2-19):
(2-21)
Using (2-4) and (2-20):
(2-22)
Using (2-5):
(2-23)
In Mode 2, (2-24)-(2-26) are the magnetic mesh equations for the equivalent
circuit:
13
Using (2-6) and (2-19):
(2-24)
Using (2-7) and (2-20):
(2-25)
Using (2-8):
(2-26)
Note that all the parameters are in terms of electrical quantities while taking into
account the magnetic connections and parameters.
In this work, state-space averaging method will be used in order to get
steady-state and dynamic model of the circuit. State variables are inductor currents
and capacitor voltages. Derivatives of capacitor voltages will be obtained in the next
chapter. For the derivatives of inductor currents, however, basic inductor equation
(2-1) cannot be used. Since a coupled-inductor is utilized, derivatives of inductor
currents should be obtained in terms of electrical quantities. Further simplification of
(2-21) and (2-22) in Mode 1 and (2-24) and (2-25) in Mode 2 gives the necessary and
ready-to-use information for the derivatives of the inductor currents. (2-23) and (2-
26) are written up to now just for mesh completeness and will not be used anymore.
Further simplification of (2-21) and (2-22) in Mode 1 gives:
(2-27)
(2-28)
where
(2-29)
(2-30)
14
(2-31)
(2-27) and (2-28) shows the representations of the derivatives of the inductor
currents in terms of the state-variables, namely
and
, in Mode 1. Note that they do not include any flux-variable. When similar
simplification is applied to (2-24) and (2-25) in Mode 2, the following equations are
obtained:
(2-32)
(2-33)
At this point (2-27)-(2-28) and (2-32)-(2-33) are left as inputs to the next chapter.
This chapter continues with the merits of using coupled-inductor in basic uk
converter.
2.2.2 Merits of Coupling
In this section, by the help of some reasonable assumptions, merits of
coupling will be proven theoretically. Those merits are either ripple-free input
current waveform or ripple-free output current waveform. (2-21)-(2-22) and (2-24)-
(2-25) are the input equations to this section. Then, some simplifications and
assumptions are needed in order to proceed. These assumptions are either mentioned
explicitly or implied in [8, 16, 26].
Assumption 2-1: Input voltage is constant.
(2-34)
Assumption 2-2: Voltage on the output capacitor,
(2-35)
15
Simplification 2-1: Mean value of the voltage on
.
[
2
]
This
fact, which is obtained in CHAPTER 3, is also explained below in another way.
When the outermost loop voltage equation starting with
in the clockwise
direction- is written, the following equation is obtained.
(2-36)
As known, the mean values of the inductors voltages are zero at steady-state.
(2-37)
(2-38)
Therefore, taking the mean values of both sides of the equation gives that:
(2-39)
Negligible ripple assumption permits this simplification be utilized in the electrical
mesh equations.
Assumption 2-3: Voltage on
capacitor,
(2-40)
As a result of these assumptions and the simplification, the following
equation sets are obtained.
In Mode 1:
Using (2-21) and (2-34):
(2-41)
16
Using (2-22), (2-39), (2-34) and (2-35):
(2-42)
In Mode 2:
Using (2-24), (2-39), (2-34) and (2-35):
(2-43)
Using (2-25) and (2-35):
(2-44)
Following two sections present the conditions for ripple-free current either at
input or output and also the simulation results that verify the theoretical analyses.
2.2.2.1 Ripple-Free Input Current
The requirement is to obtain ripple-free input current. Input current is the
same as the current of
or simply
(2-45)
Since it is constant, its derivative is zero.
(2-46)
When this condition is applied to (2-41)-(2-44), some restrictions arises. If those
restrictions turn out to be applicable, it means that zero-ripple input current
waveform is possible to obtain. Moreover, the restrictions should be applicable for
both modes.
17
In Mode 1:
Using (2-41) and (2-46):
(2-47)
Using (2-42) and (2-46):
(2-48)
In Mode 2:
Using (2-43) and (2-46):
(2-49)
Using (2-44) and (2-46):
(2-50)
For the input current to be ripple-free in Mode 1, the following condition has
to be satisfied.
Using (2-47) and (2-48):
(2-51)
The same condition is also required for the input current to be ripple-free in Mode 2.
Using (2-49) and (2-50):
(2-52)
Since the same condition makes the input current ripple-free in both modes, namely
the entire period, it can be proven that applying this condition satisfies to obtain a
ripple-free input current. (2-52) can be simplified into the following form.
18
(2-53)
Since
, equivalent reluctance of
turns out to be
a little bit smaller than
is applied to
inductors in Mode 1 and
and
generate
and
Simulation is performed with the ideal electrical components. In other words;
the switch, the diode and the capacitors are assumed to be ideal.
and
are the
ones that have been selected in the design chapter in terms of capacitance values.
: Ideal switch : Ideal diode
Above values give the following parameters.
21
As adjusted beforehand,
and
current,
. Current-axis is in ampere and time-axis is in second, like all the other current-
time graphs of the simulation results throughout this work. Note that the envelope of
waveform is very thin and negligible with respect to its mean value.
Figure 2-8
and
is seen as 26.292A and with the minimum is at 26.217A. Its ripple percent is not
zero but 0.29%. Actually, this result is expected because of the assumptions made in
the analysis. Still, the result is considered as satisfactory. Besides, an oscillation to
which the switching ripple is superimposed is observed in
waveform.
As mentioned earlier, ripple-free property is independent of current and
voltage levels. To show this fact, the load resistance is increased to its 500%.
Envelopes of the resultant current waveforms are presented in Figure 2-11. As seen
in Figure 2-12,
and
. Further information
on this issue can be found in [8, 10].
2.2.2.2 Ripple-Free Output Current
This time the requirement is to obtain a ripple-free output current. The output
current is not the same as
(2-54)
In this case,
(2-55)
The same procedure is applied for ripple-free output current case. At the end, the
following condition is obtained for a ripple-free output current waveform.
Using (2-41)-(2-44) and (2-55):
(2-56)
Further simplification gives that:
(2-57)
26
Similarly, since
, equivalent reluctance of
turns
out to be a little bit smaller than
has to be a little
bit larger than
and
in ripple-free input current case is enough. That is, the same parameters will be
used, but
and
These parameters give that:
In this situation,
is almost ripple-free
and
and
ripple percent in ripple-free input current case is 6.5%,
ripple percent in
Figure 2-15 Detailed
(2-58)
In other words, if input current is desired to be ripple-free, mutual inductance
between
and
. Similarly, one-to-one
corresponding condition of (2-57) is:
(2-59)
That is to say, if output current is desired to be ripple-free, mutual inductance
between
and
.
As a physical interpretation, another explanation is given below. Consider a
coupled-inductor model in Figure 2-17. As mentioned in [10], it is T model of a
transformer at the same time. The only difference is the direction of
. While four
Figure 2-17 Equivalent coupled-inductor model
30
electrical variables (
(2-60)
(2-61)
(2-62)
(2-63)
(2-64)
(2-65)
(2-66)
By using these basic relationships, ripple-free current conditions can be obtained.
Lets consider the ripple-free input current case. For
to be ripple-free, voltage on
current hence
is equal to
(2-67)
Then, by voltage division rule, voltage on
due to
(2-68)
Since the same flux links the both windings, voltages induced on
and
are
proportional to their number of turns.
(2-69)
Using (2-69), (2-62) and (2-66):
(2-70)
As stated in (2-58),
is adjusted to be equal to
(2-71)
Voltage applied to the negative terminal of
is also ,
and
and
and
(2-72)
(2-73)
In Mode 2:
(2-74)
33
(2-75)
2.2.3.2 Magnetic Mesh and Node Equations
Magnetic circuit already includes the parasitic elements, namely the
reluctances. Hence, (2-19) and (2-20) will directly be utilized for this analysis.
2.2.3.3 Combination of Electrical and Magnetic Equations
In Mode 1, the following equation set is obtained.
In Mode 1:
Using (2-72) and (2-19):
(2-76)
Using (2-73) and (2-20):
(2-77)
Similarly, in Mode 2, the following equation set is obtained.
In Mode 2:
Using (2-74) and (2-19):
(2-78)
34
Using (2-75) and (2-20):
(2-79)
Using (2-76)-(2-79), derivatives of the two state-variables namely
and
- can
be written in terms of state-variables and input variables. Before that, some
simplification should be done as an addition to (2-29)-(2-31).
In Mode 1:
Using (2-76) and (2-77):
(2-80)
(2-81)
In Mode 2:
Using (2-78) and (2-79):
(2-82)
35
(2-83)
These values are given with respect to ripple-free input current case
simulation. As mentioned earlier,
and
and
are also
interchanged and simulation is done with those parameters. Besides, simulations are
performed at the same operating point given in ideal electrical circuit element cases.
Figure 2-20, Figure 2-21 and Figure 2-22 belong to ripple-free input current
case. As seen in Figure 2-20, both
and
and
inductor
increases, that of
and
and
and
and
and
and
isolation transformer windings on the same magnetic core. One of the recommended
structures that can be seen in [27] is given in Figure 2-27. This structure is used in
the analysis.
Figure 2-27 Proposed integrated magnetic structure
[27]
Complete reluctance model of the integrated magnetic structure is given in Figure
2-28. As mentioned earlier, it is slightly different than the proposed reluctance model
of integrated magnetic structure in [16].
43
Figure 2-28 Complete reluctance model of the proposed integrated magnetic
structure
2.3.1.1 Electrical Mesh Equations
In Mode 1, active meshes in the equivalent circuit for IMC are shown in
Figure 2-29.
Figure 2-29 Active meshes in Mode 1 in the equivalent electrical circuit for IMC
In Mode 1:
(2-84)
(2-85)
44
(2-86)
Magnetic node equations are utilized at this step. They give the following equations.
Using Figure 2-28:
(2-87)
(2-88)
Hence, the following equation set is obtained.
In Mode 1:
Using (2-84):
(2-89)
Using (2-85) and (2-87):
(2-90)
Using (2-86) and (2-88):
(2-91)
Active meshes in Mode 2 are given in Figure 2-30.
Figure 2-30 Active meshes in Mode 2 in the equivalent electrical circuit for IMC
45
In Mode 2:
(2-92)
(2-93)
(2-94)
By the replacement of
and
(2-95)
Using (2-93) and (2-88):
(2-96)
Using (2-94):
(2-97)
Assumptions utilized in the coupled-inductor analysis are also used here.
However, there are some differences stemming from the isolation transformer. (2-39)
were valid in the previous analysis. In this analysis, however, the followings are
valid.
(2-98)
(2-99)
When these simplifications together with the previous assumptions made in coupled-
inductor analysis are applied, the following equation set is obtained.
46
In Mode 1:
Using (2-89) and (2-34):
(2-100)
Using (2-90), (2-40) and (2-98):
(2-101)
Using (2-91), (2-35), (2-40) and (2-99):
(2-102)
In Mode 2:
Using (2-95), (2-34), (2-40) and (2-98):
(2-103)
Using (2-96), (2-40) and (2-99):
(2-104)
Using (2-97) and (2-35):
(2-105)
2.3.1.2 Magnetic Mesh Equations
Using the following magnetic circuit, mesh equations can be obtained.
47
Figure 2-31 Active meshes in magnetic circuit schematic of integrated magnetic
structure
Two Major Loops:
(2-106)
(2-107)
Five Minor Loops:
(2-108)
(2-109)
(2-110)
(2-111)
(2-112)
48
and
(2-113)
Note: By using (2-113) or the assumptions such as constant capacitor voltages that
will be used later, the chance of obtaining
and
in terms of state-variables
and input variables is lost intentionally. If state-space averaging method were to be
applied to IMC in this work,
and
(2-114)
(2-115)
(2-116)
(2-117)
(2-118)
49
Two Major Loops:
Using (2-106), (2-113)-(2-115):
(2-119)
Using (2-107), (2-113)-(2-115):
(2-120)
and
and
and
terms should be
expressed in the same manner. Using (2-119) and (2-120), they can be obtained
easily.
(2-121)
(2-122)
For further simplification, Mode 1 and Mode 2 can be considered separately for
magnetic loop equations. From Figure 2-29 it is clearly seen that
during in
Mode 1.
In Mode 1:
(2-123)
Similarly,
is equal to
(2-124)
Magnetic equations can now be separately written for Mode 1 and Mode 2.
50
In Mode 1:
Using (2-121) and (2-123):
(2-125)
Using (2-122) and (2-123):
(2-126)
Using (2-114)-(2-116):
(2-127)
(2-128)
(2-129)
Using (2-117) and (2-123):
(2-130)
Using (2-118) and (2-123):
(2-131)
In Mode 2:
Using (2-121) and (2-124):
(2-132)
Using (2-122) and (2-124):
(2-133)
Using (2-114) and (2-115):
(2-134)
51
(2-135)
Using (2-116) and (2-124):
(2-136)
Using (2-117):
(2-137)
Using (2-118) and (2-124):
(2-138)
Derivatives of the fluxes are obtained in terms of the derivates of the currents. Now,
these results can be reflected to electrical mesh equations. Before that, the following
notation simplification can be useful.
(2-139)
In Mode 1:
Using (2-84), (2-125) and (2-139):
(2-140)
Using (2-101), (2-125)-(2-129), (2-131) and (2-139):
(2-141)
Using (2-102), (2-125)-(2-128), (2-130)-(2-131) and (2-139):
(2-142)
(2-142) implies that
, which is possible.
52
In Mode 2:
Using (2-103), (2-132)-(2-136) and (2-138)-(2-139):
(2-143)
(2-143) implies that
(2-144)
Using (2-105), (2-133) and (2-139):
(2-145)
When analyzed, it can be seen that all four windings have the same voltage
waveform:
in Mode 1 and
(2-146)
Using (2-143) and (2-146):
(2-147)
Since
and
(2-148)
Using (2-142) and (2-148):
(2-149)
The same argument is also valid for
(2-150)
(2-151)
(2-152)
(2-153)
Although the conditions above seem to be highly different than the coupled-
inductor case, they resembles to the results obtained in that analysis much. What the
conditions (2-150)-(2-153) say can be physically explained and verified by the
54
derivations of leakage inductance, mutual inductance, magnetizing inductance, self
inductance, coupling coefficient and effective turns ratio terms for integrated
magnetic structure. The derivations are not given in this work. However, their results
will be given in order to show the similarities between coupled-inductor and
integrated magnetic structure cases. Corresponding conditions are presented below.
(2-154) dictates that input current can be made ripple-free by balancing it with the
first winding of the transformer. Similarly, (2-156) says that output current is made
ripple-free by balancing it with the second winding of the transformer.
Corresponding equations of (2-150), (2-151), (2-152) and (2-153) respectively:
(2-154)
(2-155)
(2-156)
(2-157)
2.3.1.4 Verification by Simulation
In order to verify the conditions in (2-150)-(2-153), simulation will be done.
IMC circuit implementation will not be realized in this work. However, its
parameters are selected as if it would be implemented. In other words, its ratings are
similar to the CIC designed and implemented in this work. The following
simulation model is implemented in Simplorer. As in case of coupled-inductor,
orange lines represent the magnetic connections and black lines electrical
connections. The interactions of magnetic and electrical connections are windings
that generate
55
Figure 2-32 Simulation model of IMC in Simplorer
Note that
and
and
should
be determined such that other two conditions, (2-150) and (2-152), are also satisfied.
and
After this point, as long as the above ratios are satisfied, ripple-free input and output
waveforms can be observed. Note that turns ratio of the transformer is not mentioned
yet. According to the conditions, it seems not to affect the result. Probably it only
determines the input-to-output voltage transfer ratio. Input-to-output voltage transfer
ratio is given in [12] such that:
(2-158)
56
In order to prove this foresight, at first,
is selected as equal to
and simulation
is done with this turns ratio. Secondly, a different turns ratio is selected and the
above estimation is verified by another simulation.
In coupled-inductor case,
and
and
are selected as around 10. Then, the following parameters are obtained.
Although the numbers of turns seem to be selected randomly, they are chosen as
practical and reasonable values. IMC design is not given in this work. Just an
abstract of it can be given here. In that design,
and
and
and
and
results in lower
and
and
.
Using the above numbers of turns and reluctance values, the following
simulation results are obtained. The operating conditions are the same as coupled-
inductor case.
and
. As
understood form the thin envelopes, they are almost ripple-free.
57
Figure 2-33
and
waveforms in IMC,
, full load
Detailed waveforms of
and
waveform in IMC,
, full load
58
As seen in Figure 2-34, the maximum of input current is 23.89A and
minimum is 23.87A. Hence, its peak-to-peak ripple percent is 0.08%. Its waveform
can be considered as little bit strange because, as a known fact, inductor current
cannot change instantly. Since its magnitude is extremely low, detailed observation
seems not to be necessary.
Figure 2-35 Detailed
waveform in IMC,
, full load
Detailed
is affected by input
voltage source and
, that of
is determined by
and
and
is not large
enough,
and secondary
winding,
waveform in IMC,
, full load
Primary winding current waveform is seen above. In Mode 2, its current is
the same with
value
would result in larger ripple.
Now, the independency of turns ratio of the transformer from the ripple-free
current waveforms property will be verified by simulation. All parameters and
operating conditions except
is
60
selected as
. Then, according to
Figure 2-37 Detailed
waveform in IMC,
, full load
Simulation results are given in the following figures. As seen in Figure 2-38, input
and output current waveforms are still almost ripple-free. Since it is very obvious,
their zoomed in waveforms are not given here. By the help of this figure, the
independency of ripple-free current waveforms from turns ratio and also (2-158) is
proven. In this situation, ripple content of the transformers primary winding is
expected to be higher. This fact can even be understood by visual comparison of
Figure 2-39 with Figure 2-36. This time, current rises from 19.38A to 22.09A at the
top of the waveform. Then, its ripple percent is calculated as 13.07%. This value is
8.66% in
and
waveforms in IMC,
, full load
Figure 2-39 Detailed
waveform in IMC,
, full load
62
CHAPTER 3
STEADY-STATE AND DYNAMIC MODEL ANALYSIS
3.1 INTRODUCTION
In this chapter; steady-state model, dynamic model and transfer functions of
coupled-inductor uk converter, CIC in continuous conduction code, CCM are
obtained utilizing state-space averaging method, SSAM. SSAM is another
contribution of uk to power electronics. Detailed explanation about it can be found
in [1]. For the sake of completeness, a summary of SSAM is given in this chapter.
There, the steady-state model, dynamic model and transfer functions are derived in
matrix forms. Therefore, once the state-space model of a circuit is obtained, by
placing the matrices into the final equations give the steady-state model, dynamic
model and transfer functions of the circuit. This process is implemented for the
circuit with ideal and parasitic elements. In the previous chapter, derivatives of two
state-variables namely
and
(3-4)
Capitals such as and are constant terms and represent the dc levels at that
operating point. Terms with ~ such as and are time-dependent and
represent the small signal part at that operating point. As a first step, (3-3)-(3-4) are put
into (3-1)-(3-2).
(3-5)
(3-6)
and rearranging them gives:
(3-7)
(3-8)
At this point, the following assumption highly simplifies the analysis. Actually, it can
be considered as a reasonable assumption for small signal analysis.
65
Assumption 3-1: Magnitudes of the capitals are much larger than that of the terms
with ~.
Hence, neglecting second-order small-signal terms in (3-7)-(3-8) gives;
(3-9)
(3-10)
in which steady-state and dynamic models can be separated as;
Steady-State Model:
(3-11)
(3-12)
Dynamic Model:
(3-13)
(3-14)
As mentioned earlier, transfer functions of the system can be obtained by
applying Laplace transformation to dynamic model. Since the terms with ~
represent the small signal ac part, initial conditions can be taken as zero. Actually,
they are included in dc parts.
(3-15)
(3-16)
At this step, can be obtained by using (3-15).
66
(3-17)
Then, is obtained in terms of the known quantities and inputs by using (3-16)
and (3-17) as;
(3-18)
3.3 ANALYSIS OF COUPLED-INDUCTOR UK CONVERTER,
CIC WITH IDEAL ELEMENTS
In this section, steady-state and dynamic model of CIC in continuous
conduction mode, CCM are obtained. Using dynamic model, transfer functions are
also derived. In [2], all these derivations are present for basic uk converter. As
expected, they are highly different than those of CIC. Therefore, it is necessary to
derive them for coupled-inductor case. Besides, in [2], load namely the output
current- is not considered as a time varying input to the converter. Hence, it is not
possible to obtain a transfer function reference to the output current there. Here, the
output current is considered as a time varying input to the converter.
The circuit schematic of coupled-inductor uk converter with ideal elements
is seen in Figure 3-1. Directions of the currents are also shown.
Figure 3-1 Circuit schematic of CIC with ideal elements and directions of currents
included
67
Since there are four energy storage elements, four state-variables exist. They are
currents in inductors and voltages in capacitors. State-space equations are written in
terms of these variables, inputs and outputs. In this work, input current and output
voltage are chosen as outputs.
3.3.1 State-Space Equation Set in Mode 1
The equivalent circuit in Mode 1 is as shown in Figure 3-2, where the switch
is on.
Figure 3-2: Equivalent circuit schematic of CIC with ideal elements in Mode 1
Burrowing the derivatives of
and
(3-19)
(3-20)
Now, derivatives of capacitor voltages should be obtained.
For
(3-21)
68
(3-22)
(3-23)
.
For
(3-24)
(3-25)
(3-26)
Using (3-19), (3-20), (3-23) and (3-26); the following equation matrix is obtained.
(3-27)
As it is seen clearly, input voltage and output current are inputs to the converter.
Duty-factor of the switch is defined at this stage. Symbol represents the duty-
factor.
(3-28)
As mentioned earlier, multiplying the matrices with their durations gives their
contribution to the averaged matrices. Matrix or simply for Mode 1 is defined
as
(t) or simply
(3-29)
(3-30)
is written similarly. As mentioned earlier, outputs of the system are selected as
the output voltage and the input current, which equals
(3-31)
The following matrices are obtained.
(3-32)
(3-33)
3.3.2 State-Space Equation Set in Mode 2
The equivalent circuit in Mode 2 is shown in Figure 3-3, with the switch OFF
and the diode ON. Similar arguments followed for
and
(3-34)
(3-35)
70
Figure 3-3: Equivalent circuit schematic of CIC with ideal elements in Mode 2
For
(3-36)
(3-37)
For
:
It is the same with that in Mode 1.
(3-38)
Using (3-34), (3-35), (3-37) and (3-39); the following equation matrix is obtained.
(3-39)
This time, the normalized duration with respect to the switching period is .
Hence, we have;
71
(3-40)
(3-41)
is the same with that in Mode 1. Likewise, the following matrices are obtained.
(3-42)
(3-43)
3.3.3 Averaging of Matrices in Mode 1 and Mode 2
Averaged matrices are represented by
and
in the following manner.
(3-44)
(3-45)
(3-46)
72
(3-47)
Thus, state-space averaged model of the circuit can be represented as;
(3-48)
(3-49)
3.3.4 Decomposition of Averaged Model into Steady-State and
Dynamic Models
When (3-3), (3-4) and (3-28) are put into (3-48) and (3-49), the following
expressions are obtained as;
73
(3-50)
(3-51)
Then, the matrices in (3-11)-(3-14) can be extracted from (3-50) and (3-51),
respectively as;
(3-52)
(3-53)
74
(3-54)
(3-55)
(3-56)
(3-57)
(3-58)
(3-59)
(3-60)
3.3.5 Steady-State Model
The static model is given in (3-11) and (3-12). Utilizing the matrices above,
the following matrix equations are obtained.
75
(3-61)
(3-62)
(3-62) is just given for the sake of completeness. Required information is extracted
from (3-61) as;
First row:
(3-63)
(3-64)
Second row:
(3-65)
(3-66)
Third row:
(3-67)
(3-68)
Fourth row:
(3-69)
76
(3-70)
Some meaningful results can be inferred from (3-64), (3-66), (3-68) and (3-70).
(3-70) says that mean value of
current,
current,
is zero at
steady-state, this result is consistent and meaningful.
From (3-64) and (3-68), the following equation is obtained. Note that
is equal to
(3-71)
(3-71) gives the relationship between the mean values of the input
voltage, the output voltage and the duty-factor. It is the standard
input/output relationship met in buck-boost structures. Moreover, the
following equation is obtained from (3-64) or (3-68) by substitution of
into the equations. can be extracted from (3-71).
(3-72)
Mean value of the voltage on
can be
replaced with mean value of the input current,
. As inferred earlier,
is
equal to
(3-73)
Multiplication of (3-71) and (3-73) gives input-to-output power equality
so that,
(3-74)
77
Since all electrical elements are assumed to be ideal, the efficiency turns
out to be 100%. Again, this is a consistent result.
3.3.6 Dynamic Model and Transfer Functions
The dynamic model is given in (3-13) and (3-14). Since the aim is to get the
transfer functions, the use of (3-18) will be more proper. Simplification can be
applied before using it, because most of the matrices are constant. The simplified
version of (3-18) is given in (3-75).
(3-75)
will be obtained step by step. Substituting (3-52) and (3-60) into (3-75) gives
(3-76);
(3-76)
where
78
Then,
(3-77)
where
(3-78)
Next step is to obtain
(3-79)
Laplace transformations of and are simply and because initial
conditions are already included in and .
79
(3-80)
(3-81)
where
Using (3-75), (3-79) and (3-81); the following equations are obtained.
(3-82)
where
80
There are six transfer functions to be inferred from (3-82). The outputs are
and
and
. At this
point linearity assumption is necessary. In small signal model, perturbations around
operating point with small magnitudes are considered not to violate the linearity of
the small signal model. By the help of this assumption, superposition theorem can be
applied to (3-82). As a result, the following six transfer functions can be defined.
Duty Factor to Output Voltage Transfer Function:
(3-83)
Input Voltage to Output Voltage Transfer Function:
(3-84)
Output Current to Output Voltage Transfer Function:
(3-85)
81
Duty Factor to Input Current Transfer Function:
(3-86)
Input Voltage to Input Current Transfer Function:
(3-87)
Output Current to Input Current Transfer Function:
(3-88)
3.3.7 Verification of the Transfer Functions
In this section, the calculated transfer functions of coupled-inductor uk
converter with ideal elements through (3-83)-(3-88) are verified by simulations. Two
different simulation tools are benefited for this purpose. One of them is Matlab. In
this simulation tool, the calculated transfer functions are used as a transfer function
block. Small step disturbances in the variables appearing in the denominator of the
transfer functions are applied to the transfer functions and the dynamic responses are
observed in time domain. Second type of simulations uses Simplorer as a simulation
tool. In this simulation group, process circuit model is used instead of the transfer
functions. Only the circuit model is present. The same small step disturbances are
applied to the circuit model and the dynamic responses are observed, which are
82
naturally in time domain. If the results of both simulations turn out to be the same,
one may consider that the transfer function under investigation is valid.
At this point, it is useful to highlight one point. Since the transfer functions
are composed of labels, it is not easy to see that they include operating point
parameters such as
:
In Figure 3-4, step response of the duty-factor-to-output voltage transfer
function
It is plotted in Matlab with the following input parameters.
and
and
are the parameters taken from the simulation which is aimed to be performed at the
fourth end operating point. As it is understood, small deviation in the simulation is
also reflected to the transfer function for a better comparison.
83
Physically, Figure 3-4 is the plot of the output voltage of the circuit operating at
steady-state and disturbed with 0.0025 increase in the duty-factor at t = 0. Now, this
figure will be compared with the simulation result in Simplorer and whether the
derived transfer function is valid or not will be decided. For the time being, it can
only be said that the output voltage oscillates with the peak-to-peak magnitude of
about 0.6V and the frequency of about 75Hz. Since all the elements are ideal and
there is no damping element including the load, the oscillation continues forever.
Figure 3-4: Step response of
in Matlab,
In Figure 3-5, response of
is verified.
84
Figure 3-5: Response of
:
The second output of the circuit is the input current. What the response of the
input current is when the same step change in duty-factor occurs can be plotted by
using the following function in time domain.
Figure 3-6 says that when the duty-factor is increased 0.0025 much at t=0 at
that operating point, the input current shows this variation in time. As expected, the
horizontal axis is in seconds and the vertical axis is in Amperes. To see whether the
derived transfer function is correct or not, the same situation is formed in Simplorer
and the input current is observed in Figure 3-7. Note the perfect similarity between
the figures.
85
Figure 3-6: Step response of
in Matlab,
Figure 3-7: Response of
:
The response of the output voltage when 0.1V step change in input voltage
occurs is plotted in Figure 3-8 by using the following function in time domain.
Figure 3-8: Step response of
in Matlab,
The corresponding waveform in Simplorer is presented in Figure 3-9.
is
validated in this way.
87
Figure 3-9: Response of
Verification of
:
The response of the input current when 0.1V step change in input voltage
occurs is plotted in Figure 3-10 by using the following function in time domain.
The corresponding waveform in Simplorer is presented in Figure 3-11. Similarly,
is verified.
88
Figure 3-10: Step response of
in Matlab,
Figure 3-11: Response of
89
Verification of
:
The response of output voltage when 0.25A step change in output current
happens is plotted in Figure 3-12 by using the following function in time domain.
Figure 3-12: Step response of
in Matlab,
The corresponding waveform in Simplorer is presented in Figure 3-13. In this figure,
the output voltage seems to decrease in time at first sight. However, while its
maximum is decreasing, its minimum is increasing. That is, the envelope is being
thinner as time passes. However, the mid-sinusoidal passing the waveforms does not
change. Most probably, the increase in the thickness of the waveform is the
secondary effect of the step change. The long and the short of it,
is
validated via the similarity between the figures.
90
Figure 3-13: Response of
Verification of
:
The response of input current when 0.25A step change in output current
happens is plotted in Figure 3-14 by using the following function in time domain.
The corresponding waveform in Simplorer is presented in Figure 3-15. In a similar
manner,
is verified.
91
Figure 3-14: Step response of
in Matlab,
Figure 3-15: Response of
92
To sum up, all of 6 transfer functions derived for CIC with ideal elements
have been verified.
3.4 ANALYSIS OF COUPLED-INDUCTOR UK CONVERTER
WITH PARASITIC ELEMENTS
Different from the previous section, parasitic elements are included in the
circuit in this section. Then, the same procedure will be repeated. Explanations given
in the previous section are also valid in this section. Some steps are omitted.
The circuit schematic of CIC with parasitic elements is seen Figure 3-16.
Directions of the currents are also shown.
Figure 3-16 Circuit schematic of CIC with parasitic elements and directions of
currents included
3.4.1 State-Space Equation Set in Mode 1
Equivalent circuit schematic in Mode 1 is shown in Figure 3-17.
93
Figure 3-17 Equivalent circuit schematic of CIC with parasitic elements (directions
of currents included) in Mode 1
Burrowing the derivatives of
and
(3-89)
(3-90)
where
Note that
and
terms have been added to the equations different from the ideal
case. Besides, the coefficients of other terms have changed.
Next step is to obtain derivatives of capacitor voltages. Before proceeding
further, an important point should be considered: Equivalent series resistance, ESRs
of the capacitors are included in the element model. In this situation, capacitor
voltage
represents the voltage of ideal capacitor and not the voltage drop on ESR.
Actually, the same argument applies to other parasitic elements, as well.
For
(3-91)
where
For
(3-92)
where
Then the following equation matrix is obtained. Note that the dimensions of matrix
and vector increase. This modification is needed in Mode 2 and does not affect this
mode practically. Still, the modification is applied here for the consistency in
dimensions of the matrices. The reason of this change will be explained in Mode 2
analysis.
(3-93)
The following matrices are obtained by addition of duty-factor. Since these matrices
are valid in Mode 1, their elements are multiplied by .
95
(3-94)
(3-95)
is written in the following manner. Before that,
(3-96)
(3-97)
(3-98)
where
As a result, the following matrices are obtained as;
(3-99)
(3-100)
3.4.2 State-Space Equation Set in Mode 2
Equivalent circuit in Mode 2 is as shown in Figure 3-18.
96
Figure 3-18: Equivalent circuit schematic of CIC with parasitic elements in Mode
2
Burrowing (2-82) and (2-83) from CHAPTER 2, we have;
(3-101)
(3-102)
where
In Mode 2 operation, diode is conducting and its constant voltage drop
appears in
the equations. Where will it be placed in matrix representation? It is not a state-
variable. Actually, it appears as kind of a voltage sink in the circuit in Mode 2. That
is, it can be regarded as a constant input to the circuit. Consequently, it is going to be
added to vector. Now, other state-variables, namely capacitor voltages, can be
written in the following manner.
97
For
(3-103)
For
(3-104)
The following equation is obtained.
(3-105)
In Mode 2, the duration is . Hence,
(3-106)
(3-107)
is the same with that in Mode 1. Thus,
(3-108)
(3-109)
98
3.4.3 Averaging of Matrices in Mode 1 and Mode 2
Averaged matrices of
(3-110)
(3-111)
(3-112)
(3-113)
As a result, state-space averaged model of the circuit can be represented in
(3-114) and (3-115).
(3-114)
99
(3-115)
3.4.4 Decomposition of Averaged Model into Steady-State and
Dynamic Models
Separation is done in the same way followed in Section 3.3.4.
(3-116)
where
100
(3-117)
Individual matrices can be extracted from (3-116) and (3-117) in the following
manner. The same matrices with the ones in Section 3.3.4 are excluded in order not
to repeat.
(3-118)
(3-119)
(3-120)
(3-121)
101
Note that
(3-122)
(3-123)
3.4.5 Steady-State Model
Utilizing the matrices above, the following matrix equations are obtained.
(3-124)
(3-125)
First row of (3-124):
(3-126)
Using (3-135):
(3-127)
This result will be used later.
Second row of (3-124):
(3-128)
(3-129)
102
It is the same result with the case of ideal electrical elements. It is not affected by the
parasitic elements.
Third row of (3-124):
(3-130)
(3-131)
Fourth row of (3-124):
(3-132)
(3-133)
This result remains the same too.
First row of (3-115):
(3-134)
Using (3-133):
(3-135)
Even though ESR of the output capacitance makes the instantaneous output voltage
differs from the voltage of the output capacitance, their mean values turn out to be
the same, which is an expected result.
Second row (3-115):
(3-136)
This is just a trivial result.
(3-127) and (3-131) may be simplified as in the case of ideal elements.
However, those equations are too complex to be handled successfully. Hence, they
are utilized only for efficiency calculation. By equating both equations, the
relationship between
and
can be derived.
103
(3-137)
where
Then, the efficiency can be calculated as follows;
Using (3-129), (3-133) and (3-137):
and
Further simplification gives the final efficiency expression in (3-138). It will be used
and investigated in the design chapter. However, a clarification may be needed at this
step:
terms need not to be positive. Moreover, they should have signs such that
and
(3-138)
104
3.4.6 Dynamic Model and Transfer Functions
In this section, we will make use of (3-18) again. Simplification due to the
constant matrices can be applied before using it. The simplified version of (3-18) is
given in (3-139), as;
(3-139)
Putting the matrices into their places gives the following equation matrices step by
step. Note that superscript p (probable usage forms:
or
(3-140)
(3-141)
where
(3-142)
and
105
and
106
Using (3-122):
(3-143)
Using (3-119)-(3-121):
(3-144)
(3-145)
where
107
Using (3-143) and (3-145):
(3-146)
where
Using (3-121) and (3-123):
(3-147)
Consequently, is obtained by an addition of (3-147) to (3-146).
(3-148)
108
where
As in the case of ideal elements, again there are six transfer functions to be
inferred from (3-148).
Duty Factor to Output Voltage Transfer Function:
(3-149)
Input Voltage to Output Voltage Transfer Function:
(3-150)
Output Current to Output Voltage Transfer Function:
(3-151)
Duty Factor to Input Current Transfer Function:
(3-152)
Input Voltage to Input Current Transfer Function:
(3-153)
Output Current to Input Current Transfer Function:
(3-154)
It is not possible to express the transfer functions in explicit form due to the
crowd of the terms. However, it can be said that each transfer function is the division
109
of two polynomials in s. While the order of the numerator is at most 3, that of the
denominator is constant at 4.
3.4.7 Verification of the Transfer Functions
The same procedure followed for coupled-inductor uk converter, CIC with
ideal elements will be repeated for CIC with parasitic elements in this section.
Parameters of parasitic elements are taken from the final parameters in the design
section. That is to say, they belong to the implemented circuit. Besides, calculations
and simulations are realized at the same operating point in the ideal case.
Verification of
:
In Figure 3-19, step response of the duty-factor-to-output voltage transfer
function
It is plotted in Matlab with the following input parameters as an addition to
parameters used in the verification of the transfer functions in ideal case.
and
is
the parameter taken from the simulation which is aimed to be performed at the fourth
end operating point.
110
Physically, Figure 3-19 is the plot of the output voltage of the circuit operating at
steady-state and disturbed with 0.0025 increase in the duty-factor at t = 0. By looking
at the figure, it can be argued that the output voltage oscillates at the frequency of
about 72Hz, which is close to the frequency of ideal case. However, because of the
parasitic elements, the oscillation damps in a short time as expected.
Figure 3-19: Step response of
in Matlab,
In Figure 3-20, response of
is verified.
111
Figure 3-20: Response of
:
The second output of the circuit is the input current. What the response of the
input current is when the same step change in duty-factor occurs can be plotted by
using the following function in time domain.
Figure 3-21 says that when the duty-factor is increased 0.0025 much at t=0 at that
operating point, the input current shows this variation in time. As expected, the
horizontal axis is in seconds and the vertical axis is in Amperes. In order to see
whether the derived transfer function is correct or not, the same situation is formed in
Simplorer and the input current is observed in Figure 3-22. Note the similarity
between the figures.
112
Figure 3-21: Step response of
in Matlab,
Figure 3-22: Response of
:
The response of the output voltage when 0.1V step change in input voltage
occurs is plotted in Figure 3-23 by using the following function in time domain.
Figure 3-23: Step response of
in Matlab,
The corresponding waveform in Simplorer is presented in Figure 3-24.
is
validated in this way.
114
Figure 3-24: Response of
in Simplorer, parasitic
elements,
Verification of
:
The response of the input current when 0.1V step change in input voltage
occurs is plotted in Figure 3-25 by using the following function in time domain.
The corresponding waveform in Simplorer is introduced in Figure 3-26. In this
manner,
is also verified.
115
Figure 3-25: Step response of
in Matlab,
Figure 3-26: Response of
in Simplorer, parasitic
elements,
116
Verification of
:
The response of output voltage when 0.25A step change in output current
happens is plotted in Figure 3-27 by using the following function in time domain.
Figure 3-27: Step response of
in Matlab,
The corresponding waveform in Simplorer is given in Figure 3-28. By grounding on
the sufficient similarity between the figures,
is verified.
117
Figure 3-28: Response of
in Simplorer, parasitic
elements,
Verification of
:
The response of input current when 0.25A step change in output current
happens is plotted in Figure 3-29 by using the following function in time domain.
The corresponding waveform in Simplorer is shown in Figure 3-30. Likewise,
is verified.
In conclusion, all of 6 transfer functions derived for CIC with parasitic
elements have been validated.
118
Figure 3-29: Step response of
in Matlab,
Figure 3-30: Response of
in Simplorer, parasitic
elements,
119
CHAPTER 4
DESIGN OF A COUPLED-INDUCTOR CUK
CONVERTER
4.1 INTRODUCTION
This chapter comprises of two main topics. In the first topic, operational
requirements of the converter are defined. Where it will be used, for what purpose it
will be utilized and what are expected from the converter are given in details. Then,
based on these data, technical requirements are determined. Second topic is the
selection or design of the circuit components according to the technical requirements
in the first topic. Then, the design will be verified by simulations or experimental
tests wherever they are needed and possible.
4.2 OPERATIONAL REQUIREMENTS OF THE CONVERTER
Coupled-inductor uk converter implemented in this work will be used
mainly in one of the projects in ASELSAN. Therefore, the design criteria and
parameters are determined by the project specifications. Before determining the
design criteria, possible operating conditions must be considered at this point.
This converter will be used in the vehicle configuration of portable military
equipment or military system. This equipment has two configurations: man-portable
configuration and vehicle configuration. In man-portable configuration, the
120
equipment is powered only by lithium-ion battery pack. According to the energy
capacity of the battery pack, operating time of the equipment is determined. In
vehicle configuration, again lithium-ion battery pack is utilized but it is not the
unique energy source in this configuration. There exists another energy source,
namely the battery group of the vehicle. It is desired to draw as much power as
possible from this source to increase the operating time. Energizing the equipment
only by the battery group of the vehicle, thereby eliminating the need for lithium-ion
battery pack and providing a very long operating time compared to lithium-ion
battery pack case is aimed. At this point, an interface problem arises: Input voltage
level of the equipment does not match that of the vehicle battery. Actually, there are
vehicles having different dc buses at different voltage levels in practice. The
interface relating current levels is more sophisticated than that of voltage levels.
Power is another problematic issue. Therefore, a dc-to-dc converter having voltage,
current and/or power control modes is considered as a necessary element in the
solution of this problem. The details are given in the following paragraphs.
In vehicle configuration, the converter will be fed by a vehicle battery and
supply power to the equipment. Feeding of the converter can be realized in two
ways:
direct connection to vehicle battery,
connection to vehicle battery through vehicle cigarette-lighter adapter
(indirect connection).
This connection difference stems from the current carrying capability of the
interconnecting elements. As it is assumed, in direct connection to the battery,
current carrying capability depends on the input cable of the converter. That is,
current limit is controllable. More current can be drawn by using a thick cable.
However, cigarette lighter adapters have generally low current limits because they
are not designed to supply much power.
If direct battery connection is feasible, the converter is desired to supply its
rated power to the military equipment, or the load. The maximum power
121
consumption of the military equipment is about 400W. In the beginning, the rated
output power of the converter is thought of such that it can supply all the power
needed by the equipment alone. However, after considering the operating conditions
and restrictions, it is understood that this is not possible most of the time. The reason
will be explained in the succeeding paragraphs. For the time being lets go on with
the second scenario. If direct battery connection is not possible, converter will be fed
via vehicle cigarette-lighter adapter. In this scenario, the main power supply of the
equipment is lithium ion (li-ion) batteries, which is normally the case for man-
portable configuration of the equipment. Converter will be used as an auxiliary
power supply. The objective, in doing this, is to extend the operating time with li-ion
batteries by the support of the converter. It is clear that the more current the
converter draws from the cigarette-lighter adapter, the longer the equipment operates.
Converter can even eliminate the need for li-ion batteries if its output power meets
400W. However, vehicle cigarette lighter adapters have strict current limits. It differs
from vehicle to vehicle. One can draw safely a current between 10A to 20A from a
vehicle cigarette lighter adapter. Nevertheless, this current is not enough to feed the
equipment alone, the reason of which will be explained later.
In most of the vehicles, nominal dc bus voltage is 12V, whose source is 12V
lead-acid battery. In military vehicles, however, 2 series connected 12V lead-acid
batteries (totally 24V) are common. Since the equipment is to be used in both
military and civil vehicles, converter is required to be fed at both dc bus voltage
levels. Besides, it is generally accepted that 12V-lead-acid battery voltage varies
between 10V to 14V. When the vehicle is operating, its alternator charges the 12V-
battery group at about 14V. When the vehicle is not operating but the battery is
loaded, its voltage decreases below 12V. According to the load, it may decrease even
below 10V. However, since deep discharging of the batteries decreases the life time
of the battery, decreasing below a limit is generally not allowed. Hence, decreasing
below 10V happens in extreme conditions. Converter is designed to operate under
10V, even down to 5V. Nevertheless, the lower limit for the input voltage at rated
power is accepted as 10V. 10V-14V range for nominal 12V dc bus voltage
122
corresponds to 20V-28V range in nominal 24V dc bus voltage. When both types of
vehicles are considered, input voltage range is determined as 10V-28V.
(4-1)
Input voltage range has been determined. Now, the output voltage range must
be identified. Voltage range of li-ion battery pack gives the necessary information.
Li-ion cells have 3.6V nominal voltage. It is varying normally between 3.0V to 4.2V.
In most of the military applications, nominal 28V is used as a dc voltage level.
Hence, the li-ion battery pack used in the project, which is a standard battery pack in
military applications, is formed by 8 li-ion cells connected in series. By this way, its
nominal voltage is made to correspond to . The
maximum voltage level of the battery pack is and
minimum is . While the equipment (or the system) is
being fed by li-ion battery pack, the voltage in the equipments input may vary
between 24V and 33.6V. This corresponds to the output of the converter. The reason
of this will be explained in the following paragraphs. As a result, the output voltage
of the converter will vary between 24V to 33.6V.
(4-2)
At first sight, power demand of the system can be regarded as low. However,
after some investigation it is seen that it is relatively high for a vehicle, which is
generally not designed to supply such amount of extra power as an addition to its
inherent loads such as electronic control, lighting, air conditioning and hydraulic
systems. In a typical automobile, battery capacity is about 60AH (Ampere-Hour).
This means that if the battery of the vehicle is loaded with 10A, it can give this
amount of power for 6 hours roughly. While the engine of the vehicle is not
operating, the only energy source is battery. Battery cannot withstand probably even
an hour while supplying power only to air conditioning and lighting system.
123
Especially in 12V dc bus voltage, even 60A for 1 hour supply duration- corresponds
to a power of 720W. It means that the power consumption of a vehicle is about this
amount. As seen, energy source is very limited in this mode. Actually, vehicles are
designed to operate while its engine is active. While it is operating, an alternator
mechanically powered by the engine supplies all the electrical power needed by the
vehicle and charges the vehicle battery if power surplus exists at that instant. The
capacity of alternators also changes from vehicle to vehicle but its order is about 70A
for a typical automobile. This value is highly dependent on the rotational speed
(rev/min) of the engine and generally given at the maximum horse power rotational
speed, for example at 5000 rev/min. In other words, a moderate alternator can supply
at most 980W electrical power to the vehicle. This power will decrease to its %80
roughly due to the variation in engine rotational speed. As mentioned above and also
supported by the information gathered from the automobile producers, about 80% of
the alternator power is already consumed by the vehicle itself. The long and the short
of it, 400W-power demand of the equipment cannot be taken from the vehicle battery
-and thereby the alternator- even via direct battery connection. As a result, whereas
current carrying capability is a limiting factor in connection via cigarette lighter
adapter, power handling capability of alternator is a limitation in direct connection to
vehicle battery. At both connection modes, the converter will not meet the power
demand of the equipment alone. Hence, both of the power sources -namely the li-ion
battery pack and the converter- will be connected to the input of the equipment in
parallel. This is why the input voltage of the equipment is the same as the output
voltage of the converter. By the way, another important point should be highlighted.
The system does not always consume 400W. Rather, it has gradual power
consumption ranging from 25W to 400W. Therefore, it is highly possible that
although both of the sources are connected, only the converter can feed the system in
some situations. Note also that li-ion battery pack does not sink current in this mode
owing to its internal hardware, probably a reverse current diode or its equivalent in
terms of function. In order to make the converter the primary power source, control
method should be determined accordingly. In the other mode, again priority should
be given to the converter by the help of the control method. After all these
124
evaluations, the following two operating modes are determined and design is
implemented based on these restrictions in this work.
Direct Connection to Vehicle Battery
In direct connection to the vehicle battery, the input power and the output
voltage will be controlled. Input power control, which also brings instantaneous input
current control, is chosen as 280W in order to supply nearly 250W to the load.
Output voltage limit is determined as 34V. Being slightly greater than 33.6V, it is
selected intentionally. When the input power does not exceed 280W, the converter
can regulate its output voltage to 34V. At this voltage, even though the li-ion battery
pack is in fully charged state (33.6V), it does not supply power to the system. By this
way, it is guaranteed that the primary power source is converter in this mode. When
the load exceeds 250W, the converter delivers its rated power to the system. At the
same time, li-ion battery pack meets the rest of the power demand of the system,
namely at most 150W. In this situation, since the converter is in input power limit, its
output voltage is lower than 34V. Essentially, it is determined by the li-ion battery
pack. As it is discharged, output voltage of the converter decreases. As mentioned
Figure 4-1 Direct connection to vehicle battery,
125
earlier, this voltage range is 24V-33.6V. These two operating possibilities in direct
connection to the vehicle battery case are explained in Figure 4-1 and Figure 4-2.
Investigating the figures from right to left makes them more understandable.
In Figure 4-1, the first possible case of direct connection to the vehicle battery
is shown. In this case, power consumption of the system is lower than 250W. Hence,
the converter meets the power requirement of the system alone. Dashed lines mean
that there is a physical connection but no power flow. That is, li-ion battery pack is
connected to the system but does not supply power to it because of the fact that its
voltage is lower than the output voltage of the converter. Voltage and current level
ranges are also specified in the figure. In determining the current levels, efficiency of
the converter is assumed to be about 90%. Besides, arrows represent the direction of
the power flow.
Figure 4-2 Direct connection to vehicle battery,
In Figure 4-2, the second possible case of direct connection to the vehicle
battery is shown. In this case, power demand of the system is between 250W and
400W. Both of the sources supply power to the system. The converter supplies 250W
of the total power and the rest of the power demand, which is at most 150W, is
126
supplied by li-ion battery pack. Again, the explanations in the previous figure are
also valid here.
Connection to Vehicle Battery via Cigarette Lighter Adapter (Indirect
Connection)
In connection to the vehicle battery via cigarette lighter adapter, input current
and output voltage will be limited. Input current limit is chosen as 10A and output
voltage limit is determined as 34V. The reason behind 34V is the same with that in
the direct connection case. It enables to make the converter the primary power source
of the system. Input current limit (10A) is chosen according to the fuse and current
carrying capabilities of cigarette lighter adapters in the vehicles that are possible to
use for this purpose. As long as the load does not exceed the limit at which the input
current is less than 10A, the converter can regulate its output voltage to 34V. At this
voltage, the li-ion battery pack does not supply power to the system. In this case, one
cannot specify a definite power limit because input voltage range is very wide. At
10A, the input power may be between 100W and 280W. Therefore, in indirect
Figure 4-3 Connection to vehicle battery via cigarette lighter adapter,
and
127
connection, the input current should be taken care. When the load forces the input
current to exceed 10A in order to regulate its output voltage at 34V, the converter
limits its input current to 10A and ceases to regulate its output voltage. As a result,
the output voltage decreases and li-ion battery pack starts to supply power.
Depending on the state of the battery pack charge and the operating conditions, the
output voltage of the converter can be between 24V and 33.6V. In this case,
of the system power demand is supplied by the converter, and the rest of the
demand is supplied by the li-ion battery pack. These two operating possibilities are
explained in Figure 4-3 and Figure 4-4. Again, investigating the figures from left to
right eases to understand them.
In Figure 4-3, the first possible case of indirect connection to the vehicle
battery is shown. In this case, the power consumption of the system is lower than at
most 250W, which happens at the input voltage of 28V. This limit is 90W at the
input voltage of 10V. Whatever the input voltage is, it is assumed that the converter
meets the power requirement of the system in this case alone. Li-ion battery pack is
connected to the system but does not supply power to it. Voltage and current level
Figure 4-4 Connection to vehicle battery via cigarette lighter adapter,
and
128
ranges are also specified in the figure by the efficiency assumption of about 90%.
In Figure 4-4, the second possible case of indirect connection to the vehicle
battery is shown. In this case, power demand of the system is between 90W and
400W. Both of the sources supply power to the system. The converter supplies at
least 90W (at
), is
supplied by li-ion battery pack.
Possible operational conditions are defined so far. Based on them, the
following results are obtained in terms of the converter properties.
Input voltage range of the converter is 10-28V, and output voltage range is
24V-34V. It means that the converter is required not only to step-up but also
to step-down the input voltage. uk dc-to-dc converter, which is one of the
buck-boost converter types, is chosen for this purpose.
The converter should be designed to supply 250W maximum within the input
and output voltage ranges mentioned in (4-1) and (4-2).
Especially in controlling the input current, instead of its instantaneous or peak
value, its mean value will be controlled. However, in basic uk converter,
input current has such a waveform that it has a dc level and a triangular ac
component (ripple) superimposed on it. An example to this kind of a
waveform is shown in Figure 4-5. In this example, peak-to-peak ripple of
Figure 4-5 Input current waveform of a basic uk converter with 40% peak-to-
peak ripple
129
40% of the mean input current is selected typically just for illustration
purpose. As understood from Figure 4-5 easily, the instantaneous value of the
input current exceeds the current limit value in reality. In this situation, if the
current protection mechanism of the cigarette lighter adapter is a fast acting
one, it may be triggered by this type of current waveform. Even if the
mechanism is a slow acting one, for example a fuse, this type of a waveform
may have an adverse effect on it. Fuses normally have quite small internal
resistances and they are heated up and become open by the heat caused by the
RMS value of the current passing through them. As the ac component of the
current waveform increases while its dc level is kept constant, its RMS value
increases as well. For example, the waveform seen in Figure 4-5 has a mean
value of 10A but RMS value greater than 10A. In this situation, fuse may be
triggered as well due to the unwanted feature of the current waveform. In
buck-boost converter topology, this situation is worse because it has pulsating
input current waveform. Hence, uk converter seems as a better solution.
Moreover, coupled-inductor uk converter with ripple-free input current
version seems to be the best solution. As a result, CIC with ripple-free input
current version is decided to design in this work.
The converter is required to conduct different control strategies in direct and
indirect connections to the vehicle battery. They are mentioned in the
succeeding pages. Therefore, there must be a way for the converter to
understand which control strategy to implement at that operation. When it is
investigated, it is seen that this cannot be done automatically, for example by
sampling the input and output voltages. Only an operator can determine the
control strategy by looking at the connection type, direct or indirect. Hence, a
manual switch will be added to the converter. It will have two positions.
According to the position of it, the converter will run the corresponding
control strategy.
130
As mentioned earlier, in direct connection to the vehicle battery, output
voltage and input power will be controlled. In order to do this, feedbacks will
be taken from output voltage and input current. Most of the time, the system
is expected to consume its rated power of 400W; therefore, the converter will
be performing input power control. It will rarely control the output voltage.
Hence, this control mode is called constant input power mode in this work.
In indirect connection to the vehicle battery, output voltage and input current
will be controlled. In this control mode, again, input current and output
voltage will be taken as feedbacks. Similar to the direct connection case, the
converter will control the input current most of the time. Therefore, this
control mode is called constant input current mode.
At this point, selection and design of the circuit components can be
performed.
4.3 SELECTION AND DESIGN OF THE CIRCUIT
COMPONENTS
In this section; selection of the switch, diode, energy transferring capacitor
(
10V and
24V
Second End Operating Point:
10V and
34V
Third End Operating Point:
28V and
24V
Fourth End Operating Point:
28V and
34V
The efficiency of the coupled-inductor uk converter, CIM is assumed to be
90% in the preparation of Table 4-1. One may refer to the circuit schematic of CIM
with parasitic elements in Figure 3-16 for the labels in the tables of this section. The
instantaneous values, with small letters, are shown in that figure. In the tables,
however, steady-state (or mean) values, with capital letters, are used; because it
makes more sense. Similarly, mean value of the duty-factor, is used in the tables.
and
labels are new. When the switch or the diode is on, it carries the
sum of
and
, which is
the sum of input and output voltages. This value is labeled by
is the RMS
current of the switch. It will be used in power loss calculation.
132
Table 4-1 Voltage and current stresses on the switch
10 27,8 24 10,4 0,706 0,294 34 38,2 32,1
10 27,8 34 7,4 0,773 0,227 44 35,1 30,9
28 9,9 24 10,4 0,462 0,538 52 20,3 13,8
28 9,9 34 7,4 0,548 0,452 62 17,3 12,8
As inferred from Table 4-1, the selected switch has to withstand 62V during
will be used in conduction loss calculation for the switch. This value is
simply obtained by root mean square definition.
(4-3)
where
time interval
only. Hence, when this simplification is applied to (4-3), it yields;
(4-4)
Hence, the maximum RMS current which can be carried by the switch is given as;
133
Then, the conduction loss in the switch can be calculated by;
(4-5)
As it is seen, even 10m
Another important point in determining the suitable switch is the switching
loss. The switching loss on the switch can be expressed as
[
2
]
;
(4-6)
Hence, in order to keep
low, rise and fall times should be low for the switch
selected for use. Consequently, MOSFET with the lowest total power dissipation
among the available MOSFETs is chosen as the following one;
Company: Infineon Part Number: IPP045N10N3 G
Continuous
Continuous
Note that:
and
Thus, the conduction and the switching power losses at end operating points
for the switch are calculated and tabulated as Table 4-2.
134
Table 4-2 Power dissipation on the switch at end operating points
10 27,8 24 10,4 4,3 4,7 9,1
10 27,8 34 7,4 4,0 5,6 9,6
28 9,9 24 10,4 0,8 3,9 4,7
28 9,9 34 7,4 0,7 3,9 4,6
Table 4-2 shows that the switching losses are large. This implies that increasing the
switching frequency is not reasonable due to the increasing switching loss.
4.3.2 Selection of Diode
As a diode type, power Schottky is chosen. The reason for this is its low on-
state voltage drop and low reverse recovery time. Hence, they are generally preferred
in high frequency switched mode power supplies.
The current and voltage stresses on it are given in Table 4-3.
Table 4-3 Voltage and current stresses on the diode
10 27,8 24 10,4 0,706 0,294 34 38,2 11,2 20,7
10 27,8 34 7,4 0,773 0,227 44 35,1 8,0 16,7
28 9,9 24 10,4 0,462 0,538 52 20,3 11,0 14,9
28 9,9 34 7,4 0,548 0,452 62 17,3 7,8 11,6
Diode also carries the current
during
during
time interval.
(4-7)
and
values are given for diode power loss calculations and will be used in the
succeeding paragraphs. Before loss calculation, diode should be chosen according to
135
voltage and current stresses. As understood from Table 4-3, selected diode has to
withstand 38.2A and 62V. As a consequence, the following diode is chosen.
Company: ST Microelectronics Part Number: STPS80170C
Note that:
Actually, using 38.2A as a mean value is not exactly true. However, in many of the
diode datasheets, there is no detailed information about allowable peak repetitive
forward current. Therefore, 38.2A is used as the mean value to be on the safe side.
In the datasheet of this diode, the following equation is given for use in the
conduction loss calculation as;
(4-8)
From which, the following relationship is established;
(4-9)
It is clearly seen that
and
10 27,8 24 10,4 8,3
10 27,8 34 7,4 5,8
28 9,9 24 10,4 7,5
28 9,9 34 7,4 5,2
To sum up, power dissipation on the diode seems to be reasonable.
4.3.3 Selection of
holds
is charged by
during
and
discharged by
during
is severe. In
order to see the effect of current stress, simple capacitor model should be
investigated.
Figure 4-6 Simple capacitor model
As seen in Figure 4-6 and Figure 3-16, a capacitor has an equivalent series
resistance, ESR. This resistance has two important effects. First, the RMS current
passing through the capacitor causes power dissipation on the capacitor, which
results in heating of itself. The amount of power dissipation can be calculated in the
137
following and known manner. Note that although the mean value of the capacitor
current is zero at steady-state, its RMS value can be very large.
(4-10)
The second important effect of ESR is related to voltage ripple on the capacitor. As
mentioned earlier, capacitors have zero mean but nonzero RMS current values at
steady-state operation. It means that
(4-11)
Note that
is used instead of
. Contrary to
changes.
(4-12)
In order to clarify the situation, lets consider that the circuit is transiting form
Mode 2 to Mode 1. In Mode 2,
. At this
moment, (4-13) can be written, in which superscript
(4-13)
In Mode 1,
is
is equal to
(4-14)
Just before the switching, (4-13) is valid, but after the switching, (4-14) is valid. As a
known fact,
is equal to
.
However,
(4-15)
(4-16)
(4-17)
This sharp voltage change occurs at each switching instant. Whereas
decreases
.
When the available capacitors are searched, it is decided to use the following.
Company: Cubisic Part Number: A710160
Type: Aluminium Electrolytic
139
Aluminium electrolytic type capacitor is selected because it has higher energy
density among other capacitor types. In other words, it presents more capacitance for
unit volume.
Note that the maximum capacitor voltage is greater than the maximum
.
Table 4-5 Voltage ripple consideration of
10 27,8 24 10,4 0,649 0,025 0,674 0,680
10 27,8 34 7,4 0,597 0,019 0,616 0,880
28 9,9 24 10,4 0,346 0,016 0,362 1,040
28 9,9 34 7,4 0,294 0,013 0,307 1,240
The term
, which is
due to its ESR. In the analysis above, this voltage change is labeled by
.
From now on, it is labeled by
is replaced by
(4-18)
(4-19)
(4-20)
As it can be understood from the derivation above,
and
on the table. is 2% of
the voltage on the capacitor,
. As seen in Figure 4-7, total ripple magnitude (671mV) is very close to the
calculated one (674mV). There is a sharp change in the capacitor voltage, whose
magnitude is as expected. Lines with lower slopes at the bottom and top are where
the capacitor is being either charged or discharged. There is a small amplitude
difference between them. It is due to the fact that
RMS current and power dissipation with the selected capacitor at each end
operating points are given in Table 4-6.
Table 4-6 Power dissipation on
10 27,8 24 10,4 0,706 0,294 17,4 5,2 1,7
10 27,8 34 7,4 0,773 0,227 14,7 3,7 1,2
28 9,9 24 10,4 0,462 0,538 10,2 1,8 0,6
28 9,9 34 7,4 0,548 0,452 8,6 1,3 0,4
(4-21)
Note that the maximum RMS ripple current is less than the maximum
allowable RMS current of the capacitor bank.
142
each of
them. As it is seen in Table 4-6,
resembles to that of
(4-22)
The difference current between
and
is supplied by
, a)
waveform, b)
waveform, c)
waveform
143
dc current. In this situation, the waveforms shown in Figure 4-8 are observed. As
seen clearly,
has not severe current stress. It meets the current unbalance between
the output current and the output inductor current. Hence, it may be called output
filter capacitor. In order to determine the voltage ripple on
. Hence,
(4-23) can be written easily. Contrary to the ripple behavior in
(4-23)
The term
and
it is caused by charging and discharging of the capacitor. In order to find a
relationship between
is charged during
(4-24)
(4-25)
(4-26)
(4-27)
(4-28)
Similar derivation of (4-28) can also be found in [2]. This derivation is given
intentionally in order to show the difference between
and
in terms of voltage
ripple consideration.
(4-23) and (4-28) are ready to use in determining the output capacitance.
Before it, a simple modification may result in a better solution. 2% voltage ripple
restriction for
voltage is
inside the circuit contrary to
It is a surface mount ceramic capacitor. Since the current stress is very light and the
transferred energy is very small, the required capacitance turns out to be low. This
justifies using ceramic capacitor. Its ESR and
Using these values, Table 4-7 is formed. Note that
10 27,8 24 10,4 0,68 0,005 0,090 0,095 0,240
10 27,8 34 7,4 0,74 0,006 0,098 0,104 0,340
28 9,9 24 10,4 1,24 0,009 0,165 0,174 0,240
28 9,9 34 7,4 1,47 0,011 0,196 0,207 0,340
As it is seen in Table 4-7, over 94% of the ripple is caused by charging/discharging
process of the capacitor. The reverse is valid in
and
do not
totally coincide in time. In other words, their maximum points occur at different
times, although they are close to each other. Assuming the coincidence of the peak
values in time gives the maximum ripple voltage and this case is called as the worst
case. To be on the safe side and for simplicity, the worst case is assumed and the
design is implemented according to this assumption.
Now, the ripple voltage calculations will be verified by simulation. It is done
with all parasitic components at the same end operating point: Full load,
Calculated value is explained in Figure 2-22 and commented.
Now, lets look at the power loss on
(4-29)
By utilizing this equality, power dissipation on the output capacitor is presented in
Table 4-8. As seen in this table, total power loss on the output capacitor bank is in
1mW order. Hence, it is neglected.
147
Table 4-8 Power dissipation on
10 27,8 24 10,4 0,68 0,20 0,0003
10 27,8 34 7,4 0,74 0,21 0,0003
28 9,9 24 10,4 1,24 0,36 0,0010
28 9,9 34 7,4 1,47 0,42 0,0014
4.3.5 Design of Coupled-Inductor
A coupled-inductor providing ripple-free input current is aimed to design in
this section. First of all, its requirements such as the required inductances and the
maximum current will be determined. Then, core material selection topic will be
handled. Coupling requirement will be tried to meet. After theoretical work,
practically implemented coupled-inductor will be investigated and tested. Its
parameters will be measured and compared with the theoretical ones.
Lets start with determining the inductances of the inductors. The coupled-
inductor is aimed to provide ripple-free input current. As discussed in CHAPTER 2,
will determine
current ripple. Then, with respect to the parameters determining
must have
specific parameters in order to provide ripple-free input current waveform. In the
same chapter, it is also mentioned that
.
Further details can be found there. In short, required
? It is known that a
larger inductance results in a large volume and weight. As an advantage, core loss
and part of the copper loss caused by ac current diminish because the magnitudes of
current and flux ripples decrease. In that situation, since the magnitudes of current
ripple and flux ripple are smaller, complex calculations such as skin effect, proximity
effect and hysteresis loss can be neglected. The inverse of all these comments can be
argued for a smaller inductance. It is generally desired to have small components.
Hence, a small inductance must be preferred. What must be the lower limit for the
inductance? It is mentioned earlier that the minimum load of the system will be 25W,
148
which corresponds to the 10% of the full load of converter. It means that the current
at the output will be at least 10% of the full load. It is known that
is equal to
. In
that situation, current waveforms similar to the ones shown in Figure 4-11 appear.
Operating conditions are not given intentionally but note that the simulation is
performed with the parasitic elements included. Assume that
or
decreases
Figure 4-11
and
below zero and change the direction. This point is chosen as the limit. In other
words,
will be determined such that even at the lightest load (i.e. 10% of the full
load) at all possible operating points,
and other circuit parameters. Lets start with the known inductance
relationship.
(4-30)
Since the applied voltages are assumed to be constant, (4-30) can be simplified into
(4-31).
(4-31)
is either
equals
and equals
is 2 times 10% of the full load current. However, the full load current
differs with respect to the output voltage (24V-34V) for that operating point.
Meanwhile, duty-factor changes with the input voltage. Hence, the best way of
determining
is to calculate the required inductances at 10% load for the four end
operating points and pick the maximum of them.
Table 4-9 Required
10 27,8 24 10,4 0,706 1,04 34
10 27,8 34 7,4 0,773 0,74 53
28 9,9 24 10,4 0,462 1,04 62
28 9,9 34 7,4 0,548 0,74 104
In Table 4-9, superscript 100% means the full load values and superscript
10% means 10% of the full load values. For example, for the first end operating
point, mean value of the output current at %10 of the full load is 1.04A. It means that
is 2.08A. Also,
The same calculation is done for each end operating points. As seen in Table 4-9, the
worst case is the fourth end operating point. At this case, the required inductance
turns out to be 104uH. Hence,
(4-32)
The first specification of
is
determined to be greater than 104uH, the ripple magnitudes and the maximum
current diminish further.
Table 4-10 The maximums and the minimums of
10 27,8 24 10,42 0,706 10,75 10,08 6,49
10 27,8 34 7,35 0,773 7,72 6,98 10,06
28 9,9 24 10,42 0,462 11,04 9,80 11,88
28 9,9 34 7,35 0,548 8,09 6,62 20,00
As seen in Table 4-10, the maximum current that is carried by
instantaneously is 11.04A. Also note that the peak-to-peak ripple percent of
turn
out to be 20% at the fourth end operating point, which is not a surprise. Since the
inductance limit value criterion is based on 10% and the maximum inductance is
obtained at the fourth end operating point, superposing peak-to-peak 20% ac current
waveform on 100% dc level at full load naturally gives 20% ripple at the same end
operating point.
151
Required information in order to design
is constant and
. Besides, as
explained in CHAPTER 2, fluxes created by
and
and
occurs at the
first end operating point. Hence, these values are written in Table 4-11. Required
value (104.4uH) is another input. Diameter of the selected wire is 0.55mm. Instead of
using single but thick wire, multiple of thin wire is preferred. Thin wire is generally
used where ac current amplitude and switching frequency are high. In this way, skin
effect is considered and thin wire is used in order to decrease the ac resistance of the
wire. In this problem,
inductor were
designed high, its loss would be low because ac current magnitude is very low. Still,
skin depth may be checked. Its known formulation is given in (4-33).
(4-33)
Skin depth is calculated as:
Skin depth is defined as the distance from the surface to the inner of the conductor
where current density is
. Then, current
density limit must be set in order calculate how many wire must be paralleled in
order to carry that much of input and output currents. When many examples are
investigated, it is seen that most of them are in the range of 2A/mm
2
to 6A/mm
2
.
Midpoint of them, 4A/mm
2
, is selected and checked whether it is reasonable or not.
In this situation, it is observed that windings have moderate series resistances and
power dissipations, which will be presented in the following paragraphs. Of course,
selecting a current density below 4A/mm
2
would be better, but in 4A/mm
2
case
windings just about fit into the coil former, to which the windings are wound. Hence,
4A/mm
2
is evaluated as a satisfactory value. Using this information, required total
cross section area for
is calculated as:
How many wires can provide this total cross section area is found as the ceiling value
of
and 12 wires for
Number of turns is required now. This value comes from Table 4-14. As mentioned
earlier, the design procedure is iterative and the tables are correlated. Determination
of number of turns will be explained later. Lets take it for granted for the time being.
Then, the net occupied area by 19.5 turns of 30 wires is calculated as:
154
In Table 4-11, gross values are given. Net area is named as the area of circle,
whereas gross area is considered as the area of the square which includes that circle.
There is a correction factor of 4/ between the areas. As it is known, there exist areas
between the circles which are not utilized. In order to asses more accurately, this
correction is performed at this stage. As a result, gross occupied area of
winding is
calculated as:
Similarly, gross occupied area of
maximum 10,75
104,4
0,55
0,24
6,94
30
2,76
12
177
66
243
462
Lets go on with Table 4-12. In this table, some parameters of the selected
core are given. Before it, the identification of the core is presented as;
155
Company: Epcos Part Number: PM 74/59
Material: Ferrite N27 Air gap: 3.8mm
First of all, why gapped ferrite core is chosen can be explained briefly. As it is
known, there is a trade-off between high
and
. As a
known fact, inductance of an inductor decreases little or much as the current
increases. In other words, reluctance of an inductor increases with increasing flux
density. Characteristic of this change differs from material to material. This fact is
observed in B-H curves. For example, air has a linear B-H curve, its slope is very
low and magnetic saturation is not observed. B-H curves of ungapped ferrite
materials are comprised of two almost linear portions. First portion represents the
magnetic permeability of ferrite material and its slope is very high. Once
value
is reached, the second linear portion is observed. Because of the saturation, its slope
and hence magnetic permeability is low in this portion. This type of saturation is
called hard saturation or sharp saturation. B-H curves of gapped ferrite materials
resemble to that of ungapped ferrite materials. The only difference is the slope of the
first linear part. In gapped case,
. Due to their
non-homogeneous structure, they show different reluctance values at each operating
point. It means that powdered cores exhibit different inductance values according to
the current. Detailed explanation of this issue can be found in [29].
Based on the previous paragraph, it can be said that the leakage inductances
of the coupled-inductor (namely
and
) are independent of
and
; because
their reluctance paths are air. However,
and
would be satisfied probably at one operating point only. However, it is desired to
have ripple-free input current at all possible operating points. Hence, using gapped
ferrite material as a core and staying always at the first linear portion of it and
avoiding magnetic saturation is regarded as the best solution.
Table 4-12 Necessary core data
Core Data
315
Assumed maximum
ratio 0,05
ratio 1,05
630
462
Table 4-12 can be investigated now. In the datasheet of the selected core,
instead of the equivalent reluctance of the core and the air gap,
value is given. In
most of the core datasheets, this approach is adopted.
(4-34)
Therefore, the reluctance value of the core, namely
(4-35)
(4-36)
Using the equation above,
and
. Since
is known,
. However,
ratio should be
determined according to
) ratio or simply
and
ratio. Of
course, measuring
and
windings are
wound is possible but it does not give an accurate result. Because of the practical
differences between the groups of windings of the same inductor, their individual
inductances differ. It results in relatively high, for example 10%, difference between
the measurements taken after the first groups of wires of
and
windings are
wound and the measurements taken after all groups of wires of
and
windings
are wound. This high difference may easily disturb the balance established after the
first groups of wires of
and
ratio is assumed. In
this work, this value is assumed as 5%. According to this value, required
ratio according to (2-53) is obtained in the following manner.
If
is taken as 100,
If the measured
ratio
is planned to make 1.05 by adding a series connected adjustment inductor to
inductor as a leakage inductor. In other words, if the leakage inductance of
winding turns out to be smaller than expected, the above equation is satisfied by
increasing
winding.
Hence, there is no need to modify the coupled-inductor, which is not really practical.
However, if
ratio has
to be increased, which means rewinding of the coupled-inductor. Because of this
reason, the maximum
is calculated by using
and
. Once
is obtained, required
(4-37)
As a next step,
is calculated by using
and
values
are 18.7 and 17.8, which are not very practical. Therefore, Table 4-13 is labelled as
theoretical.
Table 4-13 Theoretical outputs with respect to core data
Theoretical Outputs with Respect to Core Data
3174603
99,4
17,8
18,7
709
223
0,35
As it is known, magnetomotive force (MMF) can be expressed as:
Since MMFs of the coupled-inductors are additive, the maximum total MMF can be
defined in the following manner.
and
In order to find the flux, the following equality can be utilized.
(4-38)
160
Hence,
This is the maximum flux present in the core. The maximum flux density can be
calculated by using the following equality.
(4-39)
When the datasheet of core material is investigated, it is seen that 0.35T remains on
the first linear portion of the B-H curve. That is, magnetic saturation is avoided and
almost constant magnetic permeability is expected. The snapshot of the B-H curve of
the used core is given in Figure 4-13.
Figure 4-13 B-H curve of the selected core material
[30]
At this point, all the parameters seem to be reasonable and iteration is
finished. Actually, two conditions are checked. First one is whether the windings fit
161
into the coil former or not. Second one is whether the maximum flux density is lower
than saturation flux density of the core or not.
As mentioned earlier,
and
18,25
19,5
1,068
104,9
738
232
0,37
For illustration purpose, photographs of the designed and implemented
coupled-inductor are presented in Figure 4-14 and Figure 4-15.
Figure 4-14 Designed and implemented coupled-inductor, photograph 1
162
Figure 4-15 Designed and implemented coupled-inductor, photograph 2
After the coupled-inductor is implemented, some measurements are taken in
order to determine magnetizing and leakage inductances. In order to ease the
visualization, the equivalent coupled-inductor model is given in Figure 4-16.
Figure 4-16 Equivalent coupled-inductor model
Similar to transformer, open circuit and short circuit test can be applied in order to
find the model parameters of the coupled-inductor. They are given as;
(4-40)
163
(4-41)
(4-42)
(4-43)
(4-40)-(4-41) or (4-42)-(4-43) can used in order to determine the parameters:
and
and
(4-44)
When the measurements belonging to the first winding side are used, the following
parameters are found.
As mentioned in CHAPTER 2,
must be equal to
From the calculations above, it is seen that there is a difference between
and
.
As mentioned at the beginning, this difference is expected. As a solution, an
adjustment inductor is used as an addition to the leakage inductance
Adding a discrete inductor with the value of 6.2uH in series with positive terminal of
winding satisfies the condition, thereby providing ripple-free input current at all
operating points. As a consequence, the following model and parameters are valid
from now on.
164
Figure 4-17 Adjusted coupled-inductor model and its parameters
Whether the adjustment inductor approach is correct or not will be checked
by simulation. For this purpose
and
As an addition to the above values, the following values are used in the simulation.
Simulation is performed with the ideal elements in order not to mix the reverse effect
of parasitic elements with the possible failure of the approach. It is performed at
and
(4-45)
As used earlier,
As given in Table 4-11, there are 30 parallel wires.
167
Parasitic resistance of
As presented in Table 4-11, there exist 12 parallel wires in
winding.
Parasitic resistance of the adjustment inductor should be considered as well.
As mentioned earlier, 6.2uH is required. In order to obtain this inductance value, 3 of
2uH inductor is decided to use. Note that rather than 6.2uH, 6uH will be used due to
the availability of the inductors. This much inconsistency may also stem from the
measurement error of the coupled-inductor parameters or the parasitic elements
easily. Therefore, trying to find exactly 6.2uH is regarded as a useless effort.
Necessary data of the selected adjustment inductor is given as;
Company: Coilcraft Part Number: SER2010-202MLB
Core Type: Ferrite
Since 3 of it will be used in series, total resistance is 3m.
As a result, equivalent coupled-inductor model with parasitic resistances can be
formed. It is seen in Figure 4-21.
168
Figure 4-21 Adjusted coupled-inductor model and its parameters with parasitic
resistances
Power dissipations of the inductors at each end operating point should be
considered at this stage. In Table 4-15, required data is given.
Table 4-15 Power dissipation on the coupled-inductor and the adjustment inductor
at end operating points
10 27,8 24 10,42 0,68 10,42 4,8 1,6 0,3 6,7
10 27,8 34 7,35 0,74 7,36 4,8 0,8 0,2 5,7
28 9,9 24 10,42 1,24 10,42 0,6 1,6 0,3 2,5
28 9,9 34 7,35 1,47 7,37 0,6 0,8 0,2 1,6
As it is known, power dissipation on a resistor depends on the RMS value of current
passing through it. Hence, RMS values of
and
is constant, its RMS value is the same with its mean value.
(4-46)
169
As shown in Figure 4-8,
(4-47)
Its dc level is equal to
. Peak-to-peak magnitude of
is labeled as
. ac
triangular waveform is not necessarily equal i.e. duty-factor is generally different
than 0.5. By using the RMS value definition in (4-3),
is expected to be very
close to its mean value. Still, the analysis seems to be necessary as a check.
(4-48)
(4-49)
Note that the integral of
is zero.
(4-50)
By using (4-29), the following data can be inferred.
(4-51)
As a consequence,
(4-52)
As seen in Table 4-15,
as expected.
Once the RMS values of currents are obtained, power dissipations on the
inductors can be calculated. As observed in Table 4-15, most of the power is
170
dissipated on
and
.
Figure 4-22 Detailed
172
The maximum flux density is also close enough to the calculated value (0.37T) in
Table 4-14. Another issue is the flux swing. Peak-to-peak flux swing percent is about
1.6%. Since it is very low, core losses are not calculated and neglected in this work.
Figure 4-24 Mutual flux
side. Voltage
is applied not only to
and the
adjustment inductor. For simplicity, this voltage will be labeled by
from now
on.
as
and
plus
adjustment inductor combination as
and
Side
Measured Voltage
on
Winding
Voltage
Ratio
0 9,36 10,20 1,09
2 9,28 9,84 1,06
4 9,44 9,68 1,03
6 9,36 9,36 1,00
As the adjustment inductance value comes close to 6.2uH, the induced
voltage on
approximates to
.
If the same voltage is applied to
and
(or
(or
) are controlled in a
cascaded manner. Since the input current can change faster than the output voltage, it
is placed the inner loop. Output voltage loop is the outer control loop and gives a
reference value for the input current according to the output voltage. Although
current mode control is adopted in a reasonable manner, voltage mode control has
also been examined. Before a more complicated current mode control, simpler
177
voltage mode control has been desired to implement as an intermediate step. Voltage
mode control simulations performed in Simplorer with many different controller
parameters did not give a stable output voltage. There existed a small amplitude
oscillation around the reference value. Oscillation frequency was the same with the
frequency of the nearest pole of the duty-factor-to-output voltage transfer function,
to the origin in s-domain. Then, voltage mode control was applied on the
implemented circuit experimentally and the same oscillation was observed. In that
situation, a large oscillation in the input current was realized. With this high
amplitude oscillation in input current, output voltage oscillation seemed to be
inevitable. Since the ultimate aim is current mode control, further elaboration in
voltage mode control was not performed and directly jumped to current mode
control. Since current mode control has been performed easily and consistent with
the calculations, voltage mode control of this circuit has then been evaluated as not
easily applicable.
In the inner and outer control loops, PI controller is preferred. It is known that
integral term lets the steady-state error be zero, which is a desired goal. Besides, in
most of the industrial applications, PI controller is seen satisfactory. Therefore, PI
controller is implemented as a first choice and gives a satisfactory result. Hence,
other controller types are not tried.
Up to this point, the controller scheme is drawn roughly and it will be
detailed in the following sections. According to the control block diagram and the
related transfer functions found earlier, necessary controller functions will be
obtained in s-domain. Then, using a bilinear transformation, controller functions will
be converted into their z-domain correspondence. By inverse z-transformation,
controller functions in discrete time domain will be obtained. Then, those functions
will be embedded to dsPIC.
178
5.2 CONTROL BLOCK DIAGRAM
Based on the discussions above, the control block diagram of the circuit can
be formed in the following manner.
Figure 5-1 Control block diagram of the circuit
The difference between the reference output voltage,
. Similarly,
the difference between
gives an input
current error,
and
produces a duty-factor, . This duty-factor is applied to the converter and the
converter produces
and
. Afterwards,
and
are taken as feedbacks and this control cycle is repeated continuously. Note
that the faster variable,
(5-1)
(5-2)
Figure 5-2 Mathematically equivalent control block diagram of the circuit
At this point, transfer functions of the control system, namely
and
or
or
terms. Therefore,
will be
determined first and then
and the
other transfer functions.
180
5.3.1 Controller Design of the Current Loop
Control block diagram of the inner loop, or current loop, is given solely in
Figure 5-3.
Figure 5-3 Control block diagram of the current loop
The closed loop transfer function of the current loop,
(5-3)
(5-4)
(5-5)
Denominator of a closed loop transfer function is defined as the characteristic
equation. Characteristic equation of the closed loop transfer function of the current
loop is labelled as
(5-6)
By definition, characteristic equation is the summation of 1 (or more accurately
) and open loop transfer function. Therefore, open loop transfer function of the
current loop,
(5-7)
181
At this point,
, or the other transfer functions derived before, has the following terms
parametrically:
and
. Note that other than the circuit parameters, also operating point parameters
take place in the transfer functions. In other words, the converters transfer function
changes as the operating point changes. As expected, transfer functions at different
operating points are calculated and it is seen that there are only minor differences
between the pole and zero locations of the transfer functions. Therefore, an operating
point can be chosen and the transfer functions can be calculated at that point. Here,
the first end operating point (
(5-8)
In order to get acquainted with
easily. Hence, zoomed-in snapshot of the red circle is given at the right figure. As it
is seen clearly,
has four poles and three zeros. Since all poles have negative
real parts,
is a minimum phase
system, bode plot approach can be applicable. In this approach, the aim is to get
positive and proper phase and gain margins in the bode plot of
. Since
.
As a starting point, the bode plot of
This plot will be corrected to a desired shape by the help of
. As mentioned
earlier, PI controller is proposed as a controller. Its transfer function can be written in
the following manner.
(5-9)
where
and
and a gain of
rad/sec.
184
Figure 5-6 Bode plot of
By adding a real zero at -1000 in s-plane, the following bode plot is obtained. This
zero placement brings the following restriction to
(5-10)
Note that the phase is almost always greater than -180 expect a frequency band
about
rad/sec in Figure 5-7. At this step, only the gain variable is left to adjust
the phase and gain margins. Figure 5-7 is actually drawn such that
. Should
When it is decreased, magnitude plot moves downward. The magnitude plot will
cross 0dB at a lower frequency and phase margin will be positive. At the same time,
gain margin increases as well. As a result, gain must be decreased.
Using the sisotool function of Matlab, this gain iteration can be performed
in an online manner. At the gain of 0.0035, reasonable gain and phase margins are
obtained. It is seen in Figure 5-8. Theoretically, positive phase and gain margins are
sufficient in order to have a stable loop. In practical, it is generally accepted that the
phase margin must be larger than 45 because of the practical reasons. It is
recommended that it should be between 45-60. Actually, this phase margin defines
a critically damped step response. While larger phase margin results in an
overdamped response, smaller phase margin leads to an underdamped response. That
is, increasing phase margin much makes the converter slower. Decreasing phase
margin much risks the stability of the converter. Hence, as seen in Figure 5-8, phase
186
Figure 5-8 Bode plot of open loop transfer function,
margin of 54.6 is obtained. There is no strict requirement on gain margin except its
sign. 6-12 dB gain margin is recommended. Gain margin of 25dB is left as it is. Gain
of 0.0035 brings another restriction to the controller transfer function.
(5-11)
Using (5-10), integral gain can be calculated as:
(5-12)
As a result, transfer function of the current controller is determined in s-domain.
(5-13)
187
It will be converted to its z-domain and discrete time domain equivalents in the
succeeding sections. Before it, transfer function of the voltage controller will be
determined.
Finally, bode plot of the closed loop transfer function is given in Figure 5-9
as an information about frequency response of the current loop.
Figure 5-9 Bode plot of
5.3.2 Controller Design of the Voltage Loop
Simplified version of control block diagram of the outer loop, or voltage
loop, is shown in Figure 5-3. Note that the current loop seen in Figure 5-2 is replaced
with its closed loop equivalent.
188
Figure 5-10 Control block diagram of the voltage loop
The closed loop transfer function of the voltage loop,
(5-14)
(5-15)
(5-16)
By using
includes all the transfer functions that exist in the cascaded loop.
(5-17)
As performed in the current loop, characteristic equation of the closed loop transfer
function of the voltage loop,
(5-18)
Then, open loop transfer function of the voltage loop,
(5-19)
189
As it is seen clearly,
here. Therefore, it is
not possible to follow the procedure applied for the design of the current controller.
Instead of that approach, different zero and gain values will be tried until a proper
phase and gain margins are obtained in Matlab. While doing it, minimum phase
property of
already has a right half-plane zero pair. Its pole-zero map is shown in
Figure 5-11. At the left figure, all poles and zeros of the related transfer function is
shown but not seen easily. Hence, zoomed-in snapshot of the red circle is given at the
right figure. Note that
.
Figure 5-11 Pole-zero map of
190
As mentioned earlier, again PI controller is proposed as a controller. Its
transfer function can be written in the following manner.
(5-20)
where
and
(5-21)
(5-22)
Therefore,
is calculated as:
(5-23)
As a result, transfer function of the voltage controller turn out to be:
(5-24)
Now, minimum phase property must be checked in order to rely on bode plot.
Pole-zero map of
Figure 5-13 Bode plot of
192
As a final graph, bode plot of the closed loop transfer function is shown in Figure
5-14 as an information about frequency response of the controlled converter.
Figure 5-14 Bode plot of
5.4 DOMAIN CONVERSIONS
Controller functions have been obtained in s-domain. In order to implement
in a digital controller, these functions will be converted to their z-domain equivalent
by using bilinear transformation. Then, by utilizing inverse z-transform, they will be
converted to their discrete time domain equivalent.
Both of the controllers are in the following form;
(5-25)
193
There are few s-domain to z-domain transformation types. Bilinear transformation
has been applied and given a satisfactory result. Hence, it will be given in this
section. According to bilinear transformation, the following conversion is performed:
(5-26)
where
(5-27)
(5-28)
At this stage, input-output relationship of the controller can be formed. Let the
controller input be and output be . While they will be named as and
in z-domain, they are called and in discrete time domain. Hence, the
following relationship exists.
(5-29)
(5-30)
(5-31)
By using inverse z-transform, (5-31) is converted to its discrete time equivalent.
(5-32)
What does (5-32) mean? is the current output of the controller and is
the previous output. is the input to the controller at that sampling instant and
194
is the previous sampled input. That is, the current output of the PI-
controller is determined based on the previous output, the input at that sampling
instant and the previous input. As a consequence, (5-32) is ready to implement in a
digital signal controller.
5.5 APPLICATION SPECIFIC POINTS
Main ideas about controller design are given so far. This section includes
some practical, complementary or application specific points, which are important as
well.
As mentioned in the design chapter, the implemented circuit has two
selectable modes: constant input current mode (CICM) and constant input power
mode (CIPM). These modes have only one difference in terms control. When Figure
5-15 is investigated, it is seen that the controller outputs are limited. Limiter at the
output of
Generally, it is not possible to use the transfer functions directly. Because of
the voltage divisions in feedback paths, output range of the current sensors, analog-
to-digital conversion output range and PWM value range; transfer function is
normalized by a meaningful coefficients. For example, while duty-factor changes
between 0-0.9 in s-domain, it corresponds to 7-9240 bit value in the used digital
signal controller. Owing to these kinds of normalizations,
coefficients are
multiplied by 613 in this application.
As it will be given in details,
waveform at
waveform at
waveform at
waveform at
is exposed to
the input voltage in Mode 1 and the inverse of the output voltage in Mode 2. This is
an expected waveform. Also, the same waveform is observed on the combination of
. Since it is same as
is not given
intentionally. Instead,
are a little
bit lower than the cursors set at
waveform at
waveform at
in magnitude.
Figure 6-9
waveform at
waveform at
waveform at
waveform at
and
and
and
208
It is very low as expected. This value has been measured as 0.49% in the simulation.
Hence, it is considered as a consistent result.
Figure 6-15
, full load
Figure 6-16
and full load
Output inductor current waveform can be observed in Figure 6-17 and Figure 6-18.
The scale of the current probe is 5A/V in this measurement. Hence, mean value of
209
the output inductor is 10.2A, peak-to-peak ripple magnitude is 0.59A and ripple
percent is 5.7%. Peak-to-peak ripple magnitude has been calculated as 0.68A at this
operating point. The reason behind is the fact that the implemented/practical
has
turned out to be greater than the targeted/theoretical
and full load
210
6.2.2.1.2
and
and 50% load
Output current waveform at the first end operating point and 50% load is
presented in Figure 6-21 and Figure 6-22. Mean value of the output inductor current
is calculated as 5.05A, peak-to-peak ripple current magnitude as 0.545A and ripple
percent as 10.8%. For the output current, peak-to-peak ripple current magnitude is
expected to be the same for the same operating point regardless of the output power.
This value turns out to be 0.59A at full load case and 0.545A at 50% load case,
Figure 6-21
and 50% load
which are very close. Besides, since the mean value decreases while the ripple
magnitude remains the same, ripple current percent increases naturally. These results
are also regarded as consistent.
6.2.2.1.3
and
and 10% load
Output current waveform at the first end operating point and 10% load is
presented in Figure 6-25. Mean value of the output inductor current is calculated as
0.96A, peak-to-peak ripple current magnitude as 0.56A and ripple percent as 58.3%.
Note that the ripple magnitude is the same with the previous load conditions at this
end operating point. The same arguments at %50 load case are also valid here. No
extraordinary situation is observed.
214
Figure 6-25
and
and
and
and full load
Output current waveform at the fourth end operating point and full load is presented
in Figure 6-28 and Figure 6-29. Mean value of the output inductor current is
calculated as 7.2A, peak-to-peak ripple current magnitude as 1.22A and ripple
percent as 16.9%. Current ripple magnitude has been calculated as 1.47A in the
design section. The difference between 1.47A and 1.22A stems from the difference
between the implemented
and full load
217
6.2.2.2.2
and
and 50% load
Output current waveform at the fourth end operating point and 50% load is presented
in Figure 6-32 and Figure 6-33. Mean value of the output inductor current is
calculated as 3.58A, peak-to-peak ripple current magnitude as 1.19A and ripple
percent as 33.2%. Note that the ripple amplitude remains the same as expected.
There is no inconsistency observed in this operating point.
Figure 6-32
and 50% load
6.2.2.2.3
and
and 10% load
bottom of the waveform is a little bit higher than zero. This is due to the higher value
of implemented
Full
First 10,00 29,20 292,0 23,95 10,39 248,8 85,2
Second 10,01 29,10 291,3 34,04 7,38 251,2 86,2
Third 27,97 9,54 266,8 23,99 10,39 249,3 93,4
Fourth 28,03 9,45 264,9 34,03 7,38 251,1 94,8
Half
First 9,99 13,68 136,7 24,07 5,18 124,7 91,2
Second 9,99 13,61 136,0 34,03 3,68 125,2 92,1
Third 27,99 4,71 131,8 24,07 5,18 124,7 94,6
Fourth 27,99 4,70 131,6 34,05 3,68 125,3 95,2
10%
First 10,01 2,68 26,8 24,02 1,03 24,7 92,2
Second 10,01 2,70 27,0 34,11 0,73 24,9 92,1
Third 28,03 0,96 26,9 24,04 1,02 24,5 91,1
Fourth 27,98 0,98 27,4 34,02 0,73 24,8 90,6
6.3 CLOSED-LOOP RESULTS
In this section, closed loop performance of the implemented circuit and the
designed digital controller in (5-33)-(5-36) will be tested. As a disturbance, load
current change is considered. In order to organize the test results in mind, the
following reminders are given:
The converter has two selectable modes: constant input current mode (CICM)
and constant input power mode (CIPM). Each mode has two regulation modes.
CICM explanation: In this mode, if no-load current is drawn by the equipment,
the converter regulates its output voltage to 34V. As the load current increases,
the input current increases as well. As long as the input current is lower than
10A, converter output voltage will be regulated to 34V. At this point, input
power of the converter varies between 100W-280W according to the input
voltage. Once the input current is started to regulate at 10A, output voltage starts
to decrease. Its level is determined by the load current, but it is guaranteed by li-
ion battery pack that it will not decrease below 24V.
223
Each modes of regulation in CICM, namely output voltage regulation and input
current regulation, has its own dynamic response characteristics. Moreover,
transition between regulation types during dynamic response is possible, which
also has different dynamic response characteristics. Therefore, the test results are
grouped into:
Output voltage regulation mode
Input current regulation mode
Transition between regulation modes
Test results are given at the transitions between the farthest possible points in
each mode. Besides, since the input voltage is an important variable here, the
results are repeated at the two end input voltages, namely 10V and 28V.
CIPM explanation: The same arguments in CICM are also valid here. The only
difference is the input current limit. While the input current limit is constant at
10A in CICM, it is determined as (280W/
and
, lasts
2.2usec. Just after ADC conversion, digital PI controller calculations are performed.
It takes about 7.8usec. At , new duty-factor is determined. Although new
PWM cycle starts 1usec earlier, new duty-factor immediately updated at that instant
and is applied at that cycle. In this way, cycle by cycle control is achieved.
Figure 6-40 Timing diagram of the control loop in CICM
6.3.1.1 Output Voltage Regulation Mode
In this section, transient response of the two output variables,
and
,
against step change in
will be given.
The maximum output currents represent the load currents at which the input current
rises to just below 10A.
From now on, red waveforms represent
.
Besides, unless otherwise noted, current probe scale is at 5A/V. It will be at 50A/V
only at two groups of graphs, where the input current is high.
The parameters defining the performance of the transient response such as
rise time, settling time, maximum percent overshoot are not used regularly in these
test results. Because the regulated variable does not deviate much from its reference
value expect the exaggerated load current changes. In order to define these
parameters, a reference band of 2% or %5 is determined beforehand. Even if this
band is determined, these parameters may not carry sufficient information in some
graphs. Hence, instead of using these parameters, presenting raw data is considered
to be more meaningful in most case.
In Figure 6-41, the first snapshot is given. In order to get acquainted with
these graphs, some explanations will be presented here. The step load current change
from 0.1A to 2.6A occurs at the second time division from the left of the figures. It is
realized at the input voltage of 10V. At
is calculated as:
Similarly, at
is at
the bottom of the voltage divisions. The reference of
is brought to
zero, which is actually 34V as explained earlier with the name of tricky solution.
When
changes,
(green) and
Decrease in
when
decreases at that instant unexpectedly. When the pole zero map of the output current-
to-input current transfer function,
The inverse of this output current step change at the same operating point is
given in Figure 6-43. This time
and
230
and zero settling time. Additionally, it takes about 30msec for
. The output
power changes 29.2 folds, between 8.5W (3.4%) and 248.2W (99.28%). Note that
this much change will not occur in the application. At most 25W to 250W output
power change will be encountered. These widest possible changes are observed just
in order to see the performance of the controller.
In Figure 6-44, it is seen that
and
Inverse operation of Figure 6-44 is presented in Figure 6-45.
follows the
symmetric attitude of the previous result except the overshoot. No overshoot is
observed at all in
and
6.3.1.2 Input Current Regulation Mode
In this section, again, transient response of
and
is
still 10A and
and
.
In Figure 6-46, it is seen that when
increases,
is expected to decrease
instantaneously (if non-minimum phase property is considered) and then increase
above 10A. However, it is decremented to 10A by control action. As shown in Figure
6-46,
decreases from 10A to 9.2A instantly, then rises to 10.4A (4% overshoot)
and stays there until
reaches 24V,
settles to 10A.
The same waveform is also observed in closed-loop simulations. Although the
232
current reference is formed correctly at 10A,
does not reach 24V. The reason may stem from the coupling between the
inductors or the control action itself.
Figure 6-46
and
The inverse action is introduced in Figure 6-47. The same arguments given in the
Figure 6-47
and
233
previous figure also apply here. Note that the output power is about 91W in these
two figures.
Previous two measurements are repeated at
changes,
and
As it is clearly seen in Figure 6-49,
and
6.3.1.3 Transition between the Regulation Modes
In this section, the output current will be changed between no-load current
and the current at which the output voltage decreases down to its lower limit, 24V.
That is,
will be changed between the farthest points among two regulation modes.
Therefore, regulation mode will change during these transitions. Different dynamic
responses are expected.
In Figure 6-50 and Figure 6-51, output power is changed stepwise between
3.4W (1.36%) and 129.2W (51.68%) at
increases to 10A at about 2 msec with a negligible overshoot. As expected,
falls
down to 24V freely. When the load is decreased,
is forced
to be 10A until
) and output (
)
capacitor during 30msec. Then,
increases
from 24V to 34V in 30msec with 500mV overshoot, which is within 2% reference
band.
235
Figure 6-50
and
Figure 6-51
and
Previous two measurements are repeated at
is 7.5msec. Since
the input voltage is high,
and
is charged faster.
Figure 6-52
and
Figure 6-53
and
6.3.2 Constant Input Power Mode
Test results are classified at the following subsections. At first, the timing
diagram of the control loop will be given. This time, the selected dsPIC makes it
237
possible to refresh the duty-factor at a rate of 3 times the switching frequency,
namely 30usec. PWM signal (red) and timing signal (blue) are shown in Figure 6-64.
Similarly, at , PWM signal goes high. After 1usec, ADC conversion starts
in synchronous with PWM signal. Then, at ADC conversion is
completed. ADC conversion of 3 channels, by the addition of
to the other 2
channels, lasts 3.5usec. Digital PI controller calculations take about 24.8usec.
Because of the division of (280W/
decreases until
decreases 3.6V
(10.5%). It returns to its reference band in 40msec and settles exactly in 50msec.
Note that the current response is nearly three times faster in CICM than in CIPM.
Voltage response is directly related to the current response as seen in Figure 6-55.
Since the way that will be taken by
to
settle increases about that much.
and
and
When the load current is decreased,
and
as a voltage
deviation.
and
Similar waveforms are observed at
is increased,
is
decreased,
and
240
Figure 6-58
and
6.3.2.2 Input Power Regulation Mode
In this mode, rather than the output voltage, input power will be regulated to
280W. Since the input voltage is constant, what is regulated actually is the input
current. Lower
values are the output currents at which the output voltage is just
below 34V. Similarly, upper
and
241
In Figure 6-59 and Figure 6-60, current probe scale is 50A/V. Input power is
regulated to 274W. When
is increased,
is decreased,
deviates 1.6A (6.1%). It has the same timing values. Note that Figure 6-59 and
Figure 6-60 are highly symmetrical.
Figure 6-60
and
In Figure 6-61 and Figure 6-62, the previous measurements are repeated at
is increased,
deviates
0.8A (8%). It returns to its reference band in 40msec and settles exactly in 50msec.
When
is decreased,
and
Figure 6-62
and
6.3.2.3 Transition between the Regulation Modes
In this section, the minimum
and
Figure 6-64
and
244
also holds here. When
is increased,
settles to 24V band at the same time. Just after
is decreased;
rises to 34V in
6msec, makes 2.2V (6.4%) overshoot, enters its reference band in 75msec and settles
exactly in 100msec.
and
Figure 6-66
and
245
When the same transitions are repeated at
is increased;
is decreased;
rises to
34V in 8.5msec, makes 0.8V (2.3%) overshoot, enters its reference band in 25msec
and settles exactly in 50msec.
and
and
is
provided by adding an adjustment inductor to