summaryrefslogtreecommitdiff
path: root/Principles_Of_Geotechnical_Engineering_by_B_M_Das/7-permiability.ipynb
diff options
context:
space:
mode:
Diffstat (limited to 'Principles_Of_Geotechnical_Engineering_by_B_M_Das/7-permiability.ipynb')
-rw-r--r--Principles_Of_Geotechnical_Engineering_by_B_M_Das/7-permiability.ipynb355
1 files changed, 355 insertions, 0 deletions
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
+}
generated by cgit (git 2.34.1) at 2025-07-07 01:31:57 +0000