防雷设计作业

研究生课程考核试卷

（适用于课程论文、提交报告）

科目：变电站直击雷防雷设计教师:杨庆

姓名：第二组__________学号:

专业：电气工程与自动化类别：（学术、专业J）

上课时间：2015年5月至2015年6月

考生成绩:

卷面成绩平时成绩课程综合成绩

阅卷评语:

阅卷教师（签名）_______________________________

重庆大学研究生院制

目录

题目：..............................................................1

课程作业2.....................................................................................................................2

1ATP仿真.........................................................2

1.1ATP仿真简介................................................2

L2ATP仿真过程.................................................3

1.3仿真结果分析................................................4

2贝杰隆编程法解题..................................................7

2.1贝杰隆法原理................................................7

2.1.1贝杰隆法概述..........................................7

2.1.2各元件等值计算电路....................................8

2.1.3本题贝杰隆等效电阻....................................9

2.2贝杰隆法程序仿真结果.......................................13

2.2.1仿真参数设置.........................................14

2.2.2贝杰隆法仿真结果与ATP仿真结果对比..................14

3作业总结.........................................................16

附件1：贝杰隆法程序：.............................................17

变电站直击雷防护设计课程作业第二组

题目:

右图为某变电站220kV系统简化接

线图.

设雷电源为-20kA。

设备入口电容（微法）

隔离开关DS6E-5

电流互感器CT6E-4

断路器DL5E-4

电压互感器PT5E-3

主变T12E-3

求:

主变、进线入口、母线中间和变压器侧断路器上的过电压？（贝杰隆法编程、

ATP)

组号：第二组

组员：

张智

刘鑫

王伍静

林翔

万佳仑

周洪宇

韦冬洪

1

变电站直击雷防护设计课程作业第二组

课程作业2

下图为某变电站220kV系统简化接线图

雷电源为-20kA。

设备入口电容（微法）

隔离开关DS6E-5

电流互感器CT6E-4

断路器DL5E-4

电压互感器PT5E-3

主变T12E-3

图1.1变电站220kV系统简化接线图以及各设备入口电容

求：

主变、进线入口、母线中间和变压器侧断路器上的过电压？（贝杰隆法编程、

ATP）同时采用ATP研究避雷器安装在哪个节点能够有效保护变电站内的设备。

1ATP仿真

1.1ATP仿真简介

EMTP是用于电力系统电磁暂态分析的仿真软件。EMTP是Electro-Magnetic

TransientProgram（电磁暂态程序）的首字母缩写。为了对高压直流输电系统仿

真，程序中增加了模拟二极管和晶闸管等开关器件的能力，像SPICE程序一样，

现在有几种EMTP版本以用于个人计算机，如MicroTran、ATP等。所有版本的

程序都具有BPA（美国邦纳维尔电力局，BonnevillePowerAdministration）的

EMTP原版的大部分功能。EMTP是一个不断发展的软件，拥有其不断发展的大

量资源，因此成为美国电力系统和电子电力仿真方面的工作标准。其中

ATP-EMTP（AlternativeTransientsProgram）是较为广泛使用的一个版本。

目前，EMTP程序能够实现的主要功能有：

①计算网络的稳态解，也可把此稳态解再作为暂态计算的初始条件。

②计算网络中交直流电源频率从/■，”加到7mg以步长4/变化时各个频率下的

稳态解，以求得网络频率特性。

③所计算的网络暂态过程包括：各种电气设备和外来因素引起的暂态过程;

由各种故障引起的故障暂态及由系统各种元件参数的相互影响引起的谐振现象;

控制系统和一次系统相互影响的暂态过程。

④此外，还可根据需要将计算中得到的各个暂态量在某一时间段进行傅里叶

级数分解,得到各次谐波分量。

基于以上功能，目前ATP主要应用于以下仿真：

2

变电站直击雷防护设计课程作业第二组

①电力系统电磁暂态计算：电力系统暂态过电压分析，暂态保护装置的综合

选择，高压并联电抗器的选择，氧化锌避雷器选择等；

②电力系统的谐振过电压计算：计算由于磁饱和元件造成的铁磁谐振问题；

③机电暂态计算:汽轮发电机的轴系扭振问题。发电机组的次同步振荡问题。

本次变电站防雷设计主要利用变电站侵入波影响的研究，以及对于利用

Matlab编译贝杰龙算法的一个验证。

1.2ATP仿真过程

首先根据220kV变电站设备布局以及接线在ATP画出基本原理图如下:

图1.2ATP仿真接线图

各个器件的具体参数设定如下:

①雷电源：电流源为标准雷电流：L2/50RS,幅值为-20kA,雷电通道波阻抗

400Q,在图中表示为并联于电流源的400Q电阻；（电源用Heidler型电源）

图1.3雷电流、线路参数设定

②线路参数设定：根据常见的220kV变电站布局，除直接连接与母线上面

的设备外，其他各个设备之间的距离为3米，母线上各个节点的距离为10m,线

3

变电站直击雷防护设计课程作业第二组

路采用分布参数，用波阻抗来等效模拟，波阻抗同样设为400,传播速度为光速；

③隔离开关DS、电流互感器CT、断路器DL、电压互感器PT、主变T1按照

题目所给的参数设定。

④避雷器：本题采用单相线性避雷器，分别加在如图所示入口处、母线处、

断路器处以及主变处测量电压的节点上，探究避雷器的安装位置对于降低雷电流

侵入波幅值的影响。避雷器的伏安特性曲线设置如下：

图1.4避雷器参数设定

⑤仿真时间设定：仿真时间设为1ms,由于侵入波在最短路路径上3nl的传

播时间为10ns,所以仿真步长设定需要比改时间短，为2ns。

1.3仿真结果分析

当不采用避雷器的时候，入口处、母线处、断路器处以及主变处过电压波形

如下图所示，由于传播速度为光速，而传播距离较短，因此各处的波形只有微小

的时间差，但是他们的幅值都达到了6.5MV,而220kV的变压器耐雷电流侵入

波幅值为1MV左右，此时对于变电站无异于会造成严重的破坏，损坏变压器等

主要设备。因此，加入避雷器是非常有必要的。

图1.5未加避雷器测量节点的过电压波形

（进线入口过电压V：/;母线过电压V：Mu；变压器侧断路器过电压V：DL；主变过电压V：T1）

4

变电站直击雷防护设计课程作业第二组

当避雷器加在入口处时，入口处、母线处、断路器处以及主变处各个过

电压波形分别如下图所示：

(a)入口处过电压波形(b)母线处过电压波形

(c)断路器处过电压波形(d)主变处过电波形

图1.6入口处加装避雷器后各个测量节点过电压波形

此时，入口处过电压降低为-530kV,母线处、断路器处以及主变处过电压的

过电压波形会产生幅值较大，持续时间较长的振荡，振荡幅值分别达到了825kV,

825kV以及920kV,持续时间为130ns,然后衰减为-530kV左右。

当避雷器加在母线处时，入口处、母线处、断路器处以及主变处各个过电压

波形分别如下图所示：

(a)入口处过电压波形(b)母线处过电压波形

变电站直击雷防护设计课程作业第二组

（C）断路器处过电压波形（d）主变处过电波形

图1.7母线处加装避雷器后各个测量节点过电压波形

此时，母线处过电压降低为-530kV,入口、断路器处以及主变处过电压的

过电压波形会产生较小的振荡,振荡幅值分别达到了750kV,690kV以及875kV,

但是很快衰减为-530kV左右，持续时间为25ns左右。

当避雷器加在断路器处时，入口处、母线处、断路器处以及主变处各个过电

压波形分别如下图所示：

图1.8断路器处加装避雷器后各个测量节点过电压波形

此时，断路器处过电压降低为-530kV,入口处、母线处以及主变处过电压的

过电压波形会产生较小的振荡,振荡幅值分别达到了690kV,690kV以及710kV,

但是很快衰减为-530kV左右，持续时间为25ns左右。

当避雷器加在主变处时，入口处、母线处、断路器处以及主变处各个过电压

波形分别如下图所示：

6

变电站直击雷防护设计课程作业第二组

（a）入口处过电压波形（b）母线处过电压

（c）断路器处过电压波形（d）主变处过电波形

图1.9主变处加装避雷器后各个测量节点过电压波形

止匕时，主变处过电压降低为-530kV,入口、断路器处以及主变处过电压的过

电压波形会产生幅值较大，持续时间较短的振荡，振荡幅值分别达到了825kV,

860kV以及875kV,但是很快衰减为-530kV左右，持续时间为25ns左右。

因此，从上述结果我们可以得出以下结论：

①安装了避雷器后，不管避雷器安装在什么位置，各个主要节点上面的过电

压均能降低10倍左右；

②避雷器安装位置在不同位置，对不同设备的保护程度不同，离某一设备电

气距离短，则产生的过电压幅值越低，振荡越短，保护效果越明显；

③通过对前面的波形分析比较，避雷器安装在断路器处除了能使该点处过电

压有效降低外，还可以使入口处、母线处以及主变处过电压的过电压波形会产生

幅值最小、持续时间最短的振荡。因此，将避雷器安装在断路器处由最好的效果。

2贝杰隆编程法解题

2.1贝杰隆法原理

2.1.1贝杰隆法概述

贝杰隆法是把求解分布参数线路波过程的特性线法和求解集中参数电路暂

态过程的梯形法二者结合起来，形成的一种数值计算方法。它需要首先把分布参

数线路和集中参数储能元件（L、C）等值成为集中参数的电阻性网络，然后应

7

变电站直击雷防护设计课程作业第二组

用求解电阻网络的通用方法，计算实际电路的波过程。其有如下两个特点：

①整个分布参数线路的等值计算电路中只包括集中参数电阻和等值电流源,

属于集中参数电路。其中阻值等于线路波阻抗或由本身参数和延迟时间决定，等

值电流源由线路两端点上的电压和电流在过去的历史记录中计算得出。

②在单根无损线路等效模型中，计算电路中线路两侧节点k和m是独立分

开的，拓扑上没有直接联系给求解带来方便。其中k与m两端点间相互的电磁

联系是通过反映历史记录的等值电流源来实现的

2.1.2各元件等值计算电路

①单根均匀无损线的贝杰隆等值计算电路

设单根均匀无损线，长度为/,传播时间片波阻抗为Z,始端k半末端

m电压、电流分别以(t)、M,a(t)、九(t)和力”(t)。计算式为无损线如图2.1所示

_7)=_Um(t一。)Uk(t__2T)

Z

<2z(2.1)

/=(」_?)=--Uk(t——I—

Iz

-T)=—Uk(t-工)_2T)一_2T)

<zZ(2.2)

I(t-r)--

Lkz

②电感的贝杰隆等值计算电路

如图2.2(a)所示线性电感3已知时刻流经电感的电流和节点电位为

ifan(t-At)和赊(t-At)、wm(t-At),&是电感L的等值计算电阻，//为时间步长，/i(t-At)

为电杆等值计算电流源。

匕»('_4)=/工(/一24)+弁[以(Z-

c(2.3)

2ri

IL(^—A^)=IL(t—2A^)H---[〃%(，一△%)—〃冽(1—A[)]

R

、L

8

变电站直击雷防护设计课程作业第二组

lL(t-At)

(b)

图2.2电感的等值计算电路

③电容的等值计算电路

如图2.3(a)所示电容元件C,其等值计算电阻和等值计算电流源为:

R=—(2.4)

c2C

2「,

Ic(t—At)=Q-2Al)-----\ukQ-Az)Az)j(2.5)

Rc

Ic(t-At)

ikm

图2.3电容的等值计算电路

④电阻的等值计算电路

电阻元件R不能储能，其暂态过程与历史记录无关，其电压电流关系如下:

(2.6)

K

ikmRimkikmimk

—m—m

UmUm

(a)(b)

图2.4电阻的等值计算电路

2.1.3本题贝杰隆等效电阻

9

变电站直击雷防护设计课程作业第二组

以原仿真图结构如图2.5为基础得到贝杰隆等效电路，以图2.5中的1点及2

点间的电路为例，其变换了贝杰隆等值计算电路如图2.6所示。

DSCTDLCTDST1

图2.5变电站系统仿真图

图2.6仿真图中1点及2点间对应贝杰隆等值计算电路图

因而可获得1点及2点构成的等值电路节点电压方程如式（2.7）所示。

(2.7)

根据变电站系统仿真图，得到全等值计算网络，列出节点电压方程

YU=1(2.8)

可获得网络中任意节点的过电压值，15个节点构成的节点导纳矩阵如下:

y="为与LuE4%%?、兀后几几］(2.9)

节点导纳矩阵为:

Y=y^E(2.10)

10

变电站直击雷防护设计课程作业第二组

其中15个节点导纳值如下所示:

221212121

Y,=--Y2=-+一；Y3=-+——=—+--=­H--------------

ZZR.Z凡Z此.Z凡，

工DS工CT工DL工C7

21v3v21v11v2I

='钎；U=7'，Y&=y~R~，Yg=y~R~，YX0=7~R~（2.H）

Ziv.Z/i\./iv.Ziv.

（DS工ZXS、PT、DS

，21,21〃21，21，11

X1=3-；乂2=3+~D-%3=3q-片4=3~~；乂5=—+

ZKZK，ZzvZZVZ

（CT、DLXCT'-DS凡

其对应的节点电流源列向量为:

品2（2.12）

15

其中根据等值计算电路图，得到节点电流源列向量值如下：

=/（）（，）一，12（，一。）

12=一42—7）—乙3Q一7）—心口§（1-X）

13=一八2«-工）-八4«-工）-Lee,«-八力）

，4=_Z）_45Q_7）_/4C"Q_A'）

15=_/54（%一乙）一46Q—。）一4”Q—A%）

16=_，65—，67"_，）_/6cmI"A'）

[7=一，76（，一了）一，78（，一7）-/7JOQ—

4=-，870-口-789（"7）-，8CDS0-At）

19=_，98«_「）_195Q_A0（2.13）

Ao=一，10,7，-7）一—7）—，ioc*（,一AQ

Al=一乙1,10Q—7）—LU2Q—7）—LlCcrQ一A。

'12=一42J1“一7）一42,13“一7）一LlC"«-A'）

/13=一/13,12Q一。）一乙3,14Q一？）一'13"（'一A'）

44二一，14,13（t-T）-/l415（t-T）-/14Cos（t-At）

，15=一/15,14«一Z）一九电（/一A。

电路图中各元件电流方程为:

11

变电站直击雷防护设计课程作业第二组

2

12cme加)-12cm('2Af)

(DS

2

晨/加)一晨/2A/)03(/AO

RS

tt

14CDs△')—14CDs2加)、t)

11RCDL

2

。5«A/)

KQT

2

16CDs(t、t)_I6CDS(t2AO。6(，A，)

(DS

2

18C*八。-Ac/2At)。8«A。(2.14)

KCDS

2

[95(/A。—kep'2Ar)U9QAr)

J

2

1IOCDS"△')—1IOCDSC2加)UJ△/)

«CDS

2

加)一LiCcT。2A，)UiT(tAO

R”

2

I

儿入。一A。=~ncDLG-2A0-——。12«一加)

2

113cB(，-")=-113%。—2加)-Ui3(tA?)

R「

CCT

2

，14CDS("加)=Tgs。-2Af)-・U"A/)

R「

CD5

2

J

15CT]Q_')=_130nq_2、t)-——C/15(r-AO

12

变电站直击雷防护设计课程作业第二组

2.2贝杰隆法程序仿真结果

根据贝杰隆法原理，在matlab中编写程序并进行仿真，程序基本流程图如

图2.7所示，程序详见附件。

图2.7贝杰隆法程序流程图

雷电流双指数波形

/a

aw

tzn

。

0.20.40.60.8

Time/s

图2.8雷电流双指数波形

13

变电站直击雷防护设计课程作业第二组

2.2.1仿真参数设置

①全系统中雷电流双指数侵入波设置如图2.8所示：

1对)=1dmyt7”⑵⑸

左=1.0409,0=1.625X1()4,&=2.456X106,〃=-20kA即雷电流波幅值设为-20kA。

②波传播时间及仿真步长设置：

3米长传输段波传播时间设为10ns,13米长传播段波传播时间设为43.3ns仿真

步长设置为10/3ns

③各设备电容值及线路波阻抗设置：

各设备入口电容值设置与题中一致，隔离开关Rds=6e-11F,电流互感器6e-10F,

断路器5e-10F,电压互感器5e-9F,主变2e-9F。线路波阻抗设置为z=400C。

2.2.2贝杰隆法仿真结果与ATP仿真结果对比

①主变过电压波形

8

x10主变电压波形

/A

g)e

o-

>-4

-5E

0.000.020.040.060.080.10[ms]0.12

(fileNoname2.pl4;x-vart)v:T1

(a)Matlab仿真结果(b)ATP仿真结果

图2.9主变过电压波形仿真结果

②进线入口过电压波形

6

x10入口电压液形

-6

0.000.020.040.060.08。)。［ms］。)2

(a)Matlab仿真结窠(fileNoname2.pl4;x-vart)v:1

(b)ATP仿真结果

图2.10入口过电压波形仿真结果

14

变电站直击雷防护设计课程作业第二组

③母线过电压波形

xio6母线电压波形

之

3

6

£

>。

'7o0.20.40.60.8-7-

0.000.020.040.060.080-10[ms]012

Time/s

(fileNoname2.pl4;x-vart)v:MU

(a)Matlab仿真结果(b)ATP仿真结果

2.11母线过电压波形仿真结果

④变压器侧断路器过电压波形

xio6断路器电压波形八

0|------------1------------->-------------!-------------1-------------!-------------10-

i：：；：[MV]

A

/

se>

_l

>o

0.000.020.040.060.08[ms]

(fileNoname2.pl4;x-vart)v:DL

(a)Matlab仿真结果(b)ATP仿真结果

2.12变压器侧断路器过电压波形仿真结果

根据仿真结果，使用两种仿真方法获得--的各部分设备值变化情况基本相同。

使用matlab仿真结果：主变处过电压峰为-6.77MV,入口、断路器处以及主变处

过电压分别为-6.73MV,-6.73MV,-6.79MV。在同样条件下，通过ATP软件仿

真的结果：主变处过电压峰值-6.65MV,入口处以及母线过电压、变压器短路侧

断路过电压期间仿真结果分别为-6.56MV,-6.55MV,-6.55MV。两种仿真方式求

得的过电压波形变化趋势及大小基本一致，总体而言，贝杰隆法编程获取的过电

压幅值相对较低，但两种方式仿真结果比较下，过电压峰值大小差距在3.5%以

内，结果有较好的一致性。

15

变电站直击雷防护设计课程作业第二组

3作业总结

通过ATP软件仿真与贝杰隆编程法获得变电站220kV系统受-20kA雷电流

影响下主变、进线入口、母线和变压器侧断路器上的过电压变化情况。

1、获得各设备过电压幅值变化波形，其中过电压峰值幅值大小为:

表3.1各位置过电压峰值幅值大小

方式

ATP仿真贝杰隆法编程

位置

主变-6.77MV-6.747MV

进线入口-6.73MV-6.746MV

母线中间-6.73MV-6.747MV

变压器侧断路器-6.79MV-6.750MV

各处位置过电压仿真达到电压峰值时间：

方式

ATP仿真贝杰隆法编程

位置

主变13,351(is13,340|is

进线入口13.376(is13.367(is

母线中间13.363(is13.350(is

变压器侧断路器13.299|is13.286gs

2、通过使用ATP软件对加避雷器后的各处过电压进行仿真，获得在不

同部位加避雷器后各处的过电压波形，发现避雷器离某一设备电气距离短，

该设备上的过电压幅值越低，振荡越短，保护效果越明显；过对前面的波形

16

变电站直击雷防护设计课程作业第二组

分析比较，避雷器安装在断路器处除了能使该点处过电压有效降低外，还可

以使入口处、母线处以及主变处过电压的过电压波形会产生幅值最小、持续

时间最短的振荡。因而选择将避雷器安装在断路器处。

附件1:贝杰隆法程序:

%贝杰隆算法仿真雷电流

clearall

clc

%globalz,tl,t2,Th,I0;0%定义全局变量z（波阻抗）,tl（3m传播时

间）,t2（13m传播时间）,Th（仿真步长）

tl=10e-9;%10ns

t2=13/3*tl;%波传播时间，tl为3m,t2为13m

传播时间

Th=tl/3;%设置仿真步长10/3ns

w=pi/(1.2*10e-6);

Im=-20e3;%侵入雷电流波幅值-20kA

C=[6e-ll,6e-10,5e-10,5e-9,2e-9];%各设备电容值

R=Th/2./C;%R=[Rds,Rct,Rdl,Rpt,Rtl]

z=400;%线路波阻抗z=400Q

Yl=[2/z,2/z+l/R(l),2/z+l/R(2),2/z+l/R(3),2/z+l/R(2),2/z+l/R(l),3/z,2/z+l/R(l

),l/z+l/R(4),2/z+l/R(l),2/z+l/R(2),2/z+l/R(3),2/z+l/R(2),2/z+l/R(l),l/z+l/R(5)];%

电导率值

Y=diag(Yl);%节点导纳矩阵

U=zeros(15,l);%定义电压向量

I=zeros(15,l);%定义电流向量

forn=l:26

I12(n)=0;

I21(n)=0;

I23(n)=0;

I32(n)=0;

I34(n)=0;

17

变电站直击雷防护设计课程作业第二组

I43(n)=0;

I45(n)=0;

I54(n)=0;

I56(n)=0;

I65(n)=0;

I67(n)=0;

I76(n)=0;

I78(n)=0;

I87(n)=0;

I89(n)=0;

I98(n)=0;

I107(n)=0;

I710(n)=0;

I1011(n)=0;

I1110(n)=0;

I1112(n)=0;

I1211(n)=0;

I1213(n)=0;

I1312(n)=0;

I1314(n)=0;

I1413(n)=0;

I1415(n)=0;

I1514(n)=0;

I2ds(n)=0;

I3ct(n)=0;

I4dl(n)=0;

I5ct(n)=0;

I6ds(n)=0;

I8ds(n)=0;

I9pt(n)=0;

I10ds(n)=0;

Illct(n)=O;

I12dl(n)=0;

I13ct(n)=0;

I14ds(n)=0;

I15tl(n)=0;

Ul(n)=O;

U2(n)=0;

U3(n)=0;

U4(n)=0;

U5(n)=0;

U6(n)=0;

U7(n)=0;

U8(n)=0;

18

变电站直击雷防护设计课程作业第二组

U9(n)=0;

U10(n)=0;

Ull(n)=O;

U12(n)=0;

U13(n)=0;

U14(n)=0;

U15(n>0;%当n<26时，所有参数初始值全部为0

end

k=1.0409;

a=1.625e4;

b=2.456e6;%双指数雷电流波形系数

m=3;%3m传播时间

p=13;%13m传播时间

forn=27:36000%仿真时间120us

aa=n;

aa

T(n)=(n-27)*Th;%对系统进行计时矩阵，最长时间120uso

I0(n)=k*Im*(exp(-a*(n-27)*Th)-exp(-b*(n-27)*Th));%注入双指数雷电流

波形

%电流系数

Il2(n-m)=-2/z*U2(n-m)-I21(n-2*m);

12l(n-m)=-2/z*Ul(n-m)-Il2(n-2*m);

I23(n-m)=-2/z*U3(n-m)-I32(n-2*m);

I32(n-m)=-2/z*U2(n-m)-I23(n-2*m);

I34(n-m)=-2/z*U4(n-m)-I43(n-2*m);

143(n-m)=-2/z*U3(n-m)-I34(n-2*m);

I45(n-m)=-2/z*U5(n-m)-I54(n-2*m);

I54(n-m)=-2/z*U4(n-m)-I45(n-2*m);

I56(n-m)=-2/z*U6(n-m)-I65(n-2*m);

I65(n-m)=-2/z*U5(n-m)-I56(n-2*m);

167(n-m)=-2/z*U7(n-m)-I76(n-2*m);

I76(n-m)=-2/z*U6(n-m)-I67(n-2*m);

I78(n-p)=2/z*U8(n-p)-I87(n-2*p);

I87(n-p)=-2/z*U7(n-p)-I78(n-2*p);

I89(n-m)=-2/z*U9(n-m)-I98(n-2*m);

I98(n-m)=-2/z*U8(n-m)-I89(n-2*m);

1710(n-p)=-2/z*U10(n-p)-Il07(n-2*p);

1107(n-p)=-2/z*U7(n-p)-I710(n-2*p);

11011(n-m)=-2/z*Ull(n-m)-Illl0(n-2*m);

11110(n-m)=-2/z*U10(n-m)-I1011(n-2*m);

11112(n-m)=-2/z*U12(n-m)-I1211(n-2*m);

H211(n-m)=-2/z*Ull(n-m)-Illl2(n-2*m);

I1213(n-m)=-2/z*U13(n-m)-I1312(n-2*m);

11312(n-m)=-2/z*U12(n-m)-I1213(n-2*m);

19

变电站直击雷防护设计课程作业第二组

11314(n-m)=-2/z*U14(n-m)-11413(n-2*m);

I1413(n-m)=-2/z*U13(n-m)-I1314(n-2*m);

11415(n-m)=-2/z*U15(n-m)-I1514(n-2*m);

11514(n-m)=-2/z*U14(n-m)-11415(n-2*m);

%元件电流参数方程

I2ds(n-l)=-I2ds(n-2)-2/R(l)*U2(n-l);

I3ct(n-l)=-I3ct(n-2)-2/R(2)*U3(n-1);

I4dl(n-l)=-I4dl(n-2)-2/R(3)*U4(n-l);

I5ct(n-l)=-I5ct(n-2)-2/R(2)*U5(n-l);

I6ds(n-l)=-I6ds(n-2)-2/R(l)*U6(n-l)