To increase the loads and generators by a factor you could do something like …
define_constants; k = 1.1; mpc = loadcase(‘case39’); mpc.bus = scale_load(k, mpc.bus); mpc.gen(:, PG) = k * mpc.gen(:, PG); Each bus can have a load, a generator or both, so there is no reason to expect 29 loads. Regarding voltage collapse, I suggest you find a good textbook discussion of voltage collapse before attempting to model it. You’ll probably want to use a continuation power flow (see the extras/cpf directory) to find the point of voltage collapse. -- Ray Zimmerman Senior Research Associate B30 Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645 On Nov 4, 2013, at 2:05 PM, Kusi, Samuel A <[email protected]> wrote: > Thank you Ray! It worked! > > I am also trying to increase the load and generators powers by a factor. This > is what I have so far. Obviously, its not working. Any comments? > > k==1.1 > define_constants; > mpc = loadcase('case39'); > newgenpwr = k* results.gen(:,PG); > % not sure about syntax to get load info > newloadpwr = k* result....(:,PD); > > > Also for a 39 bus system with 10 generators, shouldn't I expect 29 loads. The > power flow simulation shows 21 loads. Any explanation this? > Also if I want to model a scenario for voltage collapse, at what voltage > level at the generators will you describe as a collapse? 50% drop? > Please advise. Thank you > > > > ----- Original Message ----- > From: "Ray Zimmerman" <[email protected]> > To: "MATPOWER discussion forum" <[email protected]> > Sent: Monday, November 4, 2013 4:23:39 PM > Subject: Re: fname -Name of file to which print output is to be appended > > I think what you want is something like … > > mpopt = mpoption; > results = runpf('case39', mpopt, ‘flowoutput.txt'); > > The file will be a simple text file with the pretty-printed output in it. > > -- > Ray Zimmerman > Senior Research Associate > B30 Warren Hall, Cornell University, Ithaca, NY 14853 > phone: (607) 255-9645 > > > > On Nov 3, 2013, at 9:56 AM, Kusi, Samuel A <[email protected]> wrote: > >> >> How do I define a file name for print output to be appended? What format >> will this file be in? >> This is the code I am using: I keep getting error Undefined function or >> variable 'flowoutput' >> >> mpopt = mpoption; >> fname = flowoutput; >> results = runpf('case39', mpopt, flowoutput); >> >> Please advise. >> >> -- >> Samuel Kusi >> Electrical Engineering >> GTID#: 902765579 >> >> >> > > > -- > Samuel Kusi > Electrical Engineering > GTID#: 902765579 > > > >
