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