dear Dr Houssem make sure that R and X are in P.U. It'll works.
On Sat, Mar 29, 2014 at 10:30 AM, Bouchekara Houssem < [email protected]> wrote: > 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 type Pd Qd Gs Bs area Vm Va baseKV zone > Vmax Vmin > > bus = [ > > 1 3 0 0 0 0 1 1 0 13.8 1 1.1 0.9; > > 2 1 0.6 0.6 0 0 1 1 0 13.8 1 1.1 0.9; > > 3 2 0.4 0.3 0 0 1 1 0 13.8 1 1.1 0.9; > > 4 2 0.55 0.55 0 0 1 1 0 13.8 1 1.1 0.9; > > 5 2 0.3 0.3 0 0 1 1 0 13.8 1 1.1 0.9; > > 6 2 0.2 0.15 0 0 1 1 0 13.8 1 1.1 0.9; > > 7 2 0.55 0.55 0 0 1 1 0 13.8 1 1.1 0.9; > > 8 2 0.45 0.45 0 0 1 1 0 13.8 1 1.1 0.9; > > 9 2 0.4 0.4 0 0 1 1 0 13.8 1 1.1 0.9; > > 10 2 0.35 0.3 0 0 1 1 0 13.8 1 1.1 0.9; > > 11 2 0.4 0.3 0 0 1 1 0 13.8 1 1.1 0.9; > > 12 2 0.15 0.15 0 0 1 1 0 13.8 1 1.1 0.9; > > ]; > > > > %% generator data > > % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin > > gen = [ > > 12 5 0 .2 -.2 1 100 0 1.1 0; > > 11 5 0 .2 -.2 1 100 0 1.1 0; > > 10 5 0 .2 -.2 1 100 0 1.1 0; > > 9 5 0 .2 -.2 1 100 0 1.1 0; > > 8 5 0 .2 -.2 1 100 0 1.1 0; > > 7 5 0 .2 -.2 1 100 0 1.1 0; > > 6 5 0 .2 -.2 1 100 0 1.1 0; > > 5 5 0 .2 -.2 1 100 0 1.1 0; > > 4 5 0 .2 -.2 1 100 0 1.1 0; > > 3 5 0 .2 -.2 1 100 0 1.1 0; > > 1 0 0 100 -100 1 100 1 0 0; > > ]; > > > > %% branch data > > % fbus tbus r x b rateA rateB rateC ratio angle > status > > > > branch = [ > > 1 2 1.093 0.455 0 250 250 250 0 0 1; > > 2 3 1.184 0.494 0 250 250 250 0 0 1; > > 3 4 2.095 0.873 0 250 250 250 0 0 1; > > 4 5 3.188 1.329 0 250 250 250 0 0 1; > > 5 6 1.093 0.455 0 250 250 250 0 0 1; > > 6 7 1.002 0.417 0 250 250 250 0 0 1; > > 7 8 4.403 1.215 0 250 250 250 0 0 1; > > 8 9 5.642 1.597 0 250 250 250 0 0 1; > > 9 10 2.89 0.818 0 250 250 250 0 0 1; > > 10 11 1.514 0.428 0 250 250 250 0 0 1; > > 11 12 1.238 0.351 0 250 250 250 0 0 1; > > ]; > > > > return; > > >
