From 92cca121f959c6616e3da431c1e2d23c4fa5e886 Mon Sep 17 00:00:00 2001 From: hardythe1 Date: Tue, 7 Apr 2015 15:58:05 +0530 Subject: added books --- .../Chapter11.ipynb | 650 +++++++++++++++++++++ 1 file changed, 650 insertions(+) create mode 100755 Principles_Of_Foundation_Engineering/Chapter11.ipynb (limited to 'Principles_Of_Foundation_Engineering/Chapter11.ipynb') diff --git a/Principles_Of_Foundation_Engineering/Chapter11.ipynb b/Principles_Of_Foundation_Engineering/Chapter11.ipynb new file mode 100755 index 00000000..0afdb4a0 --- /dev/null +++ b/Principles_Of_Foundation_Engineering/Chapter11.ipynb @@ -0,0 +1,650 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:a4054b3d6b79658bcd763c43dc94ca6de4cd49698c12f912e4c83a9a61f47154" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter11:Pile Foundations" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.1:Pg-532" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.1\n", + "import math\n", + "\n", + "#parta\n", + "phi=30.0; # angle in degree\n", + "pa=2000.0; \n", + "q=100*50/1000.0;\n", + "Nq=55.0;\n", + "Ap=16*16/16/12; # area in ft^2\n", + "Qp=Ap*q*Nq; # in kip\n", + "qp=0.4*pa*Nq*math.tan(phi*math.pi/180)*Ap; # in lb\n", + "print round(Qp,2),\"ultimate load in lb\"\n", + "print round(qp/1000,2),\"ultimate load in kip\"\n", + "print \"there is change in answer because of calculation mistake in the book\"\n", + "\n", + "#partb\n", + "Nsigma=36;\n", + "Ap=16*16.0/12.0/12;\n", + "q=110*50.0/1000;\n", + "Qp=Ap*q*Nsigma*((1+2.0*(1-math.sin(phi*math.pi/180)))/3); # in kip\n", + "print round(Qp,2),\"ultimate load in kip\"\n", + "#partc\n", + "Nq=18.4;\n", + "Qp=Ap*q*Nq; # in kip\n", + "print round(Qp,2),\"ultimate load in kip\"\n", + "\n", + "# ANSWER IN THE BOOK IS WRONG" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "275.0 ultimate load in lb\n", + "25.4 ultimate load in kip\n", + "there is change in answer because of calculation mistake in the book\n", + "234.67 ultimate load in kip\n", + "179.91 ultimate load in kip\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.2:Pg-533" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.2\n", + "\n", + "import math\n", + "#parta\n", + "K=1.3;\n", + "f0=0;\n", + "Delta=0.8*30; # in ft\n", + "D=16.0/12; # in ft\n", + "L1=50.0;\n", + "p=4*16/12.0; # in ft\n", + "Gamma=110/1000.0; # in lb/ft^3\n", + "L=15*D; # in ft\n", + "sigma=Gamma*L; # in kip/ft^2\n", + "f20=K*sigma*math.tan(Delta*math.pi/180); # kip/ft^2\n", + "Qs=(f0+f20)/2*(p*L)+f20*p*(L1-L);\n", + "print round(Qs,2),\"ultimate load in kip\"\n", + "#partb\n", + "FS=4; # factor of safety\n", + "Qp=56.45/3+234.7/3+179.9/3; # in kip\n", + "Qu=Qs+Qp; # in kip\n", + "Qall=Qu/FS; # in kip\n", + "print round(Qall,2),\"is allowed load in kip\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "271.65 ultimate load in kip\n", + "107.17 is allowed load in kip\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.3:Pg-534" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.3\n", + "import math\n", + "K=0.25;\n", + "Ap=16*16.0/12/12; # area in ft^2\n", + "phi=30*math.pi/180;\n", + "Nq=25;\n", + "q=110*50.0/1000; # in kip\n", + "sigmao=q/2; # in kip/ft^2\n", + "p=4*16.0/12; # in ft\n", + "L=50; # in ft\n", + "FS=4; # factor of safety\n", + "Qu=q*Nq*Ap+K*sigmao*math.tan(0.8*phi)*p*L; # in kip\n", + "Qall=Qu/FS; # in kip\n", + "print round(Qall,1),\"allowed load in kip\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "81.5 allowed load in kip\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.4:Pg-535" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.4\n", + "\n", + "import math\n", + "import numpy\n", + "FS=4; # factor of safety\n", + "Ap=0.1295; # area in m^2\n", + "Nc=9;\n", + "cu2=100;\n", + "Qp=Ap*Nc*cu2; # in kN\n", + "D=[5, 10, 30]; # depth in m\n", + "avgD=[2.5, 7.5,20.0]; # average depth in m\n", + "sigma=[45.0, 110.5, 228.5]; # in KN/m^2\n", + "cu=[30, 30, 100]; # in kN/m^2\n", + "alpha=[0.6, 0.9, 0.725];\n", + "L=[5, 5, 20]; # in m\n", + "p=math.pi*0.406;\n", + "Qs=0; # in kN\n", + "cusig=numpy.zeros(3)\n", + "print round(Qp,2),\"bearing capacity in kN\"\n", + "print \"depth (m)\\t avg Depth(m)\\t avgVerticalStress(kN/m**2)\\t cu(kN/m**2)\\t cu/sigma\\t alpha\\n\"\n", + "for i in range(0,3):\n", + " cusig[i]=cu[i]/sigma[i];\n", + " Qs=Qs+alpha[i]*cu[i]*L[i]*p;\n", + " print round(D[i],2),\"\\t \\t \\t\",round(avgD[i],2),\"\\t \\t\",round(sigma[i],2),\"\\t\\t\\t \",round(cu[i],2),\"\\t \",round(cusig[i],2),\"\\t\\t \",round(alpha[i],2),\"\\n\"\n", + "print round(Qs,2),\"bearing capacity in kN\"\n", + "#part2\n", + "Lambda=0.136;\n", + "L=30;\n", + "fav=Lambda*(178.48+2*76.7);\n", + "Qs2=p*L*fav; # in kN\n", + "#part3\n", + "fav1=13;\n", + "fav2=31.9;\n", + "fav3=93.43;\n", + "Qs3=p*(fav1*5+fav2*5+fav3*20); # in kN\n", + "print round(Qs3,1),\"bearing capacity in kN\"\n", + "Qsavg=Qs/3+Qs2/3+Qs3/3; # in kN\n", + "Qu=Qp+Qsavg # in kN\n", + "Qall=Qu/FS; # in kN\n", + "print round(Qall,1),\"allowed bearing capacity in kN\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "116.55 bearing capacity in kN\n", + "depth (m)\t avg Depth(m)\t avgVerticalStress(kN/m**2)\t cu(kN/m**2)\t cu/sigma\t alpha\n", + "\n", + "5.0 \t \t \t2.5 \t \t45.0 \t\t\t 30.0 \t 0.67 \t\t 0.6 \n", + "\n", + "10.0 \t \t \t7.5 \t \t110.5 \t\t\t 30.0 \t 0.27 \t\t 0.9 \n", + "\n", + "30.0 \t \t \t20.0 \t \t228.5 \t\t\t 100.0 \t 0.44 \t\t 0.72 \n", + "\n", + "2136.44 bearing capacity in kN\n", + "2669.7 bearing capacity in kN\n", + "573.6 allowed bearing capacity in kN\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.5:Pg-538" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.5\n", + "\n", + "import numpy\n", + "D=[6, 12, 20]; # depth in m\n", + "fc=[34.34, 54.94, 70.63]; # in kN/m**2\n", + "alpha=[0.84, 0.71, 0.63];\n", + "dL=[6, 6, 8]; # in m\n", + "p=4*0.305;\n", + "Qs=0;\n", + "Q=numpy.zeros(3)\n", + "print \" depth(m)\\t fc(kN/m**2)\\t alpha \\t \\t deltaL(m)\\t Q(kN)\\n\"\n", + "for i in range (0,3):\n", + " Q[i]=alpha[i]*fc[i]*p*dL[i];\n", + " Qs=Q[i]+Qs;\n", + " print D[i],\"\\t\\t \",fc[i],\"\\t \",alpha[i],\"\\t \",dL[i],\"\\t\\t \",round(Q[i],2)\n", + "\n", + "print round(Qs),\"bearing force in kN\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + " depth(m)\t fc(kN/m**2)\t alpha \t \t deltaL(m)\t Q(kN)\n", + "\n", + "6 \t\t 34.34 \t 0.84 \t 6 \t\t 211.15\n", + "12 \t\t 54.94 \t 0.71 \t 6 \t\t 285.53\n", + "20 \t\t 70.63 \t 0.63 \t 8 \t\t 434.29\n", + "931.0 bearing force in kN\n" + ] + } + ], + "prompt_number": 17 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex5.6:Pg-545" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.6\n", + "\n", + "import math\n", + "L=21; # in m\n", + "Qwp=502-350; # in kN\n", + "Qws=350; # in kN\n", + "Ap=0.1045; # area in m^2\n", + "Ep=21e6; # in kN/m^2\n", + "epsilon=0.62;\n", + "Se1=(Qwp+epsilon*Qws)*L/Ap/Ep; # in m\n", + "#part2\n", + "Iwp=0.85;\n", + "qwp=152/Ap;\n", + "Es=25e3; # in kN/m^2\n", + "D=0.356; # in m\n", + "mus=0.35;\n", + "Se2=qwp*D/Es*Iwp*(1-mus**2); # in m\n", + "#part3\n", + "p=1.168;\n", + "Iws=2+0.35*math.sqrt(L/D);\n", + "Se3=Qws/p/L*D/Es*Iws*(1-mus**2); # in m\n", + "Se=Se1+Se2+Se3; # in m\n", + "print round(Se*1000,1),\"settlement in mm\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "19.8 settlement in mm\n" + ] + } + ], + "prompt_number": 21 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.7:Pg-560" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.7\n", + "\n", + "Ep=207e6; # in kN/m^2\n", + "Ip=123e-6; # in m^4\n", + "nh=12000; # in kN/m^3\n", + "#from table 11.13\n", + "xz=0.008;\n", + "Ax=2.435;\n", + "T=(Ep*Ip/nh)**0.2;\n", + "Qg1=xz*Ep*Ip/Ax/T**3;\n", + "#part2\n", + "Fy=248000;\n", + "d1=0.254;\n", + "Am=0.772;\n", + "Mzmax=Fy*Ip*2/d1; # in Kn-m\n", + "Qg2=Mzmax/Am/T; # in kN\n", + "if Qg2>Qg1 :\n", + " Qg=Qg1;\n", + " print round(Qg,2),\"lateral load in kN\"\n", + "# there is slight variation in answer in textbook due to approximation" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "53.27 lateral load in kN\n" + ] + } + ], + "prompt_number": 25 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.8:Pg-561" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.8\n", + "\n", + "import math\n", + "#part1\n", + "Ep=207e6; # in kN/m^2\n", + "Ip=123e-6; # in m^4\n", + "nh=12000; # in kN/m^3\n", + "#from table 11.1a\n", + "xo=0.008; # in m\n", + "L=25;\n", + "Fy=248000; # yield stress in kN/m^2\n", + "D=0.254;\n", + "Am=0.772;\n", + "Gamma=18.0; # in kN/m^3\n", + "phi=35; # in angle\n", + "Kp=(math.tan(math.pi/4+phi*math.pi/360))**2;\n", + "My=Fy*Ip*2/D; # in kN-m\n", + "Qug=140*Kp*D**3*Gamma; # in kN\n", + "\n", + "#part2\n", + "Qg1=xo*(Ep*Ip)**0.6*nh**0.4/0.15/L; # in kN\n", + "\n", + "if Qug>Qg1:\n", + " Qg=Qg1;\n", + " print round(Qg,2),\"lateral load in kN\"\n", + "\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "40.2 lateral load in kN\n" + ] + } + ], + "prompt_number": 28 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.9:Pg-567" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.9\n", + "\n", + "import math\n", + "Wrh=30*12; # in kip-ft\n", + "E=0.8;\n", + "Wr=7.5; # in kip\n", + "S=1/8.0; \n", + "C=0.1;\n", + "FS=6; # in factor of safety\n", + "n=0.4; # Coefficient of restitution\n", + "Wp=12/12.0*12/12.0*80*150+550; # in lb\n", + "Wp=Wp/1000.0;\n", + "Qu=E*Wrh/(S+C)*(Wr+n**2.0*Wp)/(Wr+Wp); # in kip\n", + "Qall=Qu/FS; # in kip\n", + "print round(Qall),\"allowed bearing capacity in kip\"\n", + "#part2\n", + "He=30*12.0;\n", + "L=80*12.0;\n", + "Ap=12*12.0; # area in in^2\n", + "Ep=3e6/1000.0; # in kip/in^2\n", + "FS=4; # factor of safety\n", + "Qu=E*He/(S+math.sqrt(E*He*L/2.0/Ap/Ep)); # in kip\n", + "Qall2=Qu/FS; # in kip\n", + "print round(Qall2),\"allowed bearing capacity in kip\"\n", + "\n", + "#partc\n", + "a=27;\n", + "b=1;\n", + "He=30;\n", + "FS=3; # factor of safety\n", + "Qu=a*math.sqrt(E*He)*(b-math.log10(S)); # in kip\n", + "Qall3=Qu/FS; # in kip\n", + "print round(Qall3),\"allowed bearing capacity in kip\"\n", + "\n", + "\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "101.0 allowed bearing capacity in kip\n", + "104.0 allowed bearing capacity in kip\n", + "84.0 allowed bearing capacity in kip\n" + ] + } + ], + "prompt_number": 35 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.10:Pg-570" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.10\n", + "\n", + "Hp=350; # in HP\n", + "vp=0.0016; # in m/s\n", + "Sl=0.762e-3; # in m/cycle\n", + "f=115; # in Hz\n", + "Qu=(0.746*Hp+98*vp)/(vp+Sl*f); # in kN\n", + "print round(Qu),\"pile load capacity in kN\"\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "2928.0 pile load capacity in kN\n" + ] + } + ], + "prompt_number": 38 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.11:Pg-578" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.11\n", + "\n", + "Lg=9.92; # in ft\n", + "Bg=7.0; # in ft\n", + "n1=3.0;\n", + "Nc=8.75;\n", + "n2=4.0/1000;\n", + "Ap=14.0**2.0/12.0**2;\n", + "cup=1775.0;\n", + "a1=0.4;#alpha1\n", + "p=4*14.0/12.0;\n", + "cu1=1050.0; # in lb/ft^2\n", + "L1=15.0;\n", + "a2=0.54;#alpha2\n", + "cu2=1775.0; # in lb/ft^2\n", + "L2=45.0;\n", + "FS=4; # factor of safety\n", + "Qu=n1*n2*(9*Ap*cup+a1*p*cu1*L1+a2*p*cu2*L2); # in kip\n", + "Qu2=Lg*Bg*cup*Nc+2*(Lg+Bg)*(cu1*L1+cu2*L2); # in kip\n", + "print round(Qu2/1000),\"load in kip\"\n", + "Qall=Qu/FS; # in kip\n", + "print round(Qall),\"allowed load in kip\"\n", + "\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "4314.0 load in kip\n", + "757.0 allowed load in kip\n" + ] + } + ], + "prompt_number": 42 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex11.12:Pg-583" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 11.12\n", + "\n", + "import math\n", + "z1=21/2.0; # in ft\n", + "Lg=9.0; # in ft\n", + "Bg=6.0;# in ft\n", + "Qg=500*1000.0; # in kip\n", + "Cc1=0.3;\n", + "Cc2=0.2;\n", + "Cc3=0.25;\n", + "H2=12;\n", + "H3=6;\n", + "H1=21;\n", + "e1=0.82;\n", + "e2=0.7;\n", + "e3=0.75;\n", + "s1=Qg/(Lg+z1)/(Bg+z1); #sigma1 in lb/ft^3\n", + "s2=500*1000/(9+27)/(6+27);#sigma2 in lb/ft^3\n", + "s3=500*1000/(9+36)/(6+36);#sigma3 in lb/ft^3\n", + "ss1=6*105+(27+21/2)*(115-62.4);#sigmadash1 in lb/ft^3\n", + "ss2=6*105+(27+21)*(115-62.4)+(120-62.4)*6;#sigmadash2 in lb/ft^3\n", + "ss3=6*105+48*(115-62.4)+12*(120-62.4)+3*(122-62.4);#sigmadash3 in lb/ft^3\n", + "sc1=Cc1*H1/(1+e1)*math.log10((ss1+s1)/ss1); # in inch\n", + "sc2=Cc2*H2/(1+e2)*math.log10((ss2+s2)/ss2); # in inch\n", + "sc3=Cc3*H3/(1+e3)*math.log10((ss3+s3)/ss3); # in inch\n", + "sc=sc1+sc2+sc3; # in inch\n", + "print round(sc*12,1),\"total settlement in inch\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "9.6 total settlement in inch\n" + ] + } + ], + "prompt_number": 45 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file -- cgit