blob: 22a9b104bad09a269cfafa5760d50422a0dae800 (
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
|
clear all; clc;
disp("Scilab Code Ex 6.16 : ")
//Given:
t1 = 15/1000; //m
t2 = 20/1000; //m
l = 250/1000; //m
b = 200/1000; //m
P = 2.4; //kN
l_a = 2; //m
l_b = 1; //m
//Internal Moment:
y1 = b/2;
y2 = t2/2;
A = (2*t1*b)+(t2*l);
y_bar = ((2*y1*t1*b)+(y2*t2*l))/A;
M = (P*l_a)+(1*y_bar);
//Section Property:
I1 = (1/12)*(l*t2^3) + (l*t2*(y_bar - y2)^2);
I2 = (1/12)*(t1*b^3) + (t1*b*(y1 - y_bar)^2);
I =I1+ 2*I2;
//Maximum Bending Stress:
c = b - y_bar;
sigma_max = (M*c)/(I*1000);
//Display:
printf('\n\nThe maximum bending stress at section a-a = %1.1f MPa',sigma_max);
//-----------------------------------------------------------------END--------------------------------------------------------------------------
|
