Hello,
Please your help, this network does not converge and I do not know where is
the problem. Is a LV grid with a trafo at node 29 and a swing generator at
node 1, the rest are demand nodes (PQ) with cosfi=0.95.
Thanks in advance.
nodo 29 trafo Prim/sek. 10/0,4 kV Rated Power 400 kVA Connection dyn
5 (150ยบ) prim 10 kV Fuse 40A sec LV1: 160A
S_base 400000 VA
U_base 400 Voltios Z_base= 0,4 I_base= 577,3502692 A Y_base= 2,5
function mpc = OLV1
%% MATPOWER Case Format : Version 2
mpc.version = '2';
%%----- Power Flow Data -----%%
%% system MVA base
mpc.baseMVA = 0.4;
%% 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 0.4 1 1.05 0.95;
29 1 0 0 0 0 1 1 0 0.4 1 1.05 0.95;
10058 1 0.001094452 0.000359672 0 0 1 1 0 0.4 1 1.05
0.95;
10059 1 0.003676261 0.001208136 0 0 1 1 0 0.4 1 1.05
0.95;
10060 1 0.000513282 0.000168681 0 0 1 1 0 0.4 1 1.05
0.95;
10062 1 0.001844909 0.000606295 0 0 1 1 0 0.4 1 1.05
0.95;
10091 1 0.000303436 9.97187E-05 0 0 1 1 0 0.4 1 1.05
0.95;
10061 1 0.001741558 0.000572331 0 0 1 1 0 0.4 1 1.05
0.95;
10063 1 0.001753704 0.000576323 0 0 1 1 0 0.4 1 1.05
0.95;
4365 1 0.001925268 0.000632704 0 0 1 1 0 0.4 1 1.05
0.95;
10064 1 0.00077464 0.000254571 0 0 1 1 0 0.4 1 1.05
0.95;
10066 1 0.001225674 0.000402795 0 0 1 1 0 0.4 1 1.05
0.95;
4366 1 0.001187717 0.000390321 0 0 1 1 0 0.4 1 1.05
0.95;
10065 1 0.000540394 0.00017759 0 0 1 1 0 0.4 1 1.05
0.95;
10067 1 0.000420017 0.000138031 0 0 1 1 0 0.4 1 1.05
0.95;
4367 1 0.00199153 0.00065448 0 0 1 1 0 0.4 1 1.05
0.95;
10068 1 0.001361558 0.000447451 0 0 1 1 0 0.4 1 1.05
0.95;
4368 1 0.000353105 0.000116041 0 0 1 1 0 0.4 1 1.05
0.95;
4369 1 0.000708487 0.000232831 0 0 1 1 0 0.4 1 1.05
0.95;
];
%% 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 20 0.575 40 0 1 0.4 1 20 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 29 0.005 0.02 0 1 0 0 1 150 1 -360 360;
29 10058 0.0321 0.0078 0 0.284056332 0 0 0 0 1 -360 360;
10058 10059 0.04494 0.01092 0 0.284056332 0 0 0 0 1 -360
360;
10059 10060 0.020865 0.00507 0 0.284056332 0 0 0 0 1
-360 360;
10059 10062 0.0377175 0.009165 0 0.284056332 0 0 0 0 1
-360 360;
10060 10091 0.0794475 0.019305 0 0.284056332 0 0 0 0 1
-360 360;
10060 10061 0.043335 0.01053 0 0.284056332 0 0 0 0 1
-360 360;
10062 10063 0.0473475 0.011505 0 0.284056332 0 0 0 0 1
-360 360;
10061 4365 0.06 0.013125 0 0.239023011 0 0 0 0 1
-360 360;
10063 10064 0.0344 0.007525 0 0.239023011 0 0 0 0 1
-360 360;
10063 10066 0.033705 0.00819 0 0.284056332 0 0 0 0 1
-360 360;
4365 4366 0.0321 0.0078 0 0.284056332 0 0 0 0 1 -360
360;
10064 10065 0.040125 0.00975 0 0.284056332 0 0 0 0 1
-360 360;
10066 10067 0.01605 0.0039 0 0.284056332 0 0 0 0 1 -360
360;
4366 4367 0.092945 0.01044 0 0.162812776 0 0 0 0 1
-360 360;
10067 10068 0.036915 0.00897 0 0.284056332 0 0 0 0 1
-360 360;
10067 4368 0.234 0.01677 0 0.114315353 0 0 0 0 1 -360
360;
10067 4369 0.0963 0.0234 0 0.284056332 0 0 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 = [
1 0 0 0;
];
