1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
|
clear;
clc;
function [Ecom,Emag]=vcurves(Prtd,Xs,Vt,Phi,K);
index=1;
if index==1;
X= input('percentage load');
Pg=X*Prtd;
Vt=Vt/sqrt(3);
pfang=Phi*%pi/180;
x=pfang:-.005:-pfang;
Pg=Pg*ones(1,length(x));
Ia=Pg./(3*abs(Vt)*cos(x));
Iacom=Ia.*(cos(x)+%i*sin(x));
Ecom=Vt+%i*Xs*Iacom;
Emag=abs(Ecom);
Ifld=Emag*1000/K;
plot(Ifld,Ia,'k');
set(gca(),"auto_clear","off");
return;
plot(Ifld,Ia,'k-');
return;
plot(Ifld,Ia,'k-.');
return;
plot(Ifld,Ia,'k.');
return;
set(gca(),"auto_clear","on");
else
end
endfunction
Prtd=36;
Vt=33;
Xs=10;
Phi=70;
K=1500;
[Ecom,Emag]=vcurves(Prtd,Xs,Vt,Phi,K);
[Ecom,Emag]=vcurves(Prtd,Xs,Vt,Phi,K);
[Ecom,Emag]=vcurves(Prtd,Xs,Vt,Phi,K);
[Ecom,Emag]=vcurves(Prtd,Xs,Vt,Phi,K);
index=0;
xlabel('Field current in amperes');
ylabel('Armature current in amp');
title('Plot of V-curves of a Synchronous machine');
|
