贝杰龙法介绍

1.2.1Bergeron特征线计算理论Bergeron特征线计算方法就是利用线路波过程的特征线方程。经过一定的转换，把分布参数的线段等值为电阻性网络。再运用求解电阻性网络的通用方法来计算整个网绘的暂态过程，也可以说Bergeron特征线计算方法是把求解分布参数线路波过程的特性线法和求解集中参数电路暂态过程的梯行法两者结合起来而形成的数值计算方法。因此它首先需要把分布参数线路和集中参数储能元件（L、C）等值成为集中参数的电阻性网络来计算实际电路的波过程［2al-2a2］。在求解网络暂态过程时，从计算开始时刻t0起，把时间离散成一系列的时间间隔，一般采用等时间步长。每一步计算时，利用t时刻以前作为历史记录存好的历史状态，计算出这一时步的电压和电流状态，并记录该状态作为历史值，为下一时步计算作准备。在讨论运特特征线方法数值求解线路波过程之前。先简要介绍一下特征线方法的基本原理。分布参数线路上任何一点的对地电压和导线中的电流是距离x和时间u的函数，是电磁波沿线路传播的过程，若先考虑到线路单位长度的电阻R、电感L、电导G和电容C均0000为常数。和频率无关。则单导线线路上的波过程可以用以下的偏微分方程来描述：dudxdudx0 0dtTOC\o"1-5"\h\z=Gu+C竺

dx 0 0dt若略去损耗；则可以有如下的无损线的偏微分方程：dx 0dtdi du——C—

dx 0dt对上述方程组进行求导换算可以得到如下二阶偏微分方程组d2uld2u

dx2V2dt2d2i_1d2i

dx2V2dt2其中1V二一.^LC~00为流动波沿线的传播速度。对无损架空线路等于光速．即电磁波在真空中的传播速度以上单根无损线波动方程的电压和电流解可以写成以下形式：u(x,t)二u(x—Vt)+u(x+Vt)fbi(x,t)二i(x—Vt)+i(x+Vt)fb其中uf和「分别代表以速度为V沿着x正方向传播的前行电压波和电流波，而你和ib分别代表以速度为V沿着x反方向传播的前行电压波和电流波。前行电压波和前行电流波之间,以及反行电压波与反行电流波之间是通过波阻抗相联系的：i(x-Vt)二u(x-Vt)/Zf f (2)i(x+vt)——u(x+vt)/Zb b■L将式(2)代入式(1)中就可以得到一下两个行波的特征方程：u(x,t)+Zi(x,t)—2u(x—vt)fu(x,t)—Zi(x,t)—2u(x+vt)b以上两个特征方程的物理意义可以描述如下：方程式中u(x,t)和i(x,t)分别表示在线路上x点在t时刻的电压和电流的瞬时值。对前行波来说，其物理意义为：因为线路均匀无损，所以电磁波沿线略向前传播时不发生畸变和衰减，当观察者沿x正方向以速度v和前行波一起运动(即x—Vt=常数)，则根据他所处的位置x在t时刻观察到的瞬时电压值u(x,t)和电流值i(x,t)所计算得到的u(x,t)+Zi(x,t)的值始终保持不变，等于两倍前行电压波的大小，这种情况从线路始端一直到末端都成立。对反行波来说，其物理意义为：若观察者沿x反方向以速度v和前行波一起运动，则他在线路上任一点x在t时刻观察到的u(x,t)—Zi(x,t)的值始终保持不变，等于两倍反行电压波的大小，这种情况从线路始端一直到末端都成立。根据以上特征方程及其物理概念可以推导出单相无损线波过程计算的等值电路及其相应的计算公式。单相无损线的波过程计算的等值电路如下所示：设一线路两端节点分别为k和m，其电压，电流为Uk(t)、\m(t)、Um(t)、Imk(t)，如图1所示：ikm 1・丁N imkuk Lim图1根据特征方程的物理概念，若观察者在t—T时刻从节点k出发，则t时刻达到节点m，由前行特性方程：

U(t—T)+Zi(t_T)=U(t)—Zi(t)TOC\o"1-5"\h\zk km m mk11可得:i(t)= U(t)— U(t—T)—i(t—T)可得:\o"CurrentDocument"mk zm zk km同时设：—.I(t—T)=——U(t—T)—i(t—T) (3)m zk km可简化为：1i(t)=u(t)+1(t—T) (4)mk zm mTOC\o"1-5"\h\z根据上式可以得到如图2所示的线路末端m在时刻t的等值计算电路，而公式(3)和(4)就是相应的计算公式，其中z是阻值等于线路波阻抗的电阻，I(t—t)是等值电流源,m他可以根据过去观察的记录，即线路始端在t—T时刻的电压U(t—T)和电流i(t—T)计算k km得到。imkimkHzUm同样，观察者可以随反行波从末瑞节点m运动到始端节点k，得到公式:I(t—T)=一-U(t—T)—i(t—T)k zm mk]⑸i(t)=—U(t)+1(t—T)km zk k由公式(5)可以得到图：综上，推广到整个线路，整个分布参数线路的等值计算电路中只包括集中参数电阻（阻值等于线路波阻抗）和等值电流源（