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 --- 3556/CH12/EX12.10/Ex12_10.sce | 107 ++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 107 insertions(+) create mode 100644 3556/CH12/EX12.10/Ex12_10.sce (limited to '3556/CH12/EX12.10/Ex12_10.sce') diff --git a/3556/CH12/EX12.10/Ex12_10.sce b/3556/CH12/EX12.10/Ex12_10.sce new file mode 100644 index 000000000..d33554222 --- /dev/null +++ b/3556/CH12/EX12.10/Ex12_10.sce @@ -0,0 +1,107 @@ +clc +// Fundamental of Electric Circuit +// Charles K. Alexander and Matthew N.O Sadiku +// Mc Graw Hill of New York +// 5th Edition + +// Part 2 : AC Circuits +// Chapter 12 : Three Phase Circuit +// Example 12 - 10 + +clear; clc; close; // Clear the work space and console. +// +// Given data +V1_mag = 120.0000 * sqrt(3); +V1_angle = 30.0000; +V2_mag = 120.0000 * sqrt(3); +V2_angle = -90.000; +Van_mag = 120.0000; +Van_angle = 0.0000; +Vbn_mag = Van_mag; +Vbn_angle = -120.0000; +Vcn_mag = Van_mag; +Vcn_angle = 120.0000; +Zan = complex(0,5); +Zbn = complex(10,0); +Zcn = complex(0,-10); +V1 = complex(V1_mag*cosd(V1_angle),V1_mag*sind(V1_angle)); +V2 = complex(V2_mag*cosd(V2_angle),V2_mag*sind(V2_angle)); +Z11 = complex(10,5) +Z12 = complex(-10,0); +Z21 = Z12; +Z22 = complex(10,-10); +// Calculations Determinants +delta = det([Z11 Z12; Z21 Z22]); +delta_mag = norm(delta); +delta_real = real(delta); +delta_imag = imag(delta); +delta_angle = atand(delta_imag,delta_real); +// Calculations Determinants For Source Voltage V1 +delta1 = det([V1 Z12; V2 Z22]); +delta1_mag = norm(delta1); +delta1_real = real(delta1); +delta1_imag = imag(delta1); +delta1_angle = atand(delta1_imag,delta1_real); +// Calculations Determinants For Source Voltage V2 +delta2 = det([Z11 V1; Z12 V2]); +delta2_mag = norm(delta2); +delta2_real = real(delta2); +delta2_imag = imag(delta2); +delta2_angle = atand(delta2_imag,delta2_real); +// Calculations Mesh Current I1 +I1_mag = delta1_mag/delta_mag; +I1_angle = delta1_angle - delta_angle +// Calculations Mesh Current I2 +I2_mag = delta2_mag/delta_mag; +I2_angle = delta2_angle - delta_angle +// Calculations Line Current a +I_a_mag = I1_mag; +I_a_angle = I1_angle +// Calculations Line Current b +I1 = complex(I1_mag*cosd(I1_angle),I1_mag*sind(I1_angle)) +I2 = complex(I2_mag*cosd(I2_angle),I2_mag*sind(I2_angle)) +I_b_mag = norm(I2 - I1); +I_b_angle = atand(imag(I2 - I1),real(I2 - I1)) +// Calculations Line Current c +I_c_mag = I2_mag +I_c_angle = -180.0000 + I2_angle; +// Calculations Power Absorbed by the Load for Phase A +Sal =(I_a_mag)^2*Zan; +// Calculations Power Absorbed by the Load for Phase B +Sbl =(I_b_mag)^2*Zbn; +// Calculations Power Absorbed by the Load for Phase C +Scl =(I_c_mag)^2*Zcn; +// Calculations Total Complex Power Absorbed by the Load +Stl = Sal + Sbl + Scl; +Stl_real = real(Stl); +Stl_imag = imag(Stl); +// Calculations Power Absorbed by the Source for Phase A +Sas =(complex(Van_mag*cosd(Van_angle),Van_mag*sind(Van_angle)))*conj(complex(I_a_mag*cosd(I_a_angle),I_a_mag*sind(I_a_angle))); +// Calculations Power Absorbed by the Load for Phase B +Sbs =(complex(Vbn_mag*cosd(Vbn_angle),Vbn_mag*sind(Vbn_angle)))*conj(complex(I_b_mag*cosd(I_b_angle),I_b_mag*sind(I_b_angle))); +// Calculations Power Absorbed by the Load for Phase C +Scs =(complex(Vcn_mag*cosd(Vcn_angle),Vcn_mag*sind(Vcn_angle)))*conj(complex(I_c_mag*cosd(I_c_angle),I_c_mag*sind(I_c_angle))); +// Calculations Total Complex Power Absorbed by The Source +Sts = Sas + Sbs + Scs; +Sts_real = -real(Sts); +Sts_imag = -imag(Sts); +// +disp("Example 12-10 Solution : "); +disp("a. The Line Currents: "); +printf(" \n I_a_mag = Magnitude of Line Currents a = %.3f A",I_a_mag) +printf(" \n I_a_angle = Angle of Line Currents a = %.3f degree",I_a_angle) +printf(" \n I_b_mag = Magnitude of Line Currents b = %.3f A",I_b_mag) +printf(" \n I_b_angle = Angle of Line Currents b = %.3f degree",I_b_angle) +printf(" \n I_c_mag = Magnitude of Line Currents c = %.3f A",I_c_mag) +printf(" \n I_c_angle = Angle of Line Currents c = %.3f degree",I_c_angle) +disp("") +disp("b Total Power Complex Absorbed By The Load"); +printf(" \n Stl_real = Real Part of Power Complex = %.3f Watt",Stl_real) +printf(" \n Stl_imag = Imaginary Part of Power Complex = %.3f Var",Stl_imag) +disp("") +disp("c. Total Power Complex Absorbed By The Source"); +printf(" \n Sts_real = Real Part of Power Complex = %.3f Watt",Sts_real) +printf(" \n Sts_imag = Imaginary Part of Power Complex = %.3f Var",Sts_imag) + + + -- cgit