diff options
Diffstat (limited to 'Principles_Of_Geotechnical_Engineering_by_B_M_Das')
15 files changed, 3086 insertions, 0 deletions
diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/10-stress_in_a_soil_mass.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/10-stress_in_a_soil_mass.ipynb new file mode 100644 index 0000000..23b49a9 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/10-stress_in_a_soil_mass.ipynb @@ -0,0 +1,96 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 10: stress in a soil mass" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"sx= 2000 //lb/ft^3\n", +"sy= 2500 //lb/ft^3\n", +"T= 800 //lb/ft^3\n", +"t= 0.348//radians\n", +"//calculations\n", +"s1= (sx+sy)/2+sqrt(((sy-sx)/2)^2+T^2)\n", +"s2= (sx+sy)/2-sqrt(((sy-sx)/2)^2+T^2)\n", +"sn= (sx+sy)/2+(sy-sx)*cos(2*t)/2-T*sin(2*t)\n", +"Tn= (sy-sx)*sin(2*t)/2+T*cos(2*t)\n", +"//results\n", +"printf ('principle stress s1 = % 2f lb/ft^3 ',s1)\n", +"printf ('principle stress s2 = % 2f lb/ft^3 ',s2)\n", +"printf ('normal stress = % 2f lb/ft^3 ',sn)\n", +"printf ('shear stress = % 2f lb/ft^3 ',Tn)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"x= 3 //m\n", +"y= 4 //m\n", +"P= 5 //kN\n", +"z= 2 //m\n", +"//calculations\n", +"r= sqrt(x^2+y^2)\n", +"k= r/z\n", +"I= 3/(2*%pi*((r/z)^2+1)^2.5)\n", +"s= P*I/z^2\n", +"//results\n", +"printf ('verticle stress increase at 2m = % 4f kN/m^3 ',s)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/11-Compressibility_of_Soil.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/11-Compressibility_of_Soil.ipynb new file mode 100644 index 0000000..3bab158 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/11-Compressibility_of_Soil.ipynb @@ -0,0 +1,439 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 11: Compressibility of Soil" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.10: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"To=3000 // lb/ft^2\n", +"eo=1.1\n", +"e1=0.9\n", +"e=eo-e1\n", +"ea=(eo+e1)/2\n", +"T1=3000 // lb/ft^2\n", +"T=1 // in\n", +"t = 2 // min\n", +"m=(e/T1)/(1+ea)\n", +"U=50\n", +"Tv=0.197\n", +"Gw=62.4 //lb/ft^3\n", +"Cv=Tv*(T/(2*12)^2)/t\n", +"k=Cv*m*Gw *10^7\n", +"printf('a)k = %f x10^-7 ft/min\n',k)\n", +"\n", +"\n", +"U=60\n", +"Tv=0.286\n", +"H=6\n", +"t60=Tv*H^2/(Cv*60*24)\n", +"printf(' b)t60 =%f days',t60)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.11: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"t50=19\n", +"Hd=2.24/2\n", +"Cv=0.197*Hd^2/t50\n", +"printf('Cv = %f cm^2/min',Cv)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.12: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"LL=40\n", +"Cc=0.009*(LL-10)\n", +"H=10*12\n", +"eo=1.0\n", +"Gss=120\n", +"Gsc=110\n", +"Gd=100\n", +"To=10*Gd +10*(Gss-62.4)+10*(Gsc-62.4)/2\n", +"\n", +"Tt=0.408\n", +"Tm=0.232\n", +"Tb=0.019\n", +"Tav= (Tt+4*Tm+Tb)/6\n", +"Sc=Cc*H*log10((To+Tav*1000)/To)/(1+eo)\n", +"printf('Sc = %f in',Sc)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.13: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"H = 6\n", +"Cc = 0.28\n", +"eo = 0.9\n", +"Cv = 0.36\n", +"To=210\n", +"Tp=115\n", +"Sc= Cc*H*log10((To+Tp)/To)/(1+eo)\n", +"t2=9\n", +"Hd=3\n", +"Tv=Cv*t2/Hd^2\n", +"U=0.67\n", +"Tf=0.677*Tp\n", +"printf('Tf = %f kN/m^2',Tf)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"Tz=150\n", +"b=1\n", +"l=2\n", +"z=5*b\n", +"Es= (10000*2 + 8000*1 +12000*2)/5\n", +"a=4\n", +"H=z\n", +"m=l/b\n", +"n=2*H/b\n", +"F1=0.641 //from tables 11.1 and 11.2\n", +"F2=0.031\n", +"u=0.3\n", +"Is= F1 + ((2-u)/(1-u))*F2\n", +"If=0.71 //from table 11.3\n", +"Sef= Tz *a*b/l *(1-u^2)*Is*If/Es\n", +"Ser=0.93*Sef\n", +"printf('The elastic settlement at the centre of foundation = %f m',Ser)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"// one value of e is done \n", +"Gs=2.75\n", +"A=30.68\n", +"Ms=128\n", +"p=1\n", +"Hs=Ms/(A*Gs*p)\n", +"H=2.540\n", +"Hv=H-Hs\n", +"e=Hv/Hs\n", +"printf('the value of e for give values = %f',e)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"e11=0.9\n", +"e21=0.8\n", +"T2=4\n", +"T1=2\n", +"Cc= (e11-e21)/log10(T2/T1) // from loading branch\n", +"e1=0.67\n", +"e2=0.655\n", +"Cs=(e1-e2)/log10(T2/T1)\n", +"k=Cs/Cc\n", +"T3=12\n", +"e3=e11-Cc*log10(T3/T1)\n", +"printf('Compression index Cc= %f\n',Cc)\n", +"printf(' Cs/Cc = %f\n',k)\n", +"printf(' e3 = %f',e3)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.4: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"Gd=14\n", +"Gss=18\n", +"Gsc=19\n", +"Gw=9.81\n", +"To= 2*Gd+4*(Gss-Gw)+2*(Gsc-Gw)\n", +"LL=40\n", +"Cc=0.009*(LL-10)\n", +"H=4\n", +"T=100\n", +"e=0.8\n", +"Sc= Cc*H*log10((To+T)/To)/(1+e)\n", +"printf('a)Primary Consolidation Sc = %f m\n',Sc)\n", +"\n", +"Tc=190\n", +"Cs=Cc/6\n", +"Sc= Cs*H*log10((To+T)/To)/(1+e)\n", +"printf(' b)Primary Consolidation Sc = %f m\n',Sc)\n", +"\n", +"Tc=170\n", +"Sc= Cc*H*log10((To+T)/Tc)/(1+e)+ Cs*H*log10(Tc/To)/(1+e)\n", +"printf(' c)Primary Consolidation Sc = %f m\n',Sc)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.5: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"Gs=18\n", +"Gw=9.81\n", +"H=10\n", +"eo=1.1\n", +"To=5*(Gs-Gw)\n", +"T1=48\n", +"T=To+T1\n", +"e1=1.045 // void ratio corresponding to T \n", +"e=eo-e1\n", +"Sc=H*e/(1+eo)\n", +"printf('The settlement in the field Sc = %f m',Sc)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.6: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"T=8.5\n", +"eo=0.8\n", +"Cc=0.28\n", +"To=2650\n", +"T1=970\n", +"C1=0.02\n", +"t2=5\n", +"t1=1.5\n", +"H=8.5*12\n", +"epr=Cc*log10((To+T1)/To)\n", +"ep=eo-epr\n", +"C2=C1/(1+ep)\n", +"Sc=epr*H/(1+eo)\n", +"Ss=C2*H*log10(t2/t1)\n", +"TS=Sc+Ss\n", +"printf('Total consolidation settlement of the clay = %f in',TS)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.7: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//T50 = Cvtlab /H^2 lab = Cvtfield?H^2 fiels\n", +"tl=140\n", +"Hf=3\n", +"Hd=0.025/2\n", +"tf=tl*Hf^2/Hd^2\n", +"k=tf/(3600*24)\n", +"printf('t field = %f days',k)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.8: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//Tv is directly proportional to U^2\n", +"t1=93.333\n", +"U2=30\n", +"U1=50\n", +"t2=t1*U2^2/U1^2\n", +"printf('t2 = %f days',t2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11.9: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"t90=75*24*60*60 // time in sec\n", +"T90=0.848\n", +"Hd=1.5*100 //in cm\n", +"Cv=T90*Hd^2/t90\n", +"printf('Cv = %f cm^2/sec',Cv)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/12-Shear_Strength_of_Soil.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/12-Shear_Strength_of_Soil.ipynb new file mode 100644 index 0000000..3921e14 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/12-Shear_Strength_of_Soil.ipynb @@ -0,0 +1,252 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 12: Shear Strength of Soil" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"D=50 // in mm\n", +"A= %pi/4 *(D/1000)^2\n", +"// solving for test 1 \n", +"N=150\n", +"Sp=157.5\n", +"Sr=44.2\n", +"Tf=Sp/A\n", +"Tr=Sr/A\n", +"// from graph\n", +"k=tand(27)\n", +"k1=tand(14.6)\n", +"\n", +"printf('Peak strength Tf = 40+ t*%f\n',k )\n", +"printf(' Residual strength Tr = t*%f',k1)\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"T3=16 // lb/in^2\n", +"Tf=25 // lb/in^2\n", +"T1=T3+Tf\n", +"a= asind((T1-T3)/(T1+T3)) // Mohr's circle\n", +"printf('a)Angle of friction,a = %f\n',a)\n", +"b= 45+ a/2\n", +"printf(' b)Angle b that the failure plane makes with the major principal plane = %f',b)\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"T1=41\n", +"T3=16\n", +"a=58\n", +"T=(T1+T3)/2 + (T1-T3)*cosd(2*a)/2\n", +"tf=(T1-T3)*sind(2*a)/2\n", +"printf('a)the normal stress T = %f lb/in^2',T)\n", +"printf(' and the shear stress tf = %f lb/in^2',tf)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12.4: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//For normally consolidated clay, c' = 0.\n", +"a=30\n", +"T3=10\n", +"T1=T3*(tand(45+a/2))^2\n", +"Tf=T1-T3\n", +"printf('The deviator stress at failure = %f lb/in^2',Tf)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12.5: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"T13=70\n", +"T1f=130\n", +"T11=T13+T1f\n", +"\n", +"T23=160\n", +"T2f=223.5\n", +"T21=T23+T2f\n", +"\n", +"a= 2*(atand(((T11-T21)/(T13-T23))^0.5)-45)\n", +"c= (T11-T13*(tand(45+a/2))^2)/(2*tand(45+a/2))\n", +"printf('the shear strength parameter c = %f kN/m^2',c)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12.6: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"T3=12\n", +"Tf=9.1\n", +"T1=T3+Tf\n", +"u=6.8\n", +"a= asind((T1-T3)/(T1+T3))\n", +"\n", +"a1= asind((T1-T3)/(T1+T3-2*u))\n", +"\n", +"printf('a)Consolidated-undrained angle of shearing resistance = %f degrees\n',a)\n", +"printf(' b)Drained friction angle = %f degrees',a1)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12.7: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"T3=12\n", +"a=27.8\n", +"T1=T3*(tand(45+a/2))^2\n", +"Tf=T1-T3\n", +"printf('the deviator stress at failure = %f lb/in^2',Tf)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12.8: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"PI=28\n", +"OCR=3.2\n", +"To=160\n", +"Kn=0.11+0.0037*PI\n", +"Ko=OCR^0.8 * Kn\n", +"Cu=Ko*To\n", +"printf('the average undrained shear strength of the clay = %f kN/m^2',Cu)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/13-Lateral_Earth_Pressure.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/13-Lateral_Earth_Pressure.ipynb new file mode 100644 index 0000000..9da380b --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/13-Lateral_Earth_Pressure.ipynb @@ -0,0 +1,350 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 13: Lateral Earth Pressure" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.10: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"H=28\n", +"C=210\n", +"b=10\n", +"G=118\n", +"c=20\n", +"kh=0.1\n", +"Ka=tand(45-c/2)\n", +"zo=2*C/(G*(Ka))\n", +"n=zo/(H-zo)\n", +"Nac=1.60\n", +"Nav=0.375\n", +"L=1.17\n", +"Pae= G*(H-zo)^2*(L*Nav)-C*(H-zo)*Nac\n", +"printf('The magnitude of the active force, Pae = %f lb/ft',Pae)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"OCR=2\n", +"a=30\n", +"Ko=(1-sind(a))*(OCR)^sind(a)\n", +"//at z=0\n", +"To1=0\n", +"Th1=0\n", +"u1=0\n", +"//at z=10\n", +"To2=10*100\n", +"Th2=Ko*To2\n", +"u2=0\n", +"//at z=15\n", +"To3= 10*100+5*(122.4-62.4)\n", +"Th3=Ko*To3\n", +"u3=5*62.4\n", +"//Lateral force Po =Area 1 + Area 2+ Area3+ Area 4\n", +"Po =(1/2)*10*707+5*707+(1/2)*5*212.1+(1/2)*5*312\n", +"z=((3535)*(5+10/3)+3535*(5/2)+530.3*(5/3)+780*(5/3))/Po\n", +"printf('z = %f ft',z)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//c=0\n", +"a=36\n", +"G=16\n", +"Ka=(1-sind(a))/(1+sind(a))\n", +"//at z=0 Tp=0\n", +"z=6\n", +"To=G*z\n", +"Ta=Ka*To\n", +"Pa=z*Ta/2\n", +"\n", +"printf('a)Rankine active force per unit length of the wall = %f kN/m',Pa)\n", +"printf(' and the location of the resultant is z = 2m\n')\n", +"\n", +"\n", +"p=36\n", +"G=16\n", +"Kp=(1+sind(a))/(1-sind(a))\n", +"//at z=0 Tp=0\n", +"z=6\n", +"To=G*z\n", +"Tp=Kp*To\n", +"Pp=z*Tp/2\n", +"\n", +"printf(' b)Rankine passive force per unit length of the wall = %f kN/m',Pp)\n", +"printf(' and the location of the resultant is z = 2m')\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"H=12\n", +"a=20\n", +"b=20\n", +"G=115\n", +"c=30\n", +"Oa= asind(sind(a)/sind(c))-a+2*b\n", +"Ka= (cosd(a-b)*sqrt(1+(sind(c))^2-2*sind(c)*cosd(Oa)))/((cosd(b))^2*(cosd(a)+sqrt((sind(c))^2-(sind(a))^2)))\n", +"Pa=G*H^2*Ka/2\n", +"B= atand((sind(c)*sind(Oa))/(1-(sind(c)*cosd(Oa))))\n", +"printf('The active force Pa per unit length of the wall = %f lb/ft\n',Pa)\n", +"printf( ' The resultant will act a distance of 12/3 = 4 ft above the bottom of the wall with B = %f degree',B)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.4: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"a=30\n", +"Ka1=(1-sind(a))/(1+sind(a))\n", +"a=35\n", +"Ka2=(1-sind(a))/(1+sind(a))\n", +"//at z=0 so T0=0\n", +"//atz=3\n", +"To=3*16\n", +"Ta1=Ka1*To\n", +"Ta2=Ka2*To\n", +"\n", +"// At z=6\n", +"To=3*16+3*(18-9.81)\n", +"Ta2=Ka2*To\n", +"\n", +"Pa =(1/2)*3*16+3*13.0+ (1/2)*3*36.1\n", +"z= (24 *(3+3/3)+39.0*(3/2)+54.15*(3/3))/Pa\n", +"printf('The force per unit length of the wall = %f kN/m\n',Pa)\n", +"printf (' The location of the resultant = %f m ',z)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.5: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"Ka= (tand(1))^2\n", +"G=16.5\n", +"cu=10\n", +"H=6\n", +"//at z=0\n", +"z=0\n", +"Ta=G*z-2*cu\n", +"//zt z=6\n", +"z=6\n", +"Ta=G*z-2*cu\n", +"\n", +"zo=2*cu/G\n", +"// Before the tensile crack occurs\n", +"Pa= G*H^2/2 - 2*cu*H\n", +"printf('Pa before the tensile crack occurs = %f kN/m\n',Pa)\n", +"//After the tensile crack occurs\n", +"Pa=(H-zo)*Ta/2\n", +"printf(' Pa after the tensile crack occurs = %f kN/m',Pa)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.6: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"H=15\n", +"a=10\n", +"G=118\n", +"b=20\n", +"C=250\n", +"Zo=2*C*sqrt((1+sind(b))/(1-sind(b)))/G\n", +"//at z=0 Ta=0\n", +"//at z=15 \n", +"z=15\n", +"K=0.3\n", +"Ta=G*z*K*cosd(a)\n", +"Pa=(H -Zo)*Ta/2\n", +"printf('The Rankine active force Pa on the retaining wall after the tensile crack occurs = %f lb/ft',Pa)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.7: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"c=30\n", +"b=15\n", +"a=10\n", +"Ka=0.3872 // from table 13.8\n", +"H=4\n", +"G=15\n", +"Pa=G*H^2*Ka/2\n", +"printf('The active force per unit length Pa = %f kN/m\n',Pa)\n", +"printf(' The resultant will act at a vertical distance equal to H/3 = 4/3 = 1.33 m above \n the bottom of the wall and will be inclined at an angle of 15to the back face of the wall.')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.9: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"kh=0.2\n", +"kv=0\n", +"H=4\n", +"a=0\n", +"b=0\n", +"c=15\n", +"d=30\n", +"G=15.5\n", +"B= atand(kh/(1-kv))\n", +"b1=b+B\n", +"a1=a+B\n", +"Ka=0.452\n", +"Pa=G*H^2*Ka/2\n", +"Pae=Pa*(1-kv)*((cosd(b1))^2/((cosd(b))^2*(cosd(B))^2))\n", +"Ka=0.3014\n", +"Pa=G*H^2*Ka/2\n", +"P1=Pae-Pa\n", +"z= ((Pa*H/3)+P1*0.6*H)/Pae\n", +"printf('Pae = %f kN/m\n',Pae)\n", +"printf(' Z = %f m',z)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/14-Lateral_Earth_Pressure_Curved_Failure_Surface.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/14-Lateral_Earth_Pressure_Curved_Failure_Surface.ipynb new file mode 100644 index 0000000..73ef301 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/14-Lateral_Earth_Pressure_Curved_Failure_Surface.ipynb @@ -0,0 +1,104 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 14: Lateral Earth Pressure Curved Failure Surface" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 14.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"G=15.7\n", +"a=0\n", +"b=15\n", +"c=30\n", +"H=3\n", +"Kp=4.977 // from table 13.9\n", +"Pp=Kp*G*H^2/2\n", +"printf('a)the passive force = %f kN/m\n',Pp)\n", +"// for part b\n", +"Kp=4.53\n", +"Pp=Kp*G*H^2/2\n", +"printf(' b)the passive force = %f kN/m\n',Pp)\n", +"// for part c\n", +"Kp=4.13\n", +"Pp=Kp*G*H^2/2\n", +"printf(' c)the passive force = %f kN/m\n',Pp)\n", +"//for part d\n", +"Kp=4.56\n", +"Pp=Kp*G*H^2/2\n", +"printf(' d)the passive force = %f kN/m\n',Pp)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 14.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"G=16\n", +"H=7\n", +"c=30\n", +"Ta=0.65*G*H*(tand(45-c/2))^2\n", +"A=Ta*3*3/4\n", +"B1=Ta*3-54.61\n", +"C=Ta*4*4/4\n", +"B2=Ta*4-97.08\n", +"s=2\n", +"As=A*s\n", +"Bs=(B1+B2)*s\n", +"Cs=C*s\n", +"printf( 'The strut loads at level A = %f kN\n',As)\n", +"printf( ' The strut loads at level B = %f kN\n',Bs)\n", +"printf( ' The strut loads at level C = %f kN',Cs)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/15-Slope_Stability.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/15-Slope_Stability.ipynb new file mode 100644 index 0000000..6968ca3 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/15-Slope_Stability.ipynb @@ -0,0 +1,333 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 15: Slope Stability" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 15.11: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"C=20\n", +"G=18.5\n", +"H=21.62\n", +"c=25\n", +"r=0.25\n", +"Fs=3.1*tand(c)\n", +"printf('Fs = %f',Fs)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 15.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"Gs=17.8\n", +"Gw=9.81\n", +"C=10\n", +"c=20\n", +"b=15\n", +"H=6\n", +"G=Gs-Gw\n", +"Fs= C/(Gs*H*cosd(b)*cosd(b)*tand(b))+G*tand(c)/(Gs*tand(b))\n", +"printf('a)The factor of safety = %f \n',Fs)\n", +"Fs=2\n", +"H=2.247/(Fs-0.61)\n", +"printf(' b)H= %f m',H)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 15.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"G=105\n", +"c=15\n", +"C=600\n", +"b=45\n", +"Fs=3\n", +"Cd=C/Fs\n", +"c1= atand(tand(c)/Fs)\n", +"\n", +"H= 4*Cd*(sind(b)*cosd(c1)/(1-cosd(b-c1)))/G\n", +"printf('The depth of the cut slope = %f ft',H)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 15.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"Cu=40\n", +"G=17.5\n", +"b=60\n", +"a=35\n", +"c=72.5\n", +"m=0.195\n", +"Hc=Cu/(G*m)\n", +"r=Hc/(2*sind(a)*sind(c/2))\n", +"BC=Hc*((1/tand(a))-(1/tand(b)))\n", +"printf('a)The maximum depth Hc = %f m\n',Hc)\n", +"printf(' b)The radius, r = %f m\n',r)\n", +"printf(' c)The distance BC.= %f m',BC)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 15.4: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"Gs=17.29\n", +"d=9.15\n", +"d1=6.1\n", +"D=d/d1\n", +"a=40\n", +"m=0.175\n", +"b=40\n", +"H=6.1\n", +"Cu=H*Gs*m\n", +"printf('a)The undrained cohesion of the clay Cu = %f kN/m^2\n',Cu)\n", +"printf(' b)The nature of the critical circle is midpointcircle\n')\n", +"d=1.5\n", +"b=40\n", +"n=0.9\n", +"D1=n*H\n", +"printf(' c)Distance = %f m',D1)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 15.5: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"Fs=1\n", +"b=56\n", +"Kh=0.25\n", +"M=3.66\n", +"Cu=500\n", +"G=100\n", +"Hc=Cu*M/G\n", +"printf('a)The maximum depth = %f ft\n',Hc)\n", +"Fs=2\n", +"H=Cu*M/(G*Fs)\n", +"printf(' b)H= %f ft',H)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 15.6: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"b=45\n", +"c=20\n", +"C=24\n", +"G=18.9\n", +"m=0.06\n", +"Hc=C/(G*m)\n", +"Cd=G*Hc*m\n", +"Fc=C/Cd\n", +"printf('a)Critical height of slope = %f \n',Hc)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 15.7: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"FSs=1\n", +"c=20\n", +"G=18.9\n", +"C=24\n", +"Hcr=C/(G*tand(c)*0.17)\n", +"printf('a)Critical height Hc = %f m\n',Hcr)\n", +"H=10\n", +"k=C/(G*H*tand(c))\n", +"Fs=4*tand(c)\n", +"printf(' b)Fs = %f',Fs)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 15.8: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"W=22.4\n", +"C=20\n", +"a=70\n", +"s=sind(a)\n", +"c=cosd(a)\n", +"l=2.924\n", +"Wn=W*s\n", +"Wn1=W*c\n", +"//doing this to all values\n", +"F1=30.501\n", +"F2=776.75\n", +"F3=1638\n", +"Fs=(F1*C+F3*tand(C))/F2\n", +"printf('Fs = %f',Fs)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 15.9: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"C=20\n", +"G=18.5\n", +"r=0.25\n", +"H=21.62\n", +"C=25\n", +"b= atand(0.5)\n", +"//from table 15.3 \n", +"m=1.624\n", +"n=1.338\n", +"Fs=m-n*r\n", +"printf(' The value of Fs for D= 1 is %f',Fs)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/16-Soil_Bearing_Capacity_for_Shallow_Foundations.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/16-Soil_Bearing_Capacity_for_Shallow_Foundations.ipynb new file mode 100644 index 0000000..a56a48e --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/16-Soil_Bearing_Capacity_for_Shallow_Foundations.ipynb @@ -0,0 +1,263 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 16: Soil Bearing Capacity for Shallow Foundations" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 16.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"c=20\n", +"// from table 16.1\n", +"Nc=17.69\n", +"Nq=7.44\n", +"Ng=3.64\n", +"\n", +"Df=3\n", +"G=110\n", +"q=G*Df\n", +"\n", +"C=200\n", +"B=4\n", +"\n", +"Qu= C*Nc+q*Nq+G*B*Ng/2\n", +"\n", +"Fs=3\n", +"Qall=Qu/Fs\n", +"printf('Qa = %f lb/ft^2',Qall)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 16.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"G=18.15\n", +"qa=30000*9.81/1000\n", +"\n", +"Nc=57.75\n", +"Nq=41.44\n", +"Ng=45.41\n", +"C=0\n", +"q=G*1\n", +"B=1\n", +"(1.3*C*Nc+q*Nq+0.4*G*B*Ng)*B^2/3 == qa\n", +"B= sqrt(294.3/(250.7+109.9))\n", +"printf(' B = %f m',B)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 16.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"B=1.2\n", +"L=1.2\n", +"c=32\n", +"C=0\n", +"Df=1\n", +"G=16\n", +"Nq=23.18\n", +"Ng=22.02\n", +"Nc=1\n", +"Lqs=1+0.1*B*(tand(45+c/2))^2/L\n", +"Lgs=Lqs\n", +"Lqd=1+0.1*Df*tand(45+c/2)/B\n", +"Lgd=Lqd\n", +"Lcs=1\n", +"Lcd=1\n", +"Gs=19.5\n", +"q=0.5*G+0.5*(Gs-9.81)\n", +"Qu= C*Lcs*Lcd*Nc+q*Lqs*Lqd*Nq+(Gs-9.81)*Lgs*Lgd*B*Ng/2\n", +"Qa=Qu/3\n", +"Q=Qa*B^2\n", +"printf('the gross load = %f kN',Q)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 16.4: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"e=0.1\n", +"B=1\n", +"X=B-2*e\n", +"Y=1.5\n", +"B1=0.8\n", +"L1=1.5\n", +"c=30\n", +"Df=1\n", +"Nq=18.4\n", +"Ng=15.668\n", +"q=1*18\n", +"G=18\n", +"Lqs=1+e*(B1/L1)*(tand(45+c/2))^2\n", +"Lgs=Lqs\n", +"Lqd=1+e*(Df/B1)*tand(45+c/2)\n", +"Lgd=Lqd\n", +"qu=q*Lqs*Lqd*Nq+Lgs*Lgd*G*B1*Ng/2\n", +"Qu=qu*B1*L1\n", +"printf('The magnitude of the gross ultimate load = %f kN',Qu)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 16.5: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"B=1.5\n", +"Df=0.75\n", +"e=0.1*B\n", +"G=17.5\n", +"c=30\n", +"C=0\n", +"q=G*Df\n", +"Nq=18.4\n", +"Ng=15.668\n", +"Lqd=1+0.1*(Df/B)*tand(45+c/2)\n", +"Lgd=Lqd\n", +"Quc=q*Nq*Lqd+Lgd*B*Ng/2\n", +"k=0.8\n", +"a=1.754\n", +"Qua=Quc*(1-a*(e/B)^k)\n", +"printf('The gross ultimate load per unit length = %f kN',Qua)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 16.6: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"Qup=280\n", +"Bp=0.7 // in m\n", +"Bf=1.5\n", +"Quf=Qup*Bf/Bp\n", +"printf('The ultimate bearing capacity = %f kN/m^2',Quf)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 16.7: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"a=2500\n", +"//doing for the first values only\n", +"Bf=4\n", +"Bp=0.305\n", +"q=a/Bf^2\n", +"Sep=4\n", +"Sef=Sep*(2*Bf/(Bf+Bp))^2\n", +"printf('Sef = %f mm',Sef)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/18-Subsoil_Exploration.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/18-Subsoil_Exploration.ipynb new file mode 100644 index 0000000..1c6c804 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/18-Subsoil_Exploration.ipynb @@ -0,0 +1,105 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 18: Subsoil Exploration" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 18.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//solving for z=5 only\n", +"To=0.275\n", +"Cn=To^(-0.5)\n", +"N60=8\n", +"N160=Cn*N60\n", +"printf('(N1)60 = %f',N160)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 18.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"z=5\n", +"To=0.275\n", +"Cn=2/(1+To)\n", +"N60=8\n", +"N160=Cn*N60\n", +"printf('(N1)60 = %f',N160)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 18.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"pa=1 // 14.7 lb/in^2 = 1ton/ft^2\n", +"To=0.275 // ton/ ft^2\n", +"N60=8\n", +"c= atand((N60/(12.2+20.3*(To/pa)))^0.34)\n", +"printf('The average soil friction angle = %f',c)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/2-weight_volume_relationships.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/2-weight_volume_relationships.ipynb new file mode 100644 index 0000000..9da5bd1 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/2-weight_volume_relationships.ipynb @@ -0,0 +1,61 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 2: weight volume relationships" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"//from graph\n", +"d= 0.15 //mm\n", +"w= 0.17 //mm\n", +"a= 0.27 //mm\n", +"//calculations\n", +"C= a/d\n", +"c= w^2/(a*d)\n", +"//results\n", +"printf ('uniformity coefficient = % f ',C)\n", +"printf ('coefficient of gradation = % 2f ',c)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/3-weight_volume_relationships.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/3-weight_volume_relationships.ipynb new file mode 100644 index 0000000..7064ae0 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/3-weight_volume_relationships.ipynb @@ -0,0 +1,164 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 3: weight volume relationships" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"V= 1.2 //m^3\n", +"M= 2350 //Kg\n", +"w= 0.086\n", +"G= 2.71\n", +"W= 1000 //kg/m^3\n", +"//calculations\n", +"R= M/V\n", +"D= M/((1+w)*V)\n", +"e= (G*W/D)-1\n", +"n= e/(e+1)\n", +"S= (w*G/e)*100\n", +"v= (M-(M/(1+w)))/W\n", +"//results\n", +"printf ('moist density = % f kg/m^3 ',R)\n", +"printf ('dry density = % f kg/m^3 ',D)\n", +"printf ('void ratio = % 3f ',e)\n", +"printf ('porosity = % 3f ',n)\n", +"printf ('Degree of saturation = % 3f ',S)\n", +"printf ('volume of water in soil sample = % 3f m^3 ',v)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"n= 0.4\n", +"G= 2.68\n", +"w= 0.12\n", +"R= 1000 //kg/m^3\n", +"V= 10 //m^3\n", +"//calculations\n", +"d= G*R*(1-n)*(1+w)\n", +"s= ((1-n)*G+n)*R\n", +"M= s-d\n", +"m= M*V\n", +"//results\n", +"printf ('mass of water to be added for full saturation = % f kg ',m)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.4: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"d= 16.19 //kN/m^3\n", +"w= 0.23\n", +"W= 9.81 //kN/m^3\n", +"//calculations\n", +"R= d*(1+w)\n", +"G= d/(W-d*w)\n", +"e= w*G\n", +"//results\n", +"printf ('satuarated unit weight = % 2f kN/m^3 ',R)\n", +"printf ('specific gravity = % 2f ',G)\n", +"printf ('void ratio = % 2f ',e)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.5: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"G= 2.68\n", +"w= 0.12\n", +"d= 1794.4 //kg/m^3\n", +"W= 1000 //kg/m^3\n", +"emax= 0.75\n", +"emin= 0.4\n", +"//calculation\n", +"e= (G*W*(1+w)/d)-1\n", +"D= ((emax-e)/(emax-emin))*100\n", +"//results\n", +"printf ('relative density of compaction in percentage = % f ',D)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/4-plasticity_and_structure_of_soil.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/4-plasticity_and_structure_of_soil.ipynb new file mode 100644 index 0000000..6fa3c99 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/4-plasticity_and_structure_of_soil.ipynb @@ -0,0 +1,61 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 4: plasticity and structure of soil" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"V1= 24.6 //cm^3\n", +"V2= 15.9 //cm^3\n", +"M1= 44 //g\n", +"M2= 30.1 //g\n", +"W= 1 //g/cm^3\n", +"//calculations\n", +"SL= (((M1-M2)/M2)*100)-(((V1-V2)/M2)*W*100)\n", +"//results\n", +"printf ('shrinkage limit of the soil = % 2f ',SL)\n", +"" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/6-soil_compaction.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/6-soil_compaction.ipynb new file mode 100644 index 0000000..ecda060 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/6-soil_compaction.ipynb @@ -0,0 +1,127 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 6: soil compaction" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"G= 2.6\n", +"LL= 20\n", +"P= 20\n", +"//calclations\n", +"R= (4804574*G-195.55*(LL)^2+156971*(P)^0.5-9527830)^0.5\n", +"n= (1.195e-4)*((LL)^2)-1.964*G-(6.617e-5)*(P)+7.651\n", +"w= %e^n\n", +"//results\n", +"printf ('maximum dry density = % f kg/m^3 ',R)\n", +"printf ('optimum moisture content = % 2f ',w)\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"do= 1570 //kg/m^3\n", +"mo= 0.545 //kg\n", +"M1= 7.59 //kg\n", +"M2= 4.78 //kg\n", +"M3= 3.007 //kg\n", +"w= 0.102 //\n", +"dmax= 19 //KN/m^3\n", +"//calculations\n", +"Ms= M1-M2\n", +"Mc= Ms-mo\n", +"Vh= Mc/do\n", +"Dc= M3/Vh\n", +"Du= Dc*9.81/1000\n", +"f= Du/(1+w)\n", +"Rc= f*100/dmax\n", +"//results\n", +"printf ('dry unit weight of compaction in the field = % 2f kN/m^3 ',f)\n", +"printf ('relative compaction in the field = % f ',Rc)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.4: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"D1= 0.36 //mm\n", +"D2= 0.52 //mm\n", +"D5= 1.42 //mm\n", +"//calculations\n", +"Sn= 1.7*(sqrt((3/(D5)^2)+(1/(D2)^2)+(1/(D1)^2)))\n", +"//results\n", +"printf ('sustainabilty number = % f ',Sn)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/7-permiability.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/7-permiability.ipynb new file mode 100644 index 0000000..b89d2d3 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/7-permiability.ipynb @@ -0,0 +1,355 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 7: permiability" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.10: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"k1= 0.302e-7 //cm/sec\n", +"k2= 0.12e-7 //cm/sec\n", +"e1= 1.1\n", +"e2= 0.9\n", +"e= 0.75\n", +"//calcualtions\n", +"n= (log10((k1/k2)*((1+e1)/(1+e2))))/log10(e1/e2)\n", +"C= k1/(e1^n/(1+e1))\n", +"k= C*(e^n/(1+e))\n", +"//results\n", +"printf ('hydraulic conductivity = %e cm/sec ',k')" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.11: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"H1= 2 //m\n", +"H2= 3 //m\n", +"H3= 4 //m\n", +"k1= 1e-4 //cm/sec\n", +"k2= 3.2e-2 //cm/sec\n", +"k3= 4.1e-5 //cm/sec\n", +"//calculations\n", +"H= H1+H2+H3\n", +"Kh= (1/H)*((k1*H1)+(k2*H2)+(k3*H3))\n", +"Kv= H/((H1/k1)+(H2/k2)+(H3/k3))\n", +"P= Kh/Kv\n", +"//results\n", +"printf ('ration of equivalent hydraulic conductivity = % 2f ',P)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.12: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"H= 450 //mm\n", +"h= 150 //mm\n", +"k1= 1e-2 //cm/sec\n", +"k2= 3e-3 //cm/sec\n", +"k3= 4.9e-4 //cm/sec\n", +"h1= 300 //mm\n", +"//calculations\n", +"Kv= H/(h*(1/k1+1/k2+1/k3))\n", +"i= h1/H\n", +"q= Kv*i*100*3600\n", +"//results\n", +"printf ('rate of water supply = % 2f cm/hr ',q)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"L= 30 //cm\n", +"A= 177 //cm^2\n", +"h= 50 //cm\n", +"Q= 350 //cm^3\n", +"t= 300 //sec\n", +"//claculations\n", +"k=Q*L/(A*h*t)\n", +"//results\n", +"printf ('hydraulic conductivity = % 3f cm/sec ',k)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"L= 203 //mm\n", +"A= 10.3 //cm^2\n", +"a= 0.39 //cm^2\n", +"h0= 508 //mm\n", +"h180= 305 //mm\n", +"t= 180 //sec\n", +"//calculations\n", +"k= 2.303*a*L*log10(h0/h180)/(A*t)\n", +"//results\n", +"printf ('hydraulic conductivity of sand = % 2f in/sec ',k)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of varilables\n", +"k= 3e-7 //cm/sec\n", +"n= 0.0911e-4 //g*sec/cm^2\n", +"dw= 1 //g/cc\n", +"//calculations\n", +"K= k*n/dw\n", +"//results\n", +"printf ('absolute premeability = % 4f cm^2 ',K)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.4: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"k= 5.3e-5 //m/sec\n", +"H= 3 //m\n", +"a= 0.139 //radians\n", +"//calculations\n", +"A= H*cos(a)\n", +"i= sin(a)\n", +"q= k*i*A*3600\n", +"//results\n", +"printf ('rate of seepage = % 4f m^3/hr/m ',q)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.5: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"L= 50 //m\n", +"k= 0.08e-2//m/sec\n", +"h= 4 //m\n", +"H1= 3 //m\n", +"H= 8 //m\n", +"a= 0.139 //radians\n", +"//calculations\n", +"i= h*cos(a)/L\n", +"A= H1*cos(a)\n", +"q= k*i*A\n", +"//results\n", +"printf ('flow rate = % 2f m^3/sec/m ',q)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.6: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"k1= 0.02 //cm/sec\n", +"e1= 0.5 \n", +"e2= 0.65\n", +"//calculations\n", +"k2= k1*(e2^3/(1+e2))/(e1^3/(1+e1))\n", +"//results\n", +"printf ('hydraulic conductivity at void ratio of 0.65 = % 2f cm/sec ',k2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.8: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"e= 0.6\n", +"D10= 0.09 //mm\n", +"//calculations\n", +"k= 2.4622*(D10^2*(e^3/(1+e)))^0.7825\n", +"//results\n", +"printf ('hydraulic conductivity = % 4f cm/sec ',k)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.9: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"e= 0.6\n", +"D10= 0.09 //mm\n", +"D60= 0.16 //mm\n", +"//calculations\n", +"Cu=D60/D10\n", +"k= 35*(e^3/(1+e))*(Cu^0.6)*(D10^2.32)\n", +"//results\n", +"printf ('hydraulic conductivity = % 3f cm/sec ',k)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/8-seepage.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/8-seepage.ipynb new file mode 100644 index 0000000..99bf4d9 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/8-seepage.ipynb @@ -0,0 +1,200 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 8: seepage" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 8.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"H1= 12 //in\n", +"H2= 20 //in\n", +"z= 8 //in\n", +"h1= 24 //in\n", +"h= 20 //in\n", +"k1= 0.026 //in/sec\n", +"D= 3 //in\n", +"//calculations\n", +"k2= H2*k1/((z/(1-h/h1))-H1)\n", +"i= h1/(H1+H2)\n", +"A= %pi/4*D^2\n", +"keq= (H1+H2)/((H1/k1)+(H2/k2))\n", +"q= keq*A*i*3600\n", +"//results\n", +"printf ('rate of water flow = % 2f in^3/hr ',q)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 8.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"Nd= 6\n", +"H1= 5.6 //m\n", +"H2= 2.2 //m\n", +"k= 5e-5 //cm/sec\n", +"dL= 4.1 //m\n", +"//calculations\n", +"H= (H1-H2)/Nd\n", +"h1= 5.61-H\n", +"h2= 5.61-5*H\n", +"q= 2.38*(H1-H2)*k/Nd\n", +"i= H/dL\n", +"//results\n", +"printf ('at point a,water will rise to height of = % 3f m ',h1)\n", +"printf ('at point b,water will rise to height of = % 3f m ',h2)\n", +"printf ('total rate of seepage per unit lenghth = %e m^3/sec/m ',q)\n", +"printf ('average hydraulic gradient at c = % 3f ',i)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 8.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"k1= 5.67 //ft/day\n", +"k2= 11.34 //ft/day\n", +"//from graph\n", +"Nd= 8\n", +"Nf= 2.5\n", +"H= 20\n", +"//calculations\n", +"q= sqrt(k1*k2)*H*Nf/Nd\n", +"//results\n", +"printf ('average rate of flow = % 2f ft^3/day/ft ',q)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 8.4: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"B= 6 //m\n", +"L= 120 //m\n", +"s= 3 //m\n", +"T= 6 //m\n", +"x= 2.4 //m\n", +"H= 5 //m\n", +"k= 0.008 //cm/sec\n", +"//calculations\n", +"b=B/2\n", +"a1= b/T\n", +"a2= s/T\n", +"a3= x/b\n", +"Q= 0.378*k*H*L*36*24\n", +"//results\n", +"printf ('seepage under the dam = % 2f m^3/day ',Q)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 8.5: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"b= %pi/4 //degrees\n", +"a= %pi/6.//degrees\n", +"B= 10 //ft\n", +"H= 20 //ft\n", +"h= 25 //ft\n", +"k= 2e-4 //ft/min\n", +"//calculations\n", +"r= H/tan(b)\n", +"d= 0.3*r+(h-H)/tan(b)+B+h/tan(a)\n", +"L= d/cos(a)-sqrt((d/cos(a))^2-(H/sin(a))^2)\n", +"q= k*L*tan(a)*sin(a)*24*60\n", +"//results\n", +"printf ('seepage rate = % 4f ft^3/day/ft ',q)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/9-in_situ_stresses.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/9-in_situ_stresses.ipynb new file mode 100644 index 0000000..5423a71 --- /dev/null +++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/9-in_situ_stresses.ipynb @@ -0,0 +1,176 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 9: in situ stresses" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.1: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"Ds= 16.5 //kN/m^3\n", +"S= 19.25 //kN/m^3\n", +"g= 9.8 //kN/m^3\n", +"h1= 6 //m\n", +"h2= 13 //m\n", +"//at point A\n", +"Sa= 0\n", +"Ua= 0\n", +"Sa1= 0\n", +"//at point B\n", +"Sb= h1*Ds\n", +"Ub= 0 \n", +"Sb1= Sb-Ub\n", +"//at point C\n", +"Sc= h1*Ds+h2*S\n", +"Uc= h2*g\n", +"Sc1= Sc-Uc\n", +"//results\n", +"printf ('total pressure at C= % 2f kN/m^3 ',Sc)\n", +"printf ('pore water pressure at C = % 2f kN/m^3 ',Uc)\n", +"printf ('effective stress at point C= % 2f kN/m^3 ',Sc1)\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.2: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"h= 20 //ft\n", +"g= 120 //kg/ft^3\n", +"h1= 12 //ft\n", +"w= 62.4 //kg/ft^3\n", +"//calculations\n", +"H= h-(h1*w/g)\n", +"//results\n", +"printf ('maximu depth that can be made in clay = % 2f ft ',H)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.3: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"G= 2.68\n", +"e= 0.52\n", +"g= 9.81 //kN/m^3\n", +"h1= 0.7 //m\n", +"h2= 1 //m\n", +"h3= 1.5 //m\n", +"h4= 2 //m\n", +"//calculations\n", +"//for soil A\n", +"sa= (G+e)*g/(1+e)\n", +"//point a\n", +"Sa= h1*g+h2*sa\n", +"u= (h2+h1+h3/2)*g\n", +"Es= Sa-u\n", +"//point b\n", +"sb= h1*g+h4*sa\n", +"ub= (h4+h1+h3)*g\n", +"Eb= sb-ub\n", +"i= h3/2\n", +"s= i*g\n", +"//results\n", +"printf ('effective stress at point a= % 2f kN/m^2 ',Es)\n", +"printf ('effective stress at point b= % 2f kN/m^2 ',Eb)\n", +"printf ('seepage force per unit voume = % 2f kN/m^3 ',s)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9.4: solved.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"//initialisation of variables\n", +"C0= 0.357\n", +"H1= 30.5 //ft\n", +"H2= 5 //ft\n", +"w= 62.4 // lb/ft^3\n", +"D= 20\n", +"g= 112 // lb/ft^3\n", +"//calculations\n", +"G= g-w\n", +"FS= D*G/(C0*w*(H1-H2))\n", +"//results\n", +"printf ('safety factor = % 2f ',FS)" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} |