blob: 7942aaeef8bf2b41804f8da68cf810446089bf5d (
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
|
clear;
clc;
funcprot(0);
//given data
//data from Exampla 10.5
Z = 3;//number of blades
D = 30;//rotor diameter in m
J = 5.0;//tip-speed ratio
l = 1.0;//blade chord in m
beta = 2;//pitch angle in deg
omega = 2.5;//in rad/s
rho = 1.2;//density in kg/m^3
cx1 = 7.5;//in m/s
sum_var1 = 6.9682;//from Table 10.3
sum_var2 = 47.509*10^-3;//from Table 10.4
//Calculations
X = sum_var1*0.5*rho*Z*l*0.5*D*cx1^2;
tau = sum_var2*0.5*rho*Z*l*(omega^2)*(0.5*D)^4;
P = tau*omega;
A2 = 0.25*%pi*D^2;
P0 = 0.5*rho*A2*cx1^3;
Cp = P/P0;
zeta = (27/16)*Cp;
//Results
printf('The total axial force = %d N.',X);
printf('\n The torque = %.3f *10^3 Nm.',tau/1000);
printf('\n The power developed = %.3f kW.',P/1000);
printf('\n The power coefficient = %.3f',Cp);
printf('\n The relative power coefficient = %.3f',zeta);
|
