I'm trying to run OPF for mpc attached. The runpf for the data was
successfully executed but the runopf('case28d') for the mpc did not
converge.
I seriously need a way out of this problem.
Thanks.
--
TIJANI MUHAMMED ADEKILEKUN
ELECTRICAL/ELECTRONICS ENGINEERING DEPARTMENT,
FEDERAL POLYTECHNIC, EDE, OSUN STATE.
function mpc = case28d
%CASE28d Power flow data for 37 bus test case.
% Please see CASEFORMAT for details on the case file format.
% (Using data obtained from Dr. Adepoju's students)
%
% 07/20/2015 UW ARCHIVE 100.0 1962 W IEEE 14 Bus Test Case
% MATPOWER
% $Id: case14.m 2394 2014-10-15 20:39:39Z ray $
%% 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 1 68.9 51.7 0 0 1 1.05 0 330 1 1.1 0.85;
2 1 0 0 0 0 1 1.05 0 330 1 1.1 0.85;
3 1 274.4 205.8 0 0 1 1 0 330 1 1.1 0.85;
4 1 344.7 258.5 0 0 1 1 0 330 1 1.1 0.85;
5 1 633.2 474.9 0 0 1 1 0 330 1 1.1 0.85;
6 1 13.8 10.3 0 0 1 1 0 330 1 1.1 0.85;
7 1 96.5 72.4 0 0 1 1 0 330 1 1.1 0.85;
8 1 383.3 287.5 0 0 1 1 0 330 1 1.1 0.85;
9 1 275.8 206.8 0 0 1 1 0 330 1 1.1 0.85;
10 1 201.2 150.9 0 0 1 1 0 330 1 1.1 0.85;
11 1 52.5 39.4 0 0 1 1.05 0 330 1 1.1 0.85;
12 1 427.0 320.2 0 0 1 1 0 330 1 1.1 0.85;
13 1 177.9 133.4 0 0 1 1 0 330 1 1.1 0.85;
14 1 184.6 138.4 0 0 1 1 0 330 1 1.1 0.85;
15 1 114.5 85.9 0 0 1 1 0 330 1 1.1 0.85;
16 1 130.6 97.9 0 0 1 1 0 330 1 1.1 0.85;
17 1 11.0 8.2 0 0 1 1 0 330 1 1.1 0.85;
18 1 0.000 0.000 0 0 1 1.05 0 330 1 1.1 0.85;
19 1 70.3 52.7 0 0 1 1 0 330 1 1.1 0.85;
20 1 193.0 144.7 0 0 1 1 0 330 1 1.1 0.85;
21 1 7.0 5.2 0 0 1 1.05 0 330 1 1.1 0.85;
22 1 199.8 149.9 0 0 1 1 0 330 1 1.1 0.85;
23 1 320.1 256.1 0 0 1 1.05 0 330 1 1.1 0.85;
24 1 68.9 51.7 0 0 1 1.05 0 330 1 1.1 0.85;
25 1 20.6 15.4 0 0 1 1 0 330 1 1.1 0.85;
26 1 0.000 0.000 0 0 1 1 0 330 1 1.1 0.85;
27 1 274.4 205.8 0 0 1 1.05 0 330 1 1.1 0.85;
28 1 344.7 258.5 0 0 1 1.05 0 330 1 1.1 0.85;
29 3 0 0 0 0 1 1.05 0 16 1 1.1 0.85;
30 2 0 0 0 0 1 1.05 0 11 1 1.1 0.85;
31 2 0 0 0 0 1 1.05 0 11 1 1.1 0.85;
32 2 0 0 0 0 1 1.05 0 16 1 1.1 0.85;
33 2 0 0 0 0 1 1.05 0 16 1 1.1 0.85;
34 2 0 0 0 0 1 1.05 0 16 1 1.1 0.85;
35 2 0 0 0 0 1 1.05 0 11 1 1.1 0.85;
36 2 0 0 0 0 1 1.05 0 16 1 1.1 0.85;
37 2 0 0 0 0 1 1.05 0 11 1 1.1 0.85;];
%% 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 = [
29 190.3 0 1000 -1000 1 100 1 1320 0 0 0 0
0 0 0 0 0 0 0 0;
30 670.0 0 1000 -1030 1.05 100 1 1000 150 0 0 0
0 0 0 0 0 0 0 0;
31 431.0 0 1000 -1000 1.05 100 1 600 0 0 0 0
0 0 0 0 0 0 0 0;
32 495.0 0 1050 -1050 1 100 1 600 0 0 0 0
0 0 0 0 0 0 0 0;
33 624.7 0 1080 -1010 1 100 1 1000 100 0 0 0
0 0 0 0 0 0 0 0;
34 388.9 0 1000 -1000 1.05 100 1 500 0 0 0 0
0 0 0 0 0 0 0 0;
35 94.9 0 1090 -1006 1 100 1 250 0 0 0 0
0 0 0 0 0 0 0 0;
36 750.0 0 1090 -1006 1 100 1 1000 0 0 0 0
0 0 0 0 0 0 0 0;
37 750.0 0 1090 -1006 1 100 1 1000 0 0 0 0
0 0 0 0 0 0 0 0;];
%mpc.gen = [
% 29 830.23 0 450 -255 1 100 1 1320 0;
% 30 200 0 450 -250 0.99 100 1 500 100;
% 31 300 0 450 -250 1 100 1 400 30;
% 32 250 0 450 -250 1 100 1 280 0;
% 33 490 0 450 -250 1.03 100 1 600 0;
% 34 350 0 450 -250 1.04 100 1 540 0;
% 35 450 0 450 -250 1.03 100 1 560 0;
% 36 300 0 450 -250 1 100 1 400 30;
% 37 300 0 450 -250 1 100 1 400 30;];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio
angle status angmin angmax
mpc.branch = [
3 1 0.0006 0.0044 0.0295 2000 2000 0.975 0 0 1
-360 360;
4 5 0.0007 0.0050 0.0333 1000 1000 1000 0 0 1
-360 360;
1 5 0.0023 0.0176 0.1176 1000 1000 1000 0 0 1
-360 360;
5 8 0.0110 0.0828 0.5500 1000 1000 1000 0 0 1
-360 360;
5 9 0.0054 0.0405 0.2669 1000 1000 1000 0 0 1
-360 360;
5 10 0.0099 0.0745 0.4949 1000 1000 1000 0 0 1
-360 360;
6 8 0.0077 0.0576 0.3830 1000 1000 1000 0 0 1
-360 360;
2 8 0.0043 0.0317 0.2101 1000 1000 1000 0 0 1
-360 360;
2 7 0.0012 0.0089 0.0589 1000 1000 1000 0 0 1
-360 360;
7 24 0.0025 0.0186 0.1237 1520 1520 1520 0 0 1
-360 360;
8 14 0.0054 0.0405 0.2691 1520 1520 1520 0 0 1
-360 360;
8 10 0.0098 0.0742 0.4930 1520 1520 1520 0 0 1
-360 360;
8 24 0.0020 0.0148 0.0982 1520 1520 1520 0 0 1
-360 360;
9 10 0.0045 0.0340 0.2257 760 760 760 0 0 1
-360 360;
15 21 0.0122 0.0916 0.6089 760 760 760 0 0 1
-360 360;
10 17 0.0061 0.0461 0.3064 760 760 760 0 0 1
-360 360;
11 12 0.0010 0.0074 0.0491 760 760 760 0 0 1
-360 360;
12 14 0.0060 0.0455 0.3025 760 760 760 0 0 1
-360 360;
13 14 0.0036 0.0272 0.1807 1520 1520 1520 0 0 1
-360 360;
16 19 0.0118 0.0887 0.5892 760 760 760 0 0 1
-360 360;
17 18 0.0002 0.0020 0.0098 760 760 760 0 0 1
-360 360;
17 23 0.0096 0.0721 0.4793 2280 2280 2280 0 0 1
-360 360;
17 21 0.0032 0.0239 0.1589 1520 1520 1520 0 0 1
-360 360;
19 20 0.0081 0.0609 0.4046 1520 1520 1520 0 0 1
-360 360;
20 22 0.0090 0.0680 0.4516 760 760 760 0 0 1
-360 360;
20 23 0.0038 0.0248 0.1886 1520 1520 1520 0 0 1
-360 360;
23 26 0.0038 0.0248 0.1886 1000 1000 1000 0 0 1
-360 360;
12 25 0.0071 0.0532 0.3800 1000 1000 1000 0 0 1
-360 360;
19 25 0.0059 0.0443 0.3060 1520 1520 1520 0 0 1
-360 360;
25 27 0.0079 0.0591 0.3900 760 760 760 0 0 1
-360 360;
5 28 0.0016 0.0118 0.0932 1000 1000 1000 0 0 1
-360 360;
1 29 0.0000 0.00025 0.0000 1000 1000 1000 1.975 0 1
-360 360;
2 30 0.0000 0.00086 0.0000 1000 1000 1000 1 0 1
-360 360;
11 31 0.0000 0.00076 0.0000 760 760 760 0.950 0 1
-360 360;
18 32 0.0000 0.00009 0.0000 760 760 760 0.950 0 1
-360 360;
21 33 0.0000 0.00080 0.0000 1520 1520 1520 1.025 0 1
-360 360;
23 34 0.0000 0.00084 0.0000 760 760 760 1 0 1
-360 360;
24 35 0.0000 0.00032 0.0000 760 760 760 1 0 1
-360 360;
27 36 0.0000 0.00032 0.0000 1520 1520 1520 0.950 0 1
-360 360;
28 37 0.0000 0.00032 0.0000 1520 1520 1520 0.950 0 1
-360 360;];
%% generator cost data
% 1 startup shutdown n x1 y1 ... xn
yn
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
2 0 0 3 0.0430 20 0;
2 0 0 3 0.25 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.25 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;];
% Warnings from cdf2matp conversion:
%
% ***** check the title format in the first line of the cdf file.
% ***** Qmax = Qmin at generator at bus 1 (Qmax set to Qmin + 10)
% ***** MVA limit of branch 1 - 2 not given, set to 0
% ***** MVA limit of branch 1 - 5 not given, set to 0
% ***** MVA limit of branch 2 - 3 not given, set to 0
% ***** MVA limit of branch 2 - 4 not given, set to 0
% ***** MVA limit of branch 2 - 5 not given, set to 0
% ***** MVA limit of branch 3 - 4 not given, set to 0
% ***** MVA limit of branch 4 - 5 not given, set to 0
% ***** MVA limit of branch 4 - 7 not given, set to 0
% ***** MVA limit of branch 4 - 9 not given, set to 0
% ***** MVA limit of branch 5 - 6 not given, set to 0
% ***** MVA limit of branch 6 - 11 not given, set to 0
% ***** MVA limit of branch 6 - 12 not given, set to 0
% ***** MVA limit of branch 6 - 13 not given, set to 0
% ***** MVA limit of branch 7 - 8 not given, set to 0
% ***** MVA limit of branch 7 - 9 not given, set to 0
% ***** MVA limit of branch 9 - 10 not given, set to 0
% ***** MVA limit of branch 9 - 14 not given, set to 0
% ***** MVA limit of branch 10 - 11 not given, set to 0
% ***** MVA limit of branch 12 - 13 not given, set to 0
% ***** MVA limit of branch 13 - 14 not given, set to 0
