Dear Power Globers,
Thank you Idris for your input I have updated the case file and it
converged.
I am working on changing topology of the network and finding out it's
effects.If suppose I want to use the tie lines and for restraining radility
opens five section switches say [3 43 48 24 18]. As I typed
>> op=[3 43 48 24 18];
>> temp=loadcase(case69);
>> for i=1:5
temp.branch(op(i),11)=0
end
>> runpf(temp)
It shows a warning as
Warning: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 1.724770e-017.
> In newtonpf at 108
In runpf at 224
I am here attaching the case file. Please help me in this.
Thanks and regards
Pinak
Kolkata,India
On Sun, Jan 11, 2015 at 8:19 PM, Idris Musa (PGR) <
[email protected]> wrote:
> Hi Panda,
>
> The real and the reactive loads (PD and QD) in column 3 and 4 of mpc.bus
> should be in MW and not kW as you currently have.
>
> Simply divide the current values representing PD and QD in columns 3 and 4
> by 1000, and the power flow will converge.
>
>
>
> Idris
>
>
>
> *From:* [email protected] [mailto:
> [email protected]] *On Behalf Of *Pinak Panda
> *Sent:* 11 January 2015 05:48
> *To:* [email protected]
> *Subject:* non convergence of 69 bus radial distribution system
>
>
>
> Hi all,
>
>
>
> I am trying to converge the power flow of 69 bus radial distribution by NR
> method here is the case file I have created.the error shows, It is not
> getting converged in 10 iterations.
>
>
>
> Please help me in this regard.
>
>
>
> Thanking you,
>
> Pinak Panda
>
> M.Tech-scholar
>
> SRIT Kolkata
>
> India
>
function mpc = case69
%% MATPOWER Case Format : Version 2
mpc.version = '2';
%%----- Power Flow Data -----%%
%% system MVA base
mpc.baseMVA = 100;
%% 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 12.66 1 1 0.90;
2 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
3 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
4 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
5 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
6 1 2.6/1000 2.2/1000 0 0 1 1 0 12.66 1
1 0.90;
7 1 40.4/1000 30/1000 0 0 1 1 0 12.66 1
1 0.90;
8 1 75/1000 54/1000 0 0 1 1 0 12.66 1
1 0.90;
9 1 30/1000 22/1000 0 0 1 1 0 12.66 1
1 0.90;
10 1 28/1000 19/1000 0 0 1 1 0 12.66 1
1 0.90;
11 1 145/1000 104/1000 0 0 1 1 0 12.66 1
1 0.90;
12 1 145/1000 104/1000 0 0 1 1 0 12.66 1
1 0.90;
13 1 8/1000 5/1000 0 0 1 1 0 12.66 1
1 0.90;
14 1 8/1000 5.5/1000 0 0 1 1 0 12.66 1
1 0.90;
15 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
16 1 45.5/1000 30/1000 0 0 1 1 0 12.66 1
1 0.90;
17 1 60/1000 35/1000 0 0 1 1 0 12.66 1
1 0.90;
18 1 60/1000 35/1000 0 0 1 1 0 12.66 1
1 0.90;
19 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
20 1 1/1000 0.6/1000 0 0 1 1 0 12.66 1
1 0.90;
21 1 114/1000 81/1000 0 0 1 1 0 12.66 1
1 0.90;
22 1 5/1000 3.5/1000 0 0 1 1 0 12.66 1
1 0.90;
23 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
24 1 28/1000 20/1000 0 0 1 1 0 12.66 1
1 0.90;
25 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
26 1 14/1000 10/1000 0 0 1 1 0 12.66 1
1 0.90;
27 1 14/1000 10/1000 0 0 1 1 0 12.66 1
1 0.90;
28 1 26/1000 18.6/1000 0 0 1 1 0 12.66 1
1 0.90;
29 1 26/1000 18.6/1000 0 0 1 1 0 12.66 1
1 0.90;
30 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
31 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
32 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
33 1 14/1000 10/1000 0 0 1 1 0 12.66 1
1 0.90;
34 1 19.5/1000 14/1000 0 0 1 1 0 12.66 1
1 0.90;
35 1 6/1000 4/1000 0 0 1 1 0 12.66 1
1 0.90;
36 1 26/1000 18.55/1000 0 0 1 1 0 12.66 1
1 0.90;
37 1 26/1000 18.55/1000 0 0 1 1 0 12.66 1
1 0.90;
38 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
39 1 24/1000 17/1000 0 0 1 1 0 12.66 1
1 0.90;
40 1 24/1000 17/1000 0 0 1 1 0 12.66 1
1 0.90;
41 1 1.2/1000 1/1000 0 0 1 1 0 12.66 1
1 0.90;
42 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
43 1 6/1000 4.3/1000 0 0 1 1 0 12.66 1
1 0.90;
44 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
45 1 39.22/1000 26.3/1000 0 0 1 1 0 12.66 1
1 0.90;
46 1 39.22/1000 26.3/1000 0 0 1 1 0 12.66 1
1 0.90;
47 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
48 1 79/1000 56.4/1000 0 0 1 1 0 12.66 1
1 0.90;
49 1 384.7/1000 274.5/1000 0 0 1 1 0 12.66 1
1 0.90;
50 1 384.7/1000 274.5/1000 0 0 1 1 0 12.66 1
1 0.90;
51 1 40.5/1000 28.3/1000 0 0 1 1 0 12.66 1
1 0.90;
52 1 3.6/1000 2.7/1000 0 0 1 1 0 12.66 1
1 0.90;
53 1 4.35/1000 3.5/1000 0 0 1 1 0 12.66 1
1 0.90;
54 1 26.4/1000 19/1000 0 0 1 1 0 12.66 1
1 0.90;
55 1 24/1000 17.2/1000 0 0 1 1 0 12.66 1
1 0.90;
56 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
57 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
58 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
59 1 100/1000 72/1000 0 0 1 1 0 12.66 1
1 0.90;
60 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
61 1 1244/1000 888/1000 0 0 1 1 0 12.66 1
1 0.90;
62 1 32/1000 23/1000 0 0 1 1 0 12.66 1
1 0.90;
63 1 0 0 0 0 1 1 0 12.66 1 1 0.90;
64 1 227/1000 162/1000 0 0 1 1 0 12.66 1
1 0.90;
65 1 59/1000 42/1000 0 0 1 1 0 12.66 1
1 0.90;
66 1 18/1000 13/1000 0 0 1 1 0 12.66 1
1 0.90;
67 1 18/1000 13/1000 0 0 1 1 0 12.66 1
1 0.90;
68 1 28/1000 20/1000 0 0 1 1 0 12.66 1
1 0.90;
69 1 28/1000 20/1000 0 0 1 1 0 12.66 1
1 0.90;
];
%% 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 0 0 0 1.00 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.0005*0.624 0.0012*0.624 0 0 0
0 0 0 1 -360 360;
2 3 0.0005*0.624 0.0012*0.624 0 0 0
0 0 0 1 -360 360;
3 4 0.0015*0.624 0.0036*0.624 0 0 0
0 0 0 1 -360 360;
4 5 0.0251*0.624 0.0294*0.624 0 0 0
0 0 0 1 -360 360;
5 6 0.366*0.624 0.1864*0.624 0 0 0
0 0 0 1 -360 360;
6 7 0.381*0.624 0.1941*0.624 0 0 0
0 0 0 1 -360 360;
7 8 0.0922*0.624 0.0470*0.624 0 0 0
0 0 0 1 -360 360;
8 9 0.0493*0.624 0.0251*0.624 0 0 0
0 0 0 1 -360 360;
9 10 0.8190*0.624 0.2707*0.624 0 0 0
0 0 0 1 -360 360;
10 11 0.1872*0.624 0.0619*0.624 0 0 0
0 0 0 1 -360 360;
11 12 0.7114*0.624 0.2351*0.624 0 0 0
0 0 0 1 -360 360;
12 13 1.03*0.624 0.34*0.624 0 0 0
0 0 0 1 -360 360;
13 14 1.044*0.624 0.345*0.624 0 0 0
0 0 0 1 -360 360;
14 15 1.0580*0.624 0.3496*0.624 0 0 0
0 0 0 1 -360 360;
15 16 0.1966*0.624 0.0650*0.624 0 0 0
0 0 0 1 -360 360;
16 17 0.3744*0.624 0.1238*0.624 0 0 0
0 0 0 1 -360 360;
17 18 0.0047*0.624 0.0016*0.624 0 0 0
0 0 0 1 -360 360;
18 19 0.3276*0.624 0.1083*0.624 0 0 0
0 0 0 1 -360 360;
19 20 0.2106*0.624 0.0690*0.624 0 0 0
0 0 0 1 -360 360;
20 21 0.3416*0.624 0.1129*0.624 0 0 0
0 0 0 1 -360 360;
21 22 0.0140*0.624 0.0046*0.624 0 0 0
0 0 0 1 -360 360;
22 23 0.1591*0.624 0.0526*0.624 0 0 0
0 0 0 1 -360 360;
23 24 0.3463*0.624 0.1145*0.624 0 0 0
0 0 0 1 -360 360;
24 25 0.7488*0.624 0.2475*0.624 0 0 0
0 0 0 1 -360 360;
25 26 0.3089*0.624 0.1021*0.624 0 0 0
0 0 0 1 -360 360;
26 27 0.1732*0.624 0.0572*0.624 0 0 0
0 0 0 1 -360 360;
3 28 0.0044*0.624 0.0108*0.624 0 0 0
0 0 0 1 -360 360;
28 29 0.0640*0.624 0.1565*0.624 0 0 0
0 0 0 1 -360 360;
29 30 0.3978*0.624 0.1315*0.624 0 0 0
0 0 0 1 -360 360;
30 31 0.0702*0.624 0.0232*0.624 0 0 0
0 0 0 1 -360 360;
31 32 0.3510*0.624 0.1160*0.624 0 0 0
0 0 0 1 -360 360;
32 33 0.8390*0.624 0.2816*0.624 0 0 0
0 0 0 1 -360 360;
33 34 1.7080*0.624 0.5646*0.624 0 0 0
0 0 0 1 -360 360;
34 35 1.4740*0.624 0.4873*0.624 0 0 0
0 0 0 1 -360 360;
3 36 0.0044*0.624 0.0108*0.624 0 0 0
0 0 0 1 -360 360;
36 37 0.0640*0.624 0.1565*0.624 0 0 0
0 0 0 1 -360 360;
37 38 0.1053*0.624 0.1230*0.624 0 0 0
0 0 0 1 -360 360;
38 39 0.0304*0.624 0.0355*0.624 0 0 0
0 0 0 1 -360 360;
39 40 0.0018*0.624 0.0021*0.624 0 0 0
0 0 0 1 -360 360;
40 41 0.7283*0.624 0.8509*0.624 0 0 0
0 0 0 1 -360 360;
41 42 0.3100*0.624 0.3623*0.624 0 0 0
0 0 0 1 -360 360;
42 43 0.0410*0.624 0.0478*0.624 0 0 0
0 0 0 1 -360 360;
43 44 0.0092*0.624 0.0116*0.624 0 0 0
0 0 0 1 -360 360;
44 45 0.1089*0.624 0.1373*0.624 0 0 0
0 0 0 1 -360 360;
45 46 0.0009*0.624 0.0012*0.624 0 0 0
0 0 0 1 -360 360;
4 47 0.0034*0.624 0.0084*0.624 0 0 0
0 0 0 1 -360 360;
47 48 0.0851*0.624 0.2083*0.624 0 0 0
0 0 0 1 -360 360;
48 49 0.2898*0.624 0.7091*0.624 0 0 0
0 0 0 1 -360 360;
49 50 0.0822*0.624 0.2011*0.624 0 0 0
0 0 0 1 -360 360;
8 51 0.0928*0.624 0.0473*0.624 0 0 0
0 0 0 1 -360 360;
51 52 0.3319*0.624 0.1114*0.624 0 0 0
0 0 0 1 -360 360;
9 53 0.1740*0.624 0.0886*0.624 0 0 0
0 0 0 1 -360 360;
53 54 0.2030*0.624 0.1034*0.624 0 0 0
0 0 0 1 -360 360;
54 55 0.2842*0.624 0.1447*0.624 0 0 0
0 0 0 1 -360 360;
55 56 0.2813*0.624 0.1433*0.624 0 0 0
0 0 0 1 -360 360;
56 57 1.5900*0.624 0.5337*0.624 0 0 0
0 0 0 1 -360 360;
57 58 0.7837*0.624 0.2630*0.624 0 0 0
0 0 0 1 -360 360;
58 59 0.3042*0.624 0.1006*0.624 0 0 0
0 0 0 1 -360 360;
59 60 0.3861*0.624 0.1172*0.624 0 0 0
0 0 0 1 -360 360;
60 61 0.5075*0.624 0.2585*0.624 0 0 0
0 0 0 1 -360 360;
61 62 0.0974*0.624 0.0496*0.624 0 0 0
0 0 0 1 -360 360;
62 63 0.1450*0.624 0.0738*0.624 0 0 0
0 0 0 1 -360 360;
63 64 0.7105*0.624 0.3619*0.624 0 0 0
0 0 0 1 -360 360;
64 65 1.0410*0.624 0.5302*0.624 0 0 0
0 0 0 1 -360 360;
11 66 0.2012*0.624 0.0611*0.624 0 0 0
0 0 0 1 -360 360;
66 67 0.0047*0.624 0.0014*0.624 0 0 0
0 0 0 1 -360 360;
12 68 0.7394*0.624 0.2444*0.624 0 0 0
0 0 0 1 -360 360;
68 69 0.0047*0.624 0.0016*0.624 0 0 0
0 0 0 1 -360 360;
52 67 0.0047*0.624 0.0016*0.624 0 0 0
0 0 0 1 -360 360;
15 69 0.0047*0.624 0.0016*0.624 0 0 0
0 0 0 1 -360 360;
10 65 0.0047*0.624 0.0016*0.624 0 0 0
0 0 0 1 -360 360;
46 50 0.0047*0.624 0.0016*0.624 0 0 0
0 0 0 1 -360 360;
35 27 0.0047*0.624 0.0016*0.624 0 0 0
0 0 0 1 -360 360;
% %
];
