Digital Compensator Design For LLC Resonant Converter
Digital Compensator Design For LLC Resonant Converter
Q1
C1
D1
Vin D3
ir isp Io
A
Ls +
Q2 Cs im
rc
R V0
ns
C2 Lm Cf
np –
D2
ns
D4
B
R1
Driver PWM Digital
Compensator ADC
Transformers Out put
Q1 R2
Q2 dsPIC33FJXXGSXXX
rs Ls ir Cs ip np: ns : ns isp Io
A
+
+ Vcr – im D3 rc
R Vo
Cf Vcf
–
Lm V'cf
+ VAB
Vin –
D4
Ts
B
The state equations are obtained in Continuous Tank EQUATION 3: TRANSFORMER PRIMARY
Current mode by using Kirchhoff’s Circuit Laws (KCL), VOLTAGE
as shown in Equation 1 through Equation 4.
di m
v' c = L m --------
EQUATION 1: RESONANT TANK VOLTAGE
f dt
di r
v AB = L s ------- + i r r s + v cr + sgn ( i p )v'c f EQUATION 4: TRANSFORMER
dt
SECONDARY CURRENT
Where:
r dv c 1
sgn(ip) = { -1, if v’cf < 0 i sp = 1 + ---c- C f ---------f + --- v c
R dt R f
+1, if v’cf ≥ 0}
CURRENT Where:
r' c = r c R
dv cr
i r = C s ----------
dt
The primary side resonant tank parameters, ir(t), vc(t) The parameters, sine component of resonant
and im(t), provided in Equation 7, are approximated to current (is), cosine component of resonant current (ic),
their fundamental harmonics, and the output filter volt- sine component of resonant capacitor voltage (vs),
age (vcf) is approximated to the DC component. The cosine component of resonant capacitor voltage (vc),
derivatives of ir(t), vcr(t) and im(t) are shown in sine component of magnetizing current (ims) and cosine
Equation 7. component of magnetizing current (imc) are slow time
varying components. Therefore, the dynamic behavior
EQUATION 7: FUNDAMENTAL of these parameters can be analyzed.
APPROXIMATION OF Figure 3 and Figure 4 illustrate the simulation wave-
PRIMARY TANK forms of the LLC resonant converter operating below
PARAMETERS the resonant frequency and continuous tank current
mode.
i r ( t ) = i s ( t ) sin ω s t – i c ( t ) cos ω s t
v c r ( t ) = v s ( t ) sin ω s t – v c ( t ) cos ω s t
i m ( t ) = i ms ( t ) sin ω s t – i mc ( t ) cos ω s t
di r di di
------- = -------s + ω s i c sin ω s t – -------c – ω s i s cos ω st
dt dt dt
dvc r dv dv
- = --------s + ω s v c sin ω s t – --------c – ω s v s cos ω s t
----------
dt dt dt
di m di m s di mc
--------- = ----------
- + ωs i mc sin ωs t – -----------
- – ω s i m s cos ω s t
dt dt dt
Where:
ωs = switching frequency in radians/second
Time (ms)
Time (ms)
Time (ms)
Time (ms)
Time (ms)
sgn ( i sp ) v' c = f 2 i ss, i sp, v'c sin ω s t – f 3 i sc, i sp, v'c cos ω s t
f f f
The switching waveform has an odd symmetry. There-
i sp = f 4 ( i ss, i sc ) fore, there is no cosine component (vec = 0, where vec
is the cosine component of the output voltage of the
Figure 5 illustrates a typical switching waveform of a half-bridge inverter) in the switching waveform and the
half-bridge inverter which is the input to the LLC sine component (ves) forms the fundamental
resonant tank (θ = Dead Time, d = Duty Cycle). component of vAB .
Vin
θ θ
dπ
0 π 2π
ωt
is rs Ls + Vs - i ps
– +
Cs i ms
Ωs Ls ic Lm Io
+ + Vps
Ves Ωs Cs Vc – +
–
– rc
Ωs Lmi mc +
R
2/p isp Vo
Ω s Lmi ms – Cf Vcf
Vec – + – –
Ω s Cs V s
+ + Vpc
Ω s Ls is Lm
ic rs Ls imc ipc
Cs
+ –
+ Vc –
Deriving Steady-State Operating Point The output filter capacitor voltage can be calculated by
substituting Equation 16 into Equation 14, as shown in
Under steady-state conditions, the state variables of Equation 17.
the modulation equations, Equation 12 through
Equation 14, do not change with time. For a chosen
EQUATION 17: FILTER CAPACITOR
operating point, the time derivatives in Equation 12
VOLTAGE
through Equation 14 are set to zero and the steady-
vc
state values are obtained (shown in upper case letters).
------f = ------ i pp = --- i s p
2n 2
The transformer currents on the primary and secondary R π π
sides are shown in Equation 16. 2n
v c = ------ i pp R
f π
EQUATION 16: TRANSFORMER CURRENTS π
v c = ------ i pp R e
Primary Current: f 4n
2 2 π
i pp = ip s + ip c V c' f = nV c = --- I pp R e
f 4
8 2
Secondary Current: R e = -----2 n R = equivalent load resistance referred
π
2 2 2 2 to primary side
i sp = i ss + i sc = n i ps + i pc = ni pp
Where:
Where: V 'c f = reflected voltage of secondary on the
n = np/ns = transformer turns ratio primary
i pc = I pc + î pc , i
pp = I pp + î pp , v c = V c + v̂ c , i m s = I m s + î m s , i mc = I mc + î mc , i s = I s + iˆs ,
f f f
i c = I c + î c , v s = V s + v̂ s , v c = V c + v̂ c , v e s = V e s + v̂ e s , v 0 = V 0 + v̂ 0
4nV c 2 4nV c
f I pc f I ps I pc- 4n I ps
v̂ ps = -------------- ---------3 î ps – -------------- -------------- î p c + ------ ------- v̂ c
π I π I pp
3 π I pp f
pp
v̂ ps = H ip î ps + H ic î p c + H vcf v̂ c
f
v̂ ps = H ip î s + H ic î c – H ip î'ims
ms – H ic î mc + H vcf v̂ c f
Where:
4nV c 2
f I pc
H ip = -------------- ---------3
π I
pp
4nV c
f I ps I pc-
H ic = – -------------- --------------
π I pp
3
4n I ps
H vcf = ------ -------
π I pp
I 2 + I 2 – ------------------------------ 1
ps pc
4nV c I I 4nV c I
2
+ I
2 I pc
f f
= – -------------- ---------------------------------3- î ps + -------------- -------------------------------------------------------------------
ps pc ps pc 4n ------------------------------ v̂
v̂ pc î + -----
-
π --- π 2 pc π cf
2 2 2 I 2
+ I
2 I ps2 + I pc2
( I ps + I pc ) ps pc
4nV c 4nV c 2
f I ps I pc- f I----------- I pc
v̂ pc = – -------------- --------------- î + -------------
- pc
- î pc + 4n
------ -------- v̂
π I pp
3 ps π I pp
3 π I pp c f
v̂ pc = G ip î s + G ic î c – G ip î ms – G ic î mc + G vcf v̂ cf
Where:
4nV c
f I ps I pc-
G ip = – -------------- ---------------
π I pp
3
4nV c I 2
f psc -
G ic = -------------- -----------
π I pp
3
4n I pc
G vcf = ------ -------
π I pp
The linearization of the input voltage (ves) is shown in Removing the steady-state terms and other higher
Equation 29. order perturbed terms in Equation 29 to get the
linearized input voltage is shown in Equation 30.
EQUATION 29: LINEARIZATION OF
HALF-BRIDGE INVERTER EQUATION 30: LINEARIZATION OF
OUTPUT VOLTAGE HALF-BRIDGE INVERTER
2v in OUTPUT VOLTAGE
π
v es = ---------- sin --- d
π 2 π π
v̂ es = --- sin --- D v̂ in + V in cos --- D dˆ
2
π 2 2
π
v̂ es = --- × ( V in + v̂ in ) sin --- ( D + dˆ )
2
π 2 v̂ es = K 1 v̂ in + K 2 dˆ
π
Expanding sin --- ( D + dˆ ) Where:
2
π
K 1 = --- sin --- D
2
π π π π
= sin --- D cos --- dˆ + cos --- D sin --- dˆ π 2
2 2 2 2
π
K 2 = V in cos --- D
2
π π π
= sin --- D + --- cos --- D dˆ
2 2 2
π π π
v̂ es = --- ( V in + v̂ in ) sin --- D + --- cos --- D dˆ
2
π 2 2 2
dî c ( Ω s L s – G ip ) ( G ic + r s ) 1 G ip G ic G vc f Ls ω0 Is
------- = ------------------------------- î s – ----------------------- iîcc – ----- v̂ c + -------- î ms + -------- î mc – ----------- v̂ c + ----------------- ω̂ sn
dt Ls Ls Ls Ls Ls Ls f Ls
K is î ps + K ic î pc
î sp = K is î s + K ic î c – K is î ms – K ic î mc
dv̂ c
rc f
1 + ----- C f × ---------- = K is î s + K ic î c – K is î ms – K ic î mc – --- v̂ c
1
R dt R f
Where:
i ps = i s – i m s and i pc = i c – i mc
2n I ps I pc
2n
K is = ------ --------------------------------- and K ic = ------ ---------------------------------
π π
I ps 2 + I pc 2 2
I ps + I pc
2
rc dv̂ c
f
------ C f + ---------- = K is î s + K ic î c – K is î ms – K ic î mc – --- v̂ c
1
r' c dt R f
dv̂ c K is r' c K is r' K is r' c K ic r' c r' c
ic c
----------f = --------------- î s + --------------- î c – --------------- î ms – --------------- î mc – --------------- v̂ c
dt Cf rc Cf rc Cf rc Cf rc RCf rc f
dî ( Ω s L s – G ip ) ( G ic + r s ) G ip G ic G vcf Ls ω0 Is
-------c = --------------------------------- î s – ------------------------- î c – ----- v̂ c + --------- î ms + --------- î mc – ------------ v̂ cf + ------------------ ω̂ sn
1
dt Ls Ls Ls Ls Ls Ls Ls
dv̂ s Cs Ω C s ω Vc
1 s 0
-------- = ------ î s – -------------- v̂ c – -------------------- ω̂ sn
dt Cs Cs Cs
dv̂ c Cs Ω C s ω Vs
1 s 0
-------- = ------ î c + -------------- v̂ s + -------------------- ω̂ sn
dt Cs Cs Cs
dî ms H ip H ic H ip H ic + L m Ω s H vcf L m I mc ω
0
- = --------- î s + --------- î c – --------- î ms – ------------------------------- î mc + ------------ v̂ cf – ------------------------- ω̂ sn
-----------
dt Lm Lm Lm Lm Lm Lm
dî mc G ip G ic ( G ip – L m Ω s ) G ic G vcf L m I ms ω
0
- = --------- î s + --------- î c – ----------------------------------- î ms – --------- î mc + ------------ v̂ cf + ------------------------- ω̂ sn
-----------
dt Lm Lm Lm Lm Lm Lm
r'c
v̂ 0 = K is r'c î s + K ic r'c î c – K is r'c î ms – K ic r'c î mc + ------ v̂ cf
rc
ŷ = ( v̂ 0 ) Output
H ip + r s ( Ω s L s + H ic ) H ip H ic H vcf
1 --------- ---------
– -------------------- – --------------------------------- – ----- 0 – ------------
Ls Ls Ls Ls Ls Ls
( Ω s L s – G ip ) G ic + r s G ip G ic G vcf
--------------------------------- 1 --------- ---------
– ------------------- - 0 – ----- – ------------
Ls Ls Ls Ls Ls Ls
1 Cs Ωs
------ 0 0 – -------------
- 0 0 0
Cs Cs
1 Cs Ωs
A = 0 ------ -------------- 0 0 0 0
Cs Cs
H ip H ic H ip H ic + L m Ω s H vcf
--------- --------- 0 0 – -------- - - ------------
– ------------------------------
Lm Lm Lm Lm Lm
G ip G ic G ip – L m Ω s G ic G vcf
--------- --------- 0 0 – ------------------------------
- – --------- ------------
Lm Lm Lm Lm Lm
K is r' c K ic r' c K is r' c K ic r' c r' c
--------------- --------------
- 0 0 – -------------- - – -------------- - – --------------
-
Cf rc Cf rc Cf rc Cf rc RC f rc
T
L s ω 0 I c L s ω 0 I s C s ω 0 V c C s ω 0 V s L m ω 0 I mc L m ω 0 I ms
B = – ----------------
- ---------------- – ------------------- ------------------- – ---------------------- ---------------------- 0
Ls Ls Cs Cs Lm Lm
r'
c
C = ( K is r' c ) ( K i c r' c ) ( 0 ) ( 0 ) ( – K is r' c ) ( – K ic r' c ) ------
rc
D = 0
For the linearized system, the required control-to-output voltage transfer function is:
v̂ 0 –1
- = C ( SI – A ) B + D = G p ( s )
---------
ω̂ sn
s s
5.5895 1 + ---------------------------- - × ---------------------------- - – 1
3.314 × 10 8.262 × 10
4 5
G p ( s ) = ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------
2 2 5
s 973.6s s 2.75 × 10 s
----------------------------- + ----------------------------- + 1 × -------------------------------- + -------------------------------- + 1
8.949 × 10 8
8.949 × 10
8 1.227 × 10 12
1.227 × 10
12
The general form of Gp(s) is shown in Equation 46. The [p, z] = pzmap (Gp (s)) command gives the poles
and zeros of the plant transfer function. Figure 7 illus-
EQUATION 46: GENERALIZED FORM OF trates the pole-zero plot for the Gp (s), which is obtained
PLANT TRANSFER from the MATLAB command, pzmap (Gp (s)). Figure 8
FUNCTION illustrates the bode plot obtained from the hardware
using the network analyzer. Figure 9 illustrates the
s s bode plot obtained using MATLAB.
G po 1 + ----------- - × ---------------- – 1
ω esr ω RHP As illustrated in Figure 8 and Figure 9, the bode plot
G p ( s ) = ------------------------------------------------------------------------------------------------------------------------ captured, using the network analyzer, matches the
2 2
--------- s
s - + ----------------------- × --------- s
s - + ----------------------- analytical bode plot obtained in MATLAB, thereby,
- + 1 - + 1
2 Q 1 × ω p1 2 Q 2 × ω p2 confirming the veracity of the mathematical model.
p1
ω p2 ω
Digital Compensator Design for LLC In this method, an analog compensator is first designed
Resonant Converter in the continuous time domain and then converted to
discrete time domain using bilinear or tustin transfor-
The plant model is derived as shown in Equation 45. In mation. Figure 10 illustrates the block diagram of the
order to attain the desired gain margin, phase margin LLC resonant converter with a digital compensator.
and crossover frequency, a digital 3P3Z compensator
is designed. The digital 3P3Z compensator is derived
using the design by emulation or digital redesign
approach.
Resonant Converter
VREF [n] + e[n] 3P3Z Vc[n] F(t) Vo[t]
DPWM
Compensator
–
As seen from Equation 46, plant transfer function con- Furthermore, in order to compensate for the effect of
sists of an ESR zero and a pair of dominant complex the complex dominant poles (reduction in system
poles. In order to compensate for the effect of ESR zero damping, and hence, increased overshoots and set-
(increased high-frequency gain, and thereby, tling time), two zeros, (s+a+jb) and (s+a-jb), are added.
increased ripple), a pole (ωp) is included in the compen- Also, to achieve sufficient attenuation at switching
sator. In order to minimize the steady-state error, an frequency, a pole is added to the compensator at half
integrator (Kc) is also added to the compensator. the switching frequency.
The discrete compensator transfer function (Gc_d) is EQUATION 50: COMPENSATOR TRANSFER
obtained using the tustin or bilinear transformation with FUNCTION IN DISCRETE
a sampling frequency of 50 kHz, as shown in DOMAIN
Equation 50.
3 2
0.2711 × z – 0.178 × z – 0.1828 × z + 0.2663
Gc_d = -----------------------------------------------------------------------------------------------------------------
3 2
z – 0.6791 × z – 0.7342 × z + 0.4133
• Microchip believes that its family of products is one of the most secure families of its kind on the market today, when used in the
intended manner and under normal conditions.
• There are dishonest and possibly illegal methods used to breach the code protection feature. All of these methods, to our
knowledge, require using the Microchip products in a manner outside the operating specifications contained in Microchip’s Data
Sheets. Most likely, the person doing so is engaged in theft of intellectual property.
• Microchip is willing to work with the customer who is concerned about the integrity of their code.
• Neither Microchip nor any other semiconductor manufacturer can guarantee the security of their code. Code protection does not
mean that we are guaranteeing the product as “unbreakable.”
Code protection is constantly evolving. We at Microchip are committed to continuously improving the code protection features of our
products. Attempts to break Microchip’s code protection feature may be a violation of the Digital Millennium Copyright Act. If such acts
allow unauthorized access to your software or other copyrighted work, you may have a right to sue for relief under that Act.
QUALITY MANAGEMENT SYSTEM Microchip received ISO/TS-16949:2009 certification for its worldwide
headquarters, design and wafer fabrication facilities in Chandler and
CERTIFIED BY DNV Tempe, Arizona; Gresham, Oregon and design centers in California
and India. The Company’s quality system processes and procedures
== ISO/TS 16949 ==
are for its PIC® MCUs and dsPIC® DSCs, KEELOQ® code hopping
devices, Serial EEPROMs, microperipherals, nonvolatile memory and
analog products. In addition, Microchip’s quality system for the design
and manufacture of development systems is ISO 9001:2000 certified.