summaryrefslogtreecommitdiff
path: root/3864/CH4/EX4.17/Ex4_17.sce
blob: 0477cf682551d68dcab3290ddecc01a8643b5bc4 (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
45
46
47
48
49
50
51
52
53
54
clear
//
//

//Initilization of Variables

F=40*10**3 //N //shear Force

//I-section

//Flanges
b=80 //mm //Width of flange
t=20 //mm //Thickness

//Web
d=200 //mm //Depth
t2=20 //mm //Thickness

//Flange-2
b2=160 //mm //Width
t3=20 //mm //Thickness

D=240 //mm //Overall Depth

//Calculations

//Distance of N-A from Top Fibre 
y=(b*t*t*2**-1+d*t2*(t+d*2**-1)+b2*t3*(t+d+t3*2**-1))*(b*t+d*t2+b2*t3)**-1 //mm

//Moment of Inertia
I=1*12**-1*b*t**3+b*t*(y-(t*2**-1))**2+1*12**-1*t2*d**3+t2*d*(y-(t+d*2**-1))**2+1*12**-1*b2*t3**3+t3*b2*((d+t+t3*2**-1)-y)**2 //mm**4

//Shear stress bottom of flange
sigma=F*b*t*(y-t*2**-1)*(b*I)**-1 //N/mm**2

//At same Level but in web
sigma2=F*b*t*(y-t*2**-1)*(t2*I)**-1 //N/mm**2

//for shear stress at N.A
X=b*t*(y-t*2**-1)+t2*(y-t)*(y-t)*2**-1 //mm**3
sigma3=F*X*(t2*I)**-1 //N/mm**2

//Shear stress at bottom of web

X=b2*t3*((D-y)-t3*2**-1) //mm**3

//Stress at bottom of web
sigma4=F*X*(t2*I)**-1 //N/mm**2

//Stress at Lower flange
sigma5=F*X*(b2*I)**-1 //N/mm**2

//Result
printf("\n The Shear Force Diagram is the result")