blob: 981db26e1323ba75b67d18edc28087eaf2589d1c (
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
|
clc
// initialization
clear
d=10 //cm
id=9.99 //cm
t=3 //mm
E=1.0*10^6 //kg/cm^2
a=2.02*10^-5 // degree/celcius
// part(a)
Tr=10 //degree C
T=(d-id)/id*1/a
printf('part(a) \n The sleeve must be heated to %.1f degree C or more for this purpose',T+Tr)
//part(b)
s_th=a*T*E
p=s_th*t*2/(d*10)
printf('\n part(b) \n The pressure developed between the rod and sleeve is %d kg/cm^2',p)
// part(c)
f=0.2
o=10 // overlap: cm
A=%pi*d*o
F=f*p*A
printf('\n part (c) \n The axial force required is %d kg',F)
//part (d)
// linked to part c
T2=20 //degree C
a2=1.17*10^-5 // /degree C
Ts=(a-a2)*(T2-Tr)*E
Ts=s_th-Ts
p2=p*Ts/s_th
F2=F*Ts/s_th
printf('\n part(d)\n The pressure developed between the rod and sleeve is %.1f kg/cm^2',p2)
printf('\n The axial force required is %d kg',F2)
//part(e)
T3=Tr+(s_th/((a-a2)*10^6))
printf('\n part(e) \n The temperature at which the sleeve comes off easily is %.1f C',T3)
// calculations in the text: rounding off errors
|
