blob: 39e1aa59864046e6cfaef66e485e8ac3d93b87f4 (
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 13.13 : ")
//Given:
sigmaB_allow = 160; //MPa
E = 200; //GPa
sigma_y = 250; //MPa
K= 1;
A = 3790; //mm^2
Ix = 17.1*10^6; //mm^4
ry = 38.2; //mm
d = 157; //mm
c= d/2;
e = 750; //mm
L = 4000; //mm
sl_ratio = (K*L)/(ry);
sl_ratio_c = sqrt((2*%pi^2*E*1000)/(sigma_y));
if(sl_ratio<sl_ratio_c)
num = (1 - (sl_ratio^2/(2*sl_ratio_c^2)))*sigma_y;
denom1 = (5/3) + ((3/8)*sl_ratio/sl_ratio_c);
denom2 = (sl_ratio^3)/(8*sl_ratio_c^3);
sigmaA_allow = num/(denom1 - denom2);
coeffP = 1/(sigmaA_allow*A) + (e*c)/(Ix*sigmaB_allow);
P = 1/coeffP;
sigA = (P/A)/(sigmaA_allow);
P = P/1000; //in kN
if(sigA < 0.15)
printf('\n\nThe maximum allowable value of eccentric load = %1.2f kN',P);
end
end
//---------------------------------------------------------------------------END--------------------------------------------------------------------------------
|
