1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
|
//Book Name:Fundamentals of Electrical Engineering
//Author:Rajendra Prasad
//Publisher: PHI Learning Private Limited
//Edition:Third ,2014
//Ex12_7.sce
clc;
clear;
z1=complex(8,6);
z2=complex(6,8);
z3=complex(4,-3);
R_YR=z1;
R_BY=z2;
R_RB=z3;
Vl=440;
printf("\n(a)Delta connected load of phase sequence RYB:")
theta1=0;
theta2=-120;
theta3=120;
V_YR=complex(Vl*cosd(theta1),Vl*sind(theta1));
V_BY=complex(Vl*cosd(theta2),Vl*sind(theta2));
V_RB=complex(Vl*cosd(theta3),Vl*sind(theta3));
I_YR=V_YR/z1;
I_BY=V_BY/z2;
I_RB=V_RB/z3;
I_YR_mag=sqrt(real(I_YR)^2+imag(I_YR)^2);
I_BY_mag=sqrt(real(I_BY)^2+imag(I_BY)^2);
I_RB_mag=sqrt(real(I_RB)^2+imag(I_RB)^2);
I_YR_angle=atand(imag(I_YR)/real(I_YR));
I_BY_angle=atand(imag(I_BY)/real(I_BY))-180;
I_RB_angle=atand(imag(I_RB)/real(I_RB))+180;
printf("\nPhase current= \tMagnitude\tAngle(deg) \n")
printf("\n\t Iyr in A \t %d \t %2.2f",I_YR_mag,I_YR_angle)
printf("\n\t Iby in A \t %d \t %2.2f",I_BY_mag,I_BY_angle)
printf("\n\t Irb in A \t %d \t %2.2f",I_RB_mag,I_RB_angle)
I_R=I_YR-I_RB;
I_Y=I_BY-I_YR;
I_B=I_RB-I_BY;
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_R_angle=atand(imag(I_R)/real(I_R));
I_Y_angle=atand(imag(I_Y)/real(I_Y))+180;
I_B_angle=atand(imag(I_B)/real(I_B))+180;
printf("\nLine current= \tMagnitude\tAngle(deg) ")
printf("\n\t Ir in A \t %2.2f %2.2f",I_R_mag,I_R_angle)
printf("\n\t Iy in A \t %2.2f \t %2.2f",I_Y_mag,I_Y_angle)
printf("\n\t Ib in A \t %2.2f \t %2.2f",I_B_mag,I_B_angle)
W_YR=(I_YR_mag)^2*real(z1);
W_BY=(I_BY_mag)^2*real(z2);
W_RB=(I_RB_mag)^2*real(z3);
printf("\n Toatal power dissipated:\n")
printf("\n\t W_YR=%d W",W_YR)
printf("\n\t W_BY=%d W",W_BY)
printf("\n\t W_RB=%d W",W_RB)
printf("\n\n(b)Delta connected load of phase sequence RBY:")
theta1=0;
theta2=120;
theta3=-120;
V_YR=complex(Vl*cosd(theta1),Vl*sind(theta1));
V_BY=complex(Vl*cosd(theta2),Vl*sind(theta2));
V_RB=complex(Vl*cosd(theta3),Vl*sind(theta3));
I_YR=V_YR/z1;
I_BY=V_BY/z2;
I_RB=V_RB/z3;
I_YR_mag=sqrt(real(I_YR)^2+imag(I_YR)^2);
I_BY_mag=sqrt(real(I_BY)^2+imag(I_BY)^2);
I_RB_mag=sqrt(real(I_RB)^2+imag(I_RB)^2);
I_YR_angle=atand(imag(I_YR)/real(I_YR));
I_BY_angle=atand(imag(I_BY)/real(I_BY));
I_RB_angle=atand(imag(I_RB)/real(I_RB));
printf("\nPhase current= \tMagnitude\tAngle(deg) \n")
printf("\n\t Iyr in A \t %d \t %2.2f",I_YR_mag,I_YR_angle)
printf("\n\t Iby in A \t %d \t %2.2f",I_BY_mag,I_BY_angle)
printf("\n\t Irb in A \t %d \t %2.2f",I_RB_mag,I_RB_angle)
I_R=I_YR-I_RB;
I_Y=I_BY-I_YR;
I_B=I_RB-I_BY;
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_R_angle=atand(imag(I_R)/real(I_R));
I_Y_angle=atand(imag(I_Y)/real(I_Y))+180;
I_B_angle=atand(imag(I_B)/real(I_B))-180;
printf("\nLine current= \tMagnitude\tAngle(deg) ")
printf("\n\t Ir in A \t %2.2f %2.2f",I_R_mag,I_R_angle)
printf("\n\t Iy in A \t %2.2f \t %2.2f",I_Y_mag,I_Y_angle)
printf("\n\t Ib in A \t %2.2f %2.2f",I_B_mag,I_B_angle)
W_YR=(I_YR_mag)^2*real(z1);
W_BY=(I_BY_mag)^2*real(z2);
W_RB=(I_RB_mag)^2*real(z3);
printf("\n Toatal power dissipated:\n")
printf("\n\t W_YR=%d W",W_YR)
printf("\n\t W_BY=%d W",W_BY)
printf("\n\t W_RB=%d W",W_RB)
|
