电力系统短路故障的计算机算法程序设计汇总.docx
《电力系统短路故障的计算机算法程序设计汇总.docx》由会员分享，可在线阅读，更多相关《电力系统短路故障的计算机算法程序设计汇总.docx（39页珍藏版）》请在冰点文库上搜索。
电力系统短路故障的计算机算法程序设计汇总
电力系统短路故障的计算机算法程序设计
学生姓名
学号
成绩
设计题目
电力系统短路故障的计算机算法程序设计
设
计
内
容
电力系统故障的计算程序设计及编制和调试。
设
计
要
求
1．在对称短路计算、简单不对称短路计算中任选一种计算作为计算作业。
2．计算机语言自选。
3．设计、编制、调试出相关的通用计算程序。
4．输入／输出数据一律以文件格式形成。
（1）对称短路计算的输入／输出数据（供参考）
输入数据可考虑为以下两个部分：
数据1（可用变量表示）：
节点数、支路数、故障节点
数据2（可用二维数组表示）：
支路参数即
支路首端号、支路末端号、支路电抗值
输出数据（可用二维数组表示）为：
节点导纳矩阵或阻抗矩阵、故障节点、故障点电流、支路电流、各节点电压。
（2）简单不对称短路计算的输入／输出数据（供参考）
输入数据可考虑为以下两个部分：
数据1（可用变量表示）：
节点数、支路数、故障节点
数据2（可用二维数组表示）：
支路参数即
支路首端号、支路末端号、支路正序电抗值、负序电抗值、零序电抗值
输出数据（可用二维数组表示）为：
故障节点、故障点电流、各支路电流、各节点电压的序分量和相分量。
时
间
安
排
2周
参
考
资
料
[1]何仰赞，温增银．电力系统分析（上册）（第三版）．武汉：
华中科技大学出版社，2002．
[2]何仰赞，温增银．电力系统分析题解（上、下）（第三版）．武汉：
华中科技大学出版社，2002．
[3]郑平安，曾大亮编．程序设计基础（C语言）第二版．北京：
清华大学出版社，2006．
目录
1课程设计说明1
2编程语言的选择及理由2
3程序设计3
3.1主体流程图3
3.2详细流程图3
3.2.1创建系统（create）流程图3
3.2.2加载系统函数（load）流程图4
3.2.3计算子函数（calculate）流程图4
3.2.4改变短路节点（change）流程图5
3.3数据及变量说明6
3.4程序代码及注释7
3.5测试算例7
4设计体会17
参考文献18
附录19
1课程设计说明
本文根据电力系统三相对称短路的特点，建立了合理的三相短路的数学模型，在此基础上，形成电力系统短路电流实用汁算方法，编制了对任意一个电力系统在任意点发生短路故障时三相短路电流及其分布的通用计算程序，该办法适用于各种复杂结构的电力系统，从另一个侧面展示了计算机应用于电力系统的广阔前景。
根据所给的电力系统，编制短路电流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。
通过设计电力系统计算程序使自己对电力系统分析有进一步理解，同时加强计算机实际应用能力的训练。
电力系统的短路故障是严重的，而又是发生几率最多的故障。
一般说来，最严重的短路是三相短路。
当发生短路时，其短路电流可达数万安甚至十几万安，它们所产生的热效应和电动力效应将使电气设备遭受严重破环。
为此，当发生短路时，继电保护装置必须迅速切除故障线路，以避免故障部分继续遭受危害，并使非故障部分从不正常运行情况下解脱出来，这要求电气设备必须有足够的机械强度和热稳定度，开关电气设备必须具备足够的开断能力，即必须经得起最大短路电流的侵扰而不致损坏。
因此，电力系统短路电流计算是电力系统运行分析和设计计算的重要环节，许多电力设计单位和个人倾注极大精力从事这一工作的研究。
由于电力系统结构复杂，随着生产发展，技术进步，电力系统日趋扩大和复杂化，短路电流计算的工作量也随之增大，采用计算机辅助计算势在必行。
2编程语言的选择及理由
一种语言之所以能存在和发展，并具有生命力，总是有其不同于（或优于）其他语言的特点。
C语言的主要特点如下。
（1）语言简洁、紧凑，使用方便、灵活。
C语言一共只有32个关键字，9种控制语句，程序书写形式自由，主要用小写字母表示，压缩了一切不必要的成分。
（2）运算符丰富。
C的运算符包含的范围很广泛，共有34种运算符。
C把括号、赋值、强制类型转换等都作为运算符处理，从而使C的运算类型极其丰富，表达式类型多样化。
灵活使用各种运算符可以实现在其他高级语言中难以实现的运算。
（3）数据结构丰富，具有现代化语言的各种数据结构。
C的数据类型有整型、实型、字符型、数组类型、指针类型、结构体类型、共用体类型等。
能用来实现各种复杂的数据结构（如链表、树、栈等）的运算。
尤其是指针类型数据，使用起来比PASCAL更为灵活、多样。
（4）具有结构化的控制语句（如if…else语句、while语句、do…while语句、switch语句、for语句）。
用函数作为程序的模块单位，便于实现程序的模块化。
C是良好的结构化语言，符合现代编程风格的要求。
（5）语法限制不太严格，程序设计自由度大。
例如对数组下标越界不做检查，由程序编写者自己保证程序的正确。
对变量的类型使用比较灵活，例如整型数据与字符型数据可以通用。
一般的高级语言语法检查比较严，能检查出几乎所有的语法错误。
而C语言允许程序编写者有较大的自由度，因此，放宽了语法检查。
程序员应当仔细检查程序，保证其正确，而不要过分依赖C编译程序去查错。
“限制”与“灵活”是一对矛盾。
限制严格，就失去灵活性;而强调灵活，就必然放松限制。
一个不熟练的编程人员，编一个正确的C程序可能会比编一个其他高级语言程序难一些。
也就是说，对用C语言的人，要求对程序设计更熟练一些。
（6）C语言能进行位（bit）操作，能实现汇编语言的大部分功能，可以直接对硬件进行操作。
因此C既具有高级语言的功能，又具有低级语言的许多功能，可用来写系统软件。
C语言的这种双重性，使它既是成功的系统描述语言，又是通用的程序设计语言。
3程序设计
3.1主体流程图
图3-1主体流程图
3.2详细流程图
3.2.1创建系统（create）流程图
图3-2create函数流程图
3.2.2加载系统函数（load）流程图
图3-3节点导纳矩阵形成流程图
3.2.3计算子函数（calculate）流程图
图3-4计算子函数流程图
3.2.4改变短路节点（change）流程图
图3-5change函数流程图
3.3数据及变量说明
本程序中将电路的支路分为三类：
电抗支路、电流支路、发电机直流，数据输入时分开输入格式如下：
电抗支路：
支路数、节点首端好、节点末端好，支路电抗值。
发电机支路：
发电机台数、发电机接入点，发电机电势、发电机之路电抗。
电容支路：
支路数、电容节点号，电容支路电抗值。
程序中的变量定义如下：
FILE*fp;//文件指针
floatIf;//短路电流
floatV0[10];//节点电压初值
floatIn[10];//注入电流
floatV[10];//节点电压
floatI[10];//支路电流
intcount[2][10];//数组用于储存支路号
floatz[10][10];//储存支路电抗值
floatyjs[10][10];//存放计算导纳矩阵中的中间值
floatye[10];//储存发电机支路导纳
floatY[10][10];//节点导纳矩阵
floatZ[10][10];//阻抗矩阵
floatdk;//电抗值
inte1[10];
floatyc[10];//节点对地电抗
floatycjs[10];
inte[10];//发电机接入点
floatE[10];
floatIg[10];
floatEmid[10];
floatze[10];//储存发电机支路导纳
intm,n,a,f,i=0,j=0,sd,md,butn,times,timesc;/*n是节点数，a是支路数，f是短路点，i为行标，j为列标*/
1、文件读入程序如下：
printf（"请输入系统的节点数，支路数，故障节点\n"）;
scanf（"%d%d%d",&n,&a,&f）;
fp=fopen（name,"w+"）;
fprintf（fp,"系统参数如下所示:
\n\n"）;
fprintf（fp,"1.系统的节点数为%d\n2.系统的支路数为%d\n3.系统的短路点为%d\n\n",n,a,f）;
fprintf（fp,"支路电抗参数如下\n"）;
2、计算主程序如下：
（1）、If=V0[f-1]/Z[f-1][f-1];（短路电流计算）
for（i=0;iV[i]=V0[i]-Z[i][f-1]*If;
（2）、for（i=0;i{
for（j=0;jY[i][i]=Y[i][i]-yjs[i][j];
Y[i][i]=Y[i][i]-ye[i]-ycjs[i];
}
for（i=0;ifor（j=i+1;j{
Y[i][j]=yjs[i][j];
Y[j][i]=yjs[i][j];
}
3.4程序代码及注释
具体源程序代码及其注释见附录。
3.5测试算例
【例6-3】在如图3-5所示的电力系统中分别在节点1和节点5接入发电机支路，其标幺值参数为：
。
在节点3发生三相短路，计算短路电流及网络中的电流分布。
线路的电阻和电容略去不计，变压器的标幺变比等于1。
各元件参数的标幺值如下：
图3-6电力系统等值网络图
图3-7三相短路时的等值网络图
由以上网络图可得计算过程如下：
（1）、形成输入文件：
系统参数如下所示:
1.系统的节点数为5
2.系统的支路数为5
3.系统的短路点为3
支路电抗参数如下
z12=0.1050z23=0.0650z24=0.0800z43=0.0500z54=0.1840
发电机接入点及其电抗如下：
共2台发电机进入系统
z1=0.1500z5=0.2200
e1=1.0000e5=1.0000
共0个节点有接地电容
（2）、计算后得到的输出文件：
计算结果如下：
1、节点导纳矩阵:
-16.19059.5238
9.5238-37.408415.384612.5000
15.3846-35.384620.0000
12.500020.0000-37.93485.4348
5.4348-9.9802
2、节点阻抗矩阵:
0.11810.09580.09020.08580.0467
0.09580.16290.15330.14590.0794
0.09020.15330.18600.16110.0877
0.08580.14590.16110.17290.0941
0.04670.07940.08770.09410.1515
3、短路节点是:
3
4、短路电流:
If=5.3767
5、各节点电压:
V1=0.5152V2=0.1758V3=0V4=0.1336V5=0.5282
6、发电机支路电流如下:
Ig[1][0]=-3.2321
Ig[5][0]=-2.1445
7、各支路电流:
I12=-3.2321I23=-2.7046I24=-0.5275I43=-2.6720I54=-2.1445
图3-8对称短路计算结果
图3-8改变短路点后计算结果
【习题6-10】如图3-9所示网络，略去负荷，试用节点阻抗矩阵法求节点5发生三相短路时，短路点的短路电流及线路L-2、L-3的电流。
已知各元件参数如下：
发电机G-1:
SN=120MVA,x”d=0.23;G-2:
SN=60MVA,x”d=0.14
变压器T-1:
SN=120MVA,Vs=10.5%;T-2:
SN=60MVA,Vs=10.5%
线路参数x1=0.4Ω/km,b1=2.8×10-6S/km
线路长度L-1:
120km,L-2:
80km,L-3:
70km
取SB=120MVA,VB=Vav
图3-9电力系统网络图
图3-10三相短路时的等值网络图
（1）、形成输入文件：
系统参数如下所示:
1.系统的节点数为：
5
2.系统的支路数为：
5
3.系统的短路点为：
5
支路电抗参数如下：
z35=0.29036z45=0.25406z34=0.43554z13=0.1050z24=0.2100
发电机接入点及其电抗如下：
共2台发电机进入系统：
z1=0.2300z2=0.2800
e1=1.0500e2=1.0500
共3个节点有接地电容：
Y3=-32.372299Y4=-34.077000Y5=-43.185200
（2）、得到的输出文件：
计算结果如下：
1、节点导纳矩阵:
-13.87169.5238
-8.33334.7619
9.5238-15.23292.29603.4440
4.76192.2960-10.96463.9361
3.44403.9361-7.3569
2、节点阻抗矩阵:
0.18190.06160.16000.10790.1326
0.06160.20780.08980.15360.1242
0.16000.08980.23300.15710.1932
0.10790.15360.15710.26890.2174
0.13260.12420.19320.21740.3427
3、短路节点是:
5
4、短路电流：
If=3.1263
5、各节点电压:
V1=0.6471V2=0.6723V3=0.4632V4=0.3890V5=0
6、发电机支路电流如下:
Ig[1][0]=-1.7516
Ig[2][0]=-1.3491
7、各支路电流:
I35=-1.5954I45=-1.5310I34=-0.1705I13=-1.7516I24=-1.3491
图3-11对称短路计算结果
【习题6-11】电力系统等值电路如图3-13所示，支路阻抗的标幺值已注明图中。
（1）形成节点导纳矩阵（或节点阻抗矩阵），并用以计算节点3的三相短路电流。
（2）另选一种方法计算短路电流，并用以验证
（1）的计算结果。
图3-12电力系统等值电路图
图3-13经Norton变换后的等值电路图
图3-13最终计算等值电路图
由以上网络可得输入文件如下：
系统参数如下所示:
1.系统的节点数为：
3
2.系统的支路数为：
2
3.系统的短路点为：
3
支路电抗参数如下：
z12=0.2000z23=0.1000
发电机接入点及其电抗如下：
共2台发电机进入系统：
z1=0.2000z3=0.130435
e1=1.0500e3=0.913043
共1个节点有接地电容：
Y2=-2.000000
（2）、得到的输出文件：
计算结果如下：
1、节点导纳矩阵:
-10.00005.0000
5.0000-14.500010.0000
10.0000-17.6667
2、节点阻抗矩阵:
0.13940.07890.0446
0.07890.15770.0893
0.04460.08930.1071
3、短路节点是:
3
4、短路电流：
If=9.1875
5、各节点电压:
V1=0.6344V2=0.2188V3=0
6、发电机支路电流如下:
Ig[1][0]=-2.0781
Ig[3][0]=-7.0000
7、各支路电流:
I12=-2.0781I23=-2.1875
图3-15对称短路计算结果
通过计算结果与算例答案的对比，计算结果与例题标准答案一样，说明了程序计算方法和思路的准确性。
说明本次课程设计基本达到了设计的要求。
4设计体会
从拿到课题到设计完成，从理论到实践，忙于查找资料，构思方案。
但是却学到了很多的东西，在巩固以前所学过的知识的同时，更是学到了很多在书本上所没有学到过的知识。
课程设计是培养我们运用所学知识，发现、提出、分析和解决实际问题、锻炼实践能力的重要环节，是对我们所学知识的具体训练和考察。
通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正服务于社会，从而提高自己的实际动手能力和独立思考的能力。
在设计的过程中遇到问题，可以说得是困难重重，难免会遇到过各种各样的问题，同时在设计的过程中发现了自己的不足之处，对以前所学过的知识理解得不够深刻，掌握得不够牢固。
通过这次设计，我懂得了学习和团队合作的重要性，了解到理论知识与实践相结合的重要意义，学会了坚持、耐心和努力，这将为自己今后的学习和工作做出了最好的榜样。
参考文献
[1]何仰赞，温增银．电力系统分析（上册）（第三版）．武汉：
华中科技大学出版社，2002．
[2]何仰赞，温增银．电力系统分析题解（上、下）（第三版）．武汉：
华中科技大学出版社，2002．
[3]郑平安，曾大亮编．程序设计基础（C语言）第二版．北京：
清华大学出版社，2006．
附录
主程序及注释（这里只给出主要计算程序）：
#include
#include
#include
声明外部函数
externvoidcreate（）;
externvoidload（）;
externvoidquit（）;
externvoidchange（charname[20]）;
externvoidcalculate（charname[20]）
{
FILE*fp;//文件指针
floatIf;//短路电流
floatV0[10];//节点电压初值
floatIn[10];//注入电流
floatV[10];//节点电压
floatI[10];//支路电流
intcount[2][10];//数组用于储存支路号
floatz[10][10];//储存支路电抗值
floatyjs[10][10];//存放计算导纳矩阵中的中间值
floatye[10];//储存发电机支路导纳
floatY[10][10];//节点导纳矩阵
floatZ[10][10];//阻抗矩阵
floatdk;//电抗值
inte1[10];
floatyc[10];//节点对地电抗
floatycjs[10];
inte[10];//发电机接入点
floatE[10];
floatIg[10];
floatEmid[10];
floatze[10];//储存发电机支路导纳
intm,n,a,f,g,i=0,j=0,sd,md,butn,times,timesc;/*n是节点数，a是支路数，f是短路点，i为行标，j为列标
for（i=0;i<10;i++）
{yc[i]=0;
ycjs[i]=0;
E[i]=0;
Emid[i]=0;}
for（i=0;i<10;i++）//初始赋
{
V0[i]=0;
for（j=0;j<10;j++）
{
z[i][j]=0;
Y[i][j]=0;
}
}
/******************************读取系统参数**************************/
fp=fopen（name,"r"）;
fscanf（fp,"系统参数如下所示:
\n\n"）;
fscanf（fp,"1.系统的节点数为%d\n2.系统的支路数为%d\n3.系统的短路点为%d\n\n",&n,&a,&f）;
/**************读电抗*************/
fscanf（fp,"支路电抗参数如下\n"）;
for（m=0;m{
fscanf（fp,"z%d%d=%f\t",&sd,&md,&dk）;
i=sd-1;
j=md-1;
count[0][m]=sd;
count[1][m]=md;
z[i][j]=dk;
z[j][i]=dk;
}
printf（"系统支路电抗如下：
\n"）;
for（m=0;m{
i=count[0][m]-1;
j=count[1][m]-1;
printf（"z%d%d=%.4f\n",count[0][m],count[1][m],z[i][j]）;
}
fscanf（fp,"\n\n发电机接入点及其电抗如下：
"）;
fscanf（fp,"共%d台发电机进入系统\n",×）;
for（m=0;mfscanf（fp,"z%d=%f\t",&e[m],&ze[m]）;
fscanf（fp,"\n"）;
for（m=0;mfscanf（fp,"e%d=%f\t",&e[m],&E[m]）;
printf（"输入的发电机支路电抗为\n"）;
for（m=0;mprintf（"z%d=%f\t",e[m],ze[m]）;
for（m=0;mfor（m=0;mprintf（"e%d=%f\t",e[m],E[m]）;
for（m=0;mprintf（"%d",e[m]）;
/****************读取节点对地电容**********************/
fscanf（fp,"\n共%d个节点有接地电容\n",×c）;
if（timesc>0）
{for（m=0;mfscanf（fp,"Y%d=%f\t",&e1[m],&yc[m]）;
printf（"\n系统有%d个节点电容\n",timesc）;
for（m=0;mprintf（"z%d0=%f",e1[m],yc[m]）;
}
fclose（fp）;
/*******************************计算节点阻抗矩阵**************************/
for（m=0;m<10;m++）
ye[m]=0;
for（m=0;mfor（j=0;j
升级会员 