diff options
Diffstat (limited to 'Strength_Of_Materials_by_B_K_Sarkar/Chapter15.ipynb')
| -rwxr-xr-x | Strength_Of_Materials_by_B_K_Sarkar/Chapter15.ipynb | 657 |
1 files changed, 657 insertions, 0 deletions
diff --git a/Strength_Of_Materials_by_B_K_Sarkar/Chapter15.ipynb b/Strength_Of_Materials_by_B_K_Sarkar/Chapter15.ipynb new file mode 100755 index 00000000..8ae049c0 --- /dev/null +++ b/Strength_Of_Materials_by_B_K_Sarkar/Chapter15.ipynb @@ -0,0 +1,657 @@ +{
+ "metadata": {
+ "name": "chapter 15.ipynb"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter no.15:Thin Cyclindrical Shell"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.1,Page no.351"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "D=0.8 #m #Diameter of Shell\n",
+ "L=3 #m #Length of shell\n",
+ "t=0.01 #m #thickness of metal\n",
+ "E=200*10**9 #Pa \n",
+ "p=2.5*10**6 #Pa #Internal Pressure\n",
+ "m=4 #Poisson's ratio\n",
+ "\n",
+ "#Calculation\n",
+ "\n",
+ "sigma_1=p*D*(2*t)**-1 #N/m**2 #Hoop stress\n",
+ "sigma_2=p*D*(4*t)**-1 #N/m**2 #Longitudinal stress\n",
+ "\n",
+ "e_1=1*E**-1*(sigma_1-sigma_2*m**-1) #Hoop strain\n",
+ "e_2=1*E**-1*(sigma_2-sigma_1*m**-1) #Hoop strain\n",
+ "\n",
+ "d=e_1*D*100 #cm #Increase in Diameter\n",
+ "l=e_2*L*100 #cm #Increase in Length\n",
+ "\n",
+ "dell_v=2*e_1+e_2 #Volumetric strain\n",
+ "V=dell_v*pi*4**-1*D**2*L*10**6 #cm**3 #Increase in Volume\n",
+ "\n",
+ "#Result\n",
+ "print\"Change in Diameter is\",round(d,3),\"cm\"\n",
+ "print\"Change in Length is\",round(l,3),\"cm\"\n",
+ "print\"Change in Volume is\",round(V,2),\"cm**3\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Change in Diameter is 0.035 cm\n",
+ "Change in Length is 0.037 cm\n",
+ "Change in Volume is 1507.96 cm**3\n"
+ ]
+ }
+ ],
+ "prompt_number": 13
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.2,Page no.352"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "D=0.8 #m #iameter of water main\n",
+ "h=100 #m #Pressure head\n",
+ "w=10*10**3 #N/m**3 #Weight of Water\n",
+ "sigma_t=20*10**6 #MPa #Permissible stress\n",
+ "\n",
+ "#Calculation\n",
+ "\n",
+ "p=w*h #N/m**2 #Pressure of inside the main\n",
+ "t=p*D*(2*sigma_t)**-1*100 #m #Thcikness of metal\n",
+ "\n",
+ "#Result\n",
+ "print\"The Thickness of metal is\",round(t,2),\"cm\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The Thickness of metal is 2.0 cm\n"
+ ]
+ }
+ ],
+ "prompt_number": 21
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.3,Page no.352"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "p=2*10**6 #MPa #Steam Pressure\n",
+ "t=0.02 #m #thickness of boiler plate\n",
+ "sigma_t=120*10**6 #MPa #Tensile stress\n",
+ "sigma_l=120*10**6 #MPa #Longitudinal stress\n",
+ "rho=0.90 #% #Efficiency of Longitudinal joint\n",
+ "rho_e=0.40 #% #Efficiency of circumferential joint\n",
+ "\n",
+ "#Calculations\n",
+ "\n",
+ "D_1=sigma_t*2*t*rho*p**-1 #Diameter of boiler \n",
+ "D_2=sigma_l*4*t*rho_e*p**-1 #Diameter of boiler \n",
+ "\n",
+ "#Max diameter of boiler is equal to minimum value of diameter\n",
+ "\n",
+ "#Result\n",
+ "print\"Maximum diameter of boiler is\",round(D_2,2),\"m\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Maximum diameter of boiler is 1.92 m\n"
+ ]
+ }
+ ],
+ "prompt_number": 29
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.4,Page no.352"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "L=0.9 #m #Length of cyclindrical shell\n",
+ "D=0.2 #m #Internal Diameter\n",
+ "t=0.008 #m #thickness of metal\n",
+ "dV=20*10**-6 #m**3 #Additional volume\n",
+ "E=200*10**9 #Pa \n",
+ "m=1*0.3**-1 #Poissoin's ratio\n",
+ "\n",
+ "#Calculations\n",
+ "\n",
+ "V=pi*4**-1*D**2*L #Volume of cyclinder\n",
+ "\n",
+ "#Let X=2*e_1+e_2\n",
+ "X=dV*V**-1 #Volumetric strain (Equation 1)\n",
+ "\n",
+ "#e_1=p*D*(2*E*t)**-1*(1-1*(2*m)**-1) #Circumferential strain\n",
+ "#e_2=p*D*(2*E*t)**-1*(1*2**-1-1*(2*m)**-1) #Circumferential strain\n",
+ "\n",
+ "#substituting above values in equation 1 we get\n",
+ "p=X*E*t*(D*((1-1*(2*m)**-1)+(1*4**-1-1*(2*m)**-1)))**-1*10**-3 #KN/m**2 #Pressure exerted by fluid\n",
+ "sigma_t=p*D*(2*t)**-1 #KN/m**2 #hoop stress\n",
+ "\n",
+ "#Result\n",
+ "print\"Pressure Exerted by Fluid on the cyclinder is\",round(p,2),\"KN/m**2\"\n",
+ "print\"Hoop stress is\",round(sigma_t,2),\"KN/m**2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Pressure Exerted by Fluid on the cyclinder is 5956.68 KN/m**2\n",
+ "Hoop stress is 74458.45 KN/m**2\n"
+ ]
+ }
+ ],
+ "prompt_number": 41
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.5,Page no.353"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "t=0.015 #m #Thickness of plate\n",
+ "sigma_t=120*10**6 #Pa #tensile stress\n",
+ "sigma_l=120*10**6 #Pa #Longitudinal stress\n",
+ "rho=0.7 #% #Efficiency of longitudinal joints\n",
+ "rho_l=0.3 #% #Efficiency of circumferential joints\n",
+ "p=2*10**6 #Pa #Internal pressure\n",
+ "D=1.5 #m #shell diameter\n",
+ "\n",
+ "#Calculations (Part-1)\n",
+ "\n",
+ "D_1=sigma_t*2*t*rho*p**-1 #m \n",
+ "D_2=sigma_l*4*t*rho_l*p**-1 #m \n",
+ "\n",
+ "#Thus max diameter of shell is min of above two cases\n",
+ "\n",
+ "#Calculations (Part-2)\n",
+ "\n",
+ "p_1=sigma_t*2*t*rho*D**-1*10**-6 #MPa\n",
+ "p_2=sigma_l*4*t*rho_l*D**-1*10**-6 #MPa\n",
+ "\n",
+ "#Thus Internal pressure is min of above two cases\n",
+ "\n",
+ "#Result\n",
+ "print\"Max Permissible Diameter of shell is\",round(D_2,2),\"m\"\n",
+ "print\"Max Permissible Internal Pressure is\",round(p_2,2),\"MPa\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Max Permissible Diameter of shell is 1.08 m\n",
+ "Max Permissible Internal Pressure is 1.44 MPa\n"
+ ]
+ }
+ ],
+ "prompt_number": 8
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.6,Page no.354"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "L=3 #m #Length\n",
+ "D=1 #m #Internal Diameter\n",
+ "t=0.015 #m #thickness\n",
+ "p=1.5*10**6 #Pa #Internal pressure\n",
+ "E=200*10**9 #Pa \n",
+ "m=1*0.3**-1 #Poissoin's ratio\n",
+ "\n",
+ "#Calculations\n",
+ "\n",
+ "sigma_t=p*D*(2*t)**-1*10**-6 #MPa #Hoop stress\n",
+ "sigma_l=p*D*(4*t)**-1*10**-6 #MPa #Longitudinal stress\n",
+ "\n",
+ "dD=(p*D**2*(2*t*E)**-1*(1-1*(2*m)**-1))*10**2 #cm #Change in Diameter\n",
+ "dL=p*D*L*(2*t*E)**-1*(1*2**-1-1*m**-1)*10**2 #cm #Change in Length\n",
+ "\n",
+ "V=pi*4**-1*D**2*L #Volume \n",
+ "dV=p*D*(2*t*E)**-1*(5*2**-1-2*(m)**-1)*V*10**6 #cm #Change in Volume\n",
+ "\n",
+ "#Result\n",
+ "print\"The circumferential stresses induced is\",round(sigma_t,2),\"MPa\"\n",
+ "print\"The Longitudinal stresses induced is\",round(sigma_l,2),\"MPa\"\n",
+ "print\"The change in dimension are:D is\",round(dD,3),\"cm\"\n",
+ "print\" :L is\",round(dL,4),\"cm\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "0.02125\n",
+ "The circumferential stresses induced is 50.0 MPa\n",
+ "The Longitudinal stresses induced is 25.0 MPa\n",
+ "The change in dimension are:D is 0.021 cm\n",
+ " :L is 0.015 cm\n"
+ ]
+ }
+ ],
+ "prompt_number": 26
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.7,Page no.355"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "L=0.9 #m #Length of cyclinder\n",
+ "D=0.4 #m #Diameter \n",
+ "t=0.006 #m #thickness\n",
+ "p=5*10**6 #Pa #Pressure\n",
+ "E=100*10**9\n",
+ "m=3 #Poissoin's ratio\n",
+ "k=2.6*10**9 #Pa #Bulk modulus\n",
+ "\n",
+ "#Calculations\n",
+ "\n",
+ "#Let X=dV_1*V_1**-1\n",
+ "X=p*(0.4-2*0.006)*(2*t*E)**-1*(5*2**-1-2*m**-1) #Volumetric strain\n",
+ "dV_1=round(X,5)*pi*4**-1*0.388**2*L #cm**3 #Increase in volume of cyclinder\n",
+ "V_1=pi*4**-1*0.388**2*L #VOlume\n",
+ "dV_2=p*k**-1*V_1 #DEcrease in volume of oil due to increase in pressure\n",
+ "\n",
+ "dV=(dV_1+dV_2)*10**6 #Resultant additional space \n",
+ "\n",
+ "#Result\n",
+ "print\"additional quantity of oil to be pumped is\",round(dV,2),\"cm**3\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "additional quantity of oil to be pumped is 519.62 cm**3\n"
+ ]
+ }
+ ],
+ "prompt_number": 23
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.8,Page no.356"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "A=1600*(3600)**-1 #Kg/sec #Amount of steam generated\n",
+ "v=0.24 #m**3/kg #specific volume of steam\n",
+ "sigma_t=4*10**6 #MPa #Tensile stress\n",
+ "V_1=30 #m/s #Velocity of steam\n",
+ "p=1*10**6 #Pa #Steam pressure\n",
+ "\n",
+ "#Calculation\n",
+ "\n",
+ "V=A*v #m**3/s #volume of steam \n",
+ "D=(V*(pi*4**-1*V_1)**-1)**0.5*100 #Diameter of pipe\n",
+ "t=p*D*(2*sigma_t)**-1 #Thicknes of pipe\n",
+ "\n",
+ "#Result\n",
+ "print\"Diameter of boiler is\",round(D,2),\"cm\"\n",
+ "print\"Thickness of steel plpe is\",round(t,2),\"cm\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Diameter of boiler is 6.73 cm\n",
+ "Thickness of steel plpe is 0.84 cm\n"
+ ]
+ }
+ ],
+ "prompt_number": 43
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.9,Page no.359"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "P=14*10**3 #N #Axial pull\n",
+ "dL=0.0084 #cm #Elongation\n",
+ "L=0.25 #m #Length\n",
+ "p=7*10**6 #Internal pressure\n",
+ "dL_2=0.0034 #cm #Longation\n",
+ "d=0.0475 #m #Internal diameter \n",
+ "D=0.05 #m #External Diameter\n",
+ "m=0.25\n",
+ "\n",
+ "#Calculation\n",
+ "\n",
+ "t=(D-d)*2**-1 #thickness od tube\n",
+ "A=pi*4**-1*(D**2-d**2) #Area of tube\n",
+ "sigma=P*A**-1 #stress\n",
+ "e=dL*(L)**-1 #strain\n",
+ "E=sigma*e**-1 #Modulus of Elasticity\n",
+ "sigma_1=p*d*(2*t)**-1 #Hoop stress\n",
+ "sigma_2=p*d*(4*t)**-1 #Longitudinal stress\n",
+ "\n",
+ "m=-(sigma_1*(dL_2*L**-1*E-sigma_2)**-1) #POissoin's ratio\\\n",
+ "\n",
+ "#Let X=1*m**-1\n",
+ "X=1*m**-1 #Poissoin's ratio\n",
+ "\n",
+ "#Result\n",
+ "print\"The value of Poissoin's ratio is\",round(X,3)\n",
+ " "
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The value of Poissoin's ratio is 0.277\n"
+ ]
+ }
+ ],
+ "prompt_number": 28
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.10,Page no.357"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "D=0.8 #m #Diameter\n",
+ "t=0.01 #m #Thickness \n",
+ "p=5*10**6 #Pa #Pressure\n",
+ "m=1*0.25**-1\n",
+ "E=200*10**9 #Pa\n",
+ "\n",
+ "#Calculations\n",
+ "\n",
+ "sigma_1=5*10**6*0.8*(4*0.01)**-1 #stress\n",
+ "sigma_2=sigma_1\n",
+ "e_1=sigma_1*E**-1-sigma_2*(m*E)**-1 #strain\n",
+ "e_v=3*e_1\n",
+ "V=4*3**-1*pi*(D*2**-1)**3 #m**3 tress\n",
+ "dell_v=e_v*V*10**6 #cm**3\n",
+ "\n",
+ "#Result\n",
+ "print\"Volume of additional Fluid\",round(dell_v,3),\"cm**3\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Volume of additional Fluid 301.593 cm**3\n"
+ ]
+ }
+ ],
+ "prompt_number": 37
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.11,Page no.358"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "d=0.3 #m #Diameter\n",
+ "D=0.003 #m #Diameter of steel wire\n",
+ "t=0.006 #m #thickness\n",
+ "sigma_w=8*10**6 #Pa #Stress\n",
+ "p=1*10**6 #Pa #Internal pressure\n",
+ "E_s=200*10**9 #Pa #Modulus of Elasticity for steel\n",
+ "E_c=100*10**9 #Pa #Modulus of Elasticity for cast iron\n",
+ "m=1*0.3**-1\n",
+ "\n",
+ "#Calculations\n",
+ "\n",
+ "sigma_p=(sigma_w*pi*2**-1*d)*(2*t)**-1 #compressive hoop stress\n",
+ "sigma_l=p*d*(4*t)**-1 #Longitudinal stress\n",
+ "\n",
+ "#when internal presure is apllied Let sigma_w_1=Tensile in wire and sigma_p_1=tensile hoop in wire\n",
+ "#sigma_p_1*2*t+sigma_w_1*2*d**-1*pi*4**-1*d**2=p*D\n",
+ "\n",
+ "#After substituting values and further simplifying we get\n",
+ "#1.2*sigma_p_1+0.471*sigma_w_1=3000 Equation 1\n",
+ "\n",
+ "#1*E_c**-1(sigma_p_1-sigma_1*m**-1+sigma_p)=1*E_s**-1(sigma_w_1-sigma_w)\n",
+ "\n",
+ "#After substituting values and further simplifying we get\n",
+ "#sigma_p_1-0.5*sigma_w_1=1.36*10**6 \n",
+ "#sigma_p_1=0.5*sigma_w_1-3.39*10**6 Equation 2\n",
+ "\n",
+ "#From Equation 2 substituting value of sigma_p_1 in Equation 1\n",
+ "\n",
+ "\n",
+ "sigma_w_1=(40.68*10**3+0.3*10**6)*(10.71238*10**-3)**-1\n",
+ "sigma_p_1=0.5*sigma_w_1-3.39*10**6\n",
+ "\n",
+ "#Let X=sigma_p_1 and Y=sigma_w_1\n",
+ "X=sigma_p_1*10**-6 #MPa #Stresses in pipe\n",
+ "Y=sigma_w_1*10**-6 #MPa #Stresses in wire\n",
+ "\n",
+ "#Result\n",
+ "print\"Stress in the pipe is\",round(X,2),\"MN/m**2\"\n",
+ "print\"Stress in the wire is\",round(Y,2),\"MN/m**2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Stress in the pipe is 12.51 MN/m**2\n",
+ "Stress in the wire is 31.8 MN/m**2\n"
+ ]
+ }
+ ],
+ "prompt_number": 29
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Problem no 15.12,Page no.359"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#Initilization of variables\n",
+ "\n",
+ "D=0.038 #m #External Diameter\n",
+ "d=0.035 #m #Internal Diameter\n",
+ "d_1=0.0008 #m #Steel wire diameter\n",
+ "p=2*10**6 #pa #Pa #Internal Pressure\n",
+ "sigma_t_1=7*10**6 #Pa #Circumferential stress\n",
+ "#E_s=1.6*E_s\n",
+ "m=0.3\n",
+ "\n",
+ "#Calculation\n",
+ "\n",
+ "t=(D-d)*2**-1 #m Thickness \n",
+ "\n",
+ "#sigma_t*2*t=pi*d*2**-1*sigma_w\n",
+ "#From Above equation we get\n",
+ "\n",
+ "#sigma_t=0.419*sigma_w (Equation 1)\n",
+ "\n",
+ "sigma_w_1=(p*d-sigma_t_1*2*t)*(2*d_1**-1*pi*4**-1*d_1**2)**-1 #stress in wire\n",
+ "sigma_l=p*d*(4*t)**-1 #Longitudinal stress in tube\n",
+ "\n",
+ "#Now Equating equations of strain in tube and wire we get\n",
+ "sigma_w=-(1.6*(sigma_t_1-sigma_l*m)-sigma_w_1)*1.67**-1*10**-6\n",
+ "\n",
+ "#Result\n",
+ "print\"The Tension at which wire must have been wound is\",round(sigma_w,2),\"MPa\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The Tension at which wire must have been wound is 20.0 MPa\n"
+ ]
+ }
+ ],
+ "prompt_number": 30
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |