Well in the given data shunt reactors are already compensated. I was able to run newton’s power flow. Thanks for help.
Raza Farid Electrical Design Engineer Green Building Technologies Contact: +92-333-5615025 [email protected] From: Fred Kanjelesa Sent: Monday, December 16, 2013 11:03 AM To: MATPOWER discussion forum Hello, I had experienced a similar problem when i tried to run a 62 bus system for our national system, until i incorporated voltage compensation equipment. owing to the size of the network, you may consider ensuring you have included the voltage compensation equipment (reactors and capacitors) at strategic load buses as it may be in the actual system. In case you need information on how to do this, for reactors +Q is used and for capacitors -Q is used, this is so because these are placed at the load bus and to include these in MATPOWER all you need is to add the (+Q for reactor OR -Q for capacitor) to the reactive part of the load and that's it. Hope it works out for you. Regards On 13 December 2013 19:02, Raza Farid <[email protected]> wrote: Hello, I have been working on an assignment and I can not figure out what the problem is. I checked my system and readings several times but still my system isn't converging. Any help will be appreciated. %%----- Power Flow Data -----%% %% system MVA base mpc.baseMVA = 100; %% bus data % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin mpc.bus = [ 1 3 0 0 0 0 1 1 0 500 1 1.06 0.94; 2 1 159 78 0 0 1 1 0 220 1 1.06 0.94; 3 1 172 86 0 0 1 1 0 220 1 1.06 0.94; 4 2 0 0 0 0 1 1 0 220 1 1.06 0.94; 5 1 170 50 0 0 1 1 0 220 1 1.06 0.94; 6 2 0 0 0 0 1 1 0 500 1 1.06 0.94; 7 1 286 52 0 0 1 1 0 220 1 1.06 0.94; 8 1 372 186 0 0 1 1 0 220 1 1.06 0.94; 9 1 206 27 0 0 1 1 0 220 1 1.06 0.94; 10 1 4.8 0 0 0 1 1 0 220 1 1.06 0.94; 11 1 -20 -20 0 0 1 1 0 220 1 1.06 0.94; 12 2 0 0 0 -222 1 1 0 500 1 1.06 0.94; 13 1 122 61 0 0 1 1 0 220 1 1.06 0.94; 14 1 148 49 0 0 1 1 0 220 1 1.06 0.94; 15 1 0 0 0 -222 1 1 0 500 1 1.06 0.94; 16 2 104 0 0 0 1 1 0 220 1 1.06 0.94; 17 1 0 0 0 -222 1 1 0 500 1 1.06 0.94; 18 1 180 80 0 0 1 1 0 220 1 1.06 0.94; 19 1 105 50 0 0 1 1 0 220 1 1.06 0.94; 20 1 191 40 0 0 1 1 0 220 1 1.06 0.94; 21 1 247 170 0 0 1 1 0 220 1 1.06 0.94; 22 2 0 0 0 0 1 1 0 220 1 1.06 0.94; 23 1 0 0 0 -282 1 1 0 500 1 1.06 0.94; 24 1 191 40 0 0 1 1 0 220 1 1.06 0.94; 25 1 165 60 0 0 1 1 0 220 1 1.06 0.94; 26 1 443 180 0 0 1 1 0 220 1 1.06 0.94; 27 1 86 100 0 0 1 1 0 220 1 1.06 0.94; 28 2 0 0 0 0 1 1 0 220 1 1.06 0.94; 29 1 0 0 0 -222 1 1 0 500 1 1.06 0.94; 30 1 0 0 0 0 1 1 0 220 1 1.06 0.94; 31 1 255 24 0 0 1 1 0 220 1 1.06 0.94; 32 1 0 0 0 -510 1 1 0 500 1 1.06 0.94; 33 1 296 114 0 0 1 1 0 220 1 1.06 0.94; 34 1 289 130 0 0 1 1 0 220 1 1.06 0.94; 35 1 137 50 0 0 1 1 0 220 1 1.06 0.94; 36 1 239 0 0 0 1 1 0 220 1 1.06 0.94; 37 1 176 120 0 0 1 1 0 220 1 1.06 0.94; 38 2 0 0 0 0 1 1 0 500 1 1.06 0.94; 39 1 0 0 0 -387 1 1 0 500 1 1.06 0.94; 40 1 194 65 0 -544 1 1 0 220 1 1.06 0.94; 41 1 110 44 0 0 1 1 0 220 1 1.06 0.94; 42 2 0 0 0 0 1 1 0 220 1 1.06 0.94; 43 2 0 0 0 0 1 1 0 220 1 1.06 0.94; 44 1 0 0 0 -222 1 1 0 500 1 1.06 0.94; 45 2 0 0 0 0 1 1 0 220 1 1.06 0.94; 46 2 0 0 0 0 1 1 0 220 1 1.06 0.94; 47 2 129 89 0 0 1 1 0 220 1 1.06 0.94; 48 1 0 0 0 -555 1 1 0 500 1 1.06 0.94; 49 1 76 8 0 0 1 1 0 220 1 1.06 0.94; 50 1 207 70 0 0 1 1 0 220 1 1.06 0.94; 51 1 0 0 0 -354 1 1 0 500 1 1.06 0.94; 52 1 10 26 0 0 1 1 0 220 1 1.06 0.94; 53 2 0 0 0 -297 1 1 0 500 1 1.06 0.94; 54 2 0 0 0 0 1 1 0 220 1 1.06 0.94; 55 2 0 0 0 0 1 1 0 220 1 1.06 0.94; 56 2 0 0 0 0 1 1 0 220 1 1.06 0.94; 57 1 550 120 0 0 1 1 0 500 1 1.06 0.94; 58 1 146 20 0 0 1 1 0 220 1 1.06 0.94; 59 1 135 -38 0 0 1 1 0 220 1 1.06 0.94; 60 1 122 60 0 0 1 1 0 220 1 1.06 0.94; 61 1 0 0 0 -264 1 1 0 500 1 1.06 0.94; 62 1 230 45 0 0 1 1 0 220 1 1.06 0.94; 63 1 157 90 0 0 1 1 0 220 1 1.06 0.94; 64 1 203 60 0 0 1 1 0 220 1 1.06 0.94; 65 1 221 79 0 0 1 1 0 220 1 1.06 0.94; 66 1 353 140 0 0 1 1 0 220 1 1.06 0.94; 67 1 0 0 0 0 1 1 0 220 1 1.06 0.94; 68 1 0 0 0 0 1 1 0 220 1 1.06 0.94; 69 2 144 200 0 0 1 1 0 220 1 1.06 0.94; 70 2 122 45 0 0 1 1 0 220 1 1.06 0.94; 71 1 24.84 5.09 0 0 1 1 0 220 1 1.06 0.94; 72 1 135 -38 0 0 1 1 0 220 1 1.06 0.94; ]; %% generator data % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf mpc.gen = [ 1 0 0 10 0 0 100 1 100 0 0 0 0 0 0 0 0 0 0 0 0; 4 607 22 10 0 0 100 1 707 0 0 0 0 0 0 0 0 0 0 0 0; 6 1530 335 10 0 0 100 1 1630 0 0 0 0 0 0 0 0 0 0 0 0; 12 870 120 10 0 0 100 1 970 0 0 0 0 0 0 0 0 0 0 0 0; 16 881 0 10 0 0 100 1 981 0 0 0 0 0 0 0 0 0 0 0 0; 22 208 50 10 0 0 100 1 308 0 0 0 0 0 0 0 0 0 0 0 0; 28 471 135 10 0 0 100 1 571 0 0 0 0 0 0 0 0 0 0 0 0; 38 416 57 10 0 0 100 1 516 0 0 0 0 0 0 0 0 0 0 0 0; 42 175 15 10 0 0 100 1 275 0 0 0 0 0 0 0 0 0 0 0 0; 43 260 110 10 0 0 100 1 360 0 0 0 0 0 0 0 0 0 0 0 0; 45 343 -144 10 0 0 100 1 443 0 0 0 0 0 0 0 0 0 0 0 0; 46 156 -132 10 0 0 100 1 256 0 0 0 0 0 0 0 0 0 0 0 0; 47 821 89 10 0 0 100 1 921 0 0 0 0 0 0 0 0 0 0 0 0; 53 855 -250 10 0 0 100 1 955 0 0 0 0 0 0 0 0 0 0 0 0; 54 174 -2 10 0 0 100 1 274 0 0 0 0 0 0 0 0 0 0 0 0; 55 0 0 10 0 0 100 1 100 0 0 0 0 0 0 0 0 0 0 0 0; 56 549 -82 10 0 0 100 1 649 0 0 0 0 0 0 0 0 0 0 0 0; 69 677 24 10 0 0 100 1 777 0 0 0 0 0 0 0 0 0 0 0 0; 70 340 70 10 0 0 100 1 440 0 0 0 0 0 0 0 0 0 0 0 0; ]; %% branch data % fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax mpc.branch = [ 1 13 0 0.00927 0 9900 0 0 0 0 1 -360 360; 6 28 0 0.00135 0 9900 0 0 0 0 1 -360 360; 12 30 0 0.01862 0 9900 0 0 0 0 1 -360 360; 15 63 0 0.00927 0 9900 0 0 0 0 1 -360 360; 17 64 0 0.00931 0 9900 0 0 0 0 1 -360 360; 23 65 0 1e-006 0 9900 0 0 0 0 1 -360 360; 32 67 0 5e-006 0 9900 0 0 0 0 1 -360 360; 29 66 0 0.009312 0 9900 0 0 0 0 1 -360 360; 39 68 0 0.01456 0 9900 0 0 0 0 1 -360 360; 44 42 0 0.01862 0 9900 0 0 0 0 1 -360 360; 44 43 0 0.01862 0 9900 0 0 0 0 1 -360 360; 48 69 0 0.0927 0 9900 0 0 0 0 1 -360 360; 51 72 0 0.02822 0 9900 0 0 0 0 1 -360 360; 61 70 0 0.02822 0 9900 0 0 0 0 1 -360 360; 6 12 0.00015 0.01451 1.3632 9900 0 0 0 0 1 -360 360; 6 1 0.00099 0.01195 1.3153 9900 0 0 0 0 1 -360 360; 6 15 0.00098 0.01131 1.248 9900 0 0 0 0 1 -360 360; 12 32 0.0034 0.03387 3.2119 9900 0 0 0 0 1 -360 360; 12 15 0.00097 0.0111 1.2255 9900 0 0 0 0 1 -360 360; 15 17 0.0017 0.01956 2.1588 9900 0 0 0 0 1 -360 360; 17 23 0.00559 0.00641 0.7072 9900 0 0 0 0 1 -360 360; 23 29 0.00135 0.01625 1.7875 9900 0 0 0 0 1 -360 360; 23 32 0.00083 0.00957 1.0569 9900 0 0 0 0 1 -360 360; 29 39 0.00144 0.0165 1.8215 9900 0 0 0 0 1 -360 360; 32 39 0.00197 0.02262 2.4961 9900 0 0 0 0 1 -360 360; 32 44 0.00249 0.02853 3.1482 9900 0 0 0 0 1 -360 360; 32 38 0.00172 0.01649 1.5491 9900 0 0 0 0 1 -360 360; 39 44 0.00056 0.00648 0.7151 9900 0 0 0 0 1 -360 360; 39 48 0.00278 0.03189 3.5193 9900 0 0 0 0 1 -360 360; 44 48 0.0028 0.02815 2.6696 9900 0 0 0 0 1 -360 360; 48 51 0.00254 0.02914 3.2157 9900 0 0 0 0 1 -360 360; 51 61 0.0017 0.0167 1.581 9900 0 0 0 0 1 -360 360; 53 61 0.00153 0.01849 2.0348 9900 0 0 0 0 1 -360 360; 53 57 0.00138 0.01589 1.754 9900 0 0 0 0 1 -360 360; 57 61 0.00132 0.01594 1.7537 9900 0 0 0 0 1 -360 360; 38 39 0.00057 0.00662 0.7308 9900 0 0 0 0 1 -360 360; 3 13 0.01894 0.10184 0.1701 9900 0 0 0 0 1 -360 360; 4 3 0.00577 0.04432 0.1441 9900 0 0 0 0 1 -360 360; 5 3 0.00654 0.05023 0.1633 9900 0 0 0 0 1 -360 360; 2 13 0.01232 0.06624 0.1106 9900 0 0 0 0 1 -360 360; 30 13 0.01389 0.10667 0.3469 9900 0 0 0 0 1 -360 360; 7 28 0.00539 0.02898 0.0484 9900 0 0 0 0 1 -360 360; 8 2 0.01386 0.07452 0.1247 9900 0 0 0 0 1 -360 360; 9 28 0.01124 0.06044 0.1009 9900 0 0 0 0 1 -360 360; 8 28 0.01047 0.0563 0.094 9900 0 0 0 0 1 -360 360; 8 30 0.01108 0.06044 0.1009 9900 0 0 0 0 1 -360 360; 9 7 0.00585 0.03146 0.0525 9900 0 0 0 0 1 -360 360; 63 14 0.00354 0.02718 0.0884 9900 0 0 0 0 1 -360 360; 63 10 0.00215 0.01654 0.0538 9900 0 0 0 0 1 -360 360; 63 9 0.00585 0.03146 0.0525 9900 0 0 0 0 1 -360 360; 16 63 0.01232 0.06624 0.1106 9900 0 0 0 0 1 -360 360; 10 9 0.00215 0.01654 0.0538 9900 0 0 0 0 1 -360 360; 11 9 0.01389 0.10667 0.3469 9900 0 0 0 0 1 -360 360; 64 16 0.00885 0.06796 0.221 9900 0 0 0 0 1 -360 360; 18 16 0.01663 0.08942 0.1493 9900 0 0 0 0 1 -360 360; 18 64 0.00462 0.03546 0.1153 9900 0 0 0 0 1 -360 360; 19 18 0.0077 0.0414 0.0691 9900 0 0 0 0 1 -360 360; 20 16 0.02556 0.13744 0.2295 9900 0 0 0 0 1 -360 360; 20 19 0.0231 0.01242 0.2074 9900 0 0 0 0 1 -360 360; 21 20 0.00184 0.01418 0.0461 9900 0 0 0 0 1 -360 360; 65 21 0.00291 0.02239 0.0728 9900 0 0 0 0 1 -360 360; 24 65 0.00261 0.01407 0.0235 9900 0 0 0 0 1 -360 360; 24 26 0.00182 0.01401 0.0455 9900 0 0 0 0 1 -360 360; 25 20 0.00415 0.02235 0.0373 9900 0 0 0 0 1 -360 360; 21 22 0.00385 0.0207 0.0345 9900 0 0 0 0 1 -360 360; 25 65 0.00415 0.02235 0.0373 9900 0 0 0 0 1 -360 360; 25 65 0.00207 0.01595 0.0518 9900 0 0 0 0 1 -360 360; 26 25 0.00261 0.01407 0.0235 9900 0 0 0 0 1 -360 360; 26 65 0.00446 0.03427 0.1114 9900 0 0 0 0 1 -360 360; 22 65 0.00061 0.00472 0.0153 9900 0 0 0 0 1 -360 360; 27 26 0.00346 0.02659 0.0864 9900 0 0 0 0 1 -360 360; 67 66 0.01447 0.07783 0.13 9900 0 0 0 0 1 -360 360; 67 31 0.0154 0.0828 0.0138 9900 0 0 0 0 1 -360 360; 33 66 0.00962 0.07387 0.2402 9900 0 0 0 0 1 -360 360; 34 67 0.00154 0.00828 0.0138 9900 0 0 0 0 1 -360 360; 35 67 0.00046 0.00248 0.0041 9900 0 0 0 0 1 -360 360; 36 35 0.00539 0.02898 0.0484 9900 0 0 0 0 1 -360 360; 67 20 0.02464 0.18912 0.615 9900 0 0 0 0 1 -360 360; 31 4 0.0154 0.1182 0.3844 9900 0 0 0 0 1 -360 360; 68 36 0.01001 0.05382 0.0898 9900 0 0 0 0 1 -360 360; 68 37 0.00215 0.01159 0.0193 9900 0 0 0 0 1 -360 360; 68 33 0.00708 0.05437 0.1768 9900 0 0 0 0 1 -360 360; 41 68 0.00426 0.03274 0.1064 9900 0 0 0 0 1 -360 360; 42 68 0.00819 0.04405 0.0735 9900 0 0 0 0 1 -360 360; 42 40 0.00831 0.06382 0.2075 9900 0 0 0 0 1 -360 360; 43 42 6e-005 0.0005 0.0016 9900 0 0 0 0 1 -360 360; 42 71 0.00369 0.01987 0.0331 9900 0 0 0 0 1 -360 360; 43 71 0.00492 0.02649 0.0442 9900 0 0 0 0 1 -360 360; 45 43 0.00138 0.01063 0.0345 9900 0 0 0 0 1 -360 360; 47 42 0.00805 0.0433 0.0723 9900 0 0 0 0 1 -360 360; 43 46 0.00138 0.01063 0.0345 9900 0 0 0 0 1 -360 360; 47 46 0.01001 0.05382 0.0898 9900 0 0 0 0 1 -360 360; 43 68 0.00431 0.03309 0.1076 9900 0 0 0 0 1 -360 360; 43 41 6e-005 0.0005 0.0016 9900 0 0 0 0 1 -360 360; 47 68 0.0154 0.0828 0.1383 9900 0 0 0 0 1 -360 360; 47 68 0.005 0.03841 0.1249 9900 0 0 0 0 1 -360 360; 50 69 0.01139 0.06127 0.1023 9900 0 0 0 0 1 -360 360; 56 50 0.02494 0.13413 0.224 9900 0 0 0 0 1 -360 360; 70 59 0.00105 0.00809 0.0263 9900 0 0 0 0 1 -360 360; 70 60 0.00639 0.03726 0.0622 9900 0 0 0 0 1 -360 360; 62 58 0.02156 0.11592 0.1936 9900 0 0 0 0 1 -360 360; 58 56 0.0215 0.13496 0.2254 9900 0 0 0 0 1 -360 360; 58 69 0.03988 0.21445 0.3581 9900 0 0 0 0 1 -360 360; 55 52 0.00562 0.03063 0.0511 9900 0 0 0 0 1 -360 360; 55 54 0.00092 0.00049 0.0082 9900 0 0 0 0 1 -360 360; 52 54 0.00184 0.00993 0.0165 9900 0 0 0 0 1 -360 360; 55 49 0.00431 0.02318 0.0387 9900 0 0 0 0 1 -360 360; 50 55 0.00631 0.03339 0.0567 9900 0 0 0 0 1 -360 360; 69 56 0.00739 0.05673 0.1845 9900 0 0 0 0 1 -360 360; ]; %%----- OPF Data -----%% %% generator cost data % 1 startup shutdown n x1 y1 ... xn yn % 2 startup shutdown n c(n-1) ... c0 mpc.gencost = [ 2 0 0 3 0.01 40 0; 2 0 0 3 0.0164744646 20 0; 2 0 0 3 0.00653594771 20 0; 2 0 0 3 0.0114942529 20 0; 2 0 0 3 0.0113507378 20 0; 2 0 0 3 0.0480769231 20 0; 2 0 0 3 0.0212314225 20 0; 2 0 0 3 0.0240384615 20 0; 2 0 0 3 0.0571428571 20 0; 2 0 0 3 0.0384615385 20 0; 2 0 0 3 0.029154519 20 0; 2 0 0 3 0.0641025641 20 0; 2 0 0 3 0.012180268 20 0; 2 0 0 3 0.0116959064 20 0; 2 0 0 3 0.0574712644 20 0; 2 0 0 3 0.01 40 0; 2 0 0 3 0.0182149362 20 0; 2 0 0 3 0.0147710487 20 0; 2 0 0 3 0.0294117647 20 0; ]; Raza Farid Electrical Design Engineer Green Building Technologies Contact: +92-333-5615025 [email protected] --_704635C2-3DD4-4E87-AC9C-1FEE07C8828C_ Content-Transfer-Encoding: base64 Content-Type: text/html; charset="utf-8" <html> <head> <meta name="generator" content="Windows Mail 17.5.9600.20315"> <style data-externalstyle="true"><!-- p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph { margin-top:0in; margin-right:0in; margin-bottom:0in; margin-left:.5in; margin-bottom:.0001pt; } p.MsoNormal, li.MsoNormal, div.MsoNormal { margin:0in; margin-bottom:.0001pt; } p.MsoListParagraphCxSpFirst, li.MsoListParagraphCxSpFirst, div.MsoListParagraphCxSpFirst, p.MsoListParagraphCxSpMiddle, li.MsoListParagraphCxSpMiddle, div.MsoListParagraphCxSpMiddle, p.MsoListParagraphCxSpLast, li.MsoListParagraphCxSpLast, div.MsoListParagraphCxSpLast { margin-top:0in; margin-right:0in; margin-bottom:0in; margin-left:.5in; margin-bottom:.0001pt; line-height:115%; } --></style></head> <body dir="ltr"> <div data-externalstyle="false" dir="ltr" style="font-family: 'Calibri', 'Segoe UI', 'Meiryo', 'Microsoft YaHei UI', 'Microsoft JhengHei UI', 'Malgun Gothic', 'sans-serif';font-size:12pt;"><div style="color: rgb(0, 0, 0);">Hello,</div><div style="color: rgb(0, 0, 0);">I have been working on an assignment and I can not figure out what the problem is. I checked my system and readings several times but still my system isn't converging. Any help will be appreciated.<br></div><div style="color: rgb(0, 0, 0);" data-signatureblock="true"><div style="color: rgb(0, 0, 0);"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><p>%%----- Power Flow Data -----%%</p> <p>%% system MVA base</p> </font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2"><p>mpc.baseMVA = 100;</p> <p> </p> </font></font><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><p>%% bus data</p> <p>% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin</p> </font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2"><p>mpc.bus = [</p> <p> 1 3 0 0 0 0 1 1 0 500 1 1.06 0.94;</p> <p> 2 1 159 78 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 3 1 172 86 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 4 2 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 5 1 170 50 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 6 2 0 0 0 0 1 1 0 500 1 1.06 0.94;</p> <p> 7 1 286 52 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 8 1 372 186 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 9 1 206 27 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 10 1 4.8 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 11 1 -20 -20 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 12 2 0 0 0 -222 1 1 0 500 1 1.06 0.94;</p> <p> 13 1 122 61 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 14 1 148 49 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 15 1 0 0 0 -222 1 1 0 500 1 1.06 0.94;</p> <p> 16 2 104 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 17 1 0 0 0 -222 1 1 0 500 1 1.06 0.94;</p> <p> 18 1 180 80 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 19 1 105 50 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 20 1 191 40 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 21 1 247 170 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 22 2 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 23 1 0 0 0 -282 1 1 0 500 1 1.06 0.94;</p> <p> 24 1 191 40 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 25 1 165 60 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 26 1 443 180 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 27 1 86 100 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 28 2 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 29 1 0 0 0 -222 1 1 0 500 1 1.06 0.94;</p> <p> 30 1 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 31 1 255 24 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 32 1 0 0 0 -510 1 1 0 500 1 1.06 0.94;</p> <p> 33 1 296 114 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 34 1 289 130 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 35 1 137 50 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 36 1 239 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 37 1 176 120 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 38 2 0 0 0 0 1 1 0 500 1 1.06 0.94;</p> <p> 39 1 0 0 0 -387 1 1 0 500 1 1.06 0.94;</p> <p> 40 1 194 65 0 -544 1 1 0 220 1 1.06 0.94;</p> <p> 41 1 110 44 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 42 2 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 43 2 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 44 1 0 0 0 -222 1 1 0 500 1 1.06 0.94;</p> <p> 45 2 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 46 2 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 47 2 129 89 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 48 1 0 0 0 -555 1 1 0 500 1 1.06 0.94;</p> <p> 49 1 76 8 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 50 1 207 70 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 51 1 0 0 0 -354 1 1 0 500 1 1.06 0.94;</p> <p> 52 1 10 26 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 53 2 0 0 0 -297 1 1 0 500 1 1.06 0.94;</p> <p> 54 2 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 55 2 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 56 2 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 57 1 550 120 0 0 1 1 0 500 1 1.06 0.94;</p> <p> 58 1 146 20 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 59 1 135 -38 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 60 1 122 60 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 61 1 0 0 0 -264 1 1 0 500 1 1.06 0.94;</p> <p> 62 1 230 45 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 63 1 157 90 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 64 1 203 60 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 65 1 221 79 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 66 1 353 140 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 67 1 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 68 1 0 0 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 69 2 144 200 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 70 2 122 45 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 71 1 24.84 5.09 0 0 1 1 0 220 1 1.06 0.94;</p> <p> 72 1 135 -38 0 0 1 1 0 220 1 1.06 0.94;</p> <p>];</p> <p> </p> </font></font><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><p>%% generator data</p> <p>% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf</p> </font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2"><p>mpc.gen = [</p> <p> 1 0 0 10 0 0 100 1 100 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 4 607 22 10 0 0 100 1 707 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 6 1530 335 10 0 0 100 1 1630 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 12 870 120 10 0 0 100 1 970 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 16 881 0 10 0 0 100 1 981 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 22 208 50 10 0 0 100 1 308 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 28 471 135 10 0 0 100 1 571 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 38 416 57 10 0 0 100 1 516 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 42 175 15 10 0 0 100 1 275 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 43 260 110 10 0 0 100 1 360 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 45 343 -144 10 0 0 100 1 443 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 46 156 -132 10 0 0 100 1 256 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 47 821 89 10 0 0 100 1 921 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 53 855 -250 10 0 0 100 1 955 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 54 174 -2 10 0 0 100 1 274 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 55 0 0 10 0 0 100 1 100 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 56 549 -82 10 0 0 100 1 649 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 69 677 24 10 0 0 100 1 777 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p> 70 340 70 10 0 0 100 1 440 0 0 0 0 0 0 0 0 0 0 0 0;</p> <p>];</p> <p> </p> </font></font><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><p>%% branch data</p> <p>% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax</p> </font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2"><p>mpc.branch = [</p> <p> 1 13 0 0.00927 0 9900 0 0 0 0 1 -360 360;</p> <p> 6 28 0 0.00135 0 9900 0 0 0 0 1 -360 360;</p> <p> 12 30 0 0.01862 0 9900 0 0 0 0 1 -360 360;</p> <p> 15 63 0 0.00927 0 9900 0 0 0 0 1 -360 360;</p> <p> 17 64 0 0.00931 0 9900 0 0 0 0 1 -360 360;</p> <p> 23 65 0 1e-006 0 9900 0 0 0 0 1 -360 360;</p> <p> 32 67 0 5e-006 0 9900 0 0 0 0 1 -360 360;</p> <p> 29 66 0 0.009312 0 9900 0 0 0 0 1 -360 360;</p> <p> 39 68 0 0.01456 0 9900 0 0 0 0 1 -360 360;</p> <p> 44 42 0 0.01862 0 9900 0 0 0 0 1 -360 360;</p> <p> 44 43 0 0.01862 0 9900 0 0 0 0 1 -360 360;</p> <p> 48 69 0 0.0927 0 9900 0 0 0 0 1 -360 360;</p> <p> 51 72 0 0.02822 0 9900 0 0 0 0 1 -360 360;</p> <p> 61 70 0 0.02822 0 9900 0 0 0 0 1 -360 360;</p> <p> 6 12 0.00015 0.01451 1.3632 9900 0 0 0 0 1 -360 360;</p> <p> 6 1 0.00099 0.01195 1.3153 9900 0 0 0 0 1 -360 360;</p> <p> 6 15 0.00098 0.01131 1.248 9900 0 0 0 0 1 -360 360;</p> <p> 12 32 0.0034 0.03387 3.2119 9900 0 0 0 0 1 -360 360;</p> <p> 12 15 0.00097 0.0111 1.2255 9900 0 0 0 0 1 -360 360;</p> <p> 15 17 0.0017 0.01956 2.1588 9900 0 0 0 0 1 -360 360;</p> <p> 17 23 0.00559 0.00641 0.7072 9900 0 0 0 0 1 -360 360;</p> <p> 23 29 0.00135 0.01625 1.7875 9900 0 0 0 0 1 -360 360;</p> <p> 23 32 0.00083 0.00957 1.0569 9900 0 0 0 0 1 -360 360;</p> <p> 29 39 0.00144 0.0165 1.8215 9900 0 0 0 0 1 -360 360;</p> <p> 32 39 0.00197 0.02262 2.4961 9900 0 0 0 0 1 -360 360;</p> <p> 32 44 0.00249 0.02853 3.1482 9900 0 0 0 0 1 -360 360;</p> <p> 32 38 0.00172 0.01649 1.5491 9900 0 0 0 0 1 -360 360;</p> <p> 39 44 0.00056 0.00648 0.7151 9900 0 0 0 0 1 -360 360;</p> <p> 39 48 0.00278 0.03189 3.5193 9900 0 0 0 0 1 -360 360;</p> <p> 44 48 0.0028 0.02815 2.6696 9900 0 0 0 0 1 -360 360;</p> <p> 48 51 0.00254 0.02914 3.2157 9900 0 0 0 0 1 -360 360;</p> <p> 51 61 0.0017 0.0167 1.581 9900 0 0 0 0 1 -360 360;</p> <p> 53 61 0.00153 0.01849 2.0348 9900 0 0 0 0 1 -360 360;</p> <p> 53 57 0.00138 0.01589 1.754 9900 0 0 0 0 1 -360 360;</p> <p> 57 61 0.00132 0.01594 1.7537 9900 0 0 0 0 1 -360 360;</p> <p> 38 39 0.00057 0.00662 0.7308 9900 0 0 0 0 1 -360 360;</p> <p> 3 13 0.01894 0.10184 0.1701 9900 0 0 0 0 1 -360 360;</p> <p> 4 3 0.00577 0.04432 0.1441 9900 0 0 0 0 1 -360 360;</p> <p> 5 3 0.00654 0.05023 0.1633 9900 0 0 0 0 1 -360 360;</p> <p> 2 13 0.01232 0.06624 0.1106 9900 0 0 0 0 1 -360 360;</p> <p> 30 13 0.01389 0.10667 0.3469 9900 0 0 0 0 1 -360 360;</p> <p> 7 28 0.00539 0.02898 0.0484 9900 0 0 0 0 1 -360 360;</p> <p> 8 2 0.01386 0.07452 0.1247 9900 0 0 0 0 1 -360 360;</p> <p> 9 28 0.01124 0.06044 0.1009 9900 0 0 0 0 1 -360 360;</p> <p> 8 28 0.01047 0.0563 0.094 9900 0 0 0 0 1 -360 360;</p> <p> 8 30 0.01108 0.06044 0.1009 9900 0 0 0 0 1 -360 360;</p> <p> 9 7 0.00585 0.03146 0.0525 9900 0 0 0 0 1 -360 360;</p> <p> 63 14 0.00354 0.02718 0.0884 9900 0 0 0 0 1 -360 360;</p> <p> 63 10 0.00215 0.01654 0.0538 9900 0 0 0 0 1 -360 360;</p> <p> 63 9 0.00585 0.03146 0.0525 9900 0 0 0 0 1 -360 360;</p> <p> 16 63 0.01232 0.06624 0.1106 9900 0 0 0 0 1 -360 360;</p> <p> 10 9 0.00215 0.01654 0.0538 9900 0 0 0 0 1 -360 360;</p> <p> 11 9 0.01389 0.10667 0.3469 9900 0 0 0 0 1 -360 360;</p> <p> 64 16 0.00885 0.06796 0.221 9900 0 0 0 0 1 -360 360;</p> <p> 18 16 0.01663 0.08942 0.1493 9900 0 0 0 0 1 -360 360;</p> <p> 18 64 0.00462 0.03546 0.1153 9900 0 0 0 0 1 -360 360;</p> <p> 19 18 0.0077 0.0414 0.0691 9900 0 0 0 0 1 -360 360;</p> <p> 20 16 0.02556 0.13744 0.2295 9900 0 0 0 0 1 -360 360;</p> <p> 20 19 0.0231 0.01242 0.2074 9900 0 0 0 0 1 -360 360;</p> <p> 21 20 0.00184 0.01418 0.0461 9900 0 0 0 0 1 -360 360;</p> <p> 65 21 0.00291 0.02239 0.0728 9900 0 0 0 0 1 -360 360;</p> <p> 24 65 0.00261 0.01407 0.0235 9900 0 0 0 0 1 -360 360;</p> <p> 24 26 0.00182 0.01401 0.0455 9900 0 0 0 0 1 -360 360;</p> <p> 25 20 0.00415 0.02235 0.0373 9900 0 0 0 0 1 -360 360;</p> <p> 21 22 0.00385 0.0207 0.0345 9900 0 0 0 0 1 -360 360;</p> <p> 25 65 0.00415 0.02235 0.0373 9900 0 0 0 0 1 -360 360;</p> <p> 25 65 0.00207 0.01595 0.0518 9900 0 0 0 0 1 -360 360;</p> <p> 26 25 0.00261 0.01407 0.0235 9900 0 0 0 0 1 -360 360;</p> <p> 26 65 0.00446 0.03427 0.1114 9900 0 0 0 0 1 -360 360;</p> <p> 22 65 0.00061 0.00472 0.0153 9900 0 0 0 0 1 -360 360;</p> <p> 27 26 0.00346 0.02659 0.0864 9900 0 0 0 0 1 -360 360;</p> <p> 67 66 0.01447 0.07783 0.13 9900 0 0 0 0 1 -360 360;</p> <p> 67 31 0.0154 0.0828 0.0138 9900 0 0 0 0 1 -360 360;</p> <p> 33 66 0.00962 0.07387 0.2402 9900 0 0 0 0 1 -360 360;</p> <p> 34 67 0.00154 0.00828 0.0138 9900 0 0 0 0 1 -360 360;</p> <p> 35 67 0.00046 0.00248 0.0041 9900 0 0 0 0 1 -360 360;</p> <p> 36 35 0.00539 0.02898 0.0484 9900 0 0 0 0 1 -360 360;</p> <p> 67 20 0.02464 0.18912 0.615 9900 0 0 0 0 1 -360 360;</p> <p> 31 4 0.0154 0.1182 0.3844 9900 0 0 0 0 1 -360 360;</p> <p> 68 36 0.01001 0.05382 0.0898 9900 0 0 0 0 1 -360 360;</p> <p> 68 37 0.00215 0.01159 0.0193 9900 0 0 0 0 1 -360 360;</p> <p> 68 33 0.00708 0.05437 0.1768 9900 0 0 0 0 1 -360 360;</p> <p> 41 68 0.00426 0.03274 0.1064 9900 0 0 0 0 1 -360 360;</p> <p> 42 68 0.00819 0.04405 0.0735 9900 0 0 0 0 1 -360 360;</p> <p> 42 40 0.00831 0.06382 0.2075 9900 0 0 0 0 1 -360 360;</p> <p> 43 42 6e-005 0.0005 0.0016 9900 0 0 0 0 1 -360 360;</p> <p> 42 71 0.00369 0.01987 0.0331 9900 0 0 0 0 1 -360 360;</p> <p> 43 71 0.00492 0.02649 0.0442 9900 0 0 0 0 1 -360 360;</p> <p> 45 43 0.00138 0.01063 0.0345 9900 0 0 0 0 1 -360 360;</p> <p> 47 42 0.00805 0.0433 0.0723 9900 0 0 0 0 1 -360 360;</p> <p> 43 46 0.00138 0.01063 0.0345 9900 0 0 0 0 1 -360 360;</p> <p> 47 46 0.01001 0.05382 0.0898 9900 0 0 0 0 1 -360 360;</p> <p> 43 68 0.00431 0.03309 0.1076 9900 0 0 0 0 1 -360 360;</p> <p> 43 41 6e-005 0.0005 0.0016 9900 0 0 0 0 1 -360 360;</p> <p> 47 68 0.0154 0.0828 0.1383 9900 0 0 0 0 1 -360 360;</p> <p> 47 68 0.005 0.03841 0.1249 9900 0 0 0 0 1 -360 360;</p> <p> 50 69 0.01139 0.06127 0.1023 9900 0 0 0 0 1 -360 360;</p> <p> 56 50 0.02494 0.13413 0.224 9900 0 0 0 0 1 -360 360;</p> <p> 70 59 0.00105 0.00809 0.0263 9900 0 0 0 0 1 -360 360;</p> <p> 70 60 0.00639 0.03726 0.0622 9900 0 0 0 0 1 -360 360;</p> <p> 62 58 0.02156 0.11592 0.1936 9900 0 0 0 0 1 -360 360;</p> <p> 58 56 0.0215 0.13496 0.2254 9900 0 0 0 0 1 -360 360;</p> <p> 58 69 0.03988 0.21445 0.3581 9900 0 0 0 0 1 -360 360;</p> <p> 55 52 0.00562 0.03063 0.0511 9900 0 0 0 0 1 -360 360;</p> <p> 55 54 0.00092 0.00049 0.0082 9900 0 0 0 0 1 -360 360;</p> <p> 52 54 0.00184 0.00993 0.0165 9900 0 0 0 0 1 -360 360;</p> <p> 55 49 0.00431 0.02318 0.0387 9900 0 0 0 0 1 -360 360;</p> <p> 50 55 0.00631 0.03339 0.0567 9900 0 0 0 0 1 -360 360;</p> <p> 69 56 0.00739 0.05673 0.1845 9900 0 0 0 0 1 -360 360;</p> <p>];</p> <p> </p> </font></font><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><p>%%----- OPF Data -----%%</p> <p>%% generator cost data</p> <p>% 1 startup shutdown n x1 y1 ... xn yn</p> <p>% 2 startup shutdown n c(n-1) ... c0</p> </font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2"><p>mpc.gencost = [</p> <p> 2 0 0 3 0.01 40 0;</p> <p> 2 0 0 3 0.0164744646 20 0;</p> <p> 2 0 0 3 0.00653594771 20 0;</p> <p> 2 0 0 3 0.0114942529 20 0;</p> <p> 2 0 0 3 0.0113507378 20 0;</p> <p> 2 0 0 3 0.0480769231 20 0;</p> <p> 2 0 0 3 0.0212314225 20 0;</p> <p> 2 0 0 3 0.0240384615 20 0;</p> <p> 2 0 0 3 0.0571428571 20 0;</p> <p> 2 0 0 3 0.0384615385 20 0;</p> <p> 2 0 0 3 0.029154519 20 0;</p> <p> 2 0 0 3 0.0641025641 20 0;</p> <p> 2 0 0 3 0.012180268 20 0;</p> <p> 2 0 0 3 0.0116959064 20 0;</p> <p> 2 0 0 3 0.0574712644 20 0;</p> <p> 2 0 0 3 0.01 40 0;</p> <p> 2 0 0 3 0.0182149362 20 0;</p> <p> 2 0 0 3 0.0147710487 20 0;</p> <p> 2 0 0 3 0.0294117647 20 0;</p> <p>];</p></font></font><br></div><div style="color: rgb(0, 0, 0);">Raza Farid<br>Electrical Design Engineer<br>Green Building Technologies<br><br>Contact:<br>+92-333-5615025<br>[email protected]</div><div style="color: rgb(0, 0, 0);"><br></div></div></div> </body> </html> -- -- Fred Kanjelesa. NCC-Trading/ZESCO Zambia/Lusaka. Mobile: +260967/955-880912 Extension: 5124 "I can do all things through Christ who strengthens me." Philippians 4:13
