Voltage stability index code
显示 更早的评论
Hi, I am doing a project which require calculation of fast voltage stability index(FVSI). Please assist on how to come up with matlab code to calculate this index in power system
12 个评论
Geoff Hayes
2018-1-10
weldon - I removed your tag since it was identical to your second sentence. See A Guide to Tags for how tags should be used.
Also, please discuss what you have tried with respect to this problem.
SOMANAGOUDA AVATI
2019-3-13
did you get the code?
NOOR habib
2020-4-29
Sir, i need the L-INDEX code too for my research work, thanks. noorhabibkhan@ncepu.edu.cn
zhibin zheng
2021-1-4
Sir, i need the L-INDEX code too for my research work, thanks.892845589@qq.com
KAMALIKA BANERJEE
2021-1-23
I need the matlab codes to calculateLine voltage stability index
KAMALIKA BANERJEE
2021-1-23
I also need L index and FVSI code for my research work
sthiti Rout
2021-3-31
yes sir....can you please help me with L index code using IEEE 30 bus test system for my research work.... sthitirout814@gmail.com
Kai Ling Hui
2021-5-11
Sir, i need the FVSI code too for my research work, thanks. kailinghui97@gmail.com
Samson Ademola
2021-6-22
Please I need assistance on how to get FVSI on power flow analysis.this is my email Sademola092@gmail.com
Amulya Goyal
2022-9-8
need help in gvsm coding please help me
Ezechukwu Ani
2023-6-7
Sir, Could you help me with FVSI code My email is ezeani201@gmail.com I await your response, sir Thank you Nelson
musa inanc
2023-12-11
Sir, i need the FVSI code too for my research work, thanks. musainanc48@gmail.com
回答（4 个）
Hafizur Rahman
2023-8-30
1 个投票
% FVSI Calculation from load flow
% Author: Hafizur Rahman
% Date: 7/08/2023
% E-mail: 123ahmedhafiz@gmail.com
% N° From To R (pu) X (pu) 1/2B line line Limits
% line bus bus (pu) code
line_data=[ 1 1 2 0.01938 0.05917 0.0264 0 200
2 1 5 0.05403 0.22304 0.0246 0 100
3 2 3 0.04699 0.19797 0.0219 0 100
4 2 4 0.05811 0.17632 0.017 0 100
5 2 5 0.05695 0.17388 0.0173 0 100
6 3 4 0.06701 0.17103 0.0064 0 50
7 4 5 0.01335 0.04211 0 0 100
8 4 7 0 0.20912 0 0.978 50
9 4 9 0 0.55618 0 0.969 50
10 5 6 0 0.25202 0 0.932 100
11 6 11 0.09498 0.1989 0 0 50
12 6 12 0.12291 0.25581 0 0 20
13 6 13 0.06615 0.13027 0 0 50
14 7 8 0 0.17615 0 1 50
15 7 9 0 0.11001 0 1 50
16 9 10 0.03181 0.0845 0 0 20
17 9 14 0.12711 0.27038 0 0 20
18 10 11 0.08205 0.19207 0 0 20
19 12 13 0.22092 0.19988 0 0 20
20 13 14 0.17093 0.34802 0 0 20];
% N° bus Type Voltage Angle Load MW Load Mvar Gen MW Gen Mvar Gen Qmin Gen Qmax Cond. Suscep.
bus_data=[
1.0000 1.0000 1.0600 0 0 0 0 0 0 0 0 0
2.0000 2.0000 1.0450 0 21.7000 12.7000 40.0000 0 -40.0000 50.0000 0 0
3.0000 2.0000 1.0100 0 94.2000 19.0000 0 0 0 40.0000 0 0
4.0000 0 1.0000 0 47.8000 -3.9000 0 0 0 0 0 0
5.0000 0 1.0000 0 7.6000 1.6000 0 0 0 0 0 0
6.0000 2.0000 1.0700 0 11.2000 7.5000 0 0 -6.0000 24.0000 0 0
7.0000 0 1.0000 0 0 0 0 0 0 0 0 0
8.0000 2.0000 1.0900 0 0 0 0 0 -6.0000 24.0000 0 0
9.0000 0 1.0000 0 29.5000 16.6000 0 0 0 0 0 0.1900
10.0000 0 1.0000 0 9.0000 5.8000 0 0 0 0 0 0
11.0000 0 1.0000 0 3.5000 1.8000 0 0 0 0 0 0
12.0000 0 1.0000 0 6.1000 1.6000 0 0 0 0 0 0
13.0000 0 1.0000 0 13.5000 5.8000 0 0 0 0 0 0
14.0000 0 1.0000 0 14.9000 5.0000 0 0 0 0 0 0];
vb=100e3; % Base Voltage
sb=100e6; % Base VA
zb=(vb^2)/sb; % Base Impedance
% Add Line 2 to 3 Colome Data
data=[line_data(:,2:3), line_data(:,4:5)*zb];
% Calculate Impedance
data(:,3)=sqrt(data(:,3).^2+data(:,4).^2);
for i = 1:length(data)
% Add Voltage and Reactive power
data(i,5:6)=[100*bus_data(data(i,1),3) bus_data(data(i,2),6)]
end
disp(' From Bus To Bus Impedance Reactance Voltage Reactive Power');
disp(' (Z) (X) (KV) (MVAR)');
disp(data)
data(:,7)=(4*(data(:,3).^2).*data(:,6))./(data(:,4).*data(:,5).^2);
fvsi_val=data;
disp(' From Bus To Bus Impedance Reactance Voltage Reactive FVSI');
disp(' (Z) (X) (KV) Power(MVAR)');
disp(fvsi_val);
1 个评论
sujono cipto sutarno
2023-10-3
Apologies in advance.
I observed that impedance and reactance values are still in p.u, while voltage values are in kV, reactive power in MVAR.
Is that really the case? Thank you very much for sharing your knowledge.
pranay dattani
2021-3-7
0 个投票
i try to calculate the Line Stability Index, voltage Stability Index, linr stability factor
clc; clear all; close all; format compact;
Vs=0.9526+j*0.55; % Sending End Voltage in pu
Vsm=abs(Vs); % Magnitude of Sending End Voltage
Vsa=angle(Vs); % Angle of Sending End Voltage
Vr=1.00+j*0; % Receiving End Voltage in pu
Vrm=abs(Vr); % Magnitude of Receiving End Voltage
Vra=angle(Vr); % Angle of Receiving End Voltage
Delta=Vsa-Vra;
R=0.1; % Line Resistance
X=1; % Line Reactance
Z=R+j*X; % Line Impedance
Zm=abs(Z);
Theta=angle(Z);
Pr=1.2; % Receiving End Active Power
Qr=0.5; % Receiving End Reactive Power
Sr=Pr+j*Qr; % Receiving End Apprant Power
fprintf('Line Stability Index \n')
Lmn=(4*X*Qr)/((Vsm*sin(Theta-Delta))^2) % Line Stability Index
fprintf('Fast Voltage Stability Index \n')
FVSI=(4*(Zm^2)*Qr)/((Vsm^2)*X) % Fast Voltage Stability Index
fprintf('Line Stability Factor \n')
LQP=4*(X/(Vsm^2))*(Qr+(((Pr^2)*X)/(Vsm^2))) % Line Stability Factor
% Line is assumed to be lossless. So,Here the sending end power is equal
% To Reaceiving End Power
2 个评论
Samson Ademola
2021-8-15
Please how can I input multiple values of R and X to calculate and give answers once instead of changing the values when getting the value of FVS and Lmn. Thanks
oukennou aziz
2022-11-18
You can put them in the vector
Aneel Kumar
2023-1-6
0 个投票
How to calculate voltage deviation of any standard IEEE Bus system by Matlab coding
Habib
2023-7-7
0 个投票
https://www.mathworks.com/matlabcentral/fileexchange/132088-fast-voltage-stability-index-fvsi
1 个评论
Ezechukwu Ani
2023-7-8
Thank you sir, I have downloaded it,let me run it sir. I am also looking fuzzy logic code for integration of FACTS devices in 39 bus transmission system. I will be glad if you can help. Thank you once again.
类别
在 帮助中心 和 File Exchange 中查找有关 Downloads 的更多信息
2023-12-11
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
