blob: bcaf8b526ab5eae0388eb693c3859a13766a5146 (
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
|
//Pump Power//
pathname=get_absolute_file_path('10.03.sce')
filename=pathname+filesep()+'10.03-data.sci'
exec(filename)
[nQ mQ]= size(Q);
[nps mps]=size(ps);
[npd mpd]= size(pd);
[nI mI]= size(I);
//Correct measured static pressures to he pump centreline p1, p2(in psig):
j=1:mps;
p1=ps(j)+px*g*zs/144
j=1:mpd;
p2=pd(j)+px*g*zd/144
//The value of Pump head(in feet):
j=1:mps;
Hp=(p2(j)-p1(j))/(px*g)*144
//Values of Hydraulic Power delivered(in hp):
j=1:mps;
Wh=Q(j).*(p2(j)-p1(j))/7.48/60*144/550
//Values of motor power output(in hp):
j=1:mI;
Pin=Effm*sqrt(3)*PF*E*I(j)/746
//Values of Pump Efficiecy:
j=1:mI;
Effp= Wh(j)./Pin(j)*100
//Plotting pump characteristics:
plot(Q,Hp,"-o")
plot(Q,Pin,"-+")
plot(Q,Effp,"-*")
xtitle('Pump Characteristics','Volume flow rate(in gpm)',['Pump Efficincy(%) ',' Pump Head(in feet) ',' Pump Power input(in hp) '])
legend('Hp','Pin','Effp')
|