电力系统潮流计算

1、潮流计算实验轮机工程学院 2014.09姚玉斌实验目的和内容目的：初步掌握潮流计算程序设计内容：第6周 潮流计算：程序各模块讲解第7周 潮流计算：数据输入模块第8周 潮流计算：主程序和数据文本文件编辑1 建立目录在D盘上建立一个目录 (以姓名的拼音命名) 在C盘已建目录：电力系统自动化实验/PowerFlow PowerFlow.cpp 源文件 PowerFlow.h 头文件PowerFlowLib.lib 库文件bus4图.doc 4节点算例接线图bus5图.doc 5节点算例接线图PowerFlow.dsw 工程文件双击，打开该工程2 使用Visaul C+6.0开始/程序/Microso

2、ft Viual Studio 6.0/ Microsoft Viual C+6.0 建立单个文件建立C+源文件选择目录输入文件名文件/新建2 使用Visaul C+6.0编写、编译、链接、执行程序在本窗口输入文件内容编译、链接、执行3 数据结构1）数据结构struct CLine / 支路数据 int i; / 首端母线号(k侧为负号)int j; / 末端母线号(k侧为负号)double r; / 电阻(在1侧)double x; / 电抗(在1侧)double yk; / B/2(或变比);3 数据结构1）数据结构struct CBus / 节点数据int i; / 母线号int Typ

3、e; / 节点类型(0:PQ，-1:PV，-2:平衡)double pg; / 有功出力double qg; / 无功出力double pl; / 有功负荷(负值)double ql; / 无功负荷(负值)double v0; / 电压给定值(PQ节点给定值为1.0);3 数据结构1）数据结构struct CCp / 电容器(或电抗器) int i; / 所在母线号double b; / 电纳值(电容器为正值，电抗器为负值);struct CYij / 导纳 double g; / 实部double b; / 虚部;3 数据结构1）数据结构struct CVp / 电压 double v; /

4、 幅值double va; / 相角;3 数据结构2）常量const int NS = 1000; / 节点数据数组大小const int ZS = 2000; / 支路数据数组大小const int NS2 = NS + NS; / 节点数据数组大小 * 2const int VA = -2; /节点类型：平衡节点const int PV = -1; /节点类型：PV节点const int PQ = 0; /节点类型：PQ节点3 数据结构3）数组CLine LineZS;CBus BusNS;CCp CpNS;CVp VpNS;CYij YijNSNS;double JacbNS2NS2;/

5、 雅可比矩阵double wNS2;/ 节点功率数组double bNS2;/ 功率不平衡量数组、解向量3 数据结构4）变量int n;/ 节点数int n2;/ int zls; / 支路数int ncp;/ 电容器个数int dsd; / 最大不平衡量的节点号int niter;/ 迭代次数int maxiter;/ 最大迭代次数double dsm; / 最大不平衡量double eps;/ 收敛精度FILE *fp1, *fp2;3 数据结构5）库函数void dataio(void);/ 从文件中读取数据void vini();/ 平启动电压初值void Pr_nod(void);/

6、 输出节点计算结果void branch(void); / 生成支路计算结果void y2(void); / 形成导纳矩阵void solution(void); / 牛顿法解方程void newtonp(void); / 极坐标牛顿法void jcbp(void);/ 雅可比矩阵void bbhl_p(void);/ 不平衡量及平衡节点功率、PV节点无功4 源文件1）包含文件4 源文件2）主函数4 源文件3）数据输入函数4 源文件3）数据输入函数4 源文件3）数据输入函数4 源文件3）数据输入函数4 源文件3）数据输入函数4 源文件3）数据输入函数4 源文件3）数据输入函数5 数据文件1)

7、例子：IEEE30节点系统的数据文件(局部) 30 41 2 0.00001 30 1 2 0.01920 0.05750 0.02640 1 3 0.04520 0.18520 0.02040 1 -2 0.00000 0.00000 0.00000 0.00000 1.05000 2 -1 0.57560 2.50000 -0.21700 -0.12700 1.03380 3 0 0.00000 0.00000 -0.02400 -0.01200 1.00000 10 0.19000 24 0.040005 数据文件2) 5节点系统图5 数据文件3) 4节点系统图6 极坐标牛顿法迭代主程序框图电压初始化计算残差P、Q修正量、V，修正、V计算雅可比矩阵P 、|Q形成导纳矩阵开始结束原始