双向DC-DC变换器广泛应用于需要双向功率传递的电能量转换设备中, 如微电网、可再生能源系统和电动汽车驱动系统等, 可作为不同电压等级的直流母线之间的接口变换器, 在能源管理系统中起着至关重要的作用。因此双向DC-DC变换器近些年来已成为一个热门的研究方向。
本文采用了一种带耦合电感的交错并联式双向DC-DC(bidirectional interleaved boost with coupled inductors, BIBCI) 变换器, 其前级为低压侧, 是交错并联式的Boost结构; 后级为高压侧, 是全桥结构, 且前后级通过耦合电感实现能量传输, 该拓扑可用于中高功率能量双向传输的电气设备上。BIBCI变换器虽然比著名的双有源全桥(dual active bridge, DAB)变换器结构稍微复杂一些, 但它拥有以下几个优点[1-4]: ①低压侧的交错并联结构使得电流在连续模式(CCM)运行; ②在2倍开关频率下可以减小电流纹波; ③通过耦合电感实现前后级电气隔离; ④能实现升压比灵活变化; ⑤所有开关管的零电压开关(ZVS)运行; ⑥确保变换器在宽功率范围内的高效率运行。
在实际应用中, 通常需要设计合适的控制策略, 来调节双向DC-DC变换器的输出电压和电流, 实现系统的控制目标。但对于设计一个性能良好的闭环控制器而言, 首先需要给出变换器的大信号和小信号的数学模型, 研究其在某一个稳态工作点的动态特性, 这是设计闭环控制器不可或缺的基础。
与DAB变换器拓扑类似, BIBCI变换器采用高频耦合电感在前后级之间实现电气隔离。由于这个基本性质, 这种结构很难通过传统的状态空间平均法来建立变换器的模型。为了解决这个问题, 有些学者采用了其他一些建模方法如: ①降阶平均模型, 忽略变压器电流动态特性, 仅考虑一个开关周期内的平均功率潮流[5-7]; ②广义平均模型, 使用傅里叶级数来表示变压器电流的状态变量, 可以较准确地预测动态性能[8]; ③离散时间模型, 采用离散状态方程表示每个开关状态区间内变压器电流动态[9-11]。
目前, 文献[1-2]比较有代表性, 基于PWM+相移调制方式, 建立了BIBCI变换器的降阶平均模型和小信号模型。然而, 在上述模型中, 都忽略了耦合电感的磁芯损耗和电路拓扑中的大量寄生参数, 而这并不能满足大带宽下对变换器性能的分析。
本文提出了一种适用于连续时间模式的BIBCI变换器的改进型降阶平均模型。所提出的模型考虑了电路中的许多主要功率损耗(如耦合电感的铜损、开关管的导通损耗、其他电感和电容等效串联电阻(ESR)以及磁芯损耗等)。并基于此种模型, 对BIBCI变换器进行大信号和小信号分析。
1 BIBCI变换器基本的BIBCI拓扑是由低压侧的交错并联Boost变换器和高压侧的全桥结构组成(如图 1所示), 前后侧通过耦合电感相连, 耦合电感通过绕组的漏电感提供能量传输的通道并实现电气隔离。值得注意的是, 图中的电感L不是一个额外的组件, 它是在高压侧等效的总漏电感, 由2个耦合电感的串联组成。详细电路如图 1所示。
前后级的功率传输由施加在等效漏电感L上的电压控制, 该电压取决于低压侧输入端电压通过耦合电感反映在高压侧的等效电压vA, 与高压侧的电压vB之差。其中vA的电压幅值为VA, 通过调节前级交错并联Boost环节的占空比D来控制; 而功率传输则通过调节BIBCI前后级之间的相移ϕ来控制, 两者结合即可实现单侧PWM+移相调制。这里, vA为三电平电压波形, vB为方波电压波形。根据D与ϕ的关系, 可以定义变换器的3个不同工作区域如下:
区域1:0 < ϕ < |D-1/2|
区域2:|D-1/2| < ϕ < (1-|D-1/2|)
区域3:(1-|D-1/2|) < ϕ < 1
其中: ϕ为相移; α为移相角(α=ϕπ); T为开关周期的一半。需要注意的是, 当ϕ=0.5时, 此时的传输功率最大[1], 所以ϕ的运行范围为[-0.5, 0.5]。而D的运行范围是[0, 1], 因此, 在实际工作中, BIBCI变换器将只会运行在区域1和区域2中ϕ < 1/2的部分。以区域2为例, 开关状态、电流和电压的波形如图 2所示, 其中ir和if分别表示低压侧和高压侧耦合电感的电流。
实际中, BIBCI拓扑的耦合电感的内阻和开关管导通损耗可以在高压侧等效为一个电阻Req
(1) |
式中, n为耦合电感的匝数比(n=np/ns)。此外, 电感Leq不是附加组件, 而是由耦合电感反映在高压侧的等效漏感
(2) |
根据(1)式和(2)式, 包含Req和Leq的BIBCI变换器等效拓扑电路如图 3所示。
2 改进型降阶模型本文针对BIBCI变换器, 提出一种改进型的降阶状态平均模型。
对于状态平均模型, 由于BIBCI变换器前级是由2个Boost结构作为输入单元, 交错运行, 等效为电流的加倍效应。同时, 2个Boost的电感电流ima和imb在一个开关周期中对称。因此, 平均模型中低压侧电流im由(3)式给出
(3) |
根据图 3所示, 建立BIBCI变换器改进型降阶模型, 选取钳位电容电压, 高压侧直流电压, 低、高压侧寄生电容电压和低压侧电流分别作为状态变量, 其状态方程用(4)式来描述:
(4) |
式中, vCL, vH, vCCL, vCH分别是施加在电容C1, C2, CCL, CH两端的电压。
下面以工作区域2为例, 由(5)式给出分别得到图 2中的t1和t2时刻为
(5) |
为了进一步说明BIBCI变换器的动态性能, 通过改进型降阶建模方法构建的BIBCI改进型降阶模型如图 4所示。从图中可以看出, 中间部分为Boost变换器模型的典型结构用来描述钳位电容电压vCL, 而ir和if作为电流源分别表示为耦合电感的前后级电流。
需要注意的是, Ts为开关周期, 而T=Ts/2。由于电路中有Req和Leq, ir和if的实际电流波形, 是分段指数函数, 可以表示为
(6) |
(7) |
式中, 0, t1, t2时刻的电流值为It0, It1, It2, 分别由以下公式计算可得
(8) |
如果考虑磁芯磁化电感LM和磁芯损耗等效电阻RM。可以推导出一个开关周期(2T)内的平均电流值If为
(9) |
式中, iM是通过LM的电流, 其平均值在开关周期Ts内为零。
(10) |
将(10)式代入(9)式, 则If可简化为
(11) |
因此, 根据(6)~(11)式, 并且考虑磁芯磁化电感LM的磁芯损耗等效电阻RM, 在一个周期T内的电流平均值Ir和If的推导过程如表达式(12)~(13)所示。
(12) |
(13) |
综合上述分析过程, 将(5)~(13)式代入(4)式, 得到了一种改进型降阶建模方法下BIBCI变换器的数学模型。接下来, 考虑稳态工作点处的小信号扰动并进行线性化处理, 可以进一步得出小信号模型。
3 小信号模型闭环控制器的设计和稳定性分析需要得到小信号, 以及控制-输出传递函数, 它反映了变换器在稳态工作点对控制扰动的小信号动态响应。
假设控制变量D, ϕ和状态变量vCL, vH, vCCL, vCDC和im的稳态值加上各自叠加的小交流扰动, 即
(14) |
将(14)式代入(4)式中, 可推导出平均模型的小信号状态方程如(15)式所示。
(15) |
为了计算模型的小信号增益, 函数
(16) |
(17) |
经过上述一系列的计算, 其中对hr, kr, gr, jr, hf, kf, gf, jf的计算过程将被省略。基于上述表达式, 可以立即推导出BIBCI变换器的小信号模型, 其状态方程可表示为
(18) |
式中
(19) |
(20) |
(21) |
(22) |
根据以上公式推导, BIBCI变换器的改进降阶模型的控制-输出传递函数, 可以由(23)式计算得到, 即
(23) |
式中:s是拉普拉斯变换符号;I是5×5的单位矩阵。
4 控制器设计考虑到应用场合为微电网或车载电力系统, BIBCI变换器的后级可以连接电压源负载或者阻性负载, 如图 5所示, 根据不同类型的负载设计出相应的控制器。
当变换器的负载为电压源负载时, 采用双闭环控制器。从图 5中可以看出, 2个比例积分(PI)调节器被放置在2个平行的控制环路中: 由占空比驱动的钳位电容电压环和由移相角驱动的低压侧电流环。根据前文所提出改进型降阶模型, 设计相应的PI调节器。通过在2个闭环调节器上设置不同的带宽, 主动解耦这2个控制环路。其中, 钳位电容电压环的带宽较大, 设为1~2 kHz; 而低压侧电流环的带宽较小, 设为几百赫兹。
对于vCL环, 令vCL*=n·VH(设为100 V), 使前级的钳位电容电压和高压侧输出电压保持匹配, 采用的大带宽调节器, 控制速度快, 不受iDCL环的影响, 无论iDCL怎样波动, vCL始终保持不变, 为耦合电感电流提供稳定的前级电压。而iDCL环带宽小, 可以认为在vCL电压恒定的情况下, 控制电流双向流通, 最终实现2个环路在频域上的解耦。
当负载为电阻负载时, 如图 5所示, 采用三闭环控制器, 在低压侧电流环的基础上, 再加一个高压侧电压环作为外环, 其带宽设为最低。vH环PI调节器的设计是一项相对简单的任务, 由于前2个控制环路已经闭环, 变换器的动态特性就像一个可控电流源, 为高压侧直流母线供电。最终, 设计vH环的带宽为100 Hz。因此, 高压端口的开环动态特性与简单的电流馈电电容器(即一阶系统)非常接近。
5 实验验证 5.1 实验平台为了对所提出的改进型降阶模型和所设计的控制器进行验证, 搭建了如图 6所示的实验平台。
其中, 控制电路将采用美国国家仪器公司生产的NI sbRIO-9606型板卡, 其控制核心是Spartan-6 LX45型FPGA芯片。控制算法是在内置FPGA芯片上实现的, 该芯片通过NI LabVIEW软件进行编程和管理。表 1列出了BIBCI变换器的关键参数, 其他一些寄生参数如(24)式所示。
(24) |
参数 | 值 |
电源电压(低压侧)VL/V | 48 |
额定高压侧电压VH/V | 400 |
额定功率P/W | 1 500 |
额定钳位电容电压VCL/V | 100 |
开关频率fsw/kHz | 20 |
匝数比n | 0.25 |
boost电感Lm/μH | 35 |
钳位电容CCL/μF | 6.8 |
高压直流电容CDC/μF | 50 |
负载电阻RH/Ω | 108 |
以下对2种BIBCI变换器的模型进行分析, 并与实际实验的测量结果进行比较:
2) 本文提出的改进型降阶平均模型。
首先, 根据2种建模方法, 分别建立了数学模型, 并计算出小信号模型和控制-输出传递函数。然后进行频域分析, 根据2种模型得到的传递函数, 在MATLAB软件中, 分别绘制出相应的波特图, 再与通过原理样机实际测量得到的结果进行对比和验证。
在实验中, 当BIBCI原理样机运行在额定工作点时, 给占空比和相移上分别施加小信号扰动, 其扰动信号
(25) |
式中:Dm和Φm分别是占空比和相移的幅值, 均设定为工作点幅值的2%。而fd和fϕ分别是扰动的频率。
通过波特图对2种模型的频域性能进行分析和验证。首先, 从图 7a)中可以看出, 对于钳位电容电压与占空比之间的传递函数GvCL-D而言, 将波特图和实验测得结果进行比较, 2种模型的理论计算值在大部分频率范围内和实测结果都符合得比较好。
然而, 对于相移与低压侧电流之间的传递函数GIm-ϕ而言, 2种模型与实测结果的差异是显而易见的。如图 7b)所示, 绿线(改进型降阶模型)的幅度比红线(降阶平均模型)在低频段高了大约10 dB, 中频段高了20 dB左右。同时, 传统降阶平均模型的没有显示出高频段的共振峰。
综上所述, 通过改进型降阶模型得到的传递函数与实际测量的传递函数之间的匹配良好, 从而证明了理论分析的准确性和有效性。而且所提出的改进型降阶模型具有良好的高阶特性, 与传统的降阶平均模型相比, 能够更准确地反映出高频响应性能。
5.3 控制器实验验证如前所述, BIBCI变换器可以作为微电网的储能系统和高压直流母线之间的接口变换器, 实现能量双向流通。
首先验证电压源负载时所设计的双闭环控制器, 选用超级电容作为能量的发出/吸收装置。在BIBCI变换器低压侧连接一个超级电容, 是由Maxwell公司生产的, 型号为BMOD0130P056B03, 其工作电压为0~56 V, 容值为130 F。高压侧并联一个100 Ω电阻和一个400 V的外部电源。
vCL环的PI调节器参数为kpvCL=0.01, kivCL=3 000 (V·s)-1; iDCL环参数为kpIm=0.001 A-1, kiIm=3 000 (A·s)-1。图 8展示了一个调制周期内vA, vB和耦合电感电流iL的波形。图 9展示了双闭环控制器的实验测试时的波形, 将参考值vCL*和设i*DCL为100 V和-5 A, 此时高压侧外部电源对超级电容进行充电。
最后, 去掉高压侧的外部电源, 对基于vCL, iDC, vH的三闭环控制器进行电阻负载测试。图 10展示了突加负载测试的实验结果, 其中参考值vCL*, vDC*和vH*分别为100 V和400 V, 高压侧负载电阻RH从360 Ω突加到180 Ω。
6 结论本文针对BIBCI变换器, 在PWM+移相双自由度调制方式下, 提出了一种改进型降阶模型, 并进行了理论分析。与传统的降阶平均模型相比, 新模型考虑了耦合电感的功率损耗、磁芯损耗, 以及电路拓扑中电感、电容和开关管等元件的寄生参数。所建立的新模型比传统降阶模型更复杂, 也更精确完善, 能够准确反映电路中元器件寄生参数和功率损耗对变换器动态特性的影响。对改进型降阶模型进一步分析, 推导出小信号模型和控制-输出的传递函数。在此基础上, 根据电压源负载和电阻负载的不同, 分别设计了双闭环和三闭环的控制器结构。
实验表明, 新模型与传统降阶模型相比较, 更加符合实际测量的结果, 能够更好地反映宽频域范围内的变换器的特性, 从而证明了新模型的准确性。然后在原理样机上对所设计的控制器进行了稳态和瞬态条件下的实验测试, 结果验证了控制器的有效性。
[1] | SPIAZZI G, BUSO S, BIADENE D, et al. High efficiency battery charger for photovoltaic inverters[C]//2017 IEEE Southern Power Electronics Conference, Puerto Varas, 2017: 1-6 |
[2] | SPIAZZI G, BUSO S. Small-signal modeling of the interleaved boost with coupled inductors converter[C]//2013 Brazilian Power Electronics Conference, Gramado, 2013: 456-461 |
[3] | TONIOLO F, LIU Q, CALDOGNETTO T, et al. Digital current control for a bidirectional interleaved boost converter with coupled inductors[C]//15th Brazilian and 5th Southern Power Electronics Conferences, 2019 |
[4] | SPIAZZI G, BUSO S. Analysis of the interleaved isolated boost converter with coupled inductors[J]. IEEE Trans on Ind Electron, 2015, 62(7): 4481-4491. DOI:10.1109/TIE.2014.2362496 |
[5] | RODRIGUE A, SEBASTIAN J, LAMAR D. An overall study of a dual active bridge for bidirectional DC/DC conversion[C]//2010 IEEE Energy Conversion Congress and Exposition, 2010: 1129-1135 |
[6] | BAI H, Mi C, WANG C, et al. The dynamic model and hybrid phase-shift control of a dual active-bridge converter[C]//2008 34th Annual Conference of IEEE Industrial Electronics, Orlando, FL, USA, 2008: 2840-2845 |
[7] | BAI H, NIE Z, MI C. Experimental comparison of traditional phase-shift, dual-phase-shift, and model-based control of isolated bidirectional DC-DC converters[J]. IEEE Trans on Power Electronics, 2010, 25(6): 1444-1449. DOI:10.1109/TPEL.2009.2039648 |
[8] | QIN H, KIMBALL J. Generalized average modeling of dual active bridge DC-DC converter[J]. IEEE Trans on Power Electronics, 2012, 27(4): 2078-2084. DOI:10.1109/TPEL.2011.2165734 |
[9] | ZHAO C, ROUND S. KOLAR J, Full-order averaging modeling of zero-voltage-switching phase-shift bidirectional DC-DC converters[J]. IET Power Electronics, 2010, 3(4): 400-419. |
[10] | SHI L, LEI W, LI Z, et al. Bilinear discrete-time modeling and stability analysis of the digitally controlled dual active bridge converter[J]. IEEE Trans on Power Electronics, 2017, 32(11): 8787-8799. DOI:10.1109/TPEL.2016.2640659 |
[11] | KRISMER F, KOLAR. Bilinear discrete-time modeling and stability analysis of the digitally controlled dual active bridge converter[J]. IEEE Trans on Power Electronics, 2017, 32(11): 8787-8799. DOI:10.1109/TPEL.2016.2640659 |