1、广东工业大学华立学院课 程 设 计（论文）课程名称 电力系统课程设计 题目名称 不对称故障分析与计算及其程序设计 学生学部（系） 电气工程系 专业班级 电气工程及其自动化（4）班 学 号 12030804035 学生姓名 覃 烽 指导教师 罗洪霞 2011 年 6 月 12 日电力系统分析课程设计（论文）1目 录摘 要 1关键词 1前 言 .11.电力系统短路故障的基本知识 .21.1 短路故障的概述 22 对称分量法在不对称短路计算中的应用 .22.1 不对称三相量的分解 32.2 对称分量法在不对称短路计算中的应用 43 简单不对称短路的分析与计算 .43.1 单相（a 相）接地短路 73
2、.3 两相（b 相和 c 相）短路 接地 .74 简单不对称短路的分析与计算计算机计算程序法 .84.1 简单故障的计算程序原理 94.2 网络节点方程的形成 .105 电力系统不对称短路计算实例 115.1 单相接地短路和两相短路不对称故障分析与计算 .115.2 两种计算方法的对比 .18结语 .19参考资料 .19附录：不对称短路电流计算程序 .20电力系统分析课程设计（论文）摘 要在电力系统运行过程中，时常会发生故障，且大多是短路故障。短路通常分为三相短路、单相接地短路、两相短路和两相接地短路。其中三相短路为对称短路，后三者为不对称短路。电力运行经验指出单相接地短路占大多数，因此分析与
3、计算不对称短路具有非常重要意义。分析计算不对称短路方法很多，目前实际最常用的方法是对称分量法。而以对称分量法为核心的计算方法又可有解析法和计算机程序算法等，本论文的主要工作即介绍这两种计算方法。解析法，是将微分方程代数化、暂态分析稳态化、不对称转化为对称并叠加完成不对称故障的分析与计算。计算机程序算法是在形成三个序网的节点导纳矩阵后，对其应用高斯消去法求得故障端点等值阻抗，根据故障类型选用相应公式计算各序电流、电压，进而合成三相电流、电压。关键词: 单相接地短路，两相短路，两相接地短路，对称分量法，节点导纳矩阵前 言电力系统分析是一门介绍电力系统稳态运行分析、故障分析和暂态过程分析的课程。电力
4、系统分析的基础为电力系统潮流计算、短路故障计算和稳定计算。在电力系统运行过程中，时常会发生故障，其中大多数是短路故障（简称短路） 。所谓短路，是指电力系统正常运行情况以外的相与相之间或相与地（或中性线）之间的连接。产生短路的主要原因是电气设备载流部分的相间绝缘或相对地绝缘损坏。此外运行人员在短路检修后未拆除地线就加电压等误操作也会引起短路故障。短路问题是电力技术方面的基本问题之一。在发电厂、变电站以及整个电力系统的设计和运行工作中，都必须事先进行短路计算，以此作为合理选择电气接线、确定限制短路电流措施等的重要依据。为此计算短路时各种运行参量（电流、电压等）是非常必要的。短路通常分为三相短路、单
5、相接地短路、两相短路和两相接地短路。其中三相短路为对称短路，后三者为不对称短路。电力运行经验指出，单相接地短路占大多数。因此分析与计算不对称短路具有非常重要意义。1.电力系统短路故障的基本知识1.1 短路故障的概述在电力系统运行过程中，时常发生故障，其中大多数是短路故障。所谓短路：是指电力系统正常运行情况以外的相与相之间或相与地（或中性线）之间的连接。除中性点外，相与相或相与地之间都是绝缘的。电力系统短路可分为三相短路，单相接地短路。两相短路和两相接地短路等。三相短路的三相回路依旧是对称的，故称为不对称短路。其他的几种短路的三相回路均不对称，故称为不对称短路。电力系统运行经念表明，单相短路电力
6、系统分析课程设计（论文）占大多数，上述短路均是指在同一地点短路，实际上也可能在不同地点同时发生短路，例如两相在不同地点接地短路。依照短路发生的地点和持续时间不同，它的后果可能使用户的供电情况部分地或全部地发生故障。当在有由多发电厂组成的电力系统发生端来了时，其后果更为严重，由于短路造成电网电压的大幅度下降，可能导致并行运行的发电机失去同步，或者导致电网枢纽点电压崩溃，所有这些可能引起电力系统瓦解而造成大面积的停电事故，这是最危险的后果。产生短路的原因很多主要有如下几个方面（1）原件损坏，例如绝缘材料的自然老化，设计，安装及维护不良所带来的设备缺陷发展成短路，(2)气象条件恶化，例如雷电造成的闪
7、络放电或避雷针动作，架空线路由于大风或导线覆冰引起电杆倒塌等。(3)违规操作，例如运行人员带负荷拉刀闸(4)其他，例如挖沟损伤电缆。在电力系统和电气设备的设计和运行中，短路计算是解决一系列技术问题所不可缺少的基本计算，比如在选择发电厂和电力系统的主接线时为啦比较不同方案接图，进行电力系统暂态稳定计算，研究短路对用户的影响。合理配置各种继电保护和自动装置并正确整定其参数都必须进行短路的计算和分析.2.对称分量法及其在计算中的应用2.1 不对称三相量的分解人们在长期的实践中发现，在三相电路中，任意一组不对称的三相相量（电压或电流） ，可以分解为三组三相对称的相量分量，式（2-1） 。在线性电路中，
8、可以用叠加原理对这三组对称分量按照三相电路去解，然后将其结果叠加起来。就是不对称三相电路的解答，这个方法就叫做对称分量法。设 ， ， 为三相系统中任意一组不对称的三相量，可以分解为三组对称的三序分量如aFbc下：120120aabbcccFF（2.1）三相序分量如图 2-1 所示：)1(aI)(b )2(aI)2(bI)(c )0()(cbaI正序分量 负序分量 零序分量电力系统分析课程设计（论文）图 2-1 三序分量正序分量： 三相的正序分量大小相等，彼此相位相差 ，与系统正常对称运行(1)(1).abcF 120对称运行方式下的相序相同，达到最大值 ，在电机内部产生正转磁场，这就是正序分a
9、bc量。此正序分量为一平衡的三相系统，因此有： (1)()(1).abcF负序分量： 三相的负序分量大小相等，彼此相位相差 ，与系统正常对称运(2)(2).abc 行对称运行方式下的相序相反，达到最大值 ，在电机内部产生反转磁场，这就是正序c分量。此正序分量为一平衡的三相系统，因此有： (2)()(2)0.abc零序分量： 三相的零序分量大小相等，相位相同，三相的零序分量同时达到最(0)(0).abcF大值，在电机内部产生漏磁，其合成磁场为零，这就是零序分量。如果以 A 相为基准相，各序分量有如下关系：正序分量： (1)a2(1)(1)baF2(1)(1)(1)cbaF负序分量： 2F 2零序
10、分量： (0)a(0)()ba (0)()ca其中： ；123jej 24032jej；1于是有：（2-2）(1)(2)(0)2(1)(2)(0)aabacaFF（2-3）(1)22(0)abcFF则有： )1(2)1(abII)1()(acI)2()2(abI )2()(acI)0()()0(acbI其逆关系式为：（2-4）2(1)2(0)31a abcFF 这样根据式（2-3）可以把三组三相对称向量合成三个不对称向量，而根据式（2-4 ）可以把三个不对称向量分解成三组对称量。2.2 对称分量法在不对称短路计算中的应用电力系统的正常运行一般是对称的，它的三相电路的参数相同，各相的电流，电压对
11、称，这就电力系统分析课程设计（论文）是说只有正序分量存在。当电力系统的某一点发生不对称故障时，三相电路的对称条件受到破坏，三相对称电路就成为不对称的了。此时，可用对称分量法，将实际的故障系统变成三个互相独立的序分量系统，而每个序分量系统本身又是三相对称的，从而就可以用进行电路计算了。图 2-2 简单系统单相接地故障图如图 2-2 所示的简单系统发生单相接地短路故障。应用对称分量法，可绘出三序网图（三序等值电路图） ，如图 2-3 所示为最简化的三序网图，三序网的参数可分为正序，负序，零序。图中分为正序阻抗，负序阻抗，零序阻抗。120,Z图 2-3 简化三序网图列出电压方程： 11aaEIZU（
12、2-22aa5）00aaIZU由此可见，应用对称分量法进行不对称故障计算时，其关键问题是先求出各序网络的等效电抗（即要求出系统中各主要原件（发电机，变压器，线路等）的各序电抗值） ，然后根据短路的类型，边界条件，把正，负，零序网连接成串，并联的形式，从而可求出电流，电压的各序分量，再应用对称分量法进而可求出各相电流和电压等.【3】3 简单不对称短路的分析与计算（解析法）电力系统发生不对称故障时，无论是单相接地短路，还是两相短路，两相短路接地，只是短路点的电压，电流出现不对称，利用对称分量法将不对称的电流电压分解为三组对称的序分量，由于每一序系统中三相对称，则在选好一相为基准后，每一序只需要计算
13、一相即可，用对称分量法计算电力系统分析课程设计（论文）电力系统的不对称故障。其大概步骤如下：（1）计算电力系统各个 原件的序阻抗；（2）制定电力系统的各序网络；（3）由各序网络和故障列出对应方程；（4）从联立方程组解出故障点电流和电压的各序分量，将相对应的各序分量相加，以求得故障点的各相电流和电压；（5）计算各序电流和各序电压在网络中的分布，进而求得各指定支路的各相电流和指定节点的各相电压。 【4】3.1 单相（a 相）接地短路单相接地短路时的系统接线图如图 3-1 所示，故障处的三相的边界条件： 0aUbI0cI图 3-1 单相接地短路时的等值接线图 用对称分量表示为： 120aaU2II（
14、3-1）120aa即有：（3-1203aaII2）根据单相接地短路时的边界条件（式 3-1，式 3-2）连接复合网，如图 3-2 所示。由复合网图可以写出各序分量：图 3-2 单相接地短路时的复合序网电力系统分析课程设计（论文）120120aaEIIZUE（3-22aI3） 00aIZ于是可以用对称分量法得到短路点的各相电流电压： 12013aaaIIIbc120aaU221 01baUIZaZ（3-4）2102ca短路点电流，电压的相量图如图 3-3 所示。这里按纯电感性电路画的，电流滞后电压 ，若不是90纯电感电路，则电流与电压角度由 的阻抗角决定，一般小于 。在相量图中，将每20Z相的序
15、分量相加，得各相电流电压的大小和相位。图 3-3 单相接地时短路处的电流，电压的相量图3.2 两相（b,c 相）短路两相短路的系统接线如图 3-4 所示，在 k 点发生 b，c 两相短路。短路点的边界条件：图 3-4 两相短路等值接线图0aIbcI0bU0bc用对称分量表示为： 120aaII电力系统分析课程设计（论文）（3-2 212010()aaaaIIII5）2 212010aaaUU于是有：(3-6)12aI0aI12aU由式（3-6）可知，故障点不与大地相连，零序电流无通路，因此无零序网络。复合网络是正负序网并联后的网络。如图 3-5 所示：图 3-5 两相短路的复合序网图从复合序网
16、中可以直接求出电流，电压的各序分量： 121aEIZ（3-7）122aaUIIZ由对称分量法可求得短路点各相电流和电压为： 120aaII2111()3bcaaaIIjI202aUZ（3-8）21bcaaaU短路点电压和电流的相量图。如图 3-6 所示，这里任然是纯电感电路。电流滞后电压 90图 3-6 两相接地短路点电流和电压相量图3.3 两相（b 相和 c 相）短路 接地 两相短路接地时系统接线图如图 3-7 所示，短路点的边界条件为：电力系统分析课程设计（论文）图 3-7 两相短路接地的等值接线图0aI0bcU用对称分量表示为： 120aa（3-9）13a由（3-9）可以画出两相短路接地
17、的复合序网图是三个序网并联，如图 3-8 所示，根据复合序网可求出电流，电压各序分量：图 3-8 两相短路接地的复合序网图123/aEIZ（3-10）0212aI2010aZI201201 1aa aUIEIZZ 用对称分量法合成各相电流电压为： 120aaII电力系统分析课程设计（论文）2021201baaZaIII（ 3-11）202101caaIIIZ12013aaaUUbc短路点流入地中的电流为：（3-12）2003bcaZII短路点电压，电流的相量图，如图 3-9 所示。这里任然是纯电感电路。电流滞后电压 90图 3-9 两相接地短路处的电压电流相量图4 简单不对称短路的分析与计算计
18、算机计算程序法4.1简单故障的计算程序原理1如果要求准确计算故障前的运行情况，则需要进行潮流计算，在近似使用计算中，对短路故障可假设节点U I0I.均为1 IGI2这里采用形成节点导纳矩阵的方法。发电机的正序电抗用X dn,可计算故障后瞬时量。发电机的负序的电抗近似等于X dn.当计算中不计负何影响时，在正，负序网络中不接负何阻值。如果记及负何影响时，负何的正序阻抗可通过其额定功率和电压计算。负何阻抗很难却确定，一般取X 2 =0.35以负何额定功率为基准。负何的中性点一般不接地，零序无通路3形成三个序网的节点导纳矩阵后，利用公式3-19和高斯消去法可求得 ,即为1UnZ1f Znf .对于短
19、路故障，只需令 =1 （其余节电电流均为零） ，分别应用高斯消去法求解一fI次所得电压，即为三序网和f点有关节点阻抗（4）根据不同短路故障，可分别应用 表3.1中 相应公式计算故障处各序电流，电压，进而合成得三相电流，电压。 （5）计算网络中任一点电压时，电力系统分析课程设计（论文）负序和零电压只须计算由故障点电流引起的电压。对于正序则还需要加上正常运行时的电压表4-1不对出故障处的各序电流电压计算公式。故障类型故障端口各序电流公式 故障端口各序电压公式单相接地短路021fffIIffff ZZU3021)(12ffff ZIUff00fffI两相短路 fffff ZI21)0(21 1fff
20、f Z22fffIU两相接地短路ffffffZUI 302)1( fffff ZI0)2(12fffffZI30)2(10 101ffff22fffZI00fff4.2 网络节点方程的形成计算机编程计算中，考虑了对地电容和变压器的实际变比 7 导纳计算公式如下：对角线元素：（ j I, Ui=1,Uj=0,i j） （4-1）iiiYIUY非对角线元素 ：（U i=1,Uj=0,i j） （4-2）jiijijyI式中 i， j-第 i , j 节点 Ii， ， Ij ，U i ，U j-第 i,j 节点电流电压变压器支路的导纳变化 （高压侧） （4-3）2KYTi（4-4） yij4） 网络
21、节点方程的形成矩阵表达式 nnnY111电力系统分析课程设计（论文）网络节点方程： （4-5）UYI（4-6） UYIInnn111.5 电力系统不对称短路计算实例5.1 例 1.已知各元件参数如下：发电机 G1 的容量为 100MVA,G2 为 200MVA,额定电压均为 10.5KV,次暂态电抗 均为 0.2；变压器 T1 的容量为 100MVA，T2 为 200MVA,变比均为 10.5/115KV，短路电“dX压百分数均为 10；三条电力线路（ 、 、 ）的参数均为 115KV,60km，电容为 0.008F/km,电1l23l抗 ；负荷 L1 为 50MW， ；负荷 L2 为 100
22、MVA, 。两台发电机10.4/kmxcos0.985cos1中性点均不接地；两台变压器为 YN,d 连接，发电机侧为三角形连接，YN 侧中性点直接接地；三条电力线路的零序电抗均为 0.20（以 50MVA 为基准值）试求 b、c 两相接地短路时故障点的短路电流（ ）0fgZ图 4.10 电力系统图方法一：解析法1 元件标幺值的计算取 ,50()BMVASBavU发电机发电机 G1 正序电抗标幺值由式（1.4）得 。发电机正序电抗标幺值等于负序1()50.2.GX电抗标幺值等于零序电抗标幺值，即 1()发电机 G2 正序电抗标幺值由式（1.4）得 。发电机正序电抗标幺值等于负2(1)005序电
23、抗标幺值等于零序电抗标幺值，即 ()2()GGXX变压器变压器 T1 正序电抗的标幺值由式（1.6）得 。变压器正序电抗标幺值等于1()50T电力系统分析课程设计（论文）负序电抗标幺值等于零序电抗标幺值，即 1()(2)1(0).5TTXX变压器 T2 正序电抗的标幺值由式（1.6）得 。变压器正序电抗标幺值等()于负序电抗标幺值等于零序电抗标幺值，即 21(2)(0).2TT线路线路 正序电抗标幺值由式（1.8）得 。1l1()250.46.1lX线路正序电抗标幺值等于负序电抗标幺值 ；线路零序电抗标幺值为 0.2（已知）()1()ll线路 正序电抗标幺值由式（1.8）得 。线路正序电抗标幺
24、值等于负2l2(1)2500.46.1lX序电抗标幺值 ；线路零序电抗标幺值为 0.2（已知）2(1)(2)0.llX线路 正序电抗标幺值由式（1.8）得 。线路正序电抗标幺值等于负3l2(1)2.0.15l序电抗标幺值 ；线路零序电抗标幺值为 0.2（已知）2(1)(2).ll负荷阻抗 2111122112508.7501.7.766LBLBLavavBLjjjjUSSQZRXPUQPS22221050.LBLBLavavBLjjjUSSQZRXPUQPS2 化简网络 本例为电力系统短路的近似计算，可作以下化简：(1)令电源电动势标幺值为 1，即 ;“21GE(2)不计负荷影响，将 略去；L
25、Z化简后的等值网络如图 4.11 所示电力系统分析课程设计（论文）图 4.11 等值网络图3 不对称短路计算正序等值电路图如图 4.12图 4.12 正序等值电路图下面通过对图 4.12 进行简化求得网络的正序阻抗如图 4.13 （a） （b） （c） （d） （e）图 4.13 电路简化流程图（a）（e）由图 4.13（e）得 10.4jZ负序等值电路图如图 4.14电力系统分析课程设计（论文）图 4.14 负序等值电路图（a） (b) (c) (d) (e)图 4.15 电路简化流程图（a）（e）由图 4.15（e）得 20.14jZ3）零序等值电路图如图 4.16图 4.16 零序等值电
26、路图下面对图 4.16 进行化简求得网络零序阻抗电力系统分析课程设计（论文）（a） （b） （c） （d）图 4.17 电路简化流程图（a）（d）由图 4.17（d）得 0.179jZ4）计算故障点各序电流 设 3U则 b、c 两相接地短路时(1.)5152016.410.4.79.aa jjjjIZ(1.)5(1.)520.79j.j3.1aajII（ ）(.) (.)205150.4j6.j2.94jZ（ ）5）计算故障点三相电流实际值b、c 两相接地短路时 (1.)(1.)(.1)(.)25250BaaIII 0.86.4.0.86)(3.42.970.5jjjjj(3.)1KA1.)1
27、.)(1(.)25550BcaaIII(8.4.4.()jj1.)2513cb AI 方法二、计算机程序法1 由式 ， , 得112(0%)/5BUK1TK12TK；/.()；101.055/2 正序导纳标幺值 Y 矩阵电力系统分析课程设计（论文）Y 矩阵的形成：正序电抗标幺值和对地电容标幺值的等值电路图如图(4.18)图 4.18 正序电抗标幺值和对地电容标幺值的等值电路图及 正序导纳矩阵如下：Y1对地导纳20BlUbS2620.134.80.51j j1. 02YY 262.0.134.80.5j j其中， 是图(4.18)节点-得对地导纳之半， 分别为节点-、-Y1.2的对地导纳之半。Y
28、 矩阵中的元素为； ；1013130.5jjy123119.050.5jjyk；45； ；12Y2240 46.0.5jjjjyY 23Y；4139.1.5.jyk2Y； ；1390320Y1524 211. 00.5j jjyyYk4838； ； ；4104239.0j3jY452 21 11. 00.5.j jjyyYk电力系统分析课程设计（论文）；0.436.28056.24jjjj；4510.jY； ； ； ；52Y531jY541jY502104jyy5034. 9.60jjy(1) 1.03.09.58.103.56.24109.jjjjjjjjjjY 3 负序阻抗 Y 矩阵同理得
29、负序阻抗标幺值和对地电纳标幺值的等值电路图如图(4.19)图 4.19 负序阻抗标幺值和对地电纳标幺值的等值电路图及 负序导纳矩阵如下Y)2(2)3019.05063.19.58.1.6.240109.jjjjjjjjjjjY 4 零序 Y 矩阵同理零序电抗标幺值的等值电路得图如图(4.20)所示电力系统分析课程设计（论文）图 4.20 零序电抗标幺值的等值电路图及 负序导纳矩阵如下Y)0(31(0)()5()34(0)2 21128.4.05jjjyyYk； ； ；34(0)34(0).j()3(0)5.jy3(0)(0)15.2jyY；2()()45(0)21116.284685jjjjj
30、Yk；45(0)45(0).jy； ；3()()5.2Y4(0)54(0)15.2jy(0)53()4(0)1jyY(0)8.1546.2510jjjY 5 求不对称电流及故障电压由式（3.19）及 3.2.1 节相关知识，运用高斯消去程序法11解得,K5 点发生两相短路时的转移阻抗值：15()23()451()0.73.64820.jjjZZ 15(2)3()452()0.4137.680.15jjjZZ15(0)23()450().814.2jjZZ对于 b、c 两相接地短路，有 5(2)0(1.)().10.8.52jjj则5()(.) 2 2(2)0.150.79331.5Zm电力系统
31、分析课程设计（论文）5(1)(1.)5.16.400.3.5aa jjjUIZ(.)1640(在此 ).(.)(1.)55.640.2.)Bcab KAmII0.25BKAI5.2两种计算方法的对比下面将解析法与计算机程序法进行简单比较。将例子结果如表5-1下：表5-1解析法与计算机程序法数据故障处短路电流 短路冲击电流（KA） 短路容量（MV.A ）比较项目短路类型1 2 1 2 1 2单相接地短路 5.803 5.868 7.4282 7.5047 508.3 568.3两相短路 2.657 2.648 3.402 3.389 265.7 264.81-代表解析法结果；2-代表计算机程序法
32、结果。由例子计算过程以及上面的表格数据及比较可得：解析法误差大，计算过程数据大多需要进行四舍五入。计算机程序法中，利用Y矩阵计算时考虑对地电容，变压器实际变比误差小，通过变换式（3-19）中的电流的取值可以将网络各个节点转移阻抗求出，使分析其他点的短路故障提供了更容易更直观的参数值，因此计算机程序法有通用性强的特点。结 语电力系统发生不对称短路后，由于短路点对地故障支路的不对称，使得整个网络电流电压三相不对称：。本论文解决不对称短路的问题核心是对称分量法。根据对称分量法采取的具体方法之一是解析法，即把该网络分解为正，负，零序三个对称序网，这三组对称序分量可分别按对称的三相电路分解。然后将其结果
33、叠加起来。求解不对称短路。首先应该计算各原件的序参数和画出等值电路。然后制定各序网络。根据不同的故障类型，确定出以相分量表示的边界条件，进而列出以序分量表示的边界条件，按边界条件将三个序网联合成复合网，由复合网求出故障处各序电流和电压，进而合成三相电流电压。方法二是计算机程序法。通过计算机形成三个序网的节点导纳矩阵，然后利用高斯消去法通过相应公式对他们进行数据运算，即可求得故障端点的等值阻抗。最后根据故障类型选取相关公式计算故障处各序电流电压，进而合成三相电流电压。通过上面两种方法放入计算进行对比可以得出，计算机程序法比较解析法具有计算过程及结果更精确，通用性更强的优点。电力系统分析课程设计（
34、论文）参考资料1 孙丽华电力工程基础北京：机械工业出版社，20062 李光琦电力系统暂态分析(第二版)北京:中国电力出版社,19953 于永源，杨绮雯电力系统分析(第二版)北京:中国电力出版社,20044 刘万顺. 电力系统故障分析北京：中国电力出版社，19995 刘从爱电力工程基础济南：山东科学技术出版社，19976 电气工程师手册第二版编辑委员会电气工程师手册北京：机械工业出版社，20007 西安交通大学，电力工业部西北电力设计院，电力工业部西北勘测设计院短路电流实用计算方法.北京：电力工业出版社，19827 何仰赞 温增银，电力系统分析(第三版).武汉 华中科技大学出版社， 20019
35、刘万顺电力系统故障分析北京：中国电力出版社，199810 张俊才电力系统故障的分析计算一对称分量法与相分量法江西工业大学11 邱关源 电路 北京 ：高等教育出版社， 199912 罗建军 杨琦. MATLAB 教程 北京： 电子工业出版社 200513 潘小辉 陈强 MATLAB 5.1 全攻略宝典. 北京： 中国水利水电出版社， 200014 Billinton Roy.ALLAN RONALD (RONALD NORMAN ).RELIABILITY EVALUATION OFPOWER SYSTEMS;NEW YORE ;PLENUM PRESS ;199615 IJ.NAGRATH.D
36、.P.K.OTHARI.MODERN.POWER SYSTEM ANALYSIS.NEW DELHI:TATA MCGRAW-HILLCO.TD.1980附录：不对称短路电流计算程序#include “stdio.h”#define SB 100#define Uav 115Main()int K,j,n,d,S Float x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,Y77,U7,I7,UK,SN,P,Q,X,b,L,Y1,Y2,Y3,K1,K2,K3,IK,IM,SD,PN;K1=K2=1.05;K3=0.96;电力系统分析课程设计（论文）Clear();
37、电抗.对地导纳Printf(“please input the patament of generater :Xd,PN,cosan”);Scanf(“%f,%f,%f”,X1=Xd*Sn/(PN/X);Scanf(“%f,%f,%f”,X10=Xd*Sn/(PN/X);Printf(“X1=%f,X10=%fn”,X0,X10);Printf(“please input the putamenter of tramsfer :UK and SNn1:”);Scanf(“%f,%f,”,X2=UK*Sn/SN;Printf(“n2:n”);Scanf(“%f,%f,”,X3=UK*Sn/SN;
38、Printf(“n3:n”);Scanf(“%f,%f,”,X4=UK*Sn/SN;Printf(“X2=%f,X3=%f,X4=%fn”,X2,X3,X4);Printf(“please input the putamenter of load:P,Qn”);Scanf(“%f,%f,”,X8=1/(Q/100)*Q/100+(P/100)*(P/100)*Q/100);X11=0.35;X12=1.2;Printf(“ please input the M of load:P,cosan”);电力系统分析课程设计（论文）Scanf(“%f,%f,”,X9=1/6.5*100/(P/X);P
39、rintf(“X8=%f,X9=%fn”,X8,X9);Printf(“please input the putamenter of line: n”);Printf(“for example:X=0.39,b=02.92e-6,i=100”);Printf(“n1:”);Scanf(“%f,%f,”,X5=X*L*Sn/(Uav*Uav);Y1=(double)12*b*L*(Uav*Uav)/100;Printf(“n2:”);Scanf(“%f,%f,”,X6=X*L*Sn/(Uav*Uav);Y2=(double)12*b*L*(Uav*Uav)/100;Printf(“n3:”);S
40、canf(“%f,%f,”,X7= X*L*Sn/(Uav*Uav);Y3=(double)12*b*L*(Uav*Uav)/100;Printf(“X5=%f,X6=%f,X7=%fn”,X5,X6,X7);Printf(“Y1=%f,Y2=%f,Y3=%fn”,Y1,Y2,Y3);Printf(“Y1=Y1,Y2=Y2,Y3=Y3n”,Y1,Y2,Y3);Printf(“please input the K(7) to see Y array: n”);Scanf(“For(n=0;nK;n+) 导纳矩阵形成电力系统分析课程设计（论文）For(j=0;jK;j+)Ynj=0;Y00=1/X
41、1+1/X2;Y01=-1/(X2*K1);Y10= Y01;Y11=1/X5+Y1+1/(X2*K1*K1);Y12=-1/X5;Y21= Y12;Y22=Y1+1/X5+Y2+1/X7+(1/2*X3*K2*K2);Y23=-1/(1/2*X3*K2);Y24=-1/X6;Y25=-1/X7;Y32= Y23;Y33=1/(1/2*X3)+(1/2*X10);Y42= Y24;Y44=Y2+1/X8 +1/X6;Y55=Y3+1/X7+1/(X4*K3*K3);Y56=-1/(X4*K3);Y65= Y56;Y66= 1/X4;For(n=0;nK;n+)For(j=0;jK;j+) Pr
42、intf(“j%,2ft”,Ynj);电力系统分析课程设计（论文）If(j=6)Printf(“n”)For(n=0;nK;n+) 短路电流和功率Un=0.0;Printf(“please input the number of short cuicit:K=I(i7)n”);Scanf(“Ud=1;I0=0.0;For(n=0;nK;n+)For(j=0;jK;j+)In= Ynj *Uj;N=0;Printf(“I0=%fn”,In+);Printf(“I1=%fn”,In+);Printf(“I2=%fn”,In+);Printf(“I3=%fn”,In+);Printf(“I4=%fn
43、”,In+);Printf(“I5=%fn”,In+);Printf(“I6=%fn”,In+);IK=Id*Sn/(sqrt(3)*Uav);IM=2.55*IK;SK=1.732*IK*Uav;Printf(“IK=%f,Sd=%f,IM=%f”,IK,Sd,IM);电力系统分析课程设计（论文）Printf(“please input the K(7) to see Y array: n”);Scanf(“ 负序导纳矩阵For(n=0;nK;n+)For(j=0;jK;j+)Ynj=0;Y00=1/X1+1/X2;Y01=-1/(X2*K1);Y10= Y01;Y11=1/X5+Y1+1/
44、(X2*K1*K1);Y12=-1/X5;Y21= Y12;Y22=Y1+1/X5+Y2+1/X7+(1/2*X3*K2*K2);Y23=-1/(1/2*X3*K2);Y24=-1/X6;Y25=-1/X7;Y32= Y23;Y33=1/(1/2*X3)+(1/2*X10);Y42= Y24;Y44=Y2+1/X8 +1/X6;Y55=Y3+1/X7+1/(X4*K3*K3);Y56=-1/(X4*K3);Y65= Y56;Y66= 1/X4;电力系统分析课程设计（论文）For(n=0;nK;n+)For(j=0;jK;j+) Printf(“j%,2ft”,Ynj);If(j=6)Print
45、f(“n”)Scanf(“零序导纳矩阵For(n=0;nK;n+)For(j=0;jK;j+)Ynj=0;Y00=1/（X2*K1*K1）+1/(3*X5);Y01=-1/(3*X5);Y10= Y01;Y11=1/(3*X5)+1/(3*X6)+1/(3*X7)+1/(1/2*X3*K2*K2);Y12=-1/(3*X6);Y21= Y12;Y22= 1/(3*X6)+Y2+1/X12;Y13= Y31;Y33=1/(X4*K3*K3)+1/(3*X7);Y31=-1/(3*X7);Y43= Y34;Y44=1/X4;For(n=0;nK;n+)For(j=0;jK;j+)电力系统分析课程设计（论文） Printf(“j%,2ft”,Yn
