High-precision Runge-Kutta method for transmission line equation
-
摘要: 提出了一种求解传输线方程的高精度龙格-库塔(RK)方法。此方法在空间上采取高阶泰勒展开,提高了对空间微分的近似精度,减少了数值色散所带来的误差。与传统的时域有限差分法(FDTD)方法相比,在每波长采样数相同时,RK方法的计算精度更高。同时,根据Taylor模型,对外界平面波激励源进行离散,成功利用RK方法对外部场激励传输线进行求解,扩大了龙格−库塔方法在求解传输线方程时的应用范围。通过编程对平面波辐照下无限大地平面上的单导体与双导体的算例分别应用FDTD方法与RK方法进行了计算,验证了RK方法的正确性。结果表明同等计算条件下RK方法的计算精度更高。Abstract: This paper presents a high-precision Runge-Kutta (RK) method for solving transmission line equations. This method adopts high-order Taylor expansion in space, which improves the approximation accuracy of spatial differentiation. Compared with the traditional finite element time-domain method, when the number of samples per wavelength is the same, RK method has higher precision. At the same time, according to the Taylor model, researchers use RK method to solve transmission line equation in the external field excitation. The correctness and high precision of the RK method are verified by numerical examples of our study.
-
Key words:
- Runge-Kutta /
- field-line coupling /
- transmission line /
- FDTD /
- Taylor expansion
-
为了实现电源系统结构的紧凑型和小型化,如何在保证功率与效率情况下减小电源的体积与重量成为了需要解决的问题。高频交流链接技术(HF AC-Link)与传统的DC-Link技术相比,减少了AC-DC-AC中的直流变换部分,省去了包括中间储能电容在内的大量元器件,从而使电源结构更加紧凑[1]。文献[2-6]对基于AC-Link技术的串联谐振充电电源进行了研究,对传统的AC-Link电路进行了改进,并使用状态平面分析法对电源的工作模式进行了分析,并进行仿真和实验验证。文献[7-8]对DC-Link技术的串联谐振的Buck-Boost电路的工作模式进行了分析,并引入状态平面分析法进行控制算法的分析。本文提出了一种基于AC-Link串联谐振的Buck-Boost电路拓扑结构,通过状态平面分析法,提出了一种控制算法,并进行仿真验证。仿真结果显示,控制算法能够实现Buck和Boost的功能,输入电流谐波含量低,并且可以实现变换器在较宽的范围内调节输出电压。
1. AC-Link串联谐振的Buck-Boost变换器结构
基于AC-Link串联谐振的Buck-Boost变换器拓扑结构如图 1所示, 该变换器拓扑结构由LC输入滤波器单元、IGBT开关矩阵单元、LC串联谐振单元、高频变压器和输出整流单元构成。LC输入滤波器单元为由滤波电容与滤波电感组成的二阶低通滤波器,用来滤除电路中的杂波,提高变换器系统的抗干扰能力。开关矩阵单元由12只IGBT组成,两个开关组成一组的双向开关构成全桥整流结构的6个桥臂。与之前提出的变换器拓扑结构整流输出单元采用4个不可控二极管作为开关器件不同的是,Buck-Boost变换器结构输出整流单元中,两个上桥臂采用不可控的二极管,下桥臂采用可控开关管IGBT。当变换器工作在Buck模式时,副边两个IGBT处于关断状态,能量可通过反并联在IGBT两端的二极管流动,此时整流结构与传统的整流结构一致。当变换器工作在Boost模式时,通过对下桥臂IGBT开关的控制,对副边变压器两端短路,原边谐振电感储存能量;当副边变压器两端导通时,谐振电感上储存的能量去补充输出电容的能量输出到负载,以实现Boost功能。
2. 状态平面图建立与分析
状态平面分析法适用于各类谐振变换器的分析, 它通过谐振变换器运行中分段线性的特点, 建立每一工作过程的特征方程, 构造每个工作过程的运行轨迹线, 并组合成闭合的稳态轨迹图, 通过图解分析其动态和暂态特性[9-10]。基于状态平面分析法可对AC-Link串联谐振的Buck-Boost变换器的工作过程进行分析,可以根据每个过程的等效电路图得到相应的轨迹方程,进而建立状态平面图。根据状态平面分析法,提出基于该变换器的控制算法。
2.1 Boost模式与Buck模式的判定
在基于AC-Link的变换器,采用在半个谐振周期内进行相间切换的方式对三相输入的电流进行控制。在串联谐振的半个谐振周期内,电网三相中的其中任一相连接在串联谐振电路一端,另外两相分别依次连接到串联谐振电路的另一端。对于串联谐振电路,由电网三相组成的2个激励电压,高线电压Vn与低线电压Vm对串联谐振电路进行激励。高线电压Vn与低线电压Vm可以通过三相输入电压Ua,Ub,Uc两两相互组合计算出。将工频输入电压的一个周期分为12个区间,每个区间的工作相为π/6。以第一个区间为例,三相线电压的大小关系为Ucb>Uab>Uca。在第一区间内选取电压值较大的两个线电压Ucb与Uab作为激励电压。将Ucb称为高线电压,记为Vn;将Uab称为低线电压,记为Vm[3]。
假设低线电压Vm激励时流入的电荷量为Q1,从高线电压Vn激励时流入的电荷量为Q2;可以根据流入的电荷量得到一个比值K。
K=Q1Q2 (1) 定义等效激励电压为Vin;在谐振电路的半个谐振周期内,由低线电压Vm和高线电压Vn提供的能量应等于等效激励电压Vin提供的能量。
Vin(Q1+Q2)=VnQ2+VmQ1 (2) 由式(1)~(2)可以得到等效激励电压的表达式为:
Vin=KVm1+K+Vn1+K (3) 根据等效激励电压Vin与输出电压Vo间的关系可以判定出变换器工作的模式。因此,当满足Vin>Vo/n时,变换器工作在Buck模式;当满足Vin < Vo/n时,变换器工作在Boost模式。其中n为高频变压器的匝数比。
2.2 Boost模式状态平面图的建立与分析
根据2.1中对Buck、Boost模式的判定方式。当Vin < Vo/n时,变换器工作在Boost模式。在Boost模式整个工作过程中采用6个过程的工作模式,具体如下:第一个工作过程采用正向高线电压Vn激励;第二个工作过程采用正向高线电压Vn激励;第三个工作过程采用正向低线电压Vm激励;第四个工作过程采用负向高线电压-Vn激励;第五个工作过程采用负向高线电压-Vn激励;第六个工作过程采用负向低线电压-Vm激励。在Boost模式下的正半周的3个工作过程等效电路图如图 2所示,负半周的3个工作过程与正半周相似。
选取谐振电感电流iLr和谐振电容电压vCr为状态变量,建立每个工作过程的状态方程。以第一个工作过程为例,可得状态方程如下
vCr=iLrCr (4) iLr=Vn−vCrLr (5) 解方程可得到对应的时域方程
vCr=Vn+ILr0Zosinωt+(VCr0−Vn)cosωt (6) iLr=ILr0cosωt−(VCr0−Vn)sinωtZo (7) 式中:ILr0为iLr起始时刻的值;VCr0为vCr起始时刻的值;ω为谐振角频率ω=1√LrCr;Zo为谐振阻抗,Zo=√LrCr。将式(6)和(7)整理可得第一个工作过程的状态轨迹方程
(vCr−Vn)2+i2LrZ2o=(VCr0−Vn)2+I2Lr0Z2o (8) 同样也可以根据上述的方法推导出第2个工作过程和第3个工作过程的状态轨迹方程
[vCr−(Vn−Vo/n)]2+i2LrZ2o=[vCr0−(Vn−Vo/n)]2+I2Lr0Z2o (9) [vCr−(Vm−Vo/n)]2+i2LrZ2o=[vCr0−(Vm−Vo/n)]2+I2Lr0Z2o (10) 根据上述的3个工作过程的状态轨迹方程可以得到稳态下正半周工作情况下的状态平面图如图 3所示。横轴坐标为谐振电容电压vCr,纵坐标为谐振电流iLr和谐振电路特征阻抗Zo的乘积。其负半周的状态平面图与正半周轨迹相同,方向相反。A1B1为第一工作过程的运行轨迹,B1B2为第二工作过程的运行轨迹,B2A2为第三工作过程的运行轨迹。B1为第一工作过程与第二工作过程切换点,B2为第二工作过程与第三工作过程切换点。可以根据状态平面图中各个参数的几何关系,得到每个工作过程的工作时间。
根据图 3的状态平面图中的几何关系可得
r21−(O1−V1)2=r22−(O2−V1)2 (11) r23−(V2−O3)2=r22−(V2−O2)2 (12) 式中:O1=Vn;O2=Vn−Von;O3=Vm−Von;r1=O1+vCr;r3=vCr−O3。
电网侧电荷分配的约束条件为
K=Q1Q2=vCr−V2V2+vCr (13) 根据r1,r3,K的关系,可以得到r2的表达式为
r2=√(O1−O2)(r23+O22−O23)+(1+K)(O2−O3)(r21+O22−O21)−2KvCr(O1−O2)(O2−O3)(O1−O3)+K(O2−O3) (14) 在Boost模式下,一个工作周期内谐振电容电压的平均值为0,在正半周第一工作过程的起始时刻(A1)与负半周第一工作过程的起始时刻(A2),谐振电容电压的绝对值相等,根据电容电压与电流关系推出其表达式为
vCr=PoTr4CrVin (15) 式中:Po为输出功率;Tr为谐振频率;Lr为谐振电感的值;Cr为谐振电感的值。
Tr=2π√LrCr (16) 由r1,r2,r3,V1和V2,可以得到每个工作过程所需要的工作时间
t1− Boost =arccos(O1−V1r1)ω (17) t2− Boost =arccos(V1−O2r2)−arccos(V2−O2r2)ω (18) t3− Boost =arccos(V2−O3r3)ω (19) 式中:ω=2π/Tr。
2.3 Buck模式状态平面图的建立与分析
当Vin>Vo/n时,变换器工作在Buck模式。Buck模式下由高线电压Vn、低线电压Vm和0电压,三个激励电压对串联谐振电路进行激励。具体如下:第一个工作过程采用高线电压Vn激励;第二个工作过程采用正向低线电压Vm激励;第三个工作过程采用0电压激励;第四个工作过程采用负向高线电压-Vn激励;第五个工作过程采用负向低线电压-Vm激励;第6个工作过程采用0电压激励。在Buck模式下的正半周的3个工作过程等效电路图如图 4所示。
采用与Boost模式相同的分析方法,根据3个工作过程等效电路图所推导的轨迹方程可以得到稳态工作情况下的状态平面图如图 5所示。
Buck模式与Boost模式相同,在正半周的起始时刻(A1)与负半周的起始时刻(A2),谐振电容上的电压绝对值相等,其表达式为
vCr=nPoTr4CrVo (20) 电网侧电荷分配的约束条件为
K=Q1Q2=V2−V1V1+vCr (21) 由图 5中状态平面图中的几何关系可得到,同样可求解出在Buck工作模式下每一个工作过程所需要的工作时间
t1− Buck =arccos(O1−V1r1)ω (22) t2− Buck =π−arccos(V2−O2r2)−arccos(O2−V1r2)ω (23) t3− Buck =arccos(V2−O3r3)ω (24) 3. 仿真验证
为了验证根据状态平面图分析法所提出的控制算法,通过之前对变换器在Buck与Boost模式下的工作状态的分析,使用Matlab/Simulink软件搭建了仿真模型并在不同的工作模式下进行了仿真分析。仿真参数设置如表 1所示。
表 1 仿真模型参数Table 1. Simulation model parametersthree-phase voltage/V resonant inductor/μH resonant capacitor/μF load/Ω transformer ratio output capacitance/μF 380 20 0.9 18.75 1:1.5 47 在三相输入电压为380 V的条件下,等效激励电压Vin的变化范围459~530 V。当设置输出电压Vo为600 V时,Vin>Vo/n成立,此时变换器工作在Buck模式。如图 6所示,变换器可以稳定的输出符合设定值的电压;输入电流谐波含量为3.13%,如图 7所示。当设置输出电压Vo为850 V时,此时变换器工作在Boost模式,与Buck模式一样可以得到稳定的输出符合设定值的电压;输入电流谐波含量为3.03%,如图 8和图 9所示。仿真结果显示,该控制算法,在Buck模式和Boost模式下,可以稳定的输出电压,具有较高的功率因数,输入电流谐波含量低。
4. 结论
本文对AC-Link串联谐振的Buck-Boost变换器拓扑结构通过状态平面图的建立和分析,提出一种控制算法。建立Matlab/Simulink建立仿真模型验证后,证明了控制算法的正确性。仿真显示,输出电压稳定,输入电流谐波含量低,可以实现变换器在较宽的范围内调节输出电压。
-
[1] 张希. 有损均匀传输线数值解的研究[D]. 重庆: 重庆大学, 2003: 42-59.Zhang Xi. Reasearch on lossy and uniform transmission line numerical solution[D]. Chongqing: Chongqing University, 2000: 21-30 [2] 张希, 刘宗行, 孙韬, 等. 传输线方程的一种数值解法[J]. 重庆大学学报(自然科学版), 2004, 27(2):116-119. (Zhang Xi, Liu Zongxing, Sun Tao, et al. A numerical method for transmission line equations[J]. Journal of Chongqing University(Natural Science Edition), 2004, 27(2): 116-119 [3] 杨阳. 龙格库塔法求模糊微分方程的数值解[D]. 哈尔滨: 哈尔滨工业大学, 2015: 25-35.Yang Yang. Numerical solution of fuzzy differential equations by Runge-Kutta method[D]. Harbin: Harbin Institute of Technology, 2015: 25-35 [4] 秦浩东, 吴志强, 张晏铭, 等. 龙格-库塔方法与差分法的比较[J]. 成都大学学报(自然科学版), 2014, 33(4):337-338. (Qin Haodong, Wu Zhiqiang, Zhang Yanming, et al. Comparison of Runge-Kutta method and difference method[J]. Journal of Chengdu University(Natural Science Edition), 2014, 33(4): 337-338 [5] 陈山. 求解波动方程的龙格-库塔型方法及其地震波传播模拟[D]. 北京: 清华大学, 2010: 8-10.Chen shan. The Runge-Kutta type method for solving wave equations and its simulation of seismic wave propagation[D]. Beijing: Tsinghua University, 2010: 8-10 [6] 梁华力, 富明慧. 一种改进的精细-龙格库塔法[J]. 中山大学学报(自然科学版), 2009, 48(5):2-5. (Liang Huali, Fu Minghui. An improved precise Runge-Kutta integration[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2009, 48(5): 2-5 [7] 关冶, 陆金甫. 数值分析基础[M]. 北京: 高等教育出版社, 2010.Guan Ye, Lu Jinpu. Fundamentals of numerical analysis[M]. Beijing: Higher Education Press, 2010 [8] Paul C R. 多导体传输线分析[M]. 3版. 北京: 中国电力出版社, 2013.Paul C R. Analysis of multiconductor transmission lines[M]. 3rd ed. Beijing: China Electric Power Press, 2013 [9] 葛德彪, 闫玉波. 电磁波时域有限差分法[M]. 3版. 西安: 西安电子科技大学出版社, 2011.Ge Debiao, Yan Yubo. Finite-difference time-domain method for electromagnetic waves[M]. 3rd ed. Xi′an: Xidian University Press, 2011 [10] 冀维林. 基于FDTD算法的多导体传输线电磁兼容的研究[D]. 北京: 北京邮电大学, 2010: 21-43.Ji Weilin. Study on electromagnetic compatibility of transmission line based on FDTD method[D]. Beijing: Beijing University of Posts and Telecommunications, 2010: 21-43 -