电力系统潮流的计算及其分析控制(下)

原创 2005年04月29日 09:33:00

第5章  PQ程序框图及程序说明

5.1   主程序流程图

初始化节点和支路各参数值

通过统计各节点静态出线数进行节点重新排序

PV节点加以标志并保存在一数组里

计算B/和B//,相应地形成因子表1和2

初始化电压幅值和相角

num=0

While(num<ITM)

num=num+1;KG[0]=1; KG[1]=1;

计算各节点有功和无功

KQO=0;

进行有功功率误差计算,并标出最大功率误差

fabs(DSM)<EPS?

false                                         true

解方程

 

KG[0]=0

相位调整

重新计算各节点有功和无功

KQO=1

进行有功功率误差计算,并标出最大功率误差

(fabs(DSM)<EPS)?

false                                          true

解方程

KG[1]=0

幅值调整

(KG[0]==0&&KG[1]==0)?

False                                          true

num>=ITM?

false               true     

对节点编号进行还原

计算节点和支路数据

 

输出提示信息

输出节点和支路数据

break

break;

5.2   子程序流程图

1.      chushihua()函数 

打开文件,

输入节点信息

关闭文件

打开文件,

输入支路信息

关闭文件

 

2.      jiedianyh()函数

根据支路信息统计各节点出线数

根据出线数设置将要排序顺序号(平衡节点排最后,)形成出线表

支路左、右节点号根据出线数表对应的进行转化

支路数据左右节点号调整(左小于右)

节点数据根据出线数表对应进行交换排序

数据输出

 

3.      daona1()函数

只用于不接地支路构成导纳矩阵,不考虑接地支路(包括变压器非标准变比)的影响,主要用来形成PQ分解法所要求的第一个因子表。

i=1步长1N

gii(i)=0;  bii(i)=0;  yds(i)=0;

k=1步长1ZLS

i=| zl[k].zli |;  j=| zl[k].zlj |;

r=| zl[k].zlr |;  x=| zl[k].zlx |;

zf=r*r+x*x;  

gij=r/zf;      bij= -x/zf;

y[k].yg= -gij;  y[k].yb= -bij;  y[k].cow=j;

gii[i]= gii[i]+gij; bii[i]= bii[i]+bij;

gii[j]= gii[j]+gij; bii[j]= bii[j]+bij;

yds[k]=yds[i]+1

ydz[i]=1

i=1步长1N

ydz[i+1]=ydz[i]+yds[i]

 

4.      daona2()函数

追加接地支路的程序框图

k=1步长1zls

i= zl[k].zli ;  j= zl[k].zlj ;

yk= zl[k].zlyk ;

i<0 V j<0

true                                        false

i<0?

true                 false

gij=0

bij=yk/2

i=abs(i)

i=abs(j)

sy(i)

sy(j)

gij=y[k].yg; bij= y[k].yb

gii[i]=gii[i]+(1-1/yk/yk)*gij

bii[i]= bii[i]+(1-1/yk/yk)*bij;

 

y[k].yg=gij/yk

y[k].yb=bij/yk

sy(i),sy(j)完成统计各节点自导计算

5.      yinzibiao()函数

(计算部分,不含前面判断及传送)

根据t值来确定对导纳阵付初值

对其他变量fd[I],di[I],u[I],m付初值

for(k=1;k<=m;k++)

对临时工作数组B[I] 付初值

x=2    im=1

while(im<=k-1)

jsq=1

jsq>fd[im]

 true                                     false

im=im+1

 

while(jsq<=fd[im])

u[x+1]!=k?

false                                      true

ai=u[x]/di[im]

jsq=jsq+1;    x=x+2;

while(jsq<=fd[im])

j=u[x+1]

b[j]=b[j]-ai*u[x]

x=x+2

jsq=jsq+1

 

im=im+1;

im>k-1?

true                                            false

di[k]=1/b[k];  ai=di[k];  jsq=0;

 

for(j=k+1;j<=m;j++)

b[j]!=0

true                false

u[x]=b[j]*ai;

 

u[x+1]=j

jsq=jsq+1

x=x+2

fd[k]=jsq

结果保存

5.2.6  其它函数其流程图不再列出;

5.3   各函数功能介绍

函数名

函数完成功能

PQjs()

函数主体,其它函数均为其子函数,完成计算

chushihua()

原始数据的输入

jiedianyh()

通过静态最少出线数优先完成节点优化排序

PVbz()

形成PV节点数组

daona1()

形成计算因子表1所需要的导纳矩阵

Daona2()

形成计算因子表1所需要的导纳矩阵

yinzibiao()

根据T值形成因子表1或2

InputV()

初始化V值

gljs()

完成节点有功和无功的计算

glwc()

计算节点功率最大误差

solution()

解方程

dyxz()

调整节点电压幅值和相角

jiedianhy()

把通过节点优化后的数据还原为以前的编号

nod()

计算节点结果数据

branch()

计算支路结果数据

shuchu()

输出节点、支路结果数据

 

5.4   程序变量说明

1.        N:系统节点的总数;WGS:系统中发电机节点的总数;

WLS:系统中负载节点的总数;EPS:迭代所要求的精度;

ITM:最大能迭代次数;

2.        struct WFL       

{

       double P; PQ节点:有功功率)

       double Q; PQ节点:无功功率,PV节点:无功功率的上限)

       int node;  (节点编号)

       double V;  (负号是PV节点的标志)

}WG[WGS],WL[WLS];

3.        struct ZL         

{

       int node1;  (一端的节点号)

       int node2;  (另一端的节点号)

       double R;  (电阻)

       double X;  (电抗)

       double K;   (输电线支路:容纳;变压器支路:非标准变比)

}pZL[ZLS]; 

4.        (*pPV)[2];      PV节点数组)  

5.        jdyh[N][2] ;     (节点优化结果数组)

6.        GII[NMAX],BII[NMAX];YDS[NMAX],YDZ[NMAX];(导纳对角元素数祖实部、虚部,各行非零非对角元个数数组,每行第一非零非对角远在Y中序号数组)

7.        DI[N],B[N];(因子表对角元素数组,计数展开用数组)

8.        Y[ZLSMAX][3]        (非对角元数组;

9.        *FD;         (因子表上三角阵各行非零元素个数)

10.    R,X,GIJ,BIJ,YK,ZF;   (电阻、抗、非对角原实部、虚部、变比容抗暂存单元)

11.    LD,PVD,GD,L,S,G;(负荷节点、PV节点、发电机节点、负荷计数、PV点计数、发电机,计数寄存单元)

12.    T (特征计数单元:  0对应B1;   1对应B2;

13.    KQO; (特征计数单元,kqo=0控制P迭代,kqo=1控制Q迭代)

14.    KG[2];         (收敛检验标识)

15.    DSD;      (迭代最大功率误差节点号)

16.    DSM;      (迭代最大功率误差)

17.    OI[NMAX];       (方程输入输出单元)

18.    W[N][2],V[N][2];  (功率、电压数组)

19.    VI,VJ,AI,BI,CI,WI;  (临时计数单元)

 

结论

 

致谢  

参考文献

1        何仰赞 温增银 旺馥瑛 周勤慧.电力系统分析(上、下册).华中理工大学出版社,1995.4

2        西安交通大学等六院校合编.电力系统计算.水利电力出版社 1995.5

3        刘浩.保留非线形的快速PQ分解状态估计程序设计.电力系统及自动化学报.1997.12

4        电力工业部电力规划设计总院编.电力系统设计手册.中国电力出版社.1998

5        王锡凡主编.电力系统优化规划.水利电力出版社.1990

6        A·R伯根 原著.刘康立译.电力系统分析.晓圆出版社.世界图书出版社.1994

7        H·E·BrownSolution of large networks by matrix methodsJohn Wiley&SonsInc1975

8        J·ArrillagaC·P·ArnoldB·J·HarkerComputer Modelling of Electrical Power SystemsJohn Wiley&SonsInc1983

9         

 

附录1(源码略)

 

 

电力系统潮流计算 迭代算法

#include #include #include #include #include #include void input();double num[3];       //结点个数,线路个数,...

期末总结-电力系统潮流上机计算

随着敲下明天要交的报告的最后一个句号,本学期的在校学习任务就算正式结束了。 已经快三周没有好好休息了吧,期末考完之后关于ACM的东西也没学几天,还错过了各种想参加练手的比赛。。。。本以为“电力系统潮...

电力系统69节点算例潮流计算

  • 2017年10月06日 16:24
  • 16KB
  • 下载

电力系统上机潮流计算

  • 2015年03月23日 13:54
  • 34KB
  • 下载

PSS E v33.40 1CD(大型电力系统仿真计算软件)

美国PTI公司产品: PSS E v33.40 1CD(大型电力系统仿真计算软件) PSS/E v32.03-ISO 1CD PSS/E v31.0-ISO 1CD PSS E v31.20 ...

基于MATLAB和VB的电力系统潮流计算

  • 2017年07月14日 14:53
  • 2.73MB
  • 下载

电力系统远程控制项目总结

很好,终于搞定了整个project,作为一个有始有终的人,咳咳,今天我来完成这个summary。 时间:2012年11月-2013年2月 人物:三个中国人+四个芬兰人+一个波兰人+巴基斯坦人...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:电力系统潮流的计算及其分析控制(下)
举报原因:
原因补充:

(最多只允许输入30个字)