the problem is still as it was. my 84 bus system is still not converging
after applying all steps directed by the Dr. Ray Zimmerman . can anyone
help me to make the system converge. i am attaching my file
On Fri, Jun 24, 2016 at 6:24 PM, Ray Zimmerman <[email protected]> wrote:
> See FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>.
>
> Ray
>
>
> On Jun 24, 2016, at 3:42 AM, anshu dahiya <[email protected]> wrote:
>
> my 84 bus system is not converging into 10 iteration. what should i take
> steps to make it converge.
>
> On Thu, Jun 23, 2016 at 8:43 PM, Ray Zimmerman <[email protected]> wrote:
>
>> Whenever you have a problem with a case file, the first thing to do is
>> always run case_info()
>> <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0b1/case_info.html>
>> to
>> check for errors in the connectivity information …
>>
>> >> case_info('casem84')
>> Bad bus numbers: 1
>> bus(1, BUS_I) = 0
>>
>> Turns out that bus numbers must be positive integers, as stated in
>> caseformat
>> <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0b1/caseformat.html>
>> .
>>
>> Ray
>>
>>
>>
>> On Jun 23, 2016, at 1:59 AM, anshu dahiya <[email protected]>
>> wrote:
>>
>> sir how to solve the following error in matpower which comes in the
>> case84 which i have attached
>> Error using sparse
>> Index into matrix must be positive.
>>
>> Error in ext2int (line 141)
>> n2i = sparse(mpc.bus(:, BUS_I), ones(nb, 1), 1:nb,
>> max(mpc.bus(:, BUS_I)), 1);
>>
>> Error in runpf (line 111)
>> mpc = ext2int(mpc);
>>
>> On Wed, Jun 22, 2016 at 6:59 PM, Ray Zimmerman <[email protected]> wrote:
>>
>>> Assuming this is for a power flow problem, you still have to decide
>>> whether the bus is voltage controlled or not. If so, make it a PV bus. If
>>> not, and the DG simply injects a fixed amount of power without controlling
>>> voltage, then make it a PQ bus.
>>>
>>> Ray
>>>
>>>
>>>
>>> > On Jun 22, 2016, at 1:27 AM, anshu dahiya <[email protected]>
>>> wrote:
>>> >
>>> > how to nominate DG bus which has load as well as generation
>>> > the nomenclature in the matpower as below
>>> > 1= PQ bus
>>> > 2= PV bus
>>> > 3= reference bus
>>> > 4= isolated bus
>>> > i have attached my file plz take a look of it. suggest me solution.
>>> > <casem84.m>
>>>
>>>
>>>
>>>
>> <casem84.m>
>>
>>
>>
>
>
function mpc = casedg84
mpc.version = '2';
mpc.baseMVA = 100;
% bus_i type Pd Qd Gs Bs area Vm Va
baseKV zone Vmax Vmin
mpc.bus = [
84 3 0 0 0 0 1 1.06 0 11 1 1.082 1;
1 1 0 0 0 0 1 0 0 11 1 1.082 1;
2 1 0.400 0.250 0 0 1 00 0 11 1 1.082 1;
3 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
4 1 0.350 .250 0 0 1 0 0 11 1 1.082 1;
5 1 0.220 .100 0 0 1 0 0 11 1 1.082 1;
6 1 1.100 .800 0 0 1 0 0 11 1 1.082 1;
7 1 .700 .580 0 0 1 0 0 11 1 1.082 1;
8 1 .300 .200 0 0 1 0 0 11 1 1.082 1;
9 1 .300 2.030 0 0 1 0 0 11 1 1.082 1;
10 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
11 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
12 1 1.200 .800 0 0 1 0 0 11 1 1.082 1;
13 1 .800 .600 0 0 1 0 0 11 1 1.082 1;
14 1 .700 .500 0 0 1 0 0 11 1 1.082 1;
15 1 0 0 0 0 1 0 0 11 1 1.082 1;
16 1 .300 .150 0 0 1 0 0 11 1 1.082 1;
17 1 .500 .350 0 0 1 0 0 11 1 1.082 1;
18 1 .700 .400 0 0 1 0 0 11 1 1.082 1;
19 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
20 1 1.500 1.300 0 0 1 0 0 11 1 1.082 1;
21 1 .400 .350 0 0 1 0 0 11 1 1.082 1;
22 1 .050 .020 0 0 1 0 0 11 1 1.082 1;
23 1 .100 .030 0 0 1 0 0 11 1 1.082 1;
24 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
25 1 .050 .030 0 0 1 0 0 11 1 1.082 1;
26 1 .100 .060 0 0 1 0 0 11 1 1.082 1;
27 1 .100 .070 0 0 1 0 0 11 1 1.082 1;
28 1 2.000 1.420 0 0 1 0 0 11 1 1.082 1;
29 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
30 1 0 0 0 0 1 0 0 11 1 1.082 1;
31 1 1.800 1.600 0 0 1 0 0 11 1 1.082 1;
32 1 .200 .100 0 0 1 0 0 11 1 1.082 1;
33 1 1.000 .700 0 0 1 0 0 11 1 1.082 1;
34 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
35 1 .100 .060 0 0 1 0 0 11 1 1.082 1;
36 1 .100 .060 0 0 1 0 0 11 1 1.082 1;
37 1 .020 .010 0 0 1 0 0 11 1 1.082 1;
38 1 .020 .010 0 0 1 0 0 11 1 1.082 1;
39 1 .040 .020 0 0 1 0 0 11 1 1.082 1;
40 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
41 1 .200 .160 0 0 1 0 0 11 1 1.082 1;
42 1 .050 .030 0 0 1 0 0 11 1 1.082 1;
43 1 0 0 0 0 1 0 0 11 1 1.082 1;
44 1 .030 .020 0 0 1 0 0 11 1 1.082 1;
45 1 .800 .700 0 0 1 0 0 11 1 1.082 1;
46 1 .200 .150 0 0 1 0 0 11 1 1.082 1;
47 1 0 0 0 0 1 0 0 11 1 1.082 1;
48 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
49 1 0 0 0 0 1 0 0 11 1 1.082 1;
50 1 .200 .160 0 0 1 0 0 11 1 1.082 1;
51 1 .800 .600 0 0 1 0 0 11 1 1.082 1;
52 1 .500 .300 0 0 1 0 0 11 1 1.082 1;
53 1 .500 .350 0 0 1 0 0 11 1 1.082 1;
54 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
55 1 .700 .380 0 0 1 0 0 11 1 1.082 1;
56 1 0 0 0 0 1 0 0 11 1 1.082 1;
57 1 .030 .020 0 0 1 0 0 11 1 1.082 1;
58 1 .600 .420 0 0 1 0 0 11 1 1.082 1;
59 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
60 1 .020 .010 0 0 1 0 0 11 1 1.082 1;
61 1 .020 .010 0 0 1 0 0 11 1 1.082 1;
62 1 .200 .130 0 0 1 0 0 11 1 1.082 1;
63 1 .600 .440 0 0 1 0 0 11 1 1.082 1;
64 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
65 1 0 0 0 0 1 0 0 11 1 1.082 1;
66 1 .050 .030 0 0 1 0 0 11 1 1.082 1;
67 1 0 0 0 0 1 0 0 11 1 1.082 1;
68 1 .400 .360 0 0 1 0 0 11 1 1.082 1;
69 1 0 0 0 0 1 0 0 11 1 1.082 1;
70 1 2.000 1.500 0 0 1 0 0 11 1 1.082 1;
71 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
72 1 .200 .150 0 0 1 0 0 11 1 1.082 1;
73 1 0 0 0 0 1 0 0 11 1 1.082 1;
74 1 1.200 .950 0 0 1 0 0 11 1 1.082 1;
75 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
76 1 .300 .180 0 0 1 0 0 11 1 1.082 1;
77 1 0 0 0 0 1 0 0 11 1 1.082 1;
78 1 .400 .360 0 0 1 0 0 11 1 1.082 1;
79 1 2.000 1.300 0 0 1 0 0 11 1 1.082 1;
80 1 .200 .140 0 0 1 0 0 11 1 1.082 1;
81 1 .500 .360 0 0 1 0 0 11 1 1.082 1;
82 1 .500 .390 0 0 1 0 0 11 1 1.082 1;
83 2 0 0 0 0 1 1.037 0 11 1 1.082 1;
];
% fbus tbus r x b rateA rateB rateC ratio
angle status angmin angmax
mpc.branch = [
84 1 0.001944 0.006624 0 0 0 0 1 0 1 -360 360;
1 2 0.002096 0.004304 0 0 0 0 1 0 1 -360 360;
2 3 0.002358 0.004842 0 0 0 0 1 0 1 -360 360;
3 4 0.000914 0.001883 0 0 0 0 1 0 1 -360 360;
4 5 0.002096 0.004304 0 0 0 0 1 0 1 -360 360;
5 6 0.000393 0.000807 0 0 0 0 1 0 1 -360 360;
6 7 0.000405 0.001380 0 0 0 0 1 0 1 -360 360;
7 8 0.001048 0.002152 0 0 0 0 1 0 1 -360 360;
7 9 0.002358 0.004842 0 0 0 0 1 0 1 -360 360;
7 10 0.001048 0.002152 0 0 0 0 1 0 1 -360 360;
84 11 0.000786 0.001614 0 0 0 0 1 0 1 -360 360;
11 12 0.003406 0.006944 0 0 0 0 1 0 1 -360 360;
12 13 0.000262 0.000538 0 0 0 0 1 0 1 -360 360;
12 14 0.000786 0.001614 0 0 0 0 1 0 1 -360 360;
84 15 0.001134 0.003864 0 0 0 0 1 0 1 -360 360;
15 16 0.000524 0.001076 0 0 0 0 1 0 1 -360 360;
16 17 0.000524 0.001076 0 0 0 0 1 0 1 -360 360;
17 18 0.001572 0.003228 0 0 0 0 1 0 1 -360 360;
18 19 0.000393 0.000807 0 0 0 0 1 0 1 -360 360;
19 20 0.001703 0.003497 0 0 0 0 1 0 1 -360 360;
20 21 0.002358 0.004842 0 0 0 0 1 0 1 -360 360;
21 22 0.001572 0.003228 0 0 0 0 1 0 1 -360 360;
21 23 0.001965 0.004035 0 0 0 0 1 0 1 -360 360;
23 24 0.001310 0.002690 0 0 0 0 1 0 1 -360 360;
84 25 0.000567 0.001932 0 0 0 0 1 0 1 -360 360;
25 26 0.001048 0.002152 0 0 0 0 1 0 1 -360 360;
26 27 0.002489 0.005111 0 0 0 0 1 0 1 -360 360;
27 28 0.000486 0.001656 0 0 0 0 1 0 1 -360 360;
28 29 0.001310 0.002690 0 0 0 0 1 0 1 -360 360;
84 30 0.001965 0.003960 0 0 0 0 1 0 1 -360 360;
30 31 0.001310 0.002690 0 0 0 0 1 0 1 -360 360;
31 32 0.001310 0.002690 0 0 0 0 1 0 1 -360 360;
32 33 0.000262 0.000538 0 0 0 0 1 0 1 -360 360;
33 34 0.001703 0.003497 0 0 0 0 1 0 1 -360 360;
34 35 0.000524 0.001076 0 0 0 0 1 0 1 -360 360;
35 36 0.004978 0.010222 0 0 0 0 1 0 1 -360 360;
36 37 0.000393 0.000807 0 0 0 0 1 0 1 -360 360;
37 38 0.000393 0.000807 0 0 0 0 1 0 1 -360 360;
38 39 0.000786 0.001614 0 0 0 0 1 0 1 -360 360;
39 40 0.002096 0.004304 0 0 0 0 1 0 1 -360 360;
38 41 0.001965 0.004035 0 0 0 0 1 0 1 -360 360;
41 42 0.002096 0.004304 0 0 0 0 1 0 1 -360 360;
84 43 0.000486 0.001656 0 0 0 0 1 0 1 -360 360;
43 44 0.000393 0.000807 0 0 0 0 1 0 1 -360 360;
44 45 0.001310 0.002690 0 0 0 0 1 0 1 -360 360;
45 46 0.002358 0.004842 0 0 0 0 1 0 1 -360 360;
84 47 0.002430 0.008280 0 0 0 0 1 0 1 -360 360;
47 48 0.000655 0.001345 0 0 0 0 1 0 1 -360 360;
48 49 0.000655 0.001345 0 0 0 0 1 0 1 -360 360;
49 50 0.000393 0.000807 0 0 0 0 1 0 1 -360 360;
50 51 0.000786 0.001614 0 0 0 0 1 0 1 -360 360;
51 52 0.000393 0.000807 0 0 0 0 1 0 1 -360 360;
52 53 0.000786 0.001614 0 0 0 0 1 0 1 -360 360;
53 54 0.000524 0.001076 0 0 0 0 1 0 1 -360 360;
54 55 0.001310 0.002690 0 0 0 0 1 0 1 -360 360;
84 56 0.002268 0.007728 0 0 0 0 1 0 1 -360 360;
56 57 0.005371 0.011029 0 0 0 0 1 0 1 -360 360;
57 58 0.000524 0.001076 0 0 0 0 1 0 1 -360 360;
58 59 0.000405 0.001380 0 0 0 0 1 0 1 -360 360;
59 60 0.000393 0.000807 0 0 0 0 1 0 1 -360 360;
60 61 0.000262 0.000538 0 0 0 0 1 0 1 -360 360;
61 62 0.001048 0.002152 0 0 0 0 1 0 1 -360 360;
62 63 0.002358 0.004842 0 0 0 0 1 0 1 -360 360;
63 64 0.000243 0.000828 0 0 0 0 1 0 1 -360 360;
84 65 0.000486 0.001656 0 0 0 0 1 0 1 -360 360;
65 66 0.001703 0.003497 0 0 0 0 1 0 1 -360 360;
66 67 0.001215 0.004140 0 0 0 0 1 0 1 -360 360;
67 68 0.002187 0.007452 0 0 0 0 1 0 1 -360 360;
68 69 0.000486 0.001656 0 0 0 0 1 0 1 -360 360;
69 70 0.000729 0.002484 0 0 0 0 1 0 1 -360 360;
70 71 0.000567 0.001932 0 0 0 0 1 0 1 -360 360;
71 72 0.000262 0.000528 0 0 0 0 1 0 1 -360 360;
84 73 0.003240 0.011040 0 0 0 0 1 0 1 -360 360;
73 74 0.000324 0.011040 0 0 0 0 1 0 1 -360 360;
74 75 0.000567 0.001932 0 0 0 0 1 0 1 -360 360;
75 76 0.000486 0.001656 0 0 0 0 1 0 1 -360 360;
84 77 0.002511 0.008556 0 0 0 0 1 0 1 -360 360;
77 78 0.001296 0.004416 0 0 0 0 1 0 1 -360 360;
78 79 0.000486 0.001656 0 0 0 0 1 0 1 -360 360;
79 80 0.001310 0.002640 0 0 0 0 1 0 1 -360 360;
80 81 0.001310 0.002690 0 0 0 0 1 0 1 -360 360;
81 82 0.000917 0.001833 0 0 0 0 1 0 1 -360 360;
82 83 0.003144 0.006456 0 0 0 0 1 0 1 -360 360;
5 55 0.001310 0.002690 0 0 0 0 1 0 1 -360 360;
7 60 0.001310 0.002690 0 0 0 0 1 0 1 -360 360;
11 43 0.001310 0.002690 0 0 0 0 1 0 1 -360 360;
12 72 0.003406 0.006994 0 0 0 0 1 0 1 -360 360;
13 76 0.004585 0.009415 0 0 0 0 1 0 1 -360 360;
14 18 0.005371 0.010824 0 0 0 0 1 0 1 -360 360;
16 26 0.000917 0.001883 0 0 0 0 1 0 1 -360 360;
20 83 0.000786 0.001614 0 0 0 0 1 0 1 -360 360;
28 32 0.000524 0.001076 0 0 0 0 1 0 1 -360 360;
29 39 0.000786 0.001614 0 0 0 0 1 0 1 -360 360;
34 46 0.000262 0.000538 0 0 0 0 1 0 1 -360 360;
40 42 0.001965 0.004035 0 0 0 0 1 0 1 -360 360;
53 64 0.000393 0.000807 0 0 0 0 1 0 1 -360 360;
];
% 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 = [
3 0.75 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
64 0.75 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
24 0.75 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
34 0.75 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
83 0.75 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
54 0.5 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
11 0.5 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
71 0.5 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
19 0.5 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
48 0.5 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
10 0.5 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
59 0.5 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
40 0.5 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
29 0.5 1 1 0 .1 100 1 1 .100 0 0 0 0 0 0 0 0 0 0 0;
75 0.5 1 1 0 .1 100 1 1 0.1 0 0 0 0 0 0 0 0 0 0 0;
84 0 0 10 1 1 100 1 7 1 0 0 0 0 0 0 0 0 0 0 0];
% 1 startup shutdown n x1 y1 ... xn
yn
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
2 11.57 0 4 0.000035 0.089 0;
2 0 0 3 0.0000058 0.021 0;
2 0 0 3 0.0000058 0.021 0;
2 0 0 3 0.0000058 0.021 0;
2 0 0 3 0.0000058 0.021 0;
2 0 0 3 0.0000058 0.021 0
2 0 0 3 0.0000053 0.02 0;
2 0 0 3 0.0000053 0.02 0;
2 0 0 3 0.0000053 0.02 0;
2 0 0 3 0.0000053 0.02 0;
2 0 0 3 0.0000053 0.02 0;
2 0 0 3 0.000005 0.02 0;
2 0 0 3 0.000005 0.02 0;
2 0 0 3 0.000005 0.02 0;
2 0 0 3 0.000005 0.02 0;
2 0 0 3 0.000005 0.02 0];
