blob: 7362b8662e7f66fbd7fc38a72bc2ca7a3a78f264 (
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
|
clc; funcprot(0);
//Example 15.6 Thrust Horsepower of a constant speed propeller at Sea level
// Initialisation of variables
V = 125*1.467; // Velocity in ft/sec
n = 1800/60; // Rotaion per second
D = 9.5; // Diameter in ft
P = 600*550;
rho = 0.002378;
// Calculations
//For Beta = 23 degree
Cp = P/(rho*n^3*D^5);
ThrustCoeff = rho*n^2*D^4;
VeloCoeff = n*60*D/88;
VBynD = [1.03 1 0.9 0.8 0.7 0.6];
Ct = [0.055 0.056 0.062 0.068 0.076 0.086]; // From figure 15.4 for corrospondiong values of V/nD at Cp = 0.066
T = Ct*ThrustCoeff;
V = VBynD*VeloCoeff;
Result = zeros(6,6);
Result(:,1) = VBynD';
Result(:,2) = Ct';
Result(:,3) = Result(:,2)/Cp;
Result(:,4) = diag(diag(Result(:,3))*diag(Result(:,1))); // Using equation 15.5.5
Result(:,5) = Result(:,4)*P/550;
Result(:,6) = V';
//Results
disp(Result,"!! V/nD Ct Ct/Cp Efficiency thp V(mph) !!") ;
|