A new computational method for threaded connection stiffness

Transcription

1 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: / 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 Motash 2 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 Abbound 9 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 method 12 and the Yamamoto method 13 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 , P.R. China. dongmei_zhang@163.com Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 3.0 License ( 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 ( open-access-at-sage).

2 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 d 1 due to bending moment, the deflection d 2 due to shear loading, the deflection d 3 due to incline at the root of the thread, the deflection d 4 due to shear at the root of the thread, and the deflection d 5 due to radial contraction and expansion on the bolt and the nut. For the ISO thread under unit force at x = l 1 as shown in Figure 2, the axial elastic deflections of unit length thread of the bolt are Figure 2. Geometrical parameters of ISO thread. d 1b = 0:034 1 y2 b E b d 2b = 1: y b E b d 3b = 0:229 1 y2 b E b d 4b = 1:18 1 y2 b E b d 5b = 0:056(1 y b ) 1 D E b p ð1þ The axial elastic deflections of unit length thread of nut are d 1n = 0:073 1 y2 n E n d 2n = 1:15 1 y n E n d 3n = 0:294 1 y2 n E n d 4n = 1:14 1 y2 n E n d 5n = 0:056 D2 0 + D2 1 D y D D2 n E n p ð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. D 0 is the outer diameter of nut.

3 Zhang et al. 3 For the engaged bolt and nut, the axial deflection is D y = D b + D n = 1 F k by y + 1 F k ny y = F k by k ny y ð13þ The stiffness of unit axial length is Figure 3. Diagram of lead angle. k y = F= y 1 = D y (d b + d n ) sin b ð14þ Therefore, for the thread on bolt and nut, the elastic deflection caused by unit force per unit width along axial direction is d b = d 1b + d 2b + d 3b + d 4b + d 5b d n = d 1n + d 2n + d 3n + d 4n + d 5n ð3þ ð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 F D b = d b ð5þ s F D n = d n ð6þ s where s is the thread length along the direction of the helix, and s = y ð7þ sin b where y is the axial length of thread and b is the lead angle as shown in Figure 3. Therefore F s = F y y s = F y sin b ð8þ Substituting equation (8) into equations (5) and (6) leads to D b = d b sin b F y = 1 F k by y D n = d n sin b F y = 1 F k ny y ð9þ ð10þ where k by and k ny represent the stiffness of unit axial length for thread on bolt and nut, respectively, and they can be expressed as k by = F= y D b = 1 d b sin b k ny = F= y D n = 1 d n sin b ð11þ ð12þ 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 e b generated by the bolt body under compressive force F(y) along y-direction is e b = F(y) ð15þ S b E b where S b is the cross-sectional area of the bolt body and E b is Young s modulus of the bolt body. Similarly, strain e n generated by the nut body under compressive force F(y) is e n = F(y) ð16þ S n E n where S n is the cross-sectional area of the nut body and E n is Young s modulus of the nut body. According to equations (9) and (10), the displacement gradients are respectively D b y = 1 2 F k by y 2 D n y = 1 2 F k ny y 2 ð17þ ð18þ In order to satisfy the deformation compatibility condition, e b, e n, D b, and D n need to meet the following relationship

4 4 Advances in Mechanical Engineering e b + e n = D b y + D n y ð19þ where Ð l 0 u(y)k bydy is the total thread stiffness of nut K b. According to equation (28), K b can be expressed as Substituting equations (15) (18) into equation (19) leads to l 2 F = 2 F ð20þ y 2 qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi (1=S where l = b E b )+(1=S n E n ) (1=k by )+(1=k ny ). The solution of equation (20) is F(y)=C 1 sinh ly + C 2 cosh ly ð21þ By the boundary conditions F(y = 0)=F B and F(y = l)=0, we can obtain C 2 = F B and C 1 = sinh ll F B. Substituting C 1 and C 2 into equation (21), the axial load of threaded connections can be expressed as F(y)=F B sinh ly + cosh ly sinh ll ð22þ where F B is the load on the first thread of the bolt as shown in Figure 1. Total stiffness of threaded connections According to the above deflections of the thread, the law of axial load distribution can be expressed as u(y)= F(y) F B = sinh ly +coshly sinh ll ð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 D b (y)=d b (0)u(y) D n (y)=d n (0)u(y) ð24þ ð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 D b (0)u(y)k by = F y ð26þ D n (0)u(y)k ny = F ð27þ y For the bolt, because the load is F B at y = 0, equation (26) can be changed into ð l D b (0) 0 u(y)k by dy = F B ð28þ K b = ð l 0 u(y)k by dy = F B D b (0) ð29þ Substituting equation (23) into equation (29), the thread stiffness of bolt can be expressed as K b = k by l = E b S b l + sinh ll + sinh ll sinh ll 1 sinh ll ð30þ Similarly, the thread stiffness of nut can be expressed as K n = k ny l = E n S n l + sinh ll + sinh ll sinh ll 1 sinh ll ð31þ Because the meshing bolt and nut is in series manner, the total stiffness of engaged thread is K = 1 1 K b + 1 K n = sinh ll l( 1) Finite element model E b S b E n S n ð32þ 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 d average 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

5 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 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 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,

7 Zhang et al. 7 Table 1. Parameters of samples. No. Size code of threads Material Parameters for material Ultimate strength for thread connections (MPa) Young s modulus (GPa) Poisson ratio Ultimate strength (MPa) 1# M Steel # M Copper # M Brass # M # M # M # M 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 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 8 Advances in Mechanical Engineering Table 2. Stiffness of thread connections (10 9 N=m). No. Size code of threads Exp. FEA Theory This study Yamamoto s method 1# M # M # M # M # M # M # M 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: 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 1. 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. 2. 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. 3. 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

9 Zhang et al. 9 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: Motash N. Determination of joint stiffness in bolted connections. J Eng Ind: T ASME 1976; 98: Wileman J, Choudhury M and Green I. Computation of member stiffness in bolted connections. J Mech Design 1991; 113: Shigley JE, Mischke CR and Budynas RG. Mechanical engineering design. New York: McGraw-Hill, Lehnhoff TF and Wistehuff WE. Nonlinear effects on the stresses and deformations of bolted joints. J Mech Design 1996; 118: Naser SA. Computation of member stiffness in bolted connections using the finite element analysis. Mech Base Des Struct Mach 2005; 33: Musto JC and Konkle NR. Computation of member stiffness in the design of bolted joints. J Mech Design 2006; 128: Nawras H. Mathematical representation of bolted-joint stiffness: a new suggested model. J Mech Sci Technol 2011; 25: Nassar SA and Abbound A. An improved stiffness model for bolted joints. J Mech Design 2009; 131: 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: Liu JH, Ouyang HJ, Peng JF, et al. Experimental and numerical studies of bolted joints subjected to axial excitation. Wear 2016; 346: Sopwith DG. The distribution of load in screw threads. Proc Inst Mech Eng 1948; 159: Yamamoto A. The theory and computation of threads connection[m]. Tokyo: Yokendo, 1980 (in Japanese).