blob: 09221472ef5b8bd313fd9efee91bc476f365f305 (
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
|
function[] = plot_format()
//Get the handle of current axes
g = gca()
//Give labels and set label properties
g.labels_font_color=5
g.font_size=3
g.grid=[1,1]
g.box="off"
endfunction
//Obtain path of solution file
path = get_absolute_file_path('solution7_18.sce')
//Obtain path of data file
datapath = path + filesep() + 'data7_18.sci'
//Obtain path of function file
funcpath = path + filesep() + 'functions7_18.sci'
//Clear all
clc
//Execute the data file
exec(datapath)
exec(funcpath,[-1])
//Calculate the endurance limit stress for bolt Sdash (N/mm2)
Sdash = (50/100)*Sut
//Calculate Ka, Kb and Kc
[Ka, Kb, Kc] = fluctuate(op, d, reliability)
//Calculate Kd
Kd = 1/Kf
//Calculate the corrected endurance limit stress Se (N/mm2)
Se = Kb * Kc * Kd * Sdash
//Calculate the stiffness of bolt kb (N/mm)
kb = ((%pi/4)*(d^2))*(E1/l)
//Calculate the area of the two plates Ac (mm2)
Ac = (%pi/4)*(((2*d)^2) - (d^2))
//Calculate the stiffness of the plates kc (N/mm)
kc = (Ac * E2)/l
//Calculate the maximum force in the bolt PMAX (N)
PMAX = (Pi * 1000) + ((kb/(kb + kc))*(Pmax * 1000))
//Calculate the minimum force in the bolt PMIN (N)
PMIN = (Pi * 1000) + ((kb/(kb + kc))*(Pmin * 1000))
//Calculate the mean force and force amplitude
Pm = (PMAX + PMIN)/2
Pa = (PMAX - PMIN)/2
//Plot modified Goodman diagram
//The common quadrilateral in the plot is the area of concern
y1 = {Se 0}
x1 = {0 Sut}
y2 = {Syt 0}
x2 = {0 Syt}
plot(x1,y1,'--*')
plot(x2,y2,'-*')
plot_format()
title('Modified Goodman diagram (Example 7.18)')
xlabel('sigmaM (N/mm2)')
ylabel('sigmaA (N/mm2)')
//Calculate the actual core cross-section area of the bolt A (mm2)
A = (Pa + ((Pi * 1000)/((1 + (Sut/Se)) * fs)))/(Sut/((1 + (Sut/Se))*fs))
//Choose proper diameter from Table 7.1
//Print results
printf('\nCore cross-section area of the bolt(A) = %f mm2\n',A)
|
