initial commit / add all books

7 files changed, 121 insertions, 0 deletions

diff --git a/746/CH3/EX3.01/3_01.jpg b/746/CH3/EX3.01/3_01.jpg
new file mode 100755
index 000000000..74d8becc0
--- /dev/null
+++ b/746/CH3/EX3.01/3_01.jpg
diff --git a/746/CH3/EX3.01/3_01.sce b/746/CH3/EX3.01/3_01.sce
new file mode 100755
index 000000000..089d4a614
--- /dev/null
+++ b/746/CH3/EX3.01/3_01.sce
@@ -0,0 +1,21 @@
+//liquid level//

+pathname=get_absolute_file_path('3.01.sce')

+filename=pathname+filesep()+'3.01-data.sci'

+exec(filename)

+//Tube diameter(in mm):

+D=1:25;

+D1=D/1000

+[m n]=size(D1)

+for i=1:n

+//Change in liquid level for water(in mm):

+dhw(i)=4*STw*cosd(thetaw)/dw/g/D1(i);

+//Change in liquid level for mercury(in mm):

+dhm(i)=4*STm*cosd(thetam)/dm/g/D1(i);

+end;

+

+//Plotting tube daimeter and water level:

+plot(D1*1000,dhw,'-o')

+//Plotting tube daimeter and mercury level:

+plot(D1*1000,dhm,'-*')

+legend(['Water';'Mercury']);

+xtitle('Liquid level vs Tube diameter','Liquid level(in mm)','Tube diameter(in mm)')

diff --git a/746/CH3/EX3.03/3_03.sce b/746/CH3/EX3.03/3_03.sce
new file mode 100755
index 000000000..92dc4e6bd
--- /dev/null
+++ b/746/CH3/EX3.03/3_03.sce
@@ -0,0 +1,8 @@
+//pressure difference//

+pathname=get_absolute_file_path('3.03.sce')

+filename=pathname+filesep()+'3.03-data.sci'

+exec(filename)

+//Pressure difference(in lbf/in^2):

+dp=g*d*(-d1+SGm*d2-SGo*d3+SGm*d4+d5)/12/144

+printf("\n\nRESULTS\n\n")

+printf("\n\nPressure difference between A and B: %f lbf/in^2\n\n",dp)

diff --git a/746/CH3/EX3.04/3_04.sce b/746/CH3/EX3.04/3_04.sce
new file mode 100755
index 000000000..4030fcd51
--- /dev/null
+++ b/746/CH3/EX3.04/3_04.sce
@@ -0,0 +1,40 @@
+//temperature and pressure//

+pathname=get_absolute_file_path('3.04.sce')

+filename=pathname+filesep()+'3.04-data.sci'

+exec(filename)

+//Assuming temperature varies linearly with altitude:

+//Temperature gradient(in F/ft):

+m=(T1-T2)/(z2-z1)

+//Value of g/(m*R):

+v=g/m/R/32.2

+//Pressure at Vail Pass(in inches of Hg):

+p12=p1*((T2+460)/(T1+460))^v

+//Percentage change in density:

+pc1=(p12/p1*(T1+460)/(T2+460)-1)*100

+//Assuming density is constant:

+//Pressure at Vail Pass(in inches of Hg):

+p22=p1*(1-(g*(z2-z1)/(R*32.2)/(T1+460)))

+//Percentage change in density:

+pc2=0;

+//Assuming temperature is constant:

+//Pressure at Vail Pass(in inches of Hg):

+p32=p1*%e^(-g*(z2-z1)/(R*32.2)/(T2+460))

+//Percentage change in density:

+pc3=(p32/p1*(T1+460)/(T1+460)-1)*100

+//For an adiabatic atmosphere:

+p42=p1*((62+460)/(80+460))^(k/(k-1))

+//Percentage change in density:

+pc4=(p42/p1*(T1+460)/(T2+460)-1)*100

+printf("\n\nRESULTS\n\n")

+printf("\n\n1) If temperature varies linearly with altitude\n\n")

+printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p12)

+printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc1)

+printf("\n\n2) If density is constant\n\n")

+printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p22)

+printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc2)

+printf("\n\n3) If temperature is constant\n\n")

+printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p32)

+printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc3)

+printf("\n\n4) For an adiabatic atmosphere\n\n")

+printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p42)

+printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc4)

diff --git a/746/CH3/EX3.05/3_05.sce b/746/CH3/EX3.05/3_05.sce
new file mode 100755
index 000000000..4cea618b9
--- /dev/null
+++ b/746/CH3/EX3.05/3_05.sce
@@ -0,0 +1,25 @@
+//force and pressure//

+pathname=get_absolute_file_path('3.05.sce')

+filename=pathname+filesep()+'3.05-data.sci'

+exec(filename)

+//Net force on the gate(in kN):

+Fr=d*g*w*(D*L+L^2/2*sind(theta))

+//Centre of pressure:

+//Calculation for y coordinate:

+   yc=D/sind(theta)+L/2

+   //Area(in m^2):

+   A=L*w

+   //Moment of inertia of rectangular gate(in m^4):

+   Ixx=w*L^3/12

+   //y coordinate(in m):

+   y=yc+Ixx/A/yc

+//Calculation for x coordinate:

+   Ixy=0

+   xc=w/2

+   //x coordinate(in m):

+   x=xc+Ixy/A/xc

+printf("\n\nRESULTS\n\n")

+printf("\n\nNet force on the gate: %f kN\n\n",Fr/1000)

+printf("\n\nCoordinate of centre of pressure:(%0.1f,%0.1f)\n\n",x,y)

+       

+      

diff --git a/746/CH3/EX3.06/3_06.sce b/746/CH3/EX3.06/3_06.sce
new file mode 100755
index 000000000..89f2daf73
--- /dev/null
+++ b/746/CH3/EX3.06/3_06.sce
@@ -0,0 +1,9 @@
+//force//

+pathname=get_absolute_file_path('3.06.sce')

+filename=pathname+filesep()+'3.06-data.sci'

+exec(filename)

+//Force required to keep the door shut(in lbf):

+function y=f(z),y=b/L*p0*z+d*b/L*(L*z-z^2),endfunction

+Ft=intg(0,L,f)

+printf("\n\nRESULTS\n\n")

+printf("\n\nForce requiredto kep the door shut: %.1f lbf\n\n",Ft)

diff --git a/746/CH3/EX3.07/3_07.sce b/746/CH3/EX3.07/3_07.sce
new file mode 100755
index 000000000..8624fca4a
--- /dev/null
+++ b/746/CH3/EX3.07/3_07.sce
@@ -0,0 +1,18 @@
+//force at equilibrium//

+pathname=get_absolute_file_path('3.07.sce')

+filename=pathname+filesep()+'3.07-data.sci'

+exec(filename)

+//Horizontal component of resultant force(in kN):

+Frh=0.5*d*g*w*D^2

+//Line of action of Frh(in m):

+y1=0.5*D+w*D^3/12/(0.5*D)/(w*D)

+//Vertical component of resultant force(in kN):

+function y=q(x), y=d*g*w*(D-sqrt(a*x)),endfunction

+Frv=intg(0,D^2/a,q)

+//Line of acion of Frv(in m):

+function k=f(x), k=d*g*w/Frv*x*(D-sqrt(a*x)),endfunction

+xa=intg(0,D^2/a,f)

+//Force required to keep the gate in equilibrium(in kN):

+Fa=1/l*(xa*Frv+(D-y1)*Frh)

+printf("\n\nRESULTS\n\n")

+printf("\n\nForce required to keep the gate at equilibrium: %f kN\n\n",Fa/1000)