blob: 0a3f795a6b7d31bff9364c239c60f1bc2c5dfcf3 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
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--------------------------------------------------------------------------
|