diff options
Diffstat (limited to '2087/CH10')
| -rwxr-xr-x | 2087/CH10/EX10.1/example10_1.sce | 33 | ||||
| -rwxr-xr-x | 2087/CH10/EX10.2/example10_2.sce | 16 | ||||
| -rwxr-xr-x | 2087/CH10/EX10.3/example10_3.sce | 25 | ||||
| -rwxr-xr-x | 2087/CH10/EX10.4/example10_4.sce | 98 | ||||
| -rwxr-xr-x | 2087/CH10/EX10.5/example10_5.jpg | bin | 0 -> 24475 bytes | |||
| -rwxr-xr-x | 2087/CH10/EX10.5/example10_5.sce | 34 |
6 files changed, 206 insertions, 0 deletions
diff --git a/2087/CH10/EX10.1/example10_1.sce b/2087/CH10/EX10.1/example10_1.sce new file mode 100755 index 000000000..216e16f07 --- /dev/null +++ b/2087/CH10/EX10.1/example10_1.sce @@ -0,0 +1,33 @@ +
+
+//example 10.1
+//calculate seepage flow per unit length of dam
+clc;funcprot(0);
+//given
+K=5D-4; //coefficient of permeability of soil
+Bt=6; //width of top of dam
+wb=146; //width of base of dam
+H=20; //heigth of dam
+hw=2; //heigth of water in reservior
+hs1=4; //slope on upstream side
+hs2=3; //slope on downstream side
+df=30; //length of drainage filter
+
+x=wb-df-72+72*0.3;
+y=18;
+s=(x^2+y^2)^0.5-x;
+
+x=[0 10 20 30 40 50 60 65.6];
+for i=1:8
+ y(i)=(4.849*x(i)+5.879)^0.5;
+ y(i)=round(y(i)*1000)/1000;
+end
+
+mprintf("\nx y");
+for i=1:8
+ mprintf("\n%f %f",x(i),y(i));
+end
+
+sf=K*s*10000;
+sf=round(sf*1000)/1000;
+mprintf("\nSeepage flow per unit length of dam=%fD-6 cumecs/metre length of dam.",sf);
diff --git a/2087/CH10/EX10.2/example10_2.sce b/2087/CH10/EX10.2/example10_2.sce new file mode 100755 index 000000000..b9bd4a9f0 --- /dev/null +++ b/2087/CH10/EX10.2/example10_2.sce @@ -0,0 +1,16 @@ +
+
+//example 10.2
+//calculate discharge per meter length of dam
+clc;funcprot(0);
+//given
+K=3D-3; //coefficient of permeability
+nd=25; //number of potential drops
+nf=4; //number of flow channels
+lf=40; //filter length
+H=52; //heigth of dam
+fb=2; //free board
+
+q=K*(H-fb)*nf/(nd*100);
+mprintf("Discharge per meter length of dam=%f cumec/metre length.",q);
+
diff --git a/2087/CH10/EX10.3/example10_3.sce b/2087/CH10/EX10.3/example10_3.sce new file mode 100755 index 000000000..e0a24bb35 --- /dev/null +++ b/2087/CH10/EX10.3/example10_3.sce @@ -0,0 +1,25 @@ +
+
+//example10.3
+//calculate factor of safety for slope
+clc;funcprot(0);
+//given
+x=4; //given scale
+An=14.4; //area of N rectangle
+At=6.4; //area of T rectangle
+Au=4.9; //area of U rectangle
+L=12.6; //length of arc;
+gamma_m=19; //unit weigth of soil
+gamma_w=9.81; //unit weigth of water
+fi=26; //effective angle(degree)
+co=19.5; //cohesion value
+
+//consider 1m length of dam
+SumN=An*x^2*gamma_m;
+SumT=At*x^2*gamma_m;
+SumU=Au*x^2*gamma_w;
+Le=x*L;
+F=((Le*co)+(SumN-SumU)*tand(fi))/SumT;
+F=round(F*100)/100;
+mprintf("Factor of safety for slope=%f.",F);
+
diff --git a/2087/CH10/EX10.4/example10_4.sce b/2087/CH10/EX10.4/example10_4.sce new file mode 100755 index 000000000..f802781fb --- /dev/null +++ b/2087/CH10/EX10.4/example10_4.sce @@ -0,0 +1,98 @@ +
+
+//example 10.4
+//check section for:
+//Stability of d/s slope against steady seepage
+//Sloughing of u/s slope against sudden drawdown
+//Stability of the foundation against shear
+//Seepage through body of dam
+clc;funcprot(0);
+//given
+//Dimensions
+H=20; //Heigth of dam
+Bt=6; //top width of dam
+s1=4; //u/s slope
+s2=3; //d/s slope
+fb=2; //free board
+//Properties of materials of dam
+gamma_d=17.27; //dry density
+wc=0.15; //optimum water content
+gamma_s=21.19; //saturated density
+gamma_w=9.81; //unit weigth of water
+wavg=19.62; //average unit weigth under seepage
+theta=26; //average angle of internal friction(degree)
+co=19.13; //average cohesion
+K=5D-4; //coefficient of permeability
+//properties of foundation materials
+gamma_f=17.27; //average unit weigth
+cof=47.87; //average cohesion
+fi=8; //average angle internal friction
+t=6; //thickness of clay
+FOSp=1.5; //permissible factor of safety of slope
+PS=8D-6; //permissible seepage
+
+
+//(a) Stability of d/s slope against steady seepage
+An=302.4; //area of N diagram
+At=91.2; //area of T diagram
+Au=98.4; //area of U diagram
+Le=60; //length of arc
+SumN=An*gamma_s;
+SumT=At*gamma_s;
+SumU=Au*gamma_w;
+F=((Le*co)+(SumN-SumU)*tand(theta))/SumT;
+F=round(F*100)/100;
+mprintf("Part(a):")
+mprintf("\nFactor of safety for slope=%f.",F);
+mprintf("\nSafe");
+
+//(b) Sloughing of u/s slope against sudden drawdown
+h1=15;
+b=80;
+P=gamma_s*H^2*tand(45-(theta/2))^2/2+gamma_w*h1^2/2;
+sav=P/b;
+smax=2*sav;
+Ne=(gamma_s-gamma_w)*b*H/2;
+R=Ne*tand(theta)+co*b;
+fs=R/P;
+fs=round(fs*100)/100;
+mprintf("\n\nPart(b):")
+mprintf("\nFactor of safety w.r.t average shear=%f.",fs);
+mprintf("\nSafe");
+sr=0.6*H*(gamma_s-gamma_w)*tand(theta)+co;
+FS=sr/smax;
+FS=round(FS*100)/100;
+mprintf("\n\nFactor of safety w.r.t maximum shear=%f.",FS);
+mprintf("\nSafe");
+
+//(c) Stability of the foundation against shear
+h1=26;
+h2=6;
+gamma_m=(wavg*(h1-h2)+gamma_f*h2)/h1;
+l=(gamma_m*h1*tand(fi)+cof)/(gamma_m*h1);
+fi1=atand(l);
+P=(h1^2-h2^2)/2*gamma_m*tand(45-(fi1/2))^2;
+sav=P/b;
+smax=2*sav;
+s1=cof+gamma_f*h2*tand(fi);
+s2=cof+gamma_m*h1*tand(fi);
+as=(s1+s2)/2;
+fs=as/sav;
+fs=round(fs*100)/100;
+mprintf("\n\nPart(c):")
+mprintf("\nFactor of safety w.r.t overall shear=%f.",fs);
+mprintf("\nSafe");
+
+gamma_av=(wavg*0.6*H+gamma_f*h2)/((0.6*H)+h2);
+s=cof+gamma_av*0.6*H*tand(fi);
+fs=s/smax;
+fs=round(fs*100)/100;
+mprintf("\n\nFactor of safety w.r.t overall shear=%f.",fs);
+mprintf("\nUnsafe");
+
+//(d) Seepage through body of dam
+s=2; //measured
+q=K*s*100000/100;
+mprintf("\n\nPart(d):")
+mprintf("\n Seepage through body of dam=%fD-5 cumecs/m length of dam",q);
+
diff --git a/2087/CH10/EX10.5/example10_5.jpg b/2087/CH10/EX10.5/example10_5.jpg Binary files differnew file mode 100755 index 000000000..a117603d4 --- /dev/null +++ b/2087/CH10/EX10.5/example10_5.jpg diff --git a/2087/CH10/EX10.5/example10_5.sce b/2087/CH10/EX10.5/example10_5.sce new file mode 100755 index 000000000..0a4aa3aae --- /dev/null +++ b/2087/CH10/EX10.5/example10_5.sce @@ -0,0 +1,34 @@ +
+
+//example 10.5
+//design upstream impervious blanket
+clc;funcprot(0);
+//given
+Zb=1.2; //thickness of blanket
+Zf=8; //distance of blanket from foundation
+kb=0.06; //coefficient of permeability of blanket material
+kf=72; //coefficient of permeability of foundation soil
+Hw=10; //heigth of water in reservior
+Xd=40;
+
+a=(kb/(kf*Zb*Zf))^0.5;
+Xo=1.414/a;
+
+//we vary value of x
+x=[0 25 50 75 100 125 151.8 300]
+for i=1:8
+ e=exp(2*a*x(i));
+ Xr(i)=(e-1)/(a*(e+1));
+ ho(i)=Xr(i)*Hw/(Xr(i)+Xd);
+ r(i)=Xr(i)*100/(Xr(i)+Xd);
+end
+mprintf("\nx Xr ho reduction q(percent)");
+for i=1:8
+ mprintf("\n%f %f %f %f",x(i),Xr(i),ho(i),r(i));
+end
+//graph is plotted between r and x.
+//after around 130m length there is only slight increase in head dissipated(ho)
+L=130;
+mprintf("\nThickness of blanket=%f m",Zb);
+mprintf("\nLength of blanket=%i m.",L);
+
|
