summaryrefslogtreecommitdiff
path: root/20/CH2/EX2.3.69/example2_3_pg69.sce
blob: 145287625847ebd5af32baca84f8c3f549833ad4 (plain)
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
// Example2_3_pg69.sce
// To find primary voltage and current supplied
// Theory of Alternating Current Machinery by Alexander Langsdorf
// First Edition 1999, Thirty Second reprint
// Tata McGraw Hill Publishing Company
// Example in Page 69


clear; clc; close;

// Given data

// Transformer A data
va_A = 100e+3; // VA rating of Transformer
v1_A = 4600; // Voltage in volts
v2_A = 230; // Voltage in volts
x_A = 0.027; // Reactance in Ohms
r_A = 0.008; // Resistance in Ohms

// Transformer B data
va_B = 200e+3; // VA rating of Transformer
v1_B = 4610; // Voltage in volts
v2_B = 225; // Voltage in volts
x_B = 0.013; // Reactance in ohms
r_B = 0.003; // Resistance in ohms

// Common Data
P_load = 150e+3; // Power in Watts
pf = +0.85; // + denotes lagging power factor
vg = 225; // Voltage in volts


// Calculations

// Transformer A
a_1 = v1_A / v2_A;
z_1 = r_A + x_A*%i;
y_1 = 1 / z_1;
y_1_HVside = y_1 / a_1;

// Transformer B
a_2 = v1_B / v2_B;
z_2 = r_B + x_B*%i;
y_2 = 1 / z_2;
y_2_HVside = y_2 / a_2;

y_K = y_1 + y_2;
y_K_HVside = y_1_HVside + y_2_HVside;

// To find the current
I = P_load / (vg * pf) ;
V2_vec = vg;
theta = acos(0.85);
I_vec = I*(cos(theta) - sin(theta)*%i); // - sign indicates I lags V

V1_vec = ((V2_vec * y_K) + I_vec) / (y_K_HVside) ;

I1_vec = (I_vec + V1_vec*((y_K / a_1) - y_K_HVside)) / (z_1 * y_K );

I2_vec = I_vec - I1_vec;

printf(' Primary Voltage of transformer = %f /_ %f Volts\n', abs(V1_vec), (atan((imag(V1_vec))/(real(V1_vec))))*180/%pi);
printf(' Current Supplied by transformer A = %f /_ %f Volts\n', abs(I1_vec), (atan((imag(I1_vec))/(real(I1_vec))))*180/%pi);
printf(' Current Supplied by transformer B = %f /_ %f Volts\n', abs(I2_vec), (atan((imag(I2_vec))/(real(I2_vec))))*180/%pi);

// Result
// Primary Voltage of transformer = 4678.867698 /_ 1.211839 Volts 
// Current Supplied by transformer A = 361.324403 /_ -44.400715 Volts
// Current Supplied by transformer B = 438.858386 /_ -21.431553 Volts