I haven’t looked at the case directly, but the first thing to do is check the binding constraints. If setting branch(:, RATE_A) to zero removes the issue, then it’s the branch limits.
Ray > On Aug 9, 2018, at 12:16 AM, Aamir Nawaz <[email protected]> wrote: > > Hi, > > In case given below, expected dispatch cannot able to go above 70% of each > generator. For example, max achievable generation is 900MW but this example > can achieve maximum of 700MW, it cant support load greater than that. Can > anyone explain why? How can i increase load more than above mentioned limit? > > casefile = 'ex_case3b'; > mpc = loadcase(casefile); > xgd = loadxgendata('ex_xgd_ramp', mpc); > profiles = getprofiles('ex_load_profile'); > nt = size(profiles(1).values, 1); % number of periods > mdi = loadmd(mpc, nt, xgd, [], [], profiles); > mpopt = mpoption('verbose', 0, 'most.dc_model', 1); > mdo = most(mdi, mpopt); > EPg = mdo.results.ExpectedDispatch > Elam = mdo.results.GenPrices > > > Regards, > Aamir Nawaz, > PhD student >
