{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Ch:31 Seals packing and gaskets"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## exa 31-1 - Page 816"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "N is 4.47      \n"
     ]
    }
   ],
   "source": [
    "from __future__ import division\n",
    "from math import sqrt, pi, log\n",
    "d=18#\n",
    "lg=25+25#\n",
    "Eb=210*10**3#\n",
    "Ecl=90*10**3#\n",
    "A=pi*d**2/4#\n",
    "kb=A*Eb/lg#\n",
    "x=(5*(lg+(0.5*d))/(lg+(2.5*d)))#\n",
    "km=pi*Ecl*d/(2*log(x))#\n",
    "C=kb/(kb+km)#\n",
    "sigp=600#\n",
    "At=192#\n",
    "Pi=0.75*sigp*At#\n",
    "F=200#\n",
    "C=0.322#\n",
    "Pb=F*C*10**3#\n",
    "FOS=2#\n",
    "W=At*sigp#\n",
    "N=Pb*FOS/(W-Pi)#\n",
    "print \"N is %0.2f      \"%(N)#"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## exa 31-2 - Page 816"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "p is 6.922 N/mm**2     \n"
     ]
    }
   ],
   "source": [
    "from math import sqrt, pi, log\n",
    "d=16#\n",
    "D=1.5*d#\n",
    "t=20#\n",
    "tg=4#\n",
    "#Let Gasket diameter in compression zone be d1\n",
    "d1=D+(2*t)+tg#\n",
    "lg=40#\n",
    "E=207*10**3#\n",
    "kb=pi*d**2*E/(lg*4)#\n",
    "Ecl=90*10**3#\n",
    "x=(5*(lg+(0.5*d))/(lg+(2.5*d)))#\n",
    "kp=pi*Ecl*d/(2*log(x))#\n",
    "Ag=pi*(d1**2-d**2)/4#\n",
    "Eg=480#\n",
    "kg=Ag*Eg/tg#\n",
    "km=kg*kp/(kg+kp)#\n",
    "C=kb/(kb+km)#\n",
    "At=157#\n",
    "sigp=600#\n",
    "Pi=0.75*At*sigp/2#\n",
    "FOS=2#\n",
    "Pf=At*sigp/FOS#\n",
    "W=Pf-Pi#\n",
    "P=W/C#\n",
    "N=5#\n",
    "F=P*N#\n",
    "p=F*4/(pi*120**2)#\n",
    "print \"p is %0.3f N/mm**2     \"%(p)#"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## exa 31-3 - Page 817"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "d is 20 mm     \n"
     ]
    }
   ],
   "source": [
    "from numpy import mat, pi, sqrt\n",
    "sigp=600#\n",
    "FOS=3#\n",
    "siga=sigp/FOS#\n",
    "d=16#\n",
    "D=1.5*d+60#\n",
    "#Let Gasket diameter in compression zone be d1\n",
    "d1=(300-160)/2#\n",
    "#Let compressive stress in gasket for leak proof joint be sigl\n",
    "sigl=12#\n",
    "At=mat([[1, 157],[2 ,192], [3, 245]])\n",
    "d=mat([[1 ,16],[2 ,18],[3, 20]])\n",
    "\n",
    "n=3#\n",
    "Pi = range(0,n)\n",
    "Pc = range(0,n)\n",
    "for i in range(0,n):\n",
    "    Pi[i]=At[i,1]*d[i,1]\n",
    "    Pc[i]=3*pi*(d1**2-d[i,1]**2)\n",
    "    if (Pi[i]>=Pc[i]):\n",
    "        print \"The Design is safe\"\n",
    "\n",
    "print \"d is %0.0f mm     \"%(d[i,1])#"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}