潮流上机课程设计报告华电

华址电力*孑

《电力系统潮流上机》课程设计报告

院系

班级：

学号：

学生姓名: 指导教师: 设计周数

成绩：

日期：年月日

q 「十•-课程课程设计报告

、课程设计的目的与要求

培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识

、设计正文（详细内容见附录）

1.手算：要求应用牛顿-拉夫逊法或P-Q分解法手算求解，要求精度为0.001MW

节点1为平衡节点，电压U, 1.0 0，节点2为PQ节点，负荷功率S20.8 j0.6，节点3 是PV 节点，P3 04U3 1.1，两条支路分别为Z13 0.01 j0.04，Z12 0.05 j0.2，对地支路y30j 0.33。

» 十•-课程课程设计报告

2•计算机计算：编写潮流计算程序，要求如下：

2.1据给定的潮流计算任务书整理潮流计算的基础数据：节点的分类，线路模型，等值变压器 模型，电压等级的

归算，标幺值的计算；

2.2基础数据的计算机存储：节点数据，支路数据（包括变压器） ；

2.3用牛顿-拉夫逊法计算；

2.4根据所选潮流计算方法画流程图，划分出功能模块，有数据输入模块，导纳阵形成模块， 解线性方程组模

块，计算不平衡功率模块，形成雅可比矩阵模块，解修正方程模块，计算 线路潮流，网损，PV 节点无功功率和平衡节点功率，数据输出模块；

2.5据上述模块编制程序并上机调试程序，得出潮流计算结果； 2.6源程序及其程序中的符号说明集、程序流图

简单系统如下图所示，支路数据如下： 乙2 0.1 j0.41，乙3

j0.3， z 14

0.12 j0.5， z 24 0.08 j0.40

y io,2 y 20,1 j 0.01528， y 10,4

y 40,1

j 0.0192， y 20,4

y 40,2 j 0.01413

k 1.1

节点数据如下：

S 1 0.30 j0.18，S 2 0.55 j0.13， S 3

0.5，U 3

1.10，U 4

1.05 0o

Z 13

Z 13

y

40,1

y

20,4

1) 节点导纳阵

#in elude #in elude #in elude #i nclude "LF.h"

y

10,4

Z

14

k(k

丄 y 20,1

Z 12

Z

13

S 4

Z 24

S 2

y

40,2

//form node conductance matrix

intMakeY( intnB, intnL, Line* sL, double** YG, double** YB ) {

in ti,j,l;

double r,x,d1,g,b,t;

for(i=0;i< nB;i++)

for(j=0;j< nB;j++)

{ YG[i][j]=0.0;

YB[i][j]=O.O;}

for(i=0;i

{r=sL[i].R;

x=sL[i].X; g=r/(r*r+x*x); b=-x/(r*r+x*x); switch(sL[i].Type)

{case 1://L ine

break;

case 2://Tra nsformer

g*=1/sL[i].K;

b*=1/sL[i].K; break;} YG[sL[i].Numl][sL[i].Numl]+=g;

YG[sL[i].NumJ][sL[i].NumJ]+=g; YG[sL[i].Numl][sL[i].NumJ]-=g;

YG[sL[i].NumJ][sL[i].Numl]-=g; YB[sL[i].Numl][sL[i].Numl]+=b+sL[i].B; YB[sL[i].NumJ][sL[i].NumJ]+=b+sL[i].B;

YB[sL[i].Numl][sL[i].NumJ]-=b;

YB[sL[i].NumJ][sL[i].Numl]-=b; }

prin tf("实部：\n");

for(i=0;i< nB;i++)

{for(j=0;j< nB;j++)

prin tf("%lf\t",YG[i][j]);

prin tf("\n");}

printf("虚部：\n”);

for(i=0;i< nB;i++)

{for(j=0;j< nB;j++)

prin tf("%lf\t",YB[i][j]);

prin tf("\n");}

/* Check the Y matrix */

ofstreamfout("out.txt");

fout<< " ------------- Y Matrix ------------------- " <

for(i=0;i< nB;i++)

{for(j=0;j< nB;j++)

fout<< YG[i][j] << "+j" << YB[i][j] << "\t";

fout<

fout.close();

return 0;}

2) 计算功率不平衡量

#in elude

#in clude

#in clude