Please see FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>.
— Ray > On Mar 9, 2015, at 8:21 AM, Xiaolin Ayon Flores <[email protected]> wrote: > > 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; > > ]; >
