From d36fc3b8f88cc3108ffff6151e376b619b9abb01 Mon Sep 17 00:00:00 2001 From: kinitrupti Date: Fri, 12 May 2017 18:40:35 +0530 Subject: Revised list of TBCs --- .../chapter10_12.ipynb | 378 --------------------- 1 file changed, 378 deletions(-) delete mode 100755 Engineering_Mechanics_by_Tayal_A.K./chapter10_12.ipynb (limited to 'Engineering_Mechanics_by_Tayal_A.K./chapter10_12.ipynb') diff --git a/Engineering_Mechanics_by_Tayal_A.K./chapter10_12.ipynb b/Engineering_Mechanics_by_Tayal_A.K./chapter10_12.ipynb deleted file mode 100755 index 3563e562..00000000 --- a/Engineering_Mechanics_by_Tayal_A.K./chapter10_12.ipynb +++ /dev/null @@ -1,378 +0,0 @@ -{ - "metadata": { - "name": "chapter10.ipynb" - }, - "nbformat": 3, - "nbformat_minor": 0, - "worksheets": [ - { - "cells": [ - { - "cell_type": "heading", - "level": 1, - "metadata": {}, - "source": [ - "Chapter 10: Uniform Flexible Suspension Cables" - ] - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 10.10-1,Page No:238" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "import math\n", - "import numpy as np\n", - "\n", - "# Initilization of variables\n", - "\n", - "W1=400 # N # vertical load at pt C\n", - "W2=600 # N # vertical load at pt D\n", - "W3=400 # N # vertical load at pt E\n", - "l=2 # m # l= Lac=Lcd=Lde=Leb\n", - "h=2.25 # m # distance of the cable from top\n", - "L=2 # m # dist of A from top\n", - "\n", - "# Calculations\n", - "\n", - "# Solving eqn's 1&2 using MATRIX for Xb & Yb\n", - "A=np.array([[-L ,4*l],[-h, 2*l]])\n", - "B=np.array([((W1*l)+(W2*2*l)+(W1*3*l)),(W1*l)])\n", - "C=np.linalg.solve(A,B)\n", - "\n", - "# Now consider the F.B.D of BE, Take moment at E\n", - "y_e=(C[1]*l)/C[0] # m / here y_e is the distance between E and the top\n", - "theta_1=arctan(y_e/l)*(180/pi) # degree # where theta_1 is the angle between BE and the horizontal\n", - "T_BE=C[0]/cos(theta_1*(pi/180)) # N (T_BE=T_max)\n", - "\n", - "# Now consider the F.B.D of portion BEDC\n", - "# Take moment at C\n", - "y_c=((C[1]*6)-(W3*4)-(W2*2))/(C[0]) # m\n", - "theta_4=arctan(((y_c)-(l))/(l))*(180/pi) # degree\n", - "T_CA=C[0]/cos(theta_4*(pi/180)) # N # Tension in CA\n", - "\n", - "# Results\n", - "\n", - "print\"(i) The horizontal reaction at B (Xb) is \",round(C[0]),\"N\"\n", - "print\"(i) The vertical reaction at B (Yb) is \",round(C[1]),\"N\"\n", - "print\"(ii) The sag at point E (y_e) is \",round(y_e,3),\"m\"\n", - "print\"(iii) The tension in portion CA (T_CA) is \",round(T_CA,1),\"N\"\n", - "print\"(iv) The max tension in the cable (T_max) is \",round(T_BE,1),\"N\" #answer varies due to decimal variance\n", - "print\"(iv) The max slope (theta_1) in the cable is \",round(theta_1,1),\"degree\"\n" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "(i) The horizontal reaction at B (Xb) is 1600.0 N\n", - "(i) The vertical reaction at B (Yb) is 1100.0 N\n", - "(ii) The sag at point E (y_e) is 1.375 m\n", - "(iii) The tension in portion CA (T_CA) is 1627.9 N\n", - "(iv) The max tension in the cable (T_max) is 1941.6 N\n", - "(iv) The max slope (theta_1) in the cable is 34.5 degree\n" - ] - } - ], - "prompt_number": 4 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 10.10-2,Page No:241" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "import math\n", - "\n", - "# Initiization of variables\n", - "\n", - "W1=100 # N # Pt load at C\n", - "W2=150 # N # Pt load at D\n", - "W3=200 # N # Pt load at E\n", - "l=1 # m # l=Lac=Lcd=Lde=Leb\n", - "h=2 # m # dist between Rb & top\n", - "Xa=200 # N\n", - "Xb=200 # N\n", - "\n", - "# Calculations\n", - "\n", - "# consider the F.B.D of entire cable\n", - "# Take moment at A\n", - "Yb=((W1*l)+(W2*2*l)+(W3*3*l)-(Xb*h))/(4*l) # N\n", - "Ya=W1+W2+W3-Yb # N # sum Fy=0\n", - "# Now consider the F.B.D of AC\n", - "\n", - "# Take moment at C,\n", - "y_c=(Ya*l)/Xa # m\n", - "theta_1=arctan(y_c/l)*(180/pi) # degree\n", - "T_AC=Xa/cos(theta_1*(pi/180)) # N # T_AC*cosd(theta_1)=horizontal component of tension in the cable\n", - "# here, T_AC=T_max\n", - "T_max=(Xa**2+Ya**2)**0.5 # N\n", - "T_AC=T_max\n", - "\n", - "# Now consider the F.B.D of portion ACD\n", - "y_d=((Ya*2*l)-(W1*l))/(Xa) # m # taking moment at D\n", - "theta_2=arctan(((y_d)-(y_c))/(l))*(180/pi) # degree\n", - "T_CD=Xa/(cos(theta_2*(pi/180))) # N \n", - "\n", - "# Results\n", - "\n", - "print\"(i) The component of support reaction at A (Ya) is \",round(Ya),\"N\"\n", - "print\"(i) The component of support reaction at B (Yb) is \",round(Yb),\"N\"\n", - "print\"(ii) The tension in portion AC (T_AC) of the cable is \",round(T_AC,1),\"N\"\n", - "print\"(ii) The tension in portion CD (T_CD) of the cable is \",round(T_CD,1),\"N\"\n", - "print\"(iii) The max tension in the cable is \",round(T_max,1),\"N\"\n" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "(i) The component of support reaction at A (Ya) is 300.0 N\n", - "(i) The component of support reaction at B (Yb) is 150.0 N\n", - "(ii) The tension in portion AC (T_AC) of the cable is 360.6 N\n", - "(ii) The tension in portion CD (T_CD) of the cable is 282.8 N\n", - "(iii) The max tension in the cable is 360.6 N\n" - ] - } - ], - "prompt_number": 14 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 10.10-3,Page No:246" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "import math\n", - "\n", - "# Initilization of variables\n", - "\n", - "w=75 # kg/m # mass per unit length of thw pipe\n", - "l=20 # m # dist between A & B\n", - "g=9.81 #m/s^2 # acc due to gravity\n", - "y=2 # m # position of C below B\n", - "\n", - "# Calculations\n", - "\n", - "# Let x_b be the distance of point C from B \n", - "# In eq'n x_b^2+32*x_b-320=0\n", - "a=1\n", - "b=32\n", - "c=-320\n", - "x_b=(-b+(b**2-(4*a*c))**0.5)/(2*a) # m # we get x_b by equating eqn's 1&2\n", - "\n", - "# Now tension T_0\n", - "T_0=((w*g*x_b**2)/(2*y))*(10**-3) #kN # from eq'n 1\n", - "\n", - "# Now the max tension occurs at point A,hence x is given as,\n", - "x=20-x_b # m\n", - "w_x=w*g*x*10**(-3) # kN \n", - "T_max=((T_0)**2+(w_x)**2)**0.5 # kN # Maximum Tension\n", - "\n", - "# Results\n", - "\n", - "print\"The lowest point C which is situated at a distance (x_b) from support B is \",round(x_b),\"m\"\n", - "print\"The maximum tension (T_max) in the cable is \",round(T_max,2),\"kN\"\n", - "print\"The minimum tension (T_0) in the cable is \",round(T_0,2),\"kN\"\n" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "The lowest point C which is situated at a distance (x_b) from support B is 8.0 m\n", - "The maximum tension (T_max) in the cable is 14.71 kN\n", - "The minimum tension (T_0) in the cable is 11.77 kN\n" - ] - } - ], - "prompt_number": 18 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 10.10-4,Page No:247" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "import math\n", - "\n", - "# Initilization of variables\n", - "\n", - "m=0.5 # kg/m # mass of the cable per unit length\n", - "g=9.81 # m/s^2\n", - "x=30#/ m # length AB\n", - "y=0.5 # m # dist between C & the horizontal\n", - "x_b=15 # m # dist of horizontal from C to B\n", - "\n", - "# Calculations\n", - "\n", - "w=m*g # N/m # weight of the cable per unit length\n", - "T_0=(w*x_b**2)/(2*y) # N # From eq'n 1\n", - "T_B=((T_0)**2+(w*x*0.5)**2)**0.5 # N # Tension in the cable at point B\n", - "W=T_B # N # As pulley is frictionless the tension in the pulley on each side is same,so W=T_B\n", - "\n", - "\n", - "# Slope of the cable at B,\n", - "theta=arccos(T_0/T_B)*(180/pi) # degree\n", - "\n", - "# Now length of the cable between C & B is,\n", - "S_cb =x_b*(1+((2/3)*(y/x_b)*0.5)) # m\n", - "\n", - "\n", - "# Now total length of the cable AB is,\n", - "S_ab=2*S_cb # m \n", - "\n", - "# Results\n", - "\n", - "print\"(i) The magnitude of load W is \",round(W),\"N\"\n", - "print\"(ii) The angle of the cable with the horizontal at B is \",round(theta,1),\"degree\"\n", - "print\"(iii) The total length of the cable AB is \",round(S_ab),\"m\"\n" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "(i) The magnitude of load W is 1106.0 N\n", - "(ii) The angle of the cable with the horizontal at B is 3.8 degree\n", - "(iii) The total length of the cable AB is 30.0 m\n" - ] - } - ], - "prompt_number": 116 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 10.10-5,Page No:249" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "import math\n", - "\n", - "# Initilization of variables\n", - "\n", - "x=30 # m # distance between two electric poles\n", - "Tmax=400 # N # Max Pull or tension\n", - "w=3 # N/m # weight per unit length of the cable\n", - "\n", - "# Calculations\n", - "\n", - "# The cable is assumed to be parabolic in shape, its eq'n is y=w*x^2/2*T_0.....(eq'n 1). Substuting the co-ordinates of point B (l/2,h), where h is the sag in the cable.This gives, T_0=(w*(l/2)^2)/(2*h)=wl^2/8*h\n", - "# Now the maximum pull or tension occurs at B,\n", - "T_B=Tmax # N \n", - "# Hence T_B=Tmax=sqrt(T_0^2+(w*l/2)^2). On simplyfyingthis eq'n we get, \n", - "h=(x**2)**0.5/((16*(((Tmax*2)**2/(w*x)**2)))+(1))**0.5 # m \n", - "\n", - "# Results \n", - "\n", - "print\"The smallest value of the sag in the cable is \",round(h,3),\"m\"\n" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "The smallest value of the sag in the cable is 0.843 m\n" - ] - } - ], - "prompt_number": 72 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 10.10-6,Page No:252" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "import math\n", - "\n", - "# Initilization of variables\n", - "\n", - "l=200 # m # length of the cable\n", - "m=1000 # kg # mass of the cable\n", - "S=50 # m # sag in the cable\n", - "s=l/2 # m\n", - "g=9.81 # m/s^2\n", - "\n", - "# Calculations\n", - "\n", - "w=(m*g)/l # N/m # mass per unit length of the cable\n", - "# Substuting the values s=l/2 & y=c+S in eq'n 1 to get the value of c,\n", - "c=7500/100 # m \n", - "Tmax=((w*c)**2+(w*s)**2)**0.5 # N # Maximum Tension\n", - "# To determine the span (2*x) let us use the eq'n of catenary, y=c*cosh(x/c), where y=c+50. On simplyfying we get y/c=cosh(x/c), here let y/c=A\n", - "y=c+50\n", - "\n", - "A=1.666 \n", - "x=c*(arccosh(A))*(pi/180)# m \n", - "L=2*x*(180/pi) # m # where L= span\n", - "\n", - "# Results\n", - "\n", - "print\"The horizontal distance between the supports and the max Tension (L) is \",round(L),\"m\"\n" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "The horizontal distance between the supports and the max Tension (L) is 165.0 m\n" - ] - } - ], - "prompt_number": 114 - } - ], - "metadata": {} - } - ] -} \ No newline at end of file -- cgit