Dr. Ray I'm so sorry as I did not see your reply on Jan 2nd. I recieve a message containing "You are out of office". Sorry again for the inconvenience. Concerning my question, Yes I started with the feasible problem you mentioned. It runs for power flow but it does not run for optimal power flow. Thanks,
Amal Elramly ________________________________ From: Ray Zimmerman <[email protected]> To: MATPOWER discussion forum <[email protected]> Cc: Amal Elramly <[email protected]> Sent: Monday, January 7, 2013 6:45 PM Subject: Re: question As I said in my reply on Jan 2nd … I'm guessing it is infeasible. Did you start with the feasible problem I mentioned in my response from Nov 28th? Please follow through on my suggestions before continuing to post the same question. -- Ray Zimmerman Senior Research Associate 419A Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645 On Jan 6, 2013, at 6:09 AM, Amal Elramly <[email protected]> wrote: > > > >Dear Dr. Ray > >I have tried to run optimal power flow for the 33 radial test bus system but >it does not converge. Knowing that the case file works for power flow (runpf). >I can not determine the source of the problem. Could you help me? > >Attached below the program warnings and results: > >>> runopf('case_amal5') >MATPOWER Version 4.1, 14-Dec-2011 -- AC Optimal Power Flow >MATLAB Interior Point Solver -- MIPS, Version 1.0, 07-Feb-2011Warning: Matrix is close to singular or badly scaled. > Results may be inaccurate. RCOND = 5.247763e-018. >> In mips at 422 > In mipsopf_solver at 145 > In opf_execute at 106 > In opf at 225 > In runopf at 96 >Warning: Matrix is close to singular or badly scaled. > Results may be inaccurate. RCOND = 3.349638e-022. >> In mips at 422 > In mipsopf_solver at 145 > In opf_execute at 106 > In opf at 225 > In runopf at 96 >Numerically Failed >Did not converge in 13 iterations. >Did not converge (0.28 seconds) >Objective Function Value = 4.40 $/hr >================================================================================ >| 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) 2.2 0.5 >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) 0.13 0.09 >Transformers 0 Branch Charging (inj) - 0.0 >Inter-ties 0 Total Inter-tie Flow 0.0 0.0 >Areas 1 > Minimum Maximum > ------------------------- -------------------------------- >Voltage Magnitude 1.506 p.u. @ bus 18 1.532 p.u. @ bus 2 >Voltage Angle -0.16 deg @ bus 18 0.09 deg @ bus 30 >P Losses (I^2*R) - 0.04 MW @ line 1-2 >Q Losses (I^2*X) - 0.02 MVAr @ line 1-2 >Lambda P -6674.13 $/MWh @ bus 1 15413474426196144000.00 $/MWh @ >bus 18 >Lambda Q 4445.12 $/MWh @ bus 1 8466416009517438000.00 $/MWh @ bus >33 >================================================================================ >| Bus Data | >================================================================================ > Bus Voltage Generation Load >Lambda($/MVA-hr) > # Mag(pu) Ang(deg) P (MW) Q (MVAr) P (MW) Q (MVAr) P Q > >----- ------- -------- -------- -------- -------- -------- ------- ------- > 1 1.526 0.000* 2.15 0.48 - - -6674.1314445.122 > 2 1.532 -0.001 - - 0.10 0.06 >482471203417005.810249602520478259.500 > 3 1.523 0.008 - - 0.09 0.04 >15072572746478563000.0008185918922953839600.000 > 4 1.522 0.017 - - 0.12 0.08 15104025286329487000.0008207876480607969300.000 > 5 1.520 0.026 - - 0.06 0.03 >15134463712170088000.0008229215530480105500.000 > 6 1.516 0.015 - - 0.06 0.02 >15206029152644973000.0008282433372989881300.000 > 7 1.515 -0.023 - - 0.20 0.10 15225286235917232000.0008290320354694113300.000 > 8 1.512 -0.060 - - 0.20 0.10 >15296769731624745000.0008314283903063420900.000 > 9 1.511 -0.079 - - 0.06 0.02 >15328004445576665000.0008322680628071578600.000 > 10 1.510 -0.096 - - 0.06 0.02 >15355208175207971000.0008329745091343159300.000 > 11 1.510 -0.097 - - 0.04 0.03 15359180129512268000.0008330590821724995600.000 > 12 1.510 -0.099 - - 0.06 0.04 >15365913344698071000.0008331654368793930800.000 > 13 1.515 -0.122 - - 0.06 0.04 >15286706390876719000.0008252011468596544500.000 > 14 1.508 -0.136 - - 0.12 0.08 15376665231265518000.0008357579409955088400.000 > 15 1.508 -0.142 - - 0.06 0.01 >15386528993294068000.0008360997889585863700.000 > 16 1.507 -0.146 - - 0.06 0.02 >15395689234238870000.0008364574717577316400.000 > 17 1.506 -0.160 - - 0.06 0.02 >15409536653083820000.0008370346040024593400.000 > 18 1.506 -0.162 - - 0.09 0.04 15413474426196144000.0008372131740058532900.000 > 19 1.526 -0.003 - - 0.09 0.04 >483409380339560.500250432571149656.120 > 20 1.525 -0.021 - - 0.09 0.04 >483931362873987.250250466952173264.060 > 21 1.525 -0.025 - - 0.09 0.04 484020789441837.310250453373003567.310 > 22 1.531 -0.029 - - 0.09 0.04 >481787513781547.190247382801145038.620 > 23 1.528 0.002 - - 0.42 0.05 >14979643414630095000.0008118178854098206700.000 > 24 1.521 0.004 - - 0.42 0.20 >15106484506421404000.0008219389687429146600.000 > 25 1.521 0.014 - - 0.06 0.20 15109715780803389000.0008229638628413411300.000 > 26 1.516 0.022 - - 0.06 0.03 >15211131578805922000.0008289737720668774400.000 > 27 1.521 0.031 - - 0.06 0.03 >15098406719421895000.0008238958818964934700.000 > 28 1.512 0.047 - - 0.12 0.02 15249229614526636000.0008384065566448683000.000 > 29 1.510 0.064 - - 0.20 0.07 >15283185053805144000.0008427418723589714900.000 > 30 1.509 0.086 - - 0.15 0.60 >15296412867248689000.0008449006203481891800.000 > 31 1.508 0.077 - - 0.21 0.07 >15317078797228669000.0008462351132649569300.000 > 32 1.507 0.077 - - 0.06 0.10 15319642074435744000.0008465341544258864100.000 > 33 1.507 0.076 - - 0.06 0.04 >15321213791940397000.0008466416009517438000.000 > -------- -------- -------- -------- > Total: 2.15 0.48 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 -11.78 -6.13 11.83 6.16 0.044 0.02 > 2 2 3 3.41 1.85 -3.39 -1.84 0.020 0.01 > 3 3 4 0.78 0.55 -0.78 -0.55 0.001 0.00 > 4 4 5 0.72 0.51 -0.72 -0.51 0.001 0.00 > 5 5 6 0.69 0.51 -0.69 -0.51 0.002 0.00 > 6 6 7 0.44 0.17 -0.44 -0.17 0.000 0.00 > 7 7 8 0.35 0.11 -0.35 -0.11 0.001 0.00 > 8 8 9 0.25 0.06 -0.25 -0.06 0.000 0.00 > 9 9 10 0.22 0.05 -0.22 -0.05 0.000 0.00 > 10 10 11 0.19 0.04 -0.19 -0.04 0.000 0.00 > 11 11 12 0.17 0.03 -0.17 -0.03 0.000 0.00 > 12 12 13 -0.45 -0.46 0.46 0.46 0.002 0.00 > 13 13 14 1.15 1.35 -1.14 -1.34 0.005 0.01 > 14 14 15 0.13 0.05 -0.13 -0.05 0.000 0.00 > 15 15 16 0.10 0.04 -0.10 -0.04 0.000 0.00 > 16 16 17 0.07 0.03 -0.07 -0.03 0.000 0.00 > 17 17 18 0.04 0.02 -0.04 -0.02 0.000 0.00 > 18 2 19 4.67 4.35 -4.65 -4.34 0.018 0.02 > 19 19 20 0.09 0.01 -0.09 -0.01 0.000 0.00 > 20 20 21 0.05 -0.01 -0.05 0.01 0.000 0.00 > 21 21 22 -0.72 -0.99 0.72 0.99 0.003 0.00 > 22 3 23 -1.85 -1.35 1.86 1.36 0.006 0.00 > 23 23 24 1.21 0.97 -1.21 -0.96 0.006 0.00 > 24 24 25 0.03 0.10 -0.03 -0.10 0.000 0.00 > 25 6 26 0.21 0.33 -0.21 -0.33 0.000 0.00 > 26 26 27 -3.76 -1.69 3.77 1.70 0.013 0.01 > 27 27 28 1.15 1.11 -1.15 -1.10 0.007 0.01 > 28 28 29 0.34 0.43 -0.34 -0.43 0.001 0.00 > 29 29 30 0.24 0.39 -0.24 -0.39 0.000 0.00 > 30 30 31 0.16 0.10 -0.16 -0.10 0.000 0.00 > 31 31 32 0.06 0.07 -0.06 -0.07 0.000 0.00 > 32 32 33 0.03 0.02 -0.03 -0.02 0.000 0.00 > -------- -------- > Total: 0.130 0.09 >================================================================================ >| Voltage Constraints | >================================================================================ >Bus # Vmin mu Vmin |V| Vmax Vmax mu >----- -------- ----- ----- ----- -------- > 1 6246.856 0.950 1.526 1.0501272740518625929200.000 > 2 6302.525 0.950 1.532 1.100355570291.153 > 3 6286.023 0.950 1.523 1.050234659548248.324 > 4 6304.242 0.950 1.522 1.05070211654568.559 > 5 6322.540 0.950 1.520 1.05029853950572.320 > 6 6369.737 0.950 1.516 1.0506097293404.070 > 7 6379.032 0.950 1.515 1.0504859694314.642 > 8 6415.832 0.950 1.512 1.0502062405473.305 > 9 6433.196 0.950 1.511 1.0501426960223.656 > 10 6449.101 0.950 1.510 1.050971125114.997 > 11 6451.543 0.950 1.510 1.050907364464.729 > 12 6455.813 0.950 1.510 1.050785716565.193 > 13 6529.122 0.950 1.515 1.10085464084.004 > 14 6480.352 0.950 1.508 1.0502838355.232 > 15 6484.456 0.950 1.508 1.050272204085.475 > 16 6488.425 0.950 1.507 1.050279697729.872 > 17 6494.291 0.950 1.506 1.050327620045.383 > 18 6496.046 0.950 1.506 1.050333944214.087 > 19 6253.843 0.950 1.526 1.0503008201144580546.000 > 20 6262.180 0.950 1.525 1.0505934815489843.966 > 21 6263.675 0.950 1.525 1.0501316284247402.528 > 22 6315.013 0.950 1.531 1.100325795191.182 > 23 6345.784 0.950 1.528 1.100251631455.982 > 24 6307.177 0.950 1.521 1.05048890166499.626 > 25 6310.559 0.950 1.521 1.05038714021061.191 > 26 6374.732 0.950 1.516 1.0502324115755.694 > 27 6434.834 0.950 1.521 1.100139848954.055 > 28 6413.029 0.950 1.512 1.0501912131557.960 > 29 6434.480 0.950 1.510 1.0501156823231.291 > 30 6442.778 0.950 1.509 1.050978580603.251 > 31 6452.337 0.950 1.508 1.050818185087.477 > 32 6453.938 0.950 1.507 1.050795031950.793 > 33 6454.696 0.950 1.507 1.050784387575.451 >================================================================================ >| Generation Constraints | >================================================================================ > Gen Bus Active Power Limits > # # Pmin mu Pmin Pg Pmax Pmax mu >---- ----- ------- -------- -------- -------- ------- > 1 1 10275.446 0.00 2.15 6.00 3600.015 >Gen Bus Reactive Power Limits > # # Qmin mu Qmin Qg Qmax Qmax mu >--- --- ------- -------- -------- -------- ------- > 1 16364.618 -0.50 0.48 1.0010809.675 >================================================================================ >| Branch Flow Constraints | >================================================================================ >Brnch From "From" End Limit "To" End To > # Bus |Sf| mu |Sf| |Smax| |St| |St| mu Bus >----- ----- ------- -------- -------- -------- ------- ----- > 1 1 4.887 13.28 0.10 13.33 4.342 2 > 2 2 0.136 3.88 0.10 3.85444868379062962690.000 3 > 3 38857013.775 0.96 0.20 0.968980790.433 4 > 4 4190975.730 0.89 0.10 0.89196199.302 5 > 5 5 95220.586 0.86 0.10 0.85100930.224 6 > 6 6 28641.260 0.47 0.30 0.47 28714.067 7 > 7 7 19558.624 0.37 0.30 0.37 19668.295 8 > 8 8 5819.117 0.26 0.10 0.26 5825.381 9 > 9 9 5534.336 0.23 0.10 0.23 5538.127 10 > 10 10 5302.956 0.20 0.10 0.20 5303.362 11 > 11 11 5150.910 0.17 0.10 0.17 5151.381 12 > 12 12 2172.596 0.64 0.10 0.64 1984.649 13 > 13 13 2250.949 1.77 0.20 1.76 2649.396 14 > 14 14 5063.615 0.14 0.10 0.14 5064.110 15 > 15 15 4945.511 0.11 0.10 0.11 4945.803 16 > 16 16 4854.824 0.08 0.10 0.08 4855.038 17 > 17 17 4795.020 0.05 0.10 0.05 4795.043 18 > 18 2 22.331 6.38 0.10 6.36 161.524 19 > 19 19 4882.909 0.10 0.10 0.10 4883.192 20 > 20 20 4782.472 0.05 0.10 0.05 4782.467 21 > 21 21 1531.158 1.22 0.10 1.23 1351.125 22 > 22 3 3636.219 2.29 0.50 2.30 3132.169 23 > 23 23 7914.348 1.55 0.50 1.54 9907.350 24 > 24 24 13363.654 0.10 0.30 0.10 13364.250 25 > 25 6 6880.566 0.39 0.10 0.39 6885.089 26 > 26 26 483.983 4.12 0.10 4.14 333.084 27 > 27 27 17981.684 1.60 0.20 1.59176578.511 28 > 28 28 45384.414 0.55 0.30 0.55 45881.876 29 > 29 29 35066.816 0.46 0.70 0.46 35143.551 30 > 30 30 14458.864 0.19 0.30 0.19 14465.501 31 > 31 31 9401.831 0.09 0.20 0.09 9401.992 32 > 32 32 4778.881 0.04 0.10 0.04 4778.887 33 > Thanks for your patience and time >Amal Elramly > > > >________________________________ > From: Amal Elramly <[email protected]> >To: MATPOWER discussion forum <[email protected]> >Sent: Sunday, December 9, 2012 11:06 AM >Subject: Fw: question > > > > > > >Dear Dr. Ray > >Thanks very much for your support. I made use of the case file that you has >sent me and I have changed the values of BR_R and BR_X to per unit values >dividing them by Zbase and then the program converges after 5 iterations and I >think the results are reasonable. I will try to run optimal power flow. > Thanks > >Amal Elramly > > > > > >________________________________ > From: Ray Zimmerman <[email protected]> >To: MATPOWER discussion forum <[email protected]> >Sent: Wednesday, November 28, 2012 6:16 PM >Subject: Re: question > > >Ok, I spent some time trying to get to something that converges … > > >- there is no reference bus (set bus type to 3 for bus 1) >- gencost format is incorrect, it should be polynomial (set first column to 2) >- the bus matrix is missing the VM column > > >If you fix these issues, then it is a valid, but infeasible problem. I believe >the infeasibility is related to voltage and reactive power. In case it is >useful, dividing BR_R by 1000, BR_X by 100 and setting RATE_A to 0 in the >branch matrix and dividing QD by 10 in the bus matrix, yields a feasible >problem. This may give you a useful starting point for this topology from >which to begin investigating where the unfeasibility comes from. > > >-- >Ray Zimmerman >Senior Research Associate >419A Warren Hall, Cornell University, Ithaca, NY 14853 >phone: (607) 255-9645 > > > > >On Nov 28, 2012, at 3:47 AM, Amal Elramly <[email protected]> wrote: > >Dr. Ray >> >>yes all lines have non-zero reactance >>The case file of the test bus system is attached below: >> system MVA basempc.baseMVA = 1; >>%% bus data >>% bus_i type Pd Qd Gs Bs Vm Va baseKV zone Vmax >>Vminmpc.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 >>apfmpc.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 angmaxmpc.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 >>> >>> >>> >> >> >> > > > > > > > > >
