Voltage stability index code

18 次查看（过去 30 天）
weldon
weldon 2018-1-10
评论： musa inanc 2023-12-11
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 个评论
显示 10更早的评论
Ezechukwu Ani
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
musa inanc 2023-12-11
Sir, i need the FVSI code too for my research work, thanks. musainanc48@gmail.com

请先登录，再进行评论。

请先登录，再回答此问题。

回答（4 个）

Hafizur Rahman
Hafizur Rahman 2023-8-30
% 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
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
pranay dattani 2021-3-7
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
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
oukennou aziz 2022-11-18
You can put them in the vector

请先登录，再进行评论。

Aneel Kumar
Aneel Kumar 2023-1-6
How to calculate voltage deviation of any standard IEEE Bus system by Matlab coding
Habib
Habib 2023-7-7
https://www.mathworks.com/matlabcentral/fileexchange/132088-fast-voltage-stability-index-fvsi
  1 个评论
Ezechukwu Ani
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.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 MATLAB 的更多信息

标签

尚未输入任何标签。

Community Treasure Hunt

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

Start Hunting!

Translated by