Dear Mr.Ray,
I have a problem with load flow solution of particular case file. Here I
am attaching the case file with the mail.
In my case file, I created some faults in the network(Branches (9,10)
(10,11) and (11,12) are turned 'off') such that the load at bus 10 is
not supplied.
when I run the powerflow....runpf('mycase1'), The solution looks like as
follows
================================================================================
| Bus Data
|
================================================================================
Bus Voltage Generation Load
# Mag(pu) Ang(deg) P (MW) Q (MVAr) P (MW) Q (MVAr)
----- ------- -------- -------- -------- -------- --------
1 1.000 0.000 22.85 7.35 - -
2 0.969 -2.611 - - 3.31 1.79
3 0.960 -2.776 - - 2.13 1.15
4 0.946 -3.844 - - 2.56 1.38
5 0.963 -2.751 - - 0.93 0.50
6 0.962 -2.776 - - - -
7 0.961 -2.784 - - - -
8 0.957 -3.023 - - 2.07 1.12
9 0.948 -3.901 - - 2.42 1.31
10 NaN NaN - - 1.26 0.68
11 NaN NaN - - - -
12 0.981 -1.623 - - 1.33 0.78
13 0.992 -0.969 - - 5.34 2.99
14 1.000 -1.524 3.00 7.99 - -
15 0.982 -1.536 - - 2.67 1.44
16 1.000 -1.524 - - - -
17 0.972 -1.819 - - - -
18 0.964 -2.192 - - 2.67 1.44
-------- -------- -------- --------
Total: 25.85 15.34 26.69 14.58
Does this mean we can't solve the power flow.?
Can't the Matpower take the connected part of the network only and do
load flow for that part?
Thanks in advance.
Sridhar Chouhan.
function [baseMVA, bus, gen, branch, areas, gencost] = mycase
%CASE9 Power flow data for 9 bus, 3 generator case.
% Please see 'help caseformat' for details on the case file format.
%
% Based on data from Joe H. Chow's book, p. 70.
% MATPOWER
% $Id: case9.m,v 1.7 2007/09/17 16:07:48 ray Exp $
%%----- Power Flow Data -----%%
%% system MVA base
baseMVA = 100;
%% bus data
% bus_i type Pd Qd Gs Bs area Vm
Va baseKV zone Vmax Vmin
bus = [
1 3 0 0 0 0 1 1 0
12.47 1 1.1 0.9;
2 1 3.3100 1.7874 0 0 1 1
0 12.47 1 1.1 0.9;
3 1 2.1300 1.1502 0 0 1 1
0 12.47 1 1.1 0.9;
4 1 2.5600 1.3800 0 0 1 1 0
12.47 1 1.1 0.9;
5 1 0.9300 0.5022 0 0 1 1 0
12.47 1 1.1 0.9;
6 1 0 0 0 0 1 1 0
12.47 1 1.1 0.9;
7 1 0 0 0 0 1 1 0
12.47 1 1.1 0.9;
8 1 2.0700 1.1178 0 0 1 1 0
12.47 1 1.1 0.9;
9 1 2.4200 1.3064 0 0 1 1 0
12.47 1 1.1 0.9;
10 1 1.2600 0.6804 0 0 1 1 0
12.47 1 1.1 0.9;
11 1 0 0 0 0 1 1 0
12.47 1 1.1 0.9;
12 1 1.3300 0.7812 0 0 1 1 0
12.47 1 1.1 0.9;
13 1 5.3400 2.9884 0 0 1 1 0
12.47 1 1.1 0.9;
14 2 0 0 0 0 1 1 0
12.47 1 1.1 0.9;
15 1 2.6700 1.4418 0 0 1 1 0
12.47 1 1.1 0.9;
16 1 0 0 0 0 1 1 0
12.47 1 1.1 0.9;
17 1 0 0 0 0 1 1 0
12.47 1 1.1 0.9;
18 1 2.6700 1.4418 0 0 1 1 0
12.47 1 1.1 0.9;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax
Pmin
gen = [
1 0 0 100 -100 1 100 1 100
0;
14 3 0 20 -20 1 100 1 10
0;
];
%% branch data
% fbus tbus r x b rateArateBrateC ratio angle
status
branch = [
1 2 0.0500 0.5760 0 250 250 250 0
0 1;
1 5 0.1270 0.9200 0 250 250 250 0
0 1;
1 13 0.1039 0.1700 0 150 150 150 0 0
1;
2 3 0.1090 0.1170 0 150 150 150 0
0 1;
2 7 0.0010 0.1100 0 300 300 300 0
0 0;%S1
3 4 0.1085 0.7200 0 150 150 150 0
0 1;
3 8 0.0010 0.1100 0 250 250 250 0
0 0;%S2
4 9 0.0010 0.1100 0 250 250 250 0
0 0;%S3
5 6 0.0132 0.0161 0 250 250 250 0
0 1;
6 7 0.0101 0.0085 0 250 250 250 0
0 1;
6 11 0.0010 0.1100 0 250 250 250 0 0
0;%S4
6 18 0.0010 0.1100 0 250 250 250 0
0 0;%S6
7 8 0.0390 0.1070 0 150 150 150 0
0 1;
8 9 0.0220 0.5860 0 300 300 300 0
0 1;
9 10 0.0119 0.1080 0 150 150 150 0
0 0;
10 11 0.1085 0.1720 0 250 250 250 0
0 0;
11 12 0.1023 0.1625 0 250 250 250 0
0 0;
12 15 0.0010 0.1100 0 250 250 250 0
0 1;%S5
13 14 0.0320 0.2161 0 250 250 250 0 0
1;
14 15 0.2032 0.1161 0 250 250 250 0 0
1;
14 16 0.2032 0.1161 0 250 250 250 0 0
1;
15 17 0.2001 0.2850 0 250 250 250 0
0 1;
16 17 0.0010 0.1100 0 250 250 250 0 0
0;%S7
17 18 0.1170 0.2920 0 250 250 250 0
0 1;
];
%----- OPF Data -----%%
%% area data
areas = [
1 5;
];
%% generator cost data
% 1 startup shutdown n x1 y1 ... xn
yn
% 2 startup shutdown n c(n-1) ... c0
% gencost = [
% 2 1500 0 3 0.11 5 150;
% 2 2000 0 3 0.085 1.2 600;
% 2 3000 0 3 0.1225 1 335;
% ];
gencost=[];
return;
