1 files changed, 52 insertions, 0 deletions
diff --git a/3772/CH14/EX14.12/Ex14_12.sce b/3772/CH14/EX14.12/Ex14_12.sce
new file mode 100644
index 000000000..a3c5c5884
--- /dev/null
+++ b/3772/CH14/EX14.12/Ex14_12.sce
@@ -0,0 +1,52 @@
+// Problem no 14.12,Page No.340
+
+clc;clear;
+close;
+H=6 //m //height of dam
+a=1.5 //m //top width
+b=3.5 //m //Bottom width
+rho_s=16 //KN/m**3 //density of soil
+rho_mason=22.5 //KN/m**3 //density of mason
+phi=30 //Degree //angle of repose
+
+//Calculations
+
+//Let Length of dam ,L=1 m
+L=1 //m
+
+//weight of dam
+W=(a+b)*2**-1*L*H*rho_mason
+
+//Rankine's coeff earth pressure
+K=((1-sin(30*%pi*180**-1))*(1+sin(phi*%pi*180**-1))**-1)
+
+//Lateral thrust
+P=rho_s*H**2*L*2**-1*K
+
+//Distance of Line of action from vertical base
+x_bar=(b**2+b*a+a**2)*(3*(b+a))**-1
+
+//distance of pt where resultant cuts the base
+x=P*W**-1*H*3**-1
+
+//Eccentricity
+e=x_bar+x-b*2**-1
+
+//Pressure at heel B
+sigma1=W*b**-1*(1-6*e*b**-1)
+
+//Pressure at heel C
+sigma2=W*b**-1*(1+6*e*b**-1)
+
+//Result
+printf("The Max Intensities of soil at the wall is %.2f",sigma2);printf(" KN/m**2")
+printf("\n The Min Intensities of soil at the wall is %.2f",sigma1);printf(" KN/m**2")
+
+//Plotting the Shear Force Diagram
+X1=[0,L,L]
+Y1=[sigma2,sigma1,0]
+Z1=[0,0,0]
+plot(X1,Y1,X1,Z1)
+xlabel("Length x in m")
+ylabel("Shear Force in kN")
+title("the Shear Force Diagram")