课程设计报告( 2014—2015年度第二学期)名称：电力系统暂态上机计算院系：电气与电子工程学院班级：电气1211学号：1121210205学生：郝阳指导教师：艳波设计周数：两周成绩：日期：2015年7月4日一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。

二、设计正文（详细容见附录,用A4纸，页数不限）1．对称短路计算过程流程图和计算结果2．不对称短路计算过程流程图和计算结果3．静态稳定性计算过程流程图和计算结果4．暂态稳定性计算过程流程图和计算结果5．思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。

通过本次课程设计，本人对电力系统故障分析有了更深刻的理解，电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。

在分析过程常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。

在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。

同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。

这与手算系统短路电流时使用的网络化简方法大大不同。

在学习中，参照潮流程序，我加深了对节点导纳矩阵建立方法的理解与学习，巩固了不同类型短路的短路电流计算方法，和序电压、序电流，相电压、相电流基于matlab软件的计算。

利用小干扰分析法判别静态稳定可以使用劳斯判据，也可以使用特征根判别的方法。

二者都能够判断，在书中的例题中使用的是劳斯判据，但是，特征根判别放法更适合在matlab软件的环境下使用，因为该软件提供了计算矩阵特征根的函数，使用起来简便易懂。

利用改进欧拉法计算最大切除角或切除时间虽然比较方便，但是由于计算机有效位数的限制而引起的舍入误差随步长的减小以及运算次数的增多而增大。

在上机编程中体会比较深的是自己犯的几个错误，比如说在算短路电流时，发电机之路的电流和电压是否要经过相位的变换，在静态稳定计算中的角度是弧度制，在计算中使用的都是标幺值等等，这些看似细微的地方，往往考验的就是自己对电力系统分析中的知识的基础，虽然小但是也非常的重要，往往就是这些错误会令程序得不到正确的结果，值得我注意和反思。

总的来说，结合我们上学期潮流上机编程的经验，虽然这次用的是matlab，但是明显感觉到对计算机处理电力系统问题的简便，程序的编写也显得更加的轻松和得心应手，体会到matlab软件的强大。

经过两个星期的回顾与学习，我巩固了电力系统暂态分析理论知识，在编程的过程中很深刻的感受到要想得出正确的结果，就必须认真地理解课本上对应的例题，知道这些原理，不能仅仅照着书抄公式；同时我也初步掌握了采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础，更好地建设我国的电力系统。

四、参考文献1.《电力系统暂态过程》，常鲜戎、书强编，机械工业，2010年1月，第一版；2.《电力系统分析基础》，庚银，机械工业，2011年，第一版；附录（设计流程图、计算结果、思考题答案）1.对称短路计算过程流程图和计算结果流程图：计算结果：（1） 导纳矩阵：正序阻抗：10.0000 -40.0000i 0.0000 +10.0000i 0.0000 +10.0000i 0.0000 +20.0000i 0.0000 + 0.0000i 0.0000 +10.0000i 0.0000 -60.0000i 0.0000 +10.0000i 0.0000 + 0.0000i 0.0000 +40.0000i 0.0000 +10.0000i0.0000 +10.0000i 0.000Y =0 -20.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 +20.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 -30.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i0.0000 +40.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 -60.0000i ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦负序阻抗：20.0000 -40.0000i0.0000 +10.0000i 0.0000 +10.0000i 0.0000 +20.0000i 0.0000 + 0.0000i 0.0000 +10.0000i 0.0000 -60.0000i 0.0000 +10.0000i 0.0000 + 0.0000i 0.0000 +40.0000i 0.0000 +10.0000i0.0000 +10.0000i 0.000Y =0 -20.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 +20.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 -30.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i0.0000 +40.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 -60.0000i ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦零序阻抗：00.0000 -30.0000i 0.0000 + 5.0000i 0.0000 + 5.0000i 0.0000 + 5.0000i 0.0000 -50.0000i 0.0000 + 5.0000i 0.0000 + 5.0000i 0.0000 + 5.0000i 0.0000 -10.0000i Y ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦（2） 节点3发生三相短路故障a. 故障点三序电流、三相电流有效值：I1 = 0.0000 - 9.8592iIa = 9.8592I2 = 0 Ib = 9.8592 I0 = 0Ic = 9.8592b. 各节点三序电压、三相电压有效值：c. 各支路三序电流、三相电流有效值：2.不对称短路计算过程流程图和计算结果（1）节点3发生A相短路接地故障：a.故障点三序电流、三相电流有效值：I1 = 0.0000 - 3.1167i Ia = 9.3500I2 = 0.0000 - 3.1167i Ib = 0.0000I0 = 0.0000 - 3.1167i Ic = 0.0000b.各节点三序电压、三相电压有效值：（2）节点3发生A相经10Ω电阻接地故障a.故障点三序电流、三相电流有效值：I1 = 0.9793 - 2.7705i Ia = 8.8154I2 = 0.9793 - 2.7705i Ib = 0.0000I0 = 0.9793 - 2.7705i Ic = 0.0000b.各节点三序电压、三相电压有效值：c.各支路三序电流、三相电流有效值：（3）节点3发生b、c两相相间短路故障a.故障点三序电流、三相电流有效值：I1 = 0.0000 - 4.9296i Ia = 0I2 = 0.0000 + 4.9296i Ib = 8.5383I0 = 0 Ic = 8.5383b.各节点三序电压、三相电压有效值：c.各支路三序电流、三相电流有效值：（4） 节点3发生b 、c 两相接地短路故障a. 故障点三序电流、三相电流有效值：I1 = 0.0000 - 6.4114i Ia = 0.0000 I2 = 0.0000 + 3.4478i Ib = 9.6262 I0 = 0.0000 + 2.9636iIc = 9.6262b. 各节点三序电压、三相电压有效值：c.各支路三序电流、三相电流有效值：3.静态稳定性计算过程流程图和计算结果流程图：计算结果：Ke=0.100000, delta=92, P=1.345052 Ke=0.200000, delta=93, P=1.363051Ke=0.300000, delta=94, P=1.380840 Ke=0.400000, delta=96, P=1.398218 Ke=0.500000, delta=97, P=1.415607 Ke=0.600000, delta=98, P=1.432731 Ke=0.700000, delta=99, P=1.449575 Ke=0.800000, delta=100, P=1.466124 Ke=0.900000, delta=101, P=1.482360 Ke=1.000000, delta=102, P=1.498265Ke=1.100000, delta=102, P=1.514127 Ke=1.200000, delta=103, P=1.529419Ke=1.300000, delta=104, P=1.544342 Ke=1.400000, delta=104, P=1.559111Ke=1.500000, delta=105, P=1.573400 Ke=1.600000, delta=106, P=1.587273Ke=1.700000, delta=106, P=1.601053 Ke=1.800000, delta=106, P=1.614344Ke=1.900000, delta=103, P=1.623638 Ke=2.000000, delta=100, P=1.626568Ke=2.100000, delta=98, P=1.62837 Ke=2.200000, delta=95, P=1.620900 Ke=2.300000, delta=93, P=1.615423 Ke=2.400000, delta=90, P=1.598429 Ke=2.500000, delta=88, P=1.586369 Ke=2.600000, delta=86, P=1.571759 Ke=2.700000, delta=84, P=1.554738 Ke=2.800000, delta=82, P=1.535446 Ke=2.900000, delta=81, P=1.527090 Ke=3.000000, delta=79, P=1.504406 Ke=3.100000, delta=77, P=1.479812 Ke=3.200000, delta=76, P=1.468345 Ke=3.300000, delta=74, P=1.440938 Ke=3.400000, delta=73, P=1.427839 Ke=3.500000, delta=71, P=1.397958 Ke=3.600000, delta=70, P=1.383436 Ke=3.700000, delta=69, P=1.368427 Ke=3.800000, delta=67, P=1.335646 Ke=3.900000, delta=66, P=1.319478 Ke=4.000000, delta=65, P=1.302918 Ke=4.100000, delta=64, P=1.285990 Ke=4.200000, delta=63, P=1.268718 Ke=4.300000, delta=62, P=1.251124 Ke=4.400000, delta=61, P=1.233229 Ke=4.500000, delta=60, P=1.215054 Ke=4.600000, delta=59, P=1.196618 Ke=4.700000, delta=58, P=1.177939 Ke=4.800000, delta=57, P=1.159036 Ke=4.900000, delta=56, P=1.139925 Ke=5.000000, delta=55, P=1.120623 Ke=5.100000, delta=54, P=1.101144 Ke=5.200000, delta=53, P=1.081504 Ke=5.300000, delta=52, P=1.061715 Ke=5.400000, delta=52, P=1.061829 Ke=5.500000, delta=51, P=1.041864 Ke=5.600000, delta=50, P=1.021780 Ke=5.700000, delta=49, P=1.001590分析结果：最终选择放大倍数Ke=2.1特性曲线：4.暂态稳定性计算过程流程图和计算结果流程图：摇摆曲线：(1)0.15s时切除故障：(2)0.25s时切除故障：分析结果：根据摇摆曲线判断：0.15s时切除故障系统暂态稳定，0.25s时切除故障系统失稳。