Dear All,
I am using Matpower for a distribution network power flow analysis.
The data I used are taken from a real network, and the network is working fine
using PowerWorld simulator. Basically, it is a radial network with 14 buses.
Bus 1 is the connection point to the main grid, which I modelled it as a slack
bus (type 3).
However, when converting the model into a Matpower case, I can only get a
failure notification "Did NOT converge (0.01 seconds)".
Does anyone have any idea about what might be the problem? Or what are the
options for me to debug?
I attached my case data in this email for your reference.
Thank you all very much.
Best regards,
Tian
*********************************************************************************************
function mpc = caseSG
%CASE SG Power flow data for a distribution network.
%% MATPOWER Case Format : Version 2
mpc.version = '2';
%%----- Power Flow Data -----%%
%% system MVA base
mpc.baseMVA = 1;
%% 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 0 22 1 1.1 0.9;
2 1 0 0 0 0 1 1 0 22 1 1.1 0.9;
3 1 0 0 0 0 1 1 0 22 1 1.1 0.9;
4 1 0 0 0 0 1 1 0 0.4 1 1.1 0.9;
5 1 0 0 0 0 1 1 0 0.4 1 1.1 0.9;
6 1 0.23 0.07 0 0 1 1 0 0.4 1 1.1 0.9;
7 1 0.2 0.06 0 0 1 1 0 0.4 1 1.1 0.9;
8 1 0.2 0.03 0 0 1 1 0 0.4 1 1.1 0.9;
9 1 0.13 0.06 0 0 1 1 0 0.4 1 1.1 0.9;
10 1 0.3 0.07 0 0 1 1 0 0.4 1 1.1 0.9;
11 1 0.2 0.05 0 0 1 1 0 0.4 1 1.1 0.9;
12 1 0.2 0.05 0 0 1 1 0 0.4 1 1.1 0.9;
13 1 0.17 0.05 0 0 1 1 0 0.4 1 1.1 0.9;
14 1 0.08 0.03 0 0 1 1 0 0.4 1 1.1 0.9;
];
%% 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 0 0 5 -5 1 1 1 10 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 2 0.000489 0.00064 0 5 5 5 0 0 1 -360 360;
1 3 0.000489 0.00064 0 5 5 5 0 0 1 -360 360;
2 4 0.094 1.34 0 5 5 5 55 0 1 -360 360;
3 5 0.094 1.34 0 5 5 5 55 0 1 -360 360;
4 6 5.876034 3.193497 0 1 1 1 0 0 1 -360 360;
4 7 5.208303 2.8306 0 1 1 1 0 0 1 -360 360;
4 8 6.343446 3.447525 0 1 1 1 0 0 1 -360 360;
4 9 8.179707 4.445493 0 1 1 1 0 0 1 -360 360;
4 10 7.678908 4.17332 0 1 1 1 0 0 1 -360 360;
5 11 6.677312 3.628974 0 1 1 1 0 0 1 -360 360;
5 12 6.343446 3.447525 0 1 1 1 0 0 1 -360 360;
5 13 3.639135 1.977791 0 1 1 1 0 0 1 -360 360;
5 14 5.876034 3.193497 0 1 1 1 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 2 1 0;
];
