On Nov 4, 2013, at 1:05 PM, Kusi, Samuel A 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);

mpc = loadcase('case39');
mpc.gen(:,PG) = k*mpc.gen(:,PG);
mpc.bus(:,PD) = k*mpc.bus(:,PD);
mpc.bus(:,QD) = k*mpc.bus(:,QD);

mpc = runpf(mpc);

> 
> Also for a 39 bus system with 10 generators, shouldn't I expect 29 loads.

This is not necessary, any bus can have generator(s) and/or load(s) incident on 
it.

> The power flow simulation shows 21 loads. Any explanation this?

This is because only 21 buses out of the 39 buses have active/reactive power 
loads incident on them.

> 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?

I don't think there is a set metric for this but I presume that 50% drop in 
voltages would certainly create (or would be past) a voltage collapse scenario. 
However, you should really run a continuation power flow to create the voltage 
collapse scenario instead of merely increasing the loading level. MATPOWER has 
support for doing a continuation power flow via increasing the load at a given 
bus. The relevant files are available in the extras/cpf folder. 

Shri

Do read cpf_intro.pdf before doing your simulations.

> 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
> 
> 



Reply via email to