From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:27:19 +0530 Subject: initial commit / add all books --- 3574/CH3/EX3.6/EX3_6.sce | 66 ++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 66 insertions(+) create mode 100644 3574/CH3/EX3.6/EX3_6.sce (limited to '3574/CH3/EX3.6/EX3_6.sce') diff --git a/3574/CH3/EX3.6/EX3_6.sce b/3574/CH3/EX3.6/EX3_6.sce new file mode 100644 index 000000000..5a57ccca1 --- /dev/null +++ b/3574/CH3/EX3.6/EX3_6.sce @@ -0,0 +1,66 @@ +// Example 3.6 +// Determine the percent of the total bank-current drawn by each transformer +// Page No. 109 + +clc; +clear; +close; + +// Given data +ZaPU_R=0.0158; // Transformer A impedance real part +ZaPU_I=0.0301; // Transformer A impedance imaginary part +ZbPU_R=0.0109; // Transformer B impedance real part +ZbPU_I=0.0398; // Transformer B impedance imaginary part +SB=200000; // Transformer B rating +VHSA=2400; // Voltage at the high side of transformer A +VHSB=2400; // Voltage at the high side of transformer B +RPUA=1.64; // Percent resistance of transformer A +XPUA=3.16; // Percent reactance of transformer A +RPUB=1.10; // Percent resistance of transformer B +XPUB=4.03; // Percent reactance of transformer B + + + +// Base impedance of transformer A +ZaPU=ZaPU_R+%i*ZaPU_I; +// Complex to Polar form... +ZaPU_Mag=sqrt(real(ZaPU)^2+imag(ZaPU)^2); // Magnitude part +ZaPU_Ang=atan(imag(ZaPU),real(ZaPU))*180/%pi; // Angle part
 + +// Base impedance of transformer B +ZbPU=ZbPU_R+%i*ZbPU_I; +// Complex to Polar form... +ZbPU_Mag=sqrt(real(ZbPU)^2+imag(ZbPU)^2); // Magnitude part +ZbPU_Ang=atan(imag(ZbPU),real(ZbPU))*180/%pi; // Angle part
 + +// Admittance of transformer A +YaPU_Mag=1/ZaPU_Mag; // Magnitude of equivalent admittance A +YaPU_Ang=0-ZaPU_Ang; // Angle of equivalent admittance A + +// Polar to Complex form + +YaPU_R=YaPU_Mag*cos(-YaPU_Ang*%pi/180); // Real part of complex number +YaPU_I=YaPU_Mag*sin(YaPU_Ang*%pi/180); //Imaginary part of complex number + +// Admittance of transformer B +YbPU_Mag=1/ZbPU_Mag; // Magnitude of equivalent admittance B +YbPU_Ang=0-ZbPU_Ang; // Angle of equivalent admittance B +// Polar to Complex form + +YbPU_R=YbPU_Mag*cos(-YbPU_Ang*%pi/180); // Real part of complex number +YbPU_I=YbPU_Mag*sin(YbPU_Ang*%pi/180); //Imaginary part of complex number + +// Parallel admittance +YP=(YaPU_R-%i*YaPU_I)+(YbPU_R-%i*YbPU_I); +// Complex to Polar form... +YP_Mag=sqrt(real(YP)^2+imag(YP)^2); // Magnitude part +YP_Ang=atan(imag(YP),real(YP))*180/%pi; // Angle part
 + +IA=YaPU_Mag/YP_Mag*100; // Percent current drawn by transformer A +IB=100-IA; + +// Display the result on the command window +printf("\n Percent of total bank current drawn by transformer A = %0.2f Percent ",IA); +printf("\n Percent of total bank current drawn by transformer B = %0.2f Percent ",IB); + + -- cgit