电气工程及其自动化AT牵引供电系统计算模板Word文档下载推荐.docx

1、尽管在长期的工程实践中，人们已经发展了许多简化方法来解决电气化铁道工程设计中的诸如变压器容量选择、电压损失计算、短路计算、防干扰计算等问题，在利用电子计算机完成牵引供电有关计算方面，人们也进行了持续的努力，但到目前为止，对牵引供电系统数学模型和电气参数的研究还不能说已经进行得很完善，特别是在利用计算机进行仿真计算方面，即使仅从满足工程设计需求角度，也还存在差距。同普通三相电力系统相关研究的水平相比，还没有达到它应该达到的深度和广度。1.3 国内外研究现状AT牵引供电系统作为一个特殊的高压输配电网络，牵引供电系统中涉及到的各式各样的问题一般总可以归为电路的稳态或电磁暂态问题。因此，从网络角度，研

2、究牵引供电系统的稳态模型和电磁暂态模型十分有必要。对三相电力系统的数学模型，经过人们长时间的大量研究，已经有广泛而深入的结果，这给建立牵引供电系统的数学模型奠定了很好的基础。实际上，电气化铁道研究者一直把三相电力系统中的有关原理、方法设法移植到牵引供电电力系统中来，便解决了大量工程实际问题。人们也针对牵引供电系统的特点，研究出了一些独特的简化计算方法，如，采用广义对称分量法进行复线AT牵引网的短路计算。在利用计算机完成牵引供电的分析计算方面，主要集中在防干扰计算、面向运行图的供电计算和谐波及负序分布计算等方向。因AT牵引网结构复杂，对其供电特性、防干扰特性的计算只有借助计算机才能得到比较准确的

3、结果，日本学者在开发AT供电技术过程中，首次把多导体传输线模型引入到牵引网的计算中，1979年日本国铁提交给国际电报电话咨询委员会（CCITT）的集中反映了这一成果。Talukdar和Koo最早在The analysis of electrified ground transportation networks中把潮流计算引入到地铁直流牵引供电系统的分析中，Cai等人则把潮流计算同直流牵引网的多导体传输线模型结合起来，随后又把这一思路运用到AT牵引供电系统2。到目前为止，还没有见到把这种建模方法应用到非AT牵引供电方式交流牵引网的文献。对当前应用较多的带负馈线的直接供电方式牵引网，尽管简化手算

4、模型能满足一般的计算要求，但要详细掌握其供电特性和防干扰特性，仍需采用类似AT牵引网的电算模型。在我国还没有见到把潮流计算引入到牵引供电系统运行仿真中的公开文献。2004年西南交通大学同铁道第四勘察设计院编制一个采用潮流算法的供电计算程序，但对牵引网采用了简化模型。在研究牵引负荷对电力系统的负序、谐波影响中，已有不少文献采用了潮流算法，这些研究中一般不对牵引供电系统详细建模，牵引变电所只作为电网的节点出现，本质上仍是三相电力系统的潮流计算。2 AT牵引供电系统计算本章节主要介绍AT牵引供电系统进行仿真计算的参数和进行初步的仿真实验。根据给定的参数，进行AT牵引供电系统的阻抗矩阵、电容矩阵、AT

5、矩阵和三绕组电源变压器矩阵仿真，且对结果进行分析。对仿真过程的总结，绘制出了相应的计算步骤图，有助于能够更清晰的了解仿真的全部过程。2.1 AT牵引供电系统计算参数图2.1 AT牵引供电系统网络图2.2 牵引网导线空间分布 计算模型依据图2.1和图2.2，基于 AT 网的主要参数如下：（1）计算导线的确定 考虑复线AT供电方式，上行和下行各有8根导线，共有2x8=16根导线，分别为: 上行导线（1.承力索 MW1；2.接触线 CW1；3.钢轨 RA1；4.钢轨 RB1；5.正馈线 PF1；6.保护线 PW1；7.通信线 T1；8.综合接地线 E1），下行导线（9.承力索 MW2；10.接触线

6、CW2；11.钢轨 RA2；12.钢轨 RB2；13.正馈线 PF2；14.保护线 PW2；15.通信线 T2；16.综合接地线 E2），牵引网的空间分布位置如图3.2所示。（2）供电臂参数表3-1 供电臂参数牵引变电所武清DK80+300分区所天津DK115+420供电臂长度（km）35.1AT所间距20.414.7电源空载电压58kV（50HZ）牵引变压器阻抗j9.60（单相VV结线，容量为31.5MVA）电力系统阻抗j1.08（系统短路容量为2800MVA）AT漏抗0.1+j0.45大地导电率10-4 1/CM钢轨泄漏电阻10保护线接地方式接钢轨接地电阻2-5牵引网导线型号接触线J、承力

7、索CMgCu-120、THJ-120正馈线F、保护线PWLGJ-240、LGJ-120钢轨R60kg/m综合接地线ETJ-70表3-2 接触悬挂类型导线型号计算半径（mm）等效半径（mm）电阻（/km）接触线 CWMgCu-1205.94.20.146承力索 MWTHJ-1207.005.310.158正馈线 PFLGJ-1859.59.030.163钢轨 RP60109.112.790.135+j0.135保护线 PWLGJ-1207.607.220.255综合地线E5.354.0550.28通信线T（3）牵引变电所采用V/X接线单相变压器 容量： 40000kVA 额定电压：220kv/2

8、27.5kv 短路阻抗（基准容量25000kVA）：高压-T绕组 10.22% 高压-F绕组 10.20% 2.2 MATLAB仿真软件简介 在20世纪70年代的中期，Cleve Moler博士和他的同事在美国国家基金会的帮助下，开发了Linpack和Eispack的Fortran语言子程序库，其中Linpack适用于线性方程求解的子程序库，Eispack适用于特征值求解的子程序库。在70年代后期，Cleve Moler博士设计了Linpack和Eispack的接口程序，并将程序取名为MATLAB。这就是现在广泛使用的MATLAB的起源。MATLAB是由美国MathWorks公司开发的大型软件

9、。在MATLAB软件中，包括了两大部分：数学计算和工程仿真。其数学计算部分提供了强大的矩阵处理和绘图功能。在工程仿真方面，MATLAB提供的软件支持几乎遍布各个领域，并且不断加以完善。MATLAB软件在电力系统建模和仿真中的应用，主要涉及到MATLAB软件中的电力系统仿真模块（SimpowerSystems）和控制系统工具箱（Control System Toolbox）。此外随着电力系统的发展，模糊逻辑工具箱（Fuzzy Logic Toolbox）、信号处理工具箱（Signal Processing Toolbox）、小波工具箱（Wavelet Toolbox）和神经网络工具箱（Neura

10、l Network Toolbox）等在电力系统中也有着广泛的应用7。本文仿真所使用的是MATLAB 7.0.1仿真软件，由其数学计算部分提供了强大的矩阵处理和绘图功能，进行本文的AT牵引供电系统的计算与分析的研究。2.3 计算步骤图 图2.3 计算步骤图2.4 计算中几个重要的函数程序函数文件是另一种形式的M文件，每一个函数文件都定义一个函数。事实上，MATLAB提供的标准函数大部分都是由函数文件定义的。在MATLAB编程中，一个完整的程序往往是由几个部分组成的，因此，把握和调试好重要的程序是非常必要的。（1）程序的基本参数 %1.承力索1；2.接触线1；3.钢轨1；4.钢轨2；5.正馈线1

11、；6.保护线1；7.通信线1；8.综合接地线1；9.承力索2；10.接触线2；3.钢轨3；4.钢轨4；5.正馈线2；6.保护线2；7.通信线2；8.综合接地线2%原始参数clear;clc;format long;Pxy=0.001*0,7500;0,6300;-755,1000;755,1000;-4400,8500;-3600,8000;-9400,6500;-4400,500;5000,7500;5000,6300;4245,1000;5755,1000;9400,8500;8600,8000;14400,6500;9400,500; %牵引网导线空间分布坐标 单位mR=0.001*7.

12、00;5.9;109.1;9.5;7.6;5.35;7.00;5.35; %16根导线计算半径 单位mRd=0.001*5.31;4.2;12.79;9.03;7.22;4.055;5.31;4.055; % 16根导线等效半径 单位mr=0.158;0.146;0.135+0.135i;0.163;0.255;0.28;0.158;0.28; %16根导线单位电阻 单位/kmf=50; %AT牵引供电系统供电频率 单位HZsigma=1e-4; %大地导电率 单位1/.cmepsilon=8.854188e-12; %空间介电系数Kat=1; %AT的变比n1/n2Zg=0.1+0.45i;

13、 %AT漏抗 单位Rg=10; %钢轨泄漏电阻 单位.M%三绕组变压器的参数Kpower=220/27.5; %三绕组变压器的变比Z1_T=0.1978592; %三绕组变压器T绕组短路后从一次侧看到的阻抗 Z1_F=0.197472; %三绕组变压器F绕组短路后从一次侧看到的阻抗Zs=1.08i; %电力系统阻抗TN=16; %导线数cs=105; %横向连接元件的参数FD=1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1.4,1,1,1,1,1,1,1,1,1,1,1,1,1,1.7; %AT牵引供电系统各分段长度构成的数组n=length（FD）; %求数组

14、元素的个数AT_Location=zeros（1,n+1）; %AT分布位置 AT_Location（20）=1;AT_Location（n+1）=1;POWER_Location=zeros（1,n+1）; %三绕组变压器的位置 POWER_Location（1）=1;lj=zeros（1,n）; %机车位置lj（25）=0;Dlcw_pf=zeros（1,n+1）; %接触线与正馈线短路位置Dlcw_pf（25）=0;Dlcw_r1=zeros（1,n+1）; %接触线与钢轨短路位置Dlcw_r1（25）=0;Dlpf_r1=zeros（1,n+1）;Dlpf_r1（25）=1; %正馈线

15、与钢轨短路位置（2）串联阻抗矩阵函数的建立 串联阻抗矩阵的函数建立是根据第2章2.3.3节串联元件的节点阻抗矩阵中的（2-115）和（2-116）两式进行编写。 串联阻抗矩阵函数：%建立函数ZK来产生16x16阶串联阻抗矩阵functionZ0=ZK（f,sigma,Pxy,r,Rd）Dg=（0.2085/sqrt（f*sigma*1e-9）/100; %导线地回路等值深度 单位 mr2=pi*pi*f*1e-4; %大地交流电阻 单位/kmfor i=1:16 for j=1: if（i=j） Z0（i,j）=r（i）+r2+0.0029i*f*log（Dg/Rd（i）; else Z0（i

16、,j）=r2+0.0029i*f*log10（Dg/sqrt（Pxy（i,1）-Pxy（j,1）.2+（Pxy（i,2）-Pxy（j,2）.2）; endend调用和调试函数，在Workspace中可以得到1616阶单位阻抗数据图形：图3.4 阻抗三维图 由图2.4可知，自阻抗的值比互阻抗的值大，符合理论计算值。（3）电容导纳矩阵函数的建立 电容导纳矩阵函数的建立是根据第2章2.3.2节含有并联导线的多导体系统的电容矩阵中（2-109）和（2-112）两式建立。 电容导纳矩阵函数：%建立函数DR来产生16x16阶的电容矩阵functionC,YDR=DR（epsilon,Pxy,R,f） if

17、（i=j） P（i,j）= 1/（2*pi*epsilon）*log（2*Pxy（i,2）/R（i）; else P（i,j）=1/（2*pi*epsilon）*log（sqrt（Pxy（i,1）-Pxy（j,1）.2+（Pxy（i,2）-Pxy（j,2）.2）/sqrt（Pxy（i,1）-Pxy（j,1）.2+（Pxy（i,2）+Pxy（j,2）.2）;C=inv（P）;YDR=2i*pi*f*C;16阶单位电容的数据图形：图3.5 电容三维图 由图2.5可知，电容的自阻抗的值比互阻抗的值大，符合理论计算值。（4）AT导纳矩阵函数的建立本文在第2章2.3.1节中推算出了AT导纳矩阵，即（2-

18、89）式。根据该式，把相应的元素代入到建立的1616阶的矩阵中。AT导纳矩阵函数：%建立函数AT来产生16x16阶矩阵functionYAT=AT（Kat,Zg）Kat1=1/（1+1/Kat）; %AT中n1/（n1+n2）Kat2=1/（1+Kat）; %AT中n2/（n1+n2）Yat=zeros（16）;Yat（2,2）=Kat12;Yat（2,3）=-Kat1;Yat（2,5）=Kat1*Kat2;Yat（3,2）=-Kat1;Yat（3,3）=1;Yat（3,5）=-Kat2;Yat（5,2）=Kat1*Kat2;Yat（5,3）=-Kat2;Yat（5,5）=Kat22;YAT=

19、Yat/Zg;（5）V/X接线电源变压器矩阵函数建立本文在第2章2.3.1节中推算出了V/X接线变压器的导纳矩阵，即（2-102）式。根据该式和V/X接线变压器在AT牵引网络的接线形式，把（2-102）式中的相应元素代入到建立的16V/X接线电源变压器矩阵函数： %建立函数POWER来产生16x16阶的三绕组变压器导纳矩阵和外界注入电流矩阵functionYPOWER,GPOWER=POWER（Kpower,Z1_T,Z1_F,Zs,TN）ZT_1=Z1_T/（Kpower2）;ZF_1=Z1_F/（Kpower2）;ZTF_1=3*ZT_1+ZF_1-4*0.24;Z1_TF=ZTF_1/（

20、Kpower2）;Z1=0.5*（Z1_T+Z1_F-Z1_TF）;Z2=0.5*（ZT_1+ZTF_1-ZF_1）;Z3=0.5*（ZTF_1+ZF_1-ZT_1）;m0=（Z1+Zs）*Z2+（Z1+Zs）*Z3+Kpower2*Z2*Z3;mpower1=（Z1+Zs+Kpower2*Z2）/m0;mpower2=-（2*（Z1+Zs）+Kpower2*Z2）/m0;mpower3=（Z1+Zs）/m0;mpower4=（4*（Z1+Zs）+Kpower2*2*Z2）/m0;Ypower=zeros（TN）;Ypower（2,2）=mpower1;Ypower（2,3）=mpower2;

21、Ypower（2,5）=mpower3;Ypower（3,2）=mpower2;Ypower（3,3）=mpower4;Ypower（3,5）=mpower2;Ypower（5,2）=mpower3;Ypower（5,3）=mpower2;Ypower（5,5）=mpower1;YPOWER=Ypower;mpower5=-Kpower2*Z2*Kpower/m0;mpower6=-mpower5;Gpower=zeros（TN,1）;Gpower（2,1）=mpower5;Gpower（5,1）=mpower6;GPOWER=-220000*Gpower;（6）M矩阵的形成在广义四端口网络表

22、征的AT供电系统中有一个系数矩阵，即（2-71）式，在该式中矩阵里的元素也是一个矩阵，是矩阵嵌套矩阵，根据这一特点本文采用了MATLAB中的for循环和cell细胞来构造矩阵。M1=cell（n+1,n+1）;n+1 if i=j M1i,j=M0（:,:,i）; elseif j-i=1 M1i,j=D（:,i+1）; elseif i-j=1,j+1）; else M1i,j=zeros（16）;M=cell2mat（M1）; cell是MATLAB中的一种数据类型，用大括号定义，括号里可以是任意类型的数据或矩阵。2.5 小结本章简要的介绍了AT牵引供电系统给定计算参数、电力系统仿真软件M

23、ATLAB的作用和对广义四端口网络表征的AT供电系统中的重要部分进行初步的计算。 （1）绘制出了相应的计算步骤图，总体把握整个计算的过程，并根据计算步骤图，计算阻抗矩阵、电容矩阵和AT导纳矩阵等。 （2）理论分析和计算结果表明，阻抗矩阵和电容导纳矩阵的对角元素的数值比非对角元素大，即自阻抗比互阻抗大，自电容比互电容大9。因此，本文推导的阻抗和电容矩阵是正确的。3 AT牵引供电系统的牵引网特性与分析本章根据第2章所给的AT牵引供电系统的参数，将对牵引供电系统正常运行及各种故障条件下进行详细的计算，运行仿真程序，得出相应的仿真图形，进行分析。通过仿真实验得出的结果，验证了多导体传输线理论在AT牵引

24、供电系统计算分析上的可行性。3.1 计算无机车运行时电压、电流分布图3.1 简单的单线AT网络示意图 根据自耦变压器的电流关系，可知图中回路中有如下关系： （3-1）由上式可知，在理想状态下，无机车运行时，供电系统处于开路，电流为零。 图3.2 无机车运行时钢轨的电流曲线图图3.3 无机车运行时接触线的电流曲线图3.4 无机车运行时馈线的电流曲线 由图3.2，3.3，3.4可知，AT牵引供电系统空载运行时，钢轨、接触线和馈线上的电流都很小，几乎为零，跟理论值接近。3.2 计算机车正常运行时电压、电流分布 为了研究AT牵引供电系统中机车正常取流时的特性，本文考虑在20km和35km处分布有AT变

25、压器，25km处机车取流1000A的情况，计算钢轨电位分布（如图3.5），钢轨、接触线和正馈线电流分布（如图3.6，37，3.8）。图3.5 机车正常运行时的钢轨电位曲线 观察图3.5，可以知道，钢轨最高电压出现在机车取流25km处，中间20kmAT变压器处也有个小峰值。正常运行时钢轨电位小于60 V，满足规定要求。图3.6 机车正常运行时的钢轨电流曲线图3.7 机车正常运行时接触线的电流曲线图3.8 机车正常运行时馈线的电流曲线观察图3.6，3.7，3.8，可以发现取流点25km处到中间 20kmAT变压器处区段钢轨、接触线和正馈线电流很大。在25km点处，机车从牵引网上取流，电流从接触线上经机车流到25km点，此时流经钢轨的电流