11 files changed, 220 insertions, 0 deletions
diff --git a/2732/CH5/EX5.1/Ex5_1.sce b/2732/CH5/EX5.1/Ex5_1.sce
new file mode 100755
index 000000000..3ae30b5a5
--- /dev/null
+++ b/2732/CH5/EX5.1/Ex5_1.sce
@@ -0,0 +1,11 @@
+clc
+//initialization of variables
+clear
+l=20 //cm
+dL=1 //m
+dl=0.004 //cm
+//calculations
+L=l*dL/dl //m
+//results
+printf('The depth of the clay bed is %d m',L)
+
diff --git a/2732/CH5/EX5.10/Ex5_10.sce b/2732/CH5/EX5.10/Ex5_10.sce
new file mode 100755
index 000000000..981db26e1
--- /dev/null
+++ b/2732/CH5/EX5.10/Ex5_10.sce
@@ -0,0 +1,41 @@
+clc
+// initialization
+clear
+d=10 //cm
+id=9.99 //cm
+t=3 //mm
+E=1.0*10^6 //kg/cm^2
+a=2.02*10^-5 // degree/celcius
+// part(a)
+Tr=10 //degree C
+T=(d-id)/id*1/a
+printf('part(a) \n The sleeve must be heated to %.1f degree C or more for this purpose',T+Tr)
+
+//part(b)
+s_th=a*T*E
+p=s_th*t*2/(d*10)
+printf('\n part(b) \n The pressure developed between the rod and sleeve is %d kg/cm^2',p)
+
+// part(c)
+f=0.2
+o=10 // overlap: cm
+A=%pi*d*o
+F=f*p*A
+printf('\n part (c) \n The axial force required is %d kg',F)
+
+//part (d)
+// linked to part c
+T2=20 //degree C
+a2=1.17*10^-5 //  /degree C
+Ts=(a-a2)*(T2-Tr)*E
+Ts=s_th-Ts
+p2=p*Ts/s_th
+F2=F*Ts/s_th
+printf('\n part(d)\n The pressure developed between the rod and sleeve is %.1f kg/cm^2',p2)
+printf('\n The axial force required is %d kg',F2)
+//part(e)
+T3=Tr+(s_th/((a-a2)*10^6))
+printf('\n part(e) \n The temperature at which the sleeve comes off easily is %.1f C',T3)
+
+// calculations in the text: rounding off errors
+
diff --git a/2732/CH5/EX5.11/Ex5_11.sce b/2732/CH5/EX5.11/Ex5_11.sce
new file mode 100755
index 000000000..9f7cab62a
--- /dev/null
+++ b/2732/CH5/EX5.11/Ex5_11.sce
@@ -0,0 +1,18 @@
+clc
+//initialization of variables
+clear
+T1=37.8 // degre C
+t=0.355 //mm
+T2=93.3 // degree C
+L=2 //cm
+m=1
+n=1.53
+a=1.86*10^-5
+//calculations
+R=2*t*(3*(1+m)^2+(1+m*n)*(m^2+(m*n)^-1))
+R=R/(6*a*(T2-T1)*(1+m^2))  // mm
+R=R/10
+def=L^2/(8*R)
+// results
+printf('The radius of curvature is %.1f cm',R)
+printf('\n The deflection is %.6f cm',def)
diff --git a/2732/CH5/EX5.12/Ex5_12.sce b/2732/CH5/EX5.12/Ex5_12.sce
new file mode 100755
index 000000000..d18002735
--- /dev/null
+++ b/2732/CH5/EX5.12/Ex5_12.sce
@@ -0,0 +1,19 @@
+clc
+// initialization of variables
+clear
+L=5 //cm
+D=1.8 //cm
+l=2.5 //cm
+d=1.5 //cm
+F=1 //tonne
+E=2.1*10^6 //kg/cm^2
+// calculations
+s1=F*1000*4/(D^2*%pi)
+s2=F*1000*4/(d^2*%pi)
+U1=1/2*s1^2/E
+U1=U1*L*D^2*%pi/4
+U2=1/2*s2^2/E
+U2=U2*l*d^2*%pi/4
+U=U1+U2
+// results
+printf('The energy stored in the bolt is %.3f kg-cm',U)
diff --git a/2732/CH5/EX5.13/Ex5_13.sce b/2732/CH5/EX5.13/Ex5_13.sce
new file mode 100755
index 000000000..cd52d987f
--- /dev/null
+++ b/2732/CH5/EX5.13/Ex5_13.sce
@@ -0,0 +1,31 @@
+clc
+// initialization of variables
+clear
+t=16 //mm
+Pt=1500 //kg/cm^2
+Ps=1025 //kg/cm^2
+Pb=2360 //kg/cm^2
+
+//part (a)
+p=6 //cm
+r=24 //mm
+d=r/10+0.15
+Ft=t*(p-d)*Pt/10
+Fs=%pi*d^2*Ps/4
+Fb=d*t*Pb
+x=min(Ft,Fs,Fb)
+effA=x*100/(p*t/10*Pt)
+
+//part (b)
+p=9 //cm
+r=30 //mm
+d=r/10+0.2
+Ft=t*(p-d)*Pt/10
+Fs=%pi*d^2*Ps/4
+Fb=d*t*Pb
+x=min(Ft,Fs,Fb)
+effB=x*100/(p*t/10*Pt)
+
+// results
+printf('The efficiencies corresponding to cases a and b are %.1f, %.1f',effA,effB)
+printf('\n Hence part b is better than part a')
diff --git a/2732/CH5/EX5.2/Ex5_2.sce b/2732/CH5/EX5.2/Ex5_2.sce
new file mode 100755
index 000000000..f54a3127a
--- /dev/null
+++ b/2732/CH5/EX5.2/Ex5_2.sce
@@ -0,0 +1,11 @@
+clc
+//initialization of variables
+clear
+A=1 //unit area
+E=2*10^6 //kg/cm^2
+// calculations
+db=3000*90/(A*E)
+dc=db+5000*60/(A*E)
+dd=dc+4000*30/(A*E)
+//results
+printf('The extension of the rod in part AB is %.2e cm in part BC is %.2e cm \n and in part CD is %.2e cm',db,dc,dd)
diff --git a/2732/CH5/EX5.3/Ex5_3.sce b/2732/CH5/EX5.3/Ex5_3.sce
new file mode 100755
index 000000000..381ecde7b
--- /dev/null
+++ b/2732/CH5/EX5.3/Ex5_3.sce
@@ -0,0 +1,11 @@
+clc
+//initialization of variables
+clear
+A=3 //cm^2
+L=18 //m
+E= 2*10^6 //kg/cm^2
+r=7833 //kg/m^3
+//calculations
+e=r*(L*100)^2/(2*E*10^6)
+// results
+printf('The elongation is %.5f cm',e)
diff --git a/2732/CH5/EX5.4/Ex5_4.sce b/2732/CH5/EX5.4/Ex5_4.sce
new file mode 100755
index 000000000..bdbe56189
--- /dev/null
+++ b/2732/CH5/EX5.4/Ex5_4.sce
@@ -0,0 +1,12 @@
+clc
+//initialization of variables
+clear
+// linked to 5_3
+P=3 //tonne
+E=2*10^6 //kg/cm^2
+d_0= 1 //cm
+d_l=2.8 //cm
+// calculations
+e=4*P*1000*d_l*10^3/(d_l^2*%pi*E*(1-((d_l-d_0)/d_l)))
+//results
+printf('The total elongation is %.2f cm',e)
diff --git a/2732/CH5/EX5.6/Ex5_6.sce b/2732/CH5/EX5.6/Ex5_6.sce
new file mode 100755
index 000000000..43f537c60
--- /dev/null
+++ b/2732/CH5/EX5.6/Ex5_6.sce
@@ -0,0 +1,20 @@
+clc
+//initialization of variables
+clear
+P=10 //tonne
+E=2*10^6 //kg/cm^2
+// calculations
+// We have to solve linear system Ax=B
+A=[1 1 1 0
+   3 1 -3 0
+   -2 2 0 -E
+   0 -1 2 -E]
+B=[P*10^3;0;0;0]
+x=inv(A)*B
+W1=x(1,1)/1000
+W2=x(2,1)/1000
+W3=x(3,1)/1000
+th=x(4,1)
+//results
+printf('The load taken by each rod is %.2f tonne, %.1f tonne, %.3f tonne',W1,W2,W3)
+printf('\n and the slope is theta = %.2e. radians',th)   
diff --git a/2732/CH5/EX5.8/Ex5_8.sce b/2732/CH5/EX5.8/Ex5_8.sce
new file mode 100755
index 000000000..75c977e50
--- /dev/null
+++ b/2732/CH5/EX5.8/Ex5_8.sce
@@ -0,0 +1,14 @@
+clc
+// initialization of variables
+clear
+b=30 // cm
+h=30 //cm
+n=6
+A=36 //cm^2
+ss_s=1500 //kg/cm^2
+ss_c=60 //kg/cm^2
+Er=15 // Elasticity ratio
+// calculations
+L=A*Er*ss_c+(b*h-A)*ss_c
+// results
+printf('The safe load is %d.kg',L)
diff --git a/2732/CH5/EX5.9/Ex5_9.sce b/2732/CH5/EX5.9/Ex5_9.sce
new file mode 100755
index 000000000..dd3570dd2
--- /dev/null
+++ b/2732/CH5/EX5.9/Ex5_9.sce
@@ -0,0 +1,32 @@
+clc
+// initiaization of variables
+clear
+gs_b=10 //cm
+gs_h=10 //cm
+d_b=2 //cm
+d_h=2 //cm
+As= 1 //cm^2
+s=10000 //kg/cm^2
+// part (a)
+Es=2*10^6 //kg/cm^2
+Ec=2*10^5 //kg/cm^2
+// calculations
+e=s/Es
+Ac=gs_b*gs_h-(d_b*d_h)
+e_c=e*Es*As/(Ec*Ac+Es*As)
+s_c=Ec*e_c
+e_s=e-e_c
+s_s=Es*e_s
+// results
+printf('part (a) \n The stress in steel and concrete are respectively %d , %.2e kg/cm^2',s_s,s_c)
+// part(b)
+P=8000 //kg
+// calculations
+e_c=(e*Es*As-P)/(Ec*Ac+Es*As)
+e_s=e-e_c
+s_c=Ec*e_c
+s_s=Es*e_s
+// results
+printf('\n part (b) \n The stress in steel and concrete are respectively %.1f , %.2f kg/cm^2',s_s,s_c)
+
+