From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 1898/CH6/EX6.17/Ex6_17.sce | 42 ++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 42 insertions(+) create mode 100755 1898/CH6/EX6.17/Ex6_17.sce (limited to '1898/CH6/EX6.17/Ex6_17.sce') diff --git a/1898/CH6/EX6.17/Ex6_17.sce b/1898/CH6/EX6.17/Ex6_17.sce new file mode 100755 index 000000000..0a3f795a6 --- /dev/null +++ b/1898/CH6/EX6.17/Ex6_17.sce @@ -0,0 +1,42 @@ +clear all; clc; + +disp("Scilab Code Ex 6.17 : ") + +//Given: +b = 60/1000; //m +h = 30/1000; //m +M = 40; //Nm +c1= h/2; +rib_t = 5/1000; //m +rib_w = 10/1000;//m + +//Without Ribs: +I1 = (1/12)*(b*h^3); +sigma_max1 = (M*c1)/(I1*10^6); + +//With Ribs: +y1 = c1; +y2 = h+(rib_t/2); +A1 = h*b; +A2 = rib_t*rib_w; +y_bar = ((y1*A1)+2*(y2*A2))/(A1 + 2*A2); + +c2 = h+rib_t - y_bar; +I2 = I1 + (b*h*(y_bar - y1)^2); +I3 = (1/12)*rib_w*rib_t^3 + (rib_w*rib_t*(y2 - y_bar)^2); +I = I2 + 2*I3; + +sigma_max2 = (M*c2)/(I*10^6); + +if(sigma_max2>sigma_max1) + + printf("\n\nThe maximum normal stress in the member without ribs = %1.2f MPa',sigma_max1); + printf("\nThe maximum normal stress in the member with ribs = %1.2f MPa',sigma_max2); + printf("\nThe ribs should be omitted."); + + end + + +//-----------------------------------------------------------------END-------------------------------------------------------------------------- + + -- cgit