summaryrefslogtreecommitdiff
path: root/503/CH8/EX8.14/ch8_14.sci
blob: e28b01d8c126cfd5e6af76dc36ae5649521e056f (plain)
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
//to calculate armature current, pf ,power angle, power , shaft torques,kVar

clc;
j=sqrt(-1);
P=8000;
Prot=500;
Pmg=P+Prot;
Pein=Pmg;
Ef=750/sqrt(3);
Vt=231;
Xs=16;
dl=asind(Xs*(Pein/3)/(Ef*Vt));
Eff=Ef*complex(cosd(-dl),sind(-dl));
Ia=(Vt-Eff)/(j*Xs);
disp(abs(Ia),'armature current(A)');
disp(cosd(atand(imag(Ia)/real(Ia))),'pf');
f=50;
p=4;
n_s=120*f/p;
w_s=2*%pi*n_s/60;
T=Pein/w_s;disp(T,'torque developed(Nm)');
T_s=P/w_s;disp(T_s,'shaft torques(Nm)');

Ef=600/sqrt(3);
Ia=(Vt-Ef)/(j*Xs);
rr=3*Vt*Ia/1000;
disp(rr,'kVar rating');
c=(abs(Ia)/Vt)/(2*%pi*f);
disp(-c,'capicator rating(F)');

Ef=300/sqrt(3);
Ia=(Vt-Ef)/(j*Xs);
rr=3*Vt*Ia/1000;
disp(-rr,'kVar rating');
L=(Vt/abs(rr))/(2*%pi*f);
disp(L,'inductor rating(H)');

Ia=j*2000/Vt;
Ef=Vt-j*Ia*Xs;
disp(abs(Ef)*sqrt(3),'excitation(V)');