Thank you very much for your help!
I have another question. I want my opf to stop as soon as it finds an answer respecting the constraints (I have a constraint on the current in the branches). Thus, I don't want it to find the optimal solution but just the first one. Is it possible to do that? Thank you for your help, Gamze ---------------------------------------- > From: abhy...@mcs.anl.gov > To: matpowe...@list.cornell.edu > Subject: Re: OPF with current limit > Date: Wed, 16 Jul 2014 10:54:55 +0000 > > > >> On Jul 16, 2014, at 2:04 AM, "Gamze Dogan" <gamgam...@hotmail.com> wrote: >> >> Hi, I want to implement the OPF with current limit of the branches, I did >> using the option 'OPF_FLOW_LIM',2 >> > This will apply current limits on "all" branches except those whose flows > have been specified as infinite. > >> >> If I understand well, this will limit my current in all branches (exept one >> involving the slack bus) to 1 pu? >> > >> Is this right? >> >> How can I modify it to have a limit at 1.2 pu? >> > Set the branch RATE_A limit. >> >> Then I wanted to know if this was the right way to calculate my current: >> >> >> Sf(d,1) = (results.branch(d, 14)/results.baseMVA) + 1j * (results.branch(d, >> 15)/results.baseMVA); >> >> >> >> Vf(d,1) = results.bus(results.branch(d,1), 8) * exp(1j * >> results.bus(results.branch(d,1), 9)*pi/180); >> >> >> If(d,1) = conj( Sf(d,1) / Vf(d,1) ); > > [Ybus, Yf, Yt] = makeYbus(mpc); % mpc is your case file struct. > > If = Yf*V; % from bus currents > It = Yt*V; % to bus currents > >> >> >> >> Thank you for your help; >> >> >> >> Gamze >> > >