Xcos examples from textbooks and for blocksHEADmaster

1 files changed, 50 insertions, 0 deletions

diff --git a/Working_Examples/293/CH18/EX18.1/eg18_1.sce b/Working_Examples/293/CH18/EX18.1/eg18_1.sce
new file mode 100755
index 0000000..068b46b
--- /dev/null
+++ b/Working_Examples/293/CH18/EX18.1/eg18_1.sce
@@ -0,0 +1,50 @@
+//a

+V1 = 440/sqrt(3);

+s = 0.025; //slip

+r1 = 0.1;

+r2 = 0.12;

+x1 = 0.35;

+x2 = 0.4;

+

+z = complex(r1 + r2/s, x1 + x2);

+i2 = V1/z; //input line current

+I2 = sqrt(real(i2)^2 + imag(i2)^2); //magnitude of input line current 

+disp("a")

+disp(i2,"input line current = ")

+

+i1 = complex(18*cos(-1.484), 18*sin(-1.484)); // magnetizing current

+I1 = sqrt(real(i1)^2 + imag(i1)^2); //magnitude of magnetizing current

+i = i1 + i2; //total current drawn from the voltage source

+I = sqrt(real(i)^2 + imag(i)^2); //magnitude of total current  

+theta = atan(imag(i)/real(i)); //phase difference between current and voltage 

+pf = cos(theta); //power factor

+disp(pf,"power factor = ")

+if theta >= 0 then

+    disp("leading")

+else disp("lagging")

+end

+

+//b

+f = 60; //hertz 

+ns = 1800; 

+ws = 2*%pi*ns/f; //stator angular velocity

+Pg = 3*I2^2*r2/s; //power 

+T = Pg/ws; //developed electromagnetic torque

+disp("b") 

+disp(T,"developed electromagneic torque (in Newton-meter) = ")

+

+//c

+Prot = 950; //rotational losses (in watts)

+Po = Pg*(1 - s) - Prot ; //output power

+HPo = Po/746; //output horse power

+disp("c")

+disp(HPo,"output horse power = ")

+

+//d

+Pc = 1200; //core losses (in W)

+SCL = 3*I^2*r1; //stator copper loss

+RCL = 3*I2^2*r2; //rotar copper loss

+loss = Pc + SCL + RCL + Prot; //total losses

+Pi = real(3*V1*i); //input power

+efficiency = 1 - (loss/Pi); 

+disp(efficiency,"efficiency = ")
\ No newline at end of file