State-Space Modelling of LLC Resonant Half-Bridge
State-Space Modelling of LLC Resonant Half-Bridge
State-Space Modelling of LLC Resonant Half-Bridge
Research Article
Abstract: Recently, LLC resonant converters have attracted significant research from industry and academia for AC–DC and
DC–DC power conversion with high efficiency and remarkable power density. They are appealing candidates for numerous
vehicular and renewable energy applications including battery chargers for electric vehicles and drivers of LED lights. This study
introduces a mathematical model of LLC resonant half-bridge DC–DC converter, which captures its steady-state behaviours for
both continuous conduction mode and discontinuous conduction mode operations. One major advantage of the proposed model
lies in accurate estimation of the switching frequency of power switches under a wide range of parametric variations. This
benefit is, however, not offered by the prevailing method based on the first harmonic approximation (FHA). The analytical
derivations of the system's state-space model, as well as equations for calculating the switching frequency by FHA, are
discussed in details. For illustration, a 340 W digitally controlled LLC resonant converter is targeted in this study. The simulation
analyses of current and voltage waveforms for light and heavy load conditions are presented. Moreover, the experimental
results, along with the comparison of switching frequency estimation for both methods, are demonstrated and discussed, which
confirms the validity and effectiveness of the proposed model.
ẋ = A2 x + B2u (5)
f s = kp k(vc t − V ref + ki∫ (kv c t − V ref )dt + f o (1)
where
In (1), k is the voltage scaler (k = Ro1 /(Ro1 + Ro2)), kp and ki are the
0 e21 0 0 0 0 0
proportional and integral gains, respectively, vc is the voltage
e22 0 0 0 0 b21 0
across the output capacitor CL, V ref is the reference voltage, f o is
the centre frequency of VCO. Let vcr (voltage across the capacitor A2 = e23 0 0 0 0 , B2 = b22 0 (6)
Cr), ir (current of the resonant tank), im (magnetising current), vc, 0 0 0 0 0 0 b23
and the feedback variable α (α = ki∫ kvc(t) − V ref dt) be the state 0 0 0 e24 0 0 b24
variables, x be the state vector, i.e. x = vcr ir im vc α T, and u be
The expressions for e21 ∼ e24 and b21 ∼ b24 are provided in the
the input vector, i.e. u = V in V ref T, then the system equation of the
Appendix. Further examination of the time-domain waveforms
converter is given by the following equation: reveals that there exists an underlying symmetric relationship
where W5 = 0 0 0 k 0 , W6 = 0 0 0 0 1 (14)
By solving (7)–(9), (11), and (12), one can find the steady-state where
solutions of t1, t2, and xt0. Therefore, the switching frequency f s of
the power converter can then be expressed by the following 0 e11 0 0 0 0 0
equation: e12 0 0 0 0 b11 −b12
A1′ = e13 0 0 0 0 , B1′ = b13 −b14 (18)
1
fs = (15) 0 −e14 −e15 0 0 0 b15
2t2
0 0 0 e16 0 0 b16
2.2 Continuous conduction mode
0 e11 0 0 0 0 0
The ideal time-domain waveforms of ir, im, id, vGS1, vGS2, and vDS2 for
e12 0 0 0 0 b11 b12
CCM operation are demonstrated in Fig. 4. Fig. 5 shows the
corresponding equivalent circuits for the periods A2′ = e13 0 0 0 0 , B2′ = b13 b14 (19)
t0 ∼ t1, t1 ∼ t2, t2 ∼ t3, and t3 ∼ t4. Similar to the case of DCM 0 e14 e15 0 0 0 b15
operation, symmetry can still be found between period t0 ∼ t2 and 0 0 0 e16 0 0 b16
t2 ∼ t4. The system equations for periods t0 ∼ t1 and t1 ∼ t2 are given
by the following equation: Then the mapping equations for xt0 to xt1 and xt1 to xt2 are given by
′ the following equation:
Period t0 ∼ t1 : ẋ = A1 x + B1′ u (16)
t1
Period t1 ∼ t2 : ẋ = A2′ x + B2′ u (17)
′
xt1 = e A1t1 xt0 + ∫ 0
′
e A1(t1 − σ)dσB1′ u (20)
Lr /Cr
3 Switching frequency estimation from FHA- Qe = (31)
Re
based method
The FHA-based method is being widely used in modelling and Consequently, (29) can be rewritten as
design process of LLC resonant converter. To compare with the
proposed model, the analytical equations for calculating the Ln f n2
switching frequency using FHA are derived here. Fig. 6 shows the Ge = (32)
Ln + 1 f n − 1 + j ( f n2 − 1) f nQeLn
2
linear equivalent circuit of LLC resonant converter. The pole
frequency f p and series resonant frequency f r are given by the
From (32), the function of f n can found as
following equations:
1 A1 f n6 + A2 f n4 + A3 f n2 + A4 = 0 (33)
fp = (22)
2π (Lr + Lm)Cr
where
1
fr = (23) A1 = Ge2Ln2 Qe2 (34)
2π LrCr
The RMS value of the input voltage vie is A2 = Ge2 −2Ln2 Qe2 + Ln2 + 2Ln + 1 − Ln2 (35)
where V out and Iout are the output voltage and current, respectively.
4 Simulation results
Then the AC equivalent load resistance Re is
To proceed with the simulation analysis, a standalone LLC
V oe, rms 8n2 resonant half-bridge DC–DC converter is utilised as the system
Re = = 2 RL (27) model in this study. The model of this converter has been
Ioe, rms π
constructed in MATLAB–Simulink environment. The design
parameters of the converter are summarised in Table 2. The
The impedances of Cr, Lr, and Lm are given by the following
Simulink model of the converter system is shown in Fig. 7. Both
equation: the heavy load and light load conditions are considered in the
analysis. For the case of heavy load operation (iL = 23 A), Figs. 8
1
XCr = , XLr = ωLr, XLm = ωLm (28) and 9 display the time-domain waveforms of the current flowing in
ωCr
the resonant tank ir, magnetising current im, overall diode current id,
output voltage vout, pole voltage vDS2, and voltage vcr across Cr
Therefore, the input-to-output voltage gain Ge is
when the input voltage V in is set to 350 and 400 V.
XLmRe The same set of figures for light load operation (iL = 6.3A) is
Ge = (29) shown in Figs. 10 and 11, respectively. It is found that both CCM
XLm + XLr − XCr Re + jXLm(XLr − XCr)
and DCM operations of the converter can be expected for different
values of V in. Specifically, when V in is varied close to the lower
limit, the switching frequency of the power switches is regulated
and stays at a low level such that the voltage drop across the
resonant tank reduces accordingly. Therefore, the nominal value of
the output voltage vout can be maintained. Since the current flowing
in the secondary side of the transformer is discontinuous, the
converter operates in DCM. On the other hand, when V in is
changed to the upper limit, the switching frequency is increased by
the feedback controller and, due to this fact, im never reaches the
same level of ir. As a result, the overall diode current at the
secondary side of the transformer remains continuous; hence the
converter operates in CCM.
Fig. 6 Linear equivalent circuit of LLC resonant DC–DC converter
Fig. 10 Simulation results when V in is set to 350 V (iL = 6.3 A, DCM operation)
(a) Current waveforms, (b) Voltage waveforms
Fig. 11 Simulation results when V in is set to 400 V (iL = 6.3 A, CCM operation)
(a) Current waveforms, (b) Voltage waveforms
the output voltage is measured by the differential probe frequency f s is conducted and the simplified estimation algorithm
THDP0200. is shown in Fig. 15. As a comparison, the evaluations by FHA-
The obtained time-domain waveforms of vout, vDS2, ir, and vcr are based method are also included. Fig. 16 demonstrates the obtained
displayed in Figs. 13 and 14, respectively, for the light load (iL = results from FHA and the proposed model along with the
6.3 A) and heavy load (iL = 23 A) conditions. Apparently, the experimental data when V in varies from 350 to 400 V. Tables 3 and
resonant converter can exhibit both DCM and CCM operations 4 summarise the selected results of both methods for light load and
under different values of the input voltage, which is consistent with heavy load operations, respectively. It can be observed that when
the results from simulation analysis. Moreover, to illustrate the V in is in the range of 350–380 V, both FHA method and the
utility of the proposed model, the calculation of switching proposed model can produce reasonably accurate values that are
Fig. 13 Experimental waveforms for light load condition (iL = 6.3 A, CH1: output voltage vout (10 V/div), CH2: pole voltage vDS2 (100 V/div), CH3: current
of resonant tank ir (1 A/div), CH4: voltage across Cr, vcr (100 V/div))
(a) DCM (Vin = 350 V), (b) CCM (Vin = 400 V)
Fig. 14 Experimental waveforms for heavy load condition (iL = 23 A, CH1: output voltage vout (10 V/div), CH2: pole voltage vDS2 (100 V/div), CH3: current
of resonant tank ir (5 A/div), CH4: voltage across Cr, vcr (100 V/div)
(a) DCM (Vin = 350 V), (b) CCM (Vin = 400 V)
close to those from test measurement. However, as V in increases values from the proposed model remain in close proximity to the
beyond this range, the evaluated results of the FHA method deviate actual results; hence decent precision can still be retained by using
considerably from targeted values. This tendency is true when the this model. Although the estimation accuracy declines slightly as
converter operates under either heavy load or light load condition, the switching frequency increases, the relative errors are still
as can be seen from Fig. 16. On the other hand, the calculated within 5% for both load conditions, as illustrated in Tables 3 and 4.
Table 3 Comparison of switching frequency estimations for the light load operation
Input voltage, V Experimental results, kHz FHA method, kHz Proposed model, kHz
350 56.95 53.76 57.04
360 62.08 59.58 62.19
370 69.36 67.81 69.27
380 79.04 80.68 79.72
390 93.37 105.1 96.19
400 125.91 177.75 121.41
Table 4 Comparison of switching frequency estimations for the heavy load operation
Input voltage, V Experimental results, kHz FHA method, kHz Proposed model, kHz
350 51.38 49.68 50.40
360 56.14 55.58 55.62
370 62.18 63.83 62.22
380 69.72 76.25 70.90
390 79.88 96.46 82.99
400 94.67 128.83 96.47
1 Ls′ 1
e12 = − + (39) b21 = b22 = (55)
Lr Lr2k3 Lr + Lm
Ls′ 1
e13 = − (40) b23 = − (56)
LmLrk3 RLCLk
n b24 = − ki (57)
e14 = − e15 = (41)
CL