串联谐振型DAB的LC等效电路的归一化推导


前言

在串联谐振型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}} \\ λ=jLC 1jω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(tt0))+C2sin(ωr(tt0))+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(tt0))+ωrCC2cos(ωr(tt0))
微分方程解的形式见之前的博客(二阶线性微分方程求解链接: 二阶线性微分方程求解
电容电压和电感电流的初值定义为
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(tt0))+(ωrCiL(t0))sin(ωr(tt0))+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(tt0))+iL(t0)cos(ωr(tt0))

归一化准备:

归一化首先需要锚定一个基准值,这里设置电压的基准值为原边的输入电压 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(tt0)=θ
代入微分方程的通解可得
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=(VC0NVN)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(VC0NVN)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=(VC0NVN)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(VC0NVN)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=(VC0NVN)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(VC0NVN)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) \\ VCNVN=(VC0NVN)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=(VC0NVN)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} (VCNVN)2+(ILN)2=((VC0NVN)cos(θ)+IL0Nsin(θ))2+((VC0NVN)sin(θ)+IL0Ncos(θ))2=(VC0NVN)2cos2(θ)+IL0N2sin2(θ)+2(VC0NVN)IL0Nsin(θ)cos(θ)+(VC0NVN)2sin2(θ)+IL0N2cos2(θ)2(VC0NVN)IL0Nsin(θ)cos(θ)=(VC0NVN)2+IL0N2R2
观察可得,归一化之后电容电压的归一化值的平方加上电感电流归一化值的平方等于定值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(IL0NVC0NVN)θ=arctan(ILNVCNVN)+θ0

二、总结和未来工作

综上所述,LC谐振腔的状态轨迹图即可求出。值得注意的是,前文中也提到了,激励 V N V_N VN是随着移相角的变化而变化的,因此实际LLC或者SR-DAB的谐振状态的轨迹并不是一个完整的圆,而是由半径和圆心都不同的原组成的。

评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

电力小圆

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值