From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 2087/CH12/EX12.2/example12_2.sce | 105 +++++++++++++++++++++++++++++++++++++++
 1 file changed, 105 insertions(+)
 create mode 100755 2087/CH12/EX12.2/example12_2.sce

(limited to '2087/CH12/EX12.2/example12_2.sce')

diff --git a/2087/CH12/EX12.2/example12_2.sce b/2087/CH12/EX12.2/example12_2.sce
new file mode 100755
index 000000000..b2956c85d
--- /dev/null
+++ b/2087/CH12/EX12.2/example12_2.sce
@@ -0,0 +1,105 @@
+
+
+//example 12.2
+//calculate uplift pressure and exit gradient
+//check whether section is safe against overturning and piping
+clc;funcprot(0);
+//given
+b=54;            //width of section
+D1D2=16;         //distance between points D1 and D2
+D2D3=37;         //distance between points D2 and D3
+
+//first pipe line
+//taking data from figure
+d=105-97;
+b1=0.5;
+alpha=b/d;
+//from the curves we get
+fic1=0.665;
+fid1=0.76;
+fie1=1;
+t=105-104;                      //floor thickness
+corec=(fid1-fic1)*100*t/d;          //correction for floor thickness
+//for pile no. 2
+D=104-97;
+d=104-97;
+bdash=16;
+C=19*(D/bdash)^0.5*(d+D)/b;    //correction for pile no. 2
+fic1=fic1*100+corec+C;         //corrected pressures
+
+//intermedite pipe line
+d=105-97;
+b1=16.5;
+alpha=b/d;
+r=b1/b;               //ratio b1/b
+//from the curves we get
+fic2=0.52;
+fie2=0.725;
+fid2=0.615;
+corec_c1=(fid2-fic2)*100*t/d;
+corec_e1=(fie2-fid2)*100/d;
+
+//for pile no. 1
+C1=C;
+d=104-97;
+bdash=37;
+D=104-95;
+C2=19*(D/bdash)^0.5*(d+D)/b;
+//correction due to slope
+corec_e2=3.3;                //from table 12.4
+//correction is negative due to upwrd slope
+l=4;               //horizontal length of slope
+corec_c2=corec_e2*l/bdash;
+
+fie2=fie2*100-corec_e1-corec_e2;
+fic2=fic2*100+corec_c1+C2-corec_c2;
+
+//pile no. 3 at d/s end
+d=103.5-95;
+alpha_=d/b;
+//for curves
+fie3=0.35;fid3=0.242;
+corec_t=(fie3-fid3)*100*(103.5-102)/d;
+
+//correction for interference at pile no. 2
+d=102-95;
+D=102-97;
+C3=19*(D/bdash)^0.5*(d+D)/b;
+fie3=fie3*100-corec_t-C3;
+
+point=['C1' 'C2' 'E2' 'E3'];            //Point
+P=[fic1 fic2 fie2 fie3];                //pressure percent
+P_=[3.55 2.78 3.39 1.58];               //pressure head
+mprintf("Points         Pressure percent        Pressure head");
+for i=1:4
+    P(i)=round(P(i)*10)/10;
+    mprintf("\n%s                %f               %f",point(i),P(i),P_(i));
+end
+
+//check for floor thickness
+Pa=P_(2)-((P_(2)-P_(4))*6.5/37);
+Pb=P_(2)-((P_(2)-P_(4))*24/37);
+Pc=P_(2)-((P_(2)-P_(4))*30/37);
+rho=2.24;                              //specific gravity of concrete
+ta=Pa/(rho-1);
+tb=Pb/(rho-1);
+tc=Pc/(rho-1);
+ta=round(ta*100)/100;
+tb=round(tb*100)/100;
+tc=round(tc*100)/100;
+mprintf("\n\nThickness required at A=%f m.",ta);
+mprintf("\nThickness required at B=%f m.",tb);
+mprintf("\nThickness required at C=%f m.",tc);
+t=103.5-102;
+mprintf("\nThickness provided=%f m.",t);
+mprintf("\nFloor thickness at B and C are adequate");
+
+//exit gradient
+H=108.5-103.5;               //seepage head
+d=103.5-95;                 //depth cut-off
+//from exit gradient curve
+alpha=6.35;
+lambda=(1+(1+alpha^2)^0.5)/2;
+Ge=H/(d*%pi*lambda^0.5);
+mprintf("\n\nexit gradient=%f.",Ge);
+mprintf("\n it is less than permissible exit gradient  < 1/6\nHence safe..");
-- 
cgit