diff options
| author | Thomas Stephen Lee | 2015-09-04 22:04:10 +0530 |
|---|---|---|
| committer | Thomas Stephen Lee | 2015-09-04 22:04:10 +0530 |
| commit | 41f1f72e9502f5c3de6ca16b303803dfcf1df594 (patch) | |
| tree | f4bf726a3e3ce5d7d9ee3781cbacfe3116115a2c /Principles_Of_Foundation_Engineering/Chapter11.ipynb | |
| parent | 9c9779ba21b9bedde88e1e8216f9e3b4f8650b0e (diff) | |
| download | Python-Textbook-Companions-41f1f72e9502f5c3de6ca16b303803dfcf1df594.tar.gz Python-Textbook-Companions-41f1f72e9502f5c3de6ca16b303803dfcf1df594.tar.bz2 Python-Textbook-Companions-41f1f72e9502f5c3de6ca16b303803dfcf1df594.zip | |
add/remove/update books
Diffstat (limited to 'Principles_Of_Foundation_Engineering/Chapter11.ipynb')
| -rwxr-xr-x | Principles_Of_Foundation_Engineering/Chapter11.ipynb | 650 |
1 files changed, 0 insertions, 650 deletions
diff --git a/Principles_Of_Foundation_Engineering/Chapter11.ipynb b/Principles_Of_Foundation_Engineering/Chapter11.ipynb deleted file mode 100755 index 0afdb4a0..00000000 --- a/Principles_Of_Foundation_Engineering/Chapter11.ipynb +++ /dev/null @@ -1,650 +0,0 @@ -{
- "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 |