Hi Ray, thanks for your reply. Yes, I took exactly the data from that
test case. Solving the case without the dclines, you get an unfeasible
case, but matpower tries to solve the case using the opf, just that it
does not converge. Adding the dc lines, I get an error before
Matpower tries to solve the case. Something like this:
MATPOWER Version 4.1, 14-Dec-2011 -- AC Optimal Power Flow
MATLAB Interior Point Solver -- MIPS, Version 1.0, 07-Feb-2011Index
exceeds matrix dimensions.
Error in opf_hessfcn (line 116)
qcost = gencost(ng+1:2*ng, :);
Error in
mipsopf_solver>@(x,lambda,cost_mult)opf_hessfcn(x,lambda,cost_mult,om,Ybus,Yf(il,:),Yt(il,:),mpopt,il)
(line 144)
hess_fcn = @(x, lambda, cost_mult)opf_hessfcn(x, lambda, cost_mult, om, Ybus,
Yf(il,:), Yt(il,:), mpopt, il);
Error in mips (line 412)
Lxx = hess_fcn(x, lambda, opt.cost_mult);
Error in mipsopf_solver (line 145)
[x, f, info, Output, Lambda] = ...
Error in opf_execute (line 106)
[results, success, raw] = feval(solver, om, mpopt);
Error in opf (line 225)
[results, success, raw] = opf_execute(om, mpopt);
Error in runopf (line 96)
[r, success] = opf(casedata, mpopt);
Error in prueba (line 7)
results=runopf(mpc, opt);
>>
I hope you can help me.
Regards,
Santiago
2013/1/28 Ray Zimmerman <[email protected]>:
> Hi Santiago,
>
> I'm not sure why this doesn't work. It looks like the DC line parameters
> were taken directly from t_case9_dcline.m, which may not yield feasible
> flows for your system. My suggestion is to start with a working case and
> make changes (e.g. single DC link with minimal flow) one at a time to see
> what causes the breakdown.
>
> I'm not sure why, but it seems even without the DC lines your case is hard
> to solve.
>
> --
> Ray Zimmerman
> Senior Research Associate
> 419A Warren Hall, Cornell University, Ithaca, NY 14853
> phone: (607) 255-9645
>
>
>
>
> On Jan 27, 2013, at 3:35 PM, Santiago Torres <[email protected]>
> wrote:
>
> Dear Ray,
>
> I am trying to work on a 6 bus test system including DC links, however
> I am having problems, I dont exactly know where.
>
> I am including AC and DC lines in the same links; In another test
> system that comes with Matpower it is possible to work with DC and AC
> lines in the same links. But
> even if I work with especial links for DC lines only, my test case
> does not work.
>
> I am attaching my test case.
>
> Please take a look if possible.
>
> Regards,
>
> Santiago
>
> mpc=loadcase('t_garver_dcline');
> mpc=toggle_dcline(mpc, 'on');
>
> opt=mpoption('OPF_ALG', 560);
>
> results=runopf(mpc, opt);
>
>
> --
> Dr.-Ing. Santiago Torres
> IEEE Senior Member
>
> Post-Doctoral Fellow
> School of Electrical and Computer Engineering
>
>
> University of Campinas, Campinas, SP, Brazil
>
> http://www.dsee.fee.unicamp.br/
>
> Albert Einstein, 400
> 13083-852, Campinas, SP, Brazil
> <t_garver_dcline.m>
>
>
--
Dr.-Ing. Santiago Torres
IEEE Senior Member
Post-Doctoral Fellow
School of Electrical and Computer Engineering
University of Campinas, Campinas, SP, Brazil
http://www.dsee.fee.unicamp.br/
Albert Einstein, 400
13083-852, Campinas, SP, Brazil
