电力系统稳定计算

1、目录一、需求分析 21.1 概述 21.2 安全稳定计算分析的目的 3二、概要设计 52.1 文档概述 52.2 计算方法 52.2.1 暂态稳定计算 5三、详细设计 63.1 定义 63.2 求解过程 73.3 系统的稳定计算框图 73.4 模块设计说明 93.4.1 模块描述 93.4.2 暂态稳定计算的过程： 93.4.3 C 语言封装 11四、使用说明 154.1 功能及性能 154.2 输入项 154.2.1 自动读取 excel 表格数据模式 154.2.2 手动键入数据模式 164.3输出项 164.4 算法 164.5 流程逻辑 16五、测试报告 175.1 测试报告简介 17

2、5.1.1 程序测试目的 175.1.2 测试参考样本 175.1.3 概述 175.1.4 测试用例设计步骤 185.2 测试摘要 185.2.1 测试用例 185.2.2 计算机计算过程 21需求分析1.1 概述电力系统的根本任务是合理开发能源， 向电力用户提供安全可靠、 质量合格和经济的电 能。在这三个基本要求中，安全可靠供电是第一位的要求。就是说，必须首先保证电力系统 的正常运行， 而正常运行所不可缺少的最基本条件是安全和稳定。 所谓安全， 是指系统中所 有的电力设备必须在不超过其允许的电压、 电流和频率条件下运行， 否则会造成设备的损坏。 所谓稳定， 是指电力系统经受扰动后能继续保持

3、向负荷正常供电的状态， 即具有承受扰动的 能力。本章的目的是进行电力系统的稳定计算。电力系统运行时， 有三种必须同时满足的稳定性要求： 同步运行稳定性、 频率稳定性和 电压稳定性。 电力系统的同步运行稳定性， 又称功角稳定性， 是指电力系统中所有发电机组 能否保持同步稳定运行的问题。 当系统在某一正常状态下运行时， 系统中所有发电机是保持 同步的，它们的电气角速度一样，它们间的角度差为一常数， 系统中每一点的电压、电流及 功率也为常数。 当受到某种扰动时， 这些运行参数会发生变化，如经过一段时间后，系统能 够回到原来的运行状态或者过渡到一个新的正常运行状态， 则系统是同步稳定的。 同步稳定 的

4、标志是各个发电机之间的功角差s?沁§?能保持有限值而不随时间无限增大。如果系统失去了同步运行稳定性， 就会发生振荡， 引起系统中各点的电压、 电流和功率大幅度周期 波动， 从而无法向负荷正常供电， 严重时将会造成电力系统大面积停电。 电力系统的频率稳 定性是指能否保持电能的频率指标正常。 频率是电能质量的一个重要指标， 它反映了电力系 统中有功功率的平衡水平。由于电能不能大量储存，所以电能的生产、 传输、 分配和消费是 同时进行的， 即系统中所有发电厂任何时刻生产的电能必须与该时刻所有负荷所需的电能及 传输、分配中损耗的电能之和相等，也就是必须满足有功平衡的等约束条件：？?沪?寸?。

5、在这种情况下，全系统的频率为一定值，且应维持在规定的允许范围内。如系统的有功 电源不足，则频率下降，当频率下降到一定程度时，将引起系统崩溃，失去频率稳定性，造 成全系统停电。 电力系统的电压稳定性是指能否保持电能的电压指标正常。电压是电能质量的又一个重要指标， 它反映电力系统中无功功率的平衡水平。 电力系统正常运行时， 全系统 包括每一地区的无功功率应处于正常的平衡水平， 从而各点电压在规定的允许范围之内。 如 全系统或某一地区的无功电源不足， 将引起全系统或该地区电压水平的降低， 当电压降到一 定程度将引起电压崩溃，失去电压稳定性，造成受影响的地区停电。随着电力系统向大机组、大电网、高电压和

6、远距离输电的发展，电压稳定问题更突出地暴露出来。近年来国际上 有多次由于电压稳定问题引起的大面积长时间的停电事故发生。因此电压稳定问题引起了电力工业界更多的关注和重视。根据2004年IEEE电力系统稳定术语、定义工作小组的建议，对电压稳定的定义和分类有如下论述。 电压稳定指的是在一个给定的初始运行状态承受扰动 后维持所有母线稳态电压的能力。电压稳定问题分为大扰动电压稳定和小扰动电压稳定。大扰动电压稳定指的是在大扰动下（如系统故障、发电机跳闸等）维持稳态电压的能力。小扰动电压稳定指的是系统承受小扰动 （如系统负荷变化） 维持稳态电压的能力。 电压稳定过程 的时间可以从几秒变化到几十秒。因此电压稳

7、定既可以是短期现象也可以是长期现象。短期电压稳定涉及到快速作用的负荷元件，如感应电动机、电子控制负荷和HVDC转换器。长期电压稳定涉及较慢作用的设备，如变压器抽头变化等。目前对电压稳定的研究内容主要包括对电压崩溃现象机理探讨、电压稳定安全分析以及预防措施、 电压稳定研究的负荷模型等。 对电压稳定分析的主要方法有， 基于物理概念的定 性分析、 基于潮流方程的静态分析方法、 基于线性化动态方程的小干扰法分析、 基于非线性 动态方程的时域仿真方法。 对以上分析方法的研究有大量的研究成果， 但对电压失稳机理的 认识并还没有达到一致，有待进一步研究。以上三种稳定性中， 同步运行稳定性是最常发生、 最受关

8、注也研究得最多的一种稳定性， 本章分析的就是这类稳定性。由于同步运行稳定性是由发电机转子运动的功率角3表征的，因而同步运行稳定性计算的目标就是求取各发电机受扰后功率角随时间的变化情况而后进 行稳定性判别。 为此必须首先了解电力系统中各个旋转元件的机械和电气特性，简称机电特性。通常它们由一组非线性微分方程和代数方程描述， 所以分析同步运行稳定性的一般方法 就是：求解一组在一定初始条件下的非线性微分方程和代数方程，得到发电机发功率角3随时间变化的曲线3 (t) ，称为摇摆曲线，从而判断其稳定性。由于电力系统的复杂性，可以 想见， 这种大量非线性微分方程和代数方程的求解， 只能借助计算机， 采用一定

9、的方法求解 数值解，这类方法称为时域仿真法，或逐步分析法( SBS,step by step )。由于其工作量很 大，所以人们一直在寻求一些简化方法： 如对小扰动， 可将非线性微分方程在运行点附近线 性化，然后按线性系统的方法，根据其特征根在复平面上的位置判断稳定，称为小扰动法； 或者更简单地， 找到一些判据直接判断稳定； 对大扰动， 设法找到一种描述系统运动的能量 函数， 从能量的角度判断其稳定性从而避免求数值解， 称为直接分析法。 这些就是分析电力 系统同步运行稳定性的方法。 由于同步运行稳定性是电力系统在受到扰动后发电机转子运 动的稳定性， 所以它和扰动的大小有关。 人们通常根据扰动的大

10、小将其分为静态稳定性和暂 态稳定性。 电力系统的静态稳定性是指： 如果在小扰动后系统达到扰动前一样或相近的稳定 运行状态， 则称系统对该特定运行情况为静态稳定， 也称为小扰动下的稳定性。 电力系统的 暂态稳定性是指： 如果在大扰动后系统达到允许的稳定运行情况， 即仍能继续保持同步运行， 则称系统对该特定运动情况和对该特定扰动为暂态稳定，也称为大扰动下的稳定性。1.2 安全稳定计算分析的目的电力系统安全稳定计算分析的目的是通过对电力系统进行详细的仿真计算和分析研究， 确定系统稳定问题的主要特征和稳定水平， 提高系统稳定水平的措施和保证系统安全稳定运 行的控制策略，庸医知道电网规划、计划、建设、生

11、产运行以及科研、试验中的相关工作。安全稳定计算分析的总体要求 电力系统安全稳定计算应根据系统的具体情况和要求， 开展对系统的静态安全分析、 静 态稳定计算、 暂态稳定计算、 动态稳定九三、 电压稳定计算、 频率稳定计算以及再同步计算， 并对计算结果进行认真、 详细的分析， 研究系统的基本稳定特性， 检验规划电网的安全稳定 水平，优化电网规划方案， 提出保证电网安全稳定运行的控制策略和提高系统稳定水平的措 施。暂态稳定计算分析暂态稳定是电力系统受到大扰动后， 各同步电机保持同步运行并过渡到新的或恢复到原 来稳态运行方式的能力， 通常保持第一、 第一第二摇摆不失步的功角稳定， 是电力系统功角 稳定

12、的一种形式。暂态稳定计算分析的目的是在规定的运行方式和故障形态下，对系统的暂态稳定性进行校验，研究保证电网安全稳定的控制策略， 并对继电保护和自动装置以及各种安全稳定措施 提出相应的要求。电力系统的暂态稳定性 电力系统的暂态稳定性是指电力系统在某一运行状态下收到大扰动后能否继续保持同 步运行。如受到扰动后任意两台发电机之间功角差3 ij 随时间一直增大， 则不稳；如3 ij 不随时间无线增大，而保持有限值，则称该运行情况对特定扰动为暂态稳定。 针对上述电力系统暂态稳定性的定义，有如下几点说明：（ 1） 定义中的大扰动指短路、 断线或线路以及符合的突然大变化等。 上述大扰动中 一般以短路挂号找你

13、噶， 特别是三相短路最严重。 各国考个电力系统咱题啊稳 定性的标准不一。 我国现对 220kV 以下电网，以三相短路为参考标准； 对 220kV 以上电网， 以单相瞬时接地、 故障后切除故障线路并重合闸成功作为系统必须 能承受的一组扰动。（ 2） 电力系统机电暂态稳定持续时间的长短既与扰动的大小有关，也与系统本上 的状况有关， 有的持续时间短， 有的持续时间长， 随所研究的问题的不同可对 不同长度的时段进行研究。 为此暂态分析时可按时间将其分为三个阶段， 在不 同阶段考虑的因素有所不同。初始阶段：一般指故障 1s 内时间的阶段。此时电力系统的继电保护和断路器 动作， 例如切除故障线路、 重合闸

14、、 切除发电机等， 但发电机的调节系统特别 是调速系统所起的作用不明显。中间阶段：指从 15s 后的时间段，此时需要考虑发电机的调节系统。 后期阶段：指 5s 以后的时间段。此时需考虑动力部分的变化产生影响、系统 频率变化的影响以及低频捡灾等自动装置的作用等。从而暂态分析可分为短期、 中期和长期三种， 本节对暂态稳定性的分析计算以 起始段为主，兼及发电机调节系统的影响。（ 3） 暂态稳定计算分析中常用假设主要有：1）不计故障电流非周期分量的作用。这一方面是因为其衰减很快，大约0.1s 左右的时间即衰减完毕；另一方面，其产生在空间不懂磁场，它 与转子绕组中的直流电流互相作用所产生的转矩以同步频率

15、周期变化， 平均值接近于 0，对发电机转子的运动影响不大。同时，其为制动性 质，故忽略它后使计算结果偏于保守。2）不计负序和零序电流的作用。由于负序电流产生的刺伤与转子绕组相互作用形成的转矩以两倍频交变，平均值也接近0；零序电流昌盛的合成磁场为零，不产生转矩，故可将其忽略。3）不计阻尼功率PD。这一方面因为其值很小，对转子运动影响不大；另一方面忽略它使计算偏于保守 （以为其制动性质） ，易于为工程界接受。（ 4） 暂态稳定的分析方法从本质上讲师求解一组描述电力系统机电特性的方程， 包括微分方程和代数方程，得到发电机功率角随时间变化的曲线摇摆曲 线，从而判断系统的稳定性。 为取微分方程的值解有多

16、重， 但均是将时间分成 小段求取各个时间的值。 因此这类方法称为逐步积分法或时域仿真法。 可以想 见，其工作量大， 十分费时。 为此人们一直在寻求一种避免求微分方程数值解 而直接判断稳定性的方法，称为暂态稳定性的直接分析法。 本文以单机无穷大系统的暂态稳定性为主要研究对象。、概要设计2.1 文档概述文档主要描述单机无穷大系统的概要设计思想。2.2 计算方法2.2.1 暂态稳定计算预防-校正法 先用欧拉法求得一个初步的近似值，称为预报值，然后用它替代梯形法右端的yi 1再直接计算 fi 1，得到校正值 yi 1，这样建立的预报 -校正系统称为改进 的欧拉格式： 预报值yi 1 yi h f xi

17、 ,yi校正值 yi 1 yi (h/2)f xi ,yif xi1,yi 1它有下列平均化形式：ypyi hf xi ,yiyc yihfxi 1,ypyi1 (ypyc)/2它的局部截断误差为O(hW),可见，改进欧拉格式较欧拉格式提高了精度，其截 断误差比欧拉格式提高了一阶。改进欧拉法改进欧拉法是对欧拉算法的改进方法。微分方程的本质特征是方程中含有导数项,U(tn 1) U(tn)f(t,U)dtg 3 f (tn，u(tn )f(tn 1,U(tn 1)2数值解法的第一步就是设法消除其导数值，这个过程称为离散化。实现离散化的基本途径是用向前差商来近似代替导数，这就是欧拉算法实现的依据。

18、欧拉(Euler) 算法是数值求解中最基本、最简单的方法，但其求解精度较低，一般不在工程中 单独进行运算。模块概述1、Eq方程系数模块：求解Eq方程的变量的系数；2>A，Eq的函数模块：求解, Eq的函数；3、系数矩阵模块：求解状态方程组的系数矩阵；4、方程的特征值模块：求解方程的特征值；5、系统稳定性判断模块：判断系统是否稳定；6摇摆曲线模块：绘制系统摇摆曲线。二、C语言模块C语言模块的实现原理和 matlab模块是相似的，考虑到 C语言封装之后无需在 matlab环境下运行，提高了程序的兼容性，所以我们尝试实现了C语言模块。三、详细设计3.1定义u：系统电压EqO:空载电势UG0机端

19、电压Xe：线路电抗Xd:同步电抗Xza n：暂态电抗Xdxt:系统同步电抗Xxtzan ：系统暂态电抗Tj ：惯性时间常数TdO：励磁绕组时间常数f ：系统频率P1：正常时最大功率P2:故障时最大功率P3:故障切除后最大功率n步长CutTime :输入故障切除时间3.2 求解过程 暂态稳定计算的求解过程大致可以分为以下步骤: 利用改进欧拉法求解转子运动方程。 静态稳定计算的求解过程大致可以分为以下步骤:1 ）列写方程并线性化；2）将?、?和?表为状态变量的函数；3）形成系数矩阵 A 并判断系统稳定性。3.3 系统的稳定计算框图图一、电力系统暂态稳定计算流程图3.4 模块设计说明3.4.1 模块

20、描述节点导纳矩阵模块 :根据所给题目要求 ,输入相应矩阵 ,运行程序得出 所需输入数据如下：系统电压 :空载电势 :机端电压 :线路电抗:同步电抗 :暂态电抗:系统同步电抗 :系统暂态电抗 :惯性时间常数 :励磁绕组时间常数 :系统频率 :正常时最大功率 :故障时最大功率 :故障切除后最大功率 :输入故障切除时间：3.4.2 暂态稳定计算的过程：预报校正法：求初始变化率?/? ?-1| = ?( ?-1，?-1)预报?= ?-1 + ?/? ?-1|?求终了变化率?/?= ?(? ?校正?= ?-1 + (?/?+ ?/?1)?/2转子运动方程d s? - ? - ?dt?取角度？的单位为度，

21、故？?>?= 360?= 18000 0/s (?= 50Hz)。Matlab程序如下：Wn=100*pi;Xe=XT1+XI+XT2;%线路电抗Xd_al匸Xe+Xd;Xdd_al匸Xe+Xdd;l=(P0-Q0*j)/(Us);Ed1= Us+Xdd_all*l*j;Ed=sqrt(reaI(Ed1)A2+imag(Ed1)A2); %计算暂态电抗后电势disp(Ed);Delta(1)=ata n( imag(Ed1)/real(Ed1); %计算初始功角disp(Delta(1);S(1)=1;%初始转速为1t(1)=0;Xl0=4*Xl;Xhe2=(X2+XT1)*(Xl+XT

22、2)/(X2+XT1+Xl+XT2);Xhe0=XT1*(Xl0+XT2)/(XT1+Xl0+XT2);Xdelta=Xhe2*Xhe0/(Xhe2+Xhe0);Xii=(Xdd+XT1)+(Xl+XT2)+(Xdd+XT1)*(XT2+Xl)/Xdelta;Pem1Pem2P2M=Ed*Us/Xii;%Xiii=Xdd+XT1+2*Xl+XT2; P3M=Ed*Us/Xiii;%h=0.001;% 步长Duration=2;%设置计算时段长度 %改进欧拉法计算故障期间相应的最大电磁功率计算故障后取相应的最大电磁功率2sfor i=1:r oun d(CutTime/h) d(i)=(S(i)

23、-1)*Wn;o(i)=(P0-P2M*si n(Delta(i)/Tj;%时段初变化率Delta(i+1)=Delta(i)+d(i)*h;S(i+1)=S(i)+o(i)*h;% 时段末估计值 dd(i)=(S(i+1)-1)*Wn;时段末变化率oo(i)=(P0-P2M*si n(Delta(i+1)/Tj;%Dav(i)=(d(i)+dd(i)/2;Oav(i)=(o(i)+oo(i)/2;% 平均变化率Delta(i+1)=Delta(i)+Dav(i)*h;S(i+1)=S(i)+Oav(i)*h;% 时段末计算值 t(i+1)=i*h;end; s=sprintf('De

24、lta=%fn',Delta(i+1);% disp(s);%故障切除后的功角变化过程for i=round(CutTime/h)+1:round(Duration/h) d(i)=(S(i)-1)*Wn;o(i)=(P0-P3M*sin(Delta(i)/Tj;% 时段 Delta(i+1)=Delta(i)+d(i)*h;S(i+1)=S(i)+o(i)*h;% 时段末估计值 dd(i)=(S(i+1)-1)*Wn;oo(i)=(P0-P3M*sin(Delta(i+1)/Tj;%此时功角时段初变化率时段末变化率Dav(i)=(d(i)+dd(i)/2;Oav(i)=(o(i)+o

25、o(i)/2;% 平均变化率Delta(i+1)=Delta(i)+Dav(i)*h;S(i+1)=S(i)+Oav(i)*h;% 时段末计算值t(i+1)=i*h;end;s=sprintf('Delta=%fn',Delta(i+1);% 此时功角disp(s);%摇摆曲线plot(t,Delta);xlabel(' 时间(秒) ');ylabel(' 功角(度数 )');title(' 发电机摇摆曲线 ');3.4.3C 语言封装exe 文件是一种可在操作系统存储空间中浮动定位的可执行程序，可执行文件分两种一种是后辍名为.C

26、OM另一种就是.EXE。 exe为自执行的文件，不需要用其他软件来打开。故 考虑将本实验程序封装成 exe 文件。具体思路：1、将MATLAB语言m文件编译成c语言c文件 该步骤为关键步骤，影响后续程序是否可执行，执行结果是否正确。使用MATLAB中的 coder 工具对 m 文件进行转译。图1、coder界面¥ jl ¥ B V- F T rLR LII I JiFntfy-PQint KI 的勺 rightmP0TtjUg“ tiQ0<kx 1)Usckbleflx ljXI如jbl叩v 1)XTLda ibled x ljWdaubls(l k L)xddaub

27、l*(l x 1)Xdddu*jblc(l t, L)丫工daubla(l > L)Tjdo jbl沁 x liTolj图2、参数类型设置界面由于MATLAB中各参数可根据后续运算确定自身数据类型，而C语言中却需要提前定义好数据类型才能进行运算，故需要对其中输入参数进行数据类型的定义。本例 中输入参数均为double类型。2、保存编译的好的c文件及其所需的头文件F-© Sourcesinain.ci 5 - xri g htx:rigHJnitilize.c rigM_teriviinate.< rt_nonfinrte.cj“"_ rtGetlnEc；”-_

28、rtGetNaM.cE B Headers' ri g hth卜-rightjnitialize.h ri g ht_type£.h rtr&nfiniteJh ；-r rtGetlrif.h i-""h3、在VS或者其他C语言环境下建立控制台工程; 本例中使用 codeblocks。11ProjectsSymbolsFile! >Q Workspace厶-叫wd0 & Sources” main.c-1 rightci i =.| righ1_initialire<| right_tenTiifr占長!rt_nonfnite.

29、c - rtGetlnf.cI 1"“_ rtGethJaM.cHeaders.|“二|righthI“ rightjnitializhJi2Jright_terminate.h=right typesih“| rt_rodinite-.1"；“mrtGetlnf.h=“rtGetlXlaN.h:rtwtypes.h4、使用主函数调用原 m文件生成的功能函数；由m生成的功能函数需要使用main.c调用才能正确地在 c语言环境下编译。main.c中的程序如下所示（其中 right函数为稳定计算的主要功能函数）：#in elude <stdio.h>#in elud

30、e <stdlib.h>int mai n()double P0;/=1;double Q0;/=0.2;double Us;/=1;double XI;/=0.235;double XT1;/=0.13;double XT2;/=0.108;double Xd;/=2.241;double Xdd;/=0.304;double X2;/=0.218;double Tj;/=8.47;double CutTime;/=0.15;double Ed,s1,s2;printf("请依次输入：初始有功功率，初始无功功率，系统电压，系统电抗，变压器 T1电抗，变压器T2电抗，同步

31、电抗，暂态电抗，负序阻抗，惯性时间常数，切除故障 时间。(数据之间使用空格分隔)：”)；scan f("%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf",&P0, &Q0,&Us, &XI, &XT1, &XT2, &Xd, &Xdd,& X2, &Tj,& CutTime);/prin tf("%f %f %f %f %f %f %f %f %f %f %f',P0,Q0,Us,XI,XT1,XT2,Xd,Xdd,X2,Tj,C

32、utTime); right(P0,Q0,Us,XI,XT1,XT2,Xd,Xdd,X2,Tj,CutTime,&Ed, &s1, &s2);printf("暂态电抗后电势：%fn故障期间功角：%fn故障后功角：%fn",Ed,s1,s2);getchar();getchar();return 0;与原m文件结果对比验证MATLAB中结果：t 1.20240. 5920Belta=0.890526J 1 I -4 2 t 0 0 J 2Exe可执行文件计算结果:n C:U sersAd m i ni &tratorDCi cu ment&am

33、p;Ter cent Files584722961Fi I e Recvwd. exeX请依坎输人初始有功功率，初始无功功率，系统电压，系统电抗，变压器T1电抗，变压器T戈电抗，同步 殿堂，暂态电抗，负序阻抗，惯性时间常数，切除散障时间。（数据之何使用空格分隔）1 0.2 1 a 235 d 13 a 108 Z 241 0. 304 0. 218 & 47 0. 15 匱态电抗后电纵L 392354 故障期间功角初贻値：0.站甜91 故障后功角初如值：6 92927从二者结果对比中可以发现，数值计算的结果并没有差异。但exe无法完成画图工作，不能够提供摇摆曲线的图像（以目前组内成员的

34、水平无法完成此工作）。在从m文件到c文件的过程中，也会有许多问题出现。下面列举部分出现的问题：1、coder只能将m文件的函数文件编译成 c语言文件，故只能以功能函数的形式进行 编译，变成c文件后，再在控制台工程中建立主程序文件去调用该功能，实现所需 功能。2、部分MATLAB中很容易实现的功能在 c语言中难以实现，故只能对源文件中的语句 进行调整，以达到可以编译的效果，比如画图。3、 C语言中的变量需提前定义，不像MATLAB中可以根据计算结果完成数据类型的定 义。虽然有如上列举以及未列举的缺点，但不可否认这样讲程序封装后的优点，即可以脱离完全脱离MATLAB环境，在任意windows环境下

35、进行使用， 而且理论上，由于C语言本 身可移植性很强的特性，以上述方法可以将本程序应用到任何平台。而且本程序的大小只有46.3k，方便传输和携带。而且可以说，通过更深入的学习，很多现在我们无法解 决的问题也能找到方法结局。四、使用说明4.1功能及性能负责接受用户输入，显示系统输出，其他模块功能的协调调用说明对该程序的全部性能要求，包括对精度、灵活性和时间特性的要求。4.2输入项程序运行分为模式：读取excel表格数据模式0，手动键入数据模式1自动读取excel表格数据模式Excel表格形成说明:1）数据表格1页，名称为：n。2）P0为有功功率，Q0为无功功率，Us为系统电压，XI为线路电抗，X

36、T1为变 压器T1电抗，XT2为变压器T2电抗，Xd为同步电抗，Xdd为暂态电抗，X2为 负序阻抗，Tj为惯性时间常数，CutTime为故障切除时间手动键入数据模式1）P0为有功功率，Q0为无功功率，Us为系统电压，Xl为线路电抗，XT1为变压器T1电抗，XT2为变压器T2电抗，Xd为同步电抗，Xdd为暂态电抗，X2为负序阻抗，Tj为惯性时间常数，CutTime为故障切除时间4.3输出项1）计算暂态电抗后电势2）初始功角3）故障期间功角4）故障切除后功角4.4算法本程序所选用的算法，具体的计算公式和计算步骤依据参考资料：电力系统暂态分析（中国电力出版社）。4.5流程逻辑本程序的逻辑流程已由暂态

37、稳定计算流程图详细介绍。五、测试报告5.1 测试报告简介5.1.1 程序测试目的本测试报告的具体编写目的为单机无穷大电力系统暂态稳定计算程序的测 试报告，目的在于总结测试阶段的测试以及测试分析结果，描述暂态稳定计算 程序是否符合要求。针对暂态稳定计算所设计的需求说明书测试用例分析指定 的相应标准过程和要求，找出程序运行过程中可能出现的问题与结果，并希望 通过测试结果了解系统所能实现的细节，进行相应调整并编写成文档便于修改 和维护。5.1.2 测试参考样本参考前期完成的需求说明书、概要设计和详细设计说明书，并以手工计算 的暂态稳定分析结果作为参考来验证暂态稳定计算程序的准确性。5.1.3 概述测

38、试设计遵循与软件设计相同的工程原则。好的软件设计包含几个对测试 设计进行精心描述的阶段。这些阶段是：测试策略 ; 测试计划; 测试描述 ; 测试过 程。这四个测试设计阶段适用于从单元测试到系统测试各个层面的测试。测试设计由软件设计说明所驱动。单元测试用于验证模块单元实现了模块 设计中定义的规格。一个完整的单元测试说明应该包含正面测试和负面的测试 正面测试验证程序应该执行的工作，负面测试验证程序不应该执行的工作。设计富有创造性的测试用例是测试设计的关键。测试用例由测试输入数据 以及与之对应的输出结果组成。 测试用例设计的好坏直接决定了测试的效果 和结果。所以说在软件测试活动中最关键的步骤就是设计

39、有效的测试用例。测试用例设计步骤步骤1 ：首先使被测程序正常运行步骤2:手工计算得到测试用例的结果步骤3:将程序运行结果与手工计算的结果进行比较步骤4:对程序进行修改，使程序运行结果接近手工计算的结果步骤5:程序测试结束5.2测试摘要测试用例为了便于手工运算，测试用例选取一个相对简单的电力系统，图5-1为电力系统的等值电路图。系统等值电路如图所示j0.304j0.130j0.235j0.108U=1.0P0=118 / 21Q0=0.2220选 Sb =250MV- A, U B(110) 115kV , Ub220 1 15 209。各元件参数121标幺值计算如下:Xd2.36250 0.8