diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /689/CH16/EX16.1 | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '689/CH16/EX16.1')
| -rw-r--r-- | 689/CH16/EX16.1/1.pdf | bin | 0 -> 20520 bytes | |||
| -rw-r--r-- | 689/CH16/EX16.1/1.sce | 39 |
2 files changed, 39 insertions, 0 deletions
diff --git a/689/CH16/EX16.1/1.pdf b/689/CH16/EX16.1/1.pdf Binary files differnew file mode 100644 index 000000000..25d86c55f --- /dev/null +++ b/689/CH16/EX16.1/1.pdf diff --git a/689/CH16/EX16.1/1.sce b/689/CH16/EX16.1/1.sce new file mode 100644 index 000000000..e4aaee889 --- /dev/null +++ b/689/CH16/EX16.1/1.sce @@ -0,0 +1,39 @@ +clc; funcprot(0);
+//Example 16.1 Horsepower required at sea level
+// Initialisation of variables
+W = 2000;
+b = 36;
+c = 6;
+Dp = 3.8; // Parasite drag equivalent
+
+// Calculations
+S = b*c;
+WingLoading = W/S;
+VeloCoeff = sqrt(WingLoading/0.00256)
+Hp_WingCoeff = 0.00256*S/375;
+Hp_parCoeff = 0.00327*Dp/375;
+
+alp = [-4 -3 -2 -1 0 4 8 12 16 18 19 20];
+Cl = [0.07 0.14 0.215 0.285 0.36 0.6455 0.93 1.19 1.435 1.545 1.560 1.540]; // Values from fig 8.8
+Cd = [0.010 0.010 0.012 0.014 0.017 0.033 0.060 0.095 0.139 0.164 0.180 0.206]; // Values from fig 8.8
+
+
+Result = zeros(12,7);
+Result(:,1) = alp';
+Result(:,2) = Cl';
+Result(:,3) = Cd';
+Result(:,4) = 60.0*diag(inv(diag(sqrt(Cl)'))); // Using equation 15.5.5
+Result(:,5) = diag(diag(Result(:,3))*diag(Result(:,4)^3))*Hp_WingCoeff;
+Result(:,6) = Result(:,4)^3*Hp_parCoeff;
+Result(:,7) = Result(:,5) + Result(:,6);
+
+//Results
+disp(Result,"!! alpha Cl Cd V HP wing HP Par HP Total !!") ;
+clf();
+plot2d(Result(:,4),[Result(:,5) Result(:,6) Result(:,7)]);
+legend(['HP Req Wing'; 'HP Req Par'; 'HP Req Total'],2);
+xlabel("Miles Per Hour");
+ylabel("HorsePower");
+title("Horsepower required for various airspeeds ");
+set(gca(),"grid",[1 1])
+
|