diff options
Diffstat (limited to '3556/CH12/EX12.10/Ex12_10.sce')
| -rw-r--r-- | 3556/CH12/EX12.10/Ex12_10.sce | 107 |
1 files changed, 107 insertions, 0 deletions
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)
+
+
+
|