% POWER ANGLE CURVE FOR SYSNCHRONOUS MACHINE

cle;
clear all;

p-4s;

pl=0.8;

vt=34.64;

xd=9,33;

xq=9.33;

viph=ve* 1000/sqrt(3);
pfang= acos(pf);

q=p* tan(pfang);
i=(p=j*q)* 10*6/(3* vtph);
delta=0:1:180;
deltarad=delta*(pi/1 80);

if xd==xq

Farin tf(\n 2 RAAB E BER SAA A ER RIE IAAI IRE III RARE HI):
fprintf(\nThe details pertains to Non Salient machine’);
ef=vtph+(j*i*xd);

excitemf=abs(ef);

fprintf(\nExcitation_volt= %g V', excitemf);

regIn=(((excitemf)}-abs(vtph))/abs(vtph))* 100;

fprintf(\nRegulation= %g ', regin);

pwrnon=(excitemf)*vtph*sin(deltarad)/xd;

netpwr=3*pwrnon/10%6;

plot(delta,netpwr);

xlabel(‘\delta(deg)--->');

ylabel(‘Three ph power in Mw--->');

title¢POWER ANGLE CURVE FOR NON SALIENT MACHINE};

lepend(‘Non salient Machine’);

grid on

end

if xd~=xq

PR (I C\n 2 OA aA aA Ia AFIS IOSISIS EIEIO IIE):

{printi(\nThe details pertains to Salient machine’);

eq=vipht+(j*i*xq);

idmag~abs(i)*sin(angle(eq)-angle(i)):
CG Scanned with OKEN Scanner
efmag=abs(eq)+((xd-xq)*idmag);
excitemf=(efmag);

fprintf(\nExcitation_volt= %g V', excitemf);

regin=(((excitemf)-abs(vtph))/abs(vtph))* 100;
fprintf(\nRegulation= %g", regIn);
pp=excitem f* vtph*sin(deltarad)/xd;
relpwr=vtph*2*(xd-xq)*sin(2*deltarad)/(2*xd*xq);
netrelpwr=3*relpwr/10%6;

pwrsal=pptrelpwr;

netpwrsal=3*pwrsal/10*6;

plot(delta,netrelpwr,'r');

hold on

plot(delta,netpwrsal,'b');

xlabel(‘\delta(deg)--->');

ylabel('Three ph power in Mw--->');

title’-POWER ANGLE CURVE FOR SALIENT MACHINE);

legend('Reluctance power','Salient power');

grid on

end

CG Scanned with OKEN Scanner

6

mt

%Transient Stability Analysis of a Synchronous Machine with Critical Clearing Time

Calculation and Power-Angle Curve Plotting.

E = input('Enter internal voltage of machine ');

V = input(Enter infinite bus voltage :');

X1 = input(Enter reactance before fault :');

X2 = input(‘Enter reactance during fault :’),

X3 = input('Enter reactance after fault :’), .

Pm = input('Enter mechanical power input : ‘5

Pmaxl ={E * V/X1) / 1000; % Max power Pre fault (MW)

Pmax2 = (E * V/ X2) / 1000; % Max power During fault (MW)
Pmax3 = (E * V/ X3) / 1000; % Max power Post fault (MW)

dO =asind(Pm / Pmax!); % Power Angle (80)

% Critical Clearing Angle (Scr)

der = acosd((pi - 2 * deg2rad(d0)) * sind(d0) - cosd(d0));

% Convert to radians
d0_rad = deg2rad(d0),
der_rad = deg2rad(dcr);

% Critical Clearing Time Calculation

H = input(‘Enter Inertia Constant H (MJ/MVA) (default 205): ');

f= 50; % Frequency in Hz

ter = sqrt((2 * H * (der_rad - d0_rad)) / (pi *f*Pm)); % Critical clearing time

ter_ms = ter * 1000; % Convert seconds to msec

% Displaying Results

disp({'Pmax1 (Pre Fault) =', num2str(Pmax!),' MW'));

disp({‘Pmax2 (During Fault) =’, num2str(Pmax2), ' MW'});

disp({'Pmax3 (Post Fault) =', num2str(Pmax3), ' MW']);

disp([‘Initial Power Angle 50 =', num2str(d0), '° = (', num2str(d0_rad), ' rad)'});
disp((‘Critical Clearing Angle 6_cr =', num2str(der), '° = (', num2str(der_rad), '
rad)']):

disp([‘Critical Clearing Time t_er =", num2str(ter), ‘sce = (', num2str(ter_ms), '
ms)']):

\o.
CG Scanned with OKEN Scanner
% Stability Condition
angle_check = input(‘Enter the clearing angle §_cl (in degrees) for stability
check: ');
if angle_check <= der
stability_status = "Stables
else
stability_status = ‘Unstable’;

end

disp([‘System Stability: ', stability_status]);

% Plot Power-Angle Characteristics

d = 0:0.1:180; % Power angle 6 from 0 to 180 degrees

Pel = Pmax! * sind(d); % Pre-fault power curve

Pe2 = Pmax2 * sind(d); % During fault power curve

Pe3 = Pmax3 * sind(d); —_% Post-fault power curve

Pm_line = Pm * ones(size(d)); % Mechanical power input line

figure;

plot(d, Pel, '-r', 'LineWidth’, }.5); hold on;

plot(d, Pe2, '-k', ‘Line Width’, 1.5);

plot(d, Pe3, '-b', 'LineWidth’, 1.5);

plot(d, Pm_line, '--k', 'LineWidth’, 1.2); ’

xlabel('Power Angle (6 in degrees)’);

ylabel(‘Power (MW)');

title(‘Power-Angle Curve with Clearing Angle’);

yline(Pm, '--k', ‘Mechanical Power P_m');

xline(der, '--g', ["\delta_{cr} =', num2str(dcr), '° (' stability_status ')')); @

legend(‘Pre Fault ', ‘During Fault ', 'Post Fault ', ‘Pm’, ‘Location’, 'best');

grid on;

CG Scanned with OKEN Scanner

% Sequence component calculations

cle;

clear all;

Vr=input(\n Enter the Phase R volatage :');

Vy=input(‘\n Enter the Phase Y volatage :');

Vb=input(‘\n Enter the Phase B volatage :');

% Vr = 230; % Phase R voltage (V)

% Vy = 220; % Phase Y voltage (V)

% Vb = 210; % Phase B voltage (V)

voltages =[Vr; Vy; Vb);

% Calculate the sequence components

V0 = (1/3) * sum(voltages); % Zero sequence component

VI = (1/3) * (voltages(1) + voltages(2)*exp(1j*2*pi/3) +

voltages(3)*exp(-1j*2*pi/3)); % Positive sequence component

V2 = (1/3) * (voltages(1) + voltages(2)* ex p(-1j*2*pi/3) +

voltages(3)*exp(1j*2*pi/3)): % Negative sequence component

% Pack the results in a vector

sequence_components = [V0;V1; V2];

fprintf('Zero sequence component (V0): %.3f + j%.3f V\n', real(V0), imag(V0));

fprintf('Positive sequence component (V1): %.3f + j%.3f Vin’, real(V 1),
imag(V1));

fprintf('Negative sequence component (V2): %.3F + j%.3f Vin’, real(V2),

imag(V2));

CG Scanned with OKEN Scanner

% Sequence components

cle;

clear all;

Ir=input('\n Enter the Phase R current ');

ly=input(\n Enter the Phase Y current :’);

Tb=input(‘\n Enter the Phase B current :');

% Construct the phase current vector

% Ir = 10; % Phase R current (A)

% ly = 9; % Phase Y current (A)

% Tb = 8; % Phase B current (A)

currents = [Ir; ly; Ib];

% Calculate the sequence components

10 = (1/3) * sum(cwrrents); % Zero sequence component

I] = (1/3) * (currents(1) + currents(2)*exp(1j*2*pi/3) +

currents(3)*exp(-1j*2*pi/3)); % Positive sequence component @
12 = (1/3) * (currents(1) + currents(2)*exp(-1j*2*pi/3) +
currents(3)*exp(1j*2*pi/3)); % Negative sequence component

% Pack the results in a vector

sequence_components = [10; 11; 12];

fprintf(‘Zero sequence component (10): %.3f + j%.3f A\n', real(10), imag(10));

fprintf('Positive sequence component (11): %.3f+ j%.3F A\n’, real(14),
imag(I1));

fprintf(‘Negative sequence component (12): %.3f + j%.3f A\n’, real(I2),

imag(I2));

CG Scanned with OKEN Scanner

% ABCD parameters of MEDIUM TRANSMISSION LINE

cle;

clear;

R=10;

X= 20;
Bc=4*10%-4;
Vr3ph=66;
Pr=10;
Qr=7.5;
Z=R+j*X;
Y=0+j*Be;

Piper ira tf 9 8 9 a IIR TY

fprintf(\nEnter |- to choose Pie Network\n');

fprintf(‘\nEnter 2- to choose T Network\n’);

Fprint l(* * #8 He RR EE EE EEE He EE EN)

Choice=input(‘\nEnter the type of network modelling : ‘);

switch Choice
case |

fprintf(‘\n You have chosen Pie type network\n');

A=1+Z*Y/2;

B=Z,;

C= Y*(14Z*Y/4);
D= 14+Z*Y/2;
ABCD=[A B; C D];

case 2

fprintf(‘\n You have chosen T type network\n’);

A=1+Z*Y/2;

B= Z*(1+Z*Y/4);

C=Y ;

D=1+Z*Y?2;

ABCD=[A B; C DJ];

otherwise

error('Your Choice is improper’)

end

fprinth(\n* ** #88 EER RRA AAA HE RE RE HN

disp(ABCD);

Tprinth(\n* A reek eee Hr

fprintf(‘\n The product AD-BC=%f,A*D-B*C);

Cpr iC \ra# #8 A eH ee eo ee a #1)

Vr=Vr3ph/sqrt(3);
Sr=Pr+j*Qr;
Ir-conj(Sr)/(3*conj(Vr));
Vs=VslIs(1):

Is=Vsls(2);
Vs3ph=sqrt(3)*abs(Vs);
Ism=1000* abs(Is);

.
+

>

CG Scanned with OKEN Scanner

plfs=cos(angle(Vs)-angle(Is));

Ss=3* Vs*conj(Is);

Reg=(Vs3ph-Vr3ph)/Vr3 ph* 100;

Effy=Pr/real(Ss)* 100;

fprintt(\n Is=%g A',Ism);

fprintf(\n Vs=%g Kv',Vs3ph);

fprintf(\n Ps=%g Mw’, real(Ss));

fprintf(\n Qs=%g Mvar’, imag(Ss));

fprintf(\n pfs=%g ', pfs);

fprintf(\n Percentage voltage regulation =%g' Reg);

fprintf(‘\n Efficiency of Transmission line=%g' Effy);

fprintf(\n* 0 OCC HBO EO a OEE oD OEE).

CG Scanned with OKEN Scanner