blob: 12ff51aff8c5ebb866a40345f02097885a4bb839 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
clc
//initialisation of variables
load= 50000 //N
torque= 1200 //Nm
d= 60//mm
t= 1.5 //mm
alpha= 60 //degrees
//CALCULATIONS
BM= load*t
axialload= (load*4)/(%pi*d^2)
bendingmoment= (BM*d*64)/(%pi*d^4*2)
Ts= axialload+bendingmoment
shearstress= (torque*10^3*d*32)/(2*%pi*d^4)
sigman= -Ts*(cosd(alpha-30))^2+shearstress*cosd(alpha-30)*sind(alpha-30)+shearstress*cosd(alpha-30)*sind(alpha-30)
T= -Ts*sind(alpha)*cosd(alpha)-shearstress*(sind(alpha))^2+shearstress*(cosd(alpha)^2)
//RESULTS
printf ('direct stress= %.1f N/mm^2',sigman)
printf (' \n Shear stress=%.1f N/mm^2',T)
|
