電力系統(tǒng)潮流計(jì)算

1、潮流計(jì)算實(shí)驗(yàn)輪機(jī)工程學(xué)院 2014.09姚玉斌實(shí)驗(yàn)?zāi)康暮蛢?nèi)容目的：初步掌握潮流計(jì)算程序設(shè)計(jì)內(nèi)容：第6周 潮流計(jì)算：程序各模塊講解第7周 潮流計(jì)算：數(shù)據(jù)輸入模塊第8周 潮流計(jì)算：主程序和數(shù)據(jù)文本文件編輯1 建立目錄在D盤上建立一個(gè)目錄 (以姓名的拼音命名) 在C盤已建目錄：電力系統(tǒng)自動化實(shí)驗(yàn)/PowerFlow PowerFlow.cpp 源文件 PowerFlow.h 頭文件PowerFlowLib.lib 庫文件bus4圖.doc 4節(jié)點(diǎn)算例接線圖bus5圖.doc 5節(jié)點(diǎn)算例接線圖PowerFlow.dsw 工程文件雙擊，打開該工程2 使用Visaul C+6.0開始/程序/Microso

2、ft Viual Studio 6.0/ Microsoft Viual C+6.0 建立單個(gè)文件建立C+源文件選擇目錄輸入文件名文件/新建2 使用Visaul C+6.0編寫、編譯、鏈接、執(zhí)行程序在本窗口輸入文件內(nèi)容編譯、鏈接、執(zhí)行3 數(shù)據(jù)結(jié)構(gòu)1）數(shù)據(jù)結(jié)構(gòu)struct CLine / 支路數(shù)據(jù) int i; / 首端母線號(k側(cè)為負(fù)號)int j; / 末端母線號(k側(cè)為負(fù)號)double r; / 電阻(在1側(cè))double x; / 電抗(在1側(cè))double yk; / B/2(或變比);3 數(shù)據(jù)結(jié)構(gòu)1）數(shù)據(jù)結(jié)構(gòu)struct CBus / 節(jié)點(diǎn)數(shù)據(jù)int i; / 母線號int Typ

3、e; / 節(jié)點(diǎn)類型(0:PQ，-1:PV，-2:平衡)double pg; / 有功出力double qg; / 無功出力double pl; / 有功負(fù)荷(負(fù)值)double ql; / 無功負(fù)荷(負(fù)值)double v0; / 電壓給定值(PQ節(jié)點(diǎn)給定值為1.0);3 數(shù)據(jù)結(jié)構(gòu)1）數(shù)據(jù)結(jié)構(gòu)struct CCp / 電容器(或電抗器) int i; / 所在母線號double b; / 電納值(電容器為正值，電抗器為負(fù)值);struct CYij / 導(dǎo)納 double g; / 實(shí)部double b; / 虛部;3 數(shù)據(jù)結(jié)構(gòu)1）數(shù)據(jù)結(jié)構(gòu)struct CVp / 電壓 double v; /

4、 幅值double va; / 相角;3 數(shù)據(jù)結(jié)構(gòu)2）常量const int NS = 1000; / 節(jié)點(diǎn)數(shù)據(jù)數(shù)組大小const int ZS = 2000; / 支路數(shù)據(jù)數(shù)組大小const int NS2 = NS + NS; / 節(jié)點(diǎn)數(shù)據(jù)數(shù)組大小 * 2const int VA = -2; /節(jié)點(diǎn)類型：平衡節(jié)點(diǎn)const int PV = -1; /節(jié)點(diǎn)類型：PV節(jié)點(diǎn)const int PQ = 0; /節(jié)點(diǎn)類型：PQ節(jié)點(diǎn)3 數(shù)據(jù)結(jié)構(gòu)3）數(shù)組CLine LineZS;CBus BusNS;CCp CpNS;CVp VpNS;CYij YijNSNS;double JacbNS2NS2;/

5、 雅可比矩陣double wNS2;/ 節(jié)點(diǎn)功率數(shù)組double bNS2;/ 功率不平衡量數(shù)組、解向量3 數(shù)據(jù)結(jié)構(gòu)4）變量int n;/ 節(jié)點(diǎn)數(shù)int n2;/ int zls; / 支路數(shù)int ncp;/ 電容器個(gè)數(shù)int dsd; / 最大不平衡量的節(jié)點(diǎn)號int niter;/ 迭代次數(shù)int maxiter;/ 最大迭代次數(shù)double dsm; / 最大不平衡量double eps;/ 收斂精度FILE *fp1, *fp2;3 數(shù)據(jù)結(jié)構(gòu)5）庫函數(shù)void dataio(void);/ 從文件中讀取數(shù)據(jù)void vini();/ 平啟動電壓初值void Pr_nod(void);/

6、 輸出節(jié)點(diǎn)計(jì)算結(jié)果void branch(void); / 生成支路計(jì)算結(jié)果void y2(void); / 形成導(dǎo)納矩陣void solution(void); / 牛頓法解方程void newtonp(void); / 極坐標(biāo)牛頓法void jcbp(void);/ 雅可比矩陣void bbhl_p(void);/ 不平衡量及平衡節(jié)點(diǎn)功率、PV節(jié)點(diǎn)無功4 源文件1）包含文件4 源文件2）主函數(shù)4 源文件3）數(shù)據(jù)輸入函數(shù)4 源文件3）數(shù)據(jù)輸入函數(shù)4 源文件3）數(shù)據(jù)輸入函數(shù)4 源文件3）數(shù)據(jù)輸入函數(shù)4 源文件3）數(shù)據(jù)輸入函數(shù)4 源文件3）數(shù)據(jù)輸入函數(shù)4 源文件3）數(shù)據(jù)輸入函數(shù)5 數(shù)據(jù)文件1)

7、例子：IEEE30節(jié)點(diǎn)系統(tǒng)的數(shù)據(jù)文件(局部) 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 數(shù)據(jù)文件2) 5節(jié)點(diǎn)系統(tǒng)圖5 數(shù)據(jù)文件3) 4節(jié)點(diǎn)系統(tǒng)圖6 極坐標(biāo)牛頓法迭代主程序框圖電壓初始化計(jì)算殘差P、Q修正量、V，修正、V計(jì)算雅可比矩陣P 、|Q形成導(dǎo)納矩陣開始結(jié)束原始