initial commit / add all books

27 files changed, 298 insertions, 0 deletions

diff --git a/2159/CH1/EX1.1/1.sce b/2159/CH1/EX1.1/1.sce
new file mode 100755
index 000000000..38f589de9
--- /dev/null
+++ b/2159/CH1/EX1.1/1.sce
@@ -0,0 +1,9 @@
+//problem 1.1

+p=343350 //pressure at any point in pa

+w=9810 //gravitational constant

+h1=(p/w)

+disp(h1 ,"pressure in term of height of water(m)")

+s1=1

+s2=13.6

+h2=h1*s1/s2

+disp(h2 ,"pressure in term of height of mercury(m)")

diff --git a/2159/CH1/EX1.10/10.sce b/2159/CH1/EX1.10/10.sce
new file mode 100755
index 000000000..9700e7b5b
--- /dev/null
+++ b/2159/CH1/EX1.10/10.sce
@@ -0,0 +1,6 @@
+// problem 1.10

+s1=1.2

+s2=1

+s3=0.7

+h=(s1-s2)*0.3/(s2-s3)

+disp(h*100,"difference in height(cm)")

diff --git a/2159/CH1/EX1.11/11.sce b/2159/CH1/EX1.11/11.sce
new file mode 100755
index 000000000..a8f79c93d
--- /dev/null
+++ b/2159/CH1/EX1.11/11.sce
@@ -0,0 +1,10 @@
+// problem 1.11

+s1=0.8

+s2=13.6

+z=0.02

+w=9810

+h2=0.2

+h1=0.1

+h=h2*s2-h1*s1+(z*h2*(s2-s1))

+p=h*w

+disp(p,"pressure of the oil in N/m2")

diff --git a/2159/CH1/EX1.12/12.sce b/2159/CH1/EX1.12/12.sce
new file mode 100755
index 000000000..95553a7bc
--- /dev/null
+++ b/2159/CH1/EX1.12/12.sce
@@ -0,0 +1,12 @@
+//problem 1.12

+l=4

+b=2

+h=3

+w=9810

+s=0.8

+p1=w*l*b*h*s

+p2=w*s*l*h*1.5

+p3=w*s*b*h*1.5

+disp(p1,"total pressure on horizontal base")

+disp(p2, "total pressure on larger vertical base")

+disp(p3,"total pressure on smaller vertical walls")

diff --git a/2159/CH1/EX1.13/13.sce b/2159/CH1/EX1.13/13.sce
new file mode 100755
index 000000000..7b97d7346
--- /dev/null
+++ b/2159/CH1/EX1.13/13.sce
@@ -0,0 +1,10 @@
+//problem 1.13

+p=490500

+w=9810

+h=p/w

+D=0.15

+A=3.142*D*D*0.25

+pt=w*A*h

+h1=(D*D)/(16*h)

+disp(pt,"total hydrostatic pressure in N")

+disp(h1,"position of centre of pressure below the centre of pipe")

diff --git a/2159/CH1/EX1.14/14.sce b/2159/CH1/EX1.14/14.sce
new file mode 100755
index 000000000..373957ac1
--- /dev/null
+++ b/2159/CH1/EX1.14/14.sce
@@ -0,0 +1,9 @@
+//problem 1.14

+w=9810

+h=4

+d=2

+a=d*d*0.25*3.142

+p=w*a*h

+h1=d*d/64

+T=p*(h1) 

+disp(T, "torque required to keep the disc in vertical position in N.m")

diff --git a/2159/CH1/EX1.15/15.sce b/2159/CH1/EX1.15/15.sce
new file mode 100755
index 000000000..5b8b07273
--- /dev/null
+++ b/2159/CH1/EX1.15/15.sce
@@ -0,0 +1,10 @@
+//problem 1.15

+w=9810

+h=2

+l=2

+b=1

+a=l*b

+p=w*a*h

+h1=h+(b*l*l*l/(12*b*l*l))

+disp(p,"total pressure")

+disp(h1,"location of its centre of pressure")

diff --git a/2159/CH1/EX1.16/16.sce b/2159/CH1/EX1.16/16.sce
new file mode 100755
index 000000000..ef63bb3f4
--- /dev/null
+++ b/2159/CH1/EX1.16/16.sce
@@ -0,0 +1,13 @@
+//problem 1.16

+h1=8

+w=9810

+wd=6

+p1=0.5*w*h1*h1*wd

+h2=4

+p2=0.5*h2*w*h2*wd

+h11=0.66666*h1

+h22=0.6666*h2

+p=p1-p2

+hf=(p1*(h1-h11)-p2*(h2-h22))/p

+disp(p,"resultant force")

+disp(hf,"position of its line of action")

diff --git a/2159/CH1/EX1.17/17.sce b/2159/CH1/EX1.17/17.sce
new file mode 100755
index 000000000..924f2716d
--- /dev/null
+++ b/2159/CH1/EX1.17/17.sce
@@ -0,0 +1,8 @@
+//problem 1.17

+z=9810

+w=10

+h=2

+p=0.5*h*h*w*z

+h1=h*0.6666

+disp(p,"total hydrostatic thrust")

+disp(h1,"its point of application")

diff --git a/2159/CH1/EX1.18/18.sce b/2159/CH1/EX1.18/18.sce
new file mode 100755
index 000000000..cc60fb70e
--- /dev/null
+++ b/2159/CH1/EX1.18/18.sce
@@ -0,0 +1,14 @@
+//problem 1.18

+a1=1.4*2.2*1.4

+x1=1.6+0.7

+x11=(1.4*1.4/(12*2.3))+x1

+x2=0.7

+x22=(1.4*1.4/(12*0.7))+x2

+z=9810

+p1=z*a1*x1

+p2=z*a1*x2/1.4

+p=p1-p2

+h=(p1*(3-x11)+p2*(1.4-x22))/p

+f=(p1*(3-x11)-p2*(1.4-x22))/1.4

+disp(p,"resultant force")

+disp(f,"force acting horizontally on the top of the gate")

diff --git a/2159/CH1/EX1.19/19.sce b/2159/CH1/EX1.19/19.sce
new file mode 100755
index 000000000..fa4d39245
--- /dev/null
+++ b/2159/CH1/EX1.19/19.sce
@@ -0,0 +1,15 @@
+//problem 1.19

+s=1.5

+s1=0.9

+w=9810

+h1=0.9

+h2=0.6

+p1=0.5*w*s*s1*h1*h1 //total pressure due to oil

+p2=w*h1*h2*s*s1 // total pressure due to oil above water

+p3=w*h2*h2*0.5*s //total pressure due to water

+p=p1+p2+p3

+h=((p1*0.6666*h1)+(p2*(h1+0.5*h2))+(p3*(0.6666*h2+h1)))/p

+disp(p,"resultant pressure on the wall in N/m2")

+disp(h,"position of centre of pressure from free surface")

+

+

diff --git a/2159/CH1/EX1.2/2.sce b/2159/CH1/EX1.2/2.sce
new file mode 100755
index 000000000..a3a4353fc
--- /dev/null
+++ b/2159/CH1/EX1.2/2.sce
@@ -0,0 +1,12 @@
+//Problem 1.2

+h1=0.75 //atm pressure in term of mercury

+w=9810

+w1=13.6*w //specific weight of mercury

+Patm=w1*h1

+w2=15000

+h2=3 //

+p=w2*h2 // gauge pressure

+Pabs=Patm+p

+disp(p,"gauge pressure(N/m2)")

+disp(Pabs,"absolute pressure(N/m2)")

+

diff --git a/2159/CH1/EX1.20/20.sce b/2159/CH1/EX1.20/20.sce
new file mode 100755
index 000000000..d0efcc38e
--- /dev/null
+++ b/2159/CH1/EX1.20/20.sce
@@ -0,0 +1,10 @@
+// problem 1.20

+d=2.4

+h=1.6

+s=1.2

+a=d*s

+w=9810

+p=w*a*h*s

+h1=((2*s*s*s*d)/(12*a*h))+h

+disp(p,"total pressure in N")

+disp(h1,"its point of application")

diff --git a/2159/CH1/EX1.21/21.sce b/2159/CH1/EX1.21/21.sce
new file mode 100755
index 000000000..59e297b00
--- /dev/null
+++ b/2159/CH1/EX1.21/21.sce
@@ -0,0 +1,14 @@
+// problem 1.21

+x=30

+d=1.2

+h=1.5

+w=9810

+z=sin(x*3.142/180)

+h1=(z*d*0.5)+h

+a=0.25*3.142*d*d

+p=a*w*h1

+h11=(d*d*z*z)/(16*h1)+h1

+disp(p,"total pressure")

+disp(h11,"position of centre of pressure")

+

+

diff --git a/2159/CH1/EX1.22/22.sce b/2159/CH1/EX1.22/22.sce
new file mode 100755
index 000000000..99850bc53
--- /dev/null
+++ b/2159/CH1/EX1.22/22.sce
@@ -0,0 +1,10 @@
+// problem 1.22

+d=2

+z=0.5

+h=z+0.5*d

+w=9810

+a=3.142*d*d*0.25

+p=a*w*h

+h11=(1/(16*1.5))+1.5

+disp(p,"total pressure on the plate")

+disp(h11,"position of centre of pressure")

diff --git a/2159/CH1/EX1.23/23.sce b/2159/CH1/EX1.23/23.sce
new file mode 100755
index 000000000..18ffef8d3
--- /dev/null
+++ b/2159/CH1/EX1.23/23.sce
@@ -0,0 +1,12 @@
+// problem 1.23

+x=30

+z=sin(x*3.142/180)

+h=6-(z*0.5)

+l=1

+b=4

+a=l*b

+w=9810

+p=w*a*h

+h11=(z*z)/(12*h)+h

+f=p*0.5072

+disp(f,"force normal to the gate at point B")

diff --git a/2159/CH1/EX1.24/24.sce b/2159/CH1/EX1.24/24.sce
new file mode 100755
index 000000000..7937a828c
--- /dev/null
+++ b/2159/CH1/EX1.24/24.sce
@@ -0,0 +1,13 @@
+// problem 1.24

+x=30

+z=sin(3.142*x/180)

+d=1.4

+h=3

+b=1.5

+h1=z+d

+a=0.5*h*b

+w=9810

+p=w*a*h1

+h11=((z*z*h*h*h*b)/(36*a*h1))+h1

+disp(p,"total pressure on the plate")

+disp(h11,"position of centre of pressure")

diff --git a/2159/CH1/EX1.25/25.sce b/2159/CH1/EX1.25/25.sce
new file mode 100755
index 000000000..218065746
--- /dev/null
+++ b/2159/CH1/EX1.25/25.sce
@@ -0,0 +1,19 @@
+// problem 1.25

+d=1.8

+h=2.4

+w=9810

+s=0.8

+p1=w*d*d*h*0.25*3.142

+h1=((d*d)/(16*h))+h

+p=w*(s*1.5+2.4)

+p2=p*3.142*d*d*0.25

+p=p2-p1

+ab=w*(s*1.5+1.5)

+de=w*(s*1.5+3.3)

+ce=de-ab

+x=((0.5*ce*d*0.3)/(0.5*(ab+de)*d))

+h2=x+h

+h12=h1-h2

+disp(p, "change in total pressure")

+disp(h2,"position of centre of pressure")

+disp(h12,"change in position of centre of pressure")

diff --git a/2159/CH1/EX1.26/26.sce b/2159/CH1/EX1.26/26.sce
new file mode 100755
index 000000000..8cbed9044
--- /dev/null
+++ b/2159/CH1/EX1.26/26.sce
@@ -0,0 +1,13 @@
+// problem 1.26

+l=5

+r=3

+a=l*r

+h=r*0.5

+w=9810

+ph=w*a*h

+pv=w*0.25*3.142*r*r*l

+p=sqrt((ph*ph)+(pv*pv))

+z=ph/pv

+theta=atand(z)

+disp(p,"resultant pressure on the gate")

+disp(theta,"angle of resultant force with vertical")

diff --git a/2159/CH1/EX1.27/27.sce b/2159/CH1/EX1.27/27.sce
new file mode 100755
index 000000000..6188b1edd
--- /dev/null
+++ b/2159/CH1/EX1.27/27.sce
@@ -0,0 +1,10 @@
+// problem 1.27

+s=5

+z=sind(45)

+a=2*s*z

+h=s*z

+w=9810

+ph=w*a*h

+pv=w*((0.25*s*s*3.142)-(0.5*a*h))

+disp(ph,"horizontal pressure")

+disp(pv,"vertical pressure")

diff --git a/2159/CH1/EX1.3/3.sce b/2159/CH1/EX1.3/3.sce
new file mode 100755
index 000000000..7e55d494b
--- /dev/null
+++ b/2159/CH1/EX1.3/3.sce
@@ -0,0 +1,11 @@
+//Problem 1.3

+h1=2.5

+h2=1.5

+s1=1

+s2=0.8

+w=9810

+p2=s2*w*h2 //Pressure intensity at interface

+p1=s1*w*h1

+p=p1+p2

+disp(p2, "pressure intensity at interface(N/m2)")

+disp(p, "pressure intensity at bottom(N/m2)")

diff --git a/2159/CH1/EX1.4/4.sce b/2159/CH1/EX1.4/4.sce
new file mode 100755
index 000000000..1b8231281
--- /dev/null
+++ b/2159/CH1/EX1.4/4.sce
@@ -0,0 +1,7 @@
+//problem 1.4

+p=71613 //gauge pressure

+w=9810

+phead=p/w

+patm=10.33

+pabs=patm+phead

+disp(pabs,"absolute pressure in term of water height in meters")

diff --git a/2159/CH1/EX1.5/5.sce b/2159/CH1/EX1.5/5.sce
new file mode 100755
index 000000000..39a6afb4b
--- /dev/null
+++ b/2159/CH1/EX1.5/5.sce
@@ -0,0 +1,10 @@
+//problem 1.5

+h1=0.05

+h2=0.1

+s1=0.8

+s2=13.6

+w=9810

+p=s2*h2*w //pressure at balance line 

+p1=s1*h1*w 

+pf=p-p1

+disp(pf, "pressure in pipe(N/m2)")

diff --git a/2159/CH1/EX1.6/6.sce b/2159/CH1/EX1.6/6.sce
new file mode 100755
index 000000000..d654f62eb
--- /dev/null
+++ b/2159/CH1/EX1.6/6.sce
@@ -0,0 +1,9 @@
+//problem 1.6

+h1=0.2

+h2=0.5

+s1=0.9

+s2=13.6

+h=-(h1*s1+h2*s2)

+w=9810

+p=h*w

+disp(p, "vacuum pressure (N/m2)")

diff --git a/2159/CH1/EX1.7/7.sce b/2159/CH1/EX1.7/7.sce
new file mode 100755
index 000000000..874b4cb8a
--- /dev/null
+++ b/2159/CH1/EX1.7/7.sce
@@ -0,0 +1,9 @@
+//problem 1.7

+s1=0.8

+s2=13.6

+dh=0.4

+h=dh*13.6-dh*0.8

+w=9810

+pd=w*h

+disp(h, "pressure difference in height of water()")

+disp(pd, "presuure difference in N/m2")

diff --git a/2159/CH1/EX1.8/8.sce b/2159/CH1/EX1.8/8.sce
new file mode 100755
index 000000000..1cf6034b5
--- /dev/null
+++ b/2159/CH1/EX1.8/8.sce
@@ -0,0 +1,12 @@
+//problem 1.8

+s1=0.8

+s2=0.7

+h1=1.5

+h2=0.3

+h3=0.7

+s3=13.6

+hd=h2*s2+h3*s3-h1*s1

+w=9810

+pd=hd*w

+disp(hd, "diffrence in pressure in term of height of water(m)")

+disp(pd, "difference in pressure (N/m2)")

diff --git a/2159/CH1/EX1.9/9.sce b/2159/CH1/EX1.9/9.sce
new file mode 100755
index 000000000..b4dac5b03
--- /dev/null
+++ b/2159/CH1/EX1.9/9.sce
@@ -0,0 +1,11 @@
+//problem 1.9

+s1=1.6

+s2=0.8

+s3=13.6

+p1=98100

+p2=176580

+w=9810

+h1=p1/w

+h2=p2/w

+h=(h2-h1+1.6*s2-4.1*s1)/(s3-s2)

+disp(h.*100 ,"difference in mercury level(cm)")