Dr. Ray yes all lines have non-zero reactance The case file of the test bus system is attached below: system MVA base mpc.baseMVA = 1; %% bus data% bus_i type Pd Qd Gs Bs Vm Va baseKV zone Vmax Vmin mpc.bus = [1 2 0 0 0 0 1 0 12.6 1 1.05 0.95; 2 1 0.1 0.06 0 0 1 0 12.6 1 1.1 0.95; 3 1 0.09 0.04 0 0 1 0 12.6 1 1.05 0.95; 4 1 0.12 0.08 0 0 1 0 12.6 1 1.05 0.95; 5 1 0.06 0.03 0 0 1 0 12.6 1 1.05 0.95; 6 1 0.06 0.02 0 0 1 0 12.6 1 1.05 0.95; 7 1 0.2 0.1 0 0 1 0 12.6 1 1.05 0.95; 8 1 0.2 0.1 0 0 1 0 12.6 1 1.05 0.95; 9 1 0.06 0.02 0 0 1 0 12.6 1 1.05 0.95; 10 1 0.06 0.02 0 0 1 0 12.6 1 1.05 0.95; 11 1 0.045 0.03 0 0 1 0 12.6 1 1.05 0.95; 12 1 0.06 0.035 0 0 1 0 12.6 1 1.05 0.95; 13 1 0.06 0.035 0 0 1 0 12.6 1 1.1 0.95; 14 1 0.12 0.08 0 0 1 0 12.6 1 1.05 0.95; 15 1 0.06 0.01 0 0 1 0 12.6 1 1.05 0.95; 16 1 0.06 0.02 0 0 1 0 12.6 1 1.05 0.95; 17 1 0.06 0.02 0 0 1 0 12.6 1 1.05 0.95; 18 1 0.09 0.04 0 0 1 0 12.6 1 1.05 0.95; 19 1 0.09 0.04 0 0 1 0 12.6 1 1.05 0.95; 20 1 0.09 0.04 0 0 1 0 12.6 1 1.05 0.95; 21 1 0.09 0.04 0 0 1 0 12.6 1 1.05 0.95; 22 1 0.09 0.04 0 0 1 0 12.6 1 1.1 0.95; 23 1 0.42 0.05 0 0 1 0 12.6 1 1.1 0.95; 24 1 0.42 0.2 0 0 1 0 12.6 1 1.05 0.95; 25 1 0.06 0.2 0 0 1 0 12.6 1 1.05 0.95; 26 1 0.06 0.025 0 0 1 0 12.6 1 1.05 0.95; 27 1 0.06 0.025 0 0 1 0 12.6 1 1.1 0.95; 28 1 0.12 0.02 0 0 1 0 12.6 1 1.05 0.95; 29 1 0.2 0.07 0 0 1 0 12.6 1 1.05 0.95; 30 1 0.15 0.6 0 0 1 0 12.6 1 1.05 0.95; 31 1 0.21 0.07 0 0 1 0 12.6 1 1.05 0.95; 32 1 0.06 0.1 0 0 1 0 12.6 1 1.05 0.95; 33 1 0.06 0.04 0 0 1 0 12.6 1 1.05 0.95; ];%% 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.5 0 1 -0.5 1 100 1 6 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 2 0.0922 0.0477 0 0.2 0.2 0.2 0 0 1 -360 360; 2 3 0.493 0.2511 0 0.1 0.1 0.1 0 0 1 -360 360; 3 4 0.3660 0.1864 0 0.2 0.2 0.2 0 0 1 -360 360; 4 5 0.3811 0.1941 0 0.1 0.1 0.1 0 0 1 -360 360; 5 6 0.8190 0.7070 0 0.1 0.1 0.1 0 0 1 -360 360; 6 7 0.1872 0.6188 0 0.3 0.3 0.3 0 0 1 -360 360; 7 8 1.7114 1.2351 0 0.3 0.3 0.3 0 0 1 -360 360; 8 9 1.0300 0.7400 0 0.1 0.1 0.1 0 0 1 -360 360; 9 10 1.0040 0.7400 0 0.1 0.1 0.1 0 0 1 -360 360; 10 11 0.1966 0.0650 0 0.1 0.1 0.1 0 0 1 -360 360; 11 12 0.3744 0.1238 0 0.1 0.1 0.1 0 0 1 -360 360; 12 13 1.4680 1.1550 0 0.1 0.1 0.1 0 0 1 -360 360; 13 14 0.5416 0.7129 0 0.2 0.2 0.2 0 0 1 -360 360; 14 15 0.5910 0.5260 0 0.1 0.1 0.1 0 0 1 -360 360; 15 16 0.7463 0.5450 0 0.1 0.1 0.1 0 0 1 -360 360; 16 17 1.2890 1.7210 0 0.1 0.1 0.1 0 0 1 -360 360; 17 18 0.7320 0.5740 0 0.1 0.1 0.1 0 0 1 -360 360; 2 19 0.1640 0.1565 0 0.1 0.1 0.1 0 0 1 -360 360; 19 20 1.5042 1.3554 0 0.1 0.1 0.1 0 0 1 -360 360; 20 21 0.4095 0.4784 0 0.1 0.1 0.1 0 0 1 -360 360; 21 22 0.7089 0.9373 0 0.1 0.1 0.1 0 0 1 -360 360; 3 23 0.4512 0.3083 0 0.5 0.5 0.5 0 0 1 -360 360; 23 24 0.8980 0.7091 0 0.5 0.5 0.5 0 0 1 -360 360; 24 25 0.8960 0.7011 0 0.3 0.3 0.3 0 0 1 -360 360; 6 26 0.2030 0.1034 0 0.1 0.1 0.1 0 0 1 -360 360; 26 27 0.2842 0.1447 0 0.1 0.1 0.1 0 0 1 -360 360; 27 28 1.0590 0.9337 0 0.2 0.2 0.2 0 0 1 -360 360; 28 29 0.8042 0.7006 0 0.3 0.3 0.3 0 0 1 -360 360; 29 30 0.5075 0.2585 0 0.7 0.7 0.7 0 0 1 -360 360; 30 31 0.9744 0.9630 0 0.3 0.3 0.3 0 0 1 -360 360; 31 32 0.3105 0.3619 0 0.2 0.2 0.2 0 0 1 -360 360; 32 33 0.3410 0.5302 0 0.1 0.1 0.1 0 0 1 -360 360; ]; %%generator cost datampc.gencost = [ 1 0 0 3 0.02 2 0; ];
Amal Elramly ________________________________ From: Ray Zimmerman <[email protected]> To: MATPOWER discussion forum <[email protected]> Sent: Tuesday, November 27, 2012 4:52 PM Subject: Re: question Do all lines have non-zero reactance? Do you have a reference bus? -- Ray Zimmerman Senior Research Associate 419A Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645 On Nov 27, 2012, at 3:20 AM, Amal Elramly <[email protected]> wrote: > >I have run power flow for 33 radial test bus system, but the it doesn't >converge and I do not know what is the problem >could you help me ? >Attached below the results of the program: > > >MATPOWER Version 4.1, 14-Dec-2011 -- AC Power Flow (Newton) >Warning: Matrix is singular to working precision. >> In newtonpf at 109 > In runpf at 220 >Warning: Matrix is singular to working precision. >> In newtonpf at 109 > In runpf at 220 >Warning: Matrix is singular to working precision. >> In newtonpf at 109 > In runpf at 220 >Warning: Matrix is singular to working precision. >> In newtonpf at 109 > In runpf at 220 >Warning: Matrix is singular to working precision. >> In newtonpf at 109 > In runpf at 220 >Warning: Matrix is singular to working precision. >> In newtonpf at 109 > In runpf at 220 >Warning: Matrix is singular to working precision. > >Newton's method power did not converge in 10 iterations. >Did not converge (0.09 seconds) >================================================================================ >| System Summary | >================================================================================ >How many? How much? P (MW) Q (MVAr) >--------------------- ------------------- ------------- ----------------- >Buses 33 Total Gen Capacity 6.0 -0.5 to 1.0 >Generators 1 On-line Capacity 6.0 -0.5 to 1.0 >Committed Gens 1 Generation (actual) NaN NaN >Loads 32 Load 3.7 2.3 > Fixed 32 Fixed 3.7 2.3 > Dispatchable 0 Dispatchable -0.0 of -0.0 -0.0 >Shunts 0 Shunt (inj) -0.0 0.0 >Branches 32 Losses (I^2 * Z) NaN NaN >Transformers 0 Branch Charging (inj) - NaN >Inter-ties 0 Total Inter-tie Flow 0.0 0.0 >Areas 1 > Minimum Maximum > ------------------------- -------------------------------- >Voltage Magnitude NaN p.u. @ bus 1 NaN p.u. @ bus 1 >Voltage Angle NaN deg @ bus 1 NaN deg @ bus 1 >P Losses (I^2*R) - NaN MW @ line 1-2 >Q Losses (I^2*X) - NaN MVAr @ line 1-2 >================================================================================ >| Bus Data | >================================================================================ > Bus Voltage Generation Load > # Mag(pu) Ang(deg) P (MW) Q (MVAr) P (MW) Q (MVAr) >----- ------- -------- -------- -------- -------- -------- > 1 NaN NaN NaN NaN - - > 2 NaN NaN - - 0.10 0.06 > 3 NaN NaN - - 0.09 0.04 > 4 NaN NaN - - 0.12 0.08 > 5 NaN NaN - - 0.06 0.03 > 6 NaN NaN - - 0.06 0.02 > 7 NaN NaN - - 0.20 0.10 > 8 NaN NaN - - 0.20 0.10 > 9 NaN NaN - - 0.06 0.02 > 10 NaN NaN - - 0.06 0.02 > 11 NaN NaN - - 0.04 0.03 > 12 NaN NaN - - 0.06 0.04 > 13 NaN NaN - - 0.06 0.04 > 14 NaN NaN - - 0.12 0.08 > 15 NaN NaN - - 0.06 0.01 > 16 NaN NaN - - 0.06 0.02 > 17 NaN NaN - - 0.06 0.02 > 18 NaN NaN - - 0.09 0.04 > 19 NaN NaN - - 0.09 0.04 > 20 NaN NaN - - 0.09 0.04 > 21 NaN NaN - - 0.09 0.04 > 22 NaN NaN - - 0.09 0.04 > 23 NaN NaN - - 0.42 0.05 > 24 NaN NaN - - 0.42 0.20 > 25 NaN NaN - - 0.06 0.20 > 26 NaN NaN - - 0.06 0.03 > 27 NaN NaN - - 0.06 0.03 > 28 NaN NaN - - 0.12 0.02 > 29 NaN NaN - - 0.20 0.07 > 30 NaN NaN - - 0.15 0.60 > 31 NaN NaN - - 0.21 0.07 > 32 NaN NaN - - 0.06 0.10 > 33 NaN NaN - - 0.06 0.04 > -------- -------- -------- -------- > Total: NaN NaN 3.69 2.30 >================================================================================ >| Branch Data >| >================================================================================ >Brnch From To From Bus Injection To Bus Injection Loss (I^2 * Z) > # Bus Bus P (MW) Q (MVAr) P (MW) Q (MVAr) P (MW) Q (MVAr) >----- ----- ----- -------- -------- -------- -------- -------- -------- > 1 1 2 NaN NaN NaN NaN NaN NaN > 2 2 3 NaN NaN NaN NaN NaN NaN > 3 3 4 NaN NaN NaN NaN NaN NaN > 4 4 5 NaN NaN NaN NaN NaN NaN > 5 5 6 NaN NaN NaN NaN NaN NaN > 6 6 7 NaN NaN NaN NaN NaN NaN > 7 7 8 NaN NaN NaN NaN NaN NaN > 8 8 9 NaN NaN NaN NaN NaN NaN > 9 9 10 NaN NaN NaN NaN NaN NaN > 10 10 11 NaN NaN NaN NaN NaN NaN > 11 11 12 NaN NaN NaN NaN NaN NaN > 12 12 13 NaN NaN NaN NaN NaN NaN > 13 13 14 NaN NaN NaN NaN NaN NaN > 14 14 15 NaN NaN NaN NaN NaN NaN > 15 15 16 NaN NaN NaN NaN NaN NaN > 16 16 17 NaN NaN NaN NaN NaN NaN > 17 17 18 NaN NaN NaN NaN NaN NaN > 18 2 19 NaN NaN NaN NaN NaN NaN > 19 19 20 NaN NaN NaN NaN NaN NaN > 20 20 21 NaN NaN NaN NaN NaN NaN > 21 21 22 NaN NaN NaN NaN NaN NaN > 22 3 23 NaN NaN NaN NaN NaN NaN > 23 23 24 NaN NaN NaN NaN NaN NaN > 24 24 25 NaN NaN NaN NaN NaN NaN > 25 6 26 NaN NaN NaN NaN NaN NaN > 26 26 27 NaN NaN NaN NaN NaN NaN > 27 27 28 NaN NaN NaN NaN NaN NaN > 28 28 29 NaN NaN NaN NaN NaN NaN > 29 29 30 NaN NaN NaN NaN NaN NaN > 30 30 31 NaN NaN NaN NaN NaN NaN > 31 31 32 NaN NaN NaN NaN NaN NaN > 32 32 33 NaN NaN NaN NaN NaN NaN > -------- -------- > Total: NaN NaN > > >Amal Elramly > > >
