1 files changed, 48 insertions, 0 deletions
diff --git a/2087/CH12/EX12.1/example12_1.sce b/2087/CH12/EX12.1/example12_1.sce
new file mode 100755
index 000000000..0a38fde25
--- /dev/null
+++ b/2087/CH12/EX12.1/example12_1.sce
@@ -0,0 +1,48 @@
+
+
+//example 12.1
+//calculate average hydraulic gradient
+//uplift presuures and thickness of floor at 6m, 12m and 18m from u/s
+clc;funcprot(0);
+//given
+rho=2.24;             //relative density of material
+gamma_w=9.81;         //unit weigth of water
+L=22;                 //total length
+lc=(2*6)+L+(2*8);     //length of creep
+hg=4/lc;              //hydraulic gradient
+mprintf("avearge hydraulic gradient=%f.",hg);
+//at 6 m from u/s
+x=6;
+lg=(6*2)+x;
+h1=4*(1-lg/50);       //unbalanced head
+up=gamma_w*h1;
+t=4*h1/(3*(rho-1));
+up=round(up*100)/100;
+t=round(t*100)/100;
+mprintf("\n\nuplift at 6 m from u/s=%f kN/square metre.",up);
+mprintf("\nthickness at 6 m from u/s=%f m.",t);
+
+//at 12 m from u/s
+x=12;
+lg=(6*2)+x;
+h1=4*(1-lg/50);       //unbalanced head
+up=gamma_w*h1;
+t=4*h1/(3*(rho-1));
+up=round(up*100)/100;
+t=round(t*100)/100;
+mprintf("\n\nuplift at 12 m from u/s=%f kN/square metre.",up);
+mprintf("\nthickness at 12 m from u/s=%f m.",t);
+
+//at 18m from u/s
+x=18;
+lg=(6*2)+x;
+h1=4*(1-lg/50);       //unbalanced head
+up=gamma_w*h1;
+t=4*h1/(3*(rho-1));
+up=round(up*10)/10;
+t=round(t*100)/100;
+mprintf("\n\nuplift at 18 m from u/s=%f kN/square metre.",up);
+mprintf("\nthickness at 18 m from u/s=%f m.",t);
+
+
+