It may be because bus 12 (including its generator) is not connected to the
network, which you can see by using case_info() in MATPOWER 5.0b1.
Ray
On Aug 5, 2014, at 6:26 AM, Simone Cochi <simoneco...@gmail.com> wrote:
> Dear all,
> I'm getting hard times to solve a power flow on a real Italian medium
> voltage network. MATPOWER does not converge with any algorithm. I'm
> wondering if the branch data must have a certain order or my problem
> is simply the R/X ratio. The casefile is the following:
>
> function mpc = LFMT
> mpc.version = 2;
> mpc.baseMVA = 10;
> mpc.bus=[
> 1 3 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 2 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9;
> 3 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9;
> 4 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9;
> 5 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9;
> 6 1 0.05400 0.02615 0 0 1 1.0 0 20 1 1.1 0.9;
> 7 1 0.08640 0.04185 0 0 1 1.0 0 20 1 1.1 0.9;
> 8 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9;
> 9 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9;
> 10 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9;
> 11 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 12 2 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 13 2 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 14 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 15 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 16 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 17 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 18 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 19 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 20 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 21 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 22 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 23 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 24 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 25 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 26 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 27 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> 28 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9;
> ];
> mpc.gen=[
> 12 0.99440 0.00000 100 -100 1.0 100.000 1 1000.0 0.0 0 0 0 0 0 0 0 0 0 0 0;
> 13 0.99440 0.00000 100 -100 1.0 100.000 1 1000.0 0.0 0 0 0 0 0 0 0 0 0 0 0;
> 1 0.00000 0.00000 100 -100 1.0 100.000 1 1000.0 0.0 0 0 0 0 0 0 0 0 0 0 0;
> ];
> mpc.branch=[
> 19 27 0.01430 0.01060 0.01240 50.0 50.0 50.0 0 0 1 360 -360;
> 27 10 0.00220 0.00120 0.00120 50.0 50.0 50.0 0 0 1 360 -360;
> 27 15 0.00320 0.00240 0.00280 50.0 50.0 50.0 0 0 1 360 -360;
> 15 24 0.00840 0.00630 0.00740 50.0 50.0 50.0 0 0 1 360 -360;
> 24 16 0.00520 0.00390 0.00450 50.0 50.0 50.0 0 0 1 360 -360;
> 24 4 0.00040 0.00010 0.00820 50.0 50.0 50.0 0 0 1 360 -360;
> 16 25 0.00310 0.00230 0.00270 50.0 50.0 50.0 0 0 1 360 -360;
> 25 28 0.00170 0.00130 0.00150 50.0 50.0 50.0 0 0 1 360 -360;
> 25 5 0.00160 0.00090 0.00090 50.0 50.0 50.0 0 0 1 360 -360;
> 28 13 0.00400 0.00290 0.04580 50.0 50.0 50.0 0 0 1 360 -360;
> 28 14 0.00070 0.00020 0.01370 50.0 50.0 50.0 0 0 1 360 -360;
> 13 6 0.01320 0.00980 0.12280 50.0 50.0 50.0 0 0 1 360 -360;
> 6 26 0.00840 0.00630 0.00740 50.0 50.0 50.0 0 0 1 360 -360;
> 26 17 0.00130 0.00100 0.00110 50.0 50.0 50.0 0 0 1 360 -360;
> 26 8 0.00740 0.00420 0.00410 50.0 50.0 50.0 0 0 1 360 -360;
> 17 11 0.00900 0.00640 0.09160 50.0 50.0 50.0 0 0 1 360 -360;
> 11 18 0.00950 0.00710 0.05630 50.0 50.0 50.0 0 0 1 360 -360;
> 18 23 0.00060 0.00050 0.00060 50.0 50.0 50.0 0 0 1 360 -360;
> 23 3 0.02320 0.01320 0.01300 50.0 50.0 50.0 0 0 1 360 -360;
> 23 20 0.00520 0.00390 0.00450 50.0 50.0 50.0 0 0 1 360 -360;
> 20 21 0.01950 0.01450 0.01700 50.0 50.0 50.0 0 0 1 360 -360;
> 21 22 0.01490 0.01110 0.01300 50.0 50.0 50.0 0 0 1 360 -360;
> 22 7 0.00700 0.00520 0.06060 50.0 50.0 50.0 0 0 1 360 -360;
> 22 9 0.02480 0.01470 0.01490 50.0 50.0 50.0 0 0 1 360 -360;
> 7 2 0.01130 0.00650 0.00630 50.0 50.0 50.0 0 0 1 360 -360;
> 1 19 0.01970 0.01460 0.02750 50.0 50.0 50.0 0 0 1 360 -360;
>
> Thank you.
>
