Note that what you want to do would give you a sub-optimal solution. You could try using a larger objective function cost tolerance to have the optimization terminate early. This can be set via mpoption using the option mips.costtol (if you are using MIPS, see the corresponding option for other solvers) in MATPOWER 5.0b1.
Shri -----Original Message----- From: Gamze Dogan <[email protected]> Reply-To: MATPOWER discussion forum <[email protected]> Date: Tue, 22 Jul 2014 12:40:13 +0000 To: MATPOWER discussion forum <[email protected]> Subject: RE: OPF with current limit >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: [email protected] >> To: [email protected] >> 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" <[email protected]> >>>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 >>> >> >> >
