blob: 214c9e24e4ad8ec097bfa09a40a4d2b8b328d0eb (
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
|
//volume and power//
pathname=get_absolute_file_path('10.02.sce')
filename=pathname+filesep()+'10.02-data.sci'
exec(filename)
U=0.5*(Dh+Dt)/2*1200*2*%pi/60
k=tand(alpha1)+cotd(betta1)
Vn1=U/k
V1=Vn1/cosd(alpha1)
Vt1=V1*sind(alpha1)
Vrb1=Vn1/sind(betta1)
//Volume flow rate (in m^3/sec):
Q=%pi/4*Vn1*(Dt^2-Dh^2)
k=(U-Vn1*cotd(betta2))/Vn1
alpha2= atand(k)
V2=Vn1/cosd(alpha2)
Vt2=V2*sind(alpha2)
//Rotor Torque (in N-m):
Tz=p*Q*(Dh+Dt)/4*(Vt2-Vt1)
//Power required (in W):
Wm=w*2*%pi/60*Tz
printf("\n\nRESULTS\n\n")
printf("\n\nVolume flow rate: %.3f m^3/sec\n\n",Q)
printf("\n\nRotor Torque: %.3f N-m\n\n",Tz)
printf("\n\nPower required: %.3f W\n\n",Wm)
|
