summaryrefslogtreecommitdiff
path: root/3864/CH7/EX7.19/Ex7_19.sce
blob: 63e01c0f39e5b42df21e4d4f328fa2501dbe9f03 (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
55
56
57
58
59
60
61
62
63
64
clear
//
//

//Initilization of Variables

d=75 //mm //diameter 
P=30*10**6 //W //Power transmitted
W=6 //N-mm/sec //Load
L=1000 //mm 
N=300 //r.p.m

//Calculations

//B.M
M=W*L*4**-1 //N-mm
T=P*60*(2*%pi*N)**-1 //Torque transmitted

//M.I
I=%pi*64**-1*d**4 //mm**4

//Bending stress
f_A=M*I**-1*(d*2**-1) //N/mm**2

//At A
p_x=f_A
p_y=0

//Polar Modulus
J=%pi*32**-1*d**4 //mm**4

//Shearing stress
q=T*J**-1*(d*2**-1) //N/mm**2

//Principal Stresses
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

//Max shear stress
q_max=(((p_x-p_y)*2**-1)**2+q**2)**0.5 //N/mm**2

//Bending stress
p_x2=0
p_y2=0

//Shearing stress
q2=T*J**-1*d*2**-1 //N/mm**2

//Principal stresses
P3=(p_x2+p_y2)*2**-1+(((p_x2-p_y2)*2**-1)**2+q2**2)**0.5 //N/mm**2
P4=(p_x2+p_y2)*2**-1-(((p_x2-p_y2)*2**-1)**2+q2**2)**0.5 //N/mm**2

//Max shear stress
q_max2=(((p_x2-p_y2)*2**-1)**2+q2**2)**0.5 //N/mm**2

//Answer for Principal Stresses P1,P2 and Max stress i.e q_max is incorrect in Book

//Result
printf("\n Principal Stresses at vertical Diameter:P1 %0.2f  N/mm**2",P1)
printf("\n                                        :P2 %0.2f  N/mm**2",P2)
printf("\n Max stress at vertical Diameter        :   %0.2f  N/mm**2",q_max)
printf("\n Principal Stresses at Horizontal Diameter:P3 %0.2f  N/mm**2",P3)
printf("\n                                        :P4 %0.2f  N/mm**2",P4)
printf("\n Max stress at Horizontal Diameter      :   %0.2f  N/mm**2",q_max2)