summaryrefslogtreecommitdiff
path: root/3864/CH7/EX7.14/Ex7_14.sce
blob: a000a0993e06d5b13014fd255bb8a5e9f0e596f6 (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
clear
//
//

//Initilization of Variables

F=40*10**3 //N //Shear Force
M=20*10**6 //Bending Moment

//Rectangular section
b=100 //mm //Width
d=200 //mm //Depth

x=20 //mm //Distance from Top surface upto point
y=80 //mm //Distance from point to Bottom

//Calculations

I=1*12**-1*b*d**3 //mm**4 //M.I

//At 20 mm Below top Fibre
f_x=M*I**-1*y //N/mm**2 //Stress

//Assuming sagging moment ,f_x is compressive p_x=f_x=-24 //N/mm**2
f_x=-24 //N/mm**2
p_x=-24

//Shearing stress
q=F*(b*I)**-1*(b*x*(b-x*2**-1)) //N/mm**2

//Direct stresses

p_y=0 //N/mm**2

p1=(p_x+p_y)*2**-1+(((p_x+p_y)*2**-1)**2+q**2)**0.5 //N/mm**2
p2=(p_x+p_y)*2**-1-(((p_x+p_y)*2**-1)**2+q**2)**0.5 //N/mm**2

//Result
printf("\n Directions of principal stresses at a point below 20mm is: %0.2f  N/mm**2",p1)
printf("\n                                                            %0.2f  N/mm**2",p2)