IEEE 14 bus Power flow code

Hi Anybody can send me the IEEE 14 bus power flow program.
Yunus

1 Comment

Chech the link. I complited it and also did fast decoupled power flow, Q limits cotrol by using Newton Raphson and using Crout's method.

Sign in to comment.

Answers (3)

function [baseMVA, bus, gen, branch, areas, gencost] = case14
%%system MVA base
baseMVA = 100;
%%bus data
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
bus = [
1 3 0 0 0 0 1 1.06 0 0 1 1.06 0.94;
2 2 21.7 12.7 0 0 1 1.045 -4.98 0 1 1.06 0.94;
3 2 94.2 19 0 0 1 1.01 -12.72 0 1 1.06 0.94;
4 1 47.8 -3.9 0 0 1 1.019 -10.33 0 1 1.06 0.94;
5 1 7.6 1.6 0 0 1 1.02 -8.78 0 1 1.06 0.94;
6 2 11.2 7.5 0 0 1 1.07 -14.22 0 1 1.06 0.94;
7 1 0 0 0 0 1 1.062 -13.37 0 1 1.06 0.94;
8 2 0 0 0 0 1 1.09 -13.36 0 1 1.06 0.94;
9 1 29.5 16.6 0 19 1 1.056 -14.94 0 1 1.06 0.94;
10 1 9 5.8 0 0 1 1.051 -15.1 0 1 1.06 0.94;
11 1 3.5 1.8 0 0 1 1.057 -14.79 0 1 1.06 0.94;
12 1 6.1 1.6 0 0 1 1.055 -15.07 0 1 1.06 0.94;
13 1 13.5 5.8 0 0 1 1.05 -15.16 0 1 1.06 0.94;
14 1 14.9 5 0 0 1 1.036 -16.04 0 1 1.06 0.94;
];
save baseMVA baseMVA;
save bus bus;
%%generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
gen = [
1 232.4 -16.9 10 0 1.06 100 1 332.4 0;
2 40 42.4 50 -40 1.045 100 1 140 0;
3 0 23.4 40 0 1.01 100 1 100 0;
6 0 12.2 24 -6 1.07 100 1 100 0;
8 0 17.4 24 -6 1.09 100 1 100 0;
];
save gen gen;
%%branch data
% fbus tbus r x b rateA rateB rateC ratio angle status
branch = [
1 2 0.01938 0.05917 0.0528 9900 0 0 0 0 1;
1 5 0.05403 0.22304 0.0492 9900 0 0 0 0 1;
2 3 0.04699 0.19797 0.0438 9900 0 0 0 0 1;
2 4 0.05811 0.17632 0.034 9900 0 0 0 0 1;
2 5 0.05695 0.17388 0.0346 9900 0 0 0 0 1;
3 4 0.06701 0.17103 0.0128 9900 0 0 0 0 1;
4 5 0.01335 0.04211 0 9900 0 0 0 0 1;
4 7 0 0.20912 0 9900 0 0 0.978 0 1;
4 9 0 0.55618 0 9900 0 0 0.969 0 1;
5 6 0 0.25202 0 9900 0 0 0.932 0 1;
6 11 0.09498 0.1989 0 9900 0 0 0 0 1;
6 12 0.12291 0.25581 0 9900 0 0 0 0 1;
6 13 0.06615 0.13027 0 9900 0 0 0 0 1;
7 8 0 0.17615 0 9900 0 0 0 0 1;
7 9 0 0.11001 0 9900 0 0 0 0 1;
9 10 0.03181 0.0845 0 9900 0 0 0 0 1;
9 14 0.12711 0.27038 0 9900 0 0 0 0 1;
10 11 0.08205 0.19207 0 9900 0 0 0 0 1;
12 13 0.22092 0.19988 0 9900 0 0 0 0 1;
13 14 0.17093 0.34802 0 9900 0 0 0 0 1;
];
save branch branch;
%%----- OPF Data -----%%
%%area data
areas = [
1 1;
];
%%generator cost data
% 1 startup shutdown n x1 y1 ... xn yn
% 2 startup shutdown n c(n-1) ... c0
gencost = [
2 0 0 3 0.0430293 20 0;
2 0 0 3 0.25 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
];
return;
this is the data file. you just try it for powerflow solution...

8 Comments

Thanks for Data file
I need also Power flow solution code. Please help me.
Was this code copied from the MATLAB File Exchange? If so then the power flow code is there.
plz send me the 14 bus system simulation model file in jasmindonga1996@yahoo.com
ezekiel gandham
ezekiel gandham on 2 Dec 2017
Edited: ezekiel gandham on 2 Dec 2017
plz send b coefficients for ieee14 bus 5units, ieee39 bus 10 units and ieee30 bus 6 unit systems to ezekielgandham123@gmail.com
plz send b coefficients for ieee14 bus 5units to nasefasmaa@yahoo.com
plz send the power flow matlab code to ecooljoe@gmail.com
https://www.mathworks.com/matlabcentral/answers/23781-ieee-14-bus-power-flow-code#comment_1771566

Sign in to comment.

POORNIMA JAYAPRAKASH
POORNIMA JAYAPRAKASH on 22 Jan 2014
Thanks for Data file I need also Power flow solution code. Please help me .....
Sonu Kumar Bairwa
Sonu Kumar Bairwa on 31 Jul 2021
any one can send me matlab code for optimal pmu placement by using integer linear programming for line outage and normal condition case. please help me. please send code on gmail enggsonubairwa@gmail.com

3 Comments

Yes, I recognize that as a project that gets assigned every couple of years, so there is someone in the world who has programmed it in the past. If you could find someone who had already programmed it, and you could convince them to share the code with you, you would be all set. It is possible.
However, in practice... it is probably going to take you longer to find someone who has already done the work and convince them to share the code with you, then it would take you to program it yourself.
Check out the link I send I completed it.
I used Newton Raphson to get it to converge and also used fast decoupled power flow analysis and Q limit control.
Let me know if you have questions.
Good luck!!

Sign in to comment.

Asked:

on 13 Dec 2011

Commented:

on 18 Oct 2021

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!