initial commit / add all books

26 files changed, 283 insertions, 0 deletions

diff --git a/629/CH6/EX6.1/ex6_1.txt b/629/CH6/EX6.1/ex6_1.txt
new file mode 100644
index 000000000..f5f532c42
--- /dev/null
+++ b/629/CH6/EX6.1/ex6_1.txt
@@ -0,0 +1,3 @@
+

+ The force acting on the beam that supports the rocket, Fb = 7.56 N.

+ 
\ No newline at end of file
diff --git a/629/CH6/EX6.1/example6_1.sce b/629/CH6/EX6.1/example6_1.sce
new file mode 100644
index 000000000..3cadf3587
--- /dev/null
+++ b/629/CH6/EX6.1/example6_1.sce
@@ -0,0 +1,14 @@
+clear

+clc

+//Example 6.1 THRUST OF ROCKET

+g=9.81; //[m/s^2]

+m=0.04; //mass[kg]

+D=0.01; //[m]

+A=%pi*D^2/4 //area[m^2]

+rho=0.5; //density[kg/m^3]

+v=450; //[m/s]

+//Sum of forces, Fz=-Fb-m.g

+Mo=-rho*A*v^2 //momentum outflow[N]

+Fz=Mo //[N]

+Fb=-Fz-m*g //Force on beam[N]

+printf("\n The force acting on the beam that supports the rocket, Fb = %.2f N.\n",Fb)
\ No newline at end of file
diff --git a/629/CH6/EX6.10/ex6_10.txt b/629/CH6/EX6.10/ex6_10.txt
new file mode 100644
index 000000000..84df14937
--- /dev/null
+++ b/629/CH6/EX6.10/ex6_10.txt
@@ -0,0 +1,3 @@
+ 

+ The frictional force acting on the block = 313 N.

+ 
\ No newline at end of file
diff --git a/629/CH6/EX6.10/example6_10.sce b/629/CH6/EX6.10/example6_10.sce
new file mode 100644
index 000000000..1d9a468f5
--- /dev/null
+++ b/629/CH6/EX6.10/example6_10.sce
@@ -0,0 +1,10 @@
+clear

+clc

+//Example 6.10 JET IMPINGING ON MOVING BLOCK

+A=5*10^-4; //area[m^2]

+rho=1000; //density[kg/m^3]

+vj=50; //[m/s]

+vb=25; //[m/s]

+//Force on cart

+F=rho*A*(vj-vb)^2 //[N]

+printf("\n The frictional force acting on the block = %.f N.\n",F)
\ No newline at end of file
diff --git a/629/CH6/EX6.11/ex6_11.txt b/629/CH6/EX6.11/ex6_11.txt
new file mode 100644
index 000000000..a24bc3c63
--- /dev/null
+++ b/629/CH6/EX6.11/ex6_11.txt
@@ -0,0 +1,2 @@
+ 

+ The ratio of propellant mass to initial mass to achieve orbital velocity = 0.907.
\ No newline at end of file
diff --git a/629/CH6/EX6.11/example6_11.sce b/629/CH6/EX6.11/example6_11.sce
new file mode 100644
index 000000000..21df5043f
--- /dev/null
+++ b/629/CH6/EX6.11/example6_11.sce
@@ -0,0 +1,10 @@
+clear

+clc

+//Example 6.11 PROPELLANT MASS RATIO FOR ACHIEVING ORBITAL VELOCITY

+Vbo=7600; //orbital velocity[m/s]

+Isp=3200; //specific impulse[m/s]

+//Vbo=Isp*log(mi/mf)

+mif=exp(Vbo/Isp) //mif=mi/mf

+//mp=mi-mf

+mpi=1-1/mif //mpi=mp/mi

+printf("\n The ratio of propellant mass to initial mass to achieve orbital velocity = %.3f.\n",mpi)
\ No newline at end of file
diff --git a/629/CH6/EX6.12/ex6_12.txt b/629/CH6/EX6.12/ex6_12.txt
new file mode 100644
index 000000000..592bbab54
--- /dev/null
+++ b/629/CH6/EX6.12/ex6_12.txt
@@ -0,0 +1,3 @@
+ 

+ The maximum pressure that develops at the downstream end = 303 psig.

+ 
\ No newline at end of file
diff --git a/629/CH6/EX6.12/example6_12.sce b/629/CH6/EX6.12/example6_12.sce
new file mode 100644
index 000000000..5c8c05bc1
--- /dev/null
+++ b/629/CH6/EX6.12/example6_12.sce
@@ -0,0 +1,16 @@
+clear

+clc

+//Example 6.12 PRESSURE RISE DUE TO WATER HAMMER EFFECT

+rho=1.94; //[slugs/ft^3]

+Ev=3.2*10^5; //[lbf/in^2]

+V=4; //[ft/s]

+//Sound speed

+c=sqrt(Ev*144/rho) //[ft/s]

+L=3000; //[ft]

+tc=2*L/c //[s]

+//Closure time of 1sec is less than tc

+//Pressure rise

+delp=rho*V*c/144 //[psi]

+pi=40; //initial pressure[psi]

+pmax=pi+delp //[psi]

+printf("\n The maximum pressure that develops at the downstream end = %.f psig.\n",pmax)
\ No newline at end of file
diff --git a/629/CH6/EX6.13/ex6_13.txt b/629/CH6/EX6.13/ex6_13.txt
new file mode 100644
index 000000000..95196410e
--- /dev/null
+++ b/629/CH6/EX6.13/ex6_13.txt
@@ -0,0 +1,2 @@
+ 

+ The moment that the support system must resist, M = 3.61 kN.m.
\ No newline at end of file
diff --git a/629/CH6/EX6.13/example6_13.sce b/629/CH6/EX6.13/example6_13.sce
new file mode 100644
index 000000000..e0eee6d1c
--- /dev/null
+++ b/629/CH6/EX6.13/example6_13.sce
@@ -0,0 +1,32 @@
+clear

+clc

+//Example 6.13 RESISTING MOMENT ON REDUCING BEND

+r1=0.15; //[m]

+r2=0.475; //[m]

+d1=0.3; //[m]

+d2=0.15; //[m]

+A1=%pi*d1^2/4 //[m^2]

+A2=%pi*d2^2/4 //[m^2]

+p1=150*10^3; //[Pa]

+p2=59.3*10^3; //[Pa]

+//Torque due to pressure

+Mp=(r1*p1*A1)+(r2*p2*A2) //[N.m]

+

+rho=998; //density[kg/m^3]

+Q=0.25; //discharge[m^3/s]

+m=rho*Q //mass flow rate[kg/s]

+v1=Q/A1 //[m/s]

+v2=Q/A2 //[m/s]

+Mi=-m*r1*v1 //moment due to inflow

+Mo=m*r2*v2 //moment due to outflow

+//Moment of momentum flow

+Mm=Mo-Mi

+

+d=0.2; //[m]

+W=1420; //weight[N]

+//Moment due to weight

+Mw=d*W

+

+//Moment exerted by support

+M=(-Mp-Mm+Mw)*10^-3 //[kN.m]

+printf("\n The moment that the support system must resist, M = %.2f kN.m.\n",-M)
\ No newline at end of file
diff --git a/629/CH6/EX6.14/ex6_14.txt b/629/CH6/EX6.14/ex6_14.txt
new file mode 100644
index 000000000..8977ab4ce
--- /dev/null
+++ b/629/CH6/EX6.14/ex6_14.txt
@@ -0,0 +1,2 @@
+ 

+ The power produced by the turbine = 343 kW.
\ No newline at end of file
diff --git a/629/CH6/EX6.14/example6_14.sce b/629/CH6/EX6.14/example6_14.sce
new file mode 100644
index 000000000..8d7746258
--- /dev/null
+++ b/629/CH6/EX6.14/example6_14.sce
@@ -0,0 +1,17 @@
+clear

+clc

+//Example 6.14 POWER DELIVERED BY A FRANCIS TURBINE

+D=1; //[m]

+l=0.04; //[m]

+Q=0.5; //discharge[m^3/s]

+rho=1000; //[kg/m^3]

+m=rho*Q //mass flow rate[kg/s]

+//Radial velocity

+Vr=Q/(%pi*D*l) //[m/s]

+theta=70; //degrees

+//Tangential velocity

+Vt=Vr*tand(theta) //[m/s]

+T=m*(D/2)*Vt //Torque[N.m]

+w=1200*2*%pi/60 //angular speed (rad/s)

+P=T*w/10^3 //Power[kW]

+printf("\n The power produced by the turbine = %.f kW.\n",P)
\ No newline at end of file
diff --git a/629/CH6/EX6.2/ex6_2.txt b/629/CH6/EX6.2/ex6_2.txt
new file mode 100644
index 000000000..92e21759e
--- /dev/null
+++ b/629/CH6/EX6.2/ex6_2.txt
@@ -0,0 +1,6 @@
+ 

+ The tension in the cable = 233 lbf.

+

+

+ The weight recorded by the scale = 1203 lbf.

+ 
\ No newline at end of file
diff --git a/629/CH6/EX6.2/example6_2.sce b/629/CH6/EX6.2/example6_2.sce
new file mode 100644
index 000000000..6c82f8cd5
--- /dev/null
+++ b/629/CH6/EX6.2/example6_2.sce
@@ -0,0 +1,18 @@
+clear

+clc

+//Example 6.2 CONCRETE FLOWING INTO CART

+A=1; //[ft^2]

+rho=150; //density[lbm/ft^3]

+v=10; //speed[ft/s]

+theta=60; //degrees

+//Momentum accumulation=0, outflow=0

+//Momentum inflow

+mi_x=rho*A*v^2*cosd(theta) //x-direction

+mi_z=rho*A*v^2*sind(theta) //z-direction

+//1slug=32.2 lbm

+//Tension in cable

+T=mi_x/32.2 //[lbf]

+printf("\n The tension in the cable = %.f lbf.\n\n",T)

+W=800; //weight[lbf]

+N=W+(mi_z/32.2) //[lbf]

+printf("\n The weight recorded by the scale = %.f lbf.\n",N)
\ No newline at end of file
diff --git a/629/CH6/EX6.3/ex6_3.txt b/629/CH6/EX6.3/ex6_3.txt
new file mode 100644
index 000000000..520c8f094
--- /dev/null
+++ b/629/CH6/EX6.3/ex6_3.txt
@@ -0,0 +1,5 @@
+ 

+ The air speed at the exit of the nozzle = 77.9 m/s.

+

+

+ The force on the flange = 9.90 N.
\ No newline at end of file
diff --git a/629/CH6/EX6.3/example6_3.sce b/629/CH6/EX6.3/example6_3.sce
new file mode 100644
index 000000000..1d9d2d663
--- /dev/null
+++ b/629/CH6/EX6.3/example6_3.sce
@@ -0,0 +1,16 @@
+clear

+clc

+//Example 6.3 FORCE ON A NOZZLE

+d1=0.06; //[m]

+A1=%pi*d1^2/4 //area[m^2]

+d2=0.01; //[m]

+rho=1.22; //density[kg/m^3]

+p1=3.7*1000; //[Pa]

+//Bernoulli equation, p1+rho*v1^2/2=rho*v2^2/2

+v2=sqrt(2*p1/(rho*(1-(d2/d1)^4))) //Exit velocity[m/s]

+printf("\n The air speed at the exit of the nozzle = %.1f m/s.\n\n",v2)

+v1=v2*(d2/d1)^2 //Inlet velocity[m/s]

+m=rho*A1*v1 //mass flow rate[kg/s]

+//Force on flange

+F=m*(v2-v1)-p1*A1 //[N]

+printf("\n The force on the flange = %.2f N.\n",-F)
\ No newline at end of file
diff --git a/629/CH6/EX6.4/ex6_4.txt b/629/CH6/EX6.4/ex6_4.txt
new file mode 100644
index 000000000..8604e0d25
--- /dev/null
+++ b/629/CH6/EX6.4/ex6_4.txt
@@ -0,0 +1,2 @@
+ 

+ The force exerted by the jet on the vane, F = (53.0 lbf)i+(91.8 lbf)j.
\ No newline at end of file
diff --git a/629/CH6/EX6.4/example6_4.sce b/629/CH6/EX6.4/example6_4.sce
new file mode 100644
index 000000000..89516873a
--- /dev/null
+++ b/629/CH6/EX6.4/example6_4.sce
@@ -0,0 +1,16 @@
+clear

+clc

+//Example 6.4 WATER DEFLECTED BY A VANE

+r=0.0417; //[ft]

+A=%pi*r^2 //area[ft^2]

+v=100; //velocty[ft/s]

+rho=1.94; //[slugs/ft^3]

+m=rho*A*v; //mass flow rate[slugs/s]

+theta=60;//degrees

+//Momentum outflow vector, mo=[mo_x mo_y]

+mo=[m*v*cosd(theta) -m*v*sind(theta)]

+//Momentum inflow vector, mi=[mi_x mi_y]

+mi=[m*v 0]

+//Force vector, F=[Fx Fy] 

+F=mi-mo //[lbf]

+printf("\n The force exerted by the jet on the vane, F = (%.1f lbf)i+(%.1f lbf)j.\n",F(1),F(2))
\ No newline at end of file
diff --git a/629/CH6/EX6.6/ex6_6.txt b/629/CH6/EX6.6/ex6_6.txt
new file mode 100644
index 000000000..1b58651f2
--- /dev/null
+++ b/629/CH6/EX6.6/ex6_6.txt
@@ -0,0 +1,3 @@
+

+The net force required to hold the bend in place, F =(-8.53 kN)i+(-31.8 kN)j+(15.1 kN)k.

+ 
\ No newline at end of file
diff --git a/629/CH6/EX6.6/example6_6.sce b/629/CH6/EX6.6/example6_6.sce
new file mode 100644
index 000000000..ea9d6a092
--- /dev/null
+++ b/629/CH6/EX6.6/example6_6.sce
@@ -0,0 +1,24 @@
+clear

+clc

+//Example 6.6 FORCES ACTING ON A PIPE BEND

+p=75*10^3 //[Pa]

+r=0.5; //[m]

+A=%pi*r^2 //area[m^2]

+S=0.94;

+rho=S*1000 //drnsity[kg/m^3]

+Gamma=S*9.81 //specific weight of oil[kN/m^3]

+V=1.2; //volume of oil[m^3]

+Q=2; //[m^3/s]

+m=rho*Q //mass flow rate[kg/s]

+v=2.55; //[m/s]

+theta=30; //degrees

+//Reaction force

+//Rx+p*A-p*A*cos(theta)=m*v*cos(theta)-m*v

+Rx=-(p*A+m*v)*(1-cosd(theta))/10^3 //[kN]

+//Ry+p*A*sin(theta)=-m*v*sin(theta)

+Ry=-(p*A+m*v)*sind(theta)/10^3 //[kN]

+We=4; //empty weight of bend[kN]

+Rz=(Gamma*V)+We //[kN]

+//Reaction force vector

+R=[Rx Ry Rz] //[kN]

+printf("\n The net force required to hold the bend in place, F =(%.2f kN)i+(%.1f kN)j+(%.1f kN)k.\n",R(1),R(2),R(3))
\ No newline at end of file
diff --git a/629/CH6/EX6.7/ex6_7.txt b/629/CH6/EX6.7/ex6_7.txt
new file mode 100644
index 000000000..295731dbc
--- /dev/null
+++ b/629/CH6/EX6.7/ex6_7.txt
@@ -0,0 +1,2 @@
+ 

+ The components of force, F required to hold the bend in place are Fx = -16.1 kN, Fz = 1.48 kN.
\ No newline at end of file
diff --git a/629/CH6/EX6.7/example6_7.sce b/629/CH6/EX6.7/example6_7.sce
new file mode 100644
index 000000000..67d05ba75
--- /dev/null
+++ b/629/CH6/EX6.7/example6_7.sce
@@ -0,0 +1,27 @@
+clear

+clc

+//Example 6.7 WATER FLOW THROUGH REDUCING BEND

+d1=0.3; //[m]

+d2=0.15; //[m]

+A1=%pi*d1^2/4 //[m^2]

+A2=%pi*d2^2/4 //[m^2]

+rho=1000; //[kg/m^3]

+Q=0.25; //[m^3/s]

+v1=Q/A1 //inlet speed[m/s]

+v2=Q/A2 //outlet speed[m/s]

+p1=150; //[kPa]

+g_w=9810; //specific weight of water

+del_z=0.325; //(del_z=z1-z2)[m]

+//Bernoulli equation, p1+(rho*v1^2/2)+(g_w*z1)=p2+(rho*v2^2/2)+(g_w*z2)

+p2=p1+(rho*(v1^2-v2^2)/2+g_w*del_z)*10^-3 //[kPa]

+//Pressure forces

+Fp=(p1*A1)+(p2*A2) //[kN]

+//Momentum flux

+Fm=rho*Q*(v1+v2)*10^-3 //[kN]

+Wb=500; //weight of metal in the bend[N]

+V=0.1; //bend volume[m^3]

+Wf=g_w*V //weight of water[N]

+//Reaction force components

+Rx=-Fp-Fm //[kN]

+Rz=(Wb+Wf)*10^-3 //[kN]

+printf("\n The components of force, F required to hold the bend in place are Fx = %.1f kN, Fz = %.2f kN.\n",Rx,Rz)
\ No newline at end of file
diff --git a/629/CH6/EX6.8/ex6_8.txt b/629/CH6/EX6.8/ex6_8.txt
new file mode 100644
index 000000000..f9985969f
--- /dev/null
+++ b/629/CH6/EX6.8/ex6_8.txt
@@ -0,0 +1,2 @@
+ 

+ The drag force on the device and support panes = 304 N.
\ No newline at end of file
diff --git a/629/CH6/EX6.8/example6_8.sce b/629/CH6/EX6.8/example6_8.sce
new file mode 100644
index 000000000..591d0f328
--- /dev/null
+++ b/629/CH6/EX6.8/example6_8.sce
@@ -0,0 +1,18 @@
+clear

+clc

+//Example 6.8 DRAG FORCE ON WIND-TUNNEL MODEL

+ro=0.5; //radius of tunnel[m]

+A1=%pi*ro^2 //[m^2]

+p1=1.5*10^3; //[Pa]

+p2=10^3; //[Pa]

+v1=30; //velocity at inlet[m/s]

+rho=1; //[kg/m^3]

+//velocity profile, v=v1*K(r/ro)

+//Q1=Q, A1*v1=A*v

+K=A1*v1*(integrate('v1*(r/ro)*2*%pi*r','r',0,ro))^-1

+F1=rho*A1*v1^2 //momentum at cross-section-1

+F2=integrate('rho*(v1*K*(r/ro))^2*2*%pi*r','r',0,ro) //at cross-section-2

+Fx=F2-F1

+//From momentum equation in x-direction, Fx=p1*A-p2*A-Fd

+Fd=(p1-p2)*A1-Fx //Drag force[N]

+printf("\n The drag force on the device and support panes = %.f N.\n",Fd)
\ No newline at end of file
diff --git a/629/CH6/EX6.9/ex6_9.txt b/629/CH6/EX6.9/ex6_9.txt
new file mode 100644
index 000000000..62e4e7387
--- /dev/null
+++ b/629/CH6/EX6.9/ex6_9.txt
@@ -0,0 +1,5 @@
+ 

+ The flow rate under the sluice gate = 2008 ft^3/s.

+

+ The force on the sluice gate = 66.6 tons.

+ 
\ No newline at end of file
diff --git a/629/CH6/EX6.9/example6_9.sce b/629/CH6/EX6.9/example6_9.sce
new file mode 100644
index 000000000..c70d35a89
--- /dev/null
+++ b/629/CH6/EX6.9/example6_9.sce
@@ -0,0 +1,25 @@
+clear

+clc

+//Example 6.9 FORCE ON A SLUICE GATE

+g=32.2; //[ft/s^2]

+d1=20; //[ft]

+d2=3; //[ft]

+w=20; //gate width[ft]

+v2=sqrt((2*g*(d1-d2))/(1-(d2/d1)^2)) //[ft/s]

+v1=d2*v2/d1 //[ft/s]

+rho=1.94; //[slugs/ft^3]

+Q=v2*d2*w //discharge[ft^3/s]

+printf("\n The flow rate under the sluice gate = %.f ft^3/s.\n",Q)

+m=rho*Q //mass flow rate[slugs/s]

+Gamma=62.4; //[lbf/ft^3]

+F1=Gamma*w*(d1^2)/2

+F2=Gamma*w*(d2^2)/2

+//momentum inflow

+mi=m*v1 //[lbf]

+//momentum outflow

+mo=m*v2 //[lbf]

+Fx=mo-mi //[lbf]

+//Sum of forces in x-direction, Fx=F1-F2-Fg

+//1ton=2000 lbf

+Fg=(F1-F2-Fx)/2000 //tons

+printf("\n The force on the sluice gate = %.1f tons.\n",Fg)
\ No newline at end of file