summaryrefslogtreecommitdiff
path: root/764/CH4/EX4.6.b/solution4_6.sce
blob: e25066af71c2af05984096d67bb5d8fe7332a9bf (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
43
44

//Function to print the designation
function[] = desig(z)
    printf('\nISLB%d is suitable for this application\n',z)
endfunction
 

//Obtain path of solution file
path = get_absolute_file_path('solution4_6.sce')
//Obtain path of data file
datapath = path + filesep() + 'data4_6.sci'
//Clear all
clc
//Execute the data file
exec(datapath)
//Load supported by each beam P (kN)
P = (W/2)
//Distance of centre of gravity of the load from the rigid support d (mm)
d  = (l/3)
//Calculate the bending moment Mb (N-mm)
Mb = (P * 1000)*d
//Calculation of ratio Ixx/y (mm3)
ratio = Mb/sigmab
//Selection of beam
//There are 5 standard beams provided 
for i = 1:1:5 
    x = Ixx(i)/(h(i)/2)
//Break the loop if value of x exceeds value of ratio
    if (x > ratio)
        break
    end
end
//Print the suitable beam designation obtained
if (i == 1)
    desig(150)
elseif (i == 2)
    desig(175)
elseif (i == 3)
    desig(200)
elseif (i == 4)
    desig(250)
else
    desig(250)
end