1 files changed, 30 insertions, 0 deletions
diff --git a/3760/CH1/EX1.16/Ex1_16.sce b/3760/CH1/EX1.16/Ex1_16.sce
new file mode 100644
index 000000000..f7e463735
--- /dev/null
+++ b/3760/CH1/EX1.16/Ex1_16.sce
@@ -0,0 +1,30 @@
+clc;
+P=200000; //rated power of transformer
+E1=4000; // primary side rated voltage 
+E2=1000; // secondary side rated voltage
+n=0.97; // efficiency
+pfn=0.25; // power factor at no load
+pff=0.8; // power factor at full load
+vr=5; // percentage voltage regulation
+Pl=((1/n)-1)*200000; // total losses at full load
+Pf=Pl*0.6; // total losses at 60% of full load
+Po=(Pl-Pf)/(1-0.36); // ohmic losses
+Pc=Pl-Po; // core losses
+re2=(Po/P)*100; // P U total resistance referred to l. v. side 
+xe2=(vr-re2*pff)/sqrt(1-pff^2); // P U total leakage reactance referred to l. v. side
+re2=(re2*E2^2)/(100*P); // total resistance in ohms 
+disp('Total resistance referred to l. v. side is ');
+printf('%f ohm',re2);
+xe2=(xe2*E2^2)/(100*P); // total leakage reactance in ohms
+disp('Total leakage reactance referred to l. v. side is ');
+printf('%f ohm',xe2);
+Rc=E2^2/Pc; 
+disp('Coreloss resistance is');
+printf('%f ohm',Rc);
+Ie2=Pc/(E2*pfn); // exciting current in Ampere
+Ic=Pc/E2; // core loss current
+Im=sqrt(Ie2^2-Ic^2); // magnetizing component of exciting current
+Xm=E2/Im; // magnetizing reactance
+disp('Magnetizing reactance is ');
+printf('%f ohm',Xm);
+disp('All parameters are known. So, equivalent circuit diagram referred to l. v. side can be drawn.');