summaryrefslogtreecommitdiff
path: root/1898/CH4/EX4.11/Ex4_11.sce
blob: 538d71672c39f03999550b6adb77aa14eed64446 (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
clear all; clc;

disp("Scilab Code Ex 4.11 : ")

//Given:
area_sleeve = 600*10^-6; //m^2
area_bolt = 400*10^-6; //m^2
T1 = 15; //degree celcius
T2 = 80; //degree celcius
alpha_bolt = 12*10^-6; //per degree celcius
alpha_sleeve = 23*10^-6; //per degree celcius
l = 0.15; //m
E_bolt = 200*10^9; //N/m^2 
E_sleeve = 73.1*10^9; //N/m^2 

//Equilibrium:
//F_s = F_b

//Compatibility:
del_T = T2 - T1; // temperature difference
delb_T = alpha_bolt*del_T*l; 
delb_F = l/(area_bolt*E_bolt);
dels_T = alpha_sleeve*del_T*l; 
dels_F = l/(area_sleeve*E_sleeve);

//delb_T + F_b*delb_F = dels_T + F_s*dels_F

F_b = (dels_T-delb_T)/(delb_F+dels_F);
F_b = F_b/1000; //in kN
F_s= F_b;

sigma_b = F_b/(area_bolt*10^3); //Average Normal Stress
sigma_s = F_s/(area_sleeve*10^3); //Average Normal Stress

//Display:


printf("\n\nThe force experienced by sleeve and bolt       = %1.2f kN",F_s);
printf('\nThe average normal stress on bolt              = %1.1f MPa',sigma_b);
printf('\nThe average normal stress on sleeve            = %1.1f MPa',sigma_s);


//-----------------------------------------------------------END-----------------------------------------------------------------------------------------------