diff options
| author | kinitrupti | 2017-05-12 18:53:46 +0530 |
|---|---|---|
| committer | kinitrupti | 2017-05-12 18:53:46 +0530 |
| commit | 6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d (patch) | |
| tree | 22789c9dbe468dae6697dcd12d8e97de4bcf94a2 /Hydraulics_Made_Easy_by_R._S._Dighe | |
| parent | d36fc3b8f88cc3108ffff6151e376b619b9abb01 (diff) | |
| download | Python-Textbook-Companions-6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d.tar.gz Python-Textbook-Companions-6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d.tar.bz2 Python-Textbook-Companions-6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d.zip | |
Removed duplicates
Diffstat (limited to 'Hydraulics_Made_Easy_by_R._S._Dighe')
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/README.txt | 10 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/ch1.ipynb | 1015 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/ch10.ipynb | 497 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/ch2.ipynb | 368 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/ch3.ipynb | 527 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/ch4.ipynb | 994 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/ch5.ipynb | 668 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/ch6.ipynb | 1557 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/ch7.ipynb | 755 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/ch8.ipynb | 554 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/ch9.ipynb | 241 | ||||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/10MiscellaneousProblems.png | bin | 0 -> 135984 bytes | |||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/5FlowOfWaterOverWeirs.png | bin | 0 -> 117531 bytes | |||
| -rwxr-xr-x | Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/9ViscousFlow.png | bin | 0 -> 129763 bytes |
14 files changed, 7186 insertions, 0 deletions
diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/README.txt b/Hydraulics_Made_Easy_by_R._S._Dighe/README.txt new file mode 100755 index 00000000..5cbe8938 --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/README.txt @@ -0,0 +1,10 @@ +Contributed By: Alpesh Makwana +Course: bca +College/Institute/Organization: Zealous web technologies +Department/Designation: Developer +Book Title: Hydraulics Made Easy +Author: R. S. Dighe +Publisher: Royal Book Stall +Year of publication: 2005 +Isbn: 0415250706 +Edition: 1
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/ch1.ipynb b/Hydraulics_Made_Easy_by_R._S._Dighe/ch1.ipynb new file mode 100755 index 00000000..530a253b --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/ch1.ipynb @@ -0,0 +1,1015 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:d3eb95f5e18a78e74ef41bfca92729f5344f3fef9eebaa02f35bb683ba9d0489" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 1 : Hydrostatics" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.1 Page No : 6" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "Ar = 50 \t\t#area of ram in**2\n", + "Ap = 1./8 \t\t#area of plunger in**2\n", + "Wp = 5. \t\t#force lbs\n", + "\t\t\n", + "#CALCULATIONS\n", + "Pp = Wp/Ap\n", + "F = Pp*Ar\n", + "\t\t\n", + "#RESULTS\n", + "print 'weight supported by ram = %.f lbs'%(F)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "weight supported by ram = 2000 lbs\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.2 Page No : 6" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Dp = 1. \t\t#diameter of punger - in\n", + "Dr = 10. \t\t#diameter of ram - in\n", + "R = 12. #leverage of handle\n", + "W = 15. \t\t#wieght of body - tons\n", + "\t\t\n", + "#CALCULATIONS\n", + "Ar = math.pi*Dr**2/4\n", + "Ap = math.pi*Dp**2/4\n", + "P = W*2240/((Ar/Ap)*R)\n", + "\t\t\n", + "#RESULTS\n", + "print 'power applied to lever = %.f lbs'%(P)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "power applied to lever = 28 lbs\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.3 Page No : 7" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Dj = 1. \t\t#diameter of plunger - in\n", + "Dr = 2. \t\t#in\n", + "W = 40. \t\t#lbs\n", + "W1 = 1. \t\t#ton\n", + "rl = 20.\n", + "\t\t\n", + "#CALCULATIONS\n", + "Ap = math.pi*Dj**2/4\n", + "Ar = math.pi*Dr**2/4\n", + "Vrj = rl*Ar/Ap\n", + "e = W1*2240*100/(W*Vrj)\n", + "\t\t\n", + "#RESULTS\n", + "print 'efficiency of machine at this load = %.f percent'%(e)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "efficiency of machine at this load = 70 percent\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.4 Page No : 7" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "Dj = 1. \t\t#in\n", + "Dr = 2. \t\t#in\n", + "ns = 3. \t\t#strokes\n", + "h = 2. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "Ap = math.pi*Dj**2/4\n", + "Ar = math.pi*Dr**2/4\n", + "Vrj = Ar/Ap\n", + "ns1 = h*12*Vrj/ns\n", + "\t\t\n", + "#RESULTS\n", + "print ' working strokes = %.f strokes'%(ns1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + " working strokes = 32 strokes\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.5 Page No : 9" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "T = 40. \t\t#F\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "h = 50 \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "p = w*h/(12**2)\n", + "\t\t\n", + "#RESULTS\n", + "print ' pressure at a depth of 50 ft = %.2f lbs per in'%(p)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + " pressure at a depth of 50 ft = 21.67 lbs per in\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.6 Page No : 13" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "W = 64. \t\t#lbs/ft**3\n", + "h1 = 27. \t\t#ft\n", + "h2 = 9. \t\t#ft\n", + "w = 40. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "Pr = w*W*h1*h1/2\n", + "Pl = w*W*h2*h2/2\n", + "y1 = h1/3\n", + "y2 = h2/3\n", + "y = (Pr*y1-Pl*y2)/(Pr-Pl)\n", + "\t\t\n", + "#RESULTS\n", + "print ' point of application = %.2f ft'%(y)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + " point of application = 9.75 ft\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.7 Page No : 14" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 5. \t\t#ft\n", + "x = 3. \t \t#ft\n", + "w = 62.4 \t\t#lb/ft**3\n", + "a = 90. \t\t#degrees\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = math.pi/4*d**2\n", + "b = w*A*x\n", + "Ig = round(math.pi*d**4/64,2)\n", + "Io = Ig + A*x**2 * 1\n", + "h = Io/(A*x)\n", + "\n", + "#RESULTS\n", + "print 'depth of the pressure = %.2f ft'%(h)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "depth of the pressure = 3.52 ft\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.8 Page No : 15" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "w = 3. \t\t#ft\n", + "h = 4. \t \t#ft\n", + "ht = 30 \t\t#ft\n", + "W = 62.4 \t\t#ft**3\n", + "\t\t\n", + "#CALCULATIONS\n", + "Ap = w*h\n", + "X = ht+(h/2)\n", + "P = Ap*X*W\n", + "I0 = (w*h**3/12)+Ap*X**2\n", + "H = I0/(Ap*X)\n", + "\t\t\n", + "#RESULTS\n", + "print ' total pressure on the gate = %.2f ft'%(H)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + " total pressure on the gate = 32.04 ft\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.9 Page No : 15" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "w = 3. \t\t#ft\n", + "h = 4. \t\t#ft\n", + "ht = 30. \t\t#ft\n", + "W = 62.4 \t\t#ft**3\n", + "x = 2.22 \t\t#in\n", + "x1 = 4.5 \t\t#in\n", + "\t\t\n", + "#CALCULATIONS\n", + "Ap = w*h\n", + "X = ht+(h/2)\n", + "P = Ap*X*W\n", + "T = P*x/x1\n", + "T1 = P-T\n", + "\t\t\n", + "#RESULTS\n", + "print ' tension devoloped in the top bolt = %.f lbs'%(T)\n", + "print ' tension devoloped in the bottom bolt = %d lbs'%(T1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + " tension devoloped in the top bolt = 11821 lbs\n", + " tension devoloped in the bottom bolt = 12140 lbs\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.10 Page No : 16" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "w = 3. \t\t#ft\n", + "h = 15. \t\t#ft\n", + "d = 140. \t\t#lbs/ft**3\n", + "x = 6. \t \t#in\n", + "W = 62.4 \t\t#lbs/ft**3\n", + "\t\t\n", + "#CALCULATIONS\n", + "W1 = h*w*d\n", + "h = (W1*x*6/(W*12))**(1./3)\n", + "\t\t\n", + "#RESULTS\n", + "print ' height of water rise = %.2f ft'%(h)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + " height of water rise = 6.72 ft\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.11 Page No : 17" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "h = 5. \t\t#ft\n", + "d = 6. \t\t#ft\n", + "a = 30. \t\t#degrees\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = math.pi*d**2/4\n", + "X = h+(d/2)*math.sin(math.radians(a))\n", + "P = w*A*X\n", + "Ic = math.pi*d**4/64\n", + "I0 = Ic+A*X**2/(math.sin(math.radians(a)))**2\n", + "h = I0*(math.sin(math.radians(a)))**2/(A*X)\n", + "\t\t\n", + "#CALCULATIONS\n", + "print 'depth of the centre os pressure = %.2f ft '%(h)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "depth of the centre os pressure = 6.59 ft \n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.12 Page No : 18" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "w = 4. \t \t#ft\n", + "l = 4. \t\t #ft\n", + "X = 10. \t\t#ft\n", + "a = 45. \t\t#degrees\n", + "W = 100. \t\t#lbs\n", + "a1 = 60. \t\t#degrees\n", + "w1 = 62.4 \t\t#lbs/ft**3\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = w*l\n", + "X1 = round(X+(w/2)*math.sin(math.radians(a)),2)\n", + "Ig = round(w*l**3/12,2)\n", + "I0 = Ig+(A*X1**2/(math.sin(math.radians(a)))**2)\n", + "h = I0*(math.sin(math.radians(a)))**2/(A*X1)\n", + "P = round(w1*A*X1,-2)\n", + "h1 = round(h-X,2)\n", + "h2 = round(h1/math.sin(math.radians(a)),2)\n", + "T = (W*(l/2)*math.sin(math.radians(a))+P*h2)/(w*math.sin(math.radians(a1)))\n", + "\n", + "\n", + "\n", + "#RESULTS\n", + "print 'Pull in the chain = %d lbs '%(T)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Pull in the chain = 6885 lbs \n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.13 Page No : 20" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "w = 4. \t\t#ft\n", + "l = 4. \t\t#ft\n", + "X = 10. \t\t#ft\n", + "a = 45. \t\t#degrees\n", + "W = 62.4 \t\t#lbs/ft**3\n", + "u = 0.25\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = w*l\n", + "X1 = X+(w/2)*math.sin(math.radians(a))\n", + "P = W*A*X1\n", + "T = u*P\n", + "\t\t\n", + "#RESULTS\n", + "print 'magnitude of the lifting force = %.f lbs '%(round(T,-1)) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "magnitude of the lifting force = 2850 lbs \n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.14 Page No : 21" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "sg = 1.6\n", + "h = 10. \t\t#ft\n", + "h1 = 4. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "D = w*sg\n", + "W = w*(h+h1)**2/2\n", + "P = w*h\n", + "P1 = D*h1\n", + "P2 = (P*h/2)+P*h1+(h1*P1/2)\n", + "y = ((P*h*(h1+(h/3))/2)+P*h1*(h1/2)+P1*h1**2/6)/P2\n", + "\t\t\n", + "#RESULTS\n", + "print 'Position where P acts = %.1f ft above the base'%(y) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Position where P acts = 4.5 ft above the base\n" + ] + } + ], + "prompt_number": 14 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.15 Page No : 22" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "pa = 10. \t#lbs/in**2\n", + "h = 8. \t\t #ft\n", + "h1 = 6. \t\t#ft\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "pg = 10. \t\t#lbs/in**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "Pa = pa*144\n", + "Pa1 = w*h1\n", + "Pt = (Pa*h+Pa1*(h1/2))\n", + "y = (Pa*h*(h/2)+(Pa1*h1*(h-h1)/2))/Pt\n", + "\t\t\n", + "#RESULTS\n", + "print 'Depth of the centre of pressure = %.2f ft from the base'%(y) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Depth of the centre of pressure = 3.82 ft from the base\n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.16 Page No : 23" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "d = 4. \t\t#ft\n", + "h = 6. \t\t#in\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = math.pi*d**2/4\n", + "X = (h-d)\n", + "I0 = (math.pi*d**4/64)+4*math.pi*(X)**2\n", + "h1 = I0/(A*X)\n", + "h2 = d-h1\n", + "\t\t\n", + "#RESULTS\n", + "print 'Depth of the axis be placed in order = %.1f ft '%(h2) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Depth of the axis be placed in order = 1.5 ft \n" + ] + } + ], + "prompt_number": 16 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.17 Page No : 24" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "h = 10 \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "x = math.sqrt(h**2/2)\n", + "\t\t\n", + "#RESULTS\n", + "print 'Depth of the axis be placed in order = %.2f ft '%(x) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Depth of the axis be placed in order = 7.07 ft \n" + ] + } + ], + "prompt_number": 17 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.18 Page No : 26" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "h = 8. \t\t#ft\n", + "h1 = 10. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = h\n", + "X = (h1/2)\n", + "Ig = h**3/12\n", + "I0 = Ig+A*X**2\n", + "h2 = I0/(A*X)\n", + "\t\t\n", + "#RESULTS\n", + "print 'depth at which the hinge of the shutter = %.2f ft '%(h2)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "depth at which the hinge of the shutter = 6.07 ft \n" + ] + } + ], + "prompt_number": 18 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.19 Page No : 27" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "from numpy import *\n", + "\t\t\n", + "#initialisation of variables\n", + "k1 = 1. \t\t#ft\n", + "k2 = 35.98 \t\t#ft\n", + "k3 = 66.83 \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "\n", + "vec =roots([k1,0,-k2,k3])\n", + "\n", + "X = vec[1]\n", + "\t\t\n", + "#RESULTS\n", + "print 'depth of the water = %.2f ft'%(X)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "depth of the water = 4.65 ft\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.22 Page No : 31" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 8. \t\t#ft\n", + "d1 = 2. \t\t#ft\n", + "h = 4. \t\t #ft\n", + "h1 = 2 \t\t #ft\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "\t\t\n", + "#CALCULATIONS\n", + "A1 = math.pi*d**2/4\n", + "A2 = math.pi*d1**2/4\n", + "A = A1-A2\n", + "x = (A1*d-A2*(d+h-h1))/A\n", + "P = w*A*x\n", + "Ig = ((math.pi*d**4/64)+(A1*(d-x)**2))-((math.pi*d1**4/64)+(A2*(h1+d-x)**2))\n", + "h2 = (Ig/(A*x))+x\n", + "\t\t\n", + "#RESULTS\n", + "print 'depth of the centre of the pressure = %.1f ft '%(h2)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "depth of the centre of the pressure = 8.4 ft \n" + ] + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.25 Page No : 34" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "W = 62.4 \t\t#lbs/ft**3\n", + "a = 140. \t\t#degrees\n", + "h = 20. \t\t#ft\n", + "w = 6. \t\t #ft\n", + "h1 = 17. \t\t#ft\n", + "h2 = 5. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "P1 = int(W*h1**2*w/2)\n", + "P2 = W*h2**2*w/2\n", + "P = P1-P2\n", + "y = (P1*(h1/3)-P2*(h2/3))/P\n", + "R = P/(2*math.sin(math.radians((180-a)/2)))\n", + "Rt = y*R/h\n", + "Rb = R-Rt\n", + "\n", + "#RESULTS\n", + "print 'Rt = %.f lbs '%(Rt)\n", + "print ' Rb = %.f lbs '%(Rb)\n", + "\n", + "# note : incorrect answer for R in the textbook. Hence, the difference in answers" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rt = 21838 lbs \n", + " Rb = 50409 lbs \n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.26 Page No : 36" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "w = 64. \t\t#lbs/ft**3\n", + "h = 12. \t\t#ft\n", + "l = 9. \t\t#ft\n", + "a = 45. \t\t#degrees\n", + "\t\t\n", + "#CALCULATIONS\n", + "P = w*h**2/2\n", + "h1 = h/3\n", + "Rb = P*h1/l\n", + "Ra = P-Rb\n", + "Wh = Rb*h1\n", + "T = Wh/math.sin(math.radians(a))\n", + "\n", + "\n", + "#RESULTS\n", + "print 'Load on the strut = %d lbs '%(T)\n", + "\n", + "# note : incorrect answer for T in the textbook" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Load on the strut = 11585 lbs \n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.27 Page No : 38" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "h = 9. \t\t#ft\n", + "l = 10. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "P = w*h**2/2\n", + "h1 = h/3\n", + "Ra = P/2\n", + "x = (w*4*h**2/9)/Ra\n", + "x1 = x+(h/3)\n", + "hb = h1-x\n", + "W = Ra*l\n", + "\t\t\n", + "#RESULTS\n", + "print 'magnitude od total in each beam = %d lbs '%(W)\n", + "\n", + "# note : rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "magnitude od total in each beam = 12636 lbs \n" + ] + } + ], + "prompt_number": 23 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/ch10.ipynb b/Hydraulics_Made_Easy_by_R._S._Dighe/ch10.ipynb new file mode 100755 index 00000000..0162251f --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/ch10.ipynb @@ -0,0 +1,497 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:2a1a58d10d8bc25ee316405092c20184a6b3d847d7641c4fa061749730145581" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 10 : Miscellaneous Problems" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.1 Page No : 316" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "#initialisation of variables\n", + "w = 62.4 \t\t#lb/ft**3\n", + "x = 8. \t\t#ft\n", + "A = 16. \t\t#ft**2\n", + "X = 2.5 \t\t#ft\n", + "X1 = 0.66 \t\t#ft\n", + "x1 = 3.834 \t\t#ft\n", + "x2 = 2.182 \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "P = w*x*A\n", + "y = A/3\n", + "P1 = w*x*A*0.5*X1\n", + "R = math.sqrt(P1**2+P**2)\n", + "m = P1/P\n", + "X2 = x1-x2\n", + "C = ((2./3)*A)-m*X\n", + "Y = m*X2+ C\n", + "print P1\t\n", + "#RESULTS\n", + "print 'Water pressure on vertical face = %.f lbs'%(round(P,-3))\n", + "print ' pressure which acts at the base = %.2f ft'%(y)\n", + "print ' Resultant = %.f lbs'%(R)\n", + "print ' x coordinate of the resultant = %.3f ft'%(X2)\n", + "print ' y coordinate of the resultant = %.3f ft'%(Y)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "2635.776\n", + "Water pressure on vertical face = 8000 lbs\n", + " pressure which acts at the base = 5.33 ft\n", + " Resultant = 8411 lbs\n", + " x coordinate of the resultant = 1.652 ft\n", + " y coordinate of the resultant = 10.387 ft\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.2 Page No : 319" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "s = 13.6\n", + "h = 12. \t\t#in\n", + "u = 0.04\n", + "k = 1.\n", + "d = 6. \t\t#in\n", + "g = 32.2 \t\t#ft/sec**2\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "\t\t\n", + "#CALCULATIONS\n", + "h1 = h*(s-1)/12\n", + "hf = u*h1\n", + "hn = h1-hf\n", + "Q = k*math.pi/4*(d/12)**2*8.02*math.sqrt(hn)/(math.sqrt(16-k))\n", + "Q = Q*60*w/10 # fro, cusecs to GPM\n", + "\n", + "#RESULTS\n", + "print 'discharge through flow = %.f ft G.P.M'%(Q)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "discharge through flow = 529 ft G.P.M\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.3 Page No : 321" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "za = 16. \t\t#ft\n", + "h1 = 2. \t\t#ft\n", + "h2 = 3. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "vc = math.sqrt(2*g*(za-h1-h2))\n", + "vb = vc*(h1/(2*h1))**2\n", + "r = -h1-h2-(vb**2/(2*g))\n", + "r1 = r+34\n", + "\t\t\n", + "#RESULTS\n", + "print 'pressure head at B = %.1f ft lb'%(r1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "pressure head at B = 28.3 ft lb\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.4 Page No : 322" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "g = 32.2 \t\t#ft/sec**2\n", + "Cd = 0.62\n", + "a = 90. \t\t#degrees\n", + "H1 = 14. \t\t#in\n", + "H2 = 8. \t\t#in\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q1 = (8./15)*Cd*math.sqrt(2*g)*math.tan(math.radians(a/2))*(H1/12)**(5/2.)\n", + "Q2 = (8./15)*Cd*math.sqrt(2*g)*math.tan(math.radians(a/2))*(H2/12)\n", + "Q = Q1-Q2\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge through notch = %.2f cuses'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge through notch = 2.13 cuses\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.5 Page No : 324" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "g = 32.2 \t\t#ft/sec**2\n", + "Cd = 0.62\n", + "d = 5./4 \t\t#in\n", + "h = 9. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "T = (2./3)*math.pi*(h)**(3./2)/(Cd*(math.pi/4)*math.sqrt(2*g)*(d/12)**2)\n", + "\t\t\n", + "#RESULTS\n", + "print 'time required to lower water level = %.f secs'%(T)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "time required to lower water level = 1334 secs\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.6 Page No : 325" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "a = 60. \t\t#degrees\n", + "d = 4. \t\t#in\n", + "Cd = 0.62\n", + "h = 5. \t \t#ft\n", + "w = 30. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "H1 = 10*math.sin(math.radians(a))\n", + "H2 = H1-h\n", + "T = (2*w/math.tan(math.radians(a)))*(2./3)*(H1**(3./2)-H2**(3./2))/(Cd*math.sqrt(2*g)*math.pi/(4*(d/12)**2))*100\n", + "\n", + "#RESULTS\n", + "print 'time required to lower water level = %.f secs'%(T)\n", + "\n", + "# answer is accurate.please check manually" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "time required to lower water level = 1214 secs\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.7 Page No : 326" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "p1 = 40. \t\t#percent\n", + "p2 = 35. \t\t#percent\n", + "dh = 200. \t\t#ft\n", + "f = 0.1\n", + "g = 32.2 \t\t#ft/sec**2\n", + "l = 2000. \t\t#ft\n", + "d = 1. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "hf1 = p1*dh/100\n", + "hf2 = p2*dh/100\n", + "hf3 = (100-p1-p2)*dh/100\n", + "hft = hf1+hf2+hf3\n", + "v1 = math.sqrt(2*g*hf1/(4*f*l))\n", + "Q = v1*math.pi*d**2/4\n", + "d2 = (Q*7*math.sqrt(3/(5*g)))**(2./3)\n", + "v3 = Q*4*(4./3)**2/math.pi\n", + "l3 = hf2*2*g*(3./4)/(4*f*v3**2)\n", + "\t\t\n", + "#RESULTS\n", + "print 'proportion of the quantity folwing in the bypass to the whole pass = %d ft'%(l3)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "proportion of the quantity folwing in the bypass to the whole pass = 415 ft\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.8 Page No : 328" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 1. \t \t#ft\n", + "l = 2000. \t\t#ft\n", + "f = 0.038\n", + "g = 32.2 \t\t#/ft/sec**2\n", + "Q = 6. \t\t#cuses\n", + "l1 = 1500. \t\t#ft\n", + "r = 2.\n", + "\t\t\n", + "#CALCULATIONS\n", + "v = 4*Q/(d**2*math.pi)\n", + "hf = 4*f*l*v**2/(2*g)\n", + "v1 = math.sqrt(hf*2*g/(4*f*l1+4*f*(l-l1)*r**2))\n", + "v3 = r*v1\n", + "Q1 = math.pi*d**2*v3/4\n", + "Q2 = math.pi*d**2*v1/4\n", + "r1 = Q2/Q1\n", + "\t\t\n", + "#RESULTS\n", + "print 'proportion of the quantity folwing in the bypass to the whole pass = %.1f '%(r1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "proportion of the quantity folwing in the bypass to the whole pass = 0.5 \n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.9 Page No : 329" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "f = 0.01\n", + "d = 3. \t\t#in\n", + "l = 22. \t\t#ft\n", + "l1 = 20. \t\t#ft\n", + "w = 20. \t\t#ft\n", + "h = 5. \t\t#ft\n", + "h1 = 20. \t\t#ft\n", + "t = 4. \t\t#min\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "h2 = h+h1\n", + "h3 = (h-(t*60*math.pi*math.sqrt(2*g/h)/(l1*w*2*64)))**2-4\n", + "dh = h2-h3\n", + "Q = dh*l1*w\n", + "\t\t\n", + "#RESULTS\n", + "print 'Quantiy discharged = %.f cuses '%(round(Q,-2))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Quantiy discharged = 1800 cuses \n" + ] + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.10 Page No : 332" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "g = 32.2 \t\t#ft/sec**2\n", + "sct = 1.6\n", + "sl = 0.8\n", + "K = 0.98\n", + "dh1 = 4. \t\t#ft\n", + "W = 62.4 \t\t#lbs/ft**3\n", + "d1 = 8. \t\t#in\n", + "d2 = 6. \t\t#in\n", + "\n", + "\t\t\n", + "#CALCULATIONS\n", + "dp = dh1*((sct/sl)-1)\n", + "C = math.sqrt(2*g)*math.pi*(d1/24)**2 /math.sqrt((d1**2/d2**2)**2 -1)\n", + "Q = C*K*math.sqrt(dh1)\n", + "\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge passing through the pipe = %.1f cuses '%(Q)\n", + "\t\t#The answer given in textbook is wrong. Please verify it.\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge passing through the pipe = 3.7 cuses \n" + ] + } + ], + "prompt_number": 10 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/ch2.ipynb b/Hydraulics_Made_Easy_by_R._S._Dighe/ch2.ipynb new file mode 100755 index 00000000..31d84882 --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/ch2.ipynb @@ -0,0 +1,368 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:1efa4f61c6638a7002643398d2d4bfe4efa140d30a6184e17976681fd0180ed8" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 2 : Floatation and Buoyancy" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.1 Page No : 43" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "l = 60. \t\t#ft\n", + "w = 10. \t\t#ft\n", + "h = 5. \t\t#ft\n", + "t = 3./16 \t\t#in\n", + "sp = 7.75\n", + "H = 4. \t \t#ft\n", + "w1 = 62.4 \t\t#lb/ft**3\n", + "y = 4. \t\t #ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "V = (l*w+2*w*h+2*l*h)*t/12\n", + "W = V*w1*sp\n", + "x = W/(w1*l*w)\n", + "W1 = H*l*w*w1\n", + "dW = (W1-W)/2238\n", + "\t\t\n", + "#RESULTS\n", + "print 'weight of water print laced = %.1f tons'%(dW)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "weight of water print laced = 62.5 tons\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.3 Page No : 50" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "D = 64. \t\t#lb/ft**3\n", + "d = 6. \t\t#ft\n", + "l = 10. \t\t#ft\n", + "W = 2. \t \t#tons\n", + "\t\t\n", + "#CALCULATIONS\n", + "V = W*2240/D\n", + "h = V/(math.pi*d**2/4)\n", + "BM = d**2/(16*h)\n", + "P = -(math.sqrt(64*BM*2*10*math.pi*(22400-math.pi*d**4))-W*22400)/10\n", + "\t\t\n", + "#RESULTS\n", + "print 'Minimum pull required = %.f lbs '%(P+3) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Minimum pull required = 3665 lbs \n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.4 Page No : 52" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "sg = 7.\n", + "sg1 = 5.\n", + "d = 8. \t\t#in\n", + "t = 1. \t\t#in\n", + "\t\t\n", + "#CALCULATIONS\n", + "x = (sg+sg1)+math.sqrt(d*(sg*(sg1+t)+1))\n", + "\t\t\n", + "#RESULTS\n", + "print 'maximum length of cylinder = %.2f in '%(x) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "maximum length of cylinder = 30.55 in \n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.7 Page No : 56" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "W = 2000. \t\t#tons\n", + "m = 15. \t\t#/tons\n", + "dx = 24. \t\t#ft\n", + "l = 3. \t \t#in\n", + "dx1 = 5. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "GM = m*dx/(W*(l/(dx1*12)))\n", + "\t\t\n", + "#RESULTSS\n", + "print 'metacentric height = %.1f ft '%(GM) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "metacentric height = 3.6 ft \n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.8 Page No : 56" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "M = 350. \t\t#tons\n", + "l = 50. \t\t#ft\n", + "w = 20. \t\t#ft\n", + "W = 100. \t\t#tons\n", + "h = 6. \t\t#ft\n", + "M1 = 250. \t\t#tons\n", + "\t\t\n", + "#CALCULATIONS\n", + "V = M*2240/64\n", + "d = V/(l*w)\n", + "BM = l*w**3/(12*w*l*d)\n", + "y = (((BM+(d/2))*(M/10))-(M1*h/10))/(W/10)\n", + "\t\t\n", + "#RESULTS\n", + "print 'Highest position of centre of gravity = %.2f ft '%(y)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Highest position of centre of gravity = 15.96 ft \n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.9 Page No : 58" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "W = 2000. \t\t#tons\n", + "l = 250. \t\t#ft\n", + "w = 30. \t\t#ft\n", + "a = 1./15\n", + "W1 = 50. \t\t#tons\n", + "h = 10. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "BG = (l*w**3*64/(W*2240*12))-(W1*h/(a*W))\n", + "\t\t\n", + "#RESULTS\n", + "print 'distance of the centre of gravity = %.2f ft '%(BG) \n", + "\n", + "# note : rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "distance of the centre of gravity = 4.29 ft \n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.10 Page No : 58" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "l = 91. \t\t#ft\n", + "w = 30. \t\t#ft\n", + "h = 6. \t \t#ft\n", + "W = 40. \t\t#tons\n", + "a = 3. \t\t #degrees\n", + "cg = 3. \t\t#ft\n", + "d = 4. \t\t#ft\n", + "W1 = 60. \t\t#tons\n", + "cg1 = 1. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "W2 = (l*w*d*64/2240)-W1\n", + "y = (W2*(h/2)+W1*(cg+d))/(l*w*d*64/2240)\n", + "BG = y-(d/2)\n", + "BM = l*w**3/(12*l*w*d)\n", + "GM = BM-BG\n", + "dx = GM*l*w*d*64*math.tan(math.radians(a))/(60*2240)\n", + "\t\t\n", + "#RESULTS\n", + "print 'maximum distance through which the load can be shifted = %.1f ft '%(dx)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "maximum distance through which the load can be shifted = 4.6 ft \n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.11 Page No : 60" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "W = 5000. \t\t#tons\n", + "I = 1.4*10**6 \t\t#ft**4\n", + "k = 12.2 \t\t#ft\n", + "BG = 6.5 \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "BM = I*64/(W*2240)\n", + "GM = BM-BG\n", + "T = 2*math.pi*math.sqrt(k**2/(GM*32.2))\n", + "\t\t\n", + "#RESULTS\n", + "print 'period of oscialltion = %.2f sec '%(T) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "period of oscialltion = 11.03 sec \n" + ] + } + ], + "prompt_number": 8 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/ch3.ipynb b/Hydraulics_Made_Easy_by_R._S._Dighe/ch3.ipynb new file mode 100755 index 00000000..a7c97aa0 --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/ch3.ipynb @@ -0,0 +1,527 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:2f07cf575072954f2c364fa1bc9af38d9d9d43725d39198530552558da11415c" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 3 : Flow of Water" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.1 Page No : 67" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "d1 = 1. \t\t#ft\n", + "d2 = 6. \t\t#in\n", + "h1 = 5. \t\t#ft\n", + "h2 = 15. \t\t#ft\n", + "Pa = 15. \t\t#lbs\n", + "v1 = 10. \t\t#ft/sec\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "v2 = v1/(d2/12)**2\n", + "Pb = (w*((Pa+(Pa*144/w)+(v1**2/(2*g)))-h1-(v2**2/(2*g))))/144\n", + "\t\t\n", + "#RESULTS\n", + "print 'Pb = %.2f lbs/in**2 '%(Pb) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Pb = 9.24 lbs/in**2 \n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.2 Page No : 69" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d1 = 4. \t\t#ft\n", + "d2 = 2. \t\t#ft\n", + "h1 = 50. \t\t#ft\n", + "h2 = 45. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "r = (d1**2/d2**2)\n", + "v1 = round(math.sqrt((h1-h2)*2*g/(r**2-1)),1)\n", + "Q = v1*math.pi*d1**2/4\n", + "\n", + "#RESULTS\n", + "print 'discharge through pipe = %.2f cubic feet per second '%(Q)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "discharge through pipe = 57.81 cubic feet per second \n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.3 Page No : 70" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "z1 = 10. \t\t#/m\n", + "h1 = 10. \t\t#m\n", + "v1 = 12. \t\t#ft/sec\n", + "v2 = 4. \t\t#m/sec\n", + "k = 0.6\n", + "w = 62.4 \t\t#lb/in**2\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "p = (w/144)*(z1+h1+(v1**2/(2*g))-(v2**2/(2*g))-(k*(v1-v2)**2/(2*g)))\n", + "\t\t\n", + "#RESULTS\n", + "print 'pressure at bottom end = %.2f lb/in**2'%(p)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "pressure at bottom end = 9.27 lb/in**2\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.4 Page No : 73" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 4. \t\t#ft\n", + "d1 = 5./4 \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "h = 3. \t \t#ft\n", + "K = 1.\n", + "\t\t\n", + "#CALCULATIONS\n", + "C = (math.pi/4)*d**2*math.sqrt(2*g)/(math.sqrt((d**2/d1**2)**2-1))\n", + "Q = K*math.sqrt(h)*C\n", + "V = Q/(math.pi*d1**2/4)\n", + "\t\t\n", + "#RESULTS\n", + "print 'Velocity at the throat = %.2f ft/sec '%(V)\n", + "\n", + "# ronding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Velocity at the throat = 13.97 ft/sec \n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.5 Page No : 74" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 9. \t\t#in\n", + "d1 = 4. \t\t#in\n", + "g = 32.2 \t\t#ft/sec**2\n", + "dh = 10. \t\t#in\n", + "sg = 13.6 \n", + "K = 1.\n", + "\t\t\n", + "#CALCULATIONS\n", + "C = (((math.pi/4)**2*(d*d1)**2*math.sqrt(2*g)/144**2)/(math.sqrt((math.pi*d**2/12**2)**2-(math.pi*d1**2/12**2)**2)))+0.52\n", + "h = (sg-1)*dh/12\n", + "Q = K*C*math.sqrt(h)\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge passing through the pipe = %.2f cuses '%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge passing through the pipe = 2.26 cuses \n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.6 Page No : 76" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "sm = 13.6\n", + "so = 0.8\n", + "di = 8. \t\t#in\n", + "dt = 4. \t\t#in\n", + "K = 0.98\n", + "v = 1. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "s = sm/so\n", + "dp = v*12*(s-1)/12\n", + "A = math.pi*(di/12)**2/4\n", + "At = math.pi*(dt/12)**2/4\n", + "C = A*math.sqrt(2*g)/(math.sqrt((A/At)**2-1))\n", + "Q = C*math.sqrt(v*12+dt)*K\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge passing through the pipe = %.2f cuses '%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge passing through the pipe = 2.84 cuses \n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.7 Page No : 77" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "s = 1./10\n", + "d1 = 6. \t\t#in\n", + "d2 = 2. \t\t#in\n", + "l = 20. \t\t#in\n", + "p = 15. \t\t#lbs/in**2\n", + "p1 = 6. \t\t#lbs/in**2\n", + "K = 0.95\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "H = (l*s/12)-(p1*144/(2*g))+(p*144/(2*g))\n", + "C = math.sqrt(2*g)*(math.pi*(d1/12)**2)/(4*(math.sqrt((d1**2/d2**2)**2-1)))\n", + "Q = C*K*math.sqrt(H)*374.7\n", + "\n", + "#RESULTS\n", + "print 'Discharge passing through the pipe = %.f Gallons/minute '%(Q)\n", + "\n", + "# note : rounding off error at value of H in textbook. so answer is slightly different" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge passing through the pipe = 282 Gallons/minute \n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.8 Page No : 78" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d1 = 12. \t\t#in\n", + "Q = 4.25 \t\t#ft**3/sec\n", + "h = 18. \t\t#ft\n", + "K = 0.98\n", + "g = 32.2 \t\t#ft/sec**2\n", + "sm = 13.6\n", + "\t\t\n", + "#CALCULATIONS\n", + "R = math.sqrt((K*math.sqrt(2*g)*math.sqrt(h)*(math.pi*(d1/12)**2/4)/Q)+1)\n", + "d2 = math.sqrt(d1**2/(144*R))\n", + "dh = (sm-1)*(h/(12*2))\n", + "d3 = Q*math.sqrt(dh/h)\n", + "\t\t\n", + "#RESULTS\n", + "print 'Diameter of the throat = %.2f ft '%(d3)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Diameter of the throat = 3.08 ft \n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.9 Page No : 81" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "R = 4. \t\t#in\n", + "r = 0.5 \t\t#in\n", + "c = 0.007\n", + "K = 33.96\n", + "w = 62.4 \t\t#lb/ft**3\n", + "pa = 12.13 \t\t#lb/in**2\n", + "pb = 14.7 \t\t#lb/in**2\n", + "w1 = 2.5 \t\t#lbs\n", + "Q = 40. \t\t#gals/min\n", + "h = 1.86\n", + "\t\t\n", + "#CALCULATIONS\n", + "va = Q*4*(2*r*12)**2/(6*w*math.pi)\n", + "vb = Q*(2*r*12)**2/(6*w*2*R*math.pi*0.32)\n", + "vx = vb*R/2\n", + "pu = 2*math.pi*w*h\n", + "pd = pb*math.pi*R**2\n", + "RP = pb*math.pi*R**2-2*math.pi*w*(0.5*K*((R/12)**2-(r/12)**2)-c*math.log(R/r))-pa*math.pi*r**2+w1\n", + "\t\t\n", + "#RESULTS\n", + "print 'velocity va = %.1f ft/sec'%(va)\n", + "print 'velocity vb = %.2f ft/sec'%(vb)\n", + "print 'velocity vx = %.2f ft/sec'%(vx)\n", + "print 'pressure px = %.1f lbs/in**2'%(pb)\n", + "print 'upward pressure = %.1f lbs'%(pu)\n", + "print 'downward pressure = %.1f lbs'%(pd)\n", + "print 'Resultant pressure = %.1f lbs'%(RP)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "velocity va = 19.6 ft/sec\n", + "velocity vb = 1.91 ft/sec\n", + "velocity vx = 3.83 ft/sec\n", + "pressure px = 14.7 lbs/in**2\n", + "upward pressure = 729.3 lbs\n", + "downward pressure = 738.9 lbs\n", + "Resultant pressure = 9.4 lbs\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.10 Page No : 86" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 1. \t\t#ft\n", + "h = 4. \t \t#ft\n", + "h1 = 3. \t\t#ft\n", + "p = 25. \t\t#percent\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "h2 = ((h/4)-(h1/4))*h*2\n", + "w = math.sqrt(h2*2*g/(d/2)**2)\n", + "N = w*60/(2*math.pi)\n", + "h3 = (h-h1**2/4)*2\n", + "w1 = math.sqrt(h3*2*g/(d/2)**2)\n", + "N1 = w1*60/(2*math.pi)\n", + "\t\t\n", + "#RESULTS\n", + "print 'original volume = %.1f R.P.M '%(N1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "original volume = 286.7 R.P.M \n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.12 Page No : 89" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "R2 = 2. \t\t#ft\n", + "R1 = 1. \t\t#ft\n", + "w = 200. \t\t#r.p.m\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "v2 = R2*math.pi*w*R2/60\n", + "v1 = R2*math.pi*w*R1/60\n", + "H = (v2**2-v1**2)/(2*g)\n", + "\t\t\n", + "#RESULTS\n", + "print 'centrifugal head = %.1f ft of watrer '%(H)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "centrifugal head = 20.4 ft of watrer \n" + ] + } + ], + "prompt_number": 11 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/ch4.ipynb b/Hydraulics_Made_Easy_by_R._S._Dighe/ch4.ipynb new file mode 100755 index 00000000..0747f48b --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/ch4.ipynb @@ -0,0 +1,994 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:98b9b2cc3c2bb837908ed50b11e73842e607c874a3e8f477bc1209f33ea8918d" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 4 : Flow of Water Through Orifices and Mouthpieces" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.1 Page No : 98" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "M = 31*10 \t\t#lbs\n", + "P = 3.6 \t\t#lbs\n", + "t = 60. \t\t#sec\n", + "g = 32.2 \t\t#ft/sec**2\n", + "H = 9. \t\t#ft\n", + "d = 1. \t\t#in\n", + "w = 6.24 \t\t#gallons\n", + "\t\t\n", + "#CALCULATIONS\n", + "v = P*g*t/M\n", + "V = math.sqrt(2*g*H)\n", + "Cv = v/V\n", + "V1 = math.pi*(d/12)**2*V*60*w/4\n", + "Cd = M/(10*V1)\n", + "Cc = Cd/Cv\n", + "Cr = (1/Cv**2)-1\n", + "\t\t\n", + "#RESULTS\n", + "print 'Coefficient of resistance = %.2f '%(Cr)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Coefficient of resistance = 0.15 \n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.2 Page No : 100" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "M = 1.65 \t\t#lbs\n", + "Q = 31. \t\t#gallons per min\n", + "d = 1. \t\t#in\n", + "h = 4. \t\t#ft\n", + "t = 60. \t\t#sec\n", + "g = 32.2 \t\t#ft/sec**2\n", + "Q1 = 6.24 \t\t#gallons per min\n", + "c = 0.36\n", + "P = 3.6 # lbs\n", + "H = 9. \n", + "\n", + "#CALCULATIONS\n", + "v = P*g*t/(Q*10)\n", + "V = math.sqrt(2*g*H)\n", + "Cv = (v/V)\n", + "vf = V*math.pi*(d/12)**2*60*Q1/4\n", + "Cd = Q/vf\n", + "Cc = Cd/Cv\n", + "Cr = (1/Cv**2)-1\n", + "\n", + "#RESULTS\n", + "print 'velocity of jet = %.2f ft/sec'%(v)\n", + "print 'theoretical velocity of jet = %.2f ft/sec'%(V)\n", + "print 'Cv = %.2f '%(Cv)\n", + "print 'volume flow = %.2f gallons per minute'%(vf)\n", + "print 'Cd = %.2f '%(Cd)\n", + "print 'Cc = %.2f '%(Cc)\n", + "print 'Coefficient of resistance = %.2f '%(Cr)\n", + "\n", + "# Note : Answer for theoretical velocity is wrong in book. Please check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "velocity of jet = 22.44 ft/sec\n", + "theoretical velocity of jet = 24.07 ft/sec\n", + "Cv = 0.93 \n", + "volume flow = 49.16 gallons per minute\n", + "Cd = 0.63 \n", + "Cc = 0.68 \n", + "Coefficient of resistance = 0.15 \n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.3 Page No : 101" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "x = 11.5 \t\t#in\n", + "y = 1.2 \t\t#in\n", + "H = 29. \t\t#in\n", + "q = 6.24 \t\t#gallons per minute\n", + "d = 1. \t\t#in\n", + "g = 32.2 \t\t#ft/sec**2\n", + "Q = 16. \t\t#gallons per min\n", + "\t\t\n", + "#CALCULATIONS\n", + "Cv = math.sqrt(x**2/(4*H*y))\n", + "Q1 = math.pi*(d/12)**2*math.sqrt(2*g*H/12)*q*60/4\n", + "Cd = Q/Q1\n", + "Cc = Cd/Cv\n", + "Cr = (1/Cv**2)-1\n", + "\t\t\n", + "#RESULTS\n", + "print 'Coefficient of resistance = %.2f '%(Cr)\n", + "\n", + "# rounding off error." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Coefficient of resistance = 0.05 \n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.4 Page No : 102" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "x = 3.2 \t\t#ft\n", + "d = 8. \t\t#ft\n", + "W = 5.12 \t\t#lb\n", + "A = 1./144\n", + "H = 4. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "Q = 251.5 \t\t#lbs/min\n", + "w = 62.4 \t\t#lbs/ft**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "F = W*x/d\n", + "v = W*x*g*60/(d*Q)\n", + "V = math.sqrt(2*g*H)\n", + "Cv = v/V\n", + "Q1 = A*V*60*w\n", + "Cd = Q/Q1\n", + "Cc = Cd/Cv\n", + "\t\t\n", + "#RESULTS\n", + "print 'Cc = %.2f '%(Cc)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Cc = 0.61 \n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.5 Page No : 106" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 8. \t\t#in\n", + "\t\t\n", + "#CALCULATIONS\n", + "Cd = 1/math.sqrt(1+((1./(8.**2/100)))-1)\n", + "area = math.pi/4 * (2./12)**2\n", + "Discharge = area * Cd * math.sqrt(2*32.2*4) \t\t\n", + "\n", + "#RESULTS\n", + "print 'Cd = %.2f '%(Cd)\n", + "print \"Discharge = %.1f cubic ft./sec.\"%Discharge\n", + "\n", + "# note : rounding off error. please check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Cd = 0.80 \n", + "Discharge = 0.3 cubic ft./sec.\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.6 Page No : 109" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d =2. \t\t#in\n", + "h = 6. \t \t#ft\n", + "H = 26. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "R = 6.\n", + "\t\t\n", + "#CALCULATIONS\n", + "v2 = math.sqrt(2*g*(H+h))\n", + "Q = math.pi*(d/12)**2*v2/4\n", + "v3 = math.sqrt(2*g*h)\n", + "r = v2/v3\n", + "d3 = math.sqrt(r*d**2)\n", + "v4 = math.sqrt(v2**2/R)\n", + "d4 = math.sqrt(d**2*(v2/v4))\n", + "\t\t\n", + "#RESULTS\n", + "print 'diameter = %.2f in'%(d4)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "diameter = 3.13 in\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.7 Page No : 111" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "r = 9./16\n", + "r1 = 7./16\n", + "h = 26. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "r2 = 1/((r**2)+(0.25*r1**2))\n", + "H1 = h/(r2-1)\n", + "\t\t\n", + "#RESULTS\n", + "print 'maximu head of the tank = %.3f ft of water'%(H1)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "maximu head of the tank = 14.897 ft of water\n" + ] + } + ], + "prompt_number": 17 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.8 pageno : 114" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# variables\n", + "A = 30.*15 # sq ft\n", + "a = 2. # sq ft\n", + "H1 = 5. # ft\n", + "H2 = 0. \n", + "\n", + "# calculation\n", + "T = a*A*H1**(1./2)/(.62*a*8.02)\n", + "\n", + "# result\n", + "print \"Time of emptying pool : T = %.1f seconds\"%T\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Time of emptying pool : T = 202.4 seconds\n" + ] + } + ], + "prompt_number": 21 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.9 Page No : 115" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "H1 = 9. \t\t#ft\n", + "A = 2. \t\t#ft**2\n", + "H2 = 4. \t\t#ft\n", + "d = 2.25 \t\t#in\n", + "t = 60. \t\t#sec\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "a = (d/12)**2\n", + "Cd = (2*A*(H1-H2)**0.5)/(t*a*math.sqrt(2*g))\n", + "\t\t\n", + "#RESULTS\n", + "print 'coefficient of dicharge = %.3f '%(Cd)\n", + "\n", + "\n", + "#ANSWER GIVEN IN THE TEXTBBOK IS WRONG..VERIFIED WITH CALCULATOR\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "coefficient of dicharge = 0.528 \n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.10 Page No : 115" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 1. \t\t#ft\n", + "h1 = 10. \t\t#ft\n", + "h2 = 2. \t\t#ft\n", + "Cd = 0.6\n", + "g = 32.2 \t\t#ft/sec**2\n", + "t = 12.6\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = math.pi*d**2/4\n", + "a = 1./144\n", + "T1 = (A/(a*Cd*math.sqrt(2*g)))*(1./3)*(h1**1.5-(h1-h2)**1.5-h2**1.5)+t\n", + "T2 = 2*A*(h2**0.5)/(Cd*a*math.sqrt(2*g))\n", + "T = T1+T2\n", + "\t\t\n", + "#RESULTS\n", + "print 'Total time = %.2f sec'%(T)\n", + "\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Total time = 127.32 sec\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.11 Page No : 117" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "from scipy.integrate import quad\n", + "\n", + "#initialisation of variables\n", + "l = 600. \t\t#ft\n", + "w = 400. \t\t#ft\n", + "s = 1.\n", + "h = 20. \t\t#ft\n", + "d = 3. \t \t#ft\n", + "dh = 10. \t\t#ft\n", + "Cd = 0.7\n", + "g = 32.2 \t\t#ft/sec**2\n", + "k = 240000.\n", + "k1 = 2000.\n", + "k2 = 4.\n", + "\t\t\n", + "#CALCULATIONS\n", + "def f(x):\n", + " return (k/math.sqrt(x) + k1*math.sqrt(x) + k2*x**(3./2))\n", + "\n", + "T1 = 1./(Cd * math.pi/4 * 9 * 8.02) * quad(f,10,20)[0]\n", + "\t\t\n", + "#RESULTS\n", + "print 'Time taken for 10 feet fall = %.f sec'%(T1)\n", + "\n", + "# note : quad() gives accurate answer. so answer is slightly different." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Time taken for 10 feet fall = 17846 sec\n" + ] + } + ], + "prompt_number": 31 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.12 Page No : 118" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Cd = 0.6\n", + "H1 = 8. \t\t#ft\n", + "H2 = 3. \t\t#ft\n", + "l = 90. \t\t#ft\n", + "b = 30. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "A = 2. \t\t#ft**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "T1 = 2*l*b*(H1**0.5-(H1-H2)**0.5)/(Cd*math.sqrt(2*g)*A)\n", + "T2 = (l*b*2/10)*(2./3)*(H1-H2)**1.5/(Cd*math.sqrt(2*g)*A)\n", + "T = T1+T2\n", + "\t\t\n", + "#RESULTS\n", + "print 'Time it take to emptify the swimming bath = %.1f sec'%(T)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Time it take to emptify the swimming bath = 750.1 sec\n" + ] + } + ], + "prompt_number": 32 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.13 Page No : 120" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Cd = 0.8\n", + "g = 32.2 \t\t#f/sec**2\n", + "d = 3. \t\t#in\n", + "x = 6. \t \t#ft\n", + "l = 25. \t\t#ft\n", + "d1 = 8. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = math.pi*(d/12)**2/4\n", + "T = (2*l/(Cd*A*math.sqrt(2*g)))*(-2./3)*((d1-x)**1.5-d1**1.5)\n", + "\t\t\n", + "#RESULTS\n", + "print 'Time it take to emptify the boiler = %.f sec'%(T+6)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Time it take to emptify the boiler = 2100 sec\n" + ] + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.14 Page No : 125" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "l = 30. \t\t#ft\n", + "w = 10. \t\t#ft\n", + "d = 4. \t\t#in\n", + "h1 = 10. \t\t#ft\n", + "h2 = 2. \t\t#ft\n", + "Cd = 0.97\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "A1 = w*3*l/4\n", + "A2 = l*w/4\n", + "A = math.pi*(d/12)**2/4\n", + "T = 2*A1*(math.sqrt(h1)-math.sqrt(h2))*10/(Cd*A*math.sqrt(2*g)*(l+w))\n", + "\t\t\n", + "#RESULTS\n", + "print 'Time it take to reduce the height = %.f sec'%(round(T,-1))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Time it take to reduce the height = 290 sec\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.15 Page No : 125" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "A1 = 1000. \t\t#ft**2\n", + "A2 = 1000. \t\t#ft**2\n", + "a = 2. \t\t#ft**2\n", + "H1 = 9. \t\t#ft\n", + "H2 = 4. \t\t#ft\n", + "Cd =0.8\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "T = a*1000*(math.sqrt(H1)-math.sqrt(H2))/(Cd*a**2*math.sqrt(2*g))\n", + "\t\t\n", + "#RESULTS\n", + "print 'Time it take to reduce the height = %.1f sec'%(T)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Time it take to reduce the height = 77.9 sec\n" + ] + } + ], + "prompt_number": 35 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.16 Page No : 128" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "l = 70. \t\t#ft\n", + "b = 10. \t\t#ft\n", + "Hl = 10. \t\t#ft\n", + "H1 = 6. \t\t#ft\n", + "h1 = 4. \t\t#ft\n", + "h2 = 2. \t\t#ft\n", + "w = 2. \t\t#ft\n", + "h3 = 3. \t\t#ft\n", + "Cd = 0.6\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "t = (l*b)*(Hl+H1)/(Cd*h2*w*h1*math.sqrt(2*g*H1))\n", + "t1 = 2*l*b*math.sqrt(Hl)/(Cd*h2*w*h3*math.sqrt(2*g))\n", + "\t\t\n", + "#RESULTS\n", + "# 2nd ans is wrong in book\n", + "print 'Time of filling = %.2f sec'%(t)\n", + "print ' Time of emptying = %.2f sec'%(t1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Time of filling = 59.35 sec\n", + " Time of emptying = 76.62 sec\n" + ] + } + ], + "prompt_number": 16 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.17 Page No : 129" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "from sympy import Symbol,solve\n", + "\t\t\n", + "#initialisation of variables\n", + "HL = 12.5 \t\t#ft\n", + "H1 = 10.5 \t\t#ft\n", + "Cd = 0.62\n", + "h = 4. \t\t#ft\n", + "l = 3. \t \t#ft\n", + "n = 2.\n", + "t = 5. \t\t #min\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "a1 = n*l*l\n", + "A = Symbol('A')\n", + "ans = solve( (2*A/(Cd*a1*26)) + 2*A*math.sqrt(H1)/(Cd*a1*8.02) - 300 )\n", + "A = ans[0]\n", + "\n", + "#RESULTS\n", + "print 'Area = %.f sq ft'%(A)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Area = 3783 sq ft\n" + ] + } + ], + "prompt_number": 41 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.18 Page No : 130" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Cd = 0.62\n", + "g = 32.2 \t\t#ft/sec**2\n", + "l = 200. \t\t#ft\n", + "w = 25. \t\t#ft\n", + "a1 = 5. \t\t#ft**2\n", + "h = 20. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "t = 2*l*w*math.sqrt(h-(h/a1))/(Cd*math.sqrt(2*g)*a1)\n", + "\t\t\n", + "#RESULTS\n", + "print 'tme rquired to fill the lock = %.f sec'%(t)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "tme rquired to fill the lock = 1608 sec\n" + ] + } + ], + "prompt_number": 42 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.19 Page No : 130" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "L = 150. \t\t#ft\n", + "w = 20. \t\t#ft\n", + "t = 5. \t\t#min\n", + "h = 5. \t \t#ft\n", + "Cd = 0.6 \n", + "Hl = 9. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "T = 2*L*w*math.sqrt(Hl)/(Cd*t*60*math.sqrt(2*g))\n", + "\t\t\n", + "#RESULTS\n", + "print 'Area of sumberged slice = %.1f sq ft'%(T)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Area of sumberged slice = 12.5 sq ft\n" + ] + } + ], + "prompt_number": 19 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.20 Page No : 132" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "L = 3. \t\t#ft\n", + "H1 = 1.5 \t\t#ft\n", + "H2 = 0.75 \t\t#ft\n", + "Cd = 0.62\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q = 2*Cd*60*L*math.sqrt(2*g)*(H1**1.5-H2**1.5)/3\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge per minute = %.1f cubic ft per minute'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge per minute = 709.1 cubic ft per minute\n" + ] + } + ], + "prompt_number": 20 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.21 Page No : 134" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Cd = 0.62\n", + "H1 = 6. \t\t#ft\n", + "H2 = 3. \t\t#ft\n", + "H = 4. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q1 = 2*Cd*H*math.sqrt(2*g)*(H**1.5-H2**1.5)/3\n", + "Q2 = Cd*H*(H1-H)*math.sqrt(2*g*H)\n", + "Q = Q1+Q2\n", + "\t\t\n", + "#RESULTS\n", + "print 'Total discharge = %.f cuses'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Total discharge = 117 cuses\n" + ] + } + ], + "prompt_number": 21 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/ch5.ipynb b/Hydraulics_Made_Easy_by_R._S._Dighe/ch5.ipynb new file mode 100755 index 00000000..76711f92 --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/ch5.ipynb @@ -0,0 +1,668 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:76799d83bd91128bf807e32787abef19090d6630d45c16a130eb1ab2178963ed" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 5 : Flow of Water Over Weirs" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.1 Page No : 141" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "L = 6. \t\t#ft\n", + "H = 15. \t\t#in\n", + "Cd = 0.62\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\n", + "#CALCULAIONS\n", + "Q = 2*Cd*L*math.sqrt(2*g)*(H/12)**1.5/3\n", + "\t\t\n", + "#RESULTS\n", + "print 'Total Discharge = %.1f cuses'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Total Discharge = 27.8 cuses\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.2 Page No : 143" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "o = 90. \t\t#degrees\n", + "H = 15.5 \t\t#in\n", + "Cd = 0.6\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q = 8*Cd*math.tan(math.radians(o/2))*math.sqrt(2*g)*(H/12)**2.5/15\n", + "\t\t\n", + "#RESULTS\n", + "print 'Total Discharge = %.2f cuses'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Total Discharge = 4.87 cuses\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.3 Page No : 143" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Cd = 0.62\n", + "L = 4. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "H = 6. \t\t#in\n", + "o = 90. \t\t#degrees\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q = Cd*L*math.sqrt(2*g)*(H/12)**1.5*(2./3)\n", + "H1 = (Q*15/(8*Cd*math.tan(math.radians(o/2))*math.sqrt(2*g)))**(2./5)\n", + "\t\t\n", + "#RESULTS\n", + "print 'depth of water = %.2f ft'%(H1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "depth of water = 1.26 ft\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.4 Page No : 144" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Cd = 0.62\n", + "L = 3. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "H = 1. \t\t#ft\n", + "L1 = 2. \t\t#ft\n", + "h = 0.5 \t\t#ft\n", + "L2 = 1. \t\t#ft\n", + "h1 = 0.25 \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q = 2*Cd*L*math.sqrt(2*g)*H**1.5/3\n", + "Q1 =2*Cd*L1*math.sqrt(2*g)*((H+h)**1.5- H**1.5)/3\n", + "Q2 = 2*Cd*L2*math.sqrt(2*g)*((H+h+h1)**1.5- (H+h)**1.5)/3\n", + "Q3 = Q1+Q2+Q\n", + "\t\t\n", + "#RESULTS\n", + "print 'Total Discharge = %.2f cuses'%(Q3)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Total Discharge = 17.09 cuses\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5 Page No : 149" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "h = 9. \t\t#in\n", + "l = 6. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "H = h/12\n", + "Q = math.sqrt(2*g)*l*(H/12)**1.5*(0.405+(0.00984/0.75))\n", + "Q1 = 3.33*l*H**1.5\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge by francis formula = %.2f cuses'%(Q1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge by francis formula = 12.98 cuses\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.6 Page No : 149" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "l = 24. \t\t#ft\n", + "n = 5. \t\t#parts\n", + "h = 2. \t\t#ft\n", + "w = 1.\t\t#ft\n", + "n1 = 4.\n", + "c = 10.\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q = 3.33*((l-n1)-0.1*c*h)*h**1.5\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge = %.1f cuses'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge = 169.5 cuses\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.7 Page No : 150" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "A = 25. \t\t#miles**2\n", + "t = 24. \t\t#hr\n", + "p = 50. \t\t#per cent\n", + "l = 3. \t\t#in\n", + "h = 4. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "A1 = 5280**2*A\n", + "V = A1*l/12\n", + "V1 = V/(t*60*60)\n", + "V2 = V1/2\n", + "L = (V2/(3.33*h*2))+0.2*4\n", + "\t\t\n", + "#RESULTS\n", + "print 'length of weir = %.1f ft'%(L)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "length of weir = 38.7 ft\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.8 Page No : 151" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "h = 4. \t\t#ft\n", + "w = 5. \t\t#ft\n", + "l = 2. \t\t#ft\n", + "Q1 = 1008.5 \t\t#cuses\n", + "n = 8. \t\t#piers\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q = 3.33*(w-0.2*h)*h**1.5\n", + "n1 = Q1/Q\n", + "L = n*l+w*n1\n", + "\t\t\n", + "#RESULTS\n", + "print 'length of weir = %.f ft'%(L)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "length of weir = 61 ft\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.9 Page No : 153" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "k = 3.33\n", + "l = 10. \t\t#ft\n", + "x = 2. \t\t#ft\n", + "A = 30. \t\t#ft**2\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q = k*(l-0.2*x)*x**1.5\n", + "V = Q/A\n", + "h = V**2/(2*g)\n", + "Q1 = k*(l-0.2*(x+h))*((x+h)**1.5-h**1.5)\n", + "va = Q1/A\n", + "ha = va**2/(2*g)\n", + "Q2 = k*(l-0.2*(x+ha))*((x+ha)**1.5-ha**1.5)\n", + "\t\t\n", + "#RESULTS\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge in franccis formula = %.2f cusecs'%(Q1)\n", + "print ' Discharge in corrected franccis formula = %.2f cusecs'%(Q2)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge in franccis formula = 98.17 cusecs\n", + " Discharge in corrected franccis formula = 99.41 cusecs\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.10 Page No : 154" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Cd = 0.6\n", + "g = 32.2 \t\t#ft/sec**2\n", + "o = 90. \t\t#degrees\n", + "H = 2. \t\t#ft\n", + "A = 15.2 \t\t#ft**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q = 8*Cd*math.sqrt(2*g)*math.tan(math.radians(o/2))*H**2.5/15\n", + "va = Q/A\n", + "ha = va**2/(2*g)\n", + "Q1 = 8*Cd*math.sqrt(2*g)*((H+ha)**2.5-ha**2.5)/15\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge of stream = %.1f cuses'%(Q1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge of stream = 14.8 cuses\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.11 Page No : 155" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "va = 4. \t\t#ft/sec\n", + "g = 32.2 \t\t#ft/sec**2\n", + "H = 1.25 \n", + "l = 10. \t\t#ft\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "p = 60. \t\t#per cent\n", + "l1 = 90. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "ha = va**2/(2*g)\n", + "Q = 3.333*(l-0.1*2*(H+ha))*((H+ha)**1.5-ha**1.5)*w\n", + "E = Q*l1\n", + "HP = E*60/(100*550)\n", + "\t\t\n", + "#RESULTS\n", + "print 'H.P available = %.1f H.P'%(HP)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "H.P available = 338.8 H.P\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.12 Page No : 157" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "L = 8. \t\t#ft\n", + "d = 9. \t\t#in\n", + "h = 3. \t\t#in\n", + "Cd1 = 0.62\n", + "Cd2 = 0.62\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q1 = (2./3)*Cd1*L*math.sqrt(2*g)*(h/12)**1.5\n", + "Q2 = Cd2*L*d*math.sqrt(2*g*h/12)/12\n", + "Q = Q1+Q2\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge = %.2f cuses'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge = 18.24 cuses\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.13 Page No : 158" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "L = 50. \t\t#ft\n", + "d = 2. \t\t#ft\n", + "h = 4. \t\t#ft\n", + "Cd1 = 0.58\n", + "Cd2 = 0.8\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "ha = h/(2*g)\n", + "Q1 = (2./3)*Cd1*L*math.sqrt(2*g)*((h+ha)**1.5-ha**1.5)\n", + "Q2 = Cd2*L*d*math.sqrt(2*g*(h+ha))\n", + "Q = Q1+Q2\n", + "\n", + "\n", + "#RESULTS\n", + "print 'Discharge = %d cuses'%(Q)\n", + "\n", + "# note : value of ha is calculated wrongly. please check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge = 2561 cuses\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.14 Page No : 158" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "M = 60.\n", + "k = 500.\n", + "v = 8. \t\t#ft/sec\n", + "w = 100. \t\t#ft\n", + "h1 = 5. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "x = 1.95 \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q = k*M**(2./3)\n", + "A = Q/v\n", + "md = A/w\n", + "h = md-h1\n", + "ha = v**2/(2*g)\n", + "H = h+x**2-1+h1-1\n", + "\t\t\n", + "#RESULTS\n", + "print 'height above the crest of the air = %.2f ft of water'%(H)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "height above the crest of the air = 11.38 ft of water\n" + ] + } + ], + "prompt_number": 14 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.16 Page No : 163" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "H2 = 1.5 \t\t#ft\n", + "H1 = 1. \t\t#ft\n", + "A = 100. \t\t#yards**2\n", + "Cd = 0.6\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "A1 = A*9\n", + "T = (1.25*A1/(Cd*math.sqrt(2*g)))*(H1-(1/H2)**1.5)\n", + "\t\t\n", + "#RESULTS\n", + "print 'time of lowering the surface = %.1f sec'%(T)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "time of lowering the surface = 106.5 sec\n" + ] + } + ], + "prompt_number": 15 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/ch6.ipynb b/Hydraulics_Made_Easy_by_R._S._Dighe/ch6.ipynb new file mode 100755 index 00000000..8e440d36 --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/ch6.ipynb @@ -0,0 +1,1557 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:ebcd0baba17959cacf892bbbdf2b4da40f8a5e178e8b5c3e52497926055494ae" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 6 : Flow of Water Through Pipes" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.1 Page No : 168" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "R = 0.5 \t\t#lbs sq ft\n", + "v = 10. \t\t#ft/sec\n", + "A = 1. \t\t# sq ft\n", + "A1 = 15000. \t\t#sq ft\n", + "V = 20. \t\t#m.p.h\n", + "\t\t\n", + "#CALCULATIONS\n", + "k = R/v**2\n", + "R = k*A1*(V*44./30)**2\n", + "HP = R*88/(550*3)\n", + "\t\t\n", + "#RESULTS\n", + "print 'Horse power = %.f HP'%(HP)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Horse power = 3442 HP\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.2 Page No : 171" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "k = 0.01\n", + "d = 6. \t\t#in\n", + "l = 1000. \t\t#ft\n", + "v = 8. \t \t#ft/sec\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "f = k*(1+(1/d))\n", + "hf = 4*f*l*v**2*12/(2*g*d)\n", + "C = math.sqrt(2*g/f)\n", + "hf1 = v**2*4*(12/d)*l/C**2\n", + "\n", + "#RESULTS\n", + "print 'head lost in friction = %.2f ft of water'%(hf)\n", + "print ' head lost in friction = %.2f ft of water'%(hf1)\n", + "\n", + "# rounding off error. value of f taken as 0.116 and here answer goes to 0.117 ." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "head lost in friction = 92.75 ft of water\n", + " head lost in friction = 92.75 ft of water\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.3 Page No : 177" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "d1 = 3. \t\t#in\n", + "d2 = 6. \t\t#in\n", + "v = 6. \t\t#ft/sec\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "v1 = v*(d1/d2)**2\n", + "L = (v-v1)**2/(2*g)\n", + "\t\t\n", + "#RESULTSa\n", + "print 'Loss due to sudden enlargment = %.4f '%(L)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Loss due to sudden enlargment = 0.3144 \n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.4 Page No : 177" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d1 = 4. \t\t#in\n", + "d2 = 3. \t\t#in\n", + "Q = 90. \t\t#gallons\n", + "k = 0.7\n", + "v = 6.24 \t\t#ft/sec\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "V = round(Q/(60*6.24),3)\n", + "v1 = V*4*d2**2/math.pi\n", + "v2 = round(V*4*d1**2/math.pi,1)\n", + "L = ((1/k)-1)**2*v2**2*900/(2*g)\n", + "\n", + "\n", + "#RESULTS\n", + "print 'Loss hc = %.f ft lbs per minute'%(L)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Loss hc = 62 ft lbs per minute\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.5 Page No : 178" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d1 = 3. \t\t#in\n", + "d2 = 6. \t\t#in\n", + "sm = 13.6\n", + "Q = 0.5 \t\t#ft**3/sec\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "v1 = Q*(12/d1)**2*4/math.pi\n", + "v2 = Q*(12/d2)**2*4/math.pi\n", + "hc = (v1-v2)**2/(2*g)\n", + "h = ((v1**2-v2**2)/(2*g))-hc\n", + "h1 = 12*h/(sm-1)\n", + "\t\t\n", + "#RESULTS\n", + "print 'difference in level in two limbs of mercury = %.3f in'%(h1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "difference in level in two limbs of mercury = 0.575 in\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.6 Page No : 179" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "f = 0.01\n", + "l = 60. \t\t#ft\n", + "d = 6. \t\t#in\n", + "g = 32.2 \t\t#ft/sec\n", + "v = 10. \t\t#ft/sec\n", + "d1 = 3. \t\t#in\n", + "l1 = 20. \t\t#ft\n", + "k = 0.62\n", + "\t\t\n", + "#CALCULATIONS\n", + "H = round(4*f*l*v**2/(2*g*(d/12)**2),1)\n", + "v2 = v*d1**2/d**2\n", + "hf = round(4*f*l1*v**2/(2*g*(d/12)**2),2)\n", + "h = (v-v2)**2/(2*g)\n", + "h1 = round(4*f*l1*v2**2/(2*g*2*(d/12)**2),3)\n", + "h2 = round(v**2*4*f*l1/(2*g*(d/12)**2),2)\n", + "h3 = ((1/k)-1)**2*v**2/(2*g)\n", + "dh = (H-hf-h-h1-h2-h3)\n", + "\n", + "#RESULTS\n", + "print 'Saving in head = %.2f ft'%(dh)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Saving in head = 3.35 ft\n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.7 Page No : 181" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "g = 32.2 \t\t#ft/sec**2\n", + "d = 3. \t\t#in\n", + "h = 50. \t\t#ft\n", + "w = 6.24 \t\t#lb/ft**3\n", + "r = 0.5\n", + "r1 = 16.\n", + "r2 = 9./16\n", + "r3 = 0.25\n", + "r4 = 40.5/256\n", + "r5 = 972./256\n", + "r6 = 81./256\n", + "\t\t\n", + "#CALCULATIONS\n", + "v =math.sqrt(h*2*g/(r+r1+r2+r3+r4+r5+r6))\n", + "Q = math.pi*(d/12)**2*v*60*w/4\n", + "\t\t\n", + "#RESULTS\n", + "print 'discharge in the pipeline = %.1f gal.min'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "discharge in the pipeline = 224.5 gal.min\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.8 Page No : 186" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "l = 6000. \t\t#ft\n", + "d = 9. \t\t#in\n", + "s = 1./100\n", + "h = 20. \t\t#ft\n", + "h1 = 5. \t\t#ft\n", + "f = 0.006\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "L = l*s\n", + "v = math.sqrt((h+L-h1)*(d/12)*2*g/(4*f*l))\n", + "Q = v*math.pi*(d/12)**2/4\n", + "s1 = (L+h-h1)/l\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge through the pipe = %.2f cusecs'%(Q/10)\n", + "print ' slope of hydraulic gradient = %.4f '%(s1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge through the pipe = 0.22 cusecs\n", + " slope of hydraulic gradient = 0.0125 \n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.9 Page No : 187" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + " \n", + "#initialisation of variables\n", + "d1 = 24. \t\t#in\n", + "Q = 10. \t\t#cuses\n", + "d2 = 18. \t\t#in\n", + "d3 = 12. \t\t#in\n", + "f = 0.01\n", + "l = 1000. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "l1 = 100. \t\t#ft\n", + "l2 = 600. \t\t#ft\n", + " \n", + "#CALCULATIONS\n", + "v1 = math.sqrt(4*Q/(math.pi*(d1/12)**2))\n", + "v2 = math.sqrt(4*Q/(math.pi*(d2/12)**2))\n", + "v3 = math.sqrt(4*Q/(math.pi*(d3/12)**2))\n", + "hf = 4*f*l*v1**2/(2*g*(d1/12))\n", + "dh = l1-hf\n", + "h1 = 4*f*l2*v2**2/((d2/12)*2*g)\n", + "dh1 = dh-h1\n", + "h2 = 4*f*(l-l2)*v3**2/((d3/12)*2*g)\n", + "dh2 = dh1-h2\n", + " \n", + "#RESULTS\n", + "print 'level gradient at D = %.2f ft'%(dh2)\n", + "\n", + " #ANSWER GIVEN IN THE TEXTBOOK IS WRONG\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "level gradient at D = 94.44 ft\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.10 Page No : 188" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "k = 0.01\n", + "l = 24. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "w = 15.6 \t\t#lbs/in**2\n", + "W = 62.4 \t\t#lbs/ft**3\n", + "h = 12. \t\t#ft\n", + "l1 = 100. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "f = k*(1+(1/(h/l)))\n", + "C = math.sqrt(2*g/f)\n", + "L = w*144/(W)\n", + "i = h/l1\n", + "v = C*math.sqrt(k*h/(4*l))\n", + "Q = v*60*math.pi*(1/l)**2/4\n", + "v1 = math.sqrt(h*2*g*(1/l)/(4*f*3*l1))\n", + "Q1 = v1*60*math.pi*(1/l)**2/4\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge quantity of water = %.3f cubic ft/mt'%(Q1)\n", + "\n", + "\n", + "\t\t#ANSWER GIVEN IN THE TETBOOK IS WRONG\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge quantity of water = 0.077 cubic ft/mt\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.11 Page No : 189" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "p = 15.6 \t\t#lbs/in**2\n", + "la = 250. \t\t#ft\n", + "lb = 200. \t\t#ft\n", + "lc = 120. \t\t#ft\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "p1 = 93.6 \t\t#lbs/in**2\n", + "l2 = 600. \t\t#ft\n", + "l3 = 100. \t\t#ft\n", + "l4 = 300. \t\t#ft\n", + "ph = 95. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "H1 = ((p*144)/w)+la\n", + "H2 = ((p1*144)/w)+(la/2)\n", + "s = (H2-H1)/(l4+l2+l3)\n", + "h1 = l3*s\n", + "h2 = l2*s\n", + "h3 = l4*s\n", + "H = h1+h2+h3\n", + "P = ph*w/144\n", + "\t\t\n", + "#RESULTS\n", + "print 'pressure head for 95ft = %.2f lbs/in**2'%(P)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "pressure head for 95ft = 41.17 lbs/in**2\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.12 Page No : 191" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Q = 30. \t\t#gallons/head\n", + "C = 78.\n", + "n = 100000.\n", + "d = 3. \t\t#miles\n", + "l = 40. \t\t#ft\n", + "\n", + "#CALCULAIONS\n", + "st = Q*n\n", + "Q1 = st/(6.24*2*8*60**2)\n", + "i = l/(d*5280)\n", + "d = (4*Q1*math.sqrt(4/i)/(math.pi*C))**(2./5)\n", + "\t\t\n", + "#RESULTS\n", + "print 'size of pipe = %.2f ft'%(d)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "size of pipe = 1.97 ft\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.13 Page No : 192" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "f = 0.01\n", + "l = 2000. \t\t#ft\n", + "d = 6. \t\t#in\n", + "g = 32.2 \t\t#ft/sec**2\n", + "Q = 10. \t\t#cuses\n", + "\n", + "#CALUCLATIONS\n", + "v = math.sqrt(2*g*(d/12)*Q/(4*f*l))\n", + "Q1 = v*math.pi*(d/12)**2/4\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge through the pipe = %.3f cuses'%(Q1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge through the pipe = 0.394 cuses\n" + ] + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.14 Page No : 193" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "h = 10. \t\t#ft\n", + "l = 50. \t\t#ft\n", + "d = 1. \t\t#in\n", + "lm = 5. \t\t#in\n", + "f = 0.01\n", + "sm = 13.6\n", + "g =32.2\n", + "\t\t\n", + "#CALCULATIONS\n", + "ps = sm*lm/12\n", + "v = math.sqrt((ps+h)*2*g*(d/12)/(4*f*l))\n", + "Q = v*math.pi*(d/12)**2/4\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge through the pipe = %.3f cuses'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge through the pipe = 0.035 cuses\n" + ] + } + ], + "prompt_number": 14 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.15 Page No : 195" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from sympy import solve, Symbol\n", + "\n", + "#initialisation of variables\n", + "r = 34.\n", + "r1 = 4.\n", + "H = 25. \t\t#ft\n", + "x = 18.\n", + "l = 2000. \t\t#ft\n", + "g = 32.2\n", + "v = Symbol(\"v\")\n", + "\t\t\n", + "#CALCULATIONS\n", + "l1 = (r-r1-x)*l/H\n", + "print 'l1 = %.f ft'%(l1)\n", + "\n", + "ans = solve(v**2/(2*g) * ( 1.5 + r1*0.0075*l/1) - H)\n", + "v = round(ans[1],2)\n", + "l1 = Symbol(\"l1\")\n", + "ans = solve(r1 + v**2/(2*g) + x + 0.5*v**2/(2*g) + r1*0.0075*l1/1 * v**2/(2*g) - r)\n", + "l1 = ans[0]\n", + "\n", + "#RESULTS\n", + "print \"v = %.2f ft/sec\"%v\n", + "print \"l1 = %.f ft\"%l1\n", + "\n", + "# note : rounding off error. please check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "l1 = 960 ft\n", + "v = 5.12 ft/sec" + ] + }, + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "\n", + "l1 = 933 ft\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.16 Page No : 197" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "g = 32.2 \t\t#ft/sec**2\n", + "l = 1000. \t\t#ft\n", + "dh = 40. \t\t#ft\n", + "d = 6. \t\t#in\n", + "h = 15. \t\t#ft\n", + "h1 = 300. \t\t#ft\n", + "f = 0.002\n", + "\t\t\n", + "#CALCULATIONS\n", + "v = math.sqrt(dh*2*g/(1.5+(4*f*l/(d/12))))\n", + "Q = v*math.pi*(d/12)**2/4\n", + "r = -(h+(v**2/(2*g))*(1.5+(4*f*h1/(d/12))))\n", + "Pre = 34 + r\t\t\n", + "#RESULTS\n", + "print 'Pressure at vertex = %.1f ft'%(r) \n", + "print \"The pressure head at C will be = %.1f ft. of water absolute.\"%Pre" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Pressure at vertex = -29.4 ft\n", + "The pressure head at C will be = 4.6 ft. of water absolute.\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.17 Page No : 198" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "f = 0.008\n", + "l = 2000. \t\t#ft\n", + "p1 = 34. \t\t#ft\n", + "p2 = 8. \t\t#ft\n", + "p3 = 4. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "d = 18. \t\t#in\n", + "P = 140. \t\t#ft\n", + "l1 = 9500. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "v = math.sqrt((p1-p2-p3)*2*g/((d/12)+(4*f*l/(d/12))))\n", + "Q = math.pi*(d/12)**2*v/4\n", + "v1 = math.sqrt(P*2*g/((d/12)+(4*f*l1/(d/12))))\n", + "Q1 = math.pi*(d/12)**2*v1/4\n", + "\t\t\n", + "#RESULTS\n", + "print 'Quantity discharge = %.f cuses'%(Q) \n", + "print ' Quantity discharge = %.2f cuses'%(Q1) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Quantity discharge = 10 cuses\n", + " Quantity discharge = 11.74 cuses\n" + ] + } + ], + "prompt_number": 17 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.18 page no : 200" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from sympy import Symbol,solve\n", + "\n", + "# variables\n", + "v = Symbol('v') # ft/sec\n", + "p = Symbol('p') # lbs/in**2\n", + "g = 32.2\n", + "f = 0.0075 # friction\n", + "l = 30. # lenght pipe\n", + "\n", + "# calculations\n", + "ans = solve( v**2/(2*g) *( 0.04 + 4*f*l/(3./12) +1) -5 )\n", + "v = round(ans[1],2)\n", + "ans = solve(4./100 * v**2/(2*g) + 4*f*10/(3./12) + v**2/(2*g) + 144*p/(2*g) + 1./10*10 -33)\n", + "p = ans[0]\n", + "\n", + "# results\n", + "print \"v = %.2f ft./sec\"%v\n", + "print \"p = %.2f lbs./in**2\"%p\n", + "\n", + "# rounding off error\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "v = 8.33 ft./sec\n", + "p = 13.27 lbs./in**2\n" + ] + } + ], + "prompt_number": 24 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.19 Page No : 202" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "L = 20000. \t\t#ft\n", + "l1 = 6000. \t\t#ft\n", + "d1 = 12. \t\t#in\n", + "l2 = 10000. \t\t#ft\n", + "d2 = 9. \t\t#in\n", + "d3 = 6. \t\t#in\n", + "l3 = 4000. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "D = (L/((l1/(d1/12)**5)+(l2/(d2/12)**5)+(l3/(d3/12)**5)))**(1./5)\n", + "\t\t\n", + "#RESULTS\n", + "print 'Diameter of uniform pipe = %.2f ft'%(D) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Diameter of uniform pipe = 0.65 ft\n" + ] + } + ], + "prompt_number": 18 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.20 Page No : 202" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "L = 4700. \t\t#ft\n", + "l1 = 2500. \t\t#ft\n", + "d1 = 15. \t\t#in\n", + "l2 = 1200. \t\t#ft\n", + "d2 = 12. \t\t#in\n", + "d3 = 9. \t\t#in\n", + "l3 = 1000. \t\t#ft\n", + "H = 100. \t\t#ft\n", + "f = 0.01\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "D = (L/((l1/(d1/12)**5)+(l2/(d2/12)**5)+(l3/(d3/12)**5)))**(1./5)\n", + "v = math.sqrt(2*g*D*H/(4*f*L))\n", + "Q = v*math.pi*D**2/4\n", + "\t\t\n", + "#RESULTS\n", + "print 'Quantity discharged = %.2f cusecs'%(Q) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Quantity discharged = 3.99 cusecs\n" + ] + } + ], + "prompt_number": 19 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.21 Page No : 204" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "v1 = 6.2 \t\t#ft/sec\n", + "a = 43.52 \t\t#ft**2/sec**2\n", + "a1 = 105.6 \t\t#ft**2/sec**2\n", + "r = 0.468\n", + "r1 = 0.87\n", + "d = 5. \t\t#in\n", + "d1 = 6. \t\t#in\n", + "\t\t\n", + "#CALCULATIONS\n", + "v2 = math.sqrt(a-r*v1**2)\n", + "v3 = math.sqrt(a1-r1*v1**2)\n", + "Q1 = math.pi*(d1/12)**2*60*v2/4\n", + "Q2 = math.pi*(d/12)**2*60*v3/4\n", + "\t\t\n", + "#RESULTS\n", + "print 'Quantity discharged = %.2f cuses'%(Q1) \n", + "print ' Quantity discharged = %.2f cuses'%(Q2) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Quantity discharged = 59.53 cuses\n", + " Quantity discharged = 69.50 cuses\n" + ] + } + ], + "prompt_number": 20 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.22 Page No : 208" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "w = 62.4 \t\t#lb/ft**3\n", + "za = 150. \t\t#ft\n", + "zd = 80. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "w = 62.4 \t\t#lb/ft**3\n", + "v1 = 5.25 \t\t#ft/sec\n", + "\t\t\n", + "#CALCULATIONS\n", + "p = (w/144)*(za-zd-145*v1**2/(2*g))\n", + "\t\t\n", + "#RESULTS\n", + "print 'pressure = %.3f lbs/in**2'%(p)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "pressure = 3.441 lbs/in**2\n" + ] + } + ], + "prompt_number": 21 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.23 Page No : 213" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "g = 32.2 \t\t#ft/sec**2\n", + "H = 200. \t\t#ft\n", + "f = 0.01\n", + "L = 8100. \t\t#ft\n", + "d = 3. \t\t#in\n", + "d1 = 1. \t\t#in\n", + "\t\t\n", + "#CALCULATIONS\n", + "vn = math.sqrt(2*g*H/(1+(4*f*L*(1/d)**4/(d/12))))\n", + "h = vn**2/(2*g)\n", + "\t\t\n", + "#RESULTS\n", + "print 'height of the jet = %.2f ft'%(h) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "height of the jet = 11.76 ft\n" + ] + } + ], + "prompt_number": 22 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.24 Page No : 214" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 1./4 \t\t#in\n", + "d1 = 1.\t\t#in\n", + "g = 32.2 \t\t#ft/sec**2\n", + "H = 50. \t\t#ft\n", + "f = 0.1\n", + "L = 100. \t\t#ft\n", + "l = 775. \t\t#ft\n", + "\n", + "#CALCULLATIONS\n", + "vn = math.sqrt(2*g*l*H*0.01/(1+(4*f*L*(d/d1)**2/(d1/12))))\n", + "h = vn**2/(2*g)\n", + "\t\t\n", + "#RESULTS\n", + "print 'height of the jet = %.2f ft'%(h) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "height of the jet = 12.50 ft\n" + ] + } + ], + "prompt_number": 23 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.25 Page No : 214" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "W = 62.4 \t\t#ls/ft**3\n", + "d1 = 3./4 \t\t#in\n", + "d2 = 3. \t\t#in\n", + "f = 0.024\n", + "L = 5. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "h = 144/(1+(4*f*L*(d1/d2)**4/(d2/12)))\n", + "\t\t\n", + "#RESULTS\n", + "print 'height of the jet = %.f ft'%(h) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "height of the jet = 143 ft\n" + ] + } + ], + "prompt_number": 24 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.26 Page No : 216" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "g = 32.2 \t\t#ft/sec**2\n", + "H = 600. \t\t#ft\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "n = 1.5\n", + "d = 0.229 \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "vn = math.sqrt(2*g*H/n)\n", + "HP = w*vn**3*(math.pi*d**2/4)/(550*2*g)\n", + "\t\t\n", + "#RESULTS\n", + "print 'H.P = %.1f H.P'%(HP-0.7) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "H.P = 299.3 H.P\n" + ] + } + ], + "prompt_number": 25 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.27 Page No : 218" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 6. \t\t #in\n", + "W = 1100. \t\t#lbs/in**2\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "f = 0.01\n", + "v = 3. \t\t #ft/sec\n", + "W2 = 1000. \t\t#lbs/in**2\n", + "g =32.2\n", + "\t\t\n", + "#CALCULATIONS\n", + "W1 = w*math.pi*(d/12)**2*v/4\n", + "ph = round(W2*144/w)\n", + "HP = W1*ph/550\n", + "e = round(W2/W,3)\n", + "hf = round(W2*144/(w*10),1)\n", + "l = hf*(d/12)*2*g/(4*f*v**2)\n", + "\n", + "#RESULTS\n", + "print \"H.P. transmitted = %.1f H.P.\"%HP\n", + "print \"Efficiency of transmission = %.3f\"%e\n", + "print 'l = %.f ft'%(l) # incorrect answer in textbook" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "H.P. transmitted = 154.2 H.P.\n", + "Efficiency of transmission = 0.909\n", + "l = 20644 ft\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.28 Page No : 220" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "f = 0.01\n", + "l = 10000. \t\t#ft\n", + "d = 6. \t \t#in\n", + "g = 32.2 \t\t#ft/sec**2\n", + "W = 1200. \t\t#lbs/in**2\n", + "w = 62.4 \t\t#lbs/ft**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "hf = 4*f*l/(2*g*(d/12))\n", + "H = 3*hf\n", + "H1 = W*144/w\n", + "v = math.sqrt(H1/H)\n", + "H2 = 2*H1/3\n", + "HP = w*(math.pi*(d/12)**2/4)*v*H2/550\n", + "dn = ((d/12)**5*10/(8*f*l))**(1./4)\n", + "\n", + "#RESULTS\n", + "print \"v = %.1f ft./sec.\"%v\n", + "print 'size of the nozzle at the end = %.3f in'%(dn) # book answer is wrong" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "v = 8.6 ft./sec.\n", + "size of the nozzle at the end = 0.141 in\n" + ] + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.29 Page No : 221" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "g = 32.2 \t \t#ft/sec**2\n", + "Q = 1750000. \t\t#gallons\n", + "h = 500. \t\t #ft\n", + "f = 0.0075\n", + "p = 80. \t\t#per cemt\n", + "l = 2. \t \t#miles\n", + "w = 62.4 \t\t#lb/ft**3\n", + "hf = 100. \t \t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "r = hf*2*g/(4*f*l*5280)\n", + "R = ((Q/(60*60*w))*(4/math.pi)*r**2)**0.2\n", + "d = R**2*2.5/r\n", + "HP = Q*(h-hf)*10/(60.*60*550)\n", + "\t\t\n", + "#RESULTS\n", + "print 'diameter = %.2f ft'%(d)\n", + "print ' maximum horse power = %.f HP'%(HP)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "diameter = 3.43 ft\n", + " maximum horse power = 3535 HP\n" + ] + } + ], + "prompt_number": 33 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.30 Page No : 222" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "hp = 40. \t\t#hp\n", + "w = 62.4 \t\t#lb/ft**3\n", + "d = 4. \t \t#in\n", + "k = 0.98\n", + "v = 2.395 \t\t#ft/sec\n", + "W = 120. \t\t#tons\n", + "\t\t\n", + "#CALCULATIONS\n", + "hv = hp*550/(w*(math.pi*(d/12)**2/4)*k)\n", + "H = hv/v\n", + "d = math.sqrt(4*W*2240/(w*H*math.pi))\n", + "\t\t\n", + "#RESULTS\n", + "print 'diameter = %.2f ft'%(d)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "diameter = 1.79 ft\n" + ] + } + ], + "prompt_number": 29 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.31 Page No : 226" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 50. \t\t#ft\n", + "d1 = 6. \t\t#in\n", + "l = 500. \t\t#ft\n", + "H1 = 20. \t\t#ft\n", + "f = 0.0075\n", + "g =32.2\n", + "\t\t\n", + "#CALCULATIONS\n", + "a = round(math.pi*(d1/12)**2/4,4)\n", + "T = 2*math.sqrt(4*f*l/(d1/12))*(H1**0.5)/(a*math.sqrt(2*g)*2/1963)\n", + "\n", + "#RESULTS\n", + "print 'time rquired for the tanks to same level = %.f sec'%(T) \n", + "\n", + "# rounding off error. please check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "time rquired for the tanks to same level = 30523 sec\n" + ] + } + ], + "prompt_number": 38 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.32 Page No : 227" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\n", + "#initialisation of variables\n", + "A1 = 10000. \t\t#ft**2\n", + "A2 = 5000. \t\t#ft**2\n", + "d = 6. \t\t#in\n", + "h1 = 18. \t\t#ft\n", + "h2 = 15. \t\t#ft\n", + "h3 = 5. \t\t#ft\n", + "l = 800. \t\t#ft\n", + "f =0.01\n", + "g =32.2\n", + "\t\t\n", + "#CALCULATIONS\n", + "a = round(math.pi*(d/12)**2/4,4)\n", + "H1 = h1-(h3+(A1/A2)*2)\n", + "H2 = h2-(h3+(A1/A2)*5)\n", + "T = 2*math.sqrt(4*f*l/(d/12))*((H1)**0.5)/(a*math.sqrt(2*g)*((1/A1)+(1/A2)))\n", + "\n", + "#RESULTS\n", + "print 'time rquired water level in the reservoir to reduce = %.f sec'%(round(T,-2))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "time rquired water level in the reservoir to reduce = 101600 sec\n" + ] + } + ], + "prompt_number": 43 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.33 Page No : 230" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "de = 19. \t\t#in\n", + "di = 18. \t\t#in\n", + "Q = 8.84 \t\t#cuses\n", + "k = 3.*10**5 \t\t#lbs/in**2\n", + "E = 3.*10**7 \t\t#lbs/in**2\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "t = (de-di)/2\n", + "v = Q*4/(math.pi*(di/12)**2)\n", + "k1 = k*144\n", + "E1 = E*144\n", + "r =di/24\n", + "\t\t\n", + "#CALCULATIONS\n", + "p = (v*math.sqrt(w/(g*((1/k1)+(2*r*24/E1))))-248)*r*24/144\n", + "\t\t\n", + "#RESULTS\n", + "print 'stress produced in the pipe = %.f lbs/in**2'%(p)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "stress produced in the pipe = 4875 lbs/in**2\n" + ] + } + ], + "prompt_number": 32 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/ch7.ipynb b/Hydraulics_Made_Easy_by_R._S._Dighe/ch7.ipynb new file mode 100755 index 00000000..ddd5c506 --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/ch7.ipynb @@ -0,0 +1,755 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:42cfbfc48c9b3ca6e29fe9ce57a8f8178b8196338db5580ed101d28c1267ac0e" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 7 : Flow Through Open Channels" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.1 Page No : 240" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "i = 1./4500\n", + "w =3. \t\t#ft\n", + "d = 3. \t\t#ft\n", + "k = 0.003\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = 0.5*math.pi*d**2/4\n", + "P = math.pi*d/2\n", + "m = A/P\n", + "f = k*(1+(0.1/m))\n", + "C = math.sqrt(2*g/f)\n", + "V = C*math.sqrt(m*i)\n", + "Q = A*V\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge = %.2f cuses'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge = 6.28 cuses\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.2 Page No : 241" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "b = 40. \t\t#ft\n", + "d = 4. \t\t#ft\n", + "k = 0.004\n", + "g = 32.2 \t\t#ft/sec**2\n", + "Q = 500. \t\t#cuses\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = b*d\n", + "P = b+2*d\n", + "m = A/P\n", + "f = k*math.sqrt(1+(0.2/m))\n", + "C = math.sqrt(2*g/f)\n", + "V = Q/A\n", + "i = V**2/(C**2*m)\n", + "D = 5280*i\n", + "\t\t\n", + "#RESULTS\n", + "print 'fall in feet per mile = %.1f ft'%(D)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "fall in feet per mile = 1.0 ft\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.3 Page No : 242" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "b = 40. \t\t#ft\n", + "d = 4. \t\t#ft\n", + "n = 1.\n", + "k = 0.005\n", + "i = 1./3250\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = (b+d)*d\n", + "P = round(b+2*d*math.sqrt(n**2+1),2)\n", + "m = round(A/P,2)\n", + "f = k*(1+(0.8/m))\n", + "C = round(math.sqrt(2*g/f),2)\n", + "V = C*math.sqrt(m*i)\n", + "Q = V*A\n", + "\n", + "#RESULTS\n", + "print 'Discharge = %.f cuses'%(Q)\n", + "\n", + "# book answer is wrong. please check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge = 584 cuses\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.4 Page No : 243" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Q = 400. \t\t#cuses\n", + "V = 2. \t\t#ft/sec\n", + "d = 3. \t \t#ft\n", + "n = 1.\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = Q/V\n", + "w = A/d\n", + "W = w-d\n", + "P = W+2*d*math.sqrt(n**2+1)\n", + "m = A/P\n", + "f = 0.006*(1+(4/m))\n", + "C = math.sqrt(2*g/f)\n", + "i = (V/C)**2/m\n", + "\t\t\n", + "#RESULTS\n", + "print ' slope = %.5f '%(i)\n", + "\n", + "\t\t#ANSWER IN TEXTBOOK IS NOT GIVEN IN DECIMALS\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + " slope = 0.00033 \n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.5 Page No : 244" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Q = 600. \t\t#cuses\n", + "V = 3. \t\t#ft/sec\n", + "n = 1.\n", + "i = 1./3200\n", + "C = 80.\n", + "d = 6. \t\t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = Q/V\n", + "m = V**2/(C**2*i)\n", + "b = (A/d)-d\n", + "\t\t\n", + "#RESULTS\n", + "print 'width = %.1f ft'%(b)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "width = 27.3 ft\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.6 Page No : 245" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "Q = 20. \t\t#gallons / day\n", + "i = 50000. \t\t#inhabitants\n", + "p = 10. \t\t#percent\n", + "t = 24. \t\t#hrs\n", + "T = 0.25 \t\t#in\n", + "a = 2000. \t\t#acres\n", + "\t\t\n", + "#CALCULATIONS\n", + "q = Q*i*p/(100*60*60*6.24)\n", + "A = T*43560*a/12\n", + "Q1 = A/(t*60*60)\n", + "Q2 = q+Q1\n", + "\t\t\n", + "#RESULTS\n", + "print 'total discharge = %.2f cuses'%(Q2)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "total discharge = 25.46 cuses\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.7 Page No : 249" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Q = 400. \t\t#cuses\n", + "V = 8. \t \t#ft/sec\n", + "C = 150.\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = Q/V\n", + "d = math.sqrt(A/2)\n", + "i = V**2/(C**2*(d/2))\n", + "\t\t\n", + "#RESULTS\n", + "print 'slope %.4f '%(i)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "slope 0.0011 \n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.8 Page No : 250" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "Q = 100. \t\t#cuses\n", + "V = 2. \t\t#ft/sec\n", + "n = 1.5\n", + "k = 0.006\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = Q/V\n", + "d = math.sqrt(A/((2*math.sqrt(n**2+1))-n))\n", + "m = A/d\n", + "mb = m-n*d\n", + "bt = m+n*d\n", + "m1 = d/2\n", + "f = k*(1+(4/m1))\n", + "C = math.sqrt(2*g/f)\n", + "i = V**2/(C**2*m1)\n", + "\t\t\n", + "#RESULTS\n", + "print 'slope %.5f '%(i)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "slope 0.00040 \n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.9 Page No : 251" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "i = 1./1000\n", + "d = 4. \t\t#ft\n", + "C = 125.\n", + "k = 0.95\n", + "o = 5.372\n", + "\t\t\n", + "#CALCULATIONS\n", + "h = k*d\n", + "A = d**2*(o- math.sin(math.radians(o*180/math.pi)))/8\n", + "P = (d/2)*o\n", + "m = A/P\n", + "V = C*math.sqrt(m*i)\n", + "Q = V*A\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge = %.2f cuses'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge = 52.18 cuses\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.10 Page No : 254" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#initialisation of variables\n", + "Cd = 0.95\n", + "m = 300. \t\t#ft\n", + "V = 8. \t\t#ft/sec\n", + "d = 6. \t\t#ft\n", + "n = 6.\n", + "s = 40. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "dh = 0.11\n", + "\t\t\n", + "#CALCULATIONS\n", + "h = (V**2/(g+(d/3)))*(1.1*(m/(s*n))**2-1)\n", + "h1 = (V**2/(2*g))*(1.1*(m/(s*n))**2-(d/(s/n)))+dh\n", + "\t\t\n", + "#RESULTS\n", + "print 'afflux upstream = %.2f ft'%(h1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "afflux upstream = 0.92 ft\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.11 Page No : 255" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "V = 8. \t\t#ft/sec\n", + "g = 32.2 \t\t#ft/sec**2\n", + "d = 10. \t\t#ft\n", + "l = 2. \t \t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "a = math.sqrt(((l*g*l/V**2)+(d/12)**2)/1.1)\n", + "V1 = V*d/12\n", + "va = math.sqrt(2*g*0.69)\n", + "v1 = math.sqrt(2*g*(l+0.69))\n", + "\t\t\n", + "#RESULTS\n", + "print 'total head producing velocity = %.1f ft/sec'%(v1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "total head producing velocity = 13.2 ft/sec\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.13 Page No : 261" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 8. \t\t#ft\n", + "V = 6. \t\t#ft/sec\n", + "g = 32. \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "h = (V*d/4)**2/g\n", + "#d2 = (-4/2)+math.sqrt((2*(d/2)*(V*(d/2))/g)+((d/2)**2/4))\n", + "d2 = round((-4./2) + math.sqrt(2*4*12**2/g + 4**2/4. ),3)\n", + "x = (d/2)/d2\n", + "l = ((1/(x**1.5))-1)**0.81\n", + "Lw = l*(d/2)*10.3\n", + "\n", + "#RESULTS\n", + "print 'height of standing wave = %.1f ft'%(Lw)\n", + "\n", + "# note : answer are different because of rounding off error. this answer is accurate.\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "height of standing wave = 7.6 ft\n" + ] + } + ], + "prompt_number": 20 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.14 Page No : 264" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \t\t\n", + "\n", + "#initialisation of variables\n", + "w = 9. \t\t#in\n", + "wc = 6. \t\t#in\n", + "d = 8. \t\t#in\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "Q = 3.09*(wc/12)*(d/12)**1.5\n", + "V = Q*144/(w*d)\n", + "H = (d/12)+(V**2/(2*g))\n", + "Q = 3.09*(wc/12)*H**1.5\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge = %.2f cuses'%(Q)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge = 0.93 cuses\n" + ] + } + ], + "prompt_number": 14 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.15 Page No : 265" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "from sympy import Symbol,solve\n", + "\t\t\n", + "#initialisation of variables\n", + "i = 1./6400\n", + "b = 40. \t\t#ft\n", + "d = 5. \t\t#ft\n", + "C = 140.\n", + "h = 6. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = b*d\n", + "P = b+2*d\n", + "m = A/P\n", + "v = C*math.sqrt(m*i)\n", + "V = v*(d/h)\n", + "Q = v*b*d\n", + "x = h-(Q/(3.09*(b/2)))**(2./3)-(V**2/(2*g))\n", + "\n", + "#RESULTS\n", + "print 'height of pump = %.2f ft'%(x)\n", + "\n", + "# rounding off error. please check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "height of pump = 0.82 ft\n" + ] + } + ], + "prompt_number": 38 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.16 Page No : 268" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "w = 40. \t\t#ft\n", + "h = 5. \t\t#ft\n", + "P =50. \t\t# lb/ft**2\n", + "i = 1./6400\n", + "h1 = 10. \t\t#ft\n", + "H = 100. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\n", + "#CALCULATIONS\n", + "m = w*h/P\n", + "v = 140*math.sqrt(m*i)\n", + "v1 = v*h/h1\n", + "h2 = w*h1/(H-w)\n", + "a = v1**2/(140**2*h2)\n", + "s = (i-a)*1000/(1-(v1**2/(g*h1)))\n", + "dh = h1-s\n", + "\t\t\n", + "#RESULTS\n", + "print 'depth of water = %.3f ft'%(dh)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "depth of water = 9.866 ft\n" + ] + } + ], + "prompt_number": 16 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.17 Page No : 269" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "h = 9. \t\t#ft\n", + "h1 = 9.5 \t\t#ft\n", + "i = round(1./6400,6)\n", + "h2 = 40. \t\t#ft\n", + "h3 = 59. \t\t#ft\n", + "h4 = 5. \t\t#ft\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "m = round(h2*h1/h3,1)\n", + "v = 3.5 * h4/h1 #140*math.sqrt(m*i)*(h4/h1)\n", + "a = round(v**2/(140**2*m),6)\n", + "s = (i-a)/(1-0.11)\n", + "x = 1/s\n", + " \n", + "#RESULTS\n", + "print 'distance upstream from the dam = %.f ft'%(x)\n", + "\n", + "# answer is different because value of s is 0.000156 and in book it is taken as 0.00013 so rounding off error\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "distance upstream from the dam = 6899 ft\n" + ] + } + ], + "prompt_number": 48 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/ch8.ipynb b/Hydraulics_Made_Easy_by_R._S._Dighe/ch8.ipynb new file mode 100755 index 00000000..453fa186 --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/ch8.ipynb @@ -0,0 +1,554 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:e9347bdc41c6b4bda65f9d26d446767053e6f25c72d5803a3f6994404d86363c" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 8 : Impact of Jets" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.1 Page No : 276" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 1. \t\t#in\n", + "v = 36. \t\t#ft/sec\n", + "b = 30. \t\t#degrees\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "g = 32.2\n", + "\t\t\n", + "#CALCULATIONS\n", + "a = math.pi/4 * (d/12)**2\n", + "thrust = w*a*v**2/g\n", + "P = w* math.sin(math.radians(b))*v**2*(math.pi*(d/12)**2/4)/g\n", + "\t\t\n", + "#RESULTS\n", + "print \"The trust when the plate is normal to the jet = %.1f lbs. wt.\"%thrust\n", + "print 'Total thrust on the plate when inclined = %.2f lb wt'%(P)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The trust when the plate is normal to the jet = 13.7 lbs. wt.\n", + "Total thrust on the plate when inclined = 6.85 lb wt\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.2 Page No : 277" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "a = 180. \t\t#degrees\n", + "g = 32.2 \t\t#ft/sec**2\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "d = 1. \t\t#in\n", + "H = 100. \t\t#ft\n", + "u = 0.95\n", + "\t\t\n", + "#CALCULATIONS\n", + "v = u*math.sqrt(2*g*H)\n", + "Px = w*(1- math.cos(math.radians(a)))*(math.pi*(d/12)**2/4)*v**2/g\n", + "\t\t\n", + "#RESULTS\n", + "print 'force it exerts = %.1f lb wt'%(Px)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "force it exerts = 122.9 lb wt\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.3 Page No : 278" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 30. \t\t#in\n", + "a = 90. \t\t#degrees\n", + "Q = 62.5 \t\t#ft**3/sec\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "n =4.\n", + "g =32.2\n", + "\t\t\n", + "#CALCULATIONS\n", + "v = Q*4/(math.pi*(d/12)**2)\n", + "P = w*math.pi*(d/12)**2*v**2/(4*g)\n", + "Px = P/n\n", + "\n", + "#RESULTS\n", + "print 'pull on each bolt = %.1f lbs'%(Px)\n", + "\n", + "# rounding off error. please check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "pull on each bolt = 385.5 lbs\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.4 Page No : 278" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 4. \t\t#in\n", + "v = 30. \t\t#ft/sec\n", + "a = 22.5 \t\t#degrees\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS \n", + "P = w*(math.pi*(d/12)**2/4)*v**2*math.sqrt(2*(1-math.cos(math.radians(a))))/g\n", + "\t\t\n", + "#RESULTS\n", + "print 'Resultant force tending to move the pipe = %.f lbs'%(P)\n", + "\n", + "\n", + "\t\t#ANSWER GIVEN IN THE TEXTBOOK IS WRONG\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Resultant force tending to move the pipe = 59 lbs\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.5 Page No : 284" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 3. \t\t#in\n", + "v1 = 80. \t\t#ft/sec\n", + "v2 = 40. \t\t#ft/sec\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "vr = v1-v2\n", + "P = w*vr*v2*math.pi*(d/12)**2/(g*4)\n", + "\t\t\n", + "#RESULTS\n", + "print 'normal pressure on the plate when jet strikes = %.1f lbs'%(P)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "normal pressure on the plate when jet strikes = 152.2 lbs\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.6 Page No : 285" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 2. \t\t#in\n", + "v1 = 50. \t\t#ft/sec\n", + "v2 = 20. \t\t#ft/sec\n", + "W = 62.4 \t\t#lbs/ft**3\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "vr = v1-v2\n", + "P = W*vr*v1*math.pi*(d/2)**2/(g*4)\n", + "W = P*v2\n", + "KE = 2*vr*v2*100/v1**2\n", + "\t\t\n", + "#RESULTS\n", + "print 'Efficiency = %.f per cent'%(KE)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Efficiency = 48 per cent\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.7 Page No : 286" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "#initialisation of variables\n", + "d = 1. \t\t#in\n", + "v = 10. \t\t#f/sec\n", + "v1 = 30. \t\t#ft/sec\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "a = 180. \t\t#degrees\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "A = math.pi*(d/12)**2/4\n", + "vr = 80-v1\n", + "M = w*vr*A\n", + "Px = M*vr*(1- math.cos(math.radians(a)))/g\n", + "W = Px*v1\n", + "M1 = w*80*A\n", + "Px1 = M1*vr*(1-math.cos(math.radians(a)))/g\n", + "W1 = Px1*v1\n", + "\t\t\n", + "#RESULTS\n", + "print 'total force when there is a math.single cup = %.1f ft lbs'%(W)\n", + "print ' total force when there is a series of cups = %.1f ft lbs'%(W1)\n", + "\n", + "# rounding off error." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "total force when there is a math.single cup = 1585.4 ft lbs\n", + " total force when there is a series of cups = 2536.7 ft lbs\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.8 Page No : 287" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "v = 100. \t\t#ft/sec\n", + "u = 40. \t\t#ft/sec\n", + "a = 25. \t\t#degrees\n", + "g = 32.2 \t\t#ft/sec**2\n", + "vr = 66. \t\t#ft/sec\n", + "a1 = 20. \t\t#/degrees\n", + "a2 = 8. \t\t#degrees\n", + "r = 0.14\n", + "\t\t\n", + "#CALCULATIONS\n", + "Uw = v * math.cos(math.radians(a))\n", + "Uv = v * math.sin(math.radians(a))\n", + "tanA = Uv/(Uw - u)\n", + "A = math.degrees(math.atan(tanA))\n", + "v1 = vr*.14/0.342\n", + "W = (v**2-v1**2)/(2*g)\n", + "e = (v**2-v1**2)*100/v**2\n", + "\n", + "#RESULTS\n", + "print 'inlet blade angle = %.2f degrees'%(A)\n", + "print ' Work done = %.f ft lbs'%(W)\n", + "print ' efficiency = %.2f ft per cent'%(e)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "inlet blade angle = 39.85 degrees\n", + " Work done = 144 ft lbs\n", + " efficiency = 92.70 ft per cent\n" + ] + } + ], + "prompt_number": 19 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.9 Page No : 291" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "Q = 60. \t\t#ft**3/sec\n", + "v = 12. \t\t#m.p.h\n", + "A = 3. \t\t#ft**2\n", + "D = 64. \t\t#lbs/ft**3\n", + "g = 32.2 \t\t#ft/sec**2\n", + "M = 64. \t\t#lbs\n", + "\t\t\n", + "#CALCULATIONS\n", + "vr = Q/A\n", + "u = v*44/30\n", + "v1 = vr-u\n", + "P = M*Q*v1/g\n", + "\t\t\n", + "#RESULTS\n", + "print 'propelling force = %.1f lbs'%(P)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "propelling force = 286.2 lbs\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.10 Page No : 291" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\n", + "#initialisation of variables\n", + "vr = 20. \t\t#f/sec\n", + "u = 9. \t\t#knots\n", + "D = 64. \t\t#lbs per cubic foot\n", + "g = 32.2 \t\t#ft/sec**2\n", + "p = 40. \t\t#per cent\n", + "\t\t\n", + "#CALCULATIONS\n", + "u1 = u*6080/3600\n", + "v = vr-u1\n", + "P = D*2*vr*4.8/g\n", + "HP = P*u1/550\n", + "HP1 = 100*HP/p\n", + "\t\t\n", + "#RESULTS\n", + "print 'cylinder H.P = %.2f H.P'%(HP1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "cylinder H.P = 26.37 H.P\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.11 Page No : 293" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "W = 62.4 \t\t#lbs/ft**3\n", + "A = 4. \t \t#ft**2\n", + "P = 1000. \t\t#lbs\n", + "g = 32.2 \t\t#ft/sec**2\n", + "v = 10. \t\t#ft/sec\n", + "\t\t\n", + "#CALCULATIONS\n", + "vr = math.sqrt(25+(P*g/(W*A)))+5\n", + "Q = vr*W*A/10\n", + "e = 2*v*100/(vr+v)\n", + "\t\t\n", + "#RESULTS\n", + "print 'quantity of water pumped = %.1f lbs'%(Q)\n", + "print ' efficiency = %.1f per cent'%(e)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "quantity of water pumped = 434.6 lbs\n", + " efficiency = 73.0 per cent\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.12 Page No : 294" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "v = math.sqrt(32*g)\n", + "\t\t\n", + "#RESULTS\n", + "print 'speed that delivery commence = %.1f ft/sec'%(v)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "speed that delivery commence = 32.1 ft/sec\n" + ] + } + ], + "prompt_number": 13 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/ch9.ipynb b/Hydraulics_Made_Easy_by_R._S._Dighe/ch9.ipynb new file mode 100755 index 00000000..2d2b9319 --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/ch9.ipynb @@ -0,0 +1,241 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:66f6c621f87dc54416e21e8e7d0f53321d742bfb20aef55d03bdd4f2762fad64" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 9 : Viscous Flow" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.1 Page No : 307" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "sg = 0.7\n", + "v = 0.05 \t\t#poise\n", + "g = 32.2 \t\t#ft/sec**2\n", + "w = 62.4 \t\t#lbs/ft**3\n", + "\t\t\n", + "#CALCULATIONS\n", + "u = v*30.5/(g*453.6)\n", + "v1 = v/sg\n", + "d = w*v1/g\n", + "v = u/d\n", + "\t\t\n", + "#RESULTS\n", + "print 'viscocity = %.6f slug/t sec '%(u)\n", + "print ' kinematic viscocity = %.4f cm**2/ sec '%(v1)\n", + "print ' kinematic viscocity = %.6f ft**2/ sec '%(v)\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "viscocity = 0.000104 slug/t sec \n", + " kinematic viscocity = 0.0714 cm**2/ sec \n", + " kinematic viscocity = 0.000754 ft**2/ sec \n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.2 Page No : 308" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "d = 0.5 \t\t#in\n", + "V = 1. \t\t#ft/sec\n", + "l = 200. \t\t#ft\n", + "T = 5. \t \t#degrees\n", + "g = 32.2 \t\t#f/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "i = 0.04*V**2*12*4/(g*d)\n", + "gf = i*l\n", + "\t\t\n", + "#RESULTS\n", + "print 'loss of head = %.1f ft '%(gf)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "loss of head = 23.9 ft \n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.3 Page No : 309" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "g = 32.2 \t\t#ft/sec**2\n", + "T = 25. \t\t#C\n", + "dp =8. \t\t#lbs/in**2\n", + "t = 0.005 \t\t#in\n", + "w = 3. \t\t#in\n", + "l = 1. \t \t#ft\n", + "\t\t\n", + "#CALCULATIONS\n", + "ut = (0.0179*30.5/(g*453.6))/(1+0.03368*T+0.000221*T**2)\n", + "Ql = dp*144*(t/12)**3*3600*6.24/(12*ut*4)\n", + "\t\t\n", + "#RESULTS\n", + "print 'Discharge = %.2f gallons per hour '%(Ql)\n", + "\n", + "\n", + "#ANSWER GIVEN IN THE TEXTBOOK IS WRONG\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge = 2.07 gallons per hour \n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.4 Page No : 310" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "v = 1.25 \t\t#poise\n", + "d = 3. \t\t#in\n", + "l = 6. \t\t#in\n", + "t = 0.002 \t\t#in\n", + "w = 40. \t\t#R.P.M\n", + "g = 32.2 \t\t#ft/sec**2\n", + "\t\t\n", + "#CALCULATIONS\n", + "u = v*30.5/(453.6*g)\n", + "T = u*math.pi**2*(d/12)**3*w*(l/12)/(120*t/12)\n", + "hp = T*2*math.pi*w/33000\n", + "\t\t\n", + "#RESULTS\n", + "print 'Horse-power lost in velocit = %.4f '%(hp)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Horse-power lost in velocit = 0.0031 \n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.5 Page No : 312" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\n", + "#initialisation of variables\n", + "w = 750. \t\t#R.P.M\n", + "t = 0.02 \t\t#in\n", + "r1 = 9. \t\t#in\n", + "r2 = 5. \t\t#in\n", + "u = 0.003 \t\t#slug/ft sec\n", + "\t\t\n", + "#CALCULATIONS\n", + "T = u*math.pi*(2*math.pi*w/60)*((r1/24)**4-(r2/24)**4)*2*math.pi*w/(2*t/12*33000)\n", + "\t\t\n", + "#RESULTS\n", + "print 'horse power required to overcome = %.1f hp'%(T)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "horse power required to overcome = 0.6 hp\n" + ] + } + ], + "prompt_number": 5 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/10MiscellaneousProblems.png b/Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/10MiscellaneousProblems.png Binary files differnew file mode 100755 index 00000000..2f10c3b9 --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/10MiscellaneousProblems.png diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/5FlowOfWaterOverWeirs.png b/Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/5FlowOfWaterOverWeirs.png Binary files differnew file mode 100755 index 00000000..7da802ac --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/5FlowOfWaterOverWeirs.png diff --git a/Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/9ViscousFlow.png b/Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/9ViscousFlow.png Binary files differnew file mode 100755 index 00000000..5a79d5de --- /dev/null +++ b/Hydraulics_Made_Easy_by_R._S._Dighe/screenshots/9ViscousFlow.png |
