Research Article
A new computational method for threaded connection stiffness
Advances in Mechanical Engineering 2016, Vol. 8(12) 1–9 Ó The Author(s) 2016 DOI: 10.1177/1687814016682653 aime.sagepub.com
Dongmei Zhang, Shiqiao Gao and Xiao Xu
Abstract For designing bolted connections in machinery applications, it is necessary to estimate the stiffness of the threaded connection. This work provides a new method for computing the stiffness of engaged screw in bolted connections according to the load distribution in screw thread. Finite element analysis is performed by building the three-dimensional model of threaded connection. A set of tensile tests are exerted to validate the accuracy of the suggested model of threaded connection. A good agreement is obtained when the analytical results are compared with finite element analysis results, experimental data, and Yamamoto method. Results reveal that the ultimate strength of thread connections is obviously lower than that of thread material. In addition, the results of calculation and finite element analysis indicated that increasing Young’s modulus of material and the engaged length or decreasing thread pitch could increase the stiffness of the thread portion of a bolt and nut. Keywords Threaded connection, stiffness, finite element analysis
Date received: 5 September 2016; accepted: 12 November 2016 Academic Editor: Yangmin Li
Introduction Threaded fasteners, mainly bolts and nuts, have played a significant role in the advancement of the design and analysis of joints and connections. In a joint, the increase in bolt tension due to external load application depends on the bolt and joint stiffness. Therefore, the correct determination of the joint stiffness is very important in the design of bolted connections. The joint stiffness mainly depends on both the stiffness of the thread connections and the effective stiffness of the clamped member material. Many researchers such as Maruyama et al.1 and Motash2 used different techniques, mainly numerical, to investigate the stiffness of the joint and the effect of different design parameters. Wileman et al.3 performed a two-dimensional (2D) finite element analysis (FEA) for joints having a range of materials and geometries. Shigley et al.4 presented an analytical solution for member stiffness, based on the work of Lehnhoff and Wistehuff.5 Naser,6
Musto and Konkle,7 Nawras,8 and Naser and Abbound9 also proposed mathematical expressions for the bolted-joint stiffness. In addition, Qin et al.10 and Liu et al.11 researched the dynamic behavior of bolted joints. Nearly all of these studies are limited to estimate the stiffness of bolted joint as a fixed unit without any consideration of engaged threaded deformation. In order to obtain the precise stiffness of a threaded connection, several authors studied the load distribution and deflection of the engaged threads. Among the analytical methods, the Sopwith method12 and the Yamamoto method13 received extensive recognition. State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing, P.R. China Corresponding author: Dongmei Zhang, State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, P.R. China. Email:
[email protected]
Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 3.0 License (http://www.creativecommons.org/licenses/by/3.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/ open-access-at-sage).
2
Advances in Mechanical Engineering
Figure 1. Force on threaded connections.
The Yamamoto method not only proposed a procedure for calculating the axial elastic deflection of thread but also provided the engaged screw stiffness and its calculation method. The assumption for Yamamoto method is that the load per unit width along the helix direction is uniformly distributed. In fact, for the threedimensional (3D) helix thread, the load distribution is not uniform. Therefore, we propose a new method which can compute the engaged screw stiffness more accurately by considering the load distribution. Meanwhile, the FEA and the tensile test for thread connections are conducted in order to verify the accuracy of the method proposed in this article. Results are compared with those obtained by Yamamoto method.
The engaged screw stiffness The thread is a critical component in connecting the bolt body and the nut body. The interaction between the bolt body and the nut body is transferred by thread body as shown in Figure 1.
The elastic deflection of thread According to the literature,13 for the ISO metric triangle thread under unit force per unit width, the axial elastic deformation of thread includes the deflection d1 due to bending moment, the deflection d2 due to shear loading, the deflection d3 due to incline at the root of the thread, the deflection d4 due to shear at the root of the thread, and the deflection d5 due to radial contraction and expansion on the bolt and the nut. For the ISO thread under unit force at x = l1 as shown in Figure 2, the axial elastic deflections of unit length thread of the bolt are
Figure 2. Geometrical parameters of ISO thread.
1 y2b d1b = 0:034 Eb 1 + yb d2b = 1:08 Eb 1 y2b d3b = 0:229 Eb 1 y2b d4b = 1:18 Eb 1 D d5b = 0:056(1 yb ) Eb p
ð1Þ
The axial elastic deflections of unit length thread of nut are 1 y2n d1n = 0:073 En 1 yn d2n = 1:15 En 1 y2n d3n = 0:294 En 1 y2n d4n = 1:14 En 2 D0 + D2 1 D + yn d5n = 0:056 2 2 En p D0 D
ð2Þ
where the sub-indices b and n represent the bolt and the nut. E is Young’s modulus, y is Poisson ratio, p is pitch, and D is the pitch diameter of screw as shown in Figure 3. D0 is the outer diameter of nut.
Zhang et al.
3 For the engaged bolt and nut, the axial deflection is Dy = Db + Dn =
1 ∂F 1 ∂F 1 1 ∂F + = + kby ∂y kny ∂y kby kny ∂y ð13Þ
The stiffness of unit axial length is ky =
Figure 3. Diagram of lead angle.
Therefore, for the thread on bolt and nut, the elastic deflection caused by unit force per unit width along axial direction is db = d1b + d2b + d3b + d4b + d5b
ð3Þ
dn = d1n + d2n + d3n + d4n + d5n
ð4Þ
For the threaded connections under compressive force F as shown in Figure 1, the axial deflection of the bolt and the nut can be expressed as Db = db
∂F ∂s
ð5Þ
Dn = dn
∂F ∂s
ð6Þ
where s is the thread length along the direction of the helix, and s=
y sin b
ð8Þ
Substituting equation (8) into equations (5) and (6) leads to Db = db sin b
∂F 1 ∂F = ∂y kby ∂y
ð9Þ
Dn = dn sin b
∂F 1 ∂F = ∂y kny ∂y
ð10Þ
where kby and kny represent the stiffness of unit axial length for thread on bolt and nut, respectively, and they can be expressed as ∂F=∂y 1 = Db db sin b
ð11Þ
∂F=∂y 1 = kny = Dn dn sin b
ð12Þ
kby =
ð14Þ
Axial load distribution in threaded connections Figure 1 shows the threaded connection structure including a fixed nut body and a bolt body. For the bolt body under compressive force F, the force will be transferred to the nut body through the threads. Because the structural deformation is presumed elastic, the bolt body, the nut body, and the thread body will exhibit deflection under force. However, the forces on the bolt body and the nut body differ from those on the thread body. Furthermore, their load distributions are also varied. Therefore, deformation modes of both the bolt body and the nut body will be different from those of the thread body. However, the bolt body and the nut body are joined by the thread body, their deformations should be compatible. The region subjected to compressive force on the thread body ranges from 0 to l, where l is the engaged length of the bolted joint. Thus, strain eb generated by the bolt body under compressive force F(y) along y-direction is
ð7Þ
where y is the axial length of thread and b is the lead angle as shown in Figure 3. Therefore ∂F ∂F ∂y ∂F = = sin b ∂s ∂y ∂s ∂y
∂F=∂y 1 = Dy (db + dn ) sin b
eb =
F(y) Sb Eb
ð15Þ
where Sb is the cross-sectional area of the bolt body and Eb is Young’s modulus of the bolt body. Similarly, strain en generated by the nut body under compressive force F(y) is en =
F(y) Sn En
ð16Þ
where Sn is the cross-sectional area of the nut body and En is Young’s modulus of the nut body. According to equations (9) and (10), the displacement gradients are respectively ∂Db 1 ∂2 F = ∂y kby ∂y2
ð17Þ
∂Dn 1 ∂2 F = kny ∂y2 ∂y
ð18Þ
In order to satisfy the deformation compatibility condition, eb , en , Db , and Dn need to meet the following relationship
4
Advances in Mechanical Engineering
eb + en =
∂Db ∂Dn + ∂y ∂y
ð19Þ
Ðl where 0 u(y)kby dy is the total thread stiffness of nut Kb . According to equation (28), Kb can be expressed as
Substituting equations (15)–(18) into equation (19) leads to l2 F = where l = (20) is
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
(1=Sb Eb ) + (1=Sn En ) (1=kby ) + (1=kny ) .
∂2 F ∂y2
The solution of equation
ð21Þ
By the boundary conditions F(y = 0) = FB and and F(y = l) = 0, we can obtain C2 = FB ll F . Substituting C and C into equation C1 = cosh B 1 2 sinh ll (21), the axial load of threaded connections can be expressed as cosh ll sinh ly + cosh ly F(y) = FB sinh ll
u(y)kby dy =
FB Db (0)
ð29Þ
ð22Þ
Total stiffness of threaded connections
ð23Þ
The deflection of thread should also satisfy the distribution law because the deflections are directly proportional to the loads. Therefore, axial deflection of the bolt and nut can be expressed as Db (y) = Db (0)u(y)
ð24Þ
Dn (y) = Dn (0)u(y)
ð25Þ
where D(0) is the axial deflection of screw at y = 0. Substituting equations (24) and (25) into equations (9) and (10) leads to Db (0)u(y)kby =
∂F ∂y
ð26Þ
Dn (0)u(y)kny =
∂F ∂y
ð27Þ
For the bolt, because the load is FB at y = 0, equation (26) can be changed into ðl Db (0) u(y)kby dy = FB
kby cosh ll cosh ll cosh ll + sinh ll + Kb = l sinh ll sinh ll cosh ll 1 ð30Þ = Eb Sb l sinh ll Similarly, the thread stiffness of nut can be expressed as kny cosh ll cosh ll cosh ll + sinh ll + l sinh ll sinh ll cosh ll 1 ð31Þ = En Sn l sinh ll
Because the meshing bolt and nut is in series manner, the total stiffness of engaged thread is K=
According to the above deflections of the thread, the law of axial load distribution can be expressed as F(y) cosh ll sinh ly + cosh ly = FB sinh ll
Substituting equation (23) into equation (29), the thread stiffness of bolt can be expressed as
Kn =
where FB is the load on the first thread of the bolt as shown in Figure 1.
0
Kb = 0
ð20Þ
F(y) = C1 sinh ly + C2 cosh ly
u(y) =
ðl
ð28Þ
1 Kb
1 +
1 Kn
sinh ll 1 1 = + ð32Þ l( cosh ll 1) Eb Sb En Sn
Finite element model Finite element models are created to study the stiffness of threaded connections and the influence of material, pitch, and engaged length on the stiffness. The FEA method is used to calculate the bolt deflection in the axial direction, which is taken as the average nodal displacement daverage at the top surface of bolt. The parameters of all the screw are listed in Table 2. Commercial finite element software ANSYS is used for modeling and analyzing the engaged bolt-nut threads. The model geometry is meshed by 8-node hexahedron elements (SOLID185). Contact and target elements are TARGE169 and CONTA172, respectively. Figure 4 shows the 3D finite element models of the standard bolt assembly. In the models, the axial loading is applied on the top surface of bolt, and the outer bottom surface of nut is assumed fixed. Convergence study is carried out on the initial finite element model by decreasing the element size near the threads. The smallest element size was 0.25 mm by 0.25 mm and there was no significant improvement in accuracy using the smaller elements. The stress distribution on the screw under external load is shown in Figure 5. Elastic material is used throughout this work. The bolt and nut materials used in this study are the same. Uniform pressure loading P is applied on the top
Zhang et al.
5
Figure 4. 3D model mesh of standard bolted joint assembly (left) and bolt itself (right).
surface of bolt. The stiffness of the thread connections KF is calculated using the following simple relation KF =
P daverage
ð33Þ
Tensile tests
Figure 5. Contour plot of the screw stress.
In order to obtain the precise values of material parameter, we conducted the tensile tests for the standard tensile bar as shown in Figure 6. The universal test machine (UTM) is used to measure the load–deflection data of samples. The tested samples are made from steel, copper, and brass material, respectively. For the thread connection samples as shown in Figure 7, the test machine can only obtain the load f. Because the nut is fixed, the axial deflection of thread
Figure 6. Tensile tests for standard bar: (a) tensile test and (b) standard bar (unit: mm).
6
Advances in Mechanical Engineering
Figure 7. Tensile tests for thread connection samples: (a) thread connection samples and (b) tensile test.
connection can be represented by the displacement variation DL between two lines as shown in Figure 7(b), which can be measured by a video gauge. The slope of f DL curve represents the joint stiffness. For trustworthy test data, 10 tests were performed for each specimen size and the average data for each case is considered.
Results and discussion The stress–strain curves for tensile standard bar are shown in Figure 8. Figure 8(a) is the whole tensile curve, by which we not only can obtain the ultimate strength of material but also can calculate Young’s modulus. In order to distinguish clearly the elastic part of tensile curve, the partially enlarged curve for Figure 8(a) is shown in Figure 8(b). Therefore, the slope of the initial linear section in the curves is Young’s modulus. By the same method, the ultimate strengths of thread connections are obtained. Table 1 presents the parameter of materials as well as the ultimate strength of material and thread connections. Comparing the
Figure 8. Stress–strain curves for tensile standard bar: (a) the whole curves and (b) the partially enlarged curves.
ultimate strength, we can know that the values of thread connections are less than those of material. This reveals that the strength of thread connections is clearly lower than fixed connections. Table 2 lists the stiffness of thread connections from experiments, FEA, the current computational method, and Yamamoto method. From this table, it is very clear that the stiffness from this study is lower than Yamamoto method because we considered the load distribution along helix. Therefore, the stiffness from this study is in good agreement with the results from experiment and FEA. An alternative way to present and compare these results is executed in Figures 9–11. Figures 9–11 are the force-deflection curves for different materials, engaged lengths and pitches, respectively. In Figure 9, all the specimens 1#, 2#, and 3# have the pitch p 4 mm, the engaged length 20 mm, and the steel material, copper,
Zhang et al.
7
Table 1. Parameters of samples. No.
1# 2# 3# 4# 5# 6# 7#
Size code of threads
M3634 20 M3634 20 M3634 20 M3634 32 M3634 12 M3633 12 M3632 12
Material
Steel Copper Brass
Parameters for material Young’s modulus (GPa)
Poisson ratio
Ultimate strength (MPa)
211 138 107
0.269 0.31 0.32
845 396 313
Figure 9. Force–deflection curves for different materials: (a) the whole curves and (b) the partially enlarged curves.
and brass, respectively. All these curves indicate that the higher the Young’s modulus of material for thread connections, the higher the stiffness. But the ultimate
Ultimate strength for thread connections (MPa)
324 149 240 328 147 132 91
Figure 10. Force–deflection curves for different engaged length: (a) the whole curves and (b) the partially enlarged curves.
strength of thread connections with brass is higher than copper due to the material properties.
8
Advances in Mechanical Engineering
Table 2. Stiffness of thread connections (109 N=m). No.
1# 2# 3# 4# 5# 6# 7#
Size code of threads
M3634 20 M3634 20 M3634 20 M3634 32 M3634 12 M3633 12 M3632 12
Exp.
5.2279 3.3964 2.6643 3.6276 1.8013 2.0853 2.3962
FEA
Theory
5.4258 3.5856 2.7616 3.6309 1.8168 2.0949 2.4080
This study
Yamamoto’s method
5.5638 3.6272 2.7765 3.6378 1.8464 2.1062 2.4332
7.9704 5.1536 3.9461 4.2821 3.1293 3.5251 4.0089
FEA: finite element analysis.
behavior can be explained by equation (32). From this equation, we can know that the total stiffness obviously increases with the length of screw when the engaged length is small, and the value will not clearly change when the engaged length is large. In Figure 11, the pitches of these specimens are 4, 3, and 2 mm. It states clearly that the stiffness of thread connections structure with the smallest pitch is the greatest, but the ultimate strength is the smallest. The reason for this behavior can be explained by equation (14). From the equation, the stiffness of screw increases with the decrease in lead angle. Because the lead angle is small for the screw with small pitch, the stiffness of screw with smaller pitch is higher. However, the component along the axial direction is small when the lead angle of screw is small. Therefore, the screw with smaller pitch has the higher ultimate strength.
Conclusion A new computational method for threaded connection stiffness was developed in this study. Effects of some geometrical parameters of the engaged threads on the estimated stiffness were investigated for three different commonly used materials. The results can be summarized as follows: 1.
2. Figure 11. Force–deflection curves for different pitches: (a) the whole curves and (b) the partially enlarged curves.
In Figure 10, material of all the specimens is brass. They have the same pitch, but the engaged lengths are 32, 20, and 12 mm. It is clear that the longer the engaged length for thread connections, the higher the stiffness and the ultimate strength. The reason for this
3.
The ultimate strength of threaded connections is lower than the corresponding material itself, which indicates that the strength of thread connections is obviously lower than fixed connection. The computed threaded stiffness obtained in this study showed good agreements when compared with FEA results and experimental data. Results obtained using Yamamoto method were larger than those obtained using the present computational model, experiment, and FEA. The stiffness was highly related to the material’s elastic properties, the engaged length, and pitch of thread. We could obtain higher stiffness through increasing Young’s modulus of the
Zhang et al. used material and the engaged length or by decreasing the pitch. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding The author(s) received no financial support for the research, authorship, and/or publication of this article.
References 1. Maruyama K, Yoshimoto I and Nakano Y. On spring constant of connected parts in bolted joints. Bull Jpn Soc Mech Eng 1975; 18: 1472–1480. 2. Motash N. Determination of joint stiffness in bolted connections. J Eng Ind: T ASME 1976; 98: 858–861. 3. Wileman J, Choudhury M and Green I. Computation of member stiffness in bolted connections. J Mech Design 1991; 113: 432–437. 4. Shigley JE, Mischke CR and Budynas RG. Mechanical engineering design. New York: McGraw-Hill, 2004.
9 5. Lehnhoff TF and Wistehuff WE. Nonlinear effects on the stresses and deformations of bolted joints. J Mech Design 1996; 118: 54–58. 6. Naser SA. Computation of member stiffness in bolted connections using the finite element analysis. Mech Base Des Struct Mach 2005; 33: 331–342. 7. Musto JC and Konkle NR. Computation of member stiffness in the design of bolted joints. J Mech Design 2006; 128: 1357–1360. 8. Nawras H. Mathematical representation of bolted-joint stiffness: a new suggested model. J Mech Sci Technol 2011; 25: 2827–2834. 9. Nassar SA and Abbound A. An improved stiffness model for bolted joints. J Mech Design 2009; 131: 121001. 10. Qin ZY, Hai QK and Chu FL. Analytical model of bolted disk-drum joints and its application to dynamic analysis of joint rotor. Proc IME C J Mech Eng Sci 2014; 228: 646–663. 11. Liu JH, Ouyang HJ, Peng JF, et al. Experimental and numerical studies of bolted joints subjected to axial excitation. Wear 2016; 346: 66–77. 12. Sopwith DG. The distribution of load in screw threads. Proc Inst Mech Eng 1948; 159: 373–383. 13. Yamamoto A. The theory and computation of threads connection[M]. Tokyo: Yokendo, 1980 (in Japanese).
