When you get errors like that, always check the case with case_info() to see if
there are any obvious errors in the case file. In this case, it appears that
you are missing bus 14.
And a clarification about power flow … by default generator limits are ignored
when solving the power flow. You can turn on the option to enforce the VAr
limits, in which case any generator exceeding the VAr limits will be converted
to PQ bus at the limit (i.e. voltage will no longer be maintained at the
setpoint).
I suspect that the main issue with using MATPOWER to solve the power flow for a
balanced distribution system is that the R/X ratios are often so high that the
Newton-Raphson solver does not converge from a flat voltage start.
Ray
>> case_info('case33_dis')
Unknown branch "from" buses: 1
branch(14, F_BUS) = 14
Unknown branch "to" buses: 1
branch(13, T_BUS) = 14
Checking connectivity ... 2 connected groups, 0 isolated buses
Elapsed time is 0.038700 seconds.
================================================================================
Full Island Island
System 1 2
Number of: ---------- ---------- ----------
buses 32 28 4
loads 31 27 4
on 31 27 4
off - - -
fixed 31 27 4
dispatchable - - -
on - - -
off - - -
generators 5 5 -
on 5 5 -
off - - -
shunt elements 1 1 -
branches 30 27 3
on 30 27 3
off - - -
ties (off) - - -
Load
active (MW)
dispatched 3595.0 3325.0 270.0
fixed 3595.0 3325.0 270.0
dispatchable - - -
nominal 3595.0 3325.0 270.0
on 3595.0 3325.0 270.0
off - - -
fixed 3595.0 3325.0 270.0
dispatchable - - -
on - - -
off - - -
reactive (MVAr)
dispatched 2220.0 2130.0 90.0
fixed 2220.0 2130.0 90.0
dispatchable - - -
nominal 2220.0 2130.0 90.0
on 2220.0 2130.0 90.0
off - - -
fixed 2220.0 2130.0 90.0
dispatchable - - -
on - - -
off - - -
Generation
active (MW)
dispatched 5500.0 5500.0 -
max capacity - - -
on - - -
off - - -
min capacity - - -
on - - -
off - - -
reactive (MVAr)
dispatched 3500.0 3500.0 -
max capacity - - -
on - - -
off - - -
min capacity - - -
on - - -
off - - -
Shunt Injections
active (MW) - - -
reactive (MVAr) 19.0 19.0 -
Branch Losses
active (MW) - - -
reactive (MVAr) - - -
DC line
export (MW)
dispatch - - -
max capacity - - -
on - - -
off - - -
min capacity - - -
on - - -
off - - -
Reference Buses
num of ref buses 1 1 0
ref bus numbers 1 1
>>
> On May 5, 2016, at 9:38 AM, Fiaz Ahmad <fiazah...@sabanciuniv.edu> wrote:
>
> I am also getting this problem even after doing the changes you suggested....
> I think now Dr. Ray should help us now
>
>
> Have a nice day
>
> Regards,
> FIAZ AHMAD
> PhD Student,
> FENS Sabanci University,
> Istanbul Turkey
>
>
>
> On Thu, May 5, 2016 at 4:31 PM, Jovan Ilic <jovan.i...@gmail.com
> <mailto:jovan.i...@gmail.com>> wrote:
>
> Fiaz,
>
> I get an error when running runpf function on your data using Matlab 2015b
> and MatPower 5.1:
>
> >> runpf(case33_dis)
> Subscript indices must either be real positive integers or logicals.
>
> Error in ext2int (line 149)
> brs = ( mpc.branch(:, BR_STATUS) & ... %% branch status
>
> Error in runpf (line 111)
> mpc = ext2int(mpc);
>
> >>
>
> Test cases that came with MatPower 5.1 run without a problem. The mpc
> structure seems to be read correctly. Maybe Ray (the developer) can
> help with this.
>
> I would suggest deleting extra generators (2-8) and corresponding cost
> functions. I would also use mBase = 1, the same as system mBase at the
> beginning. Make sure that matpower directory is in the path and just run
> runpf(case33_dis) in the directory that contains case33_dis.m file.
>
> Jovan
>
>
> On Thu, May 5, 2016 at 8:14 AM, Fiaz Ahmad <fiazah...@sabanciuniv.edu
> <mailto:fiazah...@sabanciuniv.edu>> wrote:
> %% MATPOWER Case Format : Version 2
> mpc.version = '2';
>
> %%----- Power Flow Data -----%%
> %% system MVA base
> mpc.baseMVA = 1;
>
> %% 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 0
> 11 1 1.02 0.98;
> 2 1 100 60 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 3 1 90 40 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 4 1 120 80 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 5 1 60 30 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 6 1 60 20 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 7 1 200 100 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 8 1 200 100 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 9 1 60 20 0 19 1 1.0 0
> 0 1 1.02 0.98;
> 10 1 60 20 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 11 1 45 30 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 12 1 60 35 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 13 1 60 35 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 15 1 60 10 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 16 1 60 20 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 17 1 60 20 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 18 1 90 40 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 19 1 90 40 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 20 1 90 40 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 21 1 90 40 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 22 1 90 40 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 23 1 90 50 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 24 1 420 200 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 25 1 420 200 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 26 1 60 25 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 27 1 60 25 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 28 1 60 20 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 29 1 120 70 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 30 1 200 600 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 31 1 150 70 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 32 1 210 100 0 0 1 1.0 0
> 0 1 1.02 0.98;
> 33 1 60 40 0 0 1 1.0 0
> 0 1 1.02 0.98;];
>
> %% 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 5500 3500 0 0 1.0 100 1 0
> 0 0 0 0 0 0 0 0 0
> 0 0 0;
> 2 0 0 0 0 1.0 100 1 0 0
> 0 0 0 0 0 0 0 0 0
> 0 0;
> 3 0 0 0 0 1.0 100 1 0 0
> 0 0 0 0 0 0 0 0 0
> 0 0;
> 6 0 0 0 0 1.0 100 1 0 0
> 0 0 0 0 0 0 0 0 0
> 0 0;
> 8 0 0 0 0 1.0 100 1 0 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.0470 0 0 0 0
> 0 0 1 -360 360;
> 2 3 0.4930 0.2511 0 0 0 0 0
> 0 1 -360 360;
> 3 4 0.3660 0.1864 0 0 0 0
> 0 0 1 -360 360;
> 4 5 0.3811 0.1941 0 0 0 0
> 0 0 1 -360 360;
> 5 6 0.8190 0.7070 0 0 0 0 0
> 0 1 -360 360;
> 6 7 0.1872 0.6188 0 0 0 0
> 0 0 1 -360 360;
> 7 8 0.7114 0.2351 0 0 0 0
> 0 0 1 -360 360;
> 8 9 1.0300 0.7400 0 0 0 0 0
> 0 1 -360 360;
> 9 10 1.0440 0.7400 0 0 0 0 0
> 0 1 -360 360;
> 10 11 0.1966 0.0650 0 0 0 0 0
> 0 1 -360 360;
> 11 12 0.3744 0.1238 0 0 0 0
> 0 0 1 -360 360;
> 12 13 1.4680 1.1550 0 0 0 0
> 0 0 1 -360 360;
> 13 14 0.5416 0.7129 0 0 0 0 0
> 0 1 -360 360;
> 14 15 0.5910 0.5260 0 0 0 0
> 0 0 1 -360 360;
> 15 16 0.7463 0.5450 0 0 0 0
> 0 0 1 -360 360;
> 16 17 1.2890 1.7210 0 0 0 0
> 0 0 1 -360 360;
> 17 18 0.7320 0.5740 0 0 0 0 0
> 0 1 -360 360;
> 2 19 0.1640 0.1565 0 0 0 0 0
> 0 1 -360 360;
> 19 20 1.5042 1.3554 0 0 0 0 0
> 0 1 -360 360;
> 20 21 0.4095 0.4784 0 0 0 0 0
> 0 1 -360 360;
> 21 22 0.7089 0.9373 0 0 0 0 0 0
> 1 -360 360;
> 3 23 0.4512 0.3083 0 0 0 0 0 0
> 1 -360 360;
> 23 24 0.8980 0.7091 0 0 0 0 0 0
> 1 -360 360;
> 24 25 0.8960 0.7011 0 0 0 0 0 0
> 1 -360 360;
> 6 26 0.2030 0.1034 0 0 0 0 0 0
> 1 -360 360;
> 26 27 0.2842 0.1447 0 0 0 0 0 0
> 1 -360 360;
> 27 28 1.0590 0.9337 0 0 0 0 0 0
> 1 -360 360;
> 28 29 0.8042 0.7006 0 0 0 0 0 0
> 1 -360 360;
> 29 30 0.5075 0.2585 0 0 0 0 0 0
> 1 -360 360;
> 30 31 0.9744 0.9630 0 0 0 0 0 0
> 1 -360 360;
> 31 32 0.3105 0.3619 0 0 0 0 0 0
> 1 -360 360;
> 32 33 0.3410 0.5302 0 0 0 0 0 0
> 1 -360 360;
> ];
>
>
> Have a nice day
>
> Regards,
> FIAZ AHMAD
> PhD Student,
> FENS Sabanci University,
> Istanbul Turkey
>
>
>
> On Thu, May 5, 2016 at 3:11 PM, Jovan Ilic <jovan.i...@gmail.com
> <mailto:jovan.i...@gmail.com>> wrote:
>
> Where did you add P=5000W and Q = 3500VAR in your data, it is not in the data
> you
> posted? If I remember well the default MVA base in MatPower is 100MVA; your
> P and Q data seems to be in W and VAR . Also, if I remember well, if you
> specify 0
> for Pmax and Pmin and Qmax and Qmin there will be no limits on P and Q.
>
> It seems that you need to read the MatPower Users Manual and look at the case
> files
> that come with the package. The User's Manual has very good description of
> the data
> format used.
>
> Jovan
>
>
> On Thu, May 5, 2016 at 7:54 AM, Fiaz Ahmad <fiazah...@sabanciuniv.edu
> <mailto:fiazah...@sabanciuniv.edu>> wrote:
> I am putting P=5000W and Q=3500VARs as a substation. How can I make P Q
> unlimited? Any suggestion
>
> On May 5, 2016 2:48 PM, "Jovan Ilic" <jovan.i...@gmail.com
> <mailto:jovan.i...@gmail.com>> wrote:
>
> Fiaz,
>
> If you are using the data in your original post only, there are no generators
> in
> that data set. You should put a (unlimited P and Q) generator where the
> substation
> would be in your system.
>
> Jovan
>
>
> On Thu, May 5, 2016 at 7:29 AM, Fiaz Ahmad <fiazah...@sabanciuniv.edu
> <mailto:fiazah...@sabanciuniv.edu>> wrote:
> Yes I m using balanced three phase but how can I do it with matpower. I
> changed the bus data and gen.data matrices but I couldn't get the power flows.
>
> On May 5, 2016 2:24 PM, "Jovan Ilic" <jovan.i...@gmail.com
> <mailto:jovan.i...@gmail.com>> wrote:
>
> You could use Matpower to simulate distribution systems if you are only
> interested in balanced 3 phase without mutual line coupling systems.
>
> Jovan
>
>
> On Thu, May 5, 2016 at 3:41 AM, Ji Chen <hbxia...@163.com
> <mailto:hbxia...@163.com>> wrote:
> Hi Shri,
> I am interested in power flows for distribution networks recently. Do you
> mean that Matpower is not a good chioce for distribution networks simulation
> or it just can not do distribution networks simulation?
>
> Thanks,
> Ji
>
>
>
> At 2016-05-04 23:02:50, "Abhyankar, Shrirang G." <abhy...@anl.gov
> <mailto:abhy...@anl.gov>> wrote:
> http://ewh.ieee.org/soc/pes/dsacom/testfeeders/index.html
> <http://ewh.ieee.org/soc/pes/dsacom/testfeeders/index.html>
>
> However, as Jovan mentioned, MATPOWER is for solving balanced three-phase
> transmission systems. If you are interested in solving unbalanced three-phase
> power flows for distribution networks then refer to either OpenDSS or
> GridLAB-D.
>
> Shri
>
> From: Fiaz Ahmad <fiazah...@sabanciuniv.edu
> <mailto:fiazah...@sabanciuniv.edu>>
> Reply-To: MATPOWER discussion forum <matpowe...@list.cornell.edu
> <mailto:matpowe...@list.cornell.edu>>
> Date: Wednesday, May 4, 2016 at 8:49 AM
> To: MATPOWER discussion forum <matpowe...@list.cornell.edu
> <mailto:matpowe...@list.cornell.edu>>, "jovan.i...@gmail.com
> <mailto:jovan.i...@gmail.com>" <jovan.i...@gmail.com
> <mailto:jovan.i...@gmail.com>>
> Subject: Re: modification on Matpower case.m file for power flow of
> distribution system
>
> Thanks Jovan,
>
> In Matpower there are no test cases for distribution systems... Where can i
> find IEEE distribution system test cases?
>
> Regards
>
>
> Have a nice day
>
> Regards,
> FIAZ AHMAD
> PhD Student,
> FENS Sabanci University,
> Istanbul Turkey
>
>
>
> On Wed, May 4, 2016 at 4:46 PM, Jovan Ilic <jovan.i...@gmail.com
> <mailto:jovan.i...@gmail.com>> wrote:
>
> Fiaz,
>
> You are missing lots of data here. I suggest looking at IEEE distribution
> system
> test cases to see what is involved. Matpower is used for transmission
> systems
> not distribution systems.
>
> Jovan
>
> On Wed, May 4, 2016 at 3:45 AM, Fiaz Ahmad <fiazah...@sabanciuniv.edu
> <mailto:fiazah...@sabanciuniv.edu>> wrote:
> Hi. i have the following data files for ieee 33 bus radial distribution
> system. i want to modify the case33.m file in matpower to carry out the power
> flow of the above distribution system. Please guide me how can i achieve that.
>
> linedata33
> % Bus fbus tbus R X
> 1 1 2 0.0922 0.0470
> 2 2 3 0.4930 0.2511
> 3 3 4 0.3660 0.1864
> 4 4 5 0.3811 0.1941
> 4 5 6 0.8190 0.7070
> 6 6 7 0.1872 0.6188
> 7 7 8 0.7114 0.2351
> 8 8 9 1.0300 0.7400
> 9 9 10 1.0440 0.7400
> 10 10 11 0.1966 0.0650
> 11 11 12 0.3744 0.1238
> 12 12 13 1.4680 1.1550
> 13 13 14 0.5416 0.7129
> 14 14 15 0.5910 0.5260
> 15 15 16 0.7463 0.5450
> 16 16 17 1.2890 1.7210
> 17 17 18 0.7320 0.5740
> 18 2 19 0.1640 0.1565
> 19 19 20 1.5042 1.3554
> 20 20 21 0.4095 0.4784
> 21 21 22 0.7089 0.9373
> 22 3 23 0.4512 0.3083
> 23 23 24 0.8980 0.7091
> 24 24 25 0.8960 0.7011
> 25 6 26 0.2030 0.1034
> 26 26 27 0.2842 0.1447
> 27 27 28 1.0590 0.9337
> 28 28 29 0.8042 0.7006
> 29 29 30 0.5075 0.2585
> 30 30 31 0.9744 0.9630
> 31 31 32 0.3105 0.3619
> 32 32 33 0.3410 0.5302
>
> loaddata33
> %Bus P Q
> 1 0 0 0
> 2 100 60 0
> 3 90 40 0
> 4 120 80 0
> 5 60 30 0
> 6 60 20 0
> 7 200 100 0
> 8 200 100 0
> 9 60 20 0
> 10 60 20 0
> 11 45 30 0
> 12 60 35 0
> 13 60 35 0
> 14 120 80 0
> 15 60 10 0
> 16 60 20 0
> 17 60 20 0
> 18 90 40 0
> 19 90 40 0
> 20 90 40 0
> 21 90 40 0
> 22 90 40 0
> 23 90 50 0
> 24 420 200 0
> 25 420 200 0
> 26 60 25 0
> 27 60 25 0
> 28 60 20 0
> 29 120 70 0
> 30 200 600 0
> 31 150 70 0
> 32 210 100 0
> 33 60 40 0
>
> Please help and guide.... Thanks
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
