二阶RLC电路响应与参数的关系及工程应用

李志平

(南方电网科学研究院, 广东 广州 510663)

二阶RLC串联电路的暂态过程,在电力生产实际中比较常见。高压直流输电换流器在电流不连续工况下的阀触发导通为工频激励下二阶电路过程。在控制工程中,二阶系统有很重要的地位,这不仅是因为二阶系统在数学上容易分析,更是因为可以将高阶系统降为二阶系统进行处理,仍不失其运动过程的基本性质[1]。

二阶RLC电路暂态过程的响应由一个二阶线性常系数微分方程求解,电路理论及高等数学有成熟的解法[2-3]。自动控制原理在上述基础上引入一个阻尼系数进行讨论。求解微分方程可以得出特征根,因电路参数不同其根为3种情况,对应的响应结果是过阻尼的非振荡型和欠阻尼的振荡型及临界型。文献[4-8]就暂态响应结果和电路参数的关系进行了讨论,对不同的电路参数给出了仿真结果。本文通过引入一个振荡系数,简化了二阶电路暂态响应的计算公式,对非振荡型响应提出了二阶RLC电路的响应近似等同于一阶RC电路的响应条件,对振荡型响应提出了简明的判断振荡频率周期和衰减趋势方法。同时,列举一个高压直流输电的二阶电路的暂态过程算例[9],并利用MATLAB/Simulink进行仿真[10],验证了理论推导的正确性。

1 二阶RLC串联电路响应的简述

二阶串联RLC电路如图1所示。

图1 二阶串联RLC电路示意

图1中,L,R,C分别为电感、电阻和电容,US为直流激励电压,uC和i为响应电压和电流,uC(0)=UC0为零输入响应的初值。

二阶串联RLC电路的零输入响应的微分方程为

(1)

2 非振荡型RLC电路响应与参数的关系

对于电路的响应为非振荡型,其微分方程的通解为

uC=A1e-r1t+A2e-r2t

(2)

式中:A1,A2——根据初始条件确定的常数。

A1+A2=UC0

(3)

-r1A1-r2A2=0

(4)

式(3)×r2+式(4)得

(5)

式(3)×r1+式(4)得

(6)

由式(4)得

(7)

式中:k——通解两项系数的比值,仅是m值的函数,其随m值变化如表1所示。

表1 通解两项的比值随m值变化

当t=1/r2=τ2时刻,计算m为0.1,0.3,0.5,0.8时通解第1项和第2项的比值 -(uC1/uC2)=-(A1e-r1t/A2e-r2t)=ke×e-1/k,k值较大时用级数展开-(uC1/uC2)≈ke×(1-1/k)。振荡系数m=0.5时通解第1项和第2项的比值等于13.26,忽略通解第2项造成的结果误差(小于10%),τ2时刻后的结果误差更小。

再讨论响应通解第1项的衰减系数-r1和一阶RC电路的衰减系数-1/RC的关系。m值和电感L成正比,与电容C成反比,与电阻R的平方成反比,L减少而C和R不变时导致m值变化。响应通解第1项的衰减系数为

(8)

表2 RLC电路通解的衰减系数随m值变化情况

由表2可知,m=0.3时,K(0.3)=1.090,则r1和一阶RC电路的衰减系数接近,因为响应通解第二项可以忽略,所以二阶RLC电路的零输入响应近似等同于一阶RC电路的响应。直流激励下的全响应性质相同,因此可以得出类似结论。

3 振荡型RLC电路响应与参数的关系

对于振荡欠阻尼型,其微分方程通解为

uC=e-δt[A1cos(ωt)+A2sin(ωt)]

(9)

由零输入响应的初始条件,可得A1=UC0,A2ω-δA1=0,所以,

(10)

化为

(11)

而β=arctanω/δ,回路电流为

(12)

在上述结果的基础上,利用振荡系数m对响应与元件参数关系进行分析,发现振荡频率和周期、响应的电压幅值的衰减系数e-δt随系数m的变化规律。 假定L,C值不变,R值减少使m值从1.1增大到1 000对电路响应特性的影响如表3所示。

表3 m值随R改变的二阶RLC电路的响应特性

由表3可知,衰减系数ω/ωo反映一个振荡半周期时刻电压幅值的衰减程度;m值大于10后振荡角频率和无阻尼LC电路的振荡角频率比值接近1.0。

4 实例二阶RLC电路的计算

4.1 非振荡型RLC电路的响应特性计算

改变L使m=0.5,则L=mR2C/4=500 mH,k=-A1/A2=r2/r1=5.8,计算得A1=1.21UC0,A2=-0.21UC0,则-r1τ=-K(0.5)=-1.17,uC/uRC=(A1/UC0)e-r1τ+1=1.02。

上述计算表明,振荡系数m为0.3和0.5时,τ=RC时刻二阶响应和一阶响应相近。改变L使m=1,则L=mR2C/4=1 000 mH。改变L使m为零,则L为零,即二阶RLC电路变为一阶RC电路。利用MATLAB/Simulink进行仿真,得到的非振荡二阶RLC电路仿真结果如图2所示。

图2 非振荡二阶RLC电路仿真结果

4.2 振荡型RLC电路的响应特性计算

m=100时,第1个振荡半周期时刻ωt=π时,uC=-UC0e-0.314=-73.0 V。零输入响应电流最大值β=1.47 rad,则-δtp=-0.148,所以imax=0.086 3 A。

利用MATLAB/Simulink进行仿真,得到的暂态电压、电流波形图如图3和图4所示。由图3和图4可知,振荡周期、暂态电流最大值出现时刻都与计算相符。

图3 改变R值导致m改变时响应电压仿真结果

图4 改变R值导致m改变时响应电流仿真结果

4.3 高压直流工程RLC电路的响应特性计算

换流器电流不连续工况是高压直流输电的一种,出现在启动投入运行过程的换流器单阀间歇性触发导通的暂态过程。系统交流电源经触发导通后的单阀接外部电路对换流器其他阀组件的相关RC回路电容充电。该暂态过程可等效为直流激励下的二阶RLC电路响应。对此进行计算可获得不导通阀的偏置电压。因为电路含有单向导通的晶闸管元件,其过程仅延续半个振荡周期,电流过零后阀截止即结束。

响应电流最大值β=1.26 rad,则-δtp=-0.403,所以imax=22.6 A。

本例的振荡系数m接近于10。由表3可知,当m值处于10附近时,响应结果较适中,即有一定的电流峰值和电压过冲值,电流脉冲也得到一定程度的拓宽。这是高压直流输电系统在电流不连续工况下所期望的状态。实际的激励源是下降沿的工频正弦电压波,在振荡半周内变化不大时可以近似为直流激励源。高压直流输电系统的电流不连续工况还有其他一些工作模式,也可以近似为二阶RLC电路响应进行研究。

5 结 论

(1) 振荡系数m值小于等于0.3时,二阶串联RLC电路的零输入响应近似等同于一阶RC电路的零输入响应。

(2) 利用振荡系数m可以得到振荡型暂态响应的简化计算公式。据此可以更直观地判断暂态响应结果与电路参数的关系。

(3) 利用MATLAB/Simulink进行了仿真,验证了本文所提理论的正确性。

(4) 用所推导的方法对一个典型的500 kV高压直流工程在单阀触发工况下的实例二阶RLC电路进行了分析和计算,暂态响应结果符合期望的状态。