PSS-1 Experiments
PSS-1 Experiments
T-network π-network
Vs=(1+YZ1).Vr+(Z1+Z2+YZ1Z2).Ir Vs=(1+YrZ).Vr+Z.Ir
Is=Y.Vr+(1+YZ2). Ir Is=(Ys+Yr+ZYsYr).Vr+(1+YsZ).Ir
case 2,
A=1+(y*z)/2;
B=z;
C=y*(1+y*z/4);
D=1+(y*z)/2;
otherwise
disp('wrong method')
end
case 3,
length=input('enter length of long transmission line:')
zc=sqrt(z_km/y_km);
gam=(sqrt(z_km*y_km)*length);
A=cosh(gam);
B=zc*sinh(gam);
C=(1/zc)*sinh(gam);
D=A;
otherwise
disp('wrong choice of transmission line')
end
fprintf('\n the product of AD-BC=%f',(A*D)-(B*C)')
k3=input('\n enter 1- to read vr,ir and compute vs,is and 2- to read vs,is and compute vr,ir: ')
switch(k3)
case 1,
vr_line=input('enter receiving end voltage in kV(line-line):')
vr_phase=vr_line/sqrt(3);
ir=input('enter receiving end current in kA:')
vs_phase=A*vr_phase+B*ir;
is=C*vr_phase+D*ir;
fprintf('sending end voltage in kV\n')
vs_phase
fprintf('sending end current in kA\n')
is
case 2,
vs_line=input('enter sending end voltage in kV(line-line):')
vs_phase=vs_line/sqrt(3);
is=input('enter sending end current in kA:')
vr_phase=D*vs_phase-B*is;
ir=-C*vs_phase+D*is;
fprintf('receiving end voltage in kV\n')
vr_phase
fprintf('receiving end current in kA')
ir
otherwise
disp('wrong choice')
end
receiving_power=3*real(vr_phase*conj(ir));
sending_power=3*real(vs_phase*conj(is));
eff=(receiving_power/sending_power)*100;
reg=((abs(vs_phase)-abs(vr_phase))/abs(vr_phase))*100;
fprintf('\n receiving end power=%f MW',receiving_power')
fprintf('\n Sending end power=%f MW',sending_power')
fprintf('\n Efficiency=%f%%',eff')
fprintf('\n Regulation=%f%%\n',reg')
Output:
EXPERIMENT:6
Aim: Determination of power angle characteristics for salient and non-salient pole synchronous
machines, reluctance power, excitation emf and regulation.
Machine data:
Non salient pole machine Salient pole machine
Armature resistance/phase- 0.00 Ω Armature resistance/phase-0.025pu
Synchronous reactance/phase- 5Ω Direct axis reactance/phase-1.2pu
Quadrature axis reactance/phase-0.8pu
Terminal voltage/phase-25kV Terminal voltage/phase-1pu
Output power/phase-20MW Output power/phase-0.9pu
Power factor-0.8 (lag) Power factor-0.8 (lead)
Phasor diagrams:
MATLAB program for determination of power angle characteristics for non-salient pole
synchronous machine:
clc
clear all
xs=input('enter synchronous reactance/phase in ohm or in pu: \n')
ra=input('enter armature resistance/phase in ohm or in pu: \n')
v=input('enter the terminal voltage/phase in kV or in pu: \n')
p=input('enter the output power/phase in MW or in pu: \n')
pf=input('enter the power factor: \n')
pfsign=input('enter -1 for lagging pf and 1 for leading pf\n')
ia_mag=p/(v*pf);
theta=acos(pf)*pfsign;
real_ia=ia_mag*cos(theta);
imag_ia=ia_mag*sin(theta);
ia=complex(real_ia,imag_ia);
e=v+ia*(ra+j*xs);
pmax=abs(e)*abs(v/xs);
del=0:0.01:pi;
pe=pmax*sin(del);
plot(del*180/pi,pe,'r')
title('power angle characteristics of non salint pole synchronous machine')
xlabel('delta in degrees')
ylabel('Power in MW or in pu')
reg=(abs(e)-abs(v))/abs(v)*100;
fprintf('excitation emf in kV or in pu\n')
e
fprintf('magnitude of excitation emf in kV or in pu=%f\n',abs(e))
fprintf('maximum power in MW or in pu=%f\n',pmax)
fprintf('voltage regulation=%f%%\n',reg)
output: