{
 "metadata": {
  "name": "",
  "signature": "sha256:d25e06b98904a679939aa950e605d9ba5c3f70e4cec7b065f7a8a4ce885f3659"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 1 : Fluid pressure and Its Measurement"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.1 Page No : 9"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "p = 343350. \t\t\t#pressure at any point in pa\n",
      "w = 9810. \t\t\t    #gravitational consmath.tant\n",
      "s1 = 1.\n",
      "s2 = 13.6\n",
      "\n",
      "# Calculations \n",
      "h1 = (p/w)\n",
      "h2 = h1*s1/s2\n",
      "\n",
      "# Results \n",
      "print \"pressure in term of height of water(m)\",h1\n",
      "print \"pressure in term of height of mercury(m)\",round(h2,7)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "pressure in term of height of water(m) 35.0\n",
        "pressure in term of height of mercury(m) 2.5735294\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.2 Page No : 10"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 0.75 \t\t\t#atm pressure in term of mercury\n",
      "w = 9810.\n",
      "w1 = 13.6*w \t\t\t#specific weight of mercury\n",
      "w2 = 15000.\n",
      "h2 = 3. \t\t\n",
      "\n",
      "# Calculations \n",
      "Patm = w1*h1\n",
      "p = w2*h2 \t\t\t# gauge pressure\n",
      "Pabs = Patm+p\n",
      "\n",
      "# Results \n",
      "print \"gauge pressure(N/m2)\",p\n",
      "print \"absolute pressure(N/m2)\",Pabs\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "gauge pressure(N/m2) 45000.0\n",
        "absolute pressure(N/m2) 145062.0\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.3 Page No : 10"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 2.5\n",
      "h2 = 1.5\n",
      "s1 = 1\n",
      "s2 = 0.8\n",
      "w = 9810\n",
      "\n",
      "# Calculations \n",
      "p2 = s2*w*h2 \t\t\t#Pressure intensity at interface\n",
      "p1 = s1*w*h1\n",
      "p = p1+p2\n",
      "\n",
      "# Results \n",
      "print  \"pressure intensity at interface(N/m2)\",p2\n",
      "print  \"pressure intensity at bottom(N/m2)\",p\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "pressure intensity at interface(N/m2) 11772.0\n",
        "pressure intensity at bottom(N/m2) 36297.0\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.4 Page No : 11"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "p = 71613. \t\t\t#gauge pressure\n",
      "w = 9810\n",
      "\n",
      "# Calculations \n",
      "phead = p/w\n",
      "patm = 10.33\n",
      "pabs = patm+phead\n",
      "\n",
      "# Results \n",
      "print \"absolute pressure in term of water height in meters\",pabs\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "absolute pressure in term of water height in meters 17.63\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.5 Page No : 11"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 0.05\n",
      "h2 = 0.1\n",
      "s1 = 0.8\n",
      "s2 = 13.6\n",
      "w = 9810\n",
      "\n",
      "# Calculations \n",
      "p = s2*h2*w \t\t\t#pressure at balance line \n",
      "p1 = s1*h1*w \n",
      "pf = p-p1\n",
      "\n",
      "# Results \n",
      "print  \"pressure in pipe(N/m2)\",pf\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "pressure in pipe(N/m2) 12949.2\n"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.6 Page No : 12"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 0.2\n",
      "h2 = 0.5\n",
      "s1 = 0.9\n",
      "s2 = 13.6\n",
      "\n",
      "# Calculations \n",
      "h = -(h1*s1+h2*s2)\n",
      "w = 9810\n",
      "p = h*w\n",
      "\n",
      "# Results \n",
      "print  \"vacuum pressure (N/m2)\",p\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "vacuum pressure (N/m2) -68473.8\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.7 Page No : 13"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "s1 = 0.8\n",
      "s2 = 13.6\n",
      "dh = 0.4\n",
      "\n",
      "# Calculations \n",
      "h = dh*13.6-dh*0.8\n",
      "w = 9810\n",
      "pd = w*h\n",
      "\n",
      "# Results \n",
      "print  \"pressure difference in height of water()\",h\n",
      "print  \"presuure difference in N/m2\",pd\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "pressure difference in height of water() 5.12\n",
        "presuure difference in N/m2 50227.2\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.8 Page No : 14"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "s1 = 0.8\n",
      "s2 = 0.7\n",
      "h1 = 1.5\n",
      "h2 = 0.3\n",
      "h3 = 0.7\n",
      "s3 = 13.6\n",
      "\n",
      "# Calculations \n",
      "hd = h2*s2+h3*s3-h1*s1\n",
      "w = 9810\n",
      "pd = hd*w\n",
      "\n",
      "# Results \n",
      "print  \"diffrence in pressure in term of height of water(m)\",hd\n",
      "print  \"difference in pressure (N/m2)\",pd\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "diffrence in pressure in term of height of water(m) 8.53\n",
        "difference in pressure (N/m2) 83679.3\n"
       ]
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.9 Page No : 15"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "s1 = 1.6\n",
      "s2 = 0.8\n",
      "s3 = 13.6\n",
      "p1 = 98100.\n",
      "p2 = 176580.\n",
      "w = 9810.\n",
      "\n",
      "# Calculations \n",
      "h1 = p1/w\n",
      "h2 = p2/w\n",
      "h = (h2-h1+1.6*s2-4.1*s1)/(s3-s2)\n",
      "\n",
      "# Results \n",
      "print \"difference in mercury level(cm)\",h*100\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "difference in mercury level(cm) 21.25\n"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.10 Page No : 16"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "s1 = 1.2\n",
      "s2 = 1\n",
      "s3 = 0.7\n",
      "\n",
      "# Calculations \n",
      "h = (s1-s2)*0.3/(s2-s3)\n",
      "\n",
      "# Results \n",
      "print \"difference in height(cm)\",h*100\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "difference in height(cm) 20.0\n"
       ]
      }
     ],
     "prompt_number": 10
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.11 Page No : 17"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "s1 = 0.8\n",
      "s2 = 13.6\n",
      "z = 0.02\n",
      "w = 9810\n",
      "h2 = 0.2\n",
      "h1 = 0.1\n",
      "\n",
      "# Calculations \n",
      "h = h2*s2-h1*s1+(z*h2*(s2-s1))\n",
      "p = h*w\n",
      "\n",
      "# Results \n",
      "print \"pressure of the oil in N/m2\",p\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "pressure of the oil in N/m2 26400.672\n"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.12 Page No : 18"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "l = 4\n",
      "b = 2\n",
      "h = 3\n",
      "w = 9810\n",
      "s = 0.8\n",
      "\n",
      "# Calculations \n",
      "p1 = w*l*b*h*s\n",
      "p2 = w*s*l*h*1.5\n",
      "p3 = w*s*b*h*1.5\n",
      "\n",
      "# Results \n",
      "print \"total pressure on horizontal base\",p1\n",
      "print  \"total pressure on larger vertical base\",p2\n",
      "print \"total pressure on smaller vertical walls\",p3\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "total pressure on horizontal base 188352.0\n",
        "total pressure on larger vertical base 141264.0\n",
        "total pressure on smaller vertical walls 70632.0\n"
       ]
      }
     ],
     "prompt_number": 12
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.13 Page No : 18"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "# Variables\n",
      "p = 490500.\n",
      "w = 9810.\n",
      "\n",
      "# Calculations \n",
      "h = p/w\n",
      "D = 0.15\n",
      "A = math.pi*D*D*0.25\n",
      "pt = w*A*h\n",
      "h1 = (D*D)/(16*h)\n",
      "\n",
      "# Results \n",
      "print \"total hydrostatic pressure in N\",round(pt,2)\n",
      "print \"position of centre of pressure below the centre of pipe : %f\"%h1\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "total hydrostatic pressure in N 8667.85\n",
        "position of centre of pressure below the centre of pipe : 0.000028\n"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.14 Page No : 19"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "# Variables\n",
      "w = 9810.\n",
      "h = 4.\n",
      "d = 2.\n",
      "\n",
      "# Calculations \n",
      "a = d*d*0.25*math.pi\n",
      "p = w*a*h\n",
      "h1 = d*d/64\n",
      "T = p*(h1) \n",
      "\n",
      "# Results \n",
      "print  \"torque required to keep the disc in vertical position in N.m %.4f\"%T\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "torque required to keep the disc in vertical position in N.m 7704.7560\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.15 Page No : 20"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "w = 9810.\n",
      "h = 2.\n",
      "l = 2.\n",
      "b = 1.\n",
      "\n",
      "# Calculations \n",
      "a = l*b\n",
      "p = w*a*h\n",
      "h1 = h+(b*l*l*l/(12*b*l*l))\n",
      "\n",
      "# Results \n",
      "print \"total pressure\",p\n",
      "print \"location of its centre of pressure\",round(h1,4)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "total pressure 39240.0\n",
        "location of its centre of pressure 2.1667\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.16 Page No : 21"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "h1 = 8\n",
      "w = 9810\n",
      "wd = 6\n",
      "\n",
      "# Calculations \n",
      "p1 = 0.5*w*h1*h1*wd\n",
      "h2 = 4\n",
      "p2 = 0.5*h2*w*h2*wd\n",
      "h11 = 0.66666*h1\n",
      "h22 = 0.6666*h2\n",
      "p = p1-p2\n",
      "hf = (p1*(h1-h11)-p2*(h2-h22))/p\n",
      "\n",
      "# Results \n",
      "print \"resultant force\",p\n",
      "print \"position of its line of action\",round(hf,5)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "resultant force 1412640.0\n",
        "position of its line of action 3.11109\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.17 Page No : 22"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "z = 9810\n",
      "w = 10\n",
      "h = 2\n",
      "\n",
      "# Calculations \n",
      "p = 0.5*h*h*w*z\n",
      "h1 = h*0.6666\n",
      "\n",
      "# Results \n",
      "print \"total hydrostatic thrust\",p\n",
      "print \"its point of application\",h1\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "total hydrostatic thrust 196200.0\n",
        "its point of application 1.3332\n"
       ]
      }
     ],
     "prompt_number": 17
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.18 Page No : 22"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "a1 = 1.4*2.2*1.4\n",
      "x1 = 1.6+0.7\n",
      "x11 = (1.4*1.4/(12*2.3))+x1\n",
      "x2 = 0.7\n",
      "x22 = (1.4*1.4/(12*0.7))+x2\n",
      "z = 9810\n",
      "\n",
      "# Calculations \n",
      "p1 = z*a1*x1\n",
      "p2 = z*a1*x2/1.4\n",
      "p = p1-p2\n",
      "h = (p1*(3-x11)+p2*(1.4-x22))/p\n",
      "f = (p1*(3-x11)-p2*(1.4-x22))/1.4\n",
      "\n",
      "# Results \n",
      "print \"resultant force\",p\n",
      "print \"force acting horizontally on the top of the gate\",f\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "resultant force 76141.296\n",
        "force acting horizontally on the top of the gate 36660.624\n"
       ]
      }
     ],
     "prompt_number": 18
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.19 Page No : 24"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "s = 1.5\n",
      "s1 = 0.9\n",
      "w = 9810\n",
      "h1 = 0.9\n",
      "h2 = 0.6\n",
      "\n",
      "# Calculations \n",
      "p1 = 0.5*w*s*s1*h1*h1 \t\t\t#total pressure due to oil\n",
      "p2 = w*h1*h2*s*s1 \t\t\t# total pressure due to oil above water\n",
      "p3 = w*h2*h2*0.5*s \t\t\t#total pressure due to water\n",
      "p = p1+p2+p3\n",
      "h = ((p1*0.6666*h1)+(p2*(h1+0.5*h2))+(p3*(0.6666*h2+h1)))/p\n",
      "\n",
      "# Results \n",
      "print \"resultant pressure on the wall in N/m2\",p\n",
      "print \"position of centre of pressure from free surface\",round(h,5)\n",
      "\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "resultant pressure on the wall in N/m2 15163.8075\n",
        "position of centre of pressure from free surface 1.00521\n"
       ]
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.20 Page No : 26"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\n",
      "# Variables\n",
      "d = 2.4\n",
      "h = 1.6\n",
      "s = 1.2\n",
      "\n",
      "# Calculations \n",
      "a = d*s\n",
      "w = 9810\n",
      "p = w*a*h*s\n",
      "h1 = ((2*s*s*s*d)/(12*a*h))+h\n",
      "\n",
      "# Results \n",
      "print \"total pressure in N\",p\n",
      "print \"its point of application\",h1\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "total pressure in N 54245.376\n",
        "its point of application 1.75\n"
       ]
      }
     ],
     "prompt_number": 20
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.21 Page No : 26"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "# Variables\n",
      "x = 30.\n",
      "d = 1.2\n",
      "h = 1.5\n",
      "w = 9810.\n",
      "\n",
      "# Calculations \n",
      "z = math.sin(x*math.pi/180)\n",
      "h1 = (z*d*0.5)+h\n",
      "a = 0.25*math.pi*d*d\n",
      "p = a*w*h1\n",
      "h11 = (d*d*z*z)/(16*h1)+h1\n",
      "\n",
      "# Results \n",
      "print \"total pressure\",round(p,3)\n",
      "print \"position of centre of pressure\",round(h11,4)\n",
      "\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "total pressure 19970.728\n",
        "position of centre of pressure 1.8125\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.22 Page No : 27"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "# Variables\n",
      "d = 2\n",
      "z = 0.5\n",
      "\n",
      "# Calculations \n",
      "h = z+0.5*d\n",
      "w = 9810\n",
      "a = math.pi*d*d*0.25\n",
      "p = a*w*h\n",
      "h11 = (1/(16*1.5))+1.5\n",
      "\n",
      "# Results \n",
      "print \"total pressure on the plate %.3f N\"%p\n",
      "print \"position of centre of pressure\",round(h11,5)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "total pressure on the plate 46228.536 N\n",
        "position of centre of pressure 1.54167\n"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.23 Page No : 28"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "# Variables\n",
      "x = 30.\n",
      "z = math.sin(x*math.pi/180)\n",
      "h = 6-(z*0.5)\n",
      "l = 1.\n",
      "b = 4.\n",
      "a = l*b\n",
      "w = 9810.\n",
      "\n",
      "# Calculations \n",
      "p = w*a*h\n",
      "h11 = (z*z)/(12*h)+h\n",
      "f = p*0.5072\n",
      "\n",
      "# Results \n",
      "print \"force normal to the gate at point B\",round(f,2),\"N\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "force normal to the gate at point B 114439.54 N\n"
       ]
      }
     ],
     "prompt_number": 10
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.24 Page No : 29"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "# Variables\n",
      "x = 30.\n",
      "z = math.sin(math.pi*x/180)\n",
      "d = 1.4\n",
      "h = 3.\n",
      "b = 1.5\n",
      "\n",
      "# Calculations \n",
      "h1 = z+d\n",
      "a = 0.5*h*b\n",
      "w = 9810\n",
      "p = w*a*h1\n",
      "h11 = ((z*z*h*h*h*b)/(36*a*h1))+h1\n",
      "\n",
      "# Results \n",
      "print \"total pressure on the plate\",round(p,2),\"N\"\n",
      "print \"position of centre of pressure\",round(h11,4),\"m\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "total pressure on the plate 41937.75 N\n",
        "position of centre of pressure 1.9658 m\n"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.25 Page No : 30"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "# Variables\n",
      "d = 1.8\n",
      "h = 2.4\n",
      "w = 9810\n",
      "s = 0.8\n",
      "\n",
      "# Calculations \n",
      "p1 = w*d*d*h*0.25*math.pi\n",
      "h1 = ((d*d)/(16*h))+h\n",
      "p = w*(s*1.5+2.4)\n",
      "p2 = p*math.pi*d*d*0.25\n",
      "p = p2-p1\n",
      "ab = w*(s*1.5+1.5)\n",
      "de = w*(s*1.5+3.3)\n",
      "ce = de-ab\n",
      "x = ((0.5*ce*d*0.3)/(0.5*(ab+de)*d))\n",
      "h2 = x+h\n",
      "h12 = h1-h2\n",
      "\n",
      "# Results \n",
      "print  \"change in total pressure\",round(p,3)\n",
      "print \"position of centre of pressure\",h2\n",
      "print \"change in position of centre of pressure\",round(h12,6),\"m\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "change in total pressure 29956.091\n",
        "position of centre of pressure 2.475\n",
        "change in position of centre of pressure 0.009375 m\n"
       ]
      }
     ],
     "prompt_number": 15
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.26 Page No : 31"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "\n",
      "# Variables\n",
      "l = 5.\n",
      "r = 3.\n",
      "a = l*r\n",
      "h = r*0.5\n",
      "w = 9810.\n",
      "\n",
      "# Calculations \n",
      "ph = w*a*h\n",
      "pv = w*0.25*math.pi*r*r*l\n",
      "p = math.sqrt((ph*ph)+(pv*pv))\n",
      "z = ph/pv\n",
      "theta = math.degrees(math.atan(z))\n",
      "\n",
      "# Results \n",
      "print \"resultant pressure on the gate\",round(p,2),\"N\"\n",
      "print \"angle of resultant force with vertical\",round(theta,7)\n",
      "\n",
      "# note : rounding off error."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "resultant pressure on the gate 411011.12 N\n",
        "angle of resultant force with vertical 32.4816366\n"
       ]
      }
     ],
     "prompt_number": 18
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 1.27 Page No : 32"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "# Variables\n",
      "s = 5.\n",
      "z = math.sin(math.radians(45.))\n",
      "a = 2*s*z\n",
      "h = s*z\n",
      "w = 9810.\n",
      "\n",
      "# Calculations \n",
      "ph = w*a*h\n",
      "pv = w*((0.25*s*s*math.pi)-(0.5*a*h))\n",
      "\n",
      "# Results \n",
      "print \"horizontal pressure\",ph,\"N\"\n",
      "print \"vertical pressure\",round(pv,1),\"N\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "horizontal pressure 245250.0 N\n",
        "vertical pressure 69993.9 N\n"
       ]
      }
     ],
     "prompt_number": 19
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}