summaryrefslogtreecommitdiff
path: root/1898/CH13/EX13.2/Ex13_2.sce
blob: 12909c52c82dc907d1c08d1886297eec593441e4 (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
clear all; clc;

disp("Scilab Code Ex 13.2 : ")

//Given:
E = 200; //GPa
I = 15.3*10^6; //mm^4
l= 4*1000; //mm
A = 5890; //mm^2
sigma_y = 250; //MPa

//Calculations:

Pcr = ((%pi^2)*E*10^6*I*1000^-2)/(l^2); //Pcr = (%pi^2*EI)/(l^2)

sigma_cr = (Pcr*1000)/A;

if(sigma_cr>sigma_y)
    Pcr = (sigma_y*A);
    Pcr = Pcr/1000; //in kN
end

//Display:

printf("\n\nThe maximum allowable axial load that the column can support    = %1.1f kN',Pcr);

//----------------------------------------------------------------------END--------------------------------------------------------------------------------