diff options
Diffstat (limited to '797/CH11/EX11.6s/11_06_solution.sce')
| -rw-r--r-- | 797/CH11/EX11.6s/11_06_solution.sce | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/797/CH11/EX11.6s/11_06_solution.sce b/797/CH11/EX11.6s/11_06_solution.sce new file mode 100644 index 000000000..5641cda1f --- /dev/null +++ b/797/CH11/EX11.6s/11_06_solution.sce @@ -0,0 +1,23 @@ +//Solution 11-5 +WD=get_absolute_file_path('11_06_solution.sce'); +datafile=WD+filesep()+'11_06_example.sci'; +clc; +exec(datafile) +//unit conversions +D = D / 100; //from [cm] to [m] +V = V / 3.6 //from [km/h] to [m/s] +P = P * 1.01325 * 10**5; //from [atm] to [Pa] +T = T + 273; //from [C] to [K] +rho_air = P / (R * T); //from ideal gas equation +A = %pi / 4 * D**2; //frontal area of ball +omega = 2 * %pi * N / 60; //angular velocity of ball [rad/s] +nd_rotation = omega * D /(2 * V); //non dimensional rate of rotation +//from figure 11-53 lift coefficient coefficient corresponding to nd_rotation is +C_L = 0.21; +F_L = C_L * A * rho_air * V**2 / 2; //drag force +W = m * g; +if W > F_L then + printf("Ball will drop under the combined effect of gravity and lift due to spinning with net force of %1.3f - %1.3f = %1.3f N", W, F_L, W - F_L); +else + printf("Ball will rise under under the combined effect of gravity and lift due to spinning with net force of %1.3f - %1.3f = %1.3f N", F_L, W, F_L - W); +end |