# Voltage stability index code

45 views (last 30 days)
weldon on 10 Jan 2018
Commented: musa inanc on 11 Dec 2023
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
Ezechukwu Ani on 7 Jun 2023
Sir, Could you help me with FVSI code My email is ezeani201@gmail.com I await your response, sir Thank you Nelson
musa inanc on 11 Dec 2023
Sir, i need the FVSI code too for my research work, thanks. musainanc48@gmail.com

pranay dattani on 7 Mar 2021
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
Samson Ademola on 15 Aug 2021
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 on 18 Nov 2022
You can put them in the vector

Aneel Kumar on 6 Jan 2023
How to calculate voltage deviation of any standard IEEE Bus system by Matlab coding

Habib on 7 Jul 2023
https://www.mathworks.com/matlabcentral/fileexchange/132088-fast-voltage-stability-index-fvsi
Ezechukwu Ani on 8 Jul 2023
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.

Hafizur Rahman on 30 Aug 2023
% 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);
sujono cipto sutarno on 3 Oct 2023
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.

### Categories

Find more on Power Transmission and Distribution in Help Center and File Exchange

### Tags

No tags entered yet.

### Community Treasure Hunt

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

Start Hunting!