diff options
Diffstat (limited to '689/CH9')
| -rw-r--r-- | 689/CH9/EX8.3/3.sce | 32 | ||||
| -rw-r--r-- | 689/CH9/EX9.1/1.sce | 19 | ||||
| -rw-r--r-- | 689/CH9/EX9.2/2.sce | 16 |
3 files changed, 67 insertions, 0 deletions
diff --git a/689/CH9/EX8.3/3.sce b/689/CH9/EX8.3/3.sce new file mode 100644 index 000000000..a67630ce9 --- /dev/null +++ b/689/CH9/EX8.3/3.sce @@ -0,0 +1,32 @@ +clc; funcprot(0);
+//Example 9.3 Horseshoe Vortex
+
+// Initialisation of variables
+T = 250;
+
+// Calculations
+function[y] =Velocity(R, theta1, theta2)
+ y = (T/(4*%pi*R))*(cos(theta1)-cos(theta2))
+endfunction
+//Considering Leg AB
+ theta1 = 0; //From figure
+ theta2 = %pi - atan(4/5); //From figure
+ R = 4; //From figure
+ Vt1 = Velocity(R, theta1, theta2);
+
+//Considering Leg BC
+ theta1 = atan(5/4); //From figure
+ theta2 = %pi - atan(5/6); //From figure
+ R = 5; //From figure
+ Vt2 = Velocity(R, theta1, theta2);
+
+//Considering Leg CD
+ theta1 = atan(5/6); //From figure
+ theta2 = %pi; //From figure
+ R = 6; //From figure
+ Vt3 = Velocity(R, theta1, theta2);
+
+Vt = Vt1 + Vt2 + Vt3;
+
+//Results
+disp(Vt,"Velocity at point P (ft/sec):",Vt3,"Velocity due to leg CD (ft/sec):",Vt2,"Velocity due to leg BC (ft/sec):",Vt1,"Velocity due to leg AB (ft/sec):");
diff --git a/689/CH9/EX9.1/1.sce b/689/CH9/EX9.1/1.sce new file mode 100644 index 000000000..3f764a5de --- /dev/null +++ b/689/CH9/EX9.1/1.sce @@ -0,0 +1,19 @@ +clc; funcprot(0);
+//Example 9.1 Lift due to Circulation
+
+// Initialisation of variables
+D = 4;
+L = 12;
+V = 40*1.467;
+rho = 0.002378;
+W = 100/60; // Revolution per second
+
+// Calculations
+R =D/2;
+Vt = 2*%pi*R*W;
+T = 2*%pi*R*Vt;
+Lift = rho*T*V;
+L_total = Lift*L
+
+//Results
+disp(L_total,"Total lifting force (lb):");
diff --git a/689/CH9/EX9.2/2.sce b/689/CH9/EX9.2/2.sce new file mode 100644 index 000000000..2aea7f8a2 --- /dev/null +++ b/689/CH9/EX9.2/2.sce @@ -0,0 +1,16 @@ +clc; funcprot(0);
+//Example 9.2 Biot Savart's Law
+
+// Initialisation of variables
+T = 500;
+PB = 2;
+BA1 = 3;
+
+// Calculations
+theta1 = atan(PB/BA1);
+theta2 = %pi/2;
+R = PB;
+Vt = (T/(4*%pi*R))*(cos(theta1)-cos(theta2));
+
+//Results
+disp(Vt,"Velocity at point P (ft/sec):");
\ No newline at end of file |