blob: 8f65ab20c141fab34ae4a6110d3369826cb87680 (
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
|
//New plant pf and percent decrease in line current
clc;
clear;
Pmp=5000*(10^3); // Electrical load
pfmp=0.8; // Lag
Pim=500*735;// One horse power is 735W
Effim=96/100; // Efficiency of the motor
pfim=0.9; // Lag
pfsm=0.8; // Lead
Pime=Pim/Effim;// Effective power delivered by the induction motor
deff('x=com(y,z)','x=y+(%i*y*tand(acosd(z)))')// Function to find the complex powers
//Complex Powers
Pcmp=com(Pmp,pfmp); // Manufacturing Plant Load
Pcim=com(Pime,pfim);// Induction Motor
Pcsm=com(Pime,-pfsm);// Synchronous Machine, Minus Sign indicates Lead
Pr=Pcmp-Pcim+Pcsm; // Plant Requirement after replacement
pfar=real(Pr)/abs(Pr); // New Power Factor of the plant
Pnp=abs(Pr);
Vl=poly([0 1],'Vl','c');
Io=Pmp/(pfmp*sqrt(3)*Vl);
In=Pnp/(sqrt(3)*Vl); // Improved Factor Value =1;
red=(Io-In)*100/Io; // Reduction percent in fractions
redeq=Vl-red;// Reduction percent in decimal characteristic equation
redper=roots(redeq(2));
printf('The New Power Factor of the plant = %g lag \n',pfar )
printf('The Percentage decrease in line current that will result in improved p.f = %g percent \n',redper)
|