1 files changed, 28 insertions, 0 deletions
diff --git a/3760/CH1/EX1.22/Ex1_22.sce b/3760/CH1/EX1.22/Ex1_22.sce
new file mode 100644
index 000000000..58977276d
--- /dev/null
+++ b/3760/CH1/EX1.22/Ex1_22.sce
@@ -0,0 +1,28 @@
+clc;
+P=10000; // rated power of transformer
+E1=2000; // primary side voltage
+E2=200; // secondary side voltage
+f=50; // frequency in hertz
+po=125; // no load power 
+pfo=0.15; // no load power factor
+zbh=E1^2/P; // base impedance on h v side
+k=E2/E1; // turns ratio
+zl=0.5+1*%i; // percent leakage impedance
+zlh=zl*(zbh*k^2); // percent leakage impedance referred to h v side
+Rc=E1^2/po; // coreloss resistance
+Io=po/(E1*pfo); // No load current
+Xm=E1/(Io*sqrt(1-pfo^2)); // magnetizing reactance
+p=10000; // load power
+pf=0.8; // power factor of load
+il=p/(E2*pf); // secondary load current
+ilp=il*k; // primary load current
+vp=E1+ilp*(real(zlh)*pf+imag(zlh)*sqrt(1-pf^2));
+ap=ilp^2*real(zlh); // active power loss in series resistance
+ar=ilp^2*imag(zlh); // reactive power loss in series reactance
+Ap=vp^2/Rc; // active power loss in coreloss resistance
+Ar=vp^2/Xm; // reactive power loss in magnetizing reactance
+cpl=p*(1+%i*tan(acos(0.8))); // complex power at load end in VA
+cpi=(real(cpl)+ap+Ap)+%i*(imag(cpl)+ar+Ar); // complex power input to transformer VA
+printf('real power input to transformer is %f W\n',real(cpi));
+ipf=cos(atan(imag(cpi),real(cpi)));
+printf('input power factor is %f lagging',ipf);