Open Access
 Issue Mechanics & Industry Volume 21, Number 6, 2020 616 8 https://doi.org/10.1051/meca/2020090 04 December 2020

## 1 Introduction

Different solutions can be adopted for this purpose. If the nut material is more flexible than the bolt, the force distribution in the threads becomes more uniform, and experimental tests prove that the fatigue resistance of the bolt significantly increases (between 35 and 60%, depending on the materials) [8,12,13]. Non-conventional nut designs such as tapered threads or tension nuts are also available, which have shown a better response to cyclic loading than standard nuts [7,8,13,14]. Finally, differential thread spacing (different thread pitch values in the bolt and the nut) may also be used to obtain a more uniform load distribution [6,15,16].

In this sense, this work develops a semianalytical model to estimate the thread load distribution under different thread pitch values for the bolt and the nut. The values for the parameters of the model are obtained from a Finite Element model of the bolted joint. Thus, the optimum differential thread spacing that provides the desired thread load distribution can be worked out. The results of the model are validated via FEA and compared with those predicted by the classical formulation of Sopwith .

## 2 Materials and methods

A bolted joint consisting of three threads (see Fig. 1) will be analyzed to explain the model. The pitch in the bolt and the nut are p and p + Δp, respectively. Due to the pitch difference, initially only the thread number 1 is in contact, carrying the whole preload. However, as the preload increases during the bolt tightening operation, the other threads come into contact, and therefore the preload is distributed among them. Finally, when every thread is in contact, the remaining load (until the target preload PT is reached) is distributed between all of them.

The proposed model assumes the behaviour schematized in Figure 2. The horizontal axis of the graph shows the preload (which increases as the bolt is being tightened) and the vertical axis shows how this preload is distributed among the threads. corresponds to the preload value when thread j comes into contact, is the preload value when thread j achieves its maximum percentage of load, and Mk,j is the percentage of preload that thread k is carrying at the preload value . As mentioned, due to the pitch difference, initially thread number 1 (lower thread) supports the whole preload, but when other threads get in contact, the preload is distributed among them. If the dashed line (thread 2) is observed, it can be divided into three different steps. First, when the preload reaches , the line gradually rises from zero to the maximum load percentage of thread 2, M2,2, which takes place when the preload is ; at this point, thread 2 is completely in contact. Then, the line decreases as the next threads start sharing the load. Thus, as the percentage of load in thread 3 increases (dotted line), threads 1 and 2 decrease their corresponding percentages. When the preload reaches , that is, when thread 3 comes into full contact, the load distribution gets stabilized (the percentages remains unaltered) until the tightening operation gets completed at PT because there are no more threads in this simple illustrative case.

Finite Element analyses were carried out with different configurations (bolt sizes and thread spacings), and it was observed that the response was well represented by the behaviour proposed in Figure 2. Besides, it was observed that, for any joint configuration, the and values were always proportional to the pitch difference Δp, while the Mk,j values were the same for any Δp.

Finally, two controversial aspects of the linear behaviour assumed in Figure 2 must be discussed. In first place, it could be argued that the maximum percentage of load in a thread should take place just when the next thread comes into contact, but Finite Element results have shown that this is only true for the thread number 1 in Figure 1 (for that reason, in Figure 2 but ). Second, and as a consequence of the first aspect, it can be observed that in the interval the total load percentage (i.e. the summation of the three percentages) exceeds one, which obviously is not possible. These aspects could probably be overcome by assuming a behaviour more complex than the one presented in Figure 2. Nevertheless, as it will be demonstrated in the Results and discussion section, in spite of these aspects, the model that has been developed from Figure 2 provides very accurate results. Thus, there was no need for refining the model given its simplicity and efficiency as an engineering tool.

From Figure 2, it is deducted that the total load that any thread is carrying at a given value of the preload, is equal to the area below its corresponding line in the interval between null preload and the preload under study. As illustrated in Figure 3, the load increment that a thread will experiment when the preload increases from P1 to P2 can be easily calculated as: (1)

According to equation (1) and Figure 2, the total load Fj in thread j under preload PT can be calculated as the sum of three terms: (2)

The first term Fj1 is the load in thread j when that thread is completely in contact, and therefore achieves its maximum load percentage: (3)

The second term Fj2 is the additional load that thread j carries as the rest of the threads achieve complete contact: (4)where N is the total number of threads in the bolted joint (N = 3 in the example of Fig. 2). Finally, Fj3 is the additional load that thread j carries from the point where all the threads are in complete contact up to the target preload value PT: (5)

In order to obtain the values of , and Mk,j for a particular bolted joint with a given pitch difference value Δp, a Finite Element analysis of the bolted joint with an arbitrary pitch difference Δp FE must be performed. In that model, the bolt preload is gradually increased. As a result of that analysis, the preload values at which each thread j enters into contact, the values at which it reaches its maximum load percentage and the values Mk,j FE with the percentage of preload that each thread k is carrying at the preload value are obtained for that Δp FE . As previously mentioned, the values of and are proportional to the pitch difference Δp, whereas Mk,j remains the same for any Δp. Thus, , and Mk,j for any other pitch difference Δp are calculated as: (6)

Thus, assuming that every thread is in contact at the end of the assembly (PT > PMN), equations (3)– (5) can be rewritten as: (7)

Thus, the load in each thread j for any given pitch difference value Δp under a preload PT can be calculated by replacing equations (7) in equation (2). This equation can further be used to assess the optimum pitch difference Δpoptimum that will provide a uniform load distribution among all the threads. Achieving a completely uniform load distribution is not possible because the pitch difference would have to be different for each thread. Thus, the best strategy is to establish that the first and last threads must have the same load, as proposed by Sopwith : (8)

Developing equation (8), and substituting the terms of equation (7):  where, and are the values of F12 and FN1 if Δp = Δp FE . From this equation, the differential thread spacing Δpoptimum that gives a nearly uniform load distribution along threads can be obtained: (9)

The optimum pitch difference is expressed as a function of the bolt preload. If the tightening process is not preload controlled (by using bolt tensioners, ultrasonic devices, strain gauges or similar) but torque controlled (with a torque wrench), PT in equation (9) could be substituted by any classical formula that expresses it in terms of the applied tightening torque, such as the classical one by Motosh . Fig. 1Illustrative example of a bolted joint with three threads in contact. Fig. 2Preload distribution among threads during preload operation: =preload when thread j gets in contact; = preload when thread j achieves maximum load percentage; Mk,j = preload percentage in thread k at preload . Fig. 3Calculation of the load increment in a thread under a preload increment.

## 3 Results and discussion

A bolted joint with M10 × 1 bolt under a preload of 12 kN was studied. Figure 4 shows the Finite Element model with its components and fixed boundary conditions. The location of the fixed surfaces, that is, how far they are from the bolt, will affect the stiffness of the whole assembly, and consequently it will play an important role on the load distribution among the threads under bolt preload; in this sense, the size of the material surrounding the bolt must be carefully defined in order to reproduce the stiffness and therefore the behaviour of the bolted joint in a realistic way. The components of the bolted joint were made of steel with E = 210GPa and ν = 0.3, and the friction coefficient between bolt and nut threads was µ = 0.2. The preload was applied by means of a pretension section. High order elements were used, resulting in 1.8M DoF.

As mentioned, an arbitrary pitch difference is analysed as a starting point (Δp FE ): in this case, Δp FE  = 1.5μm was used, that is, the nut pitch was defined as 1.0015 mm. Thus, preload is increasingly applied to the bolt of the FE model, obtaining the preload values at which each thread j enters into contact, the values at which it reaches its maximum load percentage and the values Mk,j with the percentage of preload that each thread k is carrying at the preload value . Table 1 lists the obtained values of these parameters , and Mk,j, to be introduced in the expressions of the previous section, and Figure 5 illustrate the load distribution obtained in the threads. As it can be appreciated, the behaviour is accurately estimated by the proposed model, the average error being 1.5%. In the table and figure, the threads are numbered as in Figure 1, from the bottom to the top of the bolt shank.

As explained, equation (9) can then be used to calculate the optimum pitch difference, which for the studied case (bolted joint and target load) provides a value of 0.49 µm. In contrast, using the formula proposed by Sopwith , a value of 0.63 µm is obtained. (10)

Where A1 and A2 are the cross sectional areas of bolt and nut respectively. As it can be observed, a difference of 30% can be found between both values. Figure 6 shows the thread load distribution obtained by means of the proposed model and the Finite Element model for these two values of the optimum pitch difference, together with the results for 0 and 1.5 µm. The results show that with the value of 0.49 µm the load in the first and last thread are the same (as imposed in equation (8), from which equation (9) was derived), according both to the proposed model and the Finite Element analysis. Using 0.63 µm this condition is not so accurately fulfilled. Besides, and what is more important, using the optimum pitch difference of equation (9) the load distribution is more uniform than using the value provided by Sopwith in (10), being the standard deviations 290 N and 330 N, respectively. The model by Sopwith is fully analytical, based on some simplifying hypotheses. Despite these simplifications, the model is very useful for its efficiency, because it provides good results with a simple formula. The methodology developed in the present work is semi-analytical, meaning that the values of the parameters that are used to reproduce the behavior of the bolted joint are tuned with the results of a Finite element model of the joint with an arbitrary thread spacing value. In this sense, it gives more accurate results than Sopwith́s model, but at the expense of a higher cost.

With this uniform load distribution, the thread stripping behaviour is much improved. In terms of fatigue response, as mentioned in the Introduction section, the stress concentration in the first engaged thread (i.e. the upper one, not to be confused with thread number 1 in Fig. 1) will presumably be reduced by a more uniform load distribution among threads. As it can be seen in Figure 7, the peak von Mises stress value in the root of the first engaged thread is 1145 MPa without differential thread spacing (0 µm) and 957 MPa with the optimum value of 0.489 µm. This means that the stress concentration factor kt is reduced almost a 20%, and consequently the fatigue response will be improved. Finally, it should be mentioned that these tolerance magnitudes are not feasible in general applications and steel bolts, but they may be possible in special applications with high precision bolts and light metal alloys. Fig. 4Finite Element model of the bolted joint.
Table 1 , and Mk,j values for the case under study (threads numbered as in Fig. 1). Fig. 5Load distribution during preload process: FE vs proposed model (threads numbered as in Fig. 1). Fig. 6Load distribution along threads at the end of the tightening operation (threads numbered as in Fig. 1). Fig. 7Von Mises stress in the first engaged thread (upper thread) under a preload of 12 kN with a differential thread spacing of: (a) 0 µm, (b) 0.489 µm.

## 4 Conclusions

A new model to estimate the thread load distribution in preloaded bolted joints with differential thread spacing between bolt and nut has been presented. Based on a preliminary Finite Element analysis and some simple analytical expressions, the load distribution for any differential thread spacing can be worked out. Furthermore, a formula to calculate the optimum pitch difference that leads to a uniform load distribution is derived, which improves the mechanical behaviour of the joint. The model results are almost coincident with Finite Element results, thus improving the accuracy of classical formulation.

## Acknowledgments

The authors want to acknowledge the financial support of the Spanish Ministry of Economy and Competitiveness [grant number DPI2017-85487-R (AEI/FEDER, UE)] and the Basque Government [grant number IT947-16].

## References

1. I. Coria, M. Abasolo, J. Aguirrebeitia, I. Heras, Study of bolt load scatter due to tightening sequence, Int. J. Press Vessels Pip 182 , 104054 (2020) [CrossRef] [Google Scholar]
2. A. Bhattacharya, A. Sen, S. Das, An investigation on the anti-loosening characteristics of threaded fasteners under vibratory conditions, Mech. Mach. Theory 45 , 1215–1225 (2010) [CrossRef] [Google Scholar]
3. I. Coria, M. Abasolo, I. Olaskoaga, A. Etxezarreta, J. Aguirrebeitia, A new methodology for the optimization of bolt tightening sequences for ring type joints. Ocean Eng. 129 , 441–450 (2017) [CrossRef] [Google Scholar]
4. I. Coria, I. Martín, A.H. Bouzid, I. Heras, M. Abasolo, Efficient assembly of bolted joints under external loads using numerical FEM. Int. J. Mech. Sci. 142–143 , 575–582 (2018) [CrossRef] [Google Scholar]
5. V. Strizhak, I. Penkov, International Symposium on History of Machines and Mechanisms Proceedings HMM 2000 Distribution of axial load on bolt and nut threads , Springer-Science+Business Media, 2000 [Google Scholar]
6. C. Brutti, Load and stress distribution in screw threads with modified washers. J. Multidiscip. Eng. Sci. Technol. 4 , 6523–6532 (2017) [Google Scholar]
7. J.H. Bickford, An introduction to the design and behavior of bolted joints, M. Dekker, 1981 [Google Scholar]
8. W.D. Pilkey, Peterson's stress concentration factors , John Wiley & Sons, Inc., New York, 1997 [CrossRef] [Google Scholar]
9. N. Motosh, Load distribution on threads of titanium tension nuts and steel bolts, J. Eng. Ind. 97 , 162–166 (1975) [CrossRef] [Google Scholar]
10. D.L. Miller, K.M. Marshek, M.R. Naji, Determination of load distribution in a threaded connection, Mech. Mach. Theory 18 , 421–430 (1983) [CrossRef] [Google Scholar]
11. N.L. Pedersen, Overall bolt stress optimization, J. Strain Anal. Eng. Des. 48 , 155–165 (2013) [CrossRef] [Google Scholar]
12. F. Kaufmann, W. Janiche, Beitrag zur Dauerhaltbarkeit von Schraubenverbindungen. Tech Mitt Krupp, Forschungsber 3 , 147 (1940) [Google Scholar]
13. H. Wiegand, Uber die Dauerfestigkeit von Schraubenwerkstoffen und Schraubenverbindungen. Thesis, Technische Hochshule Darmstadt 1933 [Google Scholar]
14. N.A. Noda, Y. Xiao, M. Kuhara, The reduction of stress concentration by tapering threads, J Solid Mech. Mater. Eng. 5 , 397–408 (2011) [CrossRef] [Google Scholar]
15. D.G. Sopwith, The distribution of load in screw threads, Proc IME 159 , 373–383 (1948) [CrossRef] [Google Scholar]
16. N.A. Noda, X. Chen, Y. Sano, M.A. Wahab, H. Maruyama, R. Fujisawa, Y. Takase, Effect of pitch difference between the bolt–nut connections upon the anti-loosening performance and fatigue life, Mater. Des. 96 , 476–489 (2016) [CrossRef] [Google Scholar]
17. N.E. Joukovsky, Collected Papers, Vol. VIII, Theory of elasticity, railways, automobiles, ONTI, Moscow, 1937, pp. 48–56 (in Russian) [Google Scholar]
18. E. Jaquet, Dber eine newartige Schraubenverbindungen. Ingenieur-Archiv, Bd.2, 1931 [Google Scholar]
19. L. Maduschka, Beansprughung von Schraubenverbindungen und Zweckmassige Gestaltung der Gewindetrager, Forsch. Geb. Ingenierwessens, 1936, 8.7 [Google Scholar]
20. I.A. Birger, Calculation of threaded joints, Oborongiz, Moscow, 1959 (in Russian) [Google Scholar]
21. I.A. Birger, G.B. Iosilevich, Threaded and flanged joints, Mechanical Engineering, Moscow, 1990 (in Russian) [Google Scholar]

Cite this article as: I. Coria, M. Abasolo, A. Gutiérrez, J. Aguirrebeitia, Achieving uniform thread load distribution in bolted joints using different pitch values, Mechanics & Industry 21, 616 (2020)

## All Tables

Table 1 , and Mk,j values for the case under study (threads numbered as in Fig. 1).

## All Figures Fig. 1Illustrative example of a bolted joint with three threads in contact. In the text Fig. 2Preload distribution among threads during preload operation: =preload when thread j gets in contact; = preload when thread j achieves maximum load percentage; Mk,j = preload percentage in thread k at preload . In the text Fig. 3Calculation of the load increment in a thread under a preload increment. In the text Fig. 4Finite Element model of the bolted joint. In the text Fig. 5Load distribution during preload process: FE vs proposed model (threads numbered as in Fig. 1). In the text Fig. 6Load distribution along threads at the end of the tightening operation (threads numbered as in Fig. 1). In the text Fig. 7Von Mises stress in the first engaged thread (upper thread) under a preload of 12 kN with a differential thread spacing of: (a) 0 µm, (b) 0.489 µm. In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.