infinite bus, 1 branch/1 gen

[Russ Patterson]

Hi,

 

I have a simple problem as shown below.  The upper bus is the infinite,
slack bus.  For the specified P & Q, I want to solve for VX.  Per my hand
calculations (iterations) the result will be VX=1.174pu @-3.3°.

 

Now, I am trying to model this in a simple load flow as my first ever
attempt to use MATPOWER (attached, which fails miserably).  

 



 

I have the manual and understand all of the bus, generator, and branch data.
What I don’t know is – what can be left out in solving this simple problem?
The machine is rated 16kV, 176MVA but I’d prefer to leave that out and just
work all in per unit (including data entry) if possible.

 

Thanks,

Russ

 

  _____  

Russell W. Patterson, P.E.

Chair – IEEE PSRC (http://pes-psrc.org) 

Office: 423-702-9981

 

 

 

 

 


function mpc = example1
%Example 1    
% Bus 1 is infinite and is 1.154 pu and zero degrees
% Bus 2 has a 176MVA generator, P=0.85pu and Q=0.279pu both flowing into bus 2
% Impedance between bus 1 & 2 is j0.0921 pu  
%
\%% MATPOWER Case Format : Version 2
mpc.version = '2';

\%%-----  Power Flow Data  -----\%%
\%% system MVA base
mpc.baseMVA = 176;

\%% bus data
%       bus_i   type    Pd      Qd      Gs      Bs      area    Vm      Va      
baseKV  zone    Vmax    Vmin
mpc.bus = [
  1 3 0 0 0 0 1 1.154 0 16 1 1.154 1.154;
  2 1 -.85 -.279 0 0 1 1 0 16 1 1.5 0.9;
];

\%% generator data
%       bus     Pg      Qg      Qmax    Qmin    Vg      mBase   status  Pmax    
Pmin    Pc1     Pc2     Qc1min  Qc1max  Qc2min  Qc2max  ramp_agc        ramp_10 
ramp_30 ramp_q  apf
mpc.gen = [
  2 149.6 49.104 49.104 49.104 1 176 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0; 
];

\%% branch data
%       fbus    tbus    r       x       b       rateA   rateB   rateC   ratio   
angle   status  angmin  angmax
mpc.branch = [
  1 2 0 .0921 0 0 0 0 0 0 1 0 0;
];