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
|
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 - 8
clear; clc; close; // Clear the work space and console.
//
// Given data
VLL = 240.0000;
f = 60.0000;
// Load 1
P1 = 30.00000;
pf1 = 0.60000;
// Load 2
Q2 = 45.00000;
pf2 = 0.80000;
// Calculations Complex, Real and Reactive Power of Load 1
S1 = P1/pf1;
Q1 = S1 * sqrt(1-((pf1)^2));
// Calculations Complex, Real and Reactive Power of Load 2
S2 = Q2/sqrt(1-((pf2)^2));
P2 = S2 * pf2;
// Calculations Total Complex, Real and Reactive Power of Load
ST = complex(P1,Q1) + complex(P2,Q2);
PT = real(ST);
QT = imag(ST);
ST_mag = norm(ST);
ST_angle = atand(QT,PT)
// Calculations Line Current For Load 1
IL1 = (50.0000 * 1000)/(sqrt(3)*240*1000);
Ia1_mag = IL1;
Ia1_angle = acosd(pf1);
// Calculations Line Current For Load 2
IL2 = (75.0000 * 1000)/(sqrt(3)*240*1000);
Ia2_mag = IL2;
Ia2_angle = acosd(pf2);
// Calculations Line Current for Load
Ia_tot_real = (Ia1_mag * cosd(Ia1_angle))+(Ia2_mag * cosd(Ia2_angle));
Ia_tot_imag = (Ia1_mag * sind(Ia1_angle))+(Ia2_mag * sind(Ia2_angle));
Ia_tot_mag = norm(complex(Ia_tot_real,Ia_tot_imag))
Ia_tot_angle = atand(Ia_tot_imag,Ia_tot_real)
// Calculations KVAR for Three Capacitors and Each Capacitor
angle_old = ST_angle;
angle_new = acosd(0.9000);
Qc = PT * (tand(angle_old)-tand(angle_new));
Qc1 = Qc/3;
// Calculations Capacitance Each Capacitor
C = Qc1/((2*(%pi)*f)*(VLL)^2)
//
disp("Example 12-8 Solution : ");
disp("a. Complex, Real and Reactive Power of Total Load : ");
disp("Complex, Real and Reactive Power of Load 1: ");
printf(" \n S1 = Complex Power of Load 1 = %.3f KVA",S1)
printf(" \n P1 = Real Power of Load 1 = %.3f KW",P1)
printf(" \n Q1 = Reactive Power of Load 1 = %.3f KVAR",Q1)
disp("")
disp("Complex, Real and Reactive Power of Load 2: ");
printf(" \n S2 = Complex Power of Load 2 = %.3f KVA",S2)
printf(" \n P2 = Real Power of Load 2 = %.3f KW",P2)
printf(" \n Q2 = Reactive Power of Load 2 = %.3f KVAR",Q2)
disp("")
disp("Complex, Real and Reactive Power of Total Load : ");
printf(" \n ST_mag = Magnitude of Complex Power of Total Load = %.3f KVA",ST_mag)
printf(" \n ST_angle = Angle of Complex Power of Total Load = %.3f Degree",ST_angle)
printf(" \n PT = Real Power of Total Load = %.3f KW",PT)
printf(" \n QT = Reactive Power of Total Load = %.3f KVAR",QT)
disp("")
disp("b. Line Current for Total Load : ");
disp("Line Current for Load 1: ");
printf(" \n Ia1_mag = Magnitude of Line Current For Load 1 = %.3f mA",Ia1_mag*1000)
printf(" \n Ia1_angle = Angle of Line Current For Load 1 = %.3f Degree",-Ia1_angle)
disp("")
disp("Line Current for Load 2: ");
printf(" \n Ia2_mag = Magnitude of Line Current For Load 2 = %.3f mA",Ia2_mag*1000)
printf(" \n Ia2_angle = Angle of Line Current For Load 2 = %.3f Degree",-Ia2_angle)
disp("")
disp("Total Line Current : ");
printf(" \n Ia2_mag = Magnitude of Total Line Current For Load = %.3f mA",Ia_tot_mag*1000)
printf(" \n Ia2_angle = Angle of Total Line Current For Load = %.3f Degree",-Ia_tot_angle)
disp("")
disp("c. Capacitance of Capasitor for Improving Power Factor : ");
printf(" \n Qc = KVAR for Three Capacitors = %.3f KVAR",Qc)
printf(" \n Qc1 = KVAR for Each Capacitor = %.3f KVAR",Qc1)
printf(" \n C = Capacitance of Capacitor = %.3f PF",C*1D+09)
|
