Dear All
I am trying to solve the power flow problem in a radial distribution
system. However, due to the high level of R/X it is not working (when I
divide R by 5 it works). I have tried the 4 PF algorithms implemented in
MATPOWER without any success. The data of this network is attached).
Please give me some advises to fix this issue.
Best Regards
Dr Houssem
% The main Programme
opt = mpoption('VERBOSE', 0, 'OUT_ALL', 1,'PF_ALG',2,'PF_MAX_IT_FD',1000);
[MVAbase, bus, gen, branch, success, et] = runpf('case3',opt);
function[baseMVA, bus, gen, branch] = case3
% CASE3 Power flow data for 12 bus, 11 gen case
%% system MVA base
baseMVA = 100;
%%bus data
%bus_i typePdQdGsBsareaVmVabaseKVzoneVmaxVmin
bus =[
13000011013.811.1 0.9;
210.60.60011013.811.1 0.9;
320.40.30011013.811.1 0.9;
420.550.550011013.811.1 0.9;
520.30.30011013.811.1 0.9;
620.20.150011013.811.1 0.9;
720.550.550011013.811.1 0.9;
820.450.450011013.811.1 0.9;
920.40.40011013.811.1 0.9;
1020.350.30011013.811.1 0.9;
1120.40.30011013.811.1 0.9;
1220.150.150011013.811.1 0.9;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
gen = [
12 5 0 .2 -.21 100 0 1.1 0;
11 5 0 .2 -.21 100 0 1.1 0;
10 5 0 .2 -.21 100 0 1.1 0;
95 0 .2 -.21 100 0 1.1 0;
85 0 .2 -.21 100 0 1.1 0;
75 0 .2 -.21 100 0 1.1 0;
65 0 .2 -.21 100 0 1.1 0;
55 0 .2 -.21 100 0 1.1 0;
45 0 .2 -.21 100 0 1.1 0;
35 0 .2 -.21 100 0 1.1 0;
10 0 100 -100 1 100 1 0 0;
];
%%branchdata
%fbustbusrxbrateArateBrateCratioanglestatus
branch=[
121.0930.4550250 250 250 001;
231.1840.4940250 250 250 001;
342.0950.8730250 250 250 001;
453.1881.3290250 250 250 001;
561.0930.4550250 250 250 001;
671.0020.4170250 250 250 001;
784.4031.2150250 250 250 001;
895.6421.5970250 250 250 001;
9102.890.8180250 250 250 001;
10111.5140.4280250 250 250 001;
11121.2380.3510250 250 250 001;
];
return;
