From 476705d693c7122d34f9b049fa79b935405c9b49 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 14 Apr 2020 10:19:27 +0530 Subject: Initial commit --- .../6-Bending.ipynb | 1099 ++++++++++++++++++++ 1 file changed, 1099 insertions(+) create mode 100644 Mechanics_of_Materials_by_R_C_Hibbeler/6-Bending.ipynb (limited to 'Mechanics_of_Materials_by_R_C_Hibbeler/6-Bending.ipynb') diff --git a/Mechanics_of_Materials_by_R_C_Hibbeler/6-Bending.ipynb b/Mechanics_of_Materials_by_R_C_Hibbeler/6-Bending.ipynb new file mode 100644 index 0000000..78011e0 --- /dev/null +++ b/Mechanics_of_Materials_by_R_C_Hibbeler/6-Bending.ipynb @@ -0,0 +1,1099 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 6: Bending" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.11: B11.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.11 : ')\n", +"\n", +"//Given:\n", +"l = 4.5; //m\n", +"R1 = 1.5; //kN\n", +"R2 = 3; //kN\n", +"uvl = 2; //kN/m\n", +"\n", +"//Shear diagram:\n", +"x = sqrt((2*R1*l)/(uvl));\n", +"M = (R1*x) - (0.5*uvl*x^3)/(3*l);\n", +"\n", +"//Display:\n", +"\n", +" \n", +" printf('\n\nV becomes zero at x = %1.1fm',x);\n", +" printf('\nThe magnitude of the maximum moment = %1.1f kNm',M);\n", +" \n", +"//-----------------------------------------------------------------END--------------------------------------------------------------------------\n", +"\n", +"\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.13: B13.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.13 : ')\n", +"\n", +"//Given:\n", +"l_ab = 4; //m\n", +"l_cd = 4; //m\n", +"l_bc = 6; //m\n", +"Rb = 8; //kN\n", +"uvl = 2; //kN/m\n", +"\n", +"//Moment diagram:\n", +"p = [-1/18 0 -3.6 17.6]\n", +"x = roots(p)\n", +"y = x(3);\n", +"\n", +"//Display:\n", +" \n", +" printf('\n\nV becomes zero at x = %1.2f m',y);\n", +"\n", +" \n", +"//-----------------------------------------------------------------END--------------------------------------------------------------------------\n", +"\n", +"\n", +"\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.14: B14.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.14 : ')\n", +"\n", +"//Given:\n", +"b = 60; //mm\n", +"h = 120; //mm\n", +"sigma_max = 20; //N/mm^2\n", +"c = b;\n", +"\n", +"//Part (a):\n", +"I = (1/12)*b*h^3;\n", +"M1 = (sigma_max*I)/(c); //sigma_max = Mc/I Flexure Formula\n", +"M1 = M1*10^-6; //in kN/m\n", +"\n", +"//Part (b):\n", +"y0=60;\n", +"y1=-60\n", +"\n", +"M2 = integrate('-(20*y^2)','y',y0,y1);\n", +"M2 = M2*10^-6;\n", +"\n", +"F = (0.5*sigma_max*b*b);\n", +"c = 2*(60 -(0.5*b)); //distance between centroids of both the volumes.\n", +"M = F*c/1000;\n", +"\n", +"//Display:\n", +" \n", +" printf('\n\nThe internal moment M calculated using : ');\n", +" printf('\na)The flexure formula = %1.2f kNm',M1);\n", +" printf('\nb)The resultant of the stress distribution using the basic principles = %1.2f kNm',M2);\n", +"\n", +" \n", +"//-----------------------------------------------------------------END--------------------------------------------------------------------------\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.15: B15.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.15 : ')\n", +"\n", +"//Given:\n", +"udl = 5; //kN/m\n", +"l1 = 3;//m\n", +"l2 = 6; //m\n", +"t = 20/1000; //mm\n", +"yb = 0.15;//m\n", +"\n", +"//Section Property:\n", +"I_bar1 = (1/12)*(0.25)*(0.02^3);\n", +"Ad2 = (0.25)*(0.02)*(yb+(t/2))^2;\n", +"I_bar2 = (1/12)*(0.02)*(0.3^3);\n", +"I = 2*(I_bar1 + Ad2) + I_bar2;\n", +"\n", +"//Bending stress:\n", +"c = 0.15 + t;\n", +"M= 22.5; //kNm\n", +"\n", +"sigma_max = (M*c)/(I*1000);\n", +"\n", +"sigma_B = (M*yb)/(I*1000);\n", +"\n", +"//Display:\n", +"\n", +" printf('\n\nThe absolute maximum bending stress is = %1.1f MPa',sigma_max);\n", +"\n", +" \n", +"//-----------------------------------------------------------------END--------------------------------------------------------------------------\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.16: B16.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.16 : ')\n", +"\n", +"//Given:\n", +"t1 = 15/1000; //m\n", +"t2 = 20/1000; //m\n", +"l = 250/1000; //m\n", +"b = 200/1000; //m\n", +"P = 2.4; //kN\n", +"l_a = 2; //m\n", +"l_b = 1; //m\n", +"\n", +"//Internal Moment:\n", +"y1 = b/2;\n", +"y2 = t2/2;\n", +"A = (2*t1*b)+(t2*l);\n", +"y_bar = ((2*y1*t1*b)+(y2*t2*l))/A;\n", +"\n", +"M = (P*l_a)+(1*y_bar);\n", +"\n", +"//Section Property:\n", +"I1 = (1/12)*(l*t2^3) + (l*t2*(y_bar - y2)^2);\n", +"I2 = (1/12)*(t1*b^3) + (t1*b*(y1 - y_bar)^2);\n", +"I =I1+ 2*I2;\n", +"\n", +"//Maximum Bending Stress:\n", +"c = b - y_bar;\n", +"sigma_max = (M*c)/(I*1000);\n", +"\n", +"//Display:\n", +" \n", +" printf('\n\nThe maximum bending stress at section a-a = %1.1f MPa',sigma_max);\n", +"\n", +" \n", +"//-----------------------------------------------------------------END--------------------------------------------------------------------------\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.17: B17.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.17 : ')\n", +"\n", +"//Given:\n", +"b = 60/1000; //m\n", +"h = 30/1000; //m\n", +"M = 40; //Nm\n", +"c1= h/2;\n", +"rib_t = 5/1000; //m\n", +"rib_w = 10/1000;//m\n", +"\n", +"//Without Ribs:\n", +"I1 = (1/12)*(b*h^3);\n", +"sigma_max1 = (M*c1)/(I1*10^6);\n", +"\n", +"//With Ribs:\n", +"y1 = c1;\n", +"y2 = h+(rib_t/2);\n", +"A1 = h*b;\n", +"A2 = rib_t*rib_w;\n", +"y_bar = ((y1*A1)+2*(y2*A2))/(A1 + 2*A2);\n", +"\n", +"c2 = h+rib_t - y_bar;\n", +"I2 = I1 + (b*h*(y_bar - y1)^2);\n", +"I3 = (1/12)*rib_w*rib_t^3 + (rib_w*rib_t*(y2 - y_bar)^2);\n", +"I = I2 + 2*I3;\n", +"\n", +"sigma_max2 = (M*c2)/(I*10^6);\n", +"\n", +"if(sigma_max2>sigma_max1)\n", +" \n", +" printf('\n\nThe maximum normal stress in the member without ribs = %1.2f MPa',sigma_max1);\n", +" printf('\nThe maximum normal stress in the member with ribs = %1.2f MPa',sigma_max2);\n", +" printf('\nThe ribs should be omitted.');\n", +" \n", +" end\n", +"\n", +" \n", +"//-----------------------------------------------------------------END--------------------------------------------------------------------------\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.18: B18.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.18 : ')\n", +"\n", +"//Given:\n", +"M = 12; //kNm\n", +"l_bc = 0.2; //m\n", +"l_be = 0.4; //m\n", +"\n", +"//Internal Moment Components:\n", +"My = (-4/5)*M;\n", +"Mz = (3/5)*M;\n", +"\n", +"Iy = (1/12)*(l_be*l_bc^3);\n", +"Iz = (1/12)*(l_bc*l_be^3); \n", +"\n", +"//Bending Stress:\n", +"sigma_B = (-Mz*1000*(l_be/2))/Iz + (My*1000*(-l_bc/2))/Iy;\n", +"sigma_B = sigma_B/10^6;\n", +"sigma_C = (-Mz*1000*(l_be/2))/Iz + (My*1000*(l_bc/2))/Iy;\n", +"sigma_C = sigma_C/10^6;\n", +"sigma_D = (-Mz*1000*(-l_be/2))/Iz + (My*1000*(l_bc/2))/Iy;\n", +"sigma_D = sigma_D/10^6;\n", +"sigma_E = (-Mz*1000*(-l_be/2))/Iz + (My*1000*(-l_bc/2))/Iy;\n", +"sigma_E = sigma_E/10^6;\n", +"\n", +"//Orientation of Nuetral Axis:\n", +"z = (0.45)/(sigma_E + sigma_B);\n", +"\n", +"//theta = -atan(4/3);\n", +"tanA = (Iz/Iy)*(-4/3);\n", +"alpha = atan(tanA);\n", +"alpha = alpha*(180/%pi);\n", +"\n", +"\n", +"//Display:\n", +"\n", +" \n", +" printf('\n\nThe normal stress at B = %1.2f MPa',sigma_B);\n", +" printf('\nThe normal stress at C = %1.2f MPa',sigma_C);\n", +" printf('\nThe normal stress at D = %1.2f MPa',sigma_D);\n", +" printf('\nThe normal stress at E = %1.2f MPa',sigma_E);\n", +" printf('\nThe orientation of the nuetral axis = %1.1f degrees',alpha);\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.19: B19.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.19 : ')\n", +"\n", +"//Given:\n", +"theta = 30*(%pi/180);\n", +"M = 15; //kNm\n", +"My = M*cos(theta); \n", +"Mz = M*sin(theta); \n", +"b = 0.1; //m\n", +"t1 = 0.04;//m\n", +"t2 = 0.03;//m\n", +"\n", +"\n", +"//Section Properties:\n", +"y1 = b/2;\n", +"y2 = b + t2/2;\n", +"A1 = (b*t1);\n", +"A2 = (b*2*t2);\n", +"z_bar = (y1*A1 + y2*A2)/(A1+A2);\n", +"\n", +"Iz = (1/12)*(b*t1^3) + (1/12)*(t2*(2*b)^3);\n", +"Iy = (1/12)*(t1*b^3) + b*t1*(z_bar - y1)^2 + (1/12)*(2*b*t2^3) + 2*b*t2*(y2 - z_bar)^2;\n", +"\n", +"//Maximum Bending Stress:\n", +"l_b = b+t2 - z_bar;\n", +"sigma_B = (-Mz*1000*(-b))/Iz + (My*1000*(l_b))/Iy;\n", +"sigma_B = sigma_B/10^6;\n", +"sigma_C = (-Mz*1000*(t1/2))/Iz + (My*1000*(-z_bar))/Iy;\n", +"sigma_C = sigma_C/10^6;\n", +"\n", +"sigma = max(abs(sigma_B),abs(sigma_C));\n", +"\n", +"//Orientation of the nuetral axis:\n", +"theta1 = 60*(%pi/180);\n", +"alpha = atan((Iz/Iy)*tan(theta1));\n", +"alpha = alpha*(180/%pi);\n", +"\n", +"//Display:\n", +"\n", +" \n", +" printf('\n\nThe maximum normal stress in the beam = %1.2f MPa',sigma);\n", +" printf('\n The orientation of the nuetral axis = %1.1f degrees',alpha);\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.20: B20.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.20 : ')\n", +"\n", +"//Given:\n", +"M =20; //kN\n", +"Iy = 0.96*10^-3; //m^4\n", +"Iz = 7.54*10^-3; //m^4\n", +"theta = 57.1*(%pi/180);\n", +"\n", +"\n", +"//Internal moment Components:\n", +"My = M*sin(theta); \n", +"Mz = M*cos(theta); \n", +"\n", +"//Bending Stress:\n", +"y_p = -0.2; //y Coordinate of P\n", +"z_p = 0.35; //z Coordinate of P\n", +"\n", +"theta1 = (%pi/2)-(theta);\n", +"yp = -z_p*sin(theta1)+ y_p*cos(theta1);\n", +"zp = z_p*cos(theta1) + y_p*sin(theta1);\n", +"\n", +"//Eq 6-17\n", +"\n", +"sigma_p = ((Mz*-yp)/Iz) + ((My*zp)/Iy) ;\n", +"sigma_p = sigma_p/10^3;\n", +"\n", +"//Orientation of the Nuetral Axis:\n", +"alpha = atan((Iz/Iy)*tan(theta));\n", +"alpha = alpha*(180/%pi);\n", +"\n", +"//Display:\n", +"\n", +" \n", +" printf('\n\nThe maximum normal stress at point P = %1.2f MPa',sigma_p);\n", +" printf('\nThe orientation of the nuetral axis = %1.1f degrees',alpha);\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.21: B21.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.21 : ')\n", +"\n", +"//Given:\n", +"M = 2; //kNm\n", +"Ew = 12; //GPa\n", +"Est = 200; //GPa\n", +"bw = 150/1000; //m\n", +"t = 20/1000; //m\n", +"rib = 9/1000;//m\n", +"\n", +"//Section Properties:\n", +"n = (Ew/Est);\n", +"bst = n*bw;\n", +"\n", +"y1 = t/2;\n", +"A1 = t*bw;\n", +"y2 = bw/2 + t;\n", +"A2 = rib*bw;\n", +"\n", +"y_bar = (y1*A1 +y2*A2)/(A1+A2);\n", +"\n", +"I1 = (1/12)*(bw)*(t^3) + A1*(y_bar - y1)^2;\n", +"I2 = (1/12)*(rib)*(bw^3) + A2*(y2-y_bar)^2;\n", +"Ina = I1+I2;\n", +"\n", +"//Normal Stress:\n", +"sigma_B = (M*(bw+t-y_bar))/(Ina*1000);\n", +"sigma_C = (M*(y_bar))/(Ina*1000);\n", +"\n", +"//Normal Stress in the wood:\n", +"sigmaB = n*sigma_B;\n", +"\n", +"//Display:\n", +"\n", +" \n", +" printf('\n\nThe normal stress at point B = %1.1f MPa',sigma_B);\n", +" printf('\nThe normal stress at point C = %1.2f MPa',sigma_C);\n", +" printf('\nThe normal stress at point B in the wood = %1.2f MPa',sigmaB);\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.22: B22.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.22 : ')\n", +"\n", +"//Given:\n", +"sigma_allow_st = 168; //MPa\n", +"sigma_allow_w = 21; //MPa\n", +"Est = 200; //GPa\n", +"Ew = 12; //GPa\n", +"Iz = 7.93*10^6; //mm^4\n", +"A1 = 5493.75; //mm^2\n", +"t = 5; //mm\n", +"h = 100; //mm\n", +"\n", +"//Without Board:\n", +"c = h+t;\n", +"M1 = (sigma_allow_st*Iz)/(c*10^6);\n", +"\n", +"//With Board:\n", +"bw = 300;//mm\n", +"n = (Ew/Est); \n", +"bst = n*bw;\n", +"\n", +"//For the transformed section:\n", +"y1 = 0;\n", +"y2 = 55;\n", +"A2 = bst*h;\n", +"\n", +"y_bar = (y1*A1 + y2*A2)/(A1+A2);\n", +"\n", +"I1 = Iz + A1*y_bar^2;\n", +"I2 = (1/12)*(bst*h^3) + (A2*(y2-y_bar)^2);\n", +"I = I1+I2;\n", +"\n", +"c = c+y_bar;\n", +"M2 = (sigma_allow_st*I)/(c*10^6);\n", +"\n", +"cw = c - y_bar;\n", +"Mw = (sigma_allow_w*I)/(n*cw*10^6);\n", +"\n", +"M = min(Mw,M2);\n", +"\n", +"//Display:\n", +"\n", +" printf('\n\nThe maximum bending moment without re-inforcement = %1.3f kNm',M1);\n", +" printf('\nThe maximum bending moment with re-inforcement = %1.2f kNm',M);\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.23: B23.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.23 : ')\n", +"\n", +"//Given:\n", +"M = 60; //kNm\n", +"Est = 200; //GPa\n", +"Econc = 25; //GPa\n", +"d = 25;//mm\n", +"r = d/2;\n", +"w = 300;//mm\n", +"ht =400; //mm\n", +"\n", +"//Section Properties:\n", +"n = Est/Econc;\n", +"Ast = 2*%pi*r^2;\n", +"A = n*Ast;\n", +"\n", +"p = [1 52.37 -20949.33]\n", +"h = roots(p)\n", +"h = h(2);\n", +"\n", +"I = (1/12)*(w*h^3) +w*h*(h/2)^2 + A*(ht - h)^2;\n", +"\n", +"//Normal Stress:\n", +"sigma_conc_max = (M*1000*h*1000)/(I);\n", +"sigma_conc = (M*1000*(ht-h)*1000)/(I);\n", +"sigma_st = n*sigma_conc;\n", +"\n", +"//Display:\n", +"\n", +" \n", +" printf('\n\nThe normal stress in each steel reinforcing rod = %1.2f MPa',sigma_st);\n", +" printf('\nThe maximum normal stress in the concrete = %1.2f MPa',sigma_conc_max);\n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.24: B24.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.24 : ')\n", +"\n", +"//Given:\n", +"sigma = 140; //Mpa\n", +"ri = 90; //mm\n", +"ro = 110; //mm\n", +"a = 20; //mm\n", +"\n", +"//Section Properties:\n", +" \n", +"y = integrate('a*(1/r)','r',ri,ro)\n", +"R = (a*a)/y;\n", +"\n", +"r_avg = (ri+ro)/2;\n", +"M1 = (-sigma*a*a*ro*(r_avg - R))/(R-ro);\n", +"M1 = M1*10^-6;\n", +"\n", +"M2 = (sigma*a*a*ri*(r_avg - R))/(R-ri);\n", +"M2 = M2*10^-6;\n", +"\n", +"M = min(M1,M2);\n", +"\n", +"sigma1 = (M*(R - ro))/(a*a*ro*(r_avg - R));\n", +"\n", +"//For a straight Bar:\n", +"I = (1/12)*(a*a^3);\n", +"c = 10; //mm\n", +"M_strt= (sigma*I)/c;\n", +"M_strt = M_strt*10^-6;\n", +"\n", +"//Display:\n", +" \n", +" printf('\n\nThe maximum bending moment that can be applied to the bar = %1.3f kNm',M);\n", +" printf('\nThe maximum bending moment that can be applied to a straight bar = %1.3f kNm',M_strt);\n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"\n", +"\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.25: B25.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.25 : ')\n", +"\n", +"//Given:\n", +"ri = 200/1000; //m\n", +"r1 = 250/1000; //m\n", +"ro = 280/1000; //m\n", +"M = 4; //kNm\n", +"a = 0.05; //m\n", +"h = 0.03; //m\n", +"\n", +"//Section Properties:\n", +"A1 = a^2 ;\n", +"A2 = (0.5*a*h);\n", +"A = A1+A2;\n", +"r_avg1 = (r1+ri)/2;\n", +"r_avg2 = r1+(h/3);\n", +"r_bar =((r_avg1*A1)+(r_avg2*A2))/A;\n", +"\n", +"int_dA_r1 = a*log(r1/ri);\n", +"int_dA_r2 = (a*ro*log(ro/r1))/(ro-r1) - a;\n", +"R = (A)/(int_dA_r1+ int_dA_r2);\n", +"k= r_bar - R;\n", +"\n", +"//Normal Stress:\n", +"sigma_B = (-M*(R-ri))/(A*ri*k*1000);\n", +"sigma_A = (-M*(R-ro))/(A*ro*k*1000);\n", +"\n", +"sigma = max(abs(sigma_B),abs(sigma_A))\n", +"\n", +"\n", +"//Display:\n", +" \n", +" printf('\n\nThe maximum normal stress in the bar = %1.0f MPa',sigma);\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.26: B26.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.26 : ')\n", +"\n", +"//Given:\n", +"M = 5; //kNm\n", +"sigma_y = 500; //MPa\n", +"r = 16; //mm\n", +"h = 80; //mm\n", +"w = 120; //mm\n", +"r_h = r/h;\n", +"w_h = w/h;\n", +"k = 1.45; \n", +"c = h/(2000);\n", +"t = 20/1000; //m\n", +"\n", +"//Calculations:\n", +"I = (1/12)*(t)*(h/1000)^3\n", +"sigma_max = (k*M*c)/(I*1000);\n", +"\n", +"//Display:\n", +" \n", +" printf('\n\nThe maximum normal stress in the steel = %1.0f MPa',sigma_max);\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.27: B27.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.27 : ')\n", +"\n", +"//Given:\n", +"sigma_y = 250; //MPa\n", +"t = 12.5; //mm\n", +"w = 200; //mm\n", +"h = 225; //mm\n", +"\n", +"//Maximum Elastic Moment:\n", +"yy = (h+t)/2;\n", +"I1 = (1/12)*(w*t^3) + (w*t*yy^2);\n", +"I = (1/12)*(t*h^3) + 2*(I1);\n", +"c = 125; //mm\n", +"\n", +"My = (sigma_y*I)/(c); //Flexure Formula\n", +"\n", +"//Plastic Moment:\n", +"C1= sigma_y*t*(h/2);\n", +"C2= sigma_y*t*(w);\n", +"Mp = (2*56.25*C1) + (2*yy*C2);\n", +"\n", +"//Shape Factor:\n", +"k = Mp/My;\n", +"\n", +"//Display:\n", +"\n", +" \n", +" printf('\n\nThe shape factor for the beam = %1.2f ',k);\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.28: B28.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.28 : ')\n", +"\n", +"//Given:\n", +"sigma_y = 250; //MPa\n", +"t = 15/1000; //m\n", +"w = 100/1000; //m\n", +"h = 120/1000; //m\n", +"c = 10/1000; //m\n", +"\n", +"//Calculations:\n", +"d = ((sigma_y*t*w)+(sigma_y*t*h))/(sigma_y*t*2);\n", +"\n", +"T = sigma_y*t*d*10^3;\n", +"C1 = sigma_y*t*c*10^3;\n", +"C2 = sigma_y*t*w*10^3;\n", +"\n", +"Mp = (T*d/2)+(C1*c/2)+(C2*(c+t/2));\n", +"\n", +"//Display:\n", +"\n", +" \n", +" printf('\n\nThe plastic moment that can be resisted by the beam = %1.1f kNm',Mp);\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.29: B29.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.29 : ')\n", +"\n", +"//Given:\n", +"ep1 = 0.01;\n", +"ep2 = 0.05;\n", +"sig1 = 1050;//N/mm^2\n", +"sig2 = 1330;//N/mm^2\n", +"sig3 = 280; //N/mm^2\n", +"y = 0.3; //cm\n", +"h = 3; //cm\n", +"w = 2; //cm\n", +"\n", +"//Calculations:\n", +"yy = (h/2)-y\n", +"T1 = (1/2)*(sig3*yy*w);\n", +"y1 = y +(2/3)*(yy);\n", +"T2 = yy*sig1*w;\n", +"y2 = y+(0.5*yy);\n", +"T3 = (0.5*y*sig1*w);\n", +"y3 = (2/3)*(y);\n", +"\n", +"M = 2*(T1*y1 + T2*y2 + T3*y3);\n", +"M = M/1000;\n", +"\n", +"//Display:\n", +"\n", +" \n", +" printf('\n\nThe bending moment applied that will cause a strain of 0.05mm/mm = %1.2f kNm',M);\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------\n", +"\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.30: B30.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.30 : ')\n", +"\n", +"//Given:\n", +"sigma_y = 250; //MPa\n", +"t = 12.5; //mm\n", +"w = 200; //mm\n", +"h = 225; //mm\n", +"c = (h/2)+t;\n", +"I = 82.44*10^6;//mm^4\n", +"Mp = 188; //kN\n", +"\n", +"//Calculations:\n", +"sigma_allow = (Mp*10^6*c)/(I);\n", +"y = (sigma_y*c)/(sigma_allow);\n", +"\n", +"//Display:\n", +" \n", +" printf('\n\nThe point of zero normal stress = %1.2f mm',y);\n", +" printf('\nThe Residual Stress distribution is shown in the text book.');\n", +" \n", +" //------------------------------------------------------------------------END---------------------------------------------------------------------------------------" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6.5: B5.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear all; clc;\n", +"\n", +"disp('Scilab Code Ex 6.5 : ')\n", +"\n", +"//Shear and Moment Diagrams:\n", +"p = [-1/9 -2 30]\n", +"x = roots(p)\n", +"y = (x(2));\n", +"\n", +" \n", +" M = (30*y) - (y^2) - (y^3)/27;\n", +"\n", +"\n", +"\n", +"//Display:\n", +" \n", +"printf('\n\nThe magnitude of the maximum moment is = %1.0f kNm', M);\n", +"printf('\nRefer to the shear and moment diagrams in the book.');\n", +"\n", +"\n", +"//---------------------------------------------------------------------------END-----------------------------------------------------------------------------\n", +"" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} -- cgit