{
 "metadata": {
  "name": "",
  "signature": "sha256:07d8c9662c407447473b639e8ef72a48c0966dd22d390e447be5846d8f6256c8"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Ch:10 Pipes and pipe joints"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "exa 10-1 - Page 295"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from math import sqrt\n",
      "sigta=140/2#\n",
      "nt=0.75#\n",
      "#Let the flow rate be Q\n",
      "Q=0.25#\n",
      "v=1.2#\n",
      "D=1.13*sqrt(Q/v)#\n",
      "D=520#\n",
      "p=0.7#\n",
      "C=9#\n",
      "t=(p*D)/(2*sigta*nt)+C#\n",
      "print \" t is %0.1f mm     \"%(t)#"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " t is 12.5 mm     \n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "exa 10-2 - Page 295"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from math import pi\n",
      "p=3*8#\n",
      "sigta=60#\n",
      "d=150#\n",
      "t=d/2*sqrt(((sigta+p)/(sigta-p))-1)#\n",
      "t=75*sqrt((84/36)-1)#\n",
      "t=40#\n",
      "do=d+(2*t)#\n",
      "D=d+(2*t)+20#\n",
      "w=10#\n",
      "Ds=d+(2*w)#\n",
      "P=pi*(Ds**2)*8/4#\n",
      "sigp=310#\n",
      "FOS=4#\n",
      "sigb=77.5#\n",
      "At=P/(sigb*2)#\n",
      "At=1300#\n",
      "D=250#\n",
      "db=45#\n",
      "b=D#\n",
      "a=1.8*b#\n",
      "CD=D+(2*db*1.2)#\n",
      "sigp=310#\n",
      "Pr=0.75*sigp*At#\n",
      "Pr=Pr*10**-3#\n",
      "t=40#\n",
      "D1=d+(2*t)+20#\n",
      "D2=D1+(4.6*31)#\n",
      "CD=D2-((3*t)+20)#\n",
      "print \" Pr is %0.2f kN     \"%(Pr)#\n",
      "print \"\\n D1 is %0.0f mm     \"%(D1)#\n",
      "print \"\\n D2 is %0.1f mm     \"%(D2)#\n",
      "print \"\\n CD is %0.1f mm     \"%(CD)#"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " Pr is 302.25 kN     \n",
        "\n",
        " D1 is 250 mm     \n",
        "\n",
        " D2 is 392.6 mm     \n",
        "\n",
        " CD is 252.6 mm     \n"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "exa 10-3 - Page 296"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "p=14#\n",
      "d=50#\n",
      "sigyp=270#\n",
      "FOS=3#\n",
      "sigta=sigyp/FOS#\n",
      "pt=2*p#\n",
      "t=d/2*sqrt(((sigta+pt)/(sigta-pt))-1)#\n",
      "t=10#\n",
      "D1=d+(2*t)#\n",
      "Ds=D1+20#\n",
      "P=pi*(Ds**2)*p/4#\n",
      "sigba=380/4#\n",
      "At=P/(4*sigba)#\n",
      "At=245#\n",
      "db=20#\n",
      "Dd=70+(2*20)+5#\n",
      "R=db+2.5#\n",
      "B=(Dd/sqrt(2))+(2*(db+2.5))#\n",
      "B=127#\n",
      "Y=Dd/(2*sqrt(2))#\n",
      "Rm=34.12#\n",
      "M=(P*Y/2)+(P*Rm/pi)#\n",
      "sigfa=250/5#\n",
      "b=127/70#\n",
      "Z=b/6#\n",
      "tf=sqrt(M/(sigfa*Z))#\n",
      "tf=44#\n",
      "print \" d is %0.0f mm     \"%(d)#\n",
      "print \"\\n t is %0.0f mm     \"%(t)#\n",
      "print \"\\n B is %0.0f mm     \"%(B)#\n",
      "print \"\\n R is %0.1f mm     \"%(R)#\n",
      "print \"\\n Y is %0.2f mm     \"%(Y)#\n",
      "print \"\\n tf is %0.0f mm     \"%(tf)#"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " d is 50 mm     \n",
        "\n",
        " t is 10 mm     \n",
        "\n",
        " B is 127 mm     \n",
        "\n",
        " R is 22.5 mm     \n",
        "\n",
        " Y is 40.66 mm     \n",
        "\n",
        " tf is 44 mm     \n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "exa 10-4 - Page 297"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from math import tan\n",
      "p=1.25#\n",
      "D=200#\n",
      "nt=0.75#\n",
      "C=9#\n",
      "sigta=20#\n",
      "t=(p*D)/(2*sigta*nt)+C#\n",
      "t=18#\n",
      "D1=D+(2*t)#\n",
      "dr=D1+10#\n",
      "sigp=310#\n",
      "sigba=sigp/4#\n",
      "db=16#\n",
      "Db=dr+32+5#\n",
      "Do=Db+(2*db)#\n",
      "P=pi*(251+db)**2*1.25/4#\n",
      "n=6#\n",
      "Y=(Db-dr)/2#\n",
      "M=P/n*Y#\n",
      "Z=dr*tan(30*pi/180)/6#\n",
      "tf=sqrt(M/(sigta*Z))#\n",
      "tf=22#\n",
      "Deff=dr+db+5#\n",
      "print \" D is %0.0f mm     \"%(D)#\n",
      "print \"\\n t is %0.0f mm     \"%(t)#\n",
      "print \"\\n Y is %0.1f mm     \"%(Y)#\n",
      "print \"\\n tf is %0.0f mm     \"%(tf)#\n",
      "print \"\\n Deff is %0.0f mm     \"%(Deff)#"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " D is 200 mm     \n",
        "\n",
        " t is 18 mm     \n",
        "\n",
        " Y is 18.5 mm     \n",
        "\n",
        " tf is 22 mm     \n",
        "\n",
        " Deff is 267 mm     \n"
       ]
      }
     ],
     "prompt_number": 4
    }
   ],
   "metadata": {}
  }
 ]
}