diff options
Diffstat (limited to 'Statics_And_Strength_Of_Materials')
19 files changed, 5628 insertions, 0 deletions
diff --git a/Statics_And_Strength_Of_Materials/README.txt b/Statics_And_Strength_Of_Materials/README.txt new file mode 100755 index 00000000..72f46bc3 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/README.txt @@ -0,0 +1,10 @@ +Contributed By: jaydeep palekar +Course: btech +College/Institute/Organization: Freelancer +Department/Designation: Computer Science +Book Title: Statics And Strength Of Materials +Author: I. J. Levinson +Publisher: Prentice Hall Inc. +Year of publication: 2002 +Isbn: 9780138445065 +Edition: 20
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch1.ipynb b/Statics_And_Strength_Of_Materials/ch1.ipynb new file mode 100755 index 00000000..94d80818 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch1.ipynb @@ -0,0 +1,256 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:e8a178cfe4f176b57d4a8a5f7fd1354a20d5cbe9cf725d0ac6feec0b0641ec3e" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 1 : Introduction" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.1 Page No : 3" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\n", + "# Variables\n", + "L = 20. \t\t\t#ft\n", + "angle = 30.\t\t\t#degrees\n", + "\n", + "# Calculations\n", + "d = L*math.sin(math.radians(angle))\n", + "\n", + "# Results\n", + "print 'Desitance from foot of Ladder = %.2f ft'%(d)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Desitance from foot of Ladder = 10.00 ft\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.2 Page No : 5" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "\n", + "# Variables\n", + "a = 5.\n", + "b = 12.\n", + "angle = 60. \t\t\t#degrees\n", + "\n", + "# Calculations\n", + "c = math.sqrt(a**2+b**2-2*a*b*math.cos(math.radians(angle)))\n", + "\n", + "# Results\n", + "print 'c = %.1f '%(c)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "c = 10.4 \n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.3 Page No : 5" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "import math \n", + "\n", + "# Variables\n", + "a = 5.\n", + "b = 12.\n", + "angle = 120. \t\t\t#degrees\n", + "\t\t\t\n", + "# Calculations\n", + "c = math.sqrt(a**2+b**2-2*a*b*math.cos(math.radians(angle)))\n", + "\t\t\t\n", + "# Results\n", + "print 'c = %.1f '%(c)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "c = 15.1 \n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.4 Page No : 6" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "b = 12.\n", + "angle1 = 35. \t\t\t#degrees\n", + "angle2 = 43. \t\t\t#degrees\n", + "\t\t\t\n", + "# Calculations\n", + "angle3 = 180-angle1-angle2\n", + "a = math.sin(math.radians(angle2))*b/math.sin(math.radians(angle3))\n", + "c = a*math.sin(math.radians(angle1))/math.sin(math.radians(angle2))\n", + "\t\t\t\n", + "# Results\n", + "print 'c = %.2f '%(c)\n", + "print 'a = %.2f.'%(a)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "c = 7.04 \n", + "a = 8.37.\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.5 Page No : 7" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Wofaninch = 0.29 \t\t\t#lb\n", + "L = 3.5 \t \t\t#ft\n", + "width = 1.75 \t\t \t#ft\n", + "t = 1. \t \t\t #in\n", + "\t\t\t\n", + "# Calculations\n", + "W = L*width*t*12*12*Wofaninch\n", + "\t\t\t\n", + "# Results\n", + "print 'W = %.f lb'%(W)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "W = 256 lb\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.6 Page No : 7" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "V = 30. \t\t\t#mph\n", + "\t\t\t\n", + "# Calculations\n", + "Vinfps = V*5280*(1./60)*(1./60)\n", + "\t\t\t\n", + "# Results\n", + "print 'v = %.f fps'%(Vinfps)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "v = 44 fps\n" + ] + } + ], + "prompt_number": 6 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch10.ipynb b/Statics_And_Strength_Of_Materials/ch10.ipynb new file mode 100755 index 00000000..b18a635e --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch10.ipynb @@ -0,0 +1,500 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:fa8306a7ccf6683f3bb43ba03db574b371457c13bfcbb27905e40f6b1eacde3b" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 10 : Torsion" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.1 Page No : 229" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \t\t\t\n", + "\n", + "# Variables\n", + "L = 50. \t\t\t#ft\n", + "Do = 2 \t\t\t#in\n", + "Di = 1.5 \t\t\t#in\n", + "Mt = 10000. \t\t\t#lb in\n", + "G = 12.*10**6\n", + "\t\t\t\n", + "# Calculations\n", + "Tmax = 16*Mt*Do/(math.pi*(Do**4-Di**4))\n", + "angle = (Mt*L*12*32)*57.3/(G*math.pi*(Do**4-Di**4))\n", + "\t\t\t\n", + "# Results\n", + "print 'Maximum shearing strees = %.f psi'%(round(Tmax,-1))\n", + "print 'twist angle = %.1f degrees'%(angle)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum shearing strees = 9310 psi\n", + "twist angle = 26.7 degrees\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.2 Page No : 229" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "d = 4. \t\t\t#ft\n", + "T = 5000. \t\t\t#psi\n", + "angle = 0.1 \t\t\t#degrees\n", + "\t\t\t\n", + "# Calculations\n", + "T1 = (math.pi*d**3)*T/16\n", + "T2 =angle*math.pi*G*math.pi*d**4/(180*12*32)\n", + "\t\t\t\n", + "# Results\n", + "if (T1<T2): \n", + " print 'Safe torque = %.2f lb in'%(T1)\n", + "else:\n", + " print 'Safe torque = %.2f lb'%(T2)\n", + "\n", + "# note : anwer is wrong in book. plz check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Safe torque = 43864.91 lb\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.3 Page No : 229" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from numpy import linalg\n", + "\t\t\t\n", + "# Variables\n", + "Ds = 1. \t\t\t#in\n", + "Db = 1.5 \t\t\t#in\n", + "Ls = 4. \t\t\t#in\n", + "Lb = 6. \t\t\t#in\n", + "Gs = 12.*10**6 \t\t\t#psi\n", + "Gb = 6.4*10**6 \t\t\t#psi\n", + "T = 10000. \t\t\t#lb in\n", + "\t\t\t\n", + "# Calculations\n", + "A = [[1,1],[(Ls*12/(Gs*Ds**4)),(-Lb*12/(Gb*Db**4))]]\n", + "b = [T,0]\n", + "c = linalg.solve(A,b)\n", + "Tab = c[0]\n", + "Tbc = c[1]\n", + "\t\t\t\n", + "# Results\n", + "print 'Torque in section AB = %.f lb in'%(Tab)\n", + "print 'Torque in section AB = %.f lb in'%(Tbc)\n", + "\n", + "# note : Answers are slightly different because of inbuilt solve function of python." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Torque in section AB = 3571 lb in\n", + "Torque in section AB = 6429 lb in\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.4 Page No : 230" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "T = 10000. \t\t\t#lb in\n", + "G = 12.*10**6\n", + "Dab = 1.5 \t\t\t#in\n", + "Lab = 4. \t\t\t#in\n", + "Dcd = 1. \t\t\t#in\n", + "Lcd = 3. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "F = T/2\n", + "Tab = F*Lab\n", + "angle = ((T*32*12*Lcd/(G*math.pi*Dcd**4))+2*(Tab*32*12*Lab/(G*math.pi*Dab**4)))*(180/math.pi)\n", + "\t\t\t\n", + "# Results\n", + "print 'angle of twist = %.0f degrees'%(angle)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "angle of twist = 36 degrees\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.5 Page No : 231" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Tallowable = 5000. \t\t\t#psi\n", + "power = 250. \t \t\t#hp\n", + "n = 1800. \t\t\t#rpm\n", + "\t\t\t\n", + "# Calculations\n", + "T = 63000*power/n\n", + "d = (16*T/(math.pi*Tallowable))**(1/3.)\n", + "\t\t\t\n", + "# Results\n", + "print 'Torque = %.2f lb in'%(T)\n", + "print 'diameter =%.2f in'%(d)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Torque = 8750.00 lb in\n", + "diameter =2.07 in\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.6 Page No : 232" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "ds = 2. \t\t\t#in\n", + "n = 315. \t\t\t#rpm\n", + "Gs = 12.*10**6\n", + "Lab = 5. \t\t\t#in\n", + "Lbc = 15. \t\t\t#in\n", + "Pa = 10. \t\t\t#hp\n", + "Pc = 40. \t\t\t#hp\n", + "Pb = 50. \t\t\t#hp\n", + "\t\t\t\n", + "# Calculations\n", + "Tab = 63000*Pa/n\n", + "Tbc = 63000*Pc/n\n", + "angle = ((32*Tbc*Lbc*12/(math.pi*ds**4*G))-(32*Tab*Lab*12/(math.pi*ds**4*G)))*(180/math.pi)\n", + "\t\t\t\n", + "# Results\n", + "print 'angle of twist of gear C releative to a = %.2f degrees'%(angle)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "angle of twist of gear C releative to a = 4.01 degrees\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.7 Page No : 234" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "k1 = 6.*10**6 \t\t\t#lb in/rad\n", + "k2 = 3.*10**6 \t\t\t#lb in/rad\n", + "k3 = 2.*10**6 \t\t\t#lb in/rad\n", + "T = 10000. \t\t\t#lb in\n", + "\t\t\t\n", + "# Calculations\n", + "ke = 1/((1/k1)+(1/k2)+(1/k3))\n", + "angle = T*180/(ke*math.pi)\n", + "\t\t\t\n", + "# Results\n", + "print 'equivalent spring constant = %.2e lb in/rad'%(ke)\n", + "print 'angle of twist d/a = %.2f degrees'%(angle)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "equivalent spring constant = 1.00e+06 lb in/rad\n", + "angle of twist d/a = 0.57 degrees\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.8 Page No : 234" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "k1 = 2.*10**6 \t\t\t#lb in/rad\n", + "k2 = 3.*10**6 \t\t\t#lb in/rad\n", + "T = 20000. \t\t\t#lb in\n", + "\t\t\t\n", + "# Calculations\n", + "ke = k1+k2\n", + "angle = T*180/(ke*math.pi)\n", + "\t\t\t\n", + "# Results\n", + "print 'equivalent spring consmath.tant = %.2e lb in/rad'%(ke)\n", + "print 'angle of twist at B = %.3f degrees'%(angle)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "equivalent spring consmath.tant = 5.00e+06 lb in/rad\n", + "angle of twist at B = 0.229 degrees\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.9 Page no : 238" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "# variables\n", + "n = 10 # coils\n", + "P = 1200. # axial load lb\n", + "R = 2. \n", + "K = 1.33 # factor\n", + "d = 1.\n", + "\n", + "# Calculations\n", + "Tmax = round(K*(16*P*R)/(math.pi*d**3),-2)\n", + "delta = 64*P*R**3*n/(12*10**6*d**4)\n", + "\n", + "# Results\n", + "print \"Stress = %d psi\"%Tmax\n", + "print \"The deflection = %.3f in\"%delta\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Stress = 16300 psi\n", + "The deflection = 0.512 in\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.10 Page No : 239" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "di = 0.2 \t\t\t#in\n", + "dm = 2. \t\t\t#in\n", + "n = 10.\n", + "F = 10. \t\t\t#lb\n", + "G = 12.*10**6\n", + "\t\t\t\n", + "# Calculations\n", + "k = G*di**4/(64*dm**3*n)\n", + "ke = 1/((1/(k+k))+(1/k)+(1/k))\n", + "delta = F/ke\n", + "\t\t\t\n", + "# Results\n", + "print 'elongation = %.2f in'%(delta)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "elongation = 6.67 in\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 10.11 Page No : 241" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "d = 0.5 \t\t\t#in\n", + "n = 315. \t\t\t#rpm\n", + "t1 = 5000. \t\t\t#psi\n", + "r1 = 8. \t\t\t#in\n", + "r2 = 4. \t\t\t#in \n", + "n1 = 6.\n", + "n2 = 4.\n", + "\t\t\t\n", + "# Calculations\n", + "t2 = r2*t1/r1\n", + "T = r1*n1*(math.pi/4)*d**2*t1+r2*n2*(math.pi/4)*d**2*t2\n", + "hp = T*n/63000\n", + "\t\t\t\n", + "# Results\n", + "print 'Premissible horsepower = %.f hp'%(hp)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Premissible horsepower = 275 hp\n" + ] + } + ], + "prompt_number": 10 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch12.ipynb b/Statics_And_Strength_Of_Materials/ch12.ipynb new file mode 100755 index 00000000..2d6e383c --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch12.ipynb @@ -0,0 +1,493 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:80d6cc3748bed03e9fd030c70419175672c3c9a4d62b5544fe9d0f908939d53e" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 12 : Stresses in Beams" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.1 Page No : 281" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "L = 20. \t\t\t#ft\n", + "b1 = 12. \t\t\t#in\n", + "h1 = 4. \t\t\t#in\n", + "b2 = 4. \t\t\t#in\n", + "h2 = 12. \t\t\t#in\n", + "Fs = 1200. \t\t\t#psi\n", + "La = 5. \t\t\t#ft\n", + "Lb = 15. \t\t\t#ft\n", + "\t\t\t\n", + "# Calculations\n", + "Ina = b1*h1**3/12\n", + "P1 = (Fs*Ina*4)/((h1/2)*12*La*3)\n", + "Ina1 = b2*h2**3/12\n", + "P2 = (Fs*Ina1*4)/((h2/2)*12*La*3)\n", + "\t\t\t\n", + "# Results\n", + "print 'P max in first case = %.0f lb'%(P1)\n", + "print 'P max in second case = %.2f lb'%(P2)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "P max in first case = 853 lb\n", + "P max in second case = 2560.00 lb\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.2 Page No : 282" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "b = 0.5 \t\t\t#in\n", + "h = 1./32 \t\t\t#in\n", + "d = 4. \t\t\t#ft\n", + "E = 30.*10**6\n", + "\t\t\t\n", + "# Calculations\n", + "stress = E*(h/2)/((d/2)*12)\n", + "Ina = b*h**3/12\n", + "M = stress*Ina/(h/2)\n", + "\t\t\t\n", + "# Results\n", + "print 'maximum stress = %.2f psi'%(round(stress,-2))\n", + "print 'internal moment = %.2f lb in'%(M) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "maximum stress = 19500.00 psi\n", + "internal moment = 1.59 lb in\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.3 Page No : 283" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "W = 1000. \t\t\t#lb/ft\n", + "L = 10. \t\t\t#in\n", + "b1 = 4. \t\t\t#in\n", + "h1 = 1. \t\t\t#in\n", + "b2 = 1. \t\t\t#in\n", + "h2 = 6. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Mmax = 12500 \t\t\t#lb ft\n", + "y = ((b1*h1*h1/2)+(b2*h2*((h2/2)+h1)))/(b1*h1+b2*h2)\n", + "Ina = round((b1*h1**3/12)+b1*h1*(y-h1/2)**2+(b2*h2**3/12)+b2*h2*(h1+h2-y-(h2/2))**2,1)\n", + "sigmat = Mmax*12*y/Ina\n", + "sigmac = Mmax*12*(h1+h2-y)/Ina\n", + "\n", + "\t\t\t\n", + "# Results\n", + "print 'maximum tensile stress = %.2f psi'%(round(sigmat,-1))\n", + "print 'maximum compressive bending stress = %.2f psi'%(round(sigmac,-2)) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "maximum tensile stress = 8180.00 psi\n", + "maximum compressive bending stress = 13800.00 psi\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.4 Page No : 287" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "st = 1200. \t\t\t#psi\n", + "sc = 100. \t\t\t#psi\n", + "h = 12. \t\t\t#in\n", + "b = 4. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "I = b*h**3/12\n", + "P1 = st*I/(b*12*(h/2))\n", + "P2 = 2*sc*b*12/3\n", + "if (P1<P2):\n", + " print 'Safe value of p = %.f lB'%(P1)\n", + "else: \n", + " print 'Safe value of p = %.f lB'%(P2)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Safe value of p = 2400 lB\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.5 Page No : 287" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "W = 600. \t\t\t#lb/ft\n", + "L1 = 8. \t\t\t#in\n", + "L2 = 4. \t\t\t#in\n", + "b = 6. \t\t\t#in\n", + "h = 8. \t\t\t#in\n", + "t = 1. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "R1 = W*(L1+L2)*((L1+L2)/2)/L1\n", + "R2 = W*(L1+L2)*(L1-(L1+L2)/2)/L1\n", + "Vmax = 3000 \t\t\t#lb\n", + "I = (b*h**3/12)-(L2*b**3/12)\n", + "Ay = b*L2*(L2/2)-L2*b/2*b/4\n", + "b = t+t\n", + "Tmax = Vmax*Ay/(I*b)\n", + "\t\t\t\n", + "# Results\n", + "print 'maximum shear stress = %.0f psi'%(Tmax) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "maximum shear stress = 245 psi\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.6 Page No : 290" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "w = 4000. \t\t\t#lb/ft\n", + "l = 20. \t\t\t#ft\n", + "y = 0.96\n", + "A = 4.18 \t\t\t#in**2\n", + "Icq = 5.6 \t\t\t#in**4\n", + "d = 28. \t\t\t#in\n", + "b = 0.5 \t\t\t#in\n", + "T = 8000. \t\t\t#psi\n", + "d1 = 0.75 \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "V = w*l/2\n", + "Ay = 2*A*((d/2)-y)\n", + "I = b*d**3/12+4*(Icq+A*((d/2)-y)**2)\n", + "p = (2*T*(math.pi/4)*d1**2*I)/(V*Ay)\n", + "\t\t\t\n", + "# Results\n", + "print 'Rivet spacing = %.2f in'%(p) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rivet spacing = 6.13 in\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.7 Page No : 293" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "Es = 30.*10**6\n", + "Ew = 1.5*10**6\n", + "w = 500. \t\t\t#lb per ft\n", + "span = 12. \t\t\t#ft\n", + "t = 0.25 \t\t\t#in\n", + "h = 12. \t\t\t#in\n", + "n = 3.\n", + "b = 5. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "bw = Es*t/Ew\n", + "Ina = n*b*h**3/12\n", + "M = (w*span*(h/2)*12)/4\n", + "S = M*(h/2)/I\n", + "Ss = Es*S/Ew\n", + "bs = Ew*bw/Es\n", + "Ina1 = n*t*h**3/12\n", + "Ss1 = M*(h/2)/Ina1\n", + "Sw = Ew*Ss1/Es\n", + "\t\t\t\n", + "# Results\n", + "print 'Maximum bending stress in steel = %.3f psi'%(Ss1)\n", + "print 'Maximum bending stress in wood = %.2f psi'%(Sw)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum bending stress in steel = 6000.000 psi\n", + "Maximum bending stress in wood = 300.00 psi\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.8 Page No : 294" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Ss = 15000. \t\t\t#psi\n", + "Sa = 6000. \t\t\t#psi\n", + "Es = 30.*10**6\n", + "Ea = 10.*10**6\n", + "Sl = 16. \t\t\t#ft\n", + "ba = 3. \t\t\t#in\n", + "ha = 8. \t\t\t#in\n", + "hs = 1. \t\t\t#in\n", + "b = 1. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations and Resuls\n", + "bs = (Ea/Es)*ba\n", + "Y = ((ba-b)*b*(hs/2)+(ha+b)*b*((ha/2)+(hs/2)))/(ba*b+ha*b)\n", + "I = (ba*hs**3/12)+ba*hs*(Y-(hs/2))**2+((b*ha**3/12)+b*ha*(ha-Y-(ha/2))**2)\n", + "w1 = Ss*I/(Y*(1./2)*ha*(ha)*12)\n", + "Ss = Es*Sa/Ea\n", + "w2 = Ss*I/((ha-Y)*(1./2)*ha*(ha)*12)\n", + "if (w1<w2):\n", + " print 'Greatest uniformly distributed load = %d lb per ft'%(w1)\n", + "else:\n", + " print 'Greatest uniformly distributed load = %d lb per ft'%(w2)\n", + "\n", + "# note : rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Greatest uniformly distributed load = 781 lb per ft\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.9 Page No : 297" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "M = 500000. \t\t\t#lb in\n", + "r = 15.\n", + "n = 3.\n", + "b = 20. \t\t\t#in\n", + "l = 12. \t\t\t#in\n", + "As = 1. \t\t\t#in**2\n", + "\t\t\t\n", + "# Calculations\n", + "At = r*As*n\n", + "x = (-2*At+math.sqrt((2*At)**2+8*At*b*l))/(2*l)\n", + "Ina = ((l*x**3)/3)+At*(b-x)**2\n", + "Scmax = M*x/Ina\n", + "Ssmax = r*M*(b-x)/Ina\n", + "\t\t\t\n", + "# Results\n", + "print 'Maximum bending stress in concrete = %.0f psi'%(Scmax)\n", + "print 'Maximum bending stress in steel = %.2f psi'%(Ssmax)\n", + "\n", + "# note : answer is accurate. please check with calculator." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum bending stress in concrete = 542 psi\n", + "Maximum bending stress in steel = 9815.22 psi\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 12.10 Page No : 298" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Sc = 800. \t\t\t#psi\n", + "Ss = 18000. \t\t\t#psi\n", + "ratio = 15.\n", + "d = 5/8. \t\t\t#in\n", + "l = 20. \t\t\t#in\n", + "b = 10. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "x = Sc*ratio*l/(Ss+Sc*ratio)\n", + "As = b*x*(x/2)/((l-x)*ratio)\n", + "Ina = (b*x**3/3)+ratio*As*(l-x)**2\n", + "M = Sc*I/x\n", + "N = As/(math.pi*(d/2)**2)\n", + "\t\t\t\n", + "# Results\n", + "print 'Number of steel bars required = %.2f'%(N)\n", + "print (\"it rounds to 6 bars\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Number of steel bars required = 5.79\n", + "it rounds to 6 bars\n" + ] + } + ], + "prompt_number": 11 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch13.ipynb b/Statics_And_Strength_Of_Materials/ch13.ipynb new file mode 100755 index 00000000..6f6d6663 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch13.ipynb @@ -0,0 +1,359 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:6b0b86bd9b0d0e04cc91f26a51f01e1dfca0b27db56af40f6330128871970df2" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 13 : Deflection of Beams" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 13.1 Page No : 313" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "E = 1.5*10**6\n", + "F1 = -100. \t\t\t#lb\n", + "F2 = -100. \t\t\t#lb\n", + "x1 = 6 \t\t\t#in\n", + "x2 = 6 \t\t\t#in\n", + "Ina = 64 \t\t\t#in**4\n", + "h1 = -600 \t\t\t#lb ft\n", + "h2 = -1200 \t\t\t#lb ft\n", + "xa1 = 10 \t\t\t# in\n", + "xa2 = 8 \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "deltamax = ((1./2)*x1*xa1*h1+(1./2)*(x1+x2)*h2*xa2)*(1728)/(E*Ina)\n", + "\t\t\t\n", + "# Results\n", + "print 'maximum deflection = %.2f in'%(deltamax)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "maximum deflection = -1.36 in\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 13.2 Page No : 314" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "E = 1.5*10**6\n", + "I = 50. \t \t\t#in**4\n", + "delta = -1. \t\t\t#in\n", + "l = 8. \t\t\t#ft\n", + "\t\t\t\n", + "# Calculations\n", + "w = -delta*8*E*I/(l**4*1728)\n", + "\t\t\t\n", + "# Results\n", + "print 'distributed weight = %.1f lb per ft'%(w)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "distributed weight = 84.8 lb per ft\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 13.3 Page No : 315" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "W = 50. \t\t\t#lb/ft\n", + "x = 5. \t\t\t#ft\n", + "x1 = 2. \t\t\t#ft\n", + "\t\t\t\n", + "# Calculations\n", + "V = W*x\n", + "M = W*((x/2)+x1)*x\n", + "M1 = W*x*(x+x1)\n", + "M2 = -M\n", + "M3 = -W*x*x/2\n", + "EIdeltamax = ((1./2)*(x+x1)*M1*((x+x1)/3))+(x+x1)*M2*((x+x1)/2)+(1./3)*x*M3*(x/4)\n", + "\t\t\t\n", + "# Results\n", + "print 'maximum value of EIdeltax = %.1f lb ft**3'%(round(EIdeltamax,-2))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "maximum value of EIdeltax = -14600.0 lb ft**3\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 13.5 Page No : 329" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "w = 180. \t\t\t#lb/ft\n", + "l = 8. \t\t\t#ft\n", + "P = 1200. \t\t\t#lb\n", + "b = 6. \t\t\t#ft\n", + "E = 3*10.**6\n", + "I = 64. \t\t\t#in**4\n", + "\t\t\t\n", + "# Calculations\n", + "delta = ((w*l**4)/(8))+((P*b**2)*(3*l-b)/(6))\n", + "\t\t\t\n", + "# Results\n", + "print 'deflection of the free end = %.1fbyEI ft'%(round(delta,-3))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "deflection of the free end = 222000.0byEI ft\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 13.6 Page No : 329" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "P = 6. \t\t\t#kips\n", + "w = 3. \t \t\t#kips/ft\n", + "L1 = 8. \t\t\t#ft\n", + "L2 = 8. \t\t\t#ft\n", + "\t\t\t\n", + "# Calculations\n", + "delta = (P*(L1+L2)**3/192)+(w*(L1+L2)**4/768)\n", + "\t\t\t\n", + "# Results\n", + "print 'midspan value of deflection = %.1f kip ft**3'%(delta)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "midspan value of deflection = 384.0 kip ft**3\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 13.7 Page No : 331" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "x1 = 3. \t\t\t#ft\n", + "x2 = 3. \t\t\t#ft\n", + "x3 = 3. \t\t\t#ft\n", + "x4 = 3. \t\t\t#ft\n", + "W1 = 4. \t\t\t#kips\n", + "W2 = 8. \t\t\t#kips\n", + "l = x1+x2+x3+x4\n", + "\t\t\t\n", + "# Calculations\n", + "b = x2+x3+x4\n", + "b1 = x4\n", + "a = x1\n", + "x = l/2\n", + "P = (((W1*b*(l/b*(x-a)**3+(l**2-b**2)*x-x**3))/(6*l))+((W2*b1*x*(l**2-x**2-b1**2))/(6*l)))*(48/l**3)\n", + "R1 = 3+2-(P/2)\n", + "R2 = P\n", + "R3 = 1+6-(P/2)\n", + "\t\t\t\n", + "# Results\n", + "print 'R1 = %.3f kips'%(R1)\n", + "print 'R2 =%.2f kips'%(R2)\n", + "print 'R3 =%.3f kips'%(R3)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "R1 = 0.875 kips\n", + "R2 =8.25 kips\n", + "R3 =2.875 kips\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 13.8 Page No : 333" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from numpy import linalg\n", + "\t\t\t\n", + "# Variables\n", + "P = 680. \t\t\t#lb\n", + "K = 1000. \t\t\t#lb/in\n", + "L = 6. \t\t\t#ft\n", + "E = 30.*10**6\n", + "Ina = 1.728 \t\t\t#in**4\n", + "\t\t\t\n", + "# Calculations\n", + "A = [[((L*12)**3/(3*E*Ina)),-(1/K)],[1,1]]\n", + "b = [0,P]\n", + "c = linalg.solve(A,b)\n", + "Pb = c[0]\n", + "Ps = c[1]\n", + "\t\t\t\n", + "# Results\n", + "print 'Force in the spring = %.2f psi'%(Ps)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Force in the spring = 480.00 psi\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 13.9 Page No : 334" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "I = 1.5 \t\t\t#in**4\n", + "Da = 0.5 \t\t\t#in\n", + "E = 30.*10**6\n", + "l = 60. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "F = 6*Da*E*I/(l**3)\n", + "\t\t\t\n", + "# Results\n", + "print 'F = %.2f lb'%(F)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "F = 625.00 lb\n" + ] + } + ], + "prompt_number": 8 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch14.ipynb b/Statics_And_Strength_Of_Materials/ch14.ipynb new file mode 100755 index 00000000..0d379010 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch14.ipynb @@ -0,0 +1,645 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:94d59c04fa299edfd91c8fbcd74f6a447f3853008eda8e5b17fed905ae2575e6" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 14 : Combined Loading" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.1 Page No : 349" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "h = 6. \t \t\t#in\n", + "x1 = 7. \t\t\t#in\n", + "x2 = 1. \t\t\t#in\n", + "x3 = 2. \t\t\t#in\n", + "P = 600. \t\t\t#lb\n", + "\t\t\t\n", + "# Calculations\n", + "By = P*(x1+x2+x3)/(x1+x2)\n", + "Bx = By*(x1+x2)/h\n", + "Fx = Bx\n", + "V = By-P\n", + "M = -P*(x2+x3)+By*x2\n", + "S1 = -Fx/(x3*h)\n", + "I = x3*h**3/12\n", + "S2 = -M*12*(h/2)/I\n", + "Scmax = S1-S2\n", + "Stmax = S1+S2\n", + "\t\t\t\n", + "# Results\n", + "print 'Maximum tensile stress at = %.1f psi'%(Scmax)\n", + "print 'Maximum compressive stress at = %.1f psi'%(Stmax)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum tensile stress at = -1133.3 psi\n", + "Maximum compressive stress at = 966.7 psi\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.2 Page No : 350" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "P = 10000. \t\t\t#lb\n", + "A = 11.77 \t\t\t#in**2\n", + "Z = 51.9 \t\t\t#in**3\n", + "x = 5. \t\t\t#ft\n", + "y = 12. \t\t\t#ft\n", + "\t\t\t\n", + "# Calculations\n", + "S1 = round(-P/A)\n", + "S2 = round(P*x*y/Z,-2)\n", + "Sc = S1-S2\n", + "St = S1+S2\n", + "\n", + "\t\t\t\n", + "# Results\n", + "print 'Axial stress at c = %.1f psi'%(Sc)\n", + "print 'Axial stress at t = %.1f psi'%(St)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Axial stress at c = -12450.0 psi\n", + "Axial stress at t = 10750.0 psi\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.3 Page No : 353" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "b = 6. \t\t \t#in\n", + "h = 12. \t\t\t#in\n", + "l = 20. \t\t\t#ft\n", + "P = 100000. \t\t#lb\n", + "\t\t\t\n", + "# Calculations\n", + "S = -P/(b*h)\n", + "S1 = l**2*6*12/(8*b*h**2)\n", + "w = -S/S1\n", + "\t\t\t\n", + "# Results\n", + "print 'Safe distributed load = %.0f lb per ft'%(w)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Safe distributed load = 333 lb per ft\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.4 Page No : 355" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "b = 4. \t\t\t#in\n", + "h = 9. \t\t\t#in\n", + "l = 6. \t\t\t#in\n", + "Mx = 600. \t\t\t#lb\n", + "My = 100. \t\t\t#lb\n", + "\t\t\t\n", + "# Calculations\n", + "Zx = b*h**3/(12*h/2)\n", + "Zy = b**3*h/(12*b/2)\n", + "S1 = Mx*l*12/Zx\n", + "S2 = My*b*12/Zy\n", + "Sb = S1+S2\n", + "Sd = -S1-S2\n", + "\t\t\t\n", + "# Results\n", + "print 'Maximum stress = %.1f psi tension)'%(Sb)\n", + "print ' Maximum stress =%.1f psi compression)'%(Sd) \n", + "\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum stress = 1000.0 psi tension)\n", + " Maximum stress =-1000.0 psi compression)\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.5 Page No : 356" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "d = 2. \t\t\t#in\n", + "Px = -600. \t\t\t#lb\n", + "Py = 1200. \t\t\t#lb\n", + "x1 = 2. \t\t\t#in\n", + "x2 = 2. \t\t\t#in\n", + "x3 = 2. \t\t\t#in\n", + "Ray = -400. \t\t\t#lb\n", + "Rax = 400. \t\t\t#lb\n", + "Rbx = 200. \t\t\t#lb\n", + "Rby = -800. \t\t\t#lb\n", + "\t\t\t\n", + "# Calculations\n", + "Mb = math.sqrt((Rax*x1)**2+(Ray*x1)**2)\n", + "Mc =math.sqrt((Rbx*x3)**2+(Rby*x3)**2)\n", + "if (Mb<Mc) :\n", + " M =Mc\n", + "else: \n", + " M = Mb \n", + "\n", + "Smax = M*12*64*(d/2)/(math.pi*d**4)\n", + "\t\t\t\n", + "# Results\n", + "print 'Maximum normal stress = %.1f psi'%(round(Smax,-2)) \n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum normal stress = 25200.0 psi\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.6 Page No : 359" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "P = 100. \t\t\t#kips\n", + "M = 400. \t\t\t#kip in\n", + "A = 14.7 \t\t\t#in**2\n", + "Z = 80.7 \t\t\t#in**3\n", + "\t\t\t\n", + "# Calculations\n", + "Smax = -(P*10**3)/A-(M*10**3)/Z\n", + "Smin = -(P*10**3)/A+(M*10**3)/Z\n", + "\t\t\t\n", + "# Results\n", + "print 'Maximum stress = %.1f psi'%(round(Smax,-2))\n", + "print 'Minimum stress =%.1f psi'%(round(Smin,-1))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum stress = -11800.0 psi\n", + "Minimum stress =-1850.0 psi\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.7 Page No : 360" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "As = 1. \t\t\t#in**2\n", + "Zs = 0.167 \t\t\t#in**3\n", + "Ah = 1. \t\t\t#in**2\n", + "Zh = 0.984 \t\t\t#in**3\n", + "es = 0.5 \t\t\t#in\n", + "eh = 0.5 \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "phbyps = (1/As+es/Zs)/(1/Ah+eh/Zh)\n", + "\t\t\t\n", + "# Results\n", + "print 'ratio = %.1f'%(phbyps)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "ratio = 2.6\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.8 Page No : 365" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Sx = 1800. \t\t\t#psi\n", + "Sy = 1000. \t\t\t#psi\n", + "angle = 30. \t\t\t#degrees\n", + "t = 0.25 \t\t\t#in\n", + "t1 = 3. \t\t\t#in\n", + "t2 = 5. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Sx1 = Sx/(t1*t)\n", + "Sy1 = Sy/(t2*t)\n", + "S = ((Sx+Sy)/2+((Sx-Sy)/2)*math.cos(math.radians(2*angle)))+(Sx-Sy)*math.cos(math.radians(2*angle))\n", + "T = (Sx-Sy)*math.sin(math.radians(2*angle))\n", + "\t\t\t\n", + "# Results\n", + "print 'S = %.1f psi'%(S)\n", + "print 'T =%.f psi'%(T)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "S = 2000.0 psi\n", + "T =693 psi\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.9 Page No : 365" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "Sx = 1800. \t\t\t#lb\n", + "Sy = 1000. \t\t\t#lb\n", + "angle = 30. \t\t\t#degrees\n", + "\t\t\t\n", + "# Calculations\n", + "Sa =-((Sx+Sy)/2+((Sx-Sy)/2)*math.cos(math.radians(2*angle)))*math.cos(math.radians(2*angle))-(Sx-Sy)\n", + "Ta = -((Sx+Sy)/2+((Sx-Sy)/2)*math.cos(math.radians(2*angle)))*math.sin(math.radians(2*angle))\n", + "Sb =((Sx+Sy)/2+((Sx-Sy)/2)*math.cos(math.radians(2*angle)))*math.cos(math.radians(2*angle))-(Sx-Sy)\n", + "Tb = ((Sx+Sy)/2+((Sx-Sy)/2)*math.cos(math.radians(2*angle)))*math.sin(math.radians(2*angle))\n", + "\t\t\t\n", + "# Results\n", + "print 'Sa = %.1f psi'%(Sa)\n", + "print 'Sb =%.1f psi'%(Sb)\n", + "print 'Ta =%.f psi'%(round(Ta,-1))\n", + "print 'Tb =%.1f psi'%(round(Tb,-1))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Sa = -1600.0 psi\n", + "Sb =0.0 psi\n", + "Ta =-1390 psi\n", + "Tb =1390.0 psi\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.10 Page No : 369" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "\t\t\n", + "# Variables\n", + "angle = 15. \t\t\t#degrees\n", + "Tyx = -1000. \t\t\t#psi\n", + "Txy = 1000. \t\t\t#psi\n", + "\t\t\t\n", + "# Calculations\n", + "Sx = Txy*math.sin(math.radians(2*angle))\n", + "Tx = Txy*math.cos(math.radians(2*angle))\n", + "Sy = Tyx*math.sin(math.radians(2*angle))\n", + "Ty = Tyx*math.cos(math.radians(2*angle))\n", + "Sx1 = Txy\n", + "Sy1 = Tyx\n", + "Txy = 0\n", + "\t\t\t\n", + "# Results\n", + "print 'Sx = %.1f psi'%(Sx)\n", + "print 'Tx = %.1f psi'%(Tx)\n", + "print 'Sy = %.1f psi'%(Sy)\n", + "print 'Ty = %.1f psi'%(Ty)\n", + "print 'Sx1 = %.1f psi'%(Sx1)\n", + "print 'Sy1 = %.1f psi'%(Sy1)\n", + "print 'Txy = %.1f psi'%(Txy)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Sx = 500.0 psi\n", + "Tx = 866.0 psi\n", + "Sy = -500.0 psi\n", + "Ty = -866.0 psi\n", + "Sx1 = 1000.0 psi\n", + "Sy1 = -1000.0 psi\n", + "Txy = 0.0 psi\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.11 Page No : 370" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "d = 4. \t\t\t#in\n", + "n = 315. \t\t\t#rpm\n", + "Ss = 8000. \t\t\t#psi\n", + "Ns = 12000. \t\t\t#psi\n", + "\t\t\t\n", + "# Calculations\n", + "T = Ss*d**4/(32*(d/2))\n", + "hp = T*math.pi*n/63000\n", + "\t\t\t\n", + "# Results\n", + "print 'T = %.1f pi lb in'%(T)\n", + "print 'horsepower rating =%.1f hp'%(round(hp,-1))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "T = 32000.0 pi lb in\n", + "horsepower rating =500.0 hp\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.12 Page No : 372" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Sx = 9. \t\t\t#ksi\n", + "Sy = -5. \t\t\t#ksi\n", + "Txy = 4. \t\t\t#ksi\n", + "\t\t\t\n", + "# Calculations\n", + "R = math.sqrt(((Sx-Sy)/2)**2+Txy**2)\n", + "Smax = ((Sx+Sy)/2)+R\n", + "Smin = ((Sx+Sy)/2)-R\n", + "ap1 = (1./2)*math.degrees(math.atan(2*Txy/(Sx-Sy)))\n", + "ap2 = 90+ap1\n", + "Sc = (Sx+Sy)/2\n", + "Tc = R\n", + "Sd = (Sx+Sy)/2\n", + "Td = -R\n", + "a1 = (90-2*ap1)/2\n", + "a2 = 90+a1\n", + "\t\t\t\n", + "# Results\n", + "print 'Smax = %.2f ksi'%(Smax)\n", + "print 'Smin = %.2f ksi'%(Smin)\n", + "print 'R = %.2f psi'%(R)\n", + "print 'palne1 = %.2f degrees'%(ap1)\n", + "print 'plane 2 = %.2f degrees'%(ap2)\n", + "print 'Sc = %.2f ksi'%(Sc)\n", + "print 'Sd = %.2f ksi'%(Sd)\n", + "print 'Tc = %.2f ksi'%(Tc)\n", + "print 'Td = %.2f ksi'%(Td)\n", + "print 'palne1 = %.2f degrees'%(a1)\n", + "print 'plane 2 = %.2f degrees'%(a2)\n", + "\n", + "# note : rounding off error." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Smax = 10.06 ksi\n", + "Smin = -6.06 ksi\n", + "R = 8.06 psi\n", + "palne1 = 14.87 degrees\n", + "plane 2 = 104.87 degrees\n", + "Sc = 2.00 ksi\n", + "Sd = 2.00 ksi\n", + "Tc = 8.06 ksi\n", + "Td = -8.06 ksi\n", + "palne1 = 30.13 degrees\n", + "plane 2 = 120.13 degrees\n" + ] + } + ], + "prompt_number": 16 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 14.13 Page No : 373" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "d = 4. \t\t \t#in\n", + "T = 40000. \t\t\t #lb in\n", + "Th = 20000. \t\t\t#lb in\n", + "\t\t\t\n", + "# Calculations\n", + "t = round(T*(d/2)*32/(math.pi*d**4),-1)\n", + "S = round(Th/(math.pi*(d/2)**2),-1)\n", + "Smax = -(S/2)-math.sqrt(t**2+(S/2)**2)\n", + "Tmax = math.sqrt(t**2+(S/2)**2)\n", + "\n", + "\t\t\t\n", + "# Results\n", + "print 'Maximum normal stress = %.f psi'%(Smax)\n", + "print 'Maximum shearing stress =%.f psi'%(Tmax)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum normal stress = -4073 psi\n", + "Maximum shearing stress =3278 psi\n" + ] + } + ], + "prompt_number": 21 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch15.ipynb b/Statics_And_Strength_Of_Materials/ch15.ipynb new file mode 100755 index 00000000..562ecba1 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch15.ipynb @@ -0,0 +1,361 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:fe60a567abd9a722f9efaa797b5de23eaf91ae1c1b70917990ddaa017e213706" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 15 : Welded, Bolted, and Riveted Connections" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 15.1 Page No : 393" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "sigma = 20000. \t\t\t#psi\n", + "b = 6. \t\t \t#in\n", + "h = 0.5 \t\t \t#in\n", + "p1 = 3750.\n", + "\t\t\t\n", + "# Calculations\n", + "P = sigma*b*h\n", + "L = (P-p1*b)/(2*p1)\n", + "\t\t\t\n", + "# Results\n", + "print 'L = %.2f in'%(L)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "L = 5.00 in\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 15.2 Page No : 294" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from numpy import linalg\n", + "\t\t\t\n", + "# Variables\n", + "P = 5000. \t\t\t#lb per in\n", + "Tl = 75. \t\t\t#kips\n", + "y1 = 2.63 \t\t\t#in\n", + "y2 = 1.37 \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "A = [[P, P],[y1*P, -y2*P]]\n", + "b = [Tl*10**3, 0]\n", + "c = linalg.solve(A,b)\n", + "L1 = c[0]\n", + "L2 = c[1]\n", + "\t\t\t\n", + "# Results\n", + "print 'L1 = %.2f in'%(L1)\n", + "print 'L2 = %.2f in'%(L2)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "L1 = 5.14 in\n", + "L2 = 9.86 in\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 15.3 Page No : 397" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "d = 3./8 \t \t\t#in\n", + "d1 = 1./8 \t\t \t#in\n", + "y = 1. \t\t \t #in\n", + "T = 15000. \t\t\t #psi\n", + "sigmab = 32000. \t\t\t#psi\n", + "sigmat = 18000. \t\t\t#psi\n", + "\t\t\t\n", + "# Calculations\n", + "Ps = math.pi*T*(d/2)**2\n", + "Pt = sigmat*d1*(y-d)\n", + "Pb = sigmab*d1*d\n", + "Pmin =Ps\n", + "sigma =T\n", + "if(Pt<Pmin):\n", + " Pmin =Pt\n", + " sigma =sigmat\n", + "else:\n", + " Pmin =Pb\n", + " sigma =sigmab\n", + "\n", + "e = Pmin*100/(sigma*d1*y)\n", + "\t\t\t\n", + "# Results\n", + "print 'e = %.2f per cent'%(e)\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "e = 62.50 per cent\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 15.4 Page No : 398" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "d = 7./8 \t\t\t#in\n", + "Ss = 15000. \t\t\t#psi\n", + "Sb = 32000. \t\t\t#psi\n", + "St = 20000. \t\t\t#psi\n", + "n = 8.\n", + "t = 3./8 \t\t\t#in\n", + "l = 10. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Ps = Ss*math.pi*n*(d/2)**2\n", + "Pb = Sb*math.pi*n*d*t\n", + "Pt1 = St*(l-d*2)*t\n", + "Pt2 = 4*St*(l-d*4)*t/3\n", + "Pt3 = 4*St*(l-d*2)*t\n", + "Pmin = Ps\n", + "sigma = Ss\n", + "if (Pb<Pmin):\n", + " Pmin =Pb\n", + " sigma =Sb\n", + "elif (Pt1<Pmin):\n", + " Pmin =Pt1\n", + " sigma =St\n", + "elif (Pt2<Pmin):\n", + " Pmin =Pt2\n", + " sigma =St\n", + "elif (Pt3<Pmin):\n", + " Pmin =Pt3\n", + " sigma =St\n", + "e = Pmin*100/(sigma*t*l)\n", + "\t\t\t\n", + "# Results\n", + "print 'e = %.1f per cent'%(e)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "e = 82.5 per cent\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 15.5 Page No : 400" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "n = 8.\n", + "shear = 15. \t\t\t#ksi\n", + "Dr = 7/8. \t\t\t#in\n", + "Ss = 32. \t\t\t#ksi\n", + "Ds = 40. \t\t\t#si\n", + "D = 3/8. \t\t\t#in\n", + "x = 0.504 \t\t\t#in\n", + "pmin=0\n", + "\t\t\t\n", + "# Calculations\n", + "Ps = shear*n*(Dr/2)**2\n", + "Pb = Ds*(n/2)*x*Dr\n", + "Pb1 = Ss*n*D*Dr\n", + "pmin = Ps\n", + "if (Pb<pmin):\n", + " Pmin = Pb\n", + "else: \n", + " Pmin = Pb1 \n", + "\t\t\t\n", + "# Results\n", + "print 'load capacity of connection = %.1f kips'%(Pb)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "load capacity of connection = 70.6 kips\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 15.6 Page No : 401" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "T = 15000. \t\t\t#psi\n", + "x1 = 3. \t\t\t#in\n", + "x2 = 3. \t\t\t#in\n", + "y1 = 3. \t\t\t#in\n", + "y2 = 3. \t\t\t#in\n", + "d = 0.5 \t\t\t#in\n", + "n = 4.\n", + "\t\t\t\n", + "# Calculations\n", + "P = T*(math.pi/4)*d**2/(math.sqrt((1/n)**2+(1/((math.sqrt(y1**2+y2**2)/y1)*n))**2+ \\\n", + "(2*(1/n)*(1/(n*(math.sqrt(y1**2+y2**2))/y1))*math.cos(math.radians(45)))))\n", + "P1 = T*(math.pi/4)*d**2/((1/n)+(y1/(n*y1)))\n", + "if (P>P1):\n", + " print 'Stornger P = %.2f lb'%(P)\n", + "else:\n", + " print 'Stornger P = %.2f lb'%(P1)\n", + "\n", + "# note : answer is different because of rounding off error. please check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Stornger P = 7450.94 lb\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 15.7 Page No : 403" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "P = 5. \t\t\t#kips\n", + "xab = 3. \t\t\t#in\n", + "xbc = 6. \t\t\t#in\n", + "xbp = 1. \t\t\t#in\n", + "y = 6. \t\t\t#in\n", + "n = 3.\n", + "\t\t\t\n", + "# Calculations\n", + "Dl = P/3\n", + "Pct = (6*P)/(((xab+xbp)*(xab+xbp)/(xbc-xbp))+(xbp/(xbc-xbp))+(xbc-xbp))\n", + "R = math.sqrt(Pct**2+Dl**2)\n", + "\t\t\t\n", + "# Results\n", + "print 'Greatest Load = %.2f kips'%(R)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Greatest Load = 3.94 kips\n" + ] + } + ], + "prompt_number": 7 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch16.ipynb b/Statics_And_Strength_Of_Materials/ch16.ipynb new file mode 100755 index 00000000..ca14e404 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch16.ipynb @@ -0,0 +1,387 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:d708156611ec89eae1db1c06f5ba02ec9d179b4f495d97ecac600c0bd4aa5921" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 16 : Columns" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 16.1 Page No : 417" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "E = 10.*10**6 \t\t\t#psi\n", + "ys = 6000. \t\t\t#psi\n", + "\t\t\t\n", + "# Calculations\n", + "lbyr = math.sqrt(math.pi**2*E/ys)\n", + "\t\t\t\n", + "# Results\n", + "print 'Slenderness Ratio = %.f '%(lbyr)\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slenderness Ratio = 128 \n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 16.2 Page No : 417" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "fs = 3.\n", + "W = 50. \t\t\t#kips\n", + "l = 20. \t\t\t#ft\n", + "E = 30.*10**6 \t\t\t#psi\n", + "\t\t\t\n", + "# Calculations\n", + "Pcr = fs*W\n", + "I = Pcr*10**3*(l*12)**2/(math.pi**2*E)\n", + "r = 2.01\n", + "lbyr = l*12/r\n", + "\t\t\t\n", + "# Results\n", + "print 'Required I = %.1f in**4'%(I)\n", + "print 'slenderness ratio =%.f '%(lbyr)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Required I = 29.2 in**4\n", + "slenderness ratio =119 \n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 16.3 Page No : 418" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "\t\t\t\n", + "# Variables\n", + "L1 = 18. \t\t\t#ft\n", + "L2 = 9. \t\t\t#ft\n", + "I1 = 12.1 \t\t\t#in**4\n", + "I2 = 1.2 \t\t\t#in**4\n", + "E = 30.*10**6 \t\t\t#psi\n", + "\t\t\t\n", + "# Calculations\n", + "r1 = 2.05\n", + "lbyr = L1*12/r1\n", + "r2 = 0.65\n", + "lbyr2 = L2*12/r2\n", + "Pcr1 = math.pi**2*E*I1/(L1*12)**2\n", + "Pcr2 = math.pi**2*E*I2/(L2*12)**2\n", + "P = Pcr1/2.5\n", + "P2 = Pcr2/2.5\n", + "\t\t\t\n", + "# Results\n", + "print 'Design load of 1 = %.2f lb'%(round(P,-2))\n", + "print 'Design load of 2 =%.2f lb'%(round(P2,-2))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Design load of 1 = 30700.00 lb\n", + "Design load of 2 =12200.00 lb\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 16.4 Page No : 419" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "E = 30.*10**6\n", + "syp = 30000. \t\t\t#psi\n", + "I = 143.5 \t\t\t#in**4\n", + "A = 7.32 \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "I1 = 2*I\n", + "A1 = 2*A\n", + "L = math.sqrt(2*math.pi**2*E*I1/(syp*A1))\n", + "\t\t\t\n", + "# Results\n", + "print 'Critical length of the column = %.0f in'%(L)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Critical length of the column = 622 in\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 16.5 Page No : 421" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "x = 30. \t\t\t#in\n", + "x1 = 10. \t\t\t#in\n", + "E = 30.*10**6\n", + "d = 0.5 \t\t\t#in\n", + "syp = 60000. \t\t\t#psi\n", + "y1 = 8. \t\t\t#in\n", + "y2 = 2. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "ratio = 0.8\n", + "l = x+x1\n", + "lr = ratio*l\n", + "I = (math.pi*(d)**4)/64\n", + "Pcr = math.pi**2*E*I/lr**2\n", + "scr = Pcr/(math.pi*(d/2)**2)\n", + "F = Pcr*y2/(y1+y2)\n", + "\t\t\t\n", + "# Results\n", + "print 'Stress in the critical load = %.0f psi'%(round(scr,-1))\n", + "print 'Critical force F =%.0f lb'%(F)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Stress in the critical load = 4520 psi\n", + "Critical force F =177 lb\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 16.6 Page No : 423" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "l = 10. \t\t\t#ft\n", + "Ys = 33000. \t\t\t#psi\n", + "E = 30.*10**6\n", + "A = 13.24 \t\t\t#in**4\n", + "\t\t\t\n", + "# Calculations\n", + "r = 2.\n", + "lbyr = l*12/r\n", + "Cc = math.sqrt(2*math.pi**2*E/Ys)\n", + "fs = 5./3+3*(lbyr)/(8*Cc)+(lbyr)**3/(5*Cc**3)\n", + "Sa =((1-((lbyr)**2/(2*Cc**2)))*(Ys))/fs\n", + "Pa = Sa*A\n", + "\n", + "\t\t\t\n", + "# Results\n", + "print 'Premissible load = %.f kips'%(Pa/1000)\n", + "\n", + "# note : rounding off error." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Premissible load = 212 kips\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 16.7 Page No : 425" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "from sympy import Symbol,solve\n", + "\t\t\t\n", + "# Variables\n", + "L = 12. \t\t\t#ft\n", + "Po = 100. \t\t\t#kips\n", + "e = 2. \t\t\t#ft\n", + "ys = 42000. \t\t\t#psi\n", + "A = 11.77 \t\t\t#in**2\n", + "rmin = 195. \t\t\t#in\n", + "Zmin = 11.0 \t\t\t#in**3\n", + "lbyr = 74.2\n", + "stress = 18. \t\t\t#ksi\n", + "\t\t\t\n", + "# Calculations\n", + "#P = (stress-(Po/A)/((1/A)+((e*12)/Zmin)))\n", + "P = Symbol(\"P\")\n", + "eq = (((Po+P)/A) + (P*(e*L)/11))/stress - 1\n", + "P = solve(eq,P)[0]\n", + "\t\t\t\n", + "# Results\n", + "print 'Additional Load = %.2f kips '%(P)\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Additional Load = 4.19 kips \n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 16.8 Page No : 425" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "l = 15. \t\t\t#ft\n", + "Al = 80. \t\t\t#kips\n", + "El = 60. \t\t\t#kips\n", + "Ys = 33. \t\t\t#ksi\n", + "e = 4. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "A = 14.4\n", + "rmin = 2.54\n", + "Zxx = 54.6\n", + "lbyr = l*12/rmin\n", + "Smax = ((Al+El)/A)+El*e/Zxx\n", + "\t\t\t\n", + "# Results\n", + "print 'Maximum stress %.1f ksi'%(Smax)\n", + "print (\"10 WF 49 is the suitable one\")\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum stress 14.1 ksi\n", + "10 WF 49 is the suitable one\n" + ] + } + ], + "prompt_number": 12 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch2.ipynb b/Statics_And_Strength_Of_Materials/ch2.ipynb new file mode 100755 index 00000000..74bc7026 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch2.ipynb @@ -0,0 +1,358 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:ddc392c3474d618cbae4adde9fb1f47f093e50fd2a55d0b1580decec813c6214" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 2 : Force System : Components, Resultants, Equivalence" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.3 Page No : 21" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "f1 = 20. \t\t\t#lb\n", + "f2 = 40. \t\t\t#lb\n", + "alpha = 30. \t\t\t#degrees\n", + "\t\t\t\n", + "# Calculations\n", + "R = math.sqrt(f1**2+f2**2+2*f1*f2*math.cos(math.radians(alpha)))\n", + "angle = math.degrees(math.asin((f2*math.sin(math.radians(180-alpha)))/(R)))\n", + "\t\t\t\n", + "# Results\n", + "print 'R = %.1f lb'%(R)\n", + "print 'angle = %.1f degrees'%(angle)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "R = 58.2 lb\n", + "angle = 20.1 degrees\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.4 Page No : 22" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "fx = 100. \t\t\t#lb\n", + "f1 = 200. \t\t\t#lb\n", + "f2 = 100. \t\t\t#lb\n", + "f3 = 50. \t\t\t#lb\n", + "a1 = 30. \t\t\t#degrees\n", + "a2 = 45. \t\t\t#degrees\n", + "a3 = 60. \t\t\t#degrees\n", + "\t\t\t\n", + "# Calculations\n", + "Rx = fx+f1*math.cos(math.radians(a1))-f2*math.cos(math.radians(a2))-f3*math.cos(math.radians(a3))\n", + "Ry = f1*math.sin(math.radians(a1))+f2*math.sin(math.radians(a2))-f3*math.sin(math.radians(a3))\n", + "R = math.sqrt(Rx**2+Ry**2)\n", + "angle = math.degrees(math.atan(Ry/Rx))\n", + "\t\t\t\n", + "# Results\n", + "print 'R = %.f lb'%(R)\n", + "print 'angle = %.1f degrees'%(angle)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "R = 218 lb\n", + "angle = 35.7 degrees\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.5 Page No : 25" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "f1 = 100. \t\t\t#lb\n", + "f2 = 200. \t\t\t#lb\n", + "x1 = 2.\n", + "x2 = -3.\n", + "y1 = 3.\n", + "y2 = 5.\n", + "z1 = 4.\n", + "z2 = -2.\n", + "\t\t\t\n", + "# Calculations\n", + "d1 = math.sqrt(x1**2+y1**2+z1**2)\n", + "d2 = math.sqrt(x2**2+y2**2+z2**2)\n", + "f1x = f1*x1/d1\n", + "f1y = f1*y1/d1\n", + "f1z = f1*z1/d1\n", + "f2x = f2*x2/d2\n", + "f2y = f2*y2/d2\n", + "f2z = f2*z2/d2\n", + "Rx = f1x+f2x\n", + "Ry = f1y+f2y\n", + "Rz = f1z+f2z\n", + "R = math.sqrt(Rx**2+Ry**2+Rz**3)\n", + "I1 = Rx/R\n", + "I2 = Ry/R\n", + "I3 = Rz/R\n", + "\t\t\t\n", + "# Results\n", + "print 'R = %d lb'%(R)\n", + "print 'I1 = %.3f '%(I1)\n", + "print 'I2 = %.3f '%(I2)\n", + "print 'I3 = %.3f '%(I3)\n", + "\n", + "# note : rounding off error would be there for R." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "R = 227 lb\n", + "I1 = -0.264 \n", + "I2 = 0.956 \n", + "I3 = 0.041 \n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.6 Page No : 27" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "F = 100. \t\t\t#lb\n", + "x1 = 6. \t\t\t#in\n", + "x2 = 8. \t\t\t#in\n", + "x3 = 2. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "xab = math.sqrt(x1**2+x2**2)\n", + "d = x3*x1/xab\n", + "M1 = F*d\n", + "Fx = F*x2/xab\n", + "Fy = F*x1/xab\n", + "M2 = Fy*xab-Fx*x1\n", + "M3 = Fy*x3\n", + "\t\t\t\n", + "# Results\n", + "print 'M1 = %.f lb.in'%(M1)\n", + "print 'M2 = %.f lb.in'%(M2)\n", + "print 'M3 = %.f lb.in'%(M3)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "M1 = 120 lb.in\n", + "M2 = 120 lb.in\n", + "M3 = 120 lb.in\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.7 Page No : 30" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Fy1 = 2. \t\t\t#kips\n", + "Fy2 = 5. \t\t\t#kips\n", + "Fy3 = 10. \t\t\t#kips\n", + "Fy4 = 3. \t\t\t#kips\n", + "L = 5. \t\t\t#ft\n", + "\n", + "# Calculations\n", + "Ry = Fy1+Fy2+Fy3+Fy4\n", + "x = (Fy1*L+Fy2*2*L+Fy3*3*L+Fy4*4*L)/Ry\n", + "\t\t\t\n", + "# Results\n", + "print 'Ry = %.2f kips'%(Ry)\n", + "print 'x = %.1f ft to the right of O'%(x)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ry = 20.00 kips\n", + "x = 13.5 ft to the right of O\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.8 Page No : 30" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \t\t\t\n", + "\n", + "# Variables\n", + "Fx1 = -15. \t\t\t#lb\n", + "Fx2 = 55. \t\t\t#lb\n", + "Fy1 = 70. \t\t\t#lb\n", + "Fy2 = -40. \t\t\t#lb\n", + "x1 = 4. \t\t\t#in\n", + "x2 = 3. \t\t\t#in\n", + "x3 = 5. \t\t\t#in\n", + "y1 = 4. \t\t\t#in\n", + "y2 = 2. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Rx = Fx1+Fx2\n", + "Ry = Fy1+Fy2\n", + "R = math.sqrt(Rx**2+Ry**2)\n", + "angle = math.degrees(math.atan(Ry/Rx))\n", + "\t\t\t\n", + "# Results\n", + "print 'R = %.2f lb'%(R)\n", + "print 'angle = %.1f degrees'%(angle)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "R = 50.00 lb\n", + "angle = 36.9 degrees\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.9 Page No : 32" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Fy = 200. \t\t\t#lb\n", + "Fx = 100. \t\t\t#lb\n", + "y = 3. \t \t\t#in\n", + "x = 6. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "M = Fy*x-Fx*y\n", + "\t\t\t\n", + "# Results\n", + "print 'Moment = %.2f lb in'%(M)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Moment = 900.00 lb in\n" + ] + } + ], + "prompt_number": 8 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch3.ipynb b/Statics_And_Strength_Of_Materials/ch3.ipynb new file mode 100755 index 00000000..b1bd4e75 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch3.ipynb @@ -0,0 +1,255 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:2df55062ff2333a87c8624b14338bf6cf25e04a49f7bd6f183057766f3ec7b9a" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 3 : Center of Gravity" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.1 Page No : 50" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "W = 3000. \t\t\t#lb\n", + "L = 10. \t\t\t#ft\n", + "Wf1 = 1200. \t\t\t#lb\n", + "Wf2 = 1500. \t\t\t#lb\n", + "angle = 30. \t\t\t#degrees\n", + "\t\t\t\n", + "# Calculations\n", + "d1 = Wf1*math.cos(angle)*L/W\n", + "d2 = Wf2*L/W\n", + "xbc = d1/math.cos(angle)\n", + "xab = d2-xbc\n", + "y = xab/math.tan(math.radians(angle))\n", + "\t\t\t\n", + "# Results\n", + "print 'x = %.2f ft'%(d2)\n", + "print 'y = %.2f ft'%(y)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "x = 5.00 ft\n", + "y = 1.73 ft\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.2 Page No : 51" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "W4 = 3. \t\t\t#lb\n", + "W3 = 5. \t\t\t#lb\n", + "W2 = 2. \t\t\t#lb\n", + "W1 = 6. \t\t\t#lb\n", + "x1 = 10. \t\t\t#in\n", + "x2 = 4. \t\t\t#in\n", + "z = 5. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "W = W1+W2+W3+W4\n", + "x = (W1*0+W2*0+W3*x2+W4*x1)/W\n", + "z = (W1*z+W2*0+W3*0+W4*0)/W\n", + "\t\t\t\n", + "# Results\n", + "print 'x = %.2f in'%(x)\n", + "print 'z = %.2f in'%(z)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "x = 3.12 in\n", + "z = 1.88 in\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.3 Page No : 52" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "import math \n", + "# Variables\n", + "W1 = 3. \t\t\t#lb\n", + "W2 = 5. \t\t\t#lb\n", + "x1 = 8. \t\t\t#in\n", + "x2 = 7. \t\t\t#in\n", + "y1 = 2. \t\t\t#in\n", + "y2 = 5. \t\t\t#in\n", + "z1 = 6. \t\t\t#in\n", + "z2 = 4. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "W = W1+W2\n", + "x = (W1*x1+W2*x2)/W\n", + "y = (W1*y1+W2*y2)/W\n", + "z = (W1*z1+W2*z2)/W\n", + "\t\t\t\n", + "# Results\n", + "print 'x = %.2f in'%(x)\n", + "print 'y = %.2f in'%(y)\n", + "print 'z = %.2f in'%(z)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "x = 7.38 in\n", + "y = 3.88 in\n", + "z = 4.75 in\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.4 Page No : 54" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "L = 9. \t\t\t#in\n", + "B = 16. \t\t\t#in\n", + "B1 = 6. \t\t\t#in\n", + "d = 2. \t \t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "x = ((L*(B-B1)*(L/2)+(1./2)*L*B1*(L/3)-(math.pi/4)*d**2*(L/2)))/(L*(B-B1)+(1./2)*L*B1-(math.pi/4)*d**2)\n", + "y = ((L*(B-B1)*((B-B1)/2)+(1./2)*L*B1*(B1/3+(B-B1))-(math.pi/4)*d**2*((B-B1)/2)))/(L*(B-B1)+(1./2)*L*B1-(math.pi/4)*d**2)\n", + "\t\t\t\n", + "# Results\n", + "print 'x = %.2f in to the right of y-axis'%(x)\n", + "print 'y = %.2f in above x axis'%(y)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "x = 4.14 in to the right of y-axis\n", + "y = 6.66 in above x axis\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.5 Page No : 57" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "Gt = 0.25 \t\t\t#in\n", + "St = 0.25 \t\t\t#in\n", + "Gw = 3.5 \t\t\t#lb/sq ft\n", + "Sw = 10. \t\t\t#lb/sq ft\n", + "Sb = 36. \t\t\t#in\n", + "Sb1 = 18. \t\t\t#in\n", + "Sb2 = 12. \t\t\t#in\n", + "Sb3 = 6. \t\t\t#in\n", + "Sy1 = 6. \t\t\t#in\n", + "Sy2 = 12. \t\t\t#in\n", + "Sy3 = 6. \t\t\t#in\n", + "Gb = 1. \t\t\t#ft\n", + "Sh = 24. \t\t\t#in\n", + "Gh = 1. \t\t\t#ft\n", + "\t\t\t\n", + "# Calculations\n", + "W = ((Sb*Sh)/(12*12)-(Gh*Gb))*Sw+(Gh*Gb)*Gw\n", + "x = ((Sb*Sh)*Sw*(Sb/24)/(12*12)-(Gh*Gb)*Sw*((Sb1+(Sb2/2))/12)+(Gh*Gb)*Gw*((Sb1+(Sb2/2))/12))/W\n", + "\t\t\t\n", + "# Results\n", + "print 'centre of gravity = %.2f ft to the right of y-axis'%(x)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "centre of gravity = 1.44 ft to the right of y-axis\n" + ] + } + ], + "prompt_number": 6 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch4.ipynb b/Statics_And_Strength_Of_Materials/ch4.ipynb new file mode 100755 index 00000000..908d7173 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch4.ipynb @@ -0,0 +1,257 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:5a553b2a442e6221dea12996ec8dc7d308a698d61fa5538169e60bbcf11198bb" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 4 : Equilibrium" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.3 Page No : 71" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "from numpy import linalg\t\t\n", + " \t\n", + "# Variables\n", + "W = 100. \t\t\t#lb\n", + "a1 = 30. \t\t\t#degrees\n", + "a2 = 45. \t\t\t#degrees\n", + "\t\t\t\n", + "# Calculations\n", + "A =[[math.cos(math.radians(a2)),-math.cos(math.radians(a1))],[math.sin(math.radians(a2)),math.sin(math.radians(a1))]]\n", + "b =[[0],[W]]\n", + "c = linalg.solve(A,b)\n", + "Tbc = c[0]\n", + "Tab = c[1]\n", + "\t\t\t\n", + "# Results\n", + "print 'Tbc = %.1f lb'%(Tbc)\n", + "print 'Tab =%.1f lb'%(Tab)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Tbc = 89.7 lb\n", + "Tab =73.2 lb\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.4 Page No : 72" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "W1 = 7000. \t\t\t#lb\n", + "W2 = 1000. \t\t\t#lb\n", + "W3 = 3000. \t\t\t#lb\n", + "x1 = 6. \t\t\t#in\n", + "x2 = 9. \t\t\t#in\n", + "x3 = 10. \t\t\t#in\n", + "x4 = 5. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Rb = (W1*x1+W2*(x1+x2)+W3*(x1+x2+x3))/(x1+x2+x3+x4)\n", + "Ra = W1+W2+W3-Rb\n", + "\t\t\t\n", + "# Results\n", + "print 'Rb = %.1f lb'%(Rb)\n", + "print 'Ra = %.1f lb'%(Ra)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rb = 4400.0 lb\n", + "Ra = 6600.0 lb\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.5 Page No : 74" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Fc = 500. \t\t\t#lb\n", + "Fd = 1000. \t\t\t#lb\n", + "xc = 2. \t\t\t#in\n", + "xd = 8. \t\t\t#in\n", + "y = 6. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Ay = Fc+Fd\n", + "Bx = (Fc*xc+Fd*xd)/y\n", + "Ax = Bx\n", + "A = math.sqrt(Ax**2+Ay**2)\n", + "\t\t\t\n", + "# Results\n", + "print 'A = %.f lb'%(A)\n", + "print 'B = %.f lb'%(Bx)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "A = 2121 lb\n", + "B = 1500 lb\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.6 Page No : 75" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "W = -300. \t\t\t#lb\n", + "r = 4. \t\t\t#in\n", + "x1 = 2. \t\t\t#ft\n", + "x2 = 3. \t\t\t#ft\n", + "x3 = 1. \t\t\t#ft\n", + "y1 = 1. \t\t\t#ft\n", + "x4 = 3. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "F = -W*r/(y1*12)\n", + "By = -W*x1/(x1+x2)\n", + "Bz = -F*(x1+x2+x3+(x4/12))/(x1+x2)\n", + "Ay = -W-By\n", + "Az = -F-Bz\n", + "\t\t\t\n", + "# Results\n", + "print 'Ay = %.2f lb'%(Ay)\n", + "print 'By = %.2f lb'%(By)\n", + "print 'Az = %.2f lb'%(Az)\n", + "print 'Bz = %.2f lb'%(Bz)\n", + "print 'F = %.2f lb'%(F)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ay = 180.00 lb\n", + "By = 120.00 lb\n", + "Az = 25.00 lb\n", + "Bz = -125.00 lb\n", + "F = 100.00 lb\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.7 Page No : 77" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "W = 500. \t\t\t#lb\n", + "r = 4. \t\t\t#in\n", + "Lx = 3. \t\t\t#in\n", + "Ly = 12. \t\t\t#in\n", + "Lz = 4. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Tbd = W*(math.sqrt((-Lx)**2+(-Ly)**2+(-Lz)**2))/Ly\n", + "Tcd = Lx*Tbd/(math.sqrt((-Lx)**2+(-Ly)**2+(-Lz)**2))\n", + "Tad = Lz*Tbd/(math.sqrt((-Lx)**2+(-Ly)**2+(-Lz)**2))\n", + "\t\t\t\n", + "# Results\n", + "print 'Tbd = %.f lb'%(Tbd)\n", + "print 'Tcd =%.f lb'%(Tcd)\n", + "print 'Tad =%.f lb'%(Tad)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Tbd = 542 lb\n", + "Tcd =125 lb\n", + "Tad =167 lb\n" + ] + } + ], + "prompt_number": 5 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch5.ipynb b/Statics_And_Strength_Of_Materials/ch5.ipynb new file mode 100755 index 00000000..997d11ff --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch5.ipynb @@ -0,0 +1,148 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:99441c9b87820b036679d4541d0d259ff9dba2feea52f8d659ae363ae4db1049" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 5 : Force Analysis of Structures" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.1 Page No : 98" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "import math \n", + "# Variables\n", + "Fc = -1000. \t\t\t#lb\n", + "A = 60. \t\t\t#degrees\n", + "E1 = 60. \t\t\t#degrees\n", + "E2 = 60. \t\t\t#degrees\n", + "D = 60. \t\t\t#degrees\n", + "L1 = 10. \t\t\t#ft\n", + "L2 = 10. \t\t\t#ft\n", + "\t\t\t\n", + "# Calculations\n", + "Ax = 0\n", + "Ay = (-Fc)*L1*math.cos(math.radians(D))/(L1+L2)\n", + "Dy = -Fc-Ay\n", + "Fab = Ay/math.sin(math.radians(A))\n", + "Fae = Fab*math.cos(math.radians(A))\n", + "Fbe = Fab*math.cos(90-E1)/math.cos(90-A)\n", + "Fbc = Fab*math.sin(math.radians(90-A))+Fbe*math.sin(math.radians(90-E1))\n", + "Fce = Fbc*math.cos(90-(180-E2-D))/math.cos(90-E2)\n", + "Fde = Fae+Fbe*math.cos(math.radians(E1))+Fce*math.cos(math.radians(E2))\n", + "Fcd = (-Fc-Fbc*math.cos(math.radians(90-E2-D)))/math.cos(math.radians(90-E1))\n", + "\t\t\t\n", + "# Results\n", + "print 'Ax = %.3f lb'%(Ax)\n", + "print 'Ay = %.2f lb'%(Ay)\n", + "print 'Dy = %.3f lb'%(Dy) \n", + "print 'Fab = %.1f lbcompression'%(Fab) \n", + "print 'Fae = %.1f lbtension'%(Fae)\n", + "print 'Fbe = %.1f lbtension'%(Fbe) \n", + "print 'Fbc = %.1f lbcompression'%(Fbc)\n", + "print 'Fce = %.1f lbcompression'%(Fce)\n", + "print 'Fde = %.1f lbtension'%(Fde)\n", + "print 'Fcd = %.1f lbcompression'%(Fcd)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ax = 0.000 lb\n", + "Ay = 250.00 lb\n", + "Dy = 750.000 lb\n", + "Fab = 288.7 lbcompression\n", + "Fae = 144.3 lbtension\n", + "Fbe = 288.7 lbtension\n", + "Fbc = 288.7 lbcompression\n", + "Fce = 288.7 lbcompression\n", + "Fde = 433.0 lbtension\n", + "Fcd = 866.0 lbcompression\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.2 Page No : 103" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "W = -100. \t\t\t#lb\n", + "angle = 45. \t\t\t#degrees\n", + "x1 = 2. \t\t\t#ft\n", + "x2 = 2. \t\t\t#ft\n", + "y1 = 2. \t\t\t#ft\n", + "y2 = 4. \t\t\t#ft\n", + "Fx = 200. \t\t\t#lb\n", + "\t\t\t\n", + "# Calculations\n", + "Cx = Fx*y1/y2\n", + "Bx = Fx+Cx\n", + "By = (y2*Bx+x1*(-W))/(x1+x2)\n", + "Cy = By\n", + "Ax = Bx\n", + "Ay = W+By\n", + "\t\t\t\n", + "# Results\n", + "print 'Ax = %.3f lb'%(Ax)\n", + "print 'Ay = %.2f lb'%(Ay)\n", + "print 'Bx = %.3f lb'%(Bx) \n", + "print 'By = %.2f lb'%(By) \n", + "print 'Cx = %.2f lb'%(Cx)\n", + "print 'Cy = %.2f lb'%(Cy) \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ax = 300.000 lb\n", + "Ay = 250.00 lb\n", + "Bx = 300.000 lb\n", + "By = 350.00 lb\n", + "Cx = 100.00 lb\n", + "Cy = 350.00 lb\n" + ] + } + ], + "prompt_number": 2 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch6.ipynb b/Statics_And_Strength_Of_Materials/ch6.ipynb new file mode 100755 index 00000000..c3a0cd02 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch6.ipynb @@ -0,0 +1,354 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:5eef56b3be475b00ee38a435583334aab5002f35094b33b8b203b9f08b619d7c" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 6 : Friction" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.1 Page No : 115" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "W = 100. \t\t\t#lb\n", + "Frictioncoefficient = 0.65\n", + "\t\t\t\n", + "# Calculations\n", + "A1 = math.degrees(math.atan(Frictioncoefficient))\n", + "\t\t\t\n", + "# Results\n", + "print 'Maximum Incliantion = %.f degrees'%(A1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum Incliantion = 33 degrees\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.2 Page No : 115" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "W = 100. \t\t\t#lb\n", + "Frictioncoefficient = 0.40\n", + "x = 3.\n", + "y = 4.\n", + "\t\t\t\n", + "# Calculations\n", + "Fmax = (W*y/(math.sqrt(x**2+y**2)))+Frictioncoefficient*W*x/(math.sqrt(x**2+y**2))\n", + "Fmin =(W*y/(math.sqrt(x**2+y**2)))-Frictioncoefficient*W*x/(math.sqrt(x**2+y**2))\n", + " \t\t\t\n", + "# Results\n", + "print 'Fmin = %.f lb'%(Fmin)\n", + "print 'Fmax =%.f lb'%(Fmax)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Fmin = 56 lb\n", + "Fmax =104 lb\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.3 Page No : 116" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "from numpy import linalg\n", + "\t\t\t\n", + "# Variables\n", + "mus = 0.25\n", + "d = 0.5 \t\t\t#in\n", + "h = 3. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "A = [[1, -1],[mus, mus]]\n", + "b = [0,1]\n", + "c = linalg.solve(A,b)\n", + "Na = c[0]\n", + "Nb = c[1]\n", + "d = -d*mus*Na+h*Na\n", + "\t\t\t\n", + "# Results\n", + "print 'minimu distance = %.2f in'%(d)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "minimu distance = 5.75 in\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.4 Page No : 118" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Ft = 1000. \t\t\t#lb\n", + "a1 = 5. \t\t\t#degrees\n", + "mu = 0.30\n", + "\t\t\t\n", + "# Calculations\n", + "R1 = Ft/math.cos(a1+math.tan(mu))\n", + "F = R1*math.sin(a1)+math.tan(mu)+math.tan(mu)/math.sin(90-math.tan(mu))\n", + "\t\t\t\n", + "# Results\n", + "print 'Forec required to start wedge = %.f lb'%(F)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Forec required to start wedge = -1705 lb\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.5 Page No : 120" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "import math \n", + "\n", + "# Variables\n", + "W = 100. \t\t\t#lb\n", + "n1 = 1/2.\n", + "n2 = 3/2.\n", + "mus = 0.40\n", + "\t\t\t\n", + "# Calculations\n", + "Ts1 = W/(math.exp(mus*n1*2*math.pi))\n", + "Ts2 = W/(math.exp(mus*n2*2*math.pi)) \n", + "\t\t\t\n", + "# Results\n", + "print 'Ts1 = %.2f lb'%(Ts1)\n", + "print 'Ts2 =%.2f lb'%(Ts2)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ts1 = 28.46 lb\n", + "Ts2 =2.31 lb\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.6 Page No : 121" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "from numpy import linalg\n", + "\t\t\t\n", + "# Variables\n", + "F = 20.\t\t\t#lb\n", + "L1 = 6. \t\t\t#in\n", + "L2 = 12. \t\t\t#in\n", + "L3 = 24. \t\t\t#in\n", + "mus = 0.60\n", + "\t\t\t\n", + "# Calculations\n", + "A =[[1,-math.exp(mus*math.pi)],[(L1+L2),(L1)]]\n", + "b =[0,F*(L1+L2+L3)]\n", + "c = linalg.solve(A,b)\n", + "TL = c[0]\n", + "Ts = c[1]\n", + "\t\t\t\n", + "# Results\n", + "print 'TL = %.2f lb'%(TL)\n", + "print 'Ts = %.2f lb'%(Ts)\n", + "\n", + "# note : answers are slightly different because of rounding off errors.\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "TL = 44.42 lb\n", + "Ts = 6.74 lb\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.7 Page No : 123" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "d = 24. \t\t\t#in\n", + "mu = 0.05\n", + "W = 2000. \t\t\t#lb\n", + "\t\t\t\n", + "# Calculations\n", + "F = W*mu*2/d\n", + "\t\t\t\n", + "# Results\n", + "print 'F = %.2f lb'%(F)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "F = 8.33 lb\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 6.8 Page No : 124" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "F = 800. \t\t\t#lb\n", + "muk = 0.10\n", + "Do = 5. \t\t\t#in\n", + "Di = 3. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "M = 2*muk*F*((Do/2)**3-(Di/2)**3)/(3*((Do/2)**2-(Di/2)**2))\n", + "\t\t\t\n", + "# Results\n", + "print 'M = %.f lb in'%(M)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "M = 163 lb in\n" + ] + } + ], + "prompt_number": 13 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch7.ipynb b/Statics_And_Strength_Of_Materials/ch7.ipynb new file mode 100755 index 00000000..5f419f92 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch7.ipynb @@ -0,0 +1,335 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:6d5c5a5b4a7fbf142c8f73269283c1a7c9727f1534378705a70188b19e0c8c3e" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 7 : Moment of Inertia" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.1 Page No : 134" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "x1 = 3. \t\t\t#in\n", + "x2 = 3. \t\t\t#in\n", + "x3 = 3. \t\t\t#in\n", + "x4 = 3. \t\t\t#in\n", + "x5 = 5. \t\t\t#in\n", + "x6 = 5. \t\t\t#in\n", + "x7 = 5. \t\t\t#in\n", + "x8 = 5. \t\t\t#in\n", + "L1 = 1. \t\t\t#in\n", + "L2 = 1. \t\t\t#in\n", + "L3 = 1. \t\t\t#in\n", + "L4 = 1. \t\t\t#in\n", + "L5 = 1. \t\t\t#in\n", + "L6 = 1. \t\t\t#in\n", + "L7 = 1. \t\t\t#in\n", + "L8 = 1. \t\t\t#in\n", + "y = 7.5 \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Ix1 = x1*L1*(y)**2\n", + "Ix2 = x2*L2*(y-L2)**2\n", + "Ix3 = x3*L3*(y-L3-L2)**2\n", + "Ix4 = x4*L4*(y-L4-L3-L2)**2\n", + "Ix5 = x5*L5*(y-L5-L4-L3-L2)**2\n", + "Ix6 = x6*L6*(y-L6-L5-L4-L3-L2)**2\n", + "Ix7 = x7*L7*(y-L7-L6-L5-L4-L3-L2)**2 \n", + "Ix8 = x8*L8*(y-L8-L7-L6-L5-L4-L3-L2)**2\n", + "Ix = Ix1+Ix2+Ix3+Ix4+Ix5+Ix6+Ix7+Ix8\n", + "\t\t\t\n", + "# Results\n", + "print 'Ix = %.f in**4'%(Ix)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ix = 552 in**4\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.2 Page No : 136" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "Iy = 60 \t\t\t#in**4\n", + "A = 25 \t\t\t #sq in\n", + "x = 10 \t\t\t #in\n", + "\t\t\t\n", + "# Calculations\n", + "Ia = Iy+ A*x**2\n", + "\t\t\t\n", + "# Results\n", + "print 'I = %.f in**4'%(Ia)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "I = 2560 in**4\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.3 Page No : 137" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "L = 5 \t\t\t#in\n", + "B = 12 \t \t\t#in\n", + "Ix = 227 \t\t\t#in**4\n", + "Iy = 10 \t\t\t#in**4\n", + "A = 10.2 \t\t\t#sq in\n", + "\t\t\t\n", + "# Calculations\n", + "Kx = math.sqrt(Ix/A)\n", + "Ky = math.sqrt(Iy/A)\n", + "\t\t\t\n", + "# Results\n", + "print 'Radius of gyration wrt x = %.2f in'%(Kx)\n", + "print 'Radius of gyration wrt y =%.2f in'%(Ky)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Radius of gyration wrt x = 4.72 in\n", + "Radius of gyration wrt y =0.99 in\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.4 Page No : 138" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "L1 = 8. \t\t\t#in\n", + "B1 = 1. \t\t\t#in\n", + "L2 = 1. \t\t\t#in\n", + "B2 = 6. \t\t\t#in\n", + "L3 = 8. \t\t\t#in\n", + "B3 = 1. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Iy = (B1*L1**3/12)+(B2*L2**3/12)+(B3*L3**3/12)\n", + "Ix = (L1*B1**3/12)+L1*B1*((B2/2)+(B1/2))**2+(L2*B2**3/12)+(L3*B3**3/12)+L3*B3*((B2/2)+(B3/2))**2\n", + "\t\t\t\n", + "# Results\n", + "print 'Ix = %.2f in**4'%(Ix)\n", + "print 'Iy =%.2f in**4'%(Iy)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ix = 215.33 in**4\n", + "Iy =85.83 in**4\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.5 Page No : 140" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "H = 8. \t\t\t#in\n", + "b = 6. \t\t\t#in\n", + "d = 4. \t\t\t#in\n", + "H1 = 5. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Ia1 = ((b*H**3)/12)+b*H*d**2\n", + "Ia2 = math.pi*((d/2)**2)+math.pi*((d/2)**2)*(H1)**2\n", + "I = Ia1-Ia2\n", + "\t\t\t\n", + "# Results\n", + "print 'I = %.2f in**4'%(I)\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "I = 697.27 in**4\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.6 Page No : 141" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "W = 64.4 \t\t\t#lb\n", + "I = 10 \t\t\t#slugft**2\n", + "g = 32.2 \t\t\t#ft/sec**2\n", + "\t\t\t\n", + "# Calculations\n", + "m = W/g\n", + "k = math.sqrt(I/m)\n", + "\t\t\t\n", + "# Results\n", + "print 'k = %.2f ft'%(k)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "k = 2.24 ft\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7.8 Page No : 142" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "d1 = 18. \t\t\t#in\n", + "d2 = 10. \t\t\t#in\n", + "d3 = 4. \t\t\t#in\n", + "Wpercuin = 0.31 \t\t\t#lb\n", + "L1 = 4. \t\t\t#in\n", + "L2 = 8. \t\t\t#in\n", + "g = 32.2 \t\t\t#ft/sec**2\n", + "\t\t\t\n", + "# Calculations\n", + "m1 = math.pi*(d1/2)**2*L1*Wpercuin/g\n", + "I1 = m1*(d1/24)**2/2\n", + "m2 = math.pi*(d2/2)**2*L2*Wpercuin/g\n", + "I2 = m2*(d2/24)**2/2\n", + "m3 = math.pi*(d3/2**2)*(L1+L2)*Wpercuin/g\n", + "I3 = m3*(d3/24)**2/2\n", + "I = I1+I2-I3\n", + "\t\t\t\n", + "# Results\n", + "print 'I = %.2f slug ft**2'%(I)\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "I = 3.28 slug ft**2\n" + ] + } + ], + "prompt_number": 7 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch8.ipynb b/Statics_And_Strength_Of_Materials/ch8.ipynb new file mode 100755 index 00000000..04957e57 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch8.ipynb @@ -0,0 +1,401 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:a09e16c002c633d5d63dac7b4f943159d57261807392bec341e8052ebf397939" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 8 : Concept of Stress" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.1 Page No : 158" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "F = -100 \t\t\t#lb\n", + "x1 = 3 \t\t\t#in\n", + "y = 6 \t\t\t#in\n", + "x2 = 24 \t\t\t#in\n", + "x3 = 12 \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Force = -F\n", + "Moment = -F*(x2+x1)\n", + "Torque = -F*y\n", + "\t\t\t\n", + "# Results\n", + "print 'Force = %.2f lb'%(Force)\n", + "print 'Moment =%.2f lb.in'%(Moment)\n", + "print 'Torque =%.2f lb.in'%(Torque)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Force = 100.00 lb\n", + "Moment =2700.00 lb.in\n", + "Torque =600.00 lb.in\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.2 Page No : 159" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "F = -5000. \t\t\t#lb\n", + "D = 250. \t\t\t#lb/ft\n", + "y1 = 4. \t\t\t#in\n", + "y2 = 2. \t\t\t#in\n", + "y3 = 5. \t\t\t#in\n", + "y4 = 3. \t\t\t#in\n", + "x = 3. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Ax = -D*y1\n", + "Ay = -F\n", + "M = (D*y1*(y2+y3+y1/2.))-F*x\n", + "\t\t\t\n", + "# Results\n", + "print 'Ax = %.2f lb'%(Ax)\n", + "print 'Ay =%.2f lb'%(Ay)\n", + "print 'M =%.2f lb.in'%(M)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ax = -1000.00 lb\n", + "Ay =5000.00 lb\n", + "M =24000.00 lb.in\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.3 Page No : 160" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "P = 5 \t\t\t#kips\n", + "angle = 30 \t\t\t#degrees\n", + "\t\t\t\n", + "# Calculations\n", + "Fn = P*math.sin(math.radians(angle))\n", + "Ft = P*math.cos(math.radians(angle))\n", + "\t\t\t\n", + "# Results\n", + "print 'Fn = %.2f lb'%(Fn)\n", + "print 'Ft =%.2f lb'%(Ft)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Fn = 2.50 lb\n", + "Ft =4.33 lb\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.4 Page No : 162" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "p = 5. \t\t\t#tons\n", + "dt = 0.75 \t\t\t#in\n", + "db = 0.5 \t\t\t#in\n", + "b = 0.5 \t\t\t#in\n", + "h = 2. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Sc = p*2000/((math.pi/4)*(dt**2))\n", + "Sr = p*2000/(b*h)\n", + "Sb = p*2000/(2*(math.pi/4)*db**2)\n", + "\t\t\t\n", + "# Results\n", + "print 'Stress in circular scetion = %.2f psi tension'%(round(Sc,-2))\n", + "print 'Stress in rectangular section = %.2f psi tension'%(Sr)\n", + "print 'Stress in bolt = %.2f psi tension'%(round(Sb,-2))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Stress in circular scetion = 22600.00 psi tension\n", + "Stress in rectangular section = 10000.00 psi tension\n", + "Stress in bolt = 25500.00 psi tension\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.5 Page No : 168" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "w = 8 \t\t\t#in\n", + "wperft = 35. \t\t\t#lb/ft\n", + "A = 10.3 \t\t\t#sq in\n", + "F1 = 3. \t\t\t#tons\n", + "F2 = 3. \t\t\t#tons\n", + "F3 = -8. \t\t\t#tons\n", + "F4 = -8. \t\t\t#tons\n", + "F5 = -5. \t\t\t#tons\n", + "F6 = -5. \t\t\t#tons\n", + "Pl = 12. \t\t\t#in\n", + "Pb = 12. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "Sa = (F1+F2)*2000/A\n", + "Sb = -(F3+F4+F1+F2)*2000/A\n", + "Sc = -(F3+F4+F1+F2+F5+F6)*2000/A\n", + "Sp = -(F3+F4+F1+F2+F5+F6)*2000/(Pl*Pb)\n", + "\t\t\t\n", + "# Results\n", + "print 'Stress in a = %.2f psi tension'%(round(Sa,-1))\n", + "print 'Stress in b = %.2f psi tension'%(round(Sb,-1))\n", + "print 'Stress in c = %.2f psi tension'%(round(Sc,-1))\n", + "print 'Stress in plate = %.f psi tension'%(Sp)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Stress in a = 1170.00 psi tension\n", + "Stress in b = 1940.00 psi tension\n", + "Stress in c = 3880.00 psi tension\n", + "Stress in plate = 278 psi tension\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.6 Page No : 166" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\t\t\t\n", + "# Variables\n", + "Ns = 8000. \t\t\t#psi\n", + "Ss = 4000. \t\t\t#psi\n", + "Ws = 25000. \t\t\t#psi\n", + "angle = 30. \t\t\t#degrees\n", + "L = 4. \t\t\t#in\n", + "b = 1. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "P = Ns*L*b/((math.cos(math.radians(2*angle)))**2)\n", + "P1 = 2*Ss*L*b/(math.sin(math.radians(2*angle)))\n", + "Pts = Ws*L*b\n", + "e = P1/Pts\n", + "\t\t\t\n", + "# Results\n", + "if (P<P1):\n", + " print 'P = %.2f lb'%(round(P,-3))\n", + "else:\n", + " print 'P1 = %.2f lb'%(round(P1,-3))\n", + "\n", + "print 'efficiency of the joint = %.2f '%(e)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "P1 = 37000.00 lb\n", + "efficiency of the joint = 0.37 \n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.7 Page No : 167" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "T = 15000 \t\t\t#psi\n", + "h1 = 3 \t\t\t#in\n", + "h2 = 2.5 \t\t\t#in\n", + "t = 0.25 \t\t\t#in\n", + "r = 5/16. \t\t\t#in\n", + "d = 1. \t\t\t#in\n", + "\t\t\t\n", + "# Calculations\n", + "P1 = T*(h1-d)*t/2.18\n", + "P2 = T*h2*t/1.7\n", + "if (P1<P2):\n", + " print 'Safe axial load = %.f lb'%(P1)\n", + "else: \n", + " print 'Safe axial load = %.f lb'%(P2)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Safe axial load = 3440 lb\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 8.8 Page No : 174" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\t\t\t\n", + "# Variables\n", + "d = 16. \t\t\t#ft\n", + "h = 24. \t \t\t#ft\n", + "P = 160. \t\t\t#lb/cu ft\n", + "hs1 = 8. \t\t\t#ft\n", + "hs2 = 8. \t\t\t#ft\n", + "hs3 = 8. \t\t\t#ft\n", + "Tsmax = 5000. \t\t\t#psi\n", + "\t\t\t\n", + "# Calculations\n", + "SW = round(P/1728,4)\n", + "P8 = round(SW*hs1*12,2)\n", + "P16 = round(SW*(hs1+hs2)*12,1)\n", + "P24 = round(SW*(hs1+hs2+hs3)*12,1)\n", + "t8 = (P8*d*12)/(2*Tsmax)\n", + "t16 = P16*d*12/(2*Tsmax)\n", + "t24 = P24*d*12/(2*Tsmax)\n", + "\t\t\t\n", + " \n", + "# Results\n", + "print 't8 = %.2f in'%(t8)\n", + "print 't16 = %.2f in'%(t16)\n", + "print 't24 = %.2f in'%(t24)\n", + "\n", + "# note : book answers are wrong. please check. " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "t8 = 0.17 in\n", + "t16 = 0.34 in\n", + "t24 = 0.51 in\n" + ] + } + ], + "prompt_number": 19 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/ch9.ipynb b/Statics_And_Strength_Of_Materials/ch9.ipynb new file mode 100755 index 00000000..ee3b4983 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/ch9.ipynb @@ -0,0 +1,509 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:9125011b4afb20e4bcb070d525e10ee7f4e22999f77038ff5e818f0e6a9e5e7f" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 9 : Concept of Strain" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.1 Page No : 193" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "length =10 \t\t\t#ft\n", + "delta =0.024 \t\t\t#in\n", + "\n", + "# Calculations\n", + "epsilon =delta/(length*12)\n", + "\n", + "# Results\n", + "print \"Axial strain =%.4f in/in\"%(epsilon)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Axial strain =0.0002 in/in\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.2 Page No : 194" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "# Variables\n", + "drop = 5. \t\t\t#in\n", + "width = 8. \t\t\t#ft\n", + "\n", + "# Calculations\n", + "deltaMB =math.sqrt((width*12/2)**2 +drop**2) - (width*12/2)\n", + "epsilon =deltaMB/(width*12/2)\n", + "\n", + "# Results\n", + "print \"Strain in the wire = %.5f in/in\"%(epsilon)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Strain in the wire = 0.00541 in/in\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.3 Page no : 198" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "# variables and Calculations\n", + "\n", + "E = 30000./0.001\n", + "deltaP1 = 66000. # psi\n", + "deltault = 116000 # psi\n", + "deltarup = 103000 # psi\n", + "\n", + "Pf = round(deltarup*math.pi/4*0.505**2,-2)\n", + "deltarup_ = round(Pf/(math.pi*0.425**2/4),-3)\n", + "percent_elongation = (2.375 - 2)/2*100\n", + "percent_reduction = ((math.pi*0.505**2/4) - (math.pi*0.425**2/4))/(math.pi*0.505**2/4) * 100\n", + "\n", + "# Results\n", + "print \"E = %.1e psi\"%E\n", + "print \"The load of failure Pf = %d lb\"%Pf\n", + "print \"Rupture strength : %.d psi\"%deltarup_\n", + "print \"Percent elongation = %.1f %%\"%percent_elongation\n", + "print \"Percent reduction in area = %.1f %%\"%percent_reduction\n", + "\n", + "\n", + "# note : last answer is wrong in book. please check." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "E = 3.0e+07 psi\n", + "The load of failure Pf = 20600 lb\n", + "Rupture strength : 145000 psi\n", + "Percent elongation = 18.8 %\n", + "Percent reduction in area = 29.2 %\n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.4 Page No : 203" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "# Variables\n", + "length =15. \t\t\t#in\n", + "tension =5000. \t\t\t#lb\n", + "UltStress =20000. \t\t\t#psi\n", + "delta =0.005 \t\t\t#in\n", + "\n", + "# Calculations\n", + "E =30*10**6 \t\t\t#psi\n", + "A1 =tension/UltStress\n", + "A2 =tension*length/(delta*E)\n", + "if A1 >= A2:\n", + " A =A1\n", + "else:\n", + " A =A2\n", + "Dia =math.sqrt(4*A/math.pi)\n", + "\n", + "# Results\n", + "print \"diameter required = %.3f in\"%( Dia)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "diameter required = 0.798 in\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.5 Page No : 204" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "# Variables\n", + "L1 =5.\n", + "L2 =10.\n", + "T1 =2.5\n", + "T2 =5.\n", + "T3 =5.\n", + "T4 =5.\n", + "T5 =2.5\n", + "E =30.*10**6 \t\t\t#psi\n", + "outDia =2. \t\t\t#in\n", + "\n", + "# Calculations # Results\n", + "inDia =1./8 \t\t\t#in\n", + "RE =(T1+T2+T3+T4+T5)/2 \t\t\t#kips\n", + "RA =RE\n", + "GH =(RA*L2-T2*L1-T1*L2)/4\n", + "print \"Stress in GH =%.1f kips\"%(GH)\n", + "A =math.pi*(outDia**2-(outDia-2*inDia)**2)/4\n", + "delta =GH*10**3 *(L1*12)/(E*A)\n", + "print \" Deformation =%.3f in\"%(delta)\n", + "sigma =GH*10**3 /A\n", + "print \" Stress =%d psi\"%(round(sigma,-3))\n", + "SF =65000/sigma\n", + "print \" Factor of safety =%.3f \"%(SF)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Stress in GH =12.5 kips\n", + " Deformation =0.034 in\n", + " Stress =17000 psi\n", + " Factor of safety =3.829 \n" + ] + } + ], + "prompt_number": 19 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.6 Page No : 205" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "from numpy import linalg\n", + "\t\t\t\n", + "# Variables\n", + "Es = 30.*10**6 \t\t\t#psi\n", + "As = 1. \t\t\t#in**2\n", + "Ea = 10.*10**6 \t\t\t#psi\n", + "Aa = 2. \t\t\t#in**2\n", + "Ls = 10. \t\t\t#ft\n", + "La = 5. \t\t\t#ft\n", + "\t\t\t\n", + "# Calculations\n", + "A =[[(Ls/(Es*As)) ,(-La/(Ea*Aa))],[1 ,1]]\n", + "b = [0,1]\n", + "c = linalg.solve(A,b)\n", + "Fa = c[0]\n", + "Fb = c[1]\n", + "d = Fb*Ls\n", + "\t\t\t\n", + "# Results\n", + "print 'distance = %.2f ft'%(d)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "distance = 5.71 ft\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.7 Page No : 206" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "# Variables\n", + "P =40000. \t\t\t#lb\n", + "L =15. \t\t\t#in\n", + "delta =0.0032 \t\t\t#in\n", + "dia =4. \t\t\t#in\n", + "axial =0.0032 \t\t\t#in\n", + "lateral =0.00022 \t\t\t#in\n", + "\n", + "# Calculations # Results\n", + "E =P*L/(delta*math.pi*(dia/2)**2)\n", + "print \"Modulus of elasticity =%.2f psi\"%(E)\n", + "Mu =lateral*L/(axial*dia)\n", + "print \"Poisson ratio = %.2f\"%(Mu)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Modulus of elasticity =14920775.91 psi\n", + "Poisson ratio = 0.26\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.8 Page No : 207" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "# Variables\n", + "alpha =11.2*10**(-6) \t\t\t#in/in/F\n", + "E =15*10**6 \t\t\t#psi\n", + "L =60. \t\t\t#in\n", + "deltaT1 =0.01 \t\t\t#in\n", + "T2 =50 \t\t\t#F\n", + "\n", + "# Calculations # Results\n", + "deltaT =deltaT1/(alpha*L)\n", + "print \"The temperature increase necessary to cause free end to touch B =%.1f F\"%(deltaT)\n", + "sigma =(alpha*L*T2-deltaT1)*E/L\n", + "print \"Stress in the rod =%d psi\"%(sigma+1)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The temperature increase necessary to cause free end to touch B =14.9 F\n", + "Stress in the rod =5900 psi\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "\n", + "Example 9.9 Page No : 208" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "weight =25000. \t\t\t#Kg\n", + "A =2 \t\t\t #sq.in\n", + "alphaS =6.5*10**(-6) \t\t\t#in/in/F\n", + "alphaB =11.2*10**(-6) \t\t\t#in/in/F\n", + "Es =30*10**6 \t\t\t#psi\n", + "Eb =15*10**6 \t\t\t#psi\n", + "\n", + "# Calculations # Results\n", + "deltaT =weight/(Es*A*(alphaB-alphaS))\n", + "print \"Net temperature drop =%.1f F\"%(deltaT)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Net temperature drop =88.7 F\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.10 Page No : 209" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math \n", + "\n", + "# Variables\n", + "S =5. \t\t\t#in\n", + "Al =6. \t\t\t#in\n", + "alphaS =6.5*10**(-6) \t\t\t#in/in/F\n", + "alphaAl =13.1*10**(-6) \t\t\t#in/in/F\n", + "Es =30.*10**6 \t\t\t#psi\n", + "EAl =10.*10**6 \t\t\t#psi\n", + "As =1. \t\t\t#in**2\n", + "AAl =2. \t\t\t#in**2\n", + "T =50. \t\t\t#F\n", + "dia =1 \t\t\t#in\n", + "\n", + "# Calculations # Results\n", + "P =(alphaS*S*12*T + alphaAl*Al*12*T)/(S*12/(Es*As) + Al*12/(EAl*AAl))\n", + "print \"Shearing force = %d lb\"%(round(P,-1))\n", + "T =P/(math.pi*(dia/2.)**2)\n", + "print \" The shear stress in the pin =%d psi\"%(round(T,-2))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Shearing force = 11900 lb\n", + " The shear stress in the pin =15200 psi\n" + ] + } + ], + "prompt_number": 21 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9.11 Page No : 211" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "edge = 2. \t\t\t#in\n", + "height =3. \t\t\t#in\n", + "F = 20000. \t\t\t#lb\n", + "deltaS = 0.00234 \t\t\t#in\n", + "deltaA = 0.00088 \t\t\t#in\n", + "\n", + "# Calculations # Results\n", + "E = F*height/(deltaA*edge*edge)\n", + "print \"Modulus of elasticity = %.1e psi\"%(E)\n", + "G =F*height/(deltaS*edge*edge)\n", + "print \" Modulus of Rigidity = %.1e psi\"%(G)\n", + "Mu =E/(2*G) -1\n", + "print \" Poisson ratio = %.1f \"%(Mu)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Modulus of elasticity = 1.7e+07 psi\n", + " Modulus of Rigidity = 6.4e+06 psi\n", + " Poisson ratio = 0.3 \n" + ] + } + ], + "prompt_number": 14 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Statics_And_Strength_Of_Materials/screenshots/ch4.png b/Statics_And_Strength_Of_Materials/screenshots/ch4.png Binary files differnew file mode 100755 index 00000000..fd954ebb --- /dev/null +++ b/Statics_And_Strength_Of_Materials/screenshots/ch4.png diff --git a/Statics_And_Strength_Of_Materials/screenshots/ch5.png b/Statics_And_Strength_Of_Materials/screenshots/ch5.png Binary files differnew file mode 100755 index 00000000..03b6529a --- /dev/null +++ b/Statics_And_Strength_Of_Materials/screenshots/ch5.png diff --git a/Statics_And_Strength_Of_Materials/screenshots/ch6.png b/Statics_And_Strength_Of_Materials/screenshots/ch6.png Binary files differnew file mode 100755 index 00000000..2b31e329 --- /dev/null +++ b/Statics_And_Strength_Of_Materials/screenshots/ch6.png |
