blob: 64726a76c17cc39eac0ee920a1ccc0f186a2d35f (
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
|
//Function to round-up a value such that it is divisible by 5
function[v] = round_five(w)
v = ceil(w)
rem = pmodulo(v,5)
if (rem ~= 0)
v = v + (5 - rem)
end
endfunction
//Obtain path of solution file
path = get_absolute_file_path('solution4_5.sce')
//Obtain path of data file
datapath = path + filesep() + 'data4_5.sci'
//Clear all
clc
//Execute the data file
exec(datapath)
//Refer Fig.4.24
//Permissible stresses (N/mm2)
//Tensile stress sigmat
sigmat = Syt/fs
//Compressive stress
sigmac = Syt/fs
//Yield strength in shear for the material Ssy (N/mm2)
Ssy = (50/100)*Syt
//Shear stress tau
tau = Ssy/fs
//Calculate diameter of the rods D (mm)
D = sqrt((4 * P)/(%pi * sigmat))
D = round_five(D)
//Calculate enlarged diameter of rods D1 (mm)
D1 = 1.1 * D
//Calculate dimensions a and b (mm)
a = 0.75 * D
a = round_five(a)
b = 1.25 * D
b = round_five(b)
//Calculate diameter of pin d
d = ((32/(%pi * sigmat))*((P/2)*((b/4) + (a/3))))^(1/3)
d = round_five(d)
//Calculate dimensions d0 and d1
d0 = 2 * d
d1 = 1.5 * d
//Print results
printf('\nDiameter of rods(D) = %f mm\n',D)
printf('\nEnlarged diameter of rods(D1) = %f mm\n',D1)
printf('\nDimensions of a and b are %f mm and %f mm respectively\n',a,b)
printf('\nDiameter of pin(d) = %f mm\n',d)
printf('\nDimensions of d0 and d1 are %f mm and %f mm respectively\n',d0,d1)
//Check for stresses in eye
s1 = (P/(b * (d0 - d)))
c1 = (P/(b * d))
t1 = (P/(b * (d0 - d)))
//Check for stresses in fork
s2 = (P/(2 * a * (d0 - d)))
c2 = (P/(2 * a * d))
t2 = (P/(2 * a * (d0 - d)))
if ((s1<sigmat) & (c1<sigmac) & (t1<tau) & (s2<sigmat) & (c2<sigmac) & (t2<tau))
printf('\nDesign of knuckle joint is safe\n')
else
printf('\nDesign of knuckle joint is not safe\n')
end
|
