If you want to model a load as a constant impedance, with the nominal MW 
occurring at 1 p.u. voltage, simply move the values from the PD and QD columns 
to the GS and BS columns and change the sign of the BS column.

In your example below, your GS and BS columns are still zero.

-- 
Ray Zimmerman
Senior Research Associate
419A Warren Hall, Cornell University, Ithaca, NY 14853
phone: (607) 255-9645




On Apr 5, 2013, at 1:56 PM, Cansu Yildirim <[email protected]> wrote:

> Dear Dr. Zimmerman, Dear all,
>  
> I want to model loads as shunt elements. How can I get the shunt conductance 
> Gs? I think I’m missing something when I use the equation
> Gs=(P/baseMVA)/(V^2)
> For the voltage I defined V=1 p.u.  Below you can see my example. I would be 
> grateful for any help. Thanks in advance!
>  
> Best regards,
>  
> Cansu Yildirim
>  
>  
>  
>  
> %% bus data
> %    bus_i    type    Pd    Qd    Gs    Bs    area    Vm    Va    baseKV    
> zone    Vmax    Vmin
> mpc.bus = [
>     1    1    0    0    0    0    1    1    0    345 1    1.1    0.9;
>     2    1    20    10    0   0    1    1    0    345    1    1.1    0.9;
>     3    1    0    0    0    0    1    1    0    345    1    1.1    0.9;
>     4    1    0    0    0    0    1    1    0    345    1    1.1    0.9;
>     5    1    40  15    0    0    1    1    0    345    1    1.1    0.9;
>     6    1    30    10    0    0    1    1    0    345    1    1.1    0.9;
>     7   1    0    0    0    0    1    1    0    345    1    1.1    0.9;
>     8   1    0    0    0    0    1    1    0    345    1    1.1    0.9;
>     9   1    0    0    0    0    1    1    0    345    1    1.1    0.9;
>     10  3    0    0    0    0    1    1    0    345    1    1.1    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 = [
>     10    0   0    1    1   1    100    1    1    1    0    0    0    0    0  
>   0    0    0    0    0    0;
>     7    10    5    1    1    1    100    1    1    1    0    0    0    0    
> 0    0    0    0    0    0    0;
>     8    30    20    1    1    1    100    1    1    1    0    0    0    0    
> 0    0    0    0    0    0    0;
>     9   20    10    1    1    1    100    1    1    1    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.05 0.2    0    0    0    0    0    0    1    -360    360;
>     2    3    0.1     0.5     0    0    0    0    0    0   1    -360    360;
>     3    4    0.2     0.8    0    0    0    0    0    0    0    -360    360;
>     4    5    0.1     0.3    0    0    0    0    0    0    1    -360    360;
>     5    6    0.2     0.4    0    0    0    0    0    0    1    -360    360;
>     6    1    0.1     0.15    0    0    0    0    0    0    1    -360    360;
>     2   5   0.2     0.5    0    0    0    0    0    0    1    -360    360;
>     10  1   0    0.004  0    0    0    0    0    0    1    -360    360;
>     7   2   0    1.56   0    0    0    0    0    0    1    -360    360;
>     8   3   0    0.5    0    0    0    0    0    0    1    -360    360;
>     9   4   0    0.4    0    0    0    0    0    0    1    -360    360;
> ];

Reply via email to