Dear Forum I am trying to run a simple power flow in the attached file but failed to converge, couldn't find what actually went wrong. If anyone could help me sort the problem out that would be really great. RegardsKawsar
function mpc = Test1 % Capstone B % Project simulation % M M Kawsar Farooki % SID: 450594762 % Using default voltage magnitude limits: VMIN = 0.9 p.u., VMAX = 1.1 p.u.
%% MATPOWER Case Format : Version 2
mpc.version = '2';
%%----- Power Flow Data -----%%
%% system MVA base
mpc.baseMVA = 100;
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va
baseKV zone Vmax Vmin
mpc.bus = [
101 1 80 22 0 0 1 1 0.8
100 1 1.1 0.9;
102 1 93 28 0 0 1 1.03505 0.6109
100 1 1.1 0.9;
103 1 90 28 0 0 1 1 0.8891
100 1 1.1 0.9;
104 1 30 8 0 0 1 1 6.5272
100 1 1.1 0.9;
105 1 63 24 0 0 1 1 1.1932
100 1 1.1 0.9;
106 2 0 0 0 0 1 1 6
100 1 1.1 0.9; % PV
107 1 105 32 0 0 1 1.045 1.7006
100 1 1.1 0.9;
108 1 34 11 0 0 1 1.0157 6.4416
100 1 1.1 0.9;
109 1 50 16 0 0 1 1.00308 -7.4573
100 1 1.1 0.9;
201 1 5 1 0 0 2 1.00744 0.0143
100 2 1.1 0.9;
202 1 5 1 0 0 2 1.00944 1.0141
100 2 1.1 0.9;
203 1 8 3 0 0 2 1.04585 1.2698
100 2 1.1 0.9;
204 1 8 2 0 0 2 1.01444 12.533 100
2 1.1 0.9;
205 1 8 1 0 0 2 1.02324 2.8034 100
2 1.1 0.9;
206 2 0 0 0 0 2 1.04484 10.6555
100 2 1.1 0.9; % PV
207 1 4 2 0 0 2 1.03741 2.901
100 2 1.1 0.9;
208 1 6 3 0 0 2 1.01499 12.4176
100 2 1.1 0.9;
209 1 4 1 0 0 2 1.03448 1.6753 100
2 1.1 0.9;
210 1 7 4 0 0 2 1.03223 17.4582
100 2 1.1 0.9;
211 1 20 6 0 0 2 1 -6.3752
100 2 1.1 0.9;
301 2 0 0 0 0 3 0.9 -1
100 3 1.1 0.9; % PV
302 1 250 72 0 0 3 0.900488
-13.9924 100 3 1.1 0.9;
303 1 250 93 0 0 3 0.90542
-22.8954 100 3 1.1 0.9;
304 1 80 28 0 0 3 1.1 3.0067
100 3 1.1 0.9;
305 1 120 32 0 0 3 0.90183
-25.5714 100 3 1.1 0.9;
306 1 90 22 0 0 3 0.902074
-25.6045 100 3 1.1 0.9;
307 2 0 0 0 0 3 1.1 3.9496
100 3 1.1 0.9; % PV
308 2 60 20 0 0 3 0.904184
-23.2471 100 3 1.1 0.9; % PV
309 1 40 15 0 0 3 0.900776
-24.165 100 3 1.1 0.9;
310 1 30 8 0 0 3 0.903746
-23.6261 100 3 1.1 0.9;
311 1 25 12 0 0 3 0.901132
-22.6716 100 3 1.1 0.9;
312 1 192 44 0 0 3 0.9015 -28.3633
100 3 1.1 0.9;
313 1 55 16 0 0 3 0.901615
-27.4327 100 3 1.1 0.9;
314 1 58 22 0 0 3 1.0904266
-34.3703 100 3 1.1 0.9;
315 1 20 10 0 0 3 0.8
-13.1037 100 3 1.1 0.9;
316 1 15 5 0 0 3 0.7 -2.2181
100 3 1.1 0.9;
317 1 90 30 0 0 3 0.6 -1.9537
100 3 1.1 0.9;
318 1 65 18 0 0 3 0.7 -4.726
100 3 1.1 0.9;
401 3 0 0 0 0 4 1 0
100 4 1.1 0.9; % Slack
];
%% 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 mu_Pmax mu_Pmin mu_Qmax mu_Qmin
mpc.gen = [
104 100 46 9999 -9999 1 100 1 9999
-9999 0 0 0 0 0 0 0 0 0
0 0 0.0000 0.0000 0.0000 0.0000;
106 300 93 9999 -9999 1 100 1 9999
-9999 0 0 0 0 0 0 0 0 0
0 0 0.0000 0.0000 0.0000 0.0000;
108 200 47 9999 -9999 1 100 1 9999
-9999 0 0 0 0 0 0 0 0 0
0 0 0.0000 0.0000 0.0000 0.0000;
204 200 70 9999 -9999 1 100 1 9999 -9999
0 0 0 0 0 0 0 0 0 0
0 0.0000 0.0000 0.0000 0.0000;
206 400 110 9999 -9999 1 100 1 9999
-9999 0 0 0 0 0 0 0 0 0
0 0 0.0000 0.0000 0.0000 0.0000;
208 350 122 9999 -9999 1 100 1 9999
-9999 0 0 0 0 0 0 0 0 0
0 0 0.0000 0.0000 0.0000 0.0000;
210 130 32 9999 -9999 1 100 1 9999
-9999 0 0 0 0 0 0 0 0 0
0 0 0.0000 0.0000 0.0000 0.0000;
304 120 37 9999 -9999 1 100 1 9999
-9999 0 0 0 0 0 0 0 0 0
0 0 0.0000 0.0000 0.0000 0.0000;
307 200 62 9999 -9999 1 100 1 9999
-9999 0 0 0 0 0 0 0 0 0
0 0 0.0000 0.0000 0.0000 0.0000;
314 160 49 9999 -9999 1 100 1 9999
-9999 0 0 0 0 0 0 0 0 0
0 0 0.0000 0.0000 0.0000 0.0000;
401 200 600 9999 -9999 1 100 1 9999 -9999
0 0 0 0 0 0 0 0 0 0
0 0.0000 0.0000 0.0000 0.0000;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio
angle status angmin angmax
mpc.branch = [
101 102 0.0084 0.6 0.0817 0 0 0 0
0 1 -360 360;
101 109 0.0084 0.6 0.0817 0 0 0 0
0 1 -360 360;
101 401 0.0078 0.62 0.076 0 0 0 0
0 1 -360 360;
102 103 0.0078 0.62 0.076 0 0 0 0
0 1 -360 360;
103 104 0.0045 0.3 0.0437 0 0 0 0
0 1 -360 360;
104 105 0.0045 0.3056 0.0437 0 0 0 0
0 1 -360 360;
105 106 0.00109 0.8068 0.076 0 0 0 0
0 1 -360 360;
106 107 0.0096 0.706 0.0931 0 0 0 0
0 1 -360 360;
107 108 0.0096 0.706 0.0931 0 0 0 0
0 1 -360 360;
108 109 0.0037 0.406 0.073 0 0 0 0
0 1 -360 360;
106 315 0.0037 0.406 0.073 0 0 0 0
0 1 -360 360;
106 206 0.0045 0.3056 0.0437 0 0 0 0
0 1 -360 360;
201 401 0.0045 0.3056 0.0437 0 0 0 0
0 1 -360 360;
201 202 0.0066 0.5027 0.0646 0 0 0 0
0 1 -360 360;
202 203 0.0066 0.5027 0.0646 0 0 0 0
0 1 -360 360;
203 204 0.0066 0.5027 0.0646 0 0 0 0
0 1 -360 360;
204 205 0.0066 0.5027 0.0646 0 0 0 0
0 1 -360 360;
205 206 0.0018 0.104 0.171 0 0 0 0
0 1 -360 360;
206 207 0.0018 0.104 0.171 0 0 0 0
0 1 -360 360;
207 208 0.0008 0.62 0.076 0 0 0 0
0 1 -360 360;
208 209 0.0031 0.2048 0.304 0 0 0 0
0 1 -360 360;
209 210 0.0031 0.2048 0.304 0 0 0 0
0 1 -360 360;
210 211 0.0031 0.2048 0.304 0 0 0 0
0 1 -360 360;
211 201 0.0045 0.3056 0.437 0 0 0 0
0 1 -360 360;
207 307 0.001 0.1405 0.54 0 0 0 0
0 1 -360 360;
301 401 0.001 0.1405 0.54 0 0 0 0
0 1 -360 360;
301 302 0.0014 0.108 0.133 0 0 0 0
0 1 -360 360;
302 303 0.0014 0.108 0.133 0 0 0 0
0 1 -360 360;
303 304 0.0019 0.155 0.19 0 0 0 0
0 1 -360 360;
304 305 0.007 0.5058 0.684 0 0 0 0
0 1 -360 360;
305 306 0.001 0.77 0.095 0 0 0 0
0 1 -360 360;
306 307 0.0049 0.3088 0.475 0 0 0 0
0 1 -360 360;
307 308 0.0051 0.403 0.494 0 0 0 0
0 1 -360 360;
308 309 0.0051 0.403 0.494 0 0 0 0
0 1 -360 360;
309 310 0.0072 0.5074 0.703 0 0 0 0
0 1 -360 360;
310 311 0.0072 0.5074 0.703 0 0 0 0
0 1 -360 360;
311 312 0.0051 0.403 0.494 0 0 0 0
0 1 -360 360;
312 301 0.001 0.4 0.48 0 0 0 0
0 1 -360 360;
301 318 0.0011 0.16 0.7 0 0 0 0
0 1 -360 360;
301 313 0.0011 0.16 0.7 0 0 0 0
0 1 -360 360;
313 314 0.0003 0.4 0.424 0 0 0 0
0 1 -360 360;
314 315 0.0002 0.3 0.32 0 0 0 0
0 1 -360 360;
314 317 0.0003 0.45 0.447 0 0 0 0
0 1 -360 360;
317 316 0.0003 0.45 0.447 0 0 0 0
0 1 -360 360;
];
%% bus names
mpc.bus_name = {
'1_101';
'1_102';
'1_103';
'1_104';
'1_105';
'1_106';
'1_107';
'1_108';
'1_109';
'2_201';
'2_202';
'2_203';
'2_204';
'2_205';
'2_206';
'2_207';
'2_208';
'2_209';
'2_210';
'2_211';
'3_301';
'3_302';
'3_303';
'3_304';
'3_305';
'3_306';
'3_307';
'3_308';
'3_309';
'3_310';
'3_311';
'3_312';
'3_313';
'3_314';
'3_315';
'3_316';
'3_317';
'3_318';
'4_401'
};
