Hi Sir
Are you planning to include 3 phase load flow for unbalanced distribution
systems in matpower?
Thanking you
On Fri, May 6, 2016 at 1:03 AM, Ray Zimmerman <r...@cornell.edu> wrote:
> 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> 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>
>> 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> 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>
>>>> 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> 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
>>>>>> > 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> 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> 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>
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>> 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>
>>>>>>>>> Reply-To: MATPOWER discussion forum <matpowe...@list.cornell.edu>
>>>>>>>>> Date: Wednesday, May 4, 2016 at 8:49 AM
>>>>>>>>> To: MATPOWER discussion forum <matpowe...@list.cornell.edu>, "
>>>>>>>>> jovan.i...@gmail.com" <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>
>>>>>>>>> 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> 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
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>
>>>>
>>>
>>
>
>
--
Thanking you,
Papiya Dutta (Gupta)
