Dear Ray,
Thanks for your reply. I have prepared a 16 bus system. The power flow
program runs OK with that but when I go for Optimal Power Flow it shows it
did not convergence. I have attached the file below here. Please would you
check the file and direct me what to change ? Waiting to hear from you.
Thanks
On Jan 6, 2017 1:33 AM, "Ray Zimmerman" <[email protected]> wrote:
> See Section 2.3.1 and Appendix B in the User’s Manual
> <http://www.pserc.cornell.edu/matpower/manual.pdf>. It’s probably easiest
> to start with an existing case file and modify it for your data. Let us
> know if you have a specific question about the case format.
>
> Ray
>
>
>
> On Jan 2, 2017, at 12:33 PM, siddique tomal <[email protected]>
> wrote:
>
> Dear all,
> i need to create a new case file of 16 bus system. i have attached the
> data table i have now. but my asking is how can i develop a new case file
> with the following data. its important to mention that the system baseMVA
> is 100 MVA and base kva of 23 kV.
> Thanks to all.
> <16 bus data.png>
>
>
>
function mpc = case_radial16
%IEEE 69 RADIAL DISTRIBUTION SYSTEM Power flow data for IEEE 69 radial distribution test case.
% Please see CASEFORMAT for details on the case file format.
% This data was converted from IEEE Common Data Format
% (ieee30cdf.txt) on 20-Sep-2004 by cdf2matp, rev. 1.11
% See end of file for warnings generated during conversion.
%
%
% MATPOWER
% $Id: case_radial69.m,v 1.8 2013/09/11 17:06:07 Koyanagi Exp $
%% MATPOWER Case Format : Version 2
mpc.version = '2'; %OK
%%----- 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.000 0 23 1 1.05 0.95;
2 1 0 0 0 0 1 1.000 0 23 1 1.05 0.95;
3 1 0 0 0 0 1 1.000 0 23 1 1.05 0.95;
4 1 2 1.6 0 0 1 0.9907 -0.3968 23 1 1.05 0.95;
5 1 3 1.5 0 0 1 0.9878 -0.5443 23 1 1.05 0.95;
6 1 2 0.8 0 0 1 0.9860 -0.6972 23 1 1.05 0.95;
7 1 1.5 1.2 0 0 1 0.9849 -0.7043 23 1 1.05 0.95;
8 1 4 2.7 0 0 1 0.9791 -0.7635 23 1 1.05 0.95;
9 1 5 3.0 0 0 1 0.9711 -1.452 23 1 1.05 0.95;
10 1 1 0.9 0 0 1 0.9769 -0.7701 23 1 1.05 0.95;
11 1 0.6 0.1 0 0 1 0.9710 -1.526 23 1 1.05 0.95;
12 1 4.5 2.0 0 0 1 0.9693 -1.837 23 1 1.05 0.95;
13 1 1 0.9 0 0 1 0.9944 -0.3293 23 1 1.05 0.95;
14 1 1 0.7 0 0 1 0.9948 -0.4562 23 1 1.05 0.95;
15 1 1 0.9 0 0 1 0.9918 -0.5228 23 1 1.05 0.95;
16 1 2.1 1.0 0 0 1 0.9913 -0.5904 23 1 1.05 0.95;
];
%% 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 500 0 1000 -100 1.0 10 1 1000 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 4 0.075 0.100 0 0 0 0 0 0 1 -360 360;
4 5 0.080 0.110 0 0 0 0 0 0 1 -360 360;
4 6 0.090 0.180 0 0 0 0 0 0 1 -360 360;
6 7 0.040 0.040 0 0 0 0 0 0 1 -360 360;
2 8 0.110 0.110 0 0 0 0 0 0 1 -360 360;
8 9 0.080 0.110 0 0 0 0 0 0 1 -360 360;
8 10 0.110 0.110 0 0 0 0 0 0 1 -360 360;
9 11 0.110 0.110 0 0 0 0 0 0 1 -360 360;
9 12 0.080 0.110 0 0 0 0 0 0 1 -360 360;
3 13 0.110 0.110 0 0 0 0 0 0 1 -360 360;
13 14 0.090 0.120 0 0 0 0 0 0 1 -360 360;
13 15 0.080 0.110 0 0 0 0 0 0 1 -360 360;
15 16 0.040 0.040 0 0 0 0 0 0 1 -360 360;
10 15 0.040 0.040 0 0 0 0 0 0 1 -360 360;
5 11 0.040 0.040 0 0 0 0 0 0 1 -360 360;
];
%%----- OPF Data -----%%
%% 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.038432 20 0;
];
% Warnings from cdf2matp conversion:
%
% ***** b not given, set to 0
% ***** MVA limit of all branch not given, set to 9900
