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

[Geoff Hayes]

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]

did you get the code?

[NOOR habib]

Sir, i need the L-INDEX code too for my research work, thanks. [e-mail]

[zhibin zheng]

Sir, i need the L-INDEX code too for my research work, [e-mail]

[KAMALIKA BANERJEE]

I need the matlab codes to calculateLine voltage stability index

[KAMALIKA BANERJEE]

I also need L index and FVSI code for my research work

[sthiti Rout]

yes sir....can you please help me with L index code using IEEE 30 bus test system for my research work.... [e-mail]

[Kai Ling Hui]

Sir, i need the FVSI code too for my research work, thanks. [e-mail]

[Samson Ademola]

Please I need assistance on how to get FVSI on power flow analysis.this is my email [e-mail]

[Amulya Goyal]

need help in gvsm coding please help me

[Ezechukwu Ani]

Sir, Could you help me with FVSI code My email is [e-mail] I await your response, sir Thank you Nelson

[musa inanc]

Sir, i need the FVSI code too for my research work, thanks. [e-mail]

Answer 1

[pranay dattani]

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]

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]

You can put them in the vector

Answer 2

[Aneel Kumar]

How to calculate voltage deviation of any standard IEEE Bus system by Matlab coding

Answer 3

[Habib]

https://www.mathworks.com/matlabcentral/fileexchange/132088-fast-voltage-stability-index-fvsi

[Ezechukwu Ani]

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.

Answer 4

[Hafizur Rahman]

% FVSI Calculation from load flow

% Author: Hafizur Rahman

% Date: 7/08/2023

% E-mail: [e-mail]

% 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]

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.