summaryrefslogtreecommitdiff
path: root/3802/CH12/EX12.10/Ex12_10.sce
diff options
context:
space:
mode:
Diffstat (limited to '3802/CH12/EX12.10/Ex12_10.sce')
-rw-r--r--3802/CH12/EX12.10/Ex12_10.sce62
1 files changed, 62 insertions, 0 deletions
diff --git a/3802/CH12/EX12.10/Ex12_10.sce b/3802/CH12/EX12.10/Ex12_10.sce
new file mode 100644
index 000000000..4396af82b
--- /dev/null
+++ b/3802/CH12/EX12.10/Ex12_10.sce
@@ -0,0 +1,62 @@
+//Book Name:Fundamentals of Electrical Engineering
+//Author:Rajendra Prasad
+//Publisher: PHI Learning Private Limited
+//Edition:Third ,2014
+
+//Ex12_10.sce
+
+clc;
+clear;
+Z_R=complex(8,6);
+Z_Y=complex(8,-6);
+Z_B=complex(5,0);
+Z_N=complex(0.5,1);
+Y_R=1/Z_R;
+Y_Y=1/Z_Y;
+Y_B=1/Z_B;
+Y_N=1/Z_N;
+E_R=220;
+E_Y=220;
+E_B=220;
+theta1=0;
+theta2=-120;
+theta3=120;
+V_R=complex(E_R*cosd(theta1),E_R*sind(theta1));
+V_Y=complex(E_Y*cosd(theta2),E_Y*sind(theta2));
+V_B=complex(E_B*cosd(theta3),E_B*sind(theta3));
+V_NN_dash=((V_R*Y_R)+(V_Y*Y_Y)+(V_B*Y_B))/(Y_R+Y_Y+Y_B+Y_N);
+
+V_R_dash=V_R-V_NN_dash;
+V_Y_dash=V_Y-V_NN_dash;
+V_B_dash=V_B-V_NN_dash;
+V_R_dash_mag=sqrt(real(V_R_dash)^2+imag(V_R_dash)^2);
+V_Y_dash_mag=sqrt(real(V_Y_dash)^2+imag(V_Y_dash)^2);
+V_B_dash_mag=sqrt(real(V_B_dash)^2+imag(V_B_dash)^2);
+V_R_dash_angle=atand(imag(V_R_dash)/real(V_R_dash));
+V_Y_dash_angle=atand(imag(V_Y_dash)/real(V_Y_dash))+180;
+V_B_dash_angle=atand(imag(V_B_dash)/real(V_B_dash))+180;
+printf("\n Load phase voltages: Magnitude\tAngle(deg)")
+printf("\n For R phase\t%3.2f\t%0.3f",V_R_dash_mag,V_R_dash_angle)
+printf("\n For Y phase\t%3.2f\t%3.2f",V_Y_dash_mag,V_Y_dash_angle)
+printf("\n For B phase\t%3.2f\t%3.2f",V_B_dash_mag,V_B_dash_angle)
+//For V_NN_dash value , the answer given in the book is wrong.So load phase voltage vary from the book answer.
+//Also V_R_dash angle is not 0.168. It is negative angle that is -0.193
+I_R=V_R_dash*Y_R;
+I_Y=V_Y_dash*Y_Y;
+I_B=V_B_dash*Y_B;
+I_N=V_NN_dash*Y_N;
+I_R_mag=sqrt(real(I_R)^2+imag(I_R)^2);
+I_Y_mag=sqrt(real(I_Y)^2+imag(I_Y)^2);
+I_B_mag=sqrt(real(I_B)^2+imag(I_B)^2);
+I_N_mag=sqrt(real(I_N)^2+imag(I_N)^2);
+I_R_angle=atand(imag(I_R)/real(I_R));
+I_Y_angle=atand(imag(I_Y)/real(I_Y))+360;
+I_B_angle=atand(imag(I_B)/real(I_B))+180;
+I_N_angle=atand(imag(I_N)/real(I_N))+180;
+printf("\n\n Load phase current: Magnitude\tAngle(deg)")
+printf("\n For R phase\t%3.2f\t%0.3f",I_R_mag,I_R_angle)
+printf("\n For Y phase\t%3.2f\t%3.2f",I_Y_mag,I_Y_angle)
+printf("\n For B phase\t%3.2f\t%3.2f",I_B_mag,I_B_angle)
+printf("\n For Neutral\t%3.2f\t%3.2f",I_N_mag,I_N_angle)
+
+