{
"metadata": {
"name": "",
"signature": "sha256:2a29c7b9c2c5d3914e37e2b72bd62a8b948e23c0720d69e0b40a353ece4dd036"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "heading",
"level": 1,
"metadata": {},
"source": [
"Chapter11-Compressibility of Soil"
]
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex1-pg303"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"#evaluvate The elastic settlement at the centre of foundation\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",
"print'%s %.3f %s'%('The elastic settlement at the centre of foundation =',Ser,'m')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"The elastic settlement at the centre of foundation = 0.012 m\n"
]
}
],
"prompt_number": 1
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex2-pg312"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"\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",
"print'%s %.3f %s'%('the value of e for give values =',e,'')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"the value of e for give values = 0.674 \n"
]
}
],
"prompt_number": 2
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex3-pg321"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math\n",
"#calculate Compression index\n",
"e11=0.9\n",
"e21=0.8\n",
"T2=4.\n",
"T1=2.\n",
"Cc= (e11-e21)/math.log10(T2/T1) ## from loading branch\n",
"e1=0.67\n",
"e2=0.655\n",
"Cs=(e1-e2)/math.log10(T2/T1)\n",
"k=Cs/Cc\n",
"T3=12.\n",
"e3=e11-Cc*math.log10(T3/T1)\n",
"print'%s %.2f %s'%('Compression index Cc= ',Cc,'')\n",
"print'%s %.2f %s'%(' Cs/Cc = ',k,'')\n",
"print'%s %.2f %s'%(' e3 = ',e3,'')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Compression index Cc= 0.33 \n",
" Cs/Cc = 0.15 \n",
" e3 = 0.64 \n"
]
}
],
"prompt_number": 3
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex4-pg323"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math\n",
"#calculate Primary Consolidation Sc in three parts\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*math.log10((To+T)/To)/(1.+e)\n",
"print'%s %.2f %s'%('a)Primary Consolidation Sc = ',Sc,' m')\n",
"\n",
"\n",
"Tc=190\n",
"Cs=Cc/6\n",
"Sc= Cs*H*math.log10((To+T)/To)/(1+e)\n",
"print'%s %.2f %s'%(' b)Primary Consolidation Sc =',Sc,'m')\n",
"\n",
"\n",
"Tc=170\n",
"Sc= Cc*H*math.log10((To+T)/Tc)/(1+e)+ Cs*H*math.log10(Tc/To)/(1+e)\n",
"print'%s %.3f %s'%(' c)Primary Consolidation Sc =',Sc,' m')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"a)Primary Consolidation Sc = 0.21 m\n",
" b)Primary Consolidation Sc = 0.04 m\n",
" c)Primary Consolidation Sc = 0.047 m\n"
]
}
],
"prompt_number": 10
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex5-pg325"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"#calculate The settlement in the field Sc\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",
"print'%s %.2f %s'%('The settlement in the field Sc = ',Sc,' m')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"The settlement in the field Sc = 0.26 m\n"
]
}
],
"prompt_number": 11
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex6-pg329"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math\n",
"#calculate Total consolidation settlement of the clay\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*math.log10((To+T1)/To)\n",
"ep=eo-epr\n",
"C2=C1/(1.+ep)\n",
"Sc=epr*H/(1.+eo)\n",
"Ss=C2*H*math.log10(t2/t1)\n",
"TS=Sc+Ss\n",
"print'%s %.1f %s'%('Total consolidation settlement of the clay =',TS,' in')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Total consolidation settlement of the clay = 2.8 in\n"
]
}
],
"prompt_number": 12
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex7-pg336"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math\n",
"#calculate t field\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",
"print'%s %.1f %s'%('t field = ',k,' days')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"t field = 93.3 days\n"
]
}
],
"prompt_number": 13
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex8-pg336"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math\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",
"print'%s %.2f %s'%('t2 =',t2,' days')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"t2 = 33.60 days\n"
]
}
],
"prompt_number": 14
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex9-pg337"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math\n",
"#evaluvate Cv\n",
"#intilization variable\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",
"print'%s %.3f %s'%('Cv =',Cv,' cm^2/sec')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Cv = 0.003 cm^2/sec\n"
]
}
],
"prompt_number": 15
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex10-pg337"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math\n",
"#calculate K and t60\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",
"print'%s %.3f %s'%('a)k = ',k,' x10^-7 ft/min')\n",
"\n",
"\n",
"U=60\n",
"Tv=0.286\n",
"H=6\n",
"t60=Tv*H**2/(Cv*60*24)\n",
"print'%s %.1f %s'%(' b)t60 =',t60,' days')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"a)k = 3.557 x10^-7 ft/min\n",
" b)t60 = 41.8 days\n"
]
}
],
"prompt_number": 5
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex11-pg344"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math\n",
"#Cv\n",
"t50=19\n",
"Hd=2.24/2\n",
"Cv=0.197*Hd**2/t50\n",
"print'%s %.3f %s'%('Cv = ',Cv,' cm^2/min')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Cv = 0.013 cm^2/min\n"
]
}
],
"prompt_number": 6
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex12-pg346"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math\n",
"#calculate SC\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*math.log10((To+Tav*1000.)/To)/(1.+eo)\n",
"print'%s %.3f %s'%('Sc =',Sc,' in')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Sc = 0.826 in\n"
]
}
],
"prompt_number": 8
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Ex13-pg356"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"#intilization variable\n",
"#Calculate total primary\n",
"import math\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*math.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",
"print'%s %.1f %s'%('Tf =',Tf,' kN/m^2')\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Tf = 77.9 kN/m^2\n"
]
}
],
"prompt_number": 7
}
],
"metadata": {}
}
]
}