From fba055ce5aa0955e22bac2413c33493b10ae6532 Mon Sep 17 00:00:00 2001 From: hardythe1 Date: Tue, 5 May 2015 14:21:39 +0530 Subject: add books --- .../Chapter12_2.ipynb | 443 +++++++++++++++++++++ 1 file changed, 443 insertions(+) create mode 100755 Principles_Of_Foundation_Engineering/Chapter12_2.ipynb (limited to 'Principles_Of_Foundation_Engineering/Chapter12_2.ipynb') diff --git a/Principles_Of_Foundation_Engineering/Chapter12_2.ipynb b/Principles_Of_Foundation_Engineering/Chapter12_2.ipynb new file mode 100755 index 00000000..fcc14387 --- /dev/null +++ b/Principles_Of_Foundation_Engineering/Chapter12_2.ipynb @@ -0,0 +1,443 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:1439e2e46e40a0d47bd1c7e06a03c417c0e1f55ad415a0af82536f86f7ecfea9" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter12:Drilled-Shaft Foundations" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex12.1:Pg-609" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 12.1\n", + "\n", + "import math\n", + "Ap=math.pi/4*1.75**2; # area in m^2\n", + "FS=4; # factor of safety\n", + "Nq=37.75; \n", + "L=8;\n", + "Es=50000.0;\n", + "mus=0.265;\n", + "pu=100.0;\n", + "Db=1.75; # in m\n", + "q=6*16.2+2*19.2;\n", + "phi=36*math.pi/180.0;\n", + "Fqs=1+math.tan(phi);\n", + "Fqd=1+2*math.tan(phi)*(1-math.sin(phi))**2*math.atan(L/Db);\n", + "Ir=Es/(2*(1+mus)*q*math.tan(phi));\n", + "delta=0.005*(1-phi/20*180/math.pi+25/20.0)*q/pu;\n", + "Irr=Ir/(1+Ir*delta);\n", + "Fqc=math.exp(-3.8*math.tan(phi)+(3.07*math.sin(phi)*math.log10(2*Irr))/(1+math.sin(phi)));\n", + "Qp=Ap*(q*(Nq-1)*Fqs*Fqd*Fqc);\n", + "Qpall=Qp/FS;\n", + "print round(Qpall,2),\"allowed load in kN\"\n", + "print \"due to rounding off error there is slight change in answer\"\n", + "\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "6648.7 allowed load in kN\n", + "due to rounding off error there is slight change in answer\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex12.2:Pg-610" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 12.2\n", + "import math\n", + "Ap=math.pi/4*1.75**2; # area in m^2\n", + "q=135.6; \n", + "w=0.83;\n", + "FS=4; # factor of safety\n", + "phi=36; # given angle\n", + "Nq=0.21*math.exp(0.17*phi);\n", + "Qp=Ap*q*(w*Nq-1); # in kN\n", + "Qpall=Qp/FS; # in kN\n", + "print round(Qpall),\"allowed load in kN\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "6383.0 allowed load in kN\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex12.3:Pg-611" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 12.3\n", + "\n", + "import math\n", + "Ap=math.pi/4*1.5**2; # area in m^2\n", + "Db=1.5; # in m\n", + "z=3.0;\n", + "p=math.pi*1;\n", + "Li=6.0;\n", + "N60=30.0;\n", + "sigmazi=16*z;\n", + "Beta=2.0-0.15*z**0.75;\n", + "fi=Beta*sigmazi; # in kN/m^2\n", + "qp=57.5*N60; # in kN/m^2\n", + "qpr=1.27/Db*qp; # in kN/m^2\n", + "Qunet=qpr*Ap+fi*p*Li; # in kN\n", + "print round(Qunet,2),\"allowed load in kN\"\n", + "#part b\n", + "k1=0.315; #from table\n", + "k2=12.0/1.5/1000*100.0;\n", + "Qunet2=qpr*Ap*k1+fi*p*Li*k2; # in kN\n", + "print round(Qunet2,2),\"allowed load in kN\"\n", + "\n", + "# the answer is slightly different in textbook due to approximation" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "4081.11 allowed load in kN\n", + "2013.14 allowed load in kN\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex12.4:Pg-617" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 12.4\n", + "\n", + "Nc=9;\n", + "Ap=math.pi/4*1.5**2; # area in m^2\n", + "cu=105; # in kN/m^2\n", + "Qpnet=Ap*cu*Nc; # in kN\n", + "print round(Qpnet),\"net ultimate bearing point capacity in kN\"\n", + "#part2\n", + "alpha=0.4;\n", + "Ds=1.5; # in m \n", + "p=math.pi*Ds;\n", + "Qs=alpha*p*(50*8+105*3); # in kN\n", + "print int(Qs),\"skin resistance in kN\"\n", + "#part3\n", + "FS=3; # factor of safety\n", + "Qu=Qpnet/FS+Qs/FS; # in kN\n", + "print round(Qu,2),\"working load in kN\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "1670.0 net ultimate bearing point capacity in kN\n", + "1347 skin resistance in kN\n", + "1005.9 working load in kN\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex12.5:Pg-618" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 12.5\n", + "\n", + "import math\n", + "cub=3000;\n", + "L=20+5; # in ft\n", + "Db=4; # in ft\n", + "Ap=math.pi/4*Db**2; # area in ft^2\n", + "alpha=0.55;\n", + "cu1=800; # in lb/ft^2\n", + "L1=7; # in ft\n", + "L2=5.5; # in ft\n", + "cu2=1200; # in lb/ft^2\n", + "p=math.pi*2.5;\n", + "k=alpha*p*(cu1*L1+cu2*L2);#f*p*deltaLi\n", + "j1=6*cub*(1+0.2*L/Db);\n", + "j2=9*cub;\n", + "qp=min(j1,j2);\n", + "Qu=k/1000+qp*Ap/1000; # in kip\n", + "print round(Qu),\"allowed load in kip\"\n", + "#part b\n", + "k1=0.57; #from table\n", + "k2=0.89;\n", + "Qunet2=qp*Ap*k1+k*k2; # in kip\n", + "print round(Qunet2/1000,2),\"allowed load in kip\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "392.0 allowed load in kip\n", + "240.3 allowed load in kip\n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex12.6:Pg-621" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 12.6\n", + "\n", + "import math\n", + "Qws=1005-250; # in kN\n", + "Qwp=250; # in kN\n", + "epsilon=0.65; \n", + "L=11; # in m\n", + "Ds=1.5; # in m\n", + "Es=14000; # in kN/m^2\n", + "Ap=math.pi/4*1.5**2; # area in m^2\n", + "Ep=21e6; # in kN/m^2\n", + "Cp=0.04; # in kN/m^2\n", + "Db=1.5;\n", + "mus=0.3;\n", + "p=math.pi*1.5;\n", + "Nc=9;\n", + "qp=105*Nc; # in kN/m^2\n", + "se1=(Qwp+epsilon*Qws)*L/(Ap*Ep); # in m \n", + "se2=Qwp*Cp/(Db*qp); # in m\n", + "Iws=2+0.35*math.sqrt(L/Ds);\n", + "se3=Qws/p/L*Ds/Es*(1-mus**2)*Iws; # in m\n", + "se=se1+se2+se3; # in m\n", + "print round(se*1000,2),\" is net settlement in mm\"\n", + "\n", + "# the answer is slightly different in textbook due to approximation" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "11.46 is net settlement in mm\n" + ] + } + ], + "prompt_number": 17 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex12.7:Pg-628" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 12.7\n", + "\n", + "import math\n", + "import numpy\n", + "from scipy.optimize import fsolve\n", + "Ds=1.0;\n", + "Ep=22e6;\n", + "Ri=1.0;\n", + "cu=100.0;\n", + "Ip=math.pi*Ds**4/64.0;\n", + "Qc=7.34*Ds**2*Ep*Ri*(cu/Ep/Ri)**0.6;\n", + "print round(Qc,2),\"bearing force in kN\"\n", + "Mc=3.86*Ds**3*Ep*Ri*(cu/Ep/Ri)**0.6;\n", + "print round(Mc,2),\"bearing moment in kNm\"\n", + "#from figure\n", + "xoQM=0.0046*1;\n", + "xoMQ=0.0041*1;\n", + "xo=0.5*(xoQM+xoMQ);\n", + "print round(xo*1000,2),\"net ground line deflection in mm\"\n", + "#partb\n", + "Ip=0.049;\n", + "Qg=150.0;\n", + "Mg=200.0;\n", + "def f(T):\n", + " return 338e-6*T**3+300.6e-6*T**2-0.00435\n", + "[x]=fsolve(f,2);\n", + "T=x;\n", + "k=[0, 0.4, 0.6, 0.8, 1.0, 1.1, 1.25];#z/T\n", + "Am=[0, 0.36, 0.52, 0.63, 0.75, 0.765, 0.75];\n", + "Bm=[1.0, 0.98, 0.95, 0.9, 0.845, 0.8, 0.73];\n", + "print \"z/T\\t Am\\t Bm\\t Mz(kN-m)\\n\"\n", + "Mz=numpy.zeros(7)\n", + "for i in range(0,7):\n", + " Mz[i]=Am[i]*Qg*T+Bm[i]*Mg;\n", + " print k[i],\"\\t\",round(Am[i],2),\"\\t\",round(Bm[i],2),\"\\t\",round(Mz[i],2)\n", + "\n", + "print round(1*T,2),\"depth in m\"\n", + "#partc\n", + "Mmax=400;\n", + "sigma=Mmax*Ds/2/Ip;\n", + "print round(sigma,2),\"tensile stress in kN/m**2\"\n", + "#partd\n", + "#from figure\n", + "k=8.5;\n", + "L=k*1;\n", + "print L,\"length in m\"\n", + "\n", + "# the answer is slightly different in textbook due to approximation\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "100615.56 bearing force in kN\n", + "52912.27 bearing moment in kNm\n", + "4.35 net ground line deflection in mm\n", + "z/T\t Am\t Bm\t Mz(kN-m)\n", + "\n", + "0 \t0.0 \t1.0 \t200.0\n", + "0.4 \t0.36 \t0.98 \t308.4\n", + "0.6 \t0.52 \t0.95 \t352.35\n", + "0.8 \t0.63 \t0.9 \t376.69\n", + "1.0 \t0.75 \t0.84 \t403.16\n", + "1.1 \t0.77 \t0.8 \t398.84\n", + "1.25 \t0.75 \t0.73 \t380.16\n", + "2.08 depth in m\n", + "4081.63 tensile stress in kN/m**2\n", + "8.5 length in m\n" + ] + } + ], + "prompt_number": 31 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Ex12.8:Pg-634" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#example 12.8\n", + "\n", + "qu=3000; # lb/in^2\n", + "Ds=3*12; #in inch\n", + "L=15*12; # in inch\n", + "FS=3; # factor of safety\n", + "Ecore=0.36e6; # in lb/in^2\n", + "f=min(2.5*qu**0.5,0.15*qu);\n", + "Qu=math.pi*Ds*L*f/1000; # in kip\n", + "Emass=Ecore*(0.266*80-1.66); # in lb/in^2\n", + "Ec=17.9*Emass; # in lb/in^2\n", + "Ac=math.pi/4*Ds**2; # area in in^2\n", + "If=0.35;\n", + "se=Qu*L/Ac/Ec+Qu*If/Ds/Emass;\n", + "Qall=Qu/FS; # in kip\n", + "print round(Qall),\"allowed load in kip\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "929.0 allowed load in kip\n" + ] + } + ], + "prompt_number": 35 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file -- cgit