summaryrefslogtreecommitdiff
path: root/1445/CH10/EX10.16/ch10_ex_16.sce
blob: f2fd3915019a9ff7f8af2cd8ab3fc1f86224d458 (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
41
42
43
44
45
46
47
48
49
50
//CHAPTER 10- THREE-PHASE INDUCTION MACHINES
//Example 16

disp("CHAPTER 10");
disp("EXAMPLE 16");

//VARIABLE INITIALIZATION
ratio1=1.5;                    //ratio of T_est and T_efl
ratio2=2.5;                    //ratio of T_em and T_efl

//SOLUTION
s=1;

//solution (a)
//directly solving the quadratic equation
a=1;
b=-3.333;
c=1;
D=(b)^2-(4*a*c);               //discriminant
sm1=(-b+sqrt(D))/(2*a);
sm2=(-b-sqrt(D))/(2*a);
if(sm1<=0 & sm2<=0) then
disp("The value of the slip at maximum torque is not valid");
else if(sm1>0 & sm1<1) 
disp(sprintf("The slip at maximum torque is %f",sm1));
else if(sm2>0 & sm2<1) 
disp(sprintf("The slip at maximum torque is %f",sm2));
end; 

//solution (b)
//directly solving the quadratic equation
a=1;
b=-1.665;
c=0.111;
D=(b)^2-(4*a*c);          
ans1=(-b+sqrt(D))/(2*a); 
ans2=(-b-sqrt(D))/(2*a);
if(ans1>0 & ans1<1)
disp(sprintf("The full load slip is %f",ans1));
sfl=ans1;
else if(ans2>0 & ans2<1)
disp(sprintf("The full load slip is %f",ans2));
sfl=ans2;
end;

//solution (c)
I=sqrt(ratio1/sfl);
disp(sprintf("The rotor current at the starting in terms of full load current is %f A",I));

//END