sir, how to connect SVC in the ckt with the matpower?
On Wed, Jul 13, 2016 at 10:27 AM, anshu dahiya <[email protected]> wrote: > sir how connect voltage regulator in matpower. is it connected in branch > data or in generator data? > > On Fri, Jul 1, 2016 at 10:12 PM, Ray Zimmerman <[email protected]> wrote: > >> mpopt = mpoption('linprog.Algorithm', 'dual-simplex'); >> >> Btw, please confirm all of the things I requested before asking any more >> questions about this topic. >> >> Ray >> >> >> On Jul 1, 2016, at 2:53 AM, anshu dahiya <[email protected]> wrote: >> >> sir i have tried very much to change solver but i couldn't understand >> how to change solver in matpower. >> >> >> On Thu, Jun 30, 2016 at 9:50 PM, Ray Zimmerman <[email protected]> wrote: >> >>> So can you confirm each of the follow? >>> - You say that the AC OPF converges but the DC OPF claims infeasibility, >>> correct? >>> - You’ve checked the system data … case_info() output looks good (no >>> islands, 1 reference bus, generation range adequate to serve total load, >>> etc.). >>> - Removing line limits (by setting RATE_A to 0) does not make that >>> problem feasible? >>> - You are using the Optimization Toolbox default solver (linprog or >>> quadprog with interior point solver). >>> >>> If so, can you try with a simplex solver (dual simplex is probably best) >>> and see if that gives a different result? >>> >>> Ray >>> >>> >>> On Jun 30, 2016, at 1:31 AM, anshu dahiya <[email protected]> >>> wrote: >>> >>> sir i have already relaxed my line limits. but it doesn't give any >>> response! >>> >>> >>> On Wed, Jun 29, 2016 at 6:29 PM, Ray Zimmerman <[email protected]> wrote: >>> >>>> It appears that your solver thinks the problem is infeasible. Try >>>> relaxing the branch flow limits to see if it solves. e.g. >>>> >>>> define_constants; >>>> mpc = loadcase(<your case file>); >>>> mpc.branch(:, RATE_A) = mpc.branch(:, RATE_A) * 2; % double the flow >>>> limits >>>> r = rundcopf(mpc); >>>> >>>> Or instead of doubling the limits, you could temporarily eliminate them >>>> completely. >>>> >>>> mpc.branch(:, RATE_A) = 0; % eliminates line flow limits completely >>>> >>>> >>>> Ray >>>> >>>> >>>> On Jun 29, 2016, at 5:19 AM, anshu dahiya <[email protected]> >>>> wrote: >>>> >>>> sir my acopf is running but not the dcopf what may be cause? The error >>>> comes in the matpower is as below: >>>> rundcopf(n84) >>>> >>>> MATPOWER Version 5.1, 20-Mar-2015 -- DC Optimal Power Flow >>>> The interior-point algorithm uses a built-in starting point; >>>> ignoring user-supplied X0. >>>> Exiting: One or more of the residuals, duality gap, or total relative >>>> error >>>> has grown 100000 times greater than its minimum value so far: >>>> the primal appears to be infeasible (and the dual unbounded). >>>> (The dual residual < TolFun=1.00e-08.) >>>> >>>> >>>>> Did NOT converge (0.06 seconds) <<<<< >>>> >>>> On Mon, Jun 27, 2016 at 5:44 PM, Saranya A <[email protected]> wrote: >>>> >>>>> Hi Anshu, >>>>> >>>>> Even I faced a similar problem with a different system. In >>>>> "mpc.branch" the 6th, 7th and 8th columns represent RATE A - MVA rating A >>>>> (long term rating) ,RATE B - MVA rating B (short term rating), RATE C - >>>>> MVA >>>>> rating C (emergency rating) respectively. Try giving a high MVA rating to >>>>> all the 3 columns and check if your solution converges. >>>>> >>>>> Regards, >>>>> Saranya.A >>>>> >>>>> On Mon, Jun 27, 2016 at 1:08 PM, anshu dahiya <[email protected] >>>>> > wrote: >>>>> >>>>>> i have double check my data it is giving the below error: >>>>>> *runopf(c84)* >>>>>> >>>>>> *MATPOWER Version 5.1, 20-Mar-2015 -- AC Optimal Power Flow* >>>>>> *MATLAB Interior Point Solver -- MIPS, Version 1.2, 20-Mar-2015* >>>>>> * (using built-in linear solver)* >>>>>> *Numerically Failed* >>>>>> >>>>>> *Did not converge in 6 iterations.* >>>>>> >>>>>> *>>>>> Did NOT converge (0.09 seconds) <<<<<* >>>>>> >>>>>> >>>>>> On Mon, Jun 27, 2016 at 11:51 AM, Fiaz Ahmad < >>>>>> [email protected]> wrote: >>>>>> >>>>>>> If there is no syntax error then check your branch matrix. It has >>>>>>> high R and X values,it will not converge. Probably you are not using the >>>>>>> correct per unit values. Your zbase=vbase*vbase/mvabase. Use zbase to >>>>>>> find >>>>>>> p.u values for R and X. >>>>>>> >>>>>>> I hope it helps >>>>>>> On Jun 27, 2016 8:53 AM, "anshu dahiya" <[email protected]> >>>>>>> wrote: >>>>>>> >>>>>>> sir my system is not converging while performing opf and dcopf , >>>>>>> what should i check in my program. Plz help me >>>>>>> >>>>>>> On Sun, Jun 26, 2016 at 5:53 AM, anshu dahiya < >>>>>>> [email protected]> wrote: >>>>>>> >>>>>>>> the problem is still as it was. my 84 bus system is still not >>>>>>>> converging after applying all steps directed by the Dr. Ray Zimmerman >>>>>>>> . can >>>>>>>> anyone help me to make the system converge. i am attaching my file >>>>>>>> >>>>>>>> On Fri, Jun 24, 2016 at 6:24 PM, Ray Zimmerman <[email protected]> >>>>>>>> wrote: >>>>>>>> >>>>>>>>> See FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>. >>>>>>>>> >>>>>>>>> Ray >>>>>>>>> >>>>>>>>> >>>>>>>>> On Jun 24, 2016, at 3:42 AM, anshu dahiya <[email protected]> >>>>>>>>> wrote: >>>>>>>>> >>>>>>>>> my 84 bus system is not converging into 10 iteration. what should >>>>>>>>> i take steps to make it converge. >>>>>>>>> >>>>>>>>> On Thu, Jun 23, 2016 at 8:43 PM, Ray Zimmerman <[email protected]> >>>>>>>>> wrote: >>>>>>>>> >>>>>>>>>> Whenever you have a problem with a case file, the first thing to >>>>>>>>>> do is always run case_info() >>>>>>>>>> <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0b1/case_info.html> >>>>>>>>>> to >>>>>>>>>> check for errors in the connectivity information … >>>>>>>>>> >>>>>>>>>> >> case_info('casem84') >>>>>>>>>> Bad bus numbers: 1 >>>>>>>>>> bus(1, BUS_I) = 0 >>>>>>>>>> >>>>>>>>>> Turns out that bus numbers must be positive integers, as stated >>>>>>>>>> in caseformat >>>>>>>>>> <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0b1/caseformat.html> >>>>>>>>>> . >>>>>>>>>> >>>>>>>>>> Ray >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> On Jun 23, 2016, at 1:59 AM, anshu dahiya < >>>>>>>>>> [email protected]> wrote: >>>>>>>>>> >>>>>>>>>> sir how to solve the following error in matpower which comes in >>>>>>>>>> the case84 which i have attached >>>>>>>>>> Error using sparse >>>>>>>>>> Index into matrix must be positive. >>>>>>>>>> >>>>>>>>>> Error in ext2int (line 141) >>>>>>>>>> n2i = sparse(mpc.bus(:, BUS_I), ones(nb, 1), 1:nb, >>>>>>>>>> max(mpc.bus(:, BUS_I)), 1); >>>>>>>>>> >>>>>>>>>> Error in runpf (line 111) >>>>>>>>>> mpc = ext2int(mpc); >>>>>>>>>> >>>>>>>>>> On Wed, Jun 22, 2016 at 6:59 PM, Ray Zimmerman <[email protected]> >>>>>>>>>> wrote: >>>>>>>>>> >>>>>>>>>>> Assuming this is for a power flow problem, you still have to >>>>>>>>>>> decide whether the bus is voltage controlled or not. If so, make it >>>>>>>>>>> a PV >>>>>>>>>>> bus. If not, and the DG simply injects a fixed amount of power >>>>>>>>>>> without >>>>>>>>>>> controlling voltage, then make it a PQ bus. >>>>>>>>>>> >>>>>>>>>>> Ray >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>>> > On Jun 22, 2016, at 1:27 AM, anshu dahiya < >>>>>>>>>>> [email protected]> wrote: >>>>>>>>>>> > >>>>>>>>>>> > how to nominate DG bus which has load as well as generation >>>>>>>>>>> > the nomenclature in the matpower as below >>>>>>>>>>> > 1= PQ bus >>>>>>>>>>> > 2= PV bus >>>>>>>>>>> > 3= reference bus >>>>>>>>>>> > 4= isolated bus >>>>>>>>>>> > i have attached my file plz take a look of it. suggest me >>>>>>>>>>> solution. >>>>>>>>>>> > <casem84.m> >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>> <casem84.m> >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>> >>>>>>> >>>>>> >>>>> >>>> >>>> >>> >>> >> >> >
