1. 手算过程
已知:
节点1:PQ节点, s(1)= -0.5000-j0.3500
节点2:PV节点, p(2)=0.4000 v(2)=1.0500
节点3:平衡节点,U(3)=1.0000∠0.0000
网络的连接图:
0.0500+j0.2000 1 0.0500+j0.2000
23
1)计算节点导纳矩阵
由 ;
;
导纳矩阵中的各元素:
;
; ;
; ;
;
; ;
;
形成导纳矩阵:
2)计算各PQ、PV节点功率的不平衡量,及PV节点电压的不平衡量:
取:
节点3是平衡节点,保持为定值。
;
;
;
;
;
;
于是:
3)计算雅可比矩阵中各元素
雅可比矩阵的各个元素分别为:
又:;
;
;
;
;
;
;
;
;
;
;
;
; ;
;
;
; ;
得到K=0时的雅可比矩阵:
4)建立修正方程组:
解得:
因为 ; ;
所以 ;
;
;
;
5)运用各节点电压的新值进行下一次迭代:
即取:
节点3时平衡节点,保持为定值。
6)计算各PQ、PV节点功率的不平衡量,及PV节点电压的不平衡量:
;
;
;
;
;
于是:
7)计算雅可比矩阵中各元素:
雅可比矩阵的各个元素分别为:
又:;
;
;
;
;
; ;
;
;
得到K=1时的雅可比矩阵:
8)建立修正方程组如下:
解得:
因为 ; ;
所以 ;
;
;
;
9)运用各节点电压的新值进行下一次迭代:………………………
2. 程序流程图
见电子版程序:余越.c
4. 输入数据
3,2,0,1,1,0.00010
1,1,-0.5000,-0.3500,1.0000,0.0000
2,2,0.4000,1.0500,1.0000,0.0000
3,3,1.0000,0.0000
1,1,2,0.0500,0.2000
2,1,3,0.0500,0.2000
5. 输出数据
潮流上机实习 华北电力大学 电自031 余越 200301030127
********** 原始数据 **********
===================================================================
节点数: 3 支路数: 2 对地支路数: 0 PQ节点数: 1 PV节点数: 1 精度:0.000100
-------------------------------------------------------------------
PQ节点 节点 1 P[1]=-0.500000 Q[1]=-0.350000
PV节点 节点 2 P[2]=0.400000 V[2]=1.050000
平衡节点 节点 3 e[3]=1.000000 f[3]=0.000000
-------------------------------------------------------------------
支路 1 相关节点: 1, 2 R=0.050000 X=0.200000
支路 2 相关节点: 1, 3 R=0.050000 X=0.200000
===================================================================
********* 计算结果 *********
节点导纳矩阵为:
2.35294+j-9.41177 -1.17647+j 4.70588 -1.17647+j 4.70588
-1.17647+j 4.70588 1.17647+j-4.70588 0.00000+j 0.00000
-1.17647+j 4.70588 0.00000+j 0.00000 1.17647+j-4.70588
==============================================================
===============================================================
--------------------输出雅可比矩阵和迭代次数---------------------
迭代的次数为 1
9.41177 2.35294 -4.70588 -1.17647 -0.50000
-2.35294 9.41177 1.17647 -4.70588 -0.35000
-4.70588 -1.17647 4.70588 1.17647 0.40000
0.00000 0.00000 0.00000 2.00000 0.10250
=============================================================
=====================================================================
--------------------------输出 df,de--------------------------------
节点为 1 df=-0.01578 de=-0.02188
节点为 2 df= 0.05094 de= 0.05125
---------------------------------------------------------------------
----------------输出迭代过程中的电压值-------------------------------
节点为 1 f=-0.01578 e= 0.97812
节点为 2 f= 0.05094 e= 1.05125
=====================================================================
===============================================================
--------------------输出雅可比矩阵和迭代次数---------------------
迭代的次数为 2
9.51875 1.95000 -4.58437 -1.22500 -0.00541
-2.95000 8.81875 1.22500 -4.58437 -0.01555
-5.00699 -0.99706 4.74136 1.39706 -0.00697
0.00000 0.00000 0.10188 2.10250 -0.00522
=============================================================
=====================================================================
--------------------------输出 df,de--------------------------------
节点为 1 df=-0.00188 de=-0.00312
节点为 2 df=-0.00343 de=-0.00232
---------------------------------------------------------------------
----------------输出迭代过程中的电压值-------------------------------
节点为 1 f=-0.01766 e= 0.97501
节点为 2 f= 0.04751 e= 1.04893
=====================================================================
===============================================================
--------------------输出雅可比矩阵和迭代次数---------------------
迭代的次数为 3
9.50304 1.95417 -4.56748 -1.23017 0.00001
-2.96650 8.76690 1.23017 -4.56748 -0.00007
-4.99205 -1.01046 4.72083 1.40412 -0.00003
0.00000 0.00000 0.09502 2.09787 -0.00002
=============================================================
=====================================================================
--------------------------输出 df,de--------------------------------
节点为 1 df=-0.00000 de=-0.00001
节点为 2 df=-0.00001 de=-0.00001
---------------------------------------------------------------------
----------------输出迭代过程中的电压值-------------------------------
节点为 1 f=-0.01766 e= 0.97499
节点为 2 f= 0.04750 e= 1.04892
=====================================================================
各节点电压为:
U[1]=0.974994 + j -0.017661
U[2]=1.048925 + j 0.047502
U[3]=1.000000 + j 0.000000
===================================================================
平衡节点功率为:
S[3]=0.112528+j 0.096897
===================================================================
线路功率如下:
线路1-2的功率: -0.388574 + j -0.257513
线路2-1的功率: 0.400000 + j 0.303216
线路1-3的功率: -0.111426 + j -0.092487
线路3-1的功率: 0.112528 + j 0.096897
==================================================================
线路上损耗的功率为:
线路1-2上的功率损耗: 0.011426 + j 0.045703
线路1-3上的功率损耗: 0.001103 + j 0.004410
====================================================================
网络总损耗: 0.012528 + j 0.050114
************************结束************************
6. 程序说明
1.输入文件为 input.txt
2.输出文件为 output.txt
3.本程序采用直角坐标系,牛顿-拉夫逊法;
4.在输入文件中:
1)节点数,总的支路数(包括对地支路),对地支路数,PQ节点数,PV节点数,精度;用逗号隔开
2)节点编号,节点类型:1表示PQ节点:有功,无功,电压初值;
2表示PV节点:有功,电压,电压初值;
3表示平衡节点:电压值;用逗号隔开;
3)支路编号,支路相关节点p1,p2;电阻和电抗;用逗号隔开;
(有平行支路要连续输入)
(对地支路要放在最后,若要输出对地支路的潮流则子程序还需适当的修改);
5.在输出文件中输出了本程序计算的原始数据、迭代过程及结果。
7.计算结果分析
1.由计算机所的结果:
各节点电压为:
U[1]=0.974994 + j -0.017661
U[2]=1.048925 + j 0.047502
U[3]=1.000000 + j 0.000000
平衡节点功率为:
S[3]=0.112528+j 0.096897
线路功率如下:
线路1-2的功率: -0.388574 + j -0.257513
线路2-1的功率: 0.400000 + j 0.303216
线路1-3的功率: -0.111426 + j -0.092487
线路3-1的功率: 0.112528 + j 0.096897
网络总损耗: 0.012528 + j 0.050114
2.结果分析
1)各节点电压均在额定电压附近,平衡功率中无功功率为0.050114,即发出无功。
2)所用线路有功都不超过1,没有过载,即转换功率比较大,运行较为安全。
3)线路损耗中线路1-2较大,须提高电压等级,以减少网损。
3.程序的优缺点
1)优点:该程序采用牛顿拉夫逊算法,收敛较快。解方程用列主元素消去法,一般不会产生溢出。采用输入输出文件,可多次进行调用,且不因源程序多次修改数据而出错,减少数据输入的麻烦,方便看运行结果。本程序考虑到PQ、PV节点,计算的节点数,可根据电脑运行容量确定的计算方程的阶数来限定,一般情况下都能满足运算要求。运算精度可由用户自定义输入。通过多道题目的试验,本程序均可使用,即通用性良好。
2)缺点:数据输入不美观,且输入格式要求严格。不能解决非标准变比变压器的潮流计算问题。没有建立良好的人机交互界面。
8. 实验总结
为期两周的电力系统潮流上机计算结束了,通过这两周的程序设计,手算题目等科目使我对运用计算机辅助计算电力系统潮流分布的原理和方法有了更深一步的了解,对电力系统潮流分布的概念有了更深刻的理解,并且通过程序来完成自己的理解。
这套程序用VC++语言编程解决潮流计算问题,采用牛顿拉夫逊算法,计算速度快,收敛性比较好。生成的数据记录包括过程和结果,查看的时候比较方便,但是输入的格式比较死板,程序解决的问题还是比较简单,由于个人水平有限,所以没有去尝试较为复杂的潮流分布问题。由于没有选数值计算这门课程,编程期间自学了一些这方面的知识,但是还是存在着缺陷,或多或少的反映在计算和编程的过程中,程序中难免有不完善的地方,有待补充完善。
最后还是要感谢各位老师的悉心指导,我们的问题各种各样,有时候自己也会被自己的程序弄的很烦,但是各位老师还是做出了耐心和富有建设性的指导,才使我们能顺利完成这次的程序设计,再次感谢!
学院专业班级学号姓名指导老师完成时间电力系统自动化报告核技术与自动化学院电气工程及其自动化20xx060505班320xx0605…
电力系统实验报告班级学号姓名实验一同步发电机准同期并列实验一实验目的1加深理解同步发电机准同期并列原理掌握准同期并列条件2掌握微机…
单机无穷大系统稳态实验一整理实验数据说明单回路送电和双回路送电对电力系统稳定运行的影响并对实验结果进行理论分析实验数据如下由实验数…
重庆邮电大学移通学院自动化系电力系统自动化及其微机保护实验专业电气工程及其自动化班级05111001学号0511100104姓名刘…
电力系统分析实验报告实验项目1单机无穷大系统稳态运行实验2电力系统暂态稳定实验学院电气信息学院专业电气工程及其自动化班级20xx级…
网页设计实验报告院部热能学院专业热能与动力工程班级112姓名范金仓学号20xx031388一实验目的及要求1确定网站主题和网站的用…
单片机实验报告班级041231学号OO姓名雷锋实验一数码管实验一实验目的1了解数码管的显示原理2掌握JXARM92440中数码管显…
重庆交通大学学生实验报告实验课程名称计算机辅助制造开课实验室重庆交通大学计算机中心学院国际学院20xx级专业机械设计制造及其自动化…
C程序实验报告实验五继承与派生实验目的1学习定义和使用类的继承关系定义派生类2熟悉不同继承方式下对基类成员的访问控制3学习利用虚基…