diff options
Diffstat (limited to '3871/CH6/EX6.20')
| -rw-r--r-- | 3871/CH6/EX6.20/Ex6_20.sce | 37 |
1 files changed, 37 insertions, 0 deletions
diff --git a/3871/CH6/EX6.20/Ex6_20.sce b/3871/CH6/EX6.20/Ex6_20.sce new file mode 100644 index 000000000..241d21cbc --- /dev/null +++ b/3871/CH6/EX6.20/Ex6_20.sce @@ -0,0 +1,37 @@ +//===========================================================================
+//chapter 6 example 20
+
+clc;
+clear all;
+
+//variable declraration
+Ts = 300; //number of turns in secondary winding
+Tp = 3; //number of turns in primary winding
+Is = 5; //current in A
+Zs = (0.583)+%i*(0.25); //secondary impedance Ω
+n1 =10;
+n2 =5;
+
+//calculations
+KT =Ts/Tp; //turn ratio
+Es = Is*Zs; //secondary voltage in volts
+Nm = n1*Es; //total magnetising amp-turns
+Ni =n2*Es; //total iron loss amp-turns
+Im =Nm/Tp; //magnetising componenet of exciting current in A
+Ie = Ni/Tp; //
+I0 = Im+%i*Ie; //exciting current on primary side in A
+I01 =sqrt(((real(I0))^2)+((imag(I0))^2))
+alpha = atan(Ie/Im); //energy component of exciting current in A
+alpha1 = (alpha*180)/%pi
+theta = atan(imag(Zs)/real(Zs));
+theta1 = (theta*180)/%pi
+x = sin(((theta1+alpha1)*%pi)/180)
+Ip = (KT*Is)+(I01*x); //primary current in A
+y = cos(((theta1+alpha1)*%pi)/180);
+b =(180/%pi)*((I01*y)/(KT*Is)); //phase angle
+
+
+//result
+mprintf("primary current = %3f. A",y);
+mprintf("\nphase angle = %3.3f ",b);
+
|