前言
在串联谐振型DAB电路拓扑中,稳态运行时,原本的直流输入和副边的直流输出往往已经达到稳定。此时去原电路等效看作一个直流源和LC串联的电路(如下图)。直流源的电压大小与原副边之间的移相比有关。此博文将详细推导该等效电路的状态轨迹( i L i_L iL- v C v_C vC)。
涉及学科知识:微分方程,电路理论,电力电子
一、基本的数学原理

如图所示,为等效的电路模型。
微分方程的列写和通解形式
其中
V
V
V的大小与原副边的电压和移相比有关。由该电路可以列写出如下微分方程组:
V
=
L
d
i
L
d
t
+
v
C
V=L\frac{d{{i}_{L}}}{dt}+{{v}_{C}} \\
V=LdtdiL+vC
i
L
=
C
d
v
C
d
t
{{i}_{L}}=C\frac{d{{v}_{C}}}{dt} \\
iL=CdtdvC
联立可得
V
=
L
C
d
2
v
C
d
t
2
+
v
C
V=LC\frac{{{d}^{2}}{{v}_{C}}}{d{{t}^{2}}}+{{v}_{C}}
V=LCdt2d2vC+vC
特征方程为
L
C
λ
2
+
1
=
0
LC{{\lambda }^{2}}+1=0 \\
LCλ2+1=0
⇒
λ
=
j
1
L
C
≡
j
ω
r
\Rightarrow \lambda =j\frac{1}{\sqrt{LC}}\equiv j{{\omega }_{r}} \\
⇒λ=jLC1≡jωr
其中,
ω
r
{{\omega }_{r}}
ωr 为谐振频率。
因此,假设初始时刻为
t
0
{{t}_{0}}
t0,原方程的解的形式为
v
C
(
t
)
=
C
1
cos
(
ω
r
(
t
−
t
0
)
)
+
C
2
sin
(
ω
r
(
t
−
t
0
)
)
+
V
{{v}_{C}}\left( t \right)={{C}_{1}}\cos \left( {{\omega }_{r}}\left( t-{{t}_{0}} \right) \right)+{{C}_{2}}\sin \left( {{\omega }_{r}}\left( t-{{t}_{0}} \right) \right)+V \\
vC(t)=C1cos(ωr(t−t0))+C2sin(ωr(t−t0))+V
i
L
(
t
)
=
−
ω
r
C
C
1
sin
(
ω
r
(
t
−
t
0
)
)
+
ω
r
C
C
2
cos
(
ω
r
(
t
−
t
0
)
)
{{i}_{L}}\left( t \right)=-{{\omega }_{r}}C{{C}_{1}}\sin \left( {{\omega }_{r}}\left( t-{{t}_{0}} \right) \right)+{{\omega }_{r}}C{{C}_{2}}\cos \left( {{\omega }_{r}}\left( t-{{t}_{0}} \right) \right) \\
iL(t)=−ωrCC1sin(ωr(t−t0))+ωrCC2cos(ωr(t−t0))
微分方程解的形式见之前的博客(二阶线性微分方程求解链接: 二阶线性微分方程求解)
电容电压和电感电流的初值定义为
i
L
(
t
0
)
=
i
L
(
t
0
)
{{i}_{L}}\left( {{t}_{0}} \right)={{i}_{L}}\left( {{t}_{0}} \right) \\
iL(t0)=iL(t0)
v
C
(
t
0
)
=
v
C
(
t
0
)
{{v}_{C}}\left( {{t}_{0}} \right)={{v}_{C}}\left( {{t}_{0}} \right) \\
vC(t0)=vC(t0)
带入上式得
v
C
(
t
0
)
=
C
1
+
V
{{v}_{C}}\left( {{t}_{0}} \right)={{C}_{1}}+V \\
vC(t0)=C1+V
i
L
(
t
0
)
=
+
ω
r
C
C
2
{{i}_{L}}\left( {{t}_{0}} \right)=+{{\omega }_{r}}C{{C}_{2}} \\
iL(t0)=+ωrCC2
⇒
C
1
=
v
C
(
t
0
)
−
V
\Rightarrow {{C}_{1}}={{v}_{C}}\left( {{t}_{0}} \right)-V \\
⇒C1=vC(t0)−V
⇒
C
2
=
i
L
(
t
0
)
ω
r
C
\Rightarrow {{C}_{2}}=\frac{{{i}_{L}}\left( {{t}_{0}} \right)}{{{\omega }_{r}}C} \\
⇒C2=ωrCiL(t0)
综上,可以解得电感电流和电容电压的通解为
v
C
(
t
)
=
(
v
C
(
t
0
)
−
V
)
cos
(
ω
r
(
t
−
t
0
)
)
+
(
i
L
(
t
0
)
ω
r
C
)
sin
(
ω
r
(
t
−
t
0
)
)
+
V
{{v}_{C}}\left( t \right)=\left( {{v}_{C}}\left( {{t}_{0}} \right)-V \right)\cos \left( {{\omega }_{r}}\left( t-{{t}_{0}} \right) \right)+\left( \frac{{{i}_{L}}\left( {{t}_{0}} \right)}{{{\omega }_{r}}C} \right)\sin \left( {{\omega }_{r}}\left( t-{{t}_{0}} \right) \right)+V \\
vC(t)=(vC(t0)−V)cos(ωr(t−t0))+(ωrCiL(t0))sin(ωr(t−t0))+V
i
L
(
t
)
=
−
ω
r
C
(
v
C
(
t
0
)
−
V
)
sin
(
ω
r
(
t
−
t
0
)
)
+
i
L
(
t
0
)
cos
(
ω
r
(
t
−
t
0
)
)
{{i}_{L}}\left( t \right)=-{{\omega }_{r}}C\left( {{v}_{C}}\left( {{t}_{0}} \right)-V \right)\sin \left( {{\omega }_{r}}\left( t-{{t}_{0}} \right) \right)+{{i}_{L}}\left( {{t}_{0}} \right)\cos \left( {{\omega }_{r}}\left( t-{{t}_{0}} \right) \right) \\
iL(t)=−ωrC(vC(t0)−V)sin(ωr(t−t0))+iL(t0)cos(ωr(t−t0))
归一化准备:
归一化首先需要锚定一个基准值,这里设置电压的基准值为原边的输入电压
V
i
n
{{V}_{in}}
Vin ,LC串联谐振电路的特征阻抗为
Z
Z
Z,
Z
=
L
C
Z=\sqrt{\frac{L}{C}}
Z=CL ,因此,电流的基准值为
V
i
n
Z
\frac{{{V}_{in}}}{Z}
ZVin
因此有以下关系:
v
C
(
t
)
=
V
C
N
V
i
n
{{v}_{C}}\left( t \right)={{V}_{CN}}{{V}_{in}} \\
vC(t)=VCNVin
v
C
(
t
0
)
=
V
C
0
N
V
i
n
{{v}_{C}}\left( {{t}_{0}} \right)={{V}_{C0N}}{{V}_{in}} \\
vC(t0)=VC0NVin
V
=
V
N
V
i
n
V={{V}_{N}}{{V}_{in}} \\
V=VNVin
i
L
(
t
)
=
I
L
N
V
i
n
Z
{{i}_{L}}\left( t \right)={{I}_{LN}}\frac{{{V}_{in}}}{Z} \\
iL(t)=ILNZVin
i
L
(
t
0
)
=
I
L
0
N
V
i
n
Z
{{i}_{L}}\left( {{t}_{0}} \right)={{I}_{L0N}}\frac{{{V}_{in}}}{Z} \\
iL(t0)=IL0NZVin
ω
r
(
t
−
t
0
)
=
θ
{{\omega }_{r}}\left( t-{{t}_{0}} \right)=\theta \\
ωr(t−t0)=θ
代入微分方程的通解可得
V
C
N
V
i
n
=
(
V
C
0
N
−
V
N
)
V
i
n
cos
(
θ
)
+
I
L
0
N
V
i
n
Z
ω
r
C
sin
(
θ
)
+
V
N
V
i
n
{{V}_{CN}}{{V}_{in}}=\left( {{V}_{C0N}}-{{V}_{N}} \right){{V}_{in}}\cos \left( \theta \right)+{{I}_{L0N}}\frac{{{V}_{in}}}{Z{{\omega }_{r}}C}\sin \left( \theta \right)+{{V}_{N}}{{V}_{in}} \\
VCNVin=(VC0N−VN)Vincos(θ)+IL0NZωrCVinsin(θ)+VNVin
I
L
N
V
i
n
Z
=
−
ω
r
C
(
V
C
0
N
−
V
N
)
V
i
n
sin
(
θ
)
+
I
L
0
N
V
i
n
Z
ω
r
C
cos
(
θ
)
{{I}_{LN}}\frac{{{V}_{in}}}{Z}=-{{\omega }_{r}}C\left( {{V}_{C0N}}-{{V}_{N}} \right){{V}_{in}}\sin \left( \theta \right)+{{I}_{L0N}}\frac{{{V}_{in}}}{Z{{\omega }_{r}}C}\cos \left( \theta \right) \\
ILNZVin=−ωrC(VC0N−VN)Vinsin(θ)+IL0NZωrCVincos(θ)
⇒
V
C
N
=
(
V
C
0
N
−
V
N
)
cos
(
θ
)
+
I
L
0
N
1
Z
ω
r
C
sin
(
θ
)
+
V
N
\Rightarrow {{V}_{CN}}=\left( {{V}_{C0N}}-{{V}_{N}} \right)\cos \left( \theta \right)+{{I}_{L0N}}\frac{1}{Z{{\omega }_{r}}C}\sin \left( \theta \right)+{{V}_{N}} \\
⇒VCN=(VC0N−VN)cos(θ)+IL0NZωrC1sin(θ)+VN
⇒
I
L
N
1
Z
=
−
ω
r
C
(
V
C
0
N
−
V
N
)
sin
(
θ
)
+
I
L
0
N
1
ω
r
C
cos
(
θ
)
\Rightarrow {{I}_{LN}}\frac{1}{Z}=-{{\omega }_{r}}C\left( {{V}_{C0N}}-{{V}_{N}} \right)\sin \left( \theta \right)+{{I}_{L0N}}\frac{1}{{{\omega }_{r}}C}\cos \left( \theta \right) \\
⇒ILNZ1=−ωrC(VC0N−VN)sin(θ)+IL0NωrC1cos(θ)
又因为
ω
r
C
=
1
L
C
=
C
L
=
1
Z
{{\omega }_{r}}C=\sqrt{\frac{1}{LC}}=\sqrt{\frac{C}{L}}=\frac{1}{Z} \\
ωrC=LC1=LC=Z1
1
Z
ω
r
C
=
1
Z
1
Z
=
1
\frac{1}{Z{{\omega }_{r}}C}=\frac{1}{Z\frac{1}{Z}}=1 \\
ZωrC1=ZZ11=1
因此,微分方程可以进一步化简为
V
C
N
=
(
V
C
0
N
−
V
N
)
cos
(
θ
)
+
I
L
0
N
1
Z
ω
r
C
sin
(
θ
)
+
V
N
{{V}_{CN}}=\left( {{V}_{C0N}}-{{V}_{N}} \right)\cos \left( \theta \right)+{{I}_{L0N}}\frac{1}{Z{{\omega }_{r}}C}\sin \left( \theta \right)+{{V}_{N}} \\
VCN=(VC0N−VN)cos(θ)+IL0NZωrC1sin(θ)+VN
I
L
N
1
Z
=
−
ω
r
C
(
V
C
0
N
−
V
N
)
sin
(
θ
)
+
I
L
0
N
cos
(
θ
)
{{I}_{LN}}\frac{1}{Z}=-{{\omega }_{r}}C\left( {{V}_{C0N}}-{{V}_{N}} \right)\sin \left( \theta \right)+{{I}_{L0N}}\cos \left( \theta \right) \\
ILNZ1=−ωrC(VC0N−VN)sin(θ)+IL0Ncos(θ)
⇒
V
C
N
−
V
N
=
(
V
C
0
N
−
V
N
)
cos
(
θ
)
+
I
L
0
N
sin
(
θ
)
\Rightarrow {{V}_{CN}}-{{V}_{N}}=\left( {{V}_{C0N}}-{{V}_{N}} \right)\cos \left( \theta \right)+{{I}_{L0N}}\sin \left( \theta \right) \\
⇒VCN−VN=(VC0N−VN)cos(θ)+IL0Nsin(θ)
⇒
I
L
N
=
−
(
V
C
0
N
−
V
N
)
sin
(
θ
)
+
I
L
0
N
cos
(
θ
)
\Rightarrow {{I}_{LN}}=-\left( {{V}_{C0N}}-{{V}_{N}} \right)\sin \left( \theta \right)+{{I}_{L0N}}\cos \left( \theta \right) \\
⇒ILN=−(VC0N−VN)sin(θ)+IL0Ncos(θ)
两式左右两边平方相加可得
(
V
C
N
−
V
N
)
2
+
(
I
L
N
)
2
=
(
(
V
C
0
N
−
V
N
)
cos
(
θ
)
+
I
L
0
N
sin
(
θ
)
)
2
+
(
−
(
V
C
0
N
−
V
N
)
sin
(
θ
)
+
I
L
0
N
cos
(
θ
)
)
2
=
(
V
C
0
N
−
V
N
)
2
cos
2
(
θ
)
+
I
L
0
N
2
sin
2
(
θ
)
+
2
(
V
C
0
N
−
V
N
)
I
L
0
N
sin
(
θ
)
cos
(
θ
)
+
(
V
C
0
N
−
V
N
)
2
sin
2
(
θ
)
+
I
L
0
N
2
cos
2
(
θ
)
−
2
(
V
C
0
N
−
V
N
)
I
L
0
N
sin
(
θ
)
cos
(
θ
)
=
(
V
C
0
N
−
V
N
)
2
+
I
L
0
N
2
≡
R
2
\begin{aligned} &{{\left( {{V}_{CN}}-{{V}_{N}} \right)}^{2}}+{{\left( {{I}_{LN}} \right)}^{2}} \\ & ={{\left( \left( {{V}_{C0N}}-{{V}_{N}} \right)\cos \left( \theta \right)+{{I}_{L0N}}\sin \left( \theta \right) \right)}^{2}}+{{\left( -\left( {{V}_{C0N}}-{{V}_{N}} \right)\sin \left( \theta \right)+{{I}_{L0N}}\cos \left( \theta \right) \right)}^{2}} \\ &={{\left( {{V}_{C0N}}-{{V}_{N}} \right)}^{2}}{{\cos }^{2}}\left( \theta \right)+{{I}_{L0N}}^{2}{{\sin }^{2}}\left( \theta \right)+2\left( {{V}_{C0N}}-{{V}_{N}} \right){{I}_{L0N}}\sin \left( \theta \right)\cos \left( \theta \right) \\ &+{{\left( {{V}_{C0N}}-{{V}_{N}} \right)}^{2}}{{\sin }^{2}}\left( \theta \right)+{{I}_{L0N}}^{2}{{\cos }^{2}}\left( \theta \right)-2\left( {{V}_{C0N}}-{{V}_{N}} \right){{I}_{L0N}}\sin \left( \theta \right)\cos \left( \theta \right) \\ &={{\left( {{V}_{C0N}}-{{V}_{N}} \right)}^{2}}+{{I}_{L0N}}^{2} \\ &\equiv {{R}^{2}} \\ \end{aligned}
(VCN−VN)2+(ILN)2=((VC0N−VN)cos(θ)+IL0Nsin(θ))2+(−(VC0N−VN)sin(θ)+IL0Ncos(θ))2=(VC0N−VN)2cos2(θ)+IL0N2sin2(θ)+2(VC0N−VN)IL0Nsin(θ)cos(θ)+(VC0N−VN)2sin2(θ)+IL0N2cos2(θ)−2(VC0N−VN)IL0Nsin(θ)cos(θ)=(VC0N−VN)2+IL0N2≡R2
观察可得,归一化之后电容电压的归一化值的平方加上电感电流归一化值的平方等于定值R的平方。即知道了初值和激励的大小,可以绘制出电感电流和电容电压的状态轨迹。

其中,初始状态的角度为
θ
0
{{\theta }_{0}}
θ0
θ
0
=
−
arctan
(
V
C
0
N
−
V
N
I
L
0
N
)
θ
=
−
arctan
(
V
C
N
−
V
N
I
L
N
)
+
θ
0
\begin{aligned} & {{\theta }_{0}}=-\arctan \left( \frac{{{V}_{C0N}}-{{V}_{N}}}{{{I}_{L0N}}} \right) \\ & \theta =-\arctan \left( \frac{{{V}_{CN}}-{{V}_{N}}}{{{I}_{LN}}} \right)+{{\theta }_{0}} \\ \end{aligned}
θ0=−arctan(IL0NVC0N−VN)θ=−arctan(ILNVCN−VN)+θ0
二、总结和未来工作
综上所述,LC谐振腔的状态轨迹图即可求出。值得注意的是,前文中也提到了,激励 V N V_N VN是随着移相角的变化而变化的,因此实际LLC或者SR-DAB的谐振状态的轨迹并不是一个完整的圆,而是由半径和圆心都不同的原组成的。
- 未来将会着重探讨傅里叶级数的非正弦电压波形的数学解析表达式。
- 可能会对LLC或者SRDAB 全周期的状态轨迹展开探讨。
往期链接:
二阶线性微分方程求解
正弦激励响应公式推导
串联谐振型DAB的LC等效电路的归一化推导
电力电子论文中方波和类方波的傅里叶级数形式推导
1241

被折叠的 条评论
为什么被折叠?



