http://www.pserc.cornell.edu/matpower/#pfconvergence
From: <bounce-121902459-73493...@list.cornell.edu> on behalf of Charalambos
Ioannou <xaralambo...@hotmail.com>
Reply-To: MATPOWER discussion forum <matpowe...@list.cornell.edu>
Date: Tuesday, October 3, 2017 at 9:12 AM
To: "MATPOWER-L@cornell.edu" <MATPOWER-L@cornell.edu>
Subject: PF and OPF did not converge
Dear Ray,
I am trying to run PF and OPF for the following case but every time it doesnt
converged the system. Please can you give me any advice??
The system consist 62 buses, 23 transformers, 8 photo voltaic plants and 1
Natural gus generator.
I consider all the buses as Type 1(PQ) and the natural gas generator as Type 2
(PV). The real power demand (Pd) for the photo voltaic plants i gave them a
negative sign.
function mpc = IsleExample3()
TCL MERGE ERROR ( 10/03/2017 10:44:16 ): "invalid command name "MATPOWER""
OutmailID: 121902634, List: 'matpower-l', MemberID: 75174736
SCRIPT: "MATPOWER Case Format : Version 2
mpc.version = '2';"
----- Power Flow Data ----- system MVA base
mpc.baseMVA = 100;
TCL MERGE ERROR ( 10/03/2017 10:44:16 ): "invalid command name "1""
OutmailID: 121902634, List: 'matpower-l', MemberID: 75174736
SCRIPT: "bus data
% bus_i type Pd Qd Gs Bs
area Vm
Va baseKV
zone Vmax
Vmin
mpc.bus = [
1 3 30.4 14.7 0 0 1 1 0 132 1 1.06 0.94;
2 1 30.4 14.7 0 0 1 1 0 132 1 1.06 0.94;
3 1 22.8 22.8 0 0 1 1 0 132 1 1.06 0.94;
4 1 22.8 22.8 0 0 1 1 0 132 1 1.06 0.94;
5 1 113 55 0 0 1 1 0 132 1 1.06 0.94;
6 1 113 55 0 0 1 1 0 132 1 1.06 0.94;
7 1 0 0 0 0 1 1 0 132 1 1.06 0.94;
8 1 113 55 0 0 1 1 0 132 1 1.06 0.94;
9 1 113 55 0 0 1 1 0 132 1 1.06 0.94;
10 1 113 55 0 0 1 1 0 33 1 1.06 0.94;
11 1 19 9 0 0 1 1 0 11 1 1.06 0.94;
12 1 113 55 0 0 1 1 0 33 1 1.06 0.94;
13 1 -6.89 0 0 0 1 1 0 33 1 1.06 0.94;
14 1 -5.00 0 0 0 1 1 0 33 1 1.06 0.94;
15 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
16 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
17 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
18 1 6 3 0 0 1 1 0 11 1 1.06 0.94;
19 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
20 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
21 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
22 1 -9.20 0 0 0 1 1 0 33 1 1.06 0.94;
23 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
24 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
25 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
26 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
27 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
28 1 17 8 0 0 1 1 0 11 1 1.06 0.94;
29 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
30 1 -5.00 0 0 0 1 1 0 33 1 1.06 0.94;
31 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
32 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
33 1 -5.00 0 0 0 1 1 0 33 1 1.06 0.94;
34 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
35 1 6 3 0 0 1 1 0 11 1 1.06 0.94;
36 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
37 1 -7.00 0 0 0 1 1 0 33 1 1.06 0.94;
38 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
39 1 6 3 0 0 1 1 0 11 1 1.06 0.94;
40 1 -1.5 0 0 0 1 1 0 33 1 1.06 0.94;
41 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
42 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
43 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
44 1 27 13 0 0 1 1 0 11 1 1.06 0.94;
45 1 27 13 0 0 1 1 0 11 1 1.06 0.94;
46 1 -5.00 0 0 0 1 1 0 33 1 1.06 0.94;
47 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
48 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
49 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
50 1 8 4 0 0 1 1 0 11 1 1.06 0.94;
51 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
52 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
53 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
54 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
55 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
56 1 12 6 0 0 1 1 0 11 1 1.06 0.94;
57 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
58 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
59 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
60 1 11 5 0 0 1 1 0 11 1 1.06 0.94;
61 1 0 0 0 0 1 1 0 33 1 1.06 0.94;
62 2 15.5 0 0 0 1 1 0 33 1 1.06 0.94;
];"
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 500 -500 1 100 1 0 0 0 0 0 0 0 0 0
0 0 0 0;
62 15.5 0 30 -30 1 100 1 30 0 0 0 0 0 0 0 0
0 0 0 0;
];
TCL MERGE ERROR ( 10/03/2017 10:44:16 ): "invalid command name "1""
OutmailID: 121902634, List: 'matpower-l', MemberID: 75174736
SCRIPT: "branch data
% fbus tbus r x b rateA
rateB rateC
ratio angle
status angmin
angmax
mpc.branch = [
1 2 0.0012 0.0221 0.0000 99 0 0 0 0 1 -360 360;
1 3 0.0019 0.0232 0.0000 138 0 0 0 0 1 -360 360;
2 4 0.0019 0.0232 0.0000 138 0 0 0 0 1 -360 360;
5 2 0.0012 0.0221 0.0000 124 0 0 0 0 1 -360 360;
3 8 0.0146 0.0303 0.0830 124 0 0 0 0 1 -360 360;
4 9 0.0145 0.0303 0.0808 124 0 0 0 0 1 -360 360;
5 6 0.0039 0.0277 0.0000 45 0 0 0 0 1 -360 360;
5 10 0 0.2771 0 45 0 0 0.95 0 1 -360 360;
6 10 0 0.2251 0 45 0 0 0.95 0 1 -360 360;
6 7 0.0019 0.0101 0.0174 145 0 0 0 0 1 -360 360;
6 7 0.0026 0.0078 0.0184 145 0 0 0 0 1 -360 360;
7 8 0.0040 0.0281 0.0000 0 0 0 0 0 1 -360 360;
7 9 0.0040 0.0281 0.0000 0 0 0 0 0 1 -360 360;
8 12 0 0.2501 0 90 0 0 0.95 0 1 -360 360;
9 12 0 0.1986 0 90 0 0 0.95 0 1 -360 360;
12 13 0.0083 0.0073 0.0006 24 0 0 0 0 1 -360 360;
12 14 0.0262 0.0231 0.0020 27 0 0 0 0 1 -360 360;
12 15 0.0389 0.1302 0.0002 44 0 0 0 0 1 -360 360;
15 16 0.0008 0.0012 0.0002 35 0 0 0 0 1 -360 360;
16 18 0 1.0800 0 15 0 0 0.955 0 1 -360 360;
15 20 0.0315 0.1054 0.0001 44 0 0 0 0 1 -360 360;
20 21 0.0584 0.0667 0.0001 19 0 0 0 0 1 -360 360;
21 22 0.0238 0.0210 0.0018 27 0 0 0 0 1 -360 360;
23 21 0.0262 0.0299 0.0000 19 0 0 0 0 1 -360 360;
27 20 0.0093 0.0310 0.0000 44 0 0 0 0 1 -360 360;
27 28 0 0.9930 0 30 0 0 0.955 0 1 -360 360;
12 19 0.0398 0.1333 0.0002 44 0 0 0 0 1 -360 360;
19 17 0.0016 0.0031 0.0000 31 0 0 0 0 1 -360 360;
17 18 0 1.0730 0 15 0 0 0.955 0 1 -360 360;
19 24 0.0398 0.1333 0.0002 44 0 0 0 0 1 -360 360;
25 24 0.0027 0.0020 0.0001 12 0 0 0 0 1 -360 360;
26 24 0.0093 0.0310 0.0000 44 0 0 0 0 1 -360 360;
26 28 0 0.9930 0 30 0 0 0.955 0 1 -360 360;
10 11 0 1.0300 0 30 0 0 0.955 0 1 -360 360;
10 11 0 1.0270 0 30 0 0 0.955 0 1 -360 360;
10 12 0.0944 0.1853 0.0002 31 0 0 0 0 1 -360 360;
10 12 0.0912 0.1790 0.0002 31 0 0 0 0 1 -360 360;
10 29 0.0087 0.0037 0.0002 31 0 0 0 0 1 -360 360;
10 29 0.0087 0.0037 0.0002 31 0 0 0 0 1 -360 360;
29 30 0.1005 0.2384 0.0010 30 0 0 0 0 1 -360 360;
31 30 0.0672 0.1601 0.0002 30 0 0 0 0 1 -360 360;
31 32 0.0320 0.0628 0.0001 31 0 0 0 0 1 -360 360;
32 33 0.0011 0.0010 0.0001 27 0 0 0 0 1 -360 360;
34 32 0.0844 0.1555 0.0012 31 0 0 0 0 1 -360 360;
34 35 0 1.0100 0 15 0 0 0.955 0 1 -360 360;
31 38 0.0542 0.2397 0.0000 44 0 0 0 0 1 -360 360;
31 39 0 0.9910 0 15 0 0 1 0 1 -360 360;
38 39 0 0.9910 0 15 0 0 1 0 1 -360 360;
38 37 0.0843 0.1655 0.0002 31 0 0 0 0 1 -360 360;
37 36 0.0261 0.0379 0.0015 31 0 0 0 0 1 -360 360;
36 35 0 1.0140 0 15 0 0 0.955 0 1 -360 360;
38 40 0.0012 0.0010 0.0001 24 0 0 0 0 1 -360 360;
38 41 0.0268 0.3173 0.0004 21 0 0 0 0 1 -360 360;
12 41 0.0445 0.0730 0.0018 31 0 0 0 0 1 -360 360;
41 44 0 0.7210 0 30 0 0 0.945 0 1 -360 360;
12 42 0.0397 0.0683 0.0035 31 0 0 0 0 1 -360 360;
42 45 0 0.7760 0 30 0 0 0.955 0 1 -360 360;
12 47 0.0481 0.1612 0.0002 44 0 0 0 0 1 -360 360;
47 46 0.0199 0.0368 0.0003 31 0 0 0 0 1 -360 360;
43 47 0.0586 0.1024 0.0029 25 0 0 0 0 1 -360 360;
43 45 0 0.7210 0 30 0 0 0.955 0 1 -360 360;
47 48 0.0667 0.1289 0.0004 31 0 0 0 0 1 -360 360;
48 49 0.1870 0.2211 0.0003 20 0 0 0 0 1 -360 360;
49 50 0 0.4670 0 30 0 0 0.955 0 1 -360 360;
48 57 0.0368 0.0722 0.0001 31 0 0 0 0 1 -360 360;
57 56 0 1.0920 0 30 0 0 0.96 0 1 -360 360;
12 54 0.1168 0.2780 0.0003 30 0 0 0 0 1 -360 360;
61 54 0.0676 0.1313 0.0006 31 0 0 0 0 1 -360 360;
61 60 0 1.0510 0 15 0 0 0.95 0 1 -360 360;
54 55 0.0278 0.0662 0.0001 30 0 0 0 0 1 -360 360;
55 56 0 1.0920 0 30 0 0 0.96 0 1 -360 360;
12 53 0.0629 0.2108 0.0003 44 0 0 0 0 1 -360 360;
53 62 0.0011 0.0022 0.0002 44 0 0 0 0 1 -360 360;
52 53 0.0111 0.0372 0.0000 44 0 0 0 0 1 -360 360;
52 51 0.1980 0.2343 0.0003 20 0 0 0 0 1 -360 360;
51 50 0 0.4770 0 30 0 0 0.955 0 1 -360 360;
52 58 0.0789 0.1688 0.0007 31 0 0 0 0 1 -360 360;
59 58 0.0016 0.0030 0.0000 27 0 0 0 0 1 -360 360;
58 60 0 1.0510 0 15 0 0 0.95 0 1 -360 360;
];"
Generator Cost Data
% model Startup Shutdown
n x1 y1 y2
mpc.gencost = [
2 0 0 3 0 1 0;
2 0 0 3 0 1 0;
];
Thanks in Advance.
Kind Regards ,
Charalampos
