diff options
Diffstat (limited to 'backup/Engineering_Mechanics_by_Tayal_A.K._version_backup/chapter12.ipynb')
| -rwxr-xr-x | backup/Engineering_Mechanics_by_Tayal_A.K._version_backup/chapter12.ipynb | 411 |
1 files changed, 411 insertions, 0 deletions
diff --git a/backup/Engineering_Mechanics_by_Tayal_A.K._version_backup/chapter12.ipynb b/backup/Engineering_Mechanics_by_Tayal_A.K._version_backup/chapter12.ipynb new file mode 100755 index 00000000..9f1a8645 --- /dev/null +++ b/backup/Engineering_Mechanics_by_Tayal_A.K._version_backup/chapter12.ipynb @@ -0,0 +1,411 @@ +{
+ "metadata": {
+ "name": "chapter12.ipynb"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter 12: Moment Of Inertia"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 12.12-7,Page No:285"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "# Initilization of variables\n",
+ "\n",
+ "A= 50 # cm^2 # area of the shaded portion\n",
+ "J_A=22.5*10**2 # cm^4 # polar moment of inertia of the shaded portion\n",
+ "d=6 # cm\n",
+ "\n",
+ "# Calculations\n",
+ "\n",
+ "J_c=J_A-(A*d**2) \n",
+ "\n",
+ "# substuting the value of I_x from eq'n 2 in eq'n 1 we get,\n",
+ "I_y=J_c*0.333 # cm^4 # M.O.I about Y-axis # 1/3=0.333\n",
+ "\n",
+ "# Now from eq'n 2,\n",
+ "I_x=2*I_y # cm^4 # M.O.I about X-axis\n",
+ "\n",
+ "# Results\n",
+ "\n",
+ "print\"The centroidal moment of inertia about X-axis (I_x) is \",round(I_x),\"cm^4\"\n",
+ "print\"The centroidal moment of inertia about Y-axis (I_y) is \",round(I_y),\"cm^4\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The centroidal moment of inertia about X-axis (I_x) is 300.0 cm^4\n",
+ "The centroidal moment of inertia about Y-axis (I_y) is 150.0 cm^4\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 12.12-8,Page No:288"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "# Initilization of variables\n",
+ "\n",
+ "b=20 # cm # width of the pate\n",
+ "d=30 # cm # depth of the plate\n",
+ "r=15 # cm # radius of the circular hole\n",
+ "h=20 # cm # distance between the centre of the circle & the x-axis\n",
+ "\n",
+ "# Calculations\n",
+ "\n",
+ "# (a) Location of the centroid of the composite area\n",
+ "\n",
+ "A_1=b*d # cm^2 # area of the plate\n",
+ "y_1=d/2 # cm # y-coordinate of the centroid\n",
+ "A_2=(pi*r**2)/4 # cm^2 # area of the circle removed (negative)\n",
+ "y_2=h # cm # y-coordinate of the centroid\n",
+ "y_c=((A_1*y_1)-(A_2*y_2))/(A_1-A_2) # cm # from the bottom edge\n",
+ "\n",
+ "# (b) Moment of Inertia of the composite area about the centroidal x-axis\n",
+ "\n",
+ "# Area (A_1) M.I of area A_1 about x-axis\n",
+ "I_x1=(b*(d**3))/12 # cm^4\n",
+ "\n",
+ "# M.I of the area A_1 about the centroidal x-axis of the composite area (By parallel-axis theorem)\n",
+ "OC_1=15 # cm # from the bottom edge\n",
+ "OC_2=20 # cm\n",
+ "OC=12.9 # cm # from the bottom edge\n",
+ "d_1=OC_1-OC # cm\n",
+ "d_2=OC_2-OC # cm \n",
+ "I_X1=(I_x1)+(A_1*d_1**2) # cm^4\n",
+ "\n",
+ "# Area(A_2) M.I of area A_2 about x-axis\n",
+ "I_x2=(pi*r**4)/64 # cm^2\n",
+ "\n",
+ "#M.I of the area A_2 about the centroidal x-axis of the composite area (By parallel-axis theorem)\n",
+ "I_X2=(I_x2)+(A_2*d_2**2) # cm^4\n",
+ "\n",
+ "# COMPOSITE AREA:M.O.I of the composite area about the centroidal x-axis\n",
+ "I_x=(I_X1)-(I_X2) # cm^4\n",
+ "\n",
+ "# Results\n",
+ "\n",
+ "print\"The M.O.I of the composite area about the centroidal x-axis is \",round(I_x),\"cm^4\" \n",
+ "#due to decimal variance answer varies by 1cm^4 from textbook"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The M.O.I of the composite area about the centroidal x-axis is 36253.0 cm^4\n"
+ ]
+ }
+ ],
+ "prompt_number": 41
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 12.12-9,Page No:289"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "# Initilization of variables\n",
+ "\n",
+ "b1=80 # mm # width of the flange pate\n",
+ "d1=20 # mm # depth of the flange plate\n",
+ "b2=40 # mm # width/thickness of the web\n",
+ "d2=60 # mm # depth of the web\n",
+ "\n",
+ "# Calculations\n",
+ "\n",
+ "# (a) Location of the centroid of the composite area\n",
+ "A_1=b1*d1 # mm^2 # area of the flange plate\n",
+ "y_1=d2+(d1/2) # mm # y-coordinate of the centroid\n",
+ "A_2=b2*d2 # mm^2 # area of the web\n",
+ "y_2=d2/2 # mm # y-coordinate of the centroid\n",
+ "y_c=((A_1*y_1)+(A_2*y_2))/(A_1+A_2) # mm # from the bottom edge\n",
+ "\n",
+ "# (b) Moment of Inertia of the composite area about the centroidal x-axis\n",
+ "\n",
+ "# Area (A_1) M.I of area A_1 about x-axis\n",
+ "I_x1=(b1*(d1**3))/12 # mm^4\n",
+ "# M.I of the area A_1 about the centroidal x-axis of the composite area (By parallel-axis theorem)\n",
+ "OC_1=70 # mm # from the bottom edge\n",
+ "OC_2=30 # mm # from the bottom edge\n",
+ "OC=y_c # mm # from the bottom edge\n",
+ "d_1=(d2-y_c)+(d1/2) # mm\n",
+ "d_2=y_c-OC_2 # mm \n",
+ "I_X1=(I_x1)+(A_1*d_1**2) # mm^4\n",
+ "\n",
+ "# Area(A_2) M.I of area A_2 about x-axis\n",
+ "I_x2=(b2*d2**3)/12 # mm^4\n",
+ "# M.I of the area A_2 about the centroidal x-axis of the composite area (By parallel-axis theorem)\n",
+ "I_X2=(I_x2)+(A_2*d_2**2) # mm^4\n",
+ "# COMPOSITE AREA:M.O.I of the composite area about the centroidal x-axis\n",
+ "I_x=(I_X1)+(I_X2) # mm^4\n",
+ "\n",
+ "# Results\n",
+ "\n",
+ "print\"The M.O.I of the composite area about the centroidal x-axis is \",round(I_x),\"mm^4\"\n",
+ "# NOTE: The answer given in the text book is 2.31*10^3 insted of 2.31*10^6.\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The M.O.I of the composite area about the centroidal x-axis is 2309333.0 mm^4\n"
+ ]
+ }
+ ],
+ "prompt_number": 42
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 12.12-10,Page No:291"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "import math\n",
+ "\n",
+ "# Initilization of variables\n",
+ "b1=120 # mm # width of the flange pate of L-section\n",
+ "d1=20 # mm # depth of the flange plate\n",
+ "b2=20 # mm # width/thickness of the web\n",
+ "d2=130 # mm # depth of the web\n",
+ "\n",
+ "# Calculations\n",
+ "\n",
+ "# (a) Location of the centroid of the composite area\n",
+ "A_1=b1*d1 # mm^2 # area of the flange plate\n",
+ "A_2=b2*d2 # mm^2 # area of the web\n",
+ "y_1=d2+(d1/2) # mm # y-coordinate of the centroid\n",
+ "y_2=d2/2 # mm # y-coordinate of the centroid\n",
+ "x_1=60 # mm # x-coordinate of the centroid\n",
+ "x_2=110 # mm # x-coordinate of the centroid\n",
+ "y_c=((A_1*y_1)+(A_2*y_2))/(A_1+A_2) # mm # from the bottom edge\n",
+ "x_c=((A_1*x_1)+(A_2*x_2))/(A_1+A_2) # mm # from the bottom edge\n",
+ "\n",
+ "# (b) Moment of Inertia of the composite area about the centroidal x-axis\n",
+ "\n",
+ "# Area (A_1) M.I of area A_1 about x-axis\n",
+ "I_x1=(b1*(d1**3))/12 # mm^4\n",
+ "# M.I of the area A_1 about the centroidal x-axis of the composite area (By parallel-axis theorem)\n",
+ "OC_1=d2+(d1/2) # mm # from the bottom edge\n",
+ "OC_2=d2/2 # mm # from the bottom edge\n",
+ "OC=y_c # mm # from the bottom edge\n",
+ "d_1=(d2-y_c)+(d1/2) # mm\n",
+ "d_2=y_c-OC_2 # mm \n",
+ "I_X1=(I_x1)+(A_1*d_1**2) # mm^4\n",
+ "\n",
+ "# Area(A_2) M.I of area A_2 about x-axis\n",
+ "I_x2=(b2*d2**3)/12 # mm^4\n",
+ "\n",
+ "# M.I of the area A_2 about the centroidal x-axis of the composite area (By parallel-axis theorem)\n",
+ "I_X2=(I_x2)+(A_2*d_2**2) # mm^4\n",
+ "# COMPOSITE AREA:M.O.I of the composite area about the centroidal x-axis\n",
+ "I_x=(I_X1)+(I_X2) # mm^4\n",
+ "\n",
+ "# (c) Moment of Inertia of the composite area about the centroidal y-axis\n",
+ "\n",
+ "# Area (A_1) M.I of area A_1 about y-axis\n",
+ "I_y1=(d1*(b1**3))/12 # mm^4\n",
+ "# M.I of the area A_1 about the centroidal y-axis of the composite area (By parallel-axis theorem)\n",
+ "d_3=x_c-(b1/2) # mm # distance between c &c1 along x axis\n",
+ "I_Y1=(I_y1)+(A_1*d_3**2) # mm^4\n",
+ "\n",
+ "# Area(A_2) M.I of area A_2 about y-axis\n",
+ "I_y2=(d2*b2**3)/12 # mm^4\n",
+ "\n",
+ "# M.I of the area A_2 about the centroidal y-axis of the composite area (By parallel-axis theorem)\n",
+ "d_4=b1-x_c-(b2/2) # mm # distance between c &c2 along x axis\n",
+ "I_Y2=(I_y2)+(A_2*d_4**2) # mm^4\n",
+ "# COMPOSITE AREA:M.O.I of the composite area about the centroidal y-axis\n",
+ "I_y=(I_Y1)+(I_Y2) # mm^4\n",
+ "\n",
+ "# Results\n",
+ "\n",
+ "print\"The M.O.I of the composite area about the centroidal x-axis is \",round(I_x),\"mm^4\"\n",
+ "print\"The M.O.I of the composite area about the centroidal Y-axis is \",round(I_y),\"mm^4\"\n",
+ "# NOTE: The answer for I_x given in text book is 0.76*10^6 insted of 10.76*10^6\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The M.O.I of the composite area about the centroidal x-axis is 10761666.0 mm^4\n",
+ "The M.O.I of the composite area about the centroidal Y-axis is 6086666.0 mm^4\n"
+ ]
+ }
+ ],
+ "prompt_number": 43
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 12.12-14,Page No:299"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "# Initilization of variables\n",
+ "\n",
+ "b=1 # cm # smaller side of the L-section\n",
+ "h=4 # cm # larger side of the L-section\n",
+ "\n",
+ "# Calculations\n",
+ "\n",
+ "# (A) RECTANGLE A_1: Using the paralel axis theorem\n",
+ "Ixy=0\n",
+ "I_xy1=(Ixy)+((h*b)*(b*0.5)*(h*0.5)) # cm**4\n",
+ "\n",
+ "# (B) RECTANGLE A_2: Using the paralel axis theorem\n",
+ "I_xy2=(Ixy)+((b*(h-1))*(1+(1.5))*(b*0.5)) # cm**4\n",
+ "\n",
+ "# Product of inertia of the total area\n",
+ "I_xy=I_xy1+I_xy2 # cm**4\n",
+ "\n",
+ "# Results\n",
+ "\n",
+ "print\"The Product of inertia of the L-section is \",round(I_xy,2),\"cm^4\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The Product of inertia of the L-section is 7.75 cm^4\n"
+ ]
+ }
+ ],
+ "prompt_number": 47
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 12.12-15,Page No:300"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "# Initilization of variables\n",
+ "\n",
+ "I_x=1548 # cm^4 # M.O.I of the Z-section about X-axis\n",
+ "I_y=2668 # cm^4 # M.O.I of the Z-section about Y-axis\n",
+ "b=12 # cm # width of flange of the Z-section\n",
+ "d=3 # cm # depth of flange of the Z-section\n",
+ "t=2 # cm # thickness of the web of the Z-section\n",
+ "h=6 # cm # depth of the web of the Z-section\n",
+ "\n",
+ "#Calculations\n",
+ "\n",
+ "A_1=b*d # cm^2 # area of top flange\n",
+ "x_1=-5 # cm # distance of the centroid from X-axis for top flange\n",
+ "y_1=4.5 # cm # distance of the centroid from Y-axis for top flange\n",
+ "A_2=t*h # cm^2 # area of web\n",
+ "x_2=0 # cm # distance of the centroid from X-axis for the web\n",
+ "y_2=0 # cm # distance of the centroid from Y-axis for the web\n",
+ "A_3=b*d # cm^2 # area of bottom flange\n",
+ "x_3=5 # cm # distance of the centroid from X-axis for top flange\n",
+ "y_3=-4.5 # cm # distance of the centroid from Y-axis for top flange\n",
+ "\n",
+ "# Product of Inertia of the total area is,\n",
+ "I_xy=((A_1*x_1*y_1)+(A_3*x_3*y_3)) # cm^4\n",
+ "# The direction of the principal axes is,\n",
+ "theta_m=(arctan((2*I_xy)/(I_y-I_x))*(180/pi))/2 # degree\n",
+ "# Principa M.O.I\n",
+ "I_max=((I_x+I_y)/2)+((((I_x-I_y)/2)**2+(I_xy)**2)**0.5) # cm**4\n",
+ "I_mini=((I_x+I_y)/2)-((((I_x-I_y)/2)**2+(I_xy)**2)**0.5) # cm**4\n",
+ "\n",
+ "# Results\n",
+ "\n",
+ "print\"The principal axes of the section about O is \",round(theta_m,2),\"degree\"\n",
+ "print\"The Maximum value of principal M.O.I is \",round(I_max),\"cm^4\"\n",
+ "print\"The Minimum value of principal M.O.I is \",round(I_mini),\"cm^4\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The principal axes of the section about O is -35.47 degree\n",
+ "The Maximum value of principal M.O.I is 3822.0 cm^4\n",
+ "The Minimum value of principal M.O.I is 394.0 cm^4\n"
+ ]
+ }
+ ],
+ "prompt_number": 50
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |