电力系统暂态上机计算课程设计报告
电力系统暂态的课程设计
电力系统暂态的课程设计一、课程目标知识目标：1. 学生能够理解电力系统暂态过程的基本概念，掌握暂态过程中的关键参数及其相互影响。
2. 学生能够描述常见电力系统暂态现象，如短路、故障、暂态稳定性问题等，并分析其原因及危害。
3. 学生能够掌握电力系统暂态稳定性分析的基本方法，并运用相关公式进行简单计算。
技能目标：1. 学生能够运用所学知识，分析实际电力系统暂态问题，提出解决策略。
2. 学生能够通过查阅资料，了解电力系统暂态过程的最新研究动态和发展趋势。
3. 学生能够利用仿真软件对电力系统暂态过程进行模拟，验证所学理论和方法。
情感态度价值观目标：1. 培养学生对电力系统暂态问题的关注和兴趣，激发其探索精神。
2. 增强学生的团队合作意识，使其在分析和解决电力系统暂态问题时，能够积极与他人交流与合作。
3. 培养学生严谨、务实的科学态度，使其在分析问题时，能够充分考虑到实际工况，注重实际应用。
课程性质：本课程为电力系统专业课程，旨在帮助学生深入了解电力系统暂态过程，掌握相关分析方法和解决策略。
学生特点：学生具备一定的电力系统基础知识，具有较强的学习能力和动手能力。
教学要求：结合学生特点，注重理论与实践相结合，提高学生的实际操作能力和解决问题的能力。
在教学过程中，注重启发式教学，引导学生主动思考，培养其创新意识。
二、教学内容1. 电力系统暂态过程基本概念：包括暂态过程的定义、分类及其在电力系统中的重要性。
- 教材章节：第二章 电力系统暂态过程概述- 内容列举：暂态过程的产生原因、暂态过程的基本类型及其特点。
2. 常见电力系统暂态现象及分析：如短路、故障、暂态稳定性问题等。
- 教材章节：第三章 常见电力系统暂态现象- 内容列举：短路故障的分类及影响、暂态稳定性问题的分析方法及其防治措施。
3. 电力系统暂态稳定性分析：介绍暂态稳定性分析的基本原理、方法和应用。
- 教材章节：第四章 电力系统暂态稳定性分析- 内容列举：暂态稳定性分析的基本理论、时域仿真法、直接法等分析方法。
电力系统暂态课程设计
电力系统暂态课程设计一、课程目标知识目标：1. 掌握电力系统暂态过程的基本概念、分类及特点；2. 理解电力系统暂态稳定性分析的基本原理和方法；3. 了解电力系统暂态过程中的主要影响因素及其作用机理。
技能目标：1. 能够运用所学知识对电力系统暂态过程进行正确分析；2. 能够运用暂态稳定性分析的方法，评估电力系统的稳定性；3. 能够根据实际案例，提出改善电力系统暂态稳定性的措施。
情感态度价值观目标：1. 培养学生对电力系统暂态过程研究的兴趣，激发学生探索精神；2. 增强学生的团队协作意识，培养合作解决问题的能力；3. 提高学生关注电力系统安全、稳定的意识，树立社会责任感。
本课程针对高年级学生，结合电力系统暂态过程的学科特点，注重理论与实践相结合。
通过本课程的学习，使学生能够深入理解电力系统暂态过程的本质，掌握相关分析方法和技能，为今后从事电力系统运行、维护和管理等工作奠定基础。
同时，培养学生关注社会热点问题，提高解决实际问题的能力。
课程目标的设定旨在使学生在知识、技能和情感态度价值观方面取得全面发展。
二、教学内容1. 电力系统暂态过程基本概念：暂态过程的定义、分类及其在电力系统中的重要性；教材章节：第一章第一节。
2. 电力系统暂态稳定性分析原理：暂态稳定性分析的基本理论、方法及其应用；教材章节：第二章。
3. 电力系统暂态过程中的影响因素：分析电力系统暂态过程中的主要因素，如发电机、变压器、线路等；教材章节：第三章。
4. 暂态稳定性分析的方法与案例：介绍暂态稳定性分析的方法，结合实际案例进行分析；教材章节：第四章。
5. 改善电力系统暂态稳定性的措施：探讨提高暂态稳定性的措施，如励磁系统、PSS、直流输电等；教材章节：第五章。
6. 电力系统暂态过程实验与仿真：组织学生进行暂态过程实验，利用仿真软件进行模拟分析；教材章节：第六章。
教学内容按照课程目标进行科学、系统地组织和安排，注重理论与实践相结合。
教学进度根据教材章节进行制定，确保学生能够在掌握基本理论和方法的基础上，通过实际案例和实验，深入理解电力系统暂态过程及其稳定性问题。
(完整word版)电力系统暂态分析课程设计
电力系统的短路计算能够帮助我们避免很多不必要的损失，随着科学技术的发展，电力已经和人们的生活密切相关，而建立结构合理的大型电力系统不仅便于电能生产与消费的集中管理、统一调度和分配，减少总装机容量，节省动力设施投资，且有利于地区能源资源的合理开发利用，更大限度地满足地区国民经济日益增长的用电需要。
电力系统的规模和技术水准已成为一个国家经济发展水平的标志之一。
电力系统的短路故障是严重的，而又是发生几率最多的故障，当发生短路时，其短路电流可达数万安以至十几万安，它们所产生的热效应和电动力效应将使电气设备遭受严重破环。
为使非故障部分从不正常运行情况下解脱出来，这要求电气设备必须有足够的机械强度和热稳定度，开关电气设备必须具备足够的开断能力。
因此，电力系统短路电流计算是电力系统运行分析，设计计算的重要环节，许多电业设计单位和个人倾注极大精力从事这一工作的研究。
由于电力系统结构复杂，随着生产发展，技术进步系统日趋扩大和复杂化，短路电流计算工作量也随之增大，采用计算机辅助计算势在并行。
我们这次课程设计使用PSCAD软件对电力系统进行线路建模，从而计算出不同短路类型时的短路电流及其波形。
关键字：短路计算；pscad软件使用；电力系统建模。
目录摘要1 课题概述 (1)1.1课题要求 (1)1.2课题内容 (5)1.3 课题目的 (6)2短路故障分析 (6)2.1不对称故障的分析 (6)2.2三相短路故障分析 (7)2.3短路电流计算步骤 (9)3方案设计 (10)3.1方案概述 (10)3.2 课题设计图 (10)3.3 课题步骤图 (11)4两相短路的仿真分析 (12)4.1 PSCAD简介 (12)4.2 两相短路故障的仿真 (13)5总结 (15)参考文献 (16)1课题概述1.1 课题要求（1）通过课程设计是学生掌握电力系统三相短路计算的基本原理和方法；（2）掌握并能熟练运用PSCAD/MATLAB仿真软件；（3）建立系统三相接线图的仿真过程；（5）编写短路计算流程图；（4）得出仿真结果。
电力系统暂态分析课程设计
内蒙古科技大学电力系统稳态分析课程设计题目：基于MATLAB的电力系统复杂潮流分析学生姓名：刘建峰学号：1167130207专业：电气工程及其自动化班级：电气2011—2班指导教师：刘景霞摘要电力系统潮流计算是电网分析的基础应用，是对复杂电力系统正常和故障条件下稳态运行状态的计算。
给定电力系统的网络结构、参数和决定电力系统运行状况的边界条件，确定电力系统运行的方法之一是朝流计算。
MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算．采用迭代法，通过建立矩阵的修正方程来依次迭代，逐步逼近真值来计算出电力网的电压，功率分布。
PQ分解法是极坐标形式牛顿-拉潮流计算的一种简化计算方法，。
P—Q分解法通过对电力系统具体特点的分析，对牛顿法修正方程式的雅可比矩阵进行了有效的简化和改进。
由于这些简化只涉及修正方程式的系数矩阵，并未改变节点功率平衡方程和收敛判据，因不会降低计算结果的精度。
用手算和计算机算法对其进行设计。
使用MATLAB软件进行编程，在很大程度上节省了内存，减少了计算量。
通过对本题计算我们了解了一些工程计算和解决工程问题的方法。
关键词：潮流计算，PQ分解法，MATLABElectrical power system complex tidal currentanalysis based on MATLABPower Flow Analysis Grid computing is the basis of applications, the complex power system under normal and fault conditions for the calculation of steady state operation. Given the power system network structure, parameters and decisions operation of the power system boundary conditions, to determine the method of operation of the power system is one of North Korea flow calculation.MATLAB is an interactive, object-oriented programming language, widely used in industry and academia, mainly for matrix calculation. Using iteration, the amendment through the establishment of matrix iterative equation to turn, gradually moving towards a true value to calculate the voltage electricity grid, power distribution.PQ decomposition method is the form of polar coordinates Newton - the widening trend of a simplified calculation method. P-Q decomposition method adopted by the specific characteristics of the power system analysis, Newton’s Law of the Jacobian matrix formula has effectively simplified and improved. As a result of these simplified formula that involves only the coefficient matrix, the balance of power has not changed node equations and the convergence criterion, because the results will not reduce the accuracy.Use MATLAB software programming, saving memory to a large extent, reduce the amount of computation. By this calculation we understand that a number of engineering calculation and solve engineering problemsKeywords : The trend, the PQ decomposition method, MATLAB目录内蒙古科技大学课程设计任务书 (5)第一章引言 (6)1.1研究背景及意义 (7)1.2潮流计算的意义 (8)1.3电力系统稳态分析潮流计算总结 (9)1.4MATLAB的概述 (9)第二章理论计算 (10)2.1P-Q法潮流计算的基本步骤 (10)2.2PQ分解法潮流计算流程图 (11)2.3两机五节点网络潮流计算 (12)第三章程序设计 (18)3.1设计流程 (18)3．2程序设计 (18)3.3程序运行结果 (22)第四章设计感想 (45)参考文献 (46)内蒙古科技大学课程设计任务书12 第一章 引言1.1 研究背景及意义电力系统在运行时，在电源电势激励作用下，电流或功率从电源通过系统各元件流入负荷，分布于电力网各处，称为潮流分布。
新版电力系统暂态上机计算课程设计报告完整无误版附程序
课程设计报告( 2014—2015年度第二学期)名称：电力系统暂态上机计算院系：电气与电子工程学院班级：电气1211学号：1121210205学生姓名：郝阳指导教师：陈艳波设计周数：两周成绩：日期：2015年7月4日一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。
二、设计正文（详细内容见附录,用A4纸，页数不限）1．对称短路计算过程流程图和计算结果2．不对称短路计算过程流程图和计算结果3．静态稳定性计算过程流程图和计算结果4．暂态稳定性计算过程流程图和计算结果5．思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。
通过本次课程设计，本人对电力系统故障分析有了更深刻的理解，电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。
在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。
在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。
同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。
这与手算系统短路电流时使用的网络化简方法大大不同。
在学习中，参照潮流程序，我加深了对节点导纳矩阵建立方法的理解与学习，巩固了不同类型短路的短路电流计算方法，和序电压、序电流，相电压、相电流基于matlab软件的计算。
利用小干扰分析法判别静态稳定可以使用劳斯判据，也可以使用特征根判别的方法。
二者都能够判断，在书中的例题中使用的是劳斯判据，但是，特征根判别放法更适合在matlab软件的环境下使用，因为该软件提供了计算矩阵特征根的函数，使用起来简便易懂。
电力系统暂态上机计算课程设计报告
-0.5489+ 0.3169i
5
0.6091
0
0
0.5977+ 0.3451i
0.0000 - 0.6901i
-0.5977+ 0.3451i
c.各支路三序电流、三相电流：
支路
I1
I2
I0
Ia
Ib
Ic
1-2
0 + 0.8451i
0
0
0 + 0.8451i
0.7319- 0.4225i
3.《电力系统故障的计算机辅助分析》 大学 米麟书等
4.《电力系统潮流计算》 大学 宋文南等
5.《电力系统故障分析》 清华大学 周荣光
6.《短路电流实用计算方法》 电力工业 交通大学等
7.《精通Matlab6.5》航空航天大学，志涌
附录（设计流程图、计算结果、思考题答案）
1、对称短路计算流程图及计算结果
-0.7319- 0.4225i
1-3
0 - 4.5070i
0
0
0 - 4.5070i
-3.9032+ 2.2535i
3.9032+ 2.2535i
2-3
0 - 5.3521i
0
0
0 - 5.3521i
-4.6351+ 2.6761i
4.6351+ 2.6761i
1-4
0 + 3.6620i
0
0
0+3.6620i
另外，通过此次课程设计，我学习并掌握了一定的Matlab知识，对简单的选择语句、循环语句、数组的使用等有了一些了解，为今后学习、工作中再次学习使用此软件打下了基础，获益匪浅。
电力系统暂态上机计算课程设计报告
课程设计报告( 2010—2011年度第二学期)名称：电力系统暂态上机院系：电气与电子工程学院班级：电气0803班日期：2011年7月6一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。
二、主要内容1．电磁暂态过程计算1.1 三相短路计算1.1.1 计算机计算三相短路流程图（程序见附录）附：式3-17 I f = U f|0| / ( Z ff + z f )≈1/( Z ff + z f ) ≈1/Z ff式3-19 U1= U1|0| +△U1 = U1|0| – Z1f I f ≈1 – Z1f I fU f = U f|0| +△U f = Z f I fU n= U n|0|+△U n = U n|0| – Z nf I f ≈1 – Z nf I f式3-20 I ij = (U i – U j) / z ij≈(△U i –△U j ) / z ij = (△U i –△U j) y ij1.1.2 计算教材《电力系统暂态分析(第三版)》P68例（3-2）例3-2：下图所示一环形网络，已知各元件参数为：发电机：G1～G3 100MW,10.5KV,cosφN=0.86，xd”=0.183。
变压器：T1～T3 120MVA,115/10.5KV,Us(%)=10.5。
线路：三条线路完全相同，长50km,电抗0.44欧姆/㎞试计算母线③三相短路后时刻的各节点电压、各支路电流以及发电机端电压。
电路图和等值电路图如下：故障分量网络：1.1.3节点发生三相短路接地计算结果得：（下面各数值均为标幺值）各节点的三相电压各支路电流发电机的端电压1.2.不对称短路计算过程流程图和计算结果1.2.1不对称短路计算过程流程图：附：式（5-39）U i(1)=U i|0|-Z if(1)I f(1)U i(2)=-Z if(2)I f(2)U i(0)=-Z if(0)I f(0)式（5-40）I ij(1) = ( U i(1) - U j(1) ) / z ij(1)I ij(2) = ( U i(2) - U j(2) ) / z ij(2)I ij(0) = ( U i(0) - U j(0) ) / z ij(0)式（5-58）U i(1)=U i|0|-(Z iq(1)-Z ik(1))I(1)U i(2)=-(Z iq(2)-Z ik(2))I(2)U i(0)=-(Zi q(0)-Z ik(0))I(0)1.2.2计算教材《电力系统暂态分析(第三版)》P127例（5-1）例5-1：在例3-2的系统中又已知三台发电机中性点均不接地；三台变压器均为YNd11接线（发电机侧为三角形）；经试验得三条输电线路的零序电抗均为0.20（以60MVA为基准值）。
华电潮流上机课程设计报告(程序最终无误版)
《电力系统潮流上机》课程设计报告院系：电气与电子工程学院班级：电气1108学号： 1111550112学生：龙日尚指导教师：宝柱设计周数：两周成绩：日期：2014年1月10日一、课程设计的目的与要求培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识二、设计正文（详细容见附录）1．手算2．计算机计算3．思考题三、课程设计总结或结论（详细容见附录）四、参考文献1.《电力系统计算：电子数字计算机的应用》，交通大学等合编。
北京：水利电力；2.《现代电力系统分析》，王锡凡主编，科学；3.《电力系统稳态分析》，珩，中国电力，1995年，第三版；附录（设计流程图、程序、表格、数据等）4.机算潮流程序及结果// dierti.cpp : 定义控制台应用程序的入口点。
//#include "stdafx.h"struct Line //线路结构体{int Num,NumI,NumJ; //线路号左节点名右节点名float R,X,B,K; //电阻电抗电纳变比（K等于1为普通支路，不等于1为变压器支路的变比） };struct Bus //节点结构体{int Num ;float Volt,Phase,GenP,GenQ,LoadP,LoadQ;int Type;};#include"stdio.h"#include"string.h"#include"math.h"#include"stdlib.h"#define NBUS 4#define NLINE 4/* Global variables */int nL,nB,nVA,nSH;float X[NBUS];int L;double def[2*NBUS];double mn[50];void Gauss(double a[50][50],double b[50], int n) /*定义高斯法 */{int JS[50];int i,j,k;float d,t,x[50];FILE *fp;int L=1;for(i=0;i<50;i++) JS[i]=0;for(k=0;k<n;k++){d=0.0;for(j=k;j<n;j++)if(fabs(a[k][j])>d){ /*在一行中找到一个最大值赋值d,并用JS[K]记住这个最大值所在的列号*/ d=fabs(a[k][j]);JS[k]=j;}if(fabs(d)<0.000001) /*如果d的数值太小，做为被除数将带来很大的误差 */L=0;else {if(JS[k]!=k)for(i=0;i<n;i++){t=a[i][k];a[i][k]=a[i][JS[k]]; /*进行列交换，让最大值始终在对角元上*/a[i][JS[k]]=t;}}if(L==0)break;for(j=k+1;j<n;j++)a[k][j]=a[k][j]/a[k][k]; /*对角元上的元素消为1*/b[k]=b[k]/a[k][k];for(i=k+1;i<n;i++){for(j=k+1;j<n;j++)a[i][j]=a[i][j]-a[i][k]*a[k][j]; /*使下三角阵的元素为0*/b[i]=b[i]-a[i][k]*b[k];}}if(fabs(a[n-1][n-1])>0.00001){ /*用追赶法，解方程组，求未知数x*/ x[n-1]=b[n-1];for(i=n-2;i>=0;i--){t=0.0;for(j=i+1;j<n;j++)t=t+a[i][j]*x[j];x[i]=(b[i]-t);}}if((fp=fopen("gauss.txt","w"))==NULL) /*将结果写到TXT文件中*/{printf("err");exit(0);}for(i=0;i<n;i++){fprintf(fp,"%lf",x[i]);mn[i]=x[i];fprintf(fp,"\n");}fclose(fp);if(fp!=NULL) fclose(fp);}int _tmain(int argc, _TCHAR* argv[]){FILE *fp;FILE *fpout;int i,j,k,l,h,n,v;int i1,i2,i3,kp,kq;float d1,d2,d3,d4,d5,d6,r,x,g,b,tt,LL,e,ps,qs,shsh,m;struct Line sL[NLINE];struct Bus sB[NBUS];float YG[NBUS+1][NBUS+1],YB[NBUS+1][NBUS+1];double u[50][2];i1=i2=i3=0;d1=d2=d3=d4=d5=d6=ps=qs=0.0;for(i=0;i<NBUS;i++)if((fp=fopen("in.txt","r"))==NULL){ printf("Can not open the file named 'in.txt' \n");exit(0);}fscanf(fp,"%d,%d,%d",&nB,&nL,&nSH);for(i=0;i<nB;i++){sB[i].Num=sB[i].Type=0;sB[i].Volt=1.0;sB[i].Phase=sB[i].GenP=sB[i].GenQ=sB[i].LoadP=sB[i].LoadQ=0.0;fscanf(fp,"%d,%f,%f,%f,%f,%f,%f,%d",&i1,&d1,&d2,&d3,&d4,&d5,&d6,&i2);sB[i].Num=i1;sB[i].Volt=d1;sB[i].Phase=d2;sB[i].GenP=d3;sB[i].GenQ=d4;sB[i].LoadP=d5,sB[i].LoadQ=d6;sB[i].T ype=i2;};for(i=0;i<nL;i++){sL[i].Num=sL[i].NumI=sL[i].NumJ=0;sL[i].R=sL[i].X=sL[i].B=0.0;sL[i].K=1.0;fscanf(fp,"%2d %3d %3d %f %f %f %f",&i1,&i2,&i3,&d1,&d2,&d3,&d4);sL[i].Num=i1;sL[i].NumI=i2;sL[i].NumJ=i3;sL[i].R=d1;sL[i].X=d2;sL[i].B=d3;sL[i].K=d4;}if(fp!=NULL) fclose(fp);/*Make Y Matrix*/for(i=1;i<nB+1;i++)for(j=1;j<nB+1;j++){YG[i][j]=0.0;YB[i][j]=0.0;};for(l=0; l<nL; l++){i=sL[l].NumI;j=sL[l].NumJ;r=sL[l].R;x=sL[l].X;d1=r*r+x*x;g=r/d1;b=-x/d1;m=sL[l].K;if(fabs(sL[l].K-1.0)<0.000001) //普通支路 {YG[i][i]=YG[i][i]+g;YG[j][j]=YG[j][j]+g;YB[i][i]=YB[i][i]+b+sL[l].B;YB[j][j]=YB[j][j]+b+sL[l].B;YG[i][j]=YG[i][j]-g;YG[j][i]=YG[j][i]-g;YB[i][j]=YB[i][j]-b;YB[j][i]=YB[j][i]-b;}else //变压器支路{YG[i][i]=YG[i][i]+g/m+g*(m-1)/m;YG[j][j]=YG[j][j]+g/m+g*(1-m)/m/m;YB[i][i]=YB[i][i]+b/m+b*(m-1)/m;YB[j][j]=YB[j][j]+b/m+b*(1-m)/m/m;YG[i][j]=YG[i][j]-g/m;YG[j][i]=YG[j][i]-g/m;YB[i][j]=YB[i][j]-b/m;YB[j][i]=YB[j][i]-b/m; }}/* Check the Y matrix */if((fp=fopen("GGBB.txt","w"))==NULL){printf("Can not open the file named 'GGBB.txt' \n");exit(0);}fprintf(fp,"---Y Matrix---\n");for(i=1;i<nB+1;i++)for(j=1;j<nB+1;j++)if(fabs(YB[i][j]-0.0)>0.000001) fprintf(fp,"Y(%3d,%-3d)=(%10.5f,%10.5f)\n",i,j,YG[i][j],YB[i][j]);if(fp!=NULL) fclose(fp);/* 节点电压附初值 */for(i=1;i<nB+1;i++){if(sB[i-1].Type==0){u[i][0]=0.0;u[i][1]=1.0;}else if(sB[i-1].Type==1){u[i][1]=sB[i-1].Volt;u[i][0]=0.0;}else if(sB[i-1].Type==2){u[i][1]=sB[i-1].Volt;u[i][0]= sB[i-1].Phase;}for(v=1;;v++)/* 迭代次数可以无限大 */{/* 节点电压附初值 */printf("迭代第%d次赋予的电压初值为e+jf:\n",v); for(i=1;i<nB+1;i++)printf("%lf,%lf\n",u[i][1],u[i][0]);printf("\n");printf("\n");/* 求偏移量 */double P_P[10];double P_Q[10];double P_UU[10];for(i=1;i<nB+1;i++){if(sB[i-1].Type==2){P_P[i]=0.0;P_Q[i]=0.0;P_UU[i]=1.05;}if(sB[i-1].Type==0){double tempP=0.0;double tempQ=0.0;for(j=1;j<nB+1;j++){tempP+=YG[i][j]*u[j][1]-YB[i][j]*u[j][0]; tempQ+=YG[i][j]*u[j][0]+YB[i][j]*u[j][1];P_P[i]=(sB[i-1].GenP-sB[i-1].LoadP)-tempP*u[i][1]-tempQ*u[i][0]; P_Q[i]=(sB[i-1].GenQ-sB[i-1].LoadQ)-tempP*u[i][0]+tempQ*u[i][1]; P_UU[i]=0.0;}if(sB[i-1].Type==1){double tempP=0.0;double tempQ=0.0;for(j=1;j<nB+1;j++){tempP+=YG[i][j]*u[j][1]-YB[i][j]*u[j][0];tempQ+=YG[i][j]*u[j][0]+YB[i][j]*u[j][1];P_P[i]=(sB[i-1].GenP-sB[i-1].LoadP)-tempP*u[i][1]-tempQ*u[i][0]; }P_UU[i]=sB[i-1].Volt*sB[i-1].Volt-u[i][1]*u[i][1]-u[i][0]*u[i][0]; P_Q[i]=0.0;}}/* 偏移量阵 */double P_PQ[6];int a=0;for(i=1;i<3;i++){P_PQ[a]=P_P[i];a=a+2;}a=1;for(i=1;i<3;i++){P_PQ[a]=P_Q[i];a=a+2;}P_PQ[4]=P_P[3];P_PQ[5]=P_UU[3];printf("迭代第%d次的偏移量为：\n",v);for(i=0;i<6;i++){printf("%f",P_PQ[i]);printf("\n");}printf("\n");printf("\n");/* 雅可比矩阵 */double H[6][6],N[6][6],J[6][6],L[6][6],R[6][6],S[6][6],aa[6],bb[6]; for(i=1;i<5;i++){ if(fabs(sB[i-1].Type-2.0)<0.000001)continue;else{for(j=1;j<5;j++)if(i!=j){H[i][j]=-YB[i][j]*u[i][1]+YG[i][j]*u[i][0];N[i][j]=YG[i][j]*u[i][1]+YB[i][j]*u[i][0];J[i][j]=-N[i][j];L[i][j]=H[i][j];R[i][j]=0;S[i][j]=0;}else{aa[i]=bb[i]=0.0;for(n=1;n<5;n++)aa[i]+=YG[i][n]*u[n][1]-YB[i][n]*u[n][0];bb[i]+=YG[i][n]*u[n][0]+YB[i][n]*u[n][1];}H[i][i]=-YB[i][i]*u[i][1]+YG[i][i]*u[i][0]+bb[i]; N[i][i]=YG[i][i]*u[i][1]+YB[i][i]*u[i][0]+aa[i]; J[i][i]=-YG[i][i]*u[i][1]-YB[i][i]*u[i][0]+aa[i]; L[i][i]=YG[i][i]*u[i][0]-YB[i][i]*u[i][1]-bb[i]; R[i][i]=2*u[i][0];S[i][i]=2*u[i][1];}}}double ss[50][50];for(i=0;i<6;i++)for(j=0;j<6;j++)ss[i][j]=0.0;for(i=1;i<3;i++)for(j=1;j<4;j++){ss[2*i-2][2*j-2]=H[i][j];ss[2*i-2][2*j-1]=N[i][j];ss[2*i-1][2*j-2]=J[i][j];ss[2*i-1][2*j-1]=L[i][j];}i=3;for(j=1;j<4;j++){ss[2*i-2][2*j-2]=H[i][j];ss[2*i-2][2*j-1]=N[i][j];ss[2*i-1][2*j-2]=R[i][j];ss[2*i-1][2*j-1]=S[i][j];}printf("迭代第%d次的雅可比矩阵为：\n",v);for(i=0;i<6;i++){for(j=0;j<6;j++)printf("%10f",ss[i][j]);printf("\n");}printf("\n");printf("\n");Gauss(ss,P_PQ,6);for(i=1;i<nB;i++){u[i][0]=u[i][0]+mn[2*(i-1)];u[i][1]=u[i][1]+mn[2*i-1];}double max;max=fabs(P_PQ[0]);for(i=0;i<=5;i++)if (max<fabs(P_PQ[i]))max=fabs(P_PQ[i]);if(fabs(max)<0.0001){printf("满足精度要求，迭代终止，迭代次数为%d\n",v); printf("\n");printf("\n");}/* 叠代循环的括号 */printf("最终求得的节点电压值为e+jf:\n");for(i=1;i<nB+1;i++)printf("%lf,%lf\n",u[i][1],u[i][0]);printf("\n");printf("\n");double uu[5],Phase[5];for(i=1;i<nB+1;i++){uu[i]=sqrt(u[i][1]*u[i][1]+u[i][0]*u[i][0]); Phase[i]=atan(u[i][0]/u[i][1]);}for(i=1;i<nB+1;i++)printf("%lf,%lf\n",uu[i],Phase[i]);*计算线路功率和平衡节点 PV节点功率*/double P[5],Q[5];double tempP=0.0;double tempQ=0.0;for(i=1;i<nB+1;i++){for(j=1;j<nB+1;j++){tempP+=YG[i][j]*u[j][1]-YB[i][j]*u[j][0]; tempQ+=YG[i][j]*u[j][0]+YB[i][j]*u[j][1];}P[i]=tempP*u[i][1]+tempQ*u[i][0];Q[i]=tempP*u[i][0]-tempQ*u[i][1];tempP=0.0;for(i=1;i<nB+1;i++)printf("节点%d注入功率为%lf,%lf\n",i,P[i],Q[i]);/* 支路功率 */double V[4][2];for(i=1;i<5;i++)for(j=0;j<3;j++)V[i][j]=u[i][j];double sP[5][5],sQ[5][5];double dsq,dsp,dp,sumgen;for(i=1;i<NBUS+1;i++){for(j=1;j<NBUS+1;j++){sP[i][j]=0.0;sQ[i][j]=0.0;}}for(l=0; l<nL; l++){i=sL[l].NumI;j=sL[l].NumJ;r=sL[l].R;x=sL[l].X;d1=r*r+x*x;g=r/d1;if(fabs(sL[l].K-1.0)<0.000001){/*Normal lines or transformers*/sP[i][j]=V[i][1]*V[i][1]*g-V[i][1]*V[j][1]*(g*cos(V[i][0]-V[j][0])+b*sin(V[i][0]-V[j][0]));sQ[i][j]=-(V[i][1]*V[i][1]*sL[l].B+V[i][1]*V[i][1]*b+V[i][1]*V[j][1]*(g*sin(V[i][0]-V[j][0])-b*cos(V[i ][0]-V[j][0])));sP[j][i]=V[j][1]*V[j][1]*g-V[i][1]*V[j][1]*(g*cos(V[j][0]-V[i][0])+b*sin(V[j][0]-V[i][0]));sQ[j][i]=-(V[j][1]*V[j][1]*sL[l].B+V[j][1]*V[j][1]*b+V[i][1]*V[j][1]*(g*sin(V[j][0]-V[i][0])-b*cos(V[j ][0]-V[i][0])));}else{/*abnormal transformer ratio*/sP[i][j]=V[i][1]*V[i][1]*g/sL[l].B/sL[l].B-V[i][1]*V[j][1]*(g*cos(V[i][0]-V[j][0])/sL[l].B+b*sin(V[i][ 0]-V[j][0])/sL[l].B);sQ[i][j]=-(V[i][1]*V[i][1]*b/sL[l].B/sL[l].B+V[i][1]*V[j][1]*(g*sin(V[i][0]-V[j][0])/sL[l].B-b*cos(V[i ][0]-V[j][0])/sL[l].B));sP[j][i]=V[j][1]*V[j][1]*g-V[i][1]*V[j][1]*(g*cos(V[j][0]-V[i][0])/sL[l].B+b*sin(V[j][0]-V[i][0])/sL[l ].B);sQ[j][i]=-(V[i][1]*V[i][1]*b+V[i][1]*V[j][1]*(g*sin(V[j][0]-V[i][0])/sL[l].B-b*cos(V[j][0]-V[i][0])/sL [l].B));}}/* 输电效率 */dsp=P[4];dsq=Q[4];for(i=0;i<NBUS;i++){dsp+=sB[i].GenP-sB[i].LoadP;dsq+=sB[i].GenQ-sB[i].LoadQ;sumgen+=sB[i].GenP;}dp=dsp/sumgen*100;/* 输出功率情况 */if((fp=fopen("功率情况.txt","w"))==NULL){printf("Can not open the file named '功率情况.txt' \n");exit(0);}fprintf(fp,"---功率情况---\n");fprintf(fp,"平衡节点功率S=%10.5f+ j%10.5f\n",P[4],Q[4]);for(i=1;i<NBUS+1;i++)for(j=1;j<NBUS+1;j++)if(fabs(sP[i][j]-0.0)>0.000001)fprintf(fp,"S(%3d,%-3d)=(%10.5f,j%10.5f)\n",i,j,sP[i][j],sQ[i][j]); fprintf(fp,"网损为%10.5f+j%10.3f,输电效率为%10.3f\n",dsp,dsq,100-dp);if(fp!=NULL) fclose(fp);return 0;}1．导纳阵Y( 1,1 )=( 1.01534, -8.19201) Y( 1,2 )=( -0.56148, 2.30208) Y( 1,3 )=( 0.00000, 3.66667) Y( 1,4 )=( -0.45386, 1.89107) Y( 2,1 )=( -0.56148, 2.30208) Y( 2,2 )=( 1.04225, -4.67651) Y( 2,4 )=( -0.48077, 2.40385) Y( 3,1 )=( 0.00000, 3.66667) Y( 3,3 )=( 0.00000, -3.33333) Y( 4,1 )=( -0.45386, 1.89107) Y( 4,2 )=( -0.48077, 2.40385) Y( 4,4 )=( 0.93463, -4.26159)2．设定电压初值01.1;01;01)0(3)0(3)0(2)0(2)0(1)0(1j jf e j jf e j jf e +=++=++=+ 3．计算功率和电压偏移；27731.0])()([41)0(11)0(1)0(141)0(11)0(1)0(11)0(11)0(1-=++--=-=∆∑∑==jj j j jj jj s s e B f G ff B e G e P P P P0.05097])()([41)0(11)0(1)0(141)0(11)0(1)0(11)0(11)0(1-=----=-=∆∑∑==j j j j j j jj s s e B f G ef B eG f Q Q Q Q同理可算出52596.0)0(22)0(2-=-=∆P P P s ,0196.0)0(22)0(2=-=∆Q Q Q s5.0)0(33)0(3=-=∆P P P s ,0.02)0(3232)0(3=-=∆U U U s4．根据求的第一次迭代时雅可比矩阵各元素的公式计算雅可比矩阵各个元素的具体值：⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤-----------⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡20000.200000066667.300003334.40052691.406629.130208.256148.00001821.182612.456148.030208.266667.3030208.256148.006298.803803.1066667.356148.030208.299265.032104.8 5．求高斯计算后的修正量：⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆-0.0000000.1276520.023566-0.108546-0.006511-0.007919-2)0(3)0(3)0(2)0(2)0(1)0(11)0()0(3)0(3)0(2)0(2)0(1)0(1U P Q P Q P e f e f e f J 6．计算各节点电压的一次近似值：12765.010855.000792.010000.197643.099349.0)0(3)0(3)1(3)0(2)0(2)1(2)0(1)0(1)1(1)0(3)0(3)1(3)0(2)0(2)1(2)0(1)0(1)1(1=∆+=-=∆+=-=∆+==∆+==∆+==∆+=f f f f f f fffe e e e e e e e e返回第三步重新迭代，并校验收敛与否，令410-=ε。
电力系统暂态上机计算报告
电力系统暂态上机计算报告系别：班级：姓名：FREEMAN学号：2011年1月三、给定题目的程序计算结果：1．电磁暂态输入文件：input.txt2 100,10.5,0.18 120,0.11,1 200,10.5,0.14 200,0.10,231,2,110,0.41,3,60,0.42,3,60,0.42,3,40100115 输出文件：output.txt------------------------------------原始数据-----------------------------------发电机G1: SGN1= 100 MV A, UGN1=10.50 KV , xd1=0.180变压器T1: STN1= 120 MV A, 变比为10.50 / 115.00 KV , Uk1%=11.00000 ,所连节点号为1电源支路数发电机的SN 、发电机的UN 、发电机的次暂态电抗 发电机相连的变压器的SN 、变压器的电抗百分数、发电机支路所接节点 节点总数 支路所接的两个节点号,线路长度,支路单位电抗 故障所在支路的两节点号，距第一个点的距离百分数 功率基准值 线路电压的基准值发电机G2: SGN2= 200 MV A, UGN2=10.50 KV, xd2=0.140变压器T2: STN2= 200 MV A, 变比为10.50 / 115.00 KV, Uk2%=10.00000 ,所连节点号为2输电线路1-2: 长110 KM, 单位电抗0.400 Ω/km输电线路1-3: 长60 KM, 单位电抗0.400 Ω/km输电线路2-3: 长60 KM, 单位电抗0.400 Ω/km故障情况：在线路2-3上距2节点40%处发生三相短路---------------------------------------------------------------------------------------------------------------本程序计算结果为:--------------------------------计算后线路电抗为：线路1-2 :0.332703, 线路1-3 :0.181474, 线路2-4 :0.072590, 线路3-4 :0.108885,节点导纳矩阵为：j(-12.197080) j( 3.005682) j( 5.510417) j( 0.000000)j( 3.005682) j(-25.115057) j( 0.000000) j(13.776042)j( 5.510417) j( 0.000000) j(-14.694444) j( 9.184028) j( 0.000000) j(13.776042) j( 9.184028) j(-22.960069)求得的R和D因子为:j(0.081987) j(-0.246426) j(-0.451782) j(0.000000) 0.00000 j(0.041027) j(-0.055711) j(-0.565185) 0.00000 0.00000 j(0.082445) j(-0.820452) 0.00000 0.00000 0.00000 j(0.142667)所求得电压、电流结果为：短路电流:相量值为if =j( -7.009324)有名值为If = 3.518987 KA各节点电压为:相量值为u1 = 0.478795, 有名值为U1 = 55.061387 KV相量值为u2 = 0.389107, 有名值为U2 = 44.747275 KV相量值为u3 = 0.179548, 有名值为U3 = 20.648020 KV各线路故障电流为:相量值i1-2 = j(-0.269573) , 有名值I1-2 = 0.135338 KA相量值i1-3 = j(-1.648974) , 有名值I1-3 = 0.827857 KA相量值i2-f = j(-5.360351) , 有名值I2-f = 2.691130 KA相量值i3-f = j(-1.648974) , 有名值I3-f = 0.827857 KA相量值△iG1 = j(-1.918547) , 有名值△IG1 = 10.549274 KA, 相量值△iG2 = j(-5.090777) , 有名值△IG2 = 27.992015 KA, -----------------------------------END----------------------------------------2．机电暂态a.输入文件：input.txt200,115240,10.5,0.80,1.1,0.9,0.3,0.44,1.54300,242,10.5,12.0280,220,121,14.0230,0.4,0.4,1.2250,0.98,1152,0.65,0容量基准值，电压基准值发电机容量，电压，功率因数，直轴电抗，交轴电抗，暂态电抗，负序电抗，惯性时间变压器T1容量，高压侧电压，低压侧电压，Uk%变压器T2容量，高压侧电压，低压侧电压，Uk%输电线路长度，单位长度正序、负序、零序电抗正常情况下向无穷大系统输出的功率、功率因数、无穷大系统电压故障类型，故障点距线路始端百分数，接地电抗b.输出文件：output.txt******************* 原始数据****************容量基准值:200.000000 电压基准值:115.000000----------------------------------------------------------------------------------- 发电机参数：额定有功功率:240.000000 额定电压:10.500000 功率因数:0.800000 直轴电抗:1.100000 交轴电抗:0.900000 直轴暂态电抗:0.300000 负序电抗:0.440000 Tj:1.540000 ----------------------------------------------------------------------------------- 变压器参数：变压器1容量:300.000000 高压侧额定电压:242.000000 低压侧额定电压电压:10.500000 电压百分数:12.000000变压器2容量:280.000000 高压侧额定电压:220.000000 低压侧额定电压电压:121.000000 电压百分数:14.000000----------------------------------------------------------------------------------- 输电线路参数：线路长度:230.000000 单位正序电抗:0.400000 单位负序电抗:0.400000 单位零序电抗=1.200000----------------------------------------------------------------------------------- 无穷大系统参数：故障前吸收有功:250.000000 功率因数:0.980000 恒定电压:115.000000----------------------------------------------------------------------------------- 故障地点：距线路始端0.650000, 故障类型：两相短路，故障点接地电抗为：j0.000000*************** 原始数据标么值归算************** 容量基准值:200.000000 电压基准值:115.000000 220KV电压等级电压基准值:209.090912----------------------------------------------------------------------------------- 发电机参数：额定有功功率:1.200000 额定电压:1.157391 功率因数:0.800000 直轴电抗:0.982340 交轴电抗:0.803733直轴暂态电抗:0.267911 负序电抗:0.392936 Tj:2.310000----------------------------------------------------------------------------------- 变压器参数：变压器T1容量:1.500000 高压侧额定电压:242.000000 低压侧额定电压电压:10.500000 电抗:0.107164 变压器T2容量:1.400000 高压侧额定电压:220.000000 低压侧额定电压电压:121.000000 电抗:0.110707----------------------------------------------------------------------------------- 输电线路参数：线路长度:230.000000 单位正序电抗:0.420870 单位负序电抗:0.420870 单位零序电抗=1.262609-----------------------------------------------------------------------------------无穷大系统参数：故障前吸收有功:1.250000 功率因数:0.980000恒定电压:1.000000-----------------------------------------------------------------------------------故障点接地电抗：0.000000****************** 正常运行时参数***************直轴暂态电抗和=0.696217交轴电抗和=1.232039虚构电动势E'=1.463569正常运行时候的功角=36.485753度。
电力系统暂态分析课内实验报告(稳定性)
实验报告
( 2006-- 2007年度第二学期)
名称：电力系统暂态分析课内实验题目：电力系统稳定性演示实验院系：电气与电子工程学院
班级：
学号：
学生姓名：
指导教师：
实验时数：2学时
成绩：
日期：2007年4月40日
电力系统暂态分析课内实验—电力系统稳定性演示
1 一、实验的目的与要求
通过该环节的学习，可以帮助学生进一步认识电力系统稳定的各种物理现象并增强对电力系统相关课程理论的理解。
二、实验正文
（１）大致了解WDT-III 型电力系统综合自动化实验台的功能
（2）了解如何绘制功角特性曲线
（3）观察电力系统失稳的现象
（4）观察不同类型的短路故障对电力系统稳定的影响
三、实验总结或结论
了解WDT-III 型电力系统综合自动化实验台
掌握电力系统机电暂态分析的原理
四、思考题
（1）WDT-III 型电力系统综合自动化实验台的功能有哪些？
（2）什么是功角特性曲线？其稳定运行点在哪儿？
（3）提高电力系统稳定性的措施有哪些？哪些措施通过本次实验仿真？
（４）实验中失去稳定的现象有哪些？
（５）不同类型的短路故障对电力系统稳定的影响有什么不同？
五、参考文献
1． 《电力系统暂态分析》，李光琦，中国电力出版社，1995年，第二版；
2． 《电力系统自动化》，李先彬，中国电力出版社，1986年，第二版；。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
课程设计报告( 2011—2012年度第二学期)名称：电力系统暂态上机院系：电气与电子工程学院班级：电气0914 班学号：1091420228学生姓名：雷刚指导教师：麻秀范老师设计周数：两周成绩：日期：2012年7月3日一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。
二、设计正文（详细内容见附录,用A4纸，页数不限）1．对称短路计算过程流程图和计算结果2．不对称短路计算过程流程图和计算结果3．静态稳定性计算过程流程图和计算结果4．暂态稳定性计算过程流程图和计算结果5．思考题三、课程设计总结或结论在这次电力系统暂态上机课程中，我们主要讨论了各种类型短路故障下（三相短路、单相接地短路、两相相间短路、两相接地短路、单相经电阻接地）系统网络中的电压、电流分布的计算以及电力系统的静态稳定、暂态稳定问题的分析。
通过本次课程设计，我对之前学习的知识有了一个巩固和提高，同时对电力系统故障分析有了更深刻的理解。
电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。
在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。
在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。
同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。
这与手算系统短路电流时使用的网络化简方法大大不同。
同时也明白了计算机编程方法对于电力系统稳态暂态计算的简便性、快速性和重要性。
最后，经过此次课程设计，对于我之前在数学建模中的学习到的Matlab知识起到了升华的作用，可以说这也算是一次比赛把，对较复杂的选择语句、循环语句的使用等有了更深的了解，为今后学习、工作中埋下一个伏笔，相信会受益匪浅。
四、参考文献1.《电力系统暂态分析》，李光琦，中国电力出版社，2006年，第三版；2.《电力系统分析》（上、下），何仰赞，华中科技大学出版社，1998年，第二版。
3.《电力系统故障的计算机辅助分析》重庆大学出版社米麟书等4.《电力系统潮流计算》天津大学出版社宋文南等5.《电力系统故障分析》清华大学出版社周荣光6.《短路电流实用计算方法》电力工业出版社西安交通大学等7.《精通Matlab6.5》北京航空航天大学出版社，张志涌附录（设计流程图、计算结果、思考题答案）1.对称短路计算过程流程图和计算结果流程图如下：计算结果如下：(1)导纳矩阵:Y1 =Y2=0 -26.6667i 0 +10.0000i 0 +10.0000i0 +10.0000i 0 -33.3333i 0 +10.0000i0 +10.0000i 0 +10.0000i 0 -20.0000iY0=0 -29.0476i 0 + 5.0000i 0 + 5.0000i0 + 5.0000i 0 -48.0228i 0 + 5.0000i0 + 5.0000i 0 + 5.0000i 0 -10.0000iYY1=YY2=0 -40.0000i 0 +10.0000i 0 +10.0000i 0 +20.0000i 00 +10.0000i 0 -60.0000i 0 +10.0000i 0 0 +40.0000i0 +10.0000i 0 +10.0000i 0 -20.0000i 0 00 +20.0000i 0 0 0 -30.0000i 00 0 +40.0000i 0 0 0 -60.0000i(2) 节点3发生三相短路故障●故障点三序电流、三相电流、三相电流有效值：I1= 0 - 9.8592i I2=0 I3=0 ;Ia= 0 - 9.8592i Ib=-8.5383 + 4.9296i Ic=8.5383 + 4.9296iIa_effective = Ib_effective = Ic_effective = 9.8592●各节点三序电压、三相电压：2．不对称短路计算过程流程图和计算结果不对称短路计算过程流程图如三相短路流程图。
(1) 节点3发生A相短路接地故障●故障点三序电流、三相电流：I1= 0 - 3.1239i I2= 0 - 3.1239i I3= 0 - 3.1239i ;Ia= 0 - 9.3718i Ib= 0 - 0.0000i Ic= 0 - 0.0000i ；(2) 节点3发生A相经10Ω电阻接地故障●故障点三序电流、三相电流：I1= 0.0333 - 0.0004i I2= 0.0333 - 0.0004i I3= 0.0333 - 0.0004i ; Ia= 0.1000 - 0.0011i Ib= 0 Ic= 0 ；(3) 节点3发生b、c两相短路故障●故障点三序电流、三相电流：I1= 0 - 4.9697i I2= 0 + 4.9697i I3= 0 ;Ia= 0 Ib= -8.6078 - 0.0000i Ic= 8.6078 + 0.0000i；●各支路三序电流、三相电流：(4) 节点3发生b、c两相短路接地故障●故障点三序电流、三相电流：I1= 0 - 6.4473i I2= 0 + 3.4921i I3= 0 + 2.9552i ; Ia= 0 Ib= -8.6078 + 4.4327i Ic= 8.6078 + 4.4327i ；3．静态稳定性计算过程流程图和计算结果静态稳定性计算过程流程图如下：计算结果：Ke=0.100000,delta=91,P=1.345117,Eq=1.999966 Ke=0.200000,delta=92,P=1.363000,Eq=2.028274 Ke=0.300000,delta=93,P=1.380693,Eq=2.056939 Ke=0.400000,delta=95,P=1.398409,Eq=2.089197 Ke=0.500000,delta=96,P=1.415747,Eq=2.119390 Ke=0.600000,delta=97,P=1.432842,Eq=2.149961 Ke=0.700000,delta=98,P=1.449679,Eq=2.180919 Ke=0.800000,delta=99,P=1.466244,Eq=2.212270 Ke=0.900000,delta=101,P=1.482360,Eq=2.251089 Ke=1.000000,delta=102,P=1.498265,Eq=2.283967 Ke=1.100000,delta=103,P=1.513819,Eq=2.317233 Ke=1.200000,delta=104,P=1.529000,Eq=2.350886 Ke=1.300000,delta=105,P=1.543786,Eq=2.384924 Ke=1.400000,delta=106,P=1.558152,Eq=2.419341 Ke=1.500000,delta=107,P=1.572075,Eq=2.454133 Ke=1.600000,delta=108,P=1.585527,Eq=2.489291 Ke=1.700000,delta=109,P=1.598483,Eq=2.524806 Ke=1.800000,delta=108,P=1.613666,Eq=2.534328 Ke=1.900000,delta=105,P=1.626492,Eq=2.515421 Ke=2.000000,delta=103,P=1.635259,Eq=2.507344 Ke=2.100000,delta=100,P=1.637025,Eq=2.483645 Ke=2.200000,delta=98,P=1.637979,Eq=2.471608 Ke=2.300000,delta=96,P=1.635815,Eq=2.457976 Ke=2.400000,delta=94,P=1.630657,Eq=2.442916 Ke=2.500000,delta=92,P=1.622630,Eq=2.426585 Ke=2.600000,delta=90,P=1.611867,Eq=2.409132 Ke=2.700000,delta=88,P=1.598501,Eq=2.390698 Ke=2.800000,delta=86,P=1.582670,Eq=2.371416 Ke=2.900000,delta=84,P=1.564509,Eq=2.351412 Ke=3.000000,delta=82,P=1.544155,Eq=2.330805 Ke=3.100000,delta=81,P=1.535174,Eq=2.323407 Ke=3.200000,delta=79,P=1.511558,Eq=2.301767 Ke=3.300000,delta=77,P=1.486103,Eq=2.279814 Ke=3.400000,delta=76,P=1.474144,Eq=2.271027 Ke=3.500000,delta=74,P=1.445989,Eq=2.248507 Ke=3.600000,delta=73,P=1.432473,Eq=2.239073 Ke=3.700000,delta=71,P=1.401945,Eq=2.216241 Ke=3.800000,delta=70,P=1.387071,Eq=2.206335 Ke=3.900000,delta=69,P=1.371732,Eq=2.196226 Ke=4.000000,delta=67,P=1.338431,Eq=2.173194 Ke=4.100000,delta=66,P=1.321987,Eq=2.162830 Ke=4.200000,delta=65,P=1.305171,Eq=2.152344 Ke=4.300000,delta=63,P=1.269626,Eq=2.129492 Ke=4.400000,delta=62,P=1.251928,Eq=2.118929 Ke=4.500000,delta=61,P=1.233937,Eq=2.108312 Ke=4.600000,delta=60,P=1.215672,Eq=2.097660 Ke=4.700000,delta=59,P=1.197153,Eq=2.086987 Ke=4.800000,delta=58,P=1.178398,Eq=2.076310 Ke=4.900000,delta=56,P=1.139925,Eq=2.054329 Ke=5.000000,delta=55,P=1.120623,Eq=2.043828 Ke=5.100000,delta=54,P=1.101144,Eq=2.033369 Ke=5.200000,delta=53,P=1.081504,Eq=2.022963 Ke=5.300000,delta=52,P=1.061715,Eq=2.012622 Ke=5.400000,delta=51,P=1.041792,Eq=2.002355 Ke=5.500000,delta=50,P=1.021746,Eq=1.992173 Ke=5.600000,delta=49,P=1.001589,Eq=1.982085 Ke=5.600000,delta=49,P=1.001589,Eq=1.972100最终选择放大倍数Ke=2.24．暂态稳定性计算过程流程图和计算结果暂态稳定性计算过程流程图如下：tDe l t a计算结果如下：0.15s 时切除故障的摇摆曲线 0.25s 时切除故障的摇摆曲线根据摇摆曲线判断：0.15s 时切除故障系统暂态稳定，0.25s 时切除故障系统的功角无限增大，系统失去稳定性。