前推回代法计算流程

1 前推回代法计算流程 要看懂前推回代法计算程序 报告叙述计算原理及计算流程 绘制计算流程框图 确定前推回 代支路次序 广度优先 或深度优先 编写前推回代计算输入文件 进行潮流计算 下列为节点配电网结构图及系统支路参数和系统负荷参数表 1 2 3 4 5 67 8 9 101112 13 图 1 2 节点配电网结构图 表 1 系统支路参数 支路R X B 2 TK S 1 23 3673 6850 0 2 32 3562 5410 0 3 41 1451 280 0 4 54 5245 040 0 2 60 8561 140 0 6 72 7452 9650 0 2 83 7434 2510 0 8 92 2372 7560 0 3 104 144 6960 0 3 111 3281 7630 0 11 122 4362 8660 0 4 133 5213 9660 0 表 2 系统负荷参数 节点编号节点类型节点初始电压 kV Pi MVA Qi MVA 2 1根节点10 400 2中间节点10 00 03420 0301 3中间节点10 00 06930 0642 4中间节点10 00 08450 0763 5叶节点10 00 02950 0261 6中间节点10 00 04740 0409 7叶节点10 00 11760 0957 8中间节点10 00 09460 0857 9叶节点10 00 09160 0859 10叶节点10 00 02710 0229 11中间节点10 00 06960 0643 12叶节点10 00 06760 0579 13叶节点10 00 02980 0242 主程序清单 PQ FT RX case114 调用数据文件 NN size PQ 1 节点数 NB size FT 1 支路数数 V PQ 1 V 初始电压相量 maxd 1 k 1 while maxd 0 0001 PQ2 PQ 每一次迭代各节点的注入有功和无功相同 PL 0 0 for i 1 NB kf FT i 1 前推始节点号 kt FT i 2 前推终节点号 x PQ2 kf 2 2 PQ2 kf 3 2 V kf V kf 计算沿线电流平方 A PQ1 i 1 PQ2 kf 2 RX i 1 x 计算支路首端有功 MW RX i 1 R PQ1 i 2 PQ2 kf 3 RX i 2 x 计算沿支路的无功损耗 Mvar RX i 2 X PQ2 kt 2 PQ2 kt 2 PQ1 i 1 用 PQ1 去修正支路末端节点的有功 P 单位 MW PQ2 kt 3 PQ2 kt 3 PQ1 i 2 用 PQ1 去修正支路末端节点的有功 Q 单位 Mvar 3 PL PL RX i 1 x end angle 1 0 0 for i NB 1 1 kf FT i 2 回代始节点号 kt FT i 1 回代终节点号 dv1 PQ1 i 1 RX i 1 PQ1 i 2 RX i 2 V kf 计算支路电压损耗的纵分量 dv1 dv2 PQ1 i 1 RX i 2 PQ1 i 2 RX i 1 V kf 计算支路电压损耗的横分量 dv2 V2 kt sqrt V kf dv1 2 dv2 2 计算支路末端电压 kV angle kt angle kf atand dv2 V kf dv1 计算支路 end maxd abs V2 2 V 2 V2 1 V 1 for i 3 1 NN if abs V2 i V i maxd maxd abs V2 i V i end end maxd k k 1 PQ1 潮流分布 即支路首端潮流 MVA V V2 节点电压模计算结果 kV angle 节点电压角度计算结果单位度 PL 网损单位 MW end clear 输入文件清单 function PQ FT RX case114 PQ 节点电压 有功 无功 10 4 0 0 4 10 0 0 0342 0 0301 10 0 0 0693 0 0642 10 0 0 0845 0 0763 10 0 0 0295 0 0261 10 0 0 0474 0 0409 10 0 0 1176 0 0957 10 0 0 0946 0 0857 10 0 0 0916 0 0859 10 0 0 0271 0 0229 10 0 0 0696 0 0643 10 0 0 0676 0 0579 10 0 0 0298 0 0242 FT 首端 末端 5 4 13 4 4 3 10 3 12 11 11 3 7 6 6 2 9 8 8 2 3 2 2 1 RX R X 4 524 5 04 3 521 3 966 1 145 1 28 4 14 4 696 2 436 2 866 1 328 1 763 2 745 2 965 0 856 1 14 2 237 2 756 3 743 4 251 2 356 2 541 3 367 3 685 计算过程 maxd 1 5 k 1 maxd 0 1780 k 2 PQ1 0 0296 0 0262 0 0299 0 0243 0 1443 0 1272 0 0272 0 0230 0 0678 0 0581 0 1378 0 1230 0 1182 0 0964 0 1660 0 1378 0 0920 0 0863 0 1890 0 1748 0 3847 0 3439 0 8099 0 7260 V Columns 1 through 8 10 4000 9 8807 9 8220 9 9672 9 9734 9 9701 9 9390 9 8550 Columns 9 through 13 9 9556 9 9780 9 9600 9 9668 9 9799 angle 6 Columns 1 through 8 0 0 3011 0 3986 0 4211 0 4387 0 3421 0 3916 0 3878 Columns 9 through 13 0 4225 0 4173 0 4444 0 4747 0 4400 PL 0 0471 maxd 0 1787 k 3 PQ1 0 0296 0 0262 0 0299 0 0243 0 1443 0 1272 0 0272 0 0230 0 0678 0 0581 0 1378 0 1230 0 1182 0 0964 0 1660 0 1378 0 0920 0 0863 0 1890 0 1749 0 3849 0 3442 0 8112 0 7274 V Columns 1 through 8 10 4000 9 8798 9 7004 9 7886 9 9405 9 8504 9 9089 9 7338 7 Columns 9 through 13 9 8100 9 7996 9 7813 9 9267 9 9470 angle Columns 1 through 8 0 0 3011 0 4011 0 4244 0 4421 0 3431 0 3929 0 3899 Columns 9 through 13 0 4257 0 4204 0 4486 0 4791 0 4434 PL 0 0484 maxd 0 1793 k 4 PQ1 0 0296 0 0262 0 0299 0 0243 0 1444 0 1272 0 0272 0 0230 0 0678 0 0581 0 1379 0 1231 0 1182 0 0964 0 1661 0 1378 0 0920 0 0864 0 1891 0 1749 0 3851 0 3444 0 8115 0 7277 8 V Columns 1 through 8 10 4000 9 8796 9 6994 9 6666 9 7614 9 8495 9 7884 9 7329 Columns 9 through 13 9 6883 9 6777 9 6591 9 7474 9 7680 angle Columns 1 through 8 0 0 3011 0 4011 0 4250 0 4433 0 3431 0 3942 0 3899 Columns 9 through 13 0 4266 0 4209 0 4498 0 4814 0 4447 PL 0 0487 maxd 0 1226 k 5 PQ1 0 0296 0 0262 0 0299 0 0243 0 1444 0 1272 0 0272 0 0230 9 0 0678 0 0581 0 1379 0 1231 0 1183 0 0964 0 1661 0 1379 0 0920 0 0864 0 1891 0 1750 0 3852 0 3444 0 8115 0 7278 V Columns 1 through 8 10 4000 9 8795 9 6991 9 6656 9 6391 9 8493 9 7875 9 7326 Columns 9 through 13 9 6873 9 6767 9 6581 9 6248 9 6457 angle Columns 1 through 8 0 0 3011 0 4011 0 4250 0 4438 0 3431 0 3942 0 3899 Columns 9 through 13 0 4266 0 4209 0 4498 0 4822 0 4452 PL 0 0487 maxd 0 0010 k 6 10 PQ1 0 0296 0 0262 0 0299 0 0243 0 1444 0 1273 0 0272 0 0230 0 0678 0 0581 0 1379 0 1231 0 1183 0 0964 0 1661 0 1379 0 0920 0 0864 0 1891 0 1750 0 3852 0 3445 0 8116 0 7278 V Columns 1 through 8 10 4000 9 8795 9 6991 9 6653 9 6381 9 8492 9 7873 9 7326 Columns 9 through 13 9 6870 9 6764 9 6579 9 6238 9 6447 angle Columns 1 through 8 0 0 3011 0 4011 0 4250 0 4438 0 3431 0 3942 0 3899 Columns 9 through 13 0 4266 0 4209 0 4498 0 4823 0 4452 PL 0 0488 11 maxd 2 6021e 004 k 7 PQ1 0 0296 0 0262 0 0299 0 0243 0 1444 0 1273 0 0272 0 0230 0 0678 0 0581 0 1379 0 1231 0 1183 0 0964 0 1661 0 1379 0 0920 0 0864 0 1891 0 1750 0 3852 0 3445 0 8116 0 7278 V Columns 1 through 8 10 4000 9 8795 9 6991 9 6652 9 6378 9 8492 9 7872 9 7326 Columns 9 through 13 9 6870 9 6764 9 6578 9 6235 9 6445 angle Columns 1 through 8 0 0 3011 0 4011 0 4250 0 4438 0 3431 0 3942 0 3899 Columns 9 through 13 12 0 4266 0 4209 0 4498 0 4823 0 4452 PL 0 0488 maxd 6 1046e 005 k 8 PQ1 0 0296 0 0262 0 0299 0 0243 0 1444 0 1273 0 0272 0 0230 0 0678 0 0581 0 1379 0 1231 0 1183 0 0964 0 1661 0 1379 0 0920 0 0864 0 1891 0 1750 0 3852 0 3445 0 8116 0 7278 V Columns 1 through 8 10 4000 9 8795 9 6991 9 6652 9 6377 9 8492 9 7872 9 7326 Columns 9 through 13 9 6870 9 6764 9 6578 9 6235 9 6444 angle 13 Columns 1 through 8 0 0 3011 0 4011 0 4250 0 4438 0 3431 0 3942 0 3899 Columns 9 through 13 0 4266 0 4209 0 4498 0 4823 0 4452 PL 0 0488 计算结果清单 maxd 1 k 1 maxd 0 1780 k 2 maxd 0 1787 k 3 maxd 0 1793 k 4 maxd 0 1226 k 5 maxd 0 0010 k 6 maxd 2 6021e 04 k 7 maxd 6 1046e 05 k 8 PQ1 0 0296 0 0262 0 0299 0 0243 0 1444 0 1273 0 0272 0 0230 0 0678 0 0581 0 1379 0 1231 0 1183 0 0964 0 1661 0 1379 0 0920 0 0864 0 1891 0 175