# Stiffness Calculation Model of Thread Connection Considering Friction Factors

Authors: Shi-kun Lu; Deng-xin Hua; Yan Li; Fang-yuan Cui; Peng-yang Li
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

In order to design a reasonable thread connection structure, it is necessary to understand the axial force distribution of threaded connections. For the application of bolted connection in mechanical design, it is necessary to estimate the stiffness of threaded connections. A calculation model for the distribution of axial force and stiffness considering the friction factor of the threaded connection is established in this paper. The method regards the thread as a tapered cantilever beam. Under the action of the thread axial force, in the consideration of friction, the two cantilever beams interact and the beam will be deformed, these deformations include bending deformation, shear deformation, inclination deformation of cantilever beam root, shear deformation of cantilever beam root, radial expansion deformation and radial shrinkage deformation, etc.; calculate each deformation of the thread, respectively, and sum them, that is, the total deformation of the thread. In this paper, on the one hand, the threaded connection stiffness was measured by experiments; on the other hand, the finite element models were established to calculate the thread stiffness; the calculation results of the method of this paper, the test results, and the finite element analysis (FEA) results were compared, respectively; the results were found to be in a reasonable range; therefore, the validity of the calculation of the method of this paper is verified.


---

## Body

## 1. Introduction

The bolts connect the equipment parts into a whole, which is used to transmit force, moment, torque, or movement. The bolt connection is widely used in various engineering fields, such as aviation machine tools, precision instruments, etc. The precision of the threaded connection affects the quality of the equipment. Especially for high-end CNC machine tools, the precision of the thread connection is very high. Therefore, it is important to study the stiffness of the threaded connection to improve the precision of the device. Many researchers have conducted research in this area.Dongmei Zhang et al. [1], propose a method which can compute the engaged screw stiffness, and the validity of the method was verified by FEA and experiments. Maruyama et al. [2] used the point matching method and the FEM, based on the experimental results of Boenick and the assumptions made by Fernlund [3] in calculating the pressure distribution between joint plates. Motash [4] assumed that the pressure distribution on any plane perpendicular to the bolt axis had zero gradient atr = dh/2 andr = ro, where it also vanishes. They mainly use numerical methods to calculate the influence of different parameters on the stiffness of bolted connections. Kenny B, Patterson E A. [5] introduced a method for measuring thread strains and stresses. Kenny B [6] et al. reviewed the distribution of loads and stresses in fastening threads. Miller D L et al. [7] established the spring model of thread force analysis and, combined with the mathematical theory, analyzed the stress of the thread and compared with the FEA results and experimental results to verify the correctness of the spring mathematical model. Wang W and Marshek K M. et al. [8] proposed an improved spring model to analyze the thread load distribution, compared the load distributions of elastic threads and yielding thread joints, and discussed the effect of the yield line on the load distribution. Wileman et al. [9] performed a two-dimensional (2D) FEA for members stiffness of joint connection. De Agostinis M et al. [10] studied the effect of lubrication on thread friction characteristics or torque. Dario Croccolo et al. [11–13] studied the effect of Engagement Ratio (ER, namely, the thread length over the thread diameter) on the tightening and untightening torque and friction coefficient of threaded joints using medium strength threaded locking devices. Zou Q. et al. [14, 15] studied the use of contact mechanics to determine the effective radius of the bolted joint and also studied the effect of lubrication on friction and torque-tension relationship in threaded fasteners. Nassar S. A. et al. [16, 17] studied the thread friction and thread friction torque in thread connection. Nassar S. A. et al. [18, 19] also investigated the effects of tightening speed and coating on the torque-tension relationship and wear pattern in threaded fastener applications in order to improve the reliability of the clamping load estimation in bolted joints. Kopfer et al. [20] believe that suitable formulations should consider contact pressure and sliding speed; based on this, the contribution shows experimental examples for main uncertainties of frictional behavior during tightening with different material combinations (results from assembly test stand). Kenny B et al. [21] reviewed the distribution of load and stress in the threads of fasteners. Shigley et al. [22] presented an analytical solution for member stiffness, based on the work of Lehnhoff and Wistehuff [23]. Nasser [24], Musto and Konkle [25], Nawras [26], and Nassar and Abbound [27] also proposed mathematical model for the bolted-joint stiffness. Qin et al. [28] established an analytical model of bolted disk-drum joints and introduced its application to dynamic analysis of joint rotor. Liu et al. [29] conducted experimental and numerical studies on axially excited bolt connections.There are also several authors that, starting from the nature of thread stiffness, from the perspective of thread deformation, established a mathematical model of the calculation of the distribution of thread axial force. The Sopwith method [30] and the Yamamoto method [31] received extensive recognition. The Sopwith method gave a method for calculating the axial force distribution of threaded connections. Yamamoto method can not only calculate the axial force distribution of threads but also calculate the stiffness of threaded connections. The assumption for Yamamoto method is that the load per unit width along the helix direction is uniformly distributed. In fact, for the three-dimensional (3D) helix thread, the load distribution is not uniform. Therefore, based on the Yamamoto method, Dongmei Zhang et al. [1] propose a method which can compute the engaged screw stiffness by considering the load distribution, and the validity of the method was verified by FEA and experiments. The method of Zhang Dongmei does not consider the influence of the friction coefficient of the thread contact surface. In fact, the friction coefficient of the contact surface of the thread connection has an influence on the distribution of the axial force of the thread and the stiffness of the thread. Therefore, we propose a new method which can compute the engaged screw stiffness more accurately by considering the effects of friction and the load distribution. The accuracy of the method was verified by the FEA and bolt tensile test. The flow chart of the article is shown in Figure 2.


## 2. Mathematical Model

### 2.1. Axial Load Distribution

According to Yamamoto [31], the thread is regarded as a cantilever beam, and the thread is deformed under axial force and preload. These deformations include the following (shown in Figure 3): thread bending deformation, thread shear deformation, thread root inclination deformation, thread root shear deformation, radial direction extended deformation, or radial shrinkage deformation.For theISO thread, the axial deformation of the thread at x at the axial unit width force wz is thread bending deformation δ1, thread shear deformation δ2, thread root inclination deformation δ3, thread root shear deformation δ4, and radial direction extended deformation δ5 (nut) or radial shrinkage deformation δ5 (screw), and calculate these deformations of the thread, respectively, and then sum them, that is, the total deformation.


#### 2.1.1. Bending Deformation

In the threaded connection, under the action of the load, the contact surface friction coefficient isμ, when the sliding force along the inclined plane is greater than the friction force along the inclined plane, the relative sliding occurs between the two inclined planes, and the axial unit width force (shown in Figure 3) is wz; if the influence of the lead angle is ignored, the force per unit width perpendicular to the thread surface can be expressed as(1)w=wzμsin⁡α+cos⁡αThe force per unit width perpendicular to the thread surface can be decomposed into thex-direction component force and they-direction component force, respectively(2)wcos⁡α=wzcos⁡αμsin⁡α+cos⁡αand(3)wsin⁡α=wzsin⁡αμsin⁡α+cos⁡αThe friction generated along the slope iswμ; i.e.,(4)wμ=wzμμsin⁡α+cos⁡αThe forcewμ is also decomposed intox-direction force andy-direction force, which are μsin⁡α and μcos⁡α, respectively.(5)wμsin⁡α=wzμsin⁡αμsin⁡α+cos⁡α(6)wμcos⁡α=wzμcos⁡αμsin⁡α+cos⁡αIn the unit width, the thread is regarded as a rectangular variable-section cantilever beam. Under the action of the above-mentioned force, the thread undergoes bending deformation, and the virtual work done by the bending momentE- on the beam dy section is(7)E-dθ=E-MwEbIydyAccording to the principle of virtual work, the deflectionδ1 (see Figure 3(a)) of the beam subjected to the load is(8)δ1=∫0cE-MwEbIydywhereE- is the bending moment of the unit load beam. Mw is the bending moment of the beam under the actual load.I(y) is the area moment of inertia of the beam at y. Eb is Young’s Modulus of the material.c is the length of the beam. Here, the forces are assumed as acting on the mean diameter of the thread.As shown in Figure5, the height h(y) of the beam section per unit width and the area moment of inertia I(y) of the section can be expressed by using the function interpolation.(9)hy=1+β1-1c-ych,0≤y≤cIy=112bh31+β1-1c-yc3,0≤y≤cβ1=ahwhereh is the beam end section height;b is the beam section width; β1 is the beam root section height and the beam end section height ratio; see Figure 5.From Figure5, the bending moment of the beam is related to they-axis component of w andwμ, and these components cause the beam to bend; therefore, the analytical solution shows that the bending moment of the unit width beam subjected to the friction force and the vertical load of the thread surface is(10)Mw=wzμsin⁡α+cos⁡α·cos⁡α·c+μsin⁡α·c-sin⁡αa2-ctan⁡α+μcos⁡αa2-ctan⁡αSubstituting (10) and (9) to (8) and integrating to obtain the analytical expression of the deflection δ1 (shown in Figure 3(a)) of the cantilever beam with variable cross-section under load one has(11)δ1=12Mwc2Ebbh3·121β12+1-1β1·1β1-12


#### 2.1.2. Shear Deformation

Assume that the distribution of shear stress on any section is distributed according to the parabola [31] and the deformation δ2 (see Figure 3(b)) caused by the shear force within the width of unit 1 is(12)δ2=logeab·61+vcos⁡α+μsin⁡αcot⁡α5Eb·wzμsin⁡α+cos⁡α


#### 2.1.3. Inclination Deformation of the Thread Root

Under the action of the load, the thread surface is subjected to a bending moment, and the root of the thread is tilted, as shown in Figure3(c). Due to the inclination of the thread, axial displacement occurs at the point of action of the thread surface force, and the axial displacement can be expressed as [31](13)δ3=wzμsin⁡α+cos⁡α·12c1-v2πEba2·cos⁡α·c+μsin⁡α·c-sin⁡αa2-ctan⁡α+μcos⁡αa2-ctan⁡α


#### 2.1.4. Deformation due to Radial Expansion and Radial Shrinkage

According to the static analysis, the thread is subjected to radial forcewsin⁡α-wμcos⁡α (shown in Figure 4), and it is known from the literature [31] that the internal and external thread radial deformation (shown in Figure 3(d)) are(14)δ4b=1-νbdp2EbPtan⁡α·wzsin⁡α-μcos⁡αμsin⁡α+cos⁡αand(15)δ4n=D02+dp2D02-dp2+νndp2EnP·tan⁡α·wzsin⁡α-μcos⁡αμsin⁡α+cos⁡α


#### 2.1.5. Shear Deformation of the Root

Assuming that the shear stress of the root section is evenly distributed, the displacement of theO point in the x direction caused by the shear deformation (shown in Figure 3(e)) is the same as the displacement of the thread in the x direction; this displacement can be expressed as [31](16)δ5=wzμsin⁡α+cos⁡α·21-v2·cos⁡α+μsin⁡απEb·PalogeP+a/2P-a/2+12loge4P2a2-1For ISO internal threads, the relationship between a, b, c, and pitchP is(17)a=0.833Pb=0.5Pc=0.289PSubstituting (17) into (9), (10), (11), (12), (13), (14), and (16) one gets the relation(18)δ1b=3.468MwbEbP2βb3·121βb2+1-1βb·1βb-12,βb=1.6684(19)Mwb=wzμsin⁡α+cos⁡α0.289P·cos⁡α+μ·sin⁡α-0.4165P-0.289P·tan⁡αsin⁡α-μ·cos⁡α(20)δ2b=0.51·61+vcos⁡α+μsin⁡αcot⁡α5Eb·wzμsin⁡α+cos⁡α(21)δ3b=wzμsin⁡α+cos⁡α·12c1-v2πEba2·0.289P·cos⁡α+μ·sin⁡α-0.4165P-0.289P·tan⁡αsin⁡α-μ·cos⁡α(22)δ4b=1-νbdp2EbPsin⁡α-μcos⁡αtan⁡α·wzμsin⁡α+cos⁡α(23)δ5b=1.8449wzμsin⁡α+cos⁡α·21-v2·cos⁡α+μsin⁡απEbFor ISO internal threads, the relationship betweena,b,c, and pitch P is(24)a=0.875Pb=0.5Pc=0.325PSubstituting (24) into (9), (10), (11), (12), (13), (15), and (16) type one gets the relation(25)δ1n=3.784MwnEnP2βn2·121βn2+1-1βn·1βn-12,βn=1.751(26)Mwn=wzμsin⁡α+cos⁡α0.325P·cos⁡α+μ·sin⁡α-0.4375P-0.325P·tan⁡αsin⁡α-μ·cos⁡α(27)δ2n=0.55962·61+vcos⁡α+μsin⁡αcot⁡α5Eb·wzμsin⁡α+cos⁡α(28)δ3n=wzμsin⁡α+cos⁡α·12c1-v2πEba2·0.325P·cos⁡α+μ·sin⁡α-0.4375P-0.325P·tan⁡αsin⁡α-μ·cos⁡α(29)δ4n=D02+dp2D02-dp2+νndp2EnP·sin⁡α-μcos⁡αtan⁡α·wzμsin⁡α+cos⁡α(30)δ5n=1.7928wzμsin⁡α+cos⁡α·21-v2·cos⁡α+μsin⁡απEbBy adding these deformations separately, the total deformation (shown in Figure4) of screw thread and nut thread can be obtained under the action of force wz.(31)δb=δ1b+δ2b+δ3b+δ4b+δ5b(32)δn=δ1n+δ2n+δ3n+δ4n+δ5nThe unit force per unit width of the axial direction can be expressed as(33)fΔ=wzwz=1Under the action of unit force of axial unit width, the total deformation of external thread and internal thread is(34)δb1=δ1b+δ2b+δ3b+δ4b+δ5bwzand(35)δn1=δ1n+δ2n+δ3n+δ4n+δ5nwzFor threaded connections, at thex-axis of the loadF, the axial deformation of screws and nuts can be expressed as(36)zbr=δb1·∂F∂r(37)znr=δn1·∂F∂rHerer is the length along the helical direction, and the relation between the axial height x and the length along the helix direction can be represented by the following formula according to the geometric relation shown in Figure 6.(38)r=xsin⁡βHere,β is the lead angle of the thread shown in Figure 6, and then(39)znx=δn1·∂F∂r=δn1·∂F∂x·∂x∂r=δn1·sin⁡β·∂F∂x(40)zbx=δb1·∂F∂r=δb1·∂F∂x·∂x∂r=δb1·sin⁡β·∂F∂xAssume(41)∂F∂x·zbx=1δb1·sin⁡β=kubxx(42)∂F∂x·znx=1δn1·sin⁡β=kunxxHere,kbx(x) and knx(x) represent the stiffness of the unit axial length of the nut and the screw, respectively, for the unit force.The axial total deformation of the threaded connection atx is denoted as(43)zxx=zbx+znxThe stiffness of the unit axial length of the threaded connection is expressed as(44)kuxx=∂F∂x·zxx=∂F∂x÷δb1·sin⁡β·∂F∂x+δn1·sin⁡β·∂F∂x=1δn1+δb1·sin⁡βAs shown in Figure1(a), the threaded connection structure includes a nut body and a screw body. The nut is fixed, the screw is subjected to pulling force, the total axial force is Fb, and the axial force at the threaded connection screw x isF(x). If the position of the bottom end face of the nut is the origin 0 and, at the x position, the axial force isF(x), the screw elongation amount εb and the nut compression εn can be obtained from the following:(45)εbx=FxAbxEb(46)εnx=FxAnxEnFigure 1
Thread force.
(a) (b)Figure 2
Schematic diagram of the article.Figure 3
Thread deformation caused by various reasons.
(a) (b) (c) (d) (e)Figure 4
The comprehensive deformation of the thread.
(a) (b)Figure 5
The force on the thread.Figure 6
The lead angle of the thread.whereAb(x) and An(x) are the vertical cross-sectional areas of screws and nuts at the x position. Eb and En are, respectively, Young’s modulus of the screw body and Young’s modulus of the nut body. Find the displacement gradient for the expression, which is, respectively, expressed as(47)∂zbx∂x=1kubxx·∂2F∂x2(48)∂znx∂x=1kunxx·∂2F∂x2Here,kubxx=1/δb1·sin⁡β, and kunxx=1/δn1·sin⁡β.As shown in Figure1(a), the screw is subjected to the tensile force Fb, with the bottom of the nut as the coordinate origin, and the force at the x position is Fx, and then the elongation of the screw at x is ∫xLεbxdx=wb, and the compressed shortening amount of the nut at x is ∫xLεnxdx=wn. The relationship between wb, wn, zb, and zn is ∫xLεbxdx+∫xLεnxdx=zbx+znxx=L-zbx+znxx=x (see Figures 7 and 1(a),), and the partial derivative of this relation can be obtained by the following formula:(49)εbx+εnx=∂zbx∂x+∂znx∂xFigure 7
Illustration of the elastic deformation of the screwed portion of the threaded connection.Substituting (45), (46), (47), and (48) into (49) and simplifying it(50)AnEn+AbEbAbEbAnEn·kubxkunxkubx+kunxFx=∂2Fx∂x2Let(51)n=AnEn+AbEbAbEbAnEn·kubxkunxkubx+kunxThen(52)n2Fx=∂2Fx∂x2From mathematical knowledge, the equation is a differential equation. The general solution of the equation can be expressed as(53)Fx=C1sinh⁡nx+C2cosh⁡nxAs can be seen from Figure1, the axial force at the first thread at the connection surface of the nut and the screw is Fb, and the axial force at the last thread at the lower end of the thread joint surface of the nut and screw is 0; that is, the boundary condition isF(x=0)=Fb andF(x=L)=0. Taking these boundary conditions into the equation will give C1=-Fb(cosh⁡(nL))/sinh⁡(nL) and C2=Fb, so we get the expression of the threaded connection axial load as(54)Fx=Fbcosh⁡nx-cosh⁡nLsinh⁡nLsinh⁡nxTherefore, the axial force distribution density of the thread connection along thex direction can be expressed as(55)fx=FxFb=cosh⁡nx-cosh⁡nLsinh⁡nLsinh⁡nx


### 2.2. Thread Connection Stiffness

The stiffness in the axial directionx of the bolted connection is equal to the axial force distribution of the threaded connection multiplied by the unit stiffness; i.e.,(56)kxx=fxkuxxThe overall stiffness of the bolt connection can be expressed as(57)Kc=∫0Lkxx=∫0LfxkuxxdxSubstituting (44) and (55) into (57), the stiffness of the bolt connection is expressed as(58)Kc=∫0Lkxx=∫0Lfxkuxxdx=1nδb1+δn1·sin⁡β·cosh⁡nL-1sinh⁡nL


## 2.1. Axial Load Distribution

According to Yamamoto [31], the thread is regarded as a cantilever beam, and the thread is deformed under axial force and preload. These deformations include the following (shown in Figure 3): thread bending deformation, thread shear deformation, thread root inclination deformation, thread root shear deformation, radial direction extended deformation, or radial shrinkage deformation.For theISO thread, the axial deformation of the thread at x at the axial unit width force wz is thread bending deformation δ1, thread shear deformation δ2, thread root inclination deformation δ3, thread root shear deformation δ4, and radial direction extended deformation δ5 (nut) or radial shrinkage deformation δ5 (screw), and calculate these deformations of the thread, respectively, and then sum them, that is, the total deformation.


### 2.1.1. Bending Deformation

In the threaded connection, under the action of the load, the contact surface friction coefficient isμ, when the sliding force along the inclined plane is greater than the friction force along the inclined plane, the relative sliding occurs between the two inclined planes, and the axial unit width force (shown in Figure 3) is wz; if the influence of the lead angle is ignored, the force per unit width perpendicular to the thread surface can be expressed as(1)w=wzμsin⁡α+cos⁡αThe force per unit width perpendicular to the thread surface can be decomposed into thex-direction component force and they-direction component force, respectively(2)wcos⁡α=wzcos⁡αμsin⁡α+cos⁡αand(3)wsin⁡α=wzsin⁡αμsin⁡α+cos⁡αThe friction generated along the slope iswμ; i.e.,(4)wμ=wzμμsin⁡α+cos⁡αThe forcewμ is also decomposed intox-direction force andy-direction force, which are μsin⁡α and μcos⁡α, respectively.(5)wμsin⁡α=wzμsin⁡αμsin⁡α+cos⁡α(6)wμcos⁡α=wzμcos⁡αμsin⁡α+cos⁡αIn the unit width, the thread is regarded as a rectangular variable-section cantilever beam. Under the action of the above-mentioned force, the thread undergoes bending deformation, and the virtual work done by the bending momentE- on the beam dy section is(7)E-dθ=E-MwEbIydyAccording to the principle of virtual work, the deflectionδ1 (see Figure 3(a)) of the beam subjected to the load is(8)δ1=∫0cE-MwEbIydywhereE- is the bending moment of the unit load beam. Mw is the bending moment of the beam under the actual load.I(y) is the area moment of inertia of the beam at y. Eb is Young’s Modulus of the material.c is the length of the beam. Here, the forces are assumed as acting on the mean diameter of the thread.As shown in Figure5, the height h(y) of the beam section per unit width and the area moment of inertia I(y) of the section can be expressed by using the function interpolation.(9)hy=1+β1-1c-ych,0≤y≤cIy=112bh31+β1-1c-yc3,0≤y≤cβ1=ahwhereh is the beam end section height;b is the beam section width; β1 is the beam root section height and the beam end section height ratio; see Figure 5.From Figure5, the bending moment of the beam is related to they-axis component of w andwμ, and these components cause the beam to bend; therefore, the analytical solution shows that the bending moment of the unit width beam subjected to the friction force and the vertical load of the thread surface is(10)Mw=wzμsin⁡α+cos⁡α·cos⁡α·c+μsin⁡α·c-sin⁡αa2-ctan⁡α+μcos⁡αa2-ctan⁡αSubstituting (10) and (9) to (8) and integrating to obtain the analytical expression of the deflection δ1 (shown in Figure 3(a)) of the cantilever beam with variable cross-section under load one has(11)δ1=12Mwc2Ebbh3·121β12+1-1β1·1β1-12


### 2.1.2. Shear Deformation

Assume that the distribution of shear stress on any section is distributed according to the parabola [31] and the deformation δ2 (see Figure 3(b)) caused by the shear force within the width of unit 1 is(12)δ2=logeab·61+vcos⁡α+μsin⁡αcot⁡α5Eb·wzμsin⁡α+cos⁡α


### 2.1.3. Inclination Deformation of the Thread Root

Under the action of the load, the thread surface is subjected to a bending moment, and the root of the thread is tilted, as shown in Figure3(c). Due to the inclination of the thread, axial displacement occurs at the point of action of the thread surface force, and the axial displacement can be expressed as [31](13)δ3=wzμsin⁡α+cos⁡α·12c1-v2πEba2·cos⁡α·c+μsin⁡α·c-sin⁡αa2-ctan⁡α+μcos⁡αa2-ctan⁡α


### 2.1.4. Deformation due to Radial Expansion and Radial Shrinkage

According to the static analysis, the thread is subjected to radial forcewsin⁡α-wμcos⁡α (shown in Figure 4), and it is known from the literature [31] that the internal and external thread radial deformation (shown in Figure 3(d)) are(14)δ4b=1-νbdp2EbPtan⁡α·wzsin⁡α-μcos⁡αμsin⁡α+cos⁡αand(15)δ4n=D02+dp2D02-dp2+νndp2EnP·tan⁡α·wzsin⁡α-μcos⁡αμsin⁡α+cos⁡α


### 2.1.5. Shear Deformation of the Root

Assuming that the shear stress of the root section is evenly distributed, the displacement of theO point in the x direction caused by the shear deformation (shown in Figure 3(e)) is the same as the displacement of the thread in the x direction; this displacement can be expressed as [31](16)δ5=wzμsin⁡α+cos⁡α·21-v2·cos⁡α+μsin⁡απEb·PalogeP+a/2P-a/2+12loge4P2a2-1For ISO internal threads, the relationship between a, b, c, and pitchP is(17)a=0.833Pb=0.5Pc=0.289PSubstituting (17) into (9), (10), (11), (12), (13), (14), and (16) one gets the relation(18)δ1b=3.468MwbEbP2βb3·121βb2+1-1βb·1βb-12,βb=1.6684(19)Mwb=wzμsin⁡α+cos⁡α0.289P·cos⁡α+μ·sin⁡α-0.4165P-0.289P·tan⁡αsin⁡α-μ·cos⁡α(20)δ2b=0.51·61+vcos⁡α+μsin⁡αcot⁡α5Eb·wzμsin⁡α+cos⁡α(21)δ3b=wzμsin⁡α+cos⁡α·12c1-v2πEba2·0.289P·cos⁡α+μ·sin⁡α-0.4165P-0.289P·tan⁡αsin⁡α-μ·cos⁡α(22)δ4b=1-νbdp2EbPsin⁡α-μcos⁡αtan⁡α·wzμsin⁡α+cos⁡α(23)δ5b=1.8449wzμsin⁡α+cos⁡α·21-v2·cos⁡α+μsin⁡απEbFor ISO internal threads, the relationship betweena,b,c, and pitch P is(24)a=0.875Pb=0.5Pc=0.325PSubstituting (24) into (9), (10), (11), (12), (13), (15), and (16) type one gets the relation(25)δ1n=3.784MwnEnP2βn2·121βn2+1-1βn·1βn-12,βn=1.751(26)Mwn=wzμsin⁡α+cos⁡α0.325P·cos⁡α+μ·sin⁡α-0.4375P-0.325P·tan⁡αsin⁡α-μ·cos⁡α(27)δ2n=0.55962·61+vcos⁡α+μsin⁡αcot⁡α5Eb·wzμsin⁡α+cos⁡α(28)δ3n=wzμsin⁡α+cos⁡α·12c1-v2πEba2·0.325P·cos⁡α+μ·sin⁡α-0.4375P-0.325P·tan⁡αsin⁡α-μ·cos⁡α(29)δ4n=D02+dp2D02-dp2+νndp2EnP·sin⁡α-μcos⁡αtan⁡α·wzμsin⁡α+cos⁡α(30)δ5n=1.7928wzμsin⁡α+cos⁡α·21-v2·cos⁡α+μsin⁡απEbBy adding these deformations separately, the total deformation (shown in Figure4) of screw thread and nut thread can be obtained under the action of force wz.(31)δb=δ1b+δ2b+δ3b+δ4b+δ5b(32)δn=δ1n+δ2n+δ3n+δ4n+δ5nThe unit force per unit width of the axial direction can be expressed as(33)fΔ=wzwz=1Under the action of unit force of axial unit width, the total deformation of external thread and internal thread is(34)δb1=δ1b+δ2b+δ3b+δ4b+δ5bwzand(35)δn1=δ1n+δ2n+δ3n+δ4n+δ5nwzFor threaded connections, at thex-axis of the loadF, the axial deformation of screws and nuts can be expressed as(36)zbr=δb1·∂F∂r(37)znr=δn1·∂F∂rHerer is the length along the helical direction, and the relation between the axial height x and the length along the helix direction can be represented by the following formula according to the geometric relation shown in Figure 6.(38)r=xsin⁡βHere,β is the lead angle of the thread shown in Figure 6, and then(39)znx=δn1·∂F∂r=δn1·∂F∂x·∂x∂r=δn1·sin⁡β·∂F∂x(40)zbx=δb1·∂F∂r=δb1·∂F∂x·∂x∂r=δb1·sin⁡β·∂F∂xAssume(41)∂F∂x·zbx=1δb1·sin⁡β=kubxx(42)∂F∂x·znx=1δn1·sin⁡β=kunxxHere,kbx(x) and knx(x) represent the stiffness of the unit axial length of the nut and the screw, respectively, for the unit force.The axial total deformation of the threaded connection atx is denoted as(43)zxx=zbx+znxThe stiffness of the unit axial length of the threaded connection is expressed as(44)kuxx=∂F∂x·zxx=∂F∂x÷δb1·sin⁡β·∂F∂x+δn1·sin⁡β·∂F∂x=1δn1+δb1·sin⁡βAs shown in Figure1(a), the threaded connection structure includes a nut body and a screw body. The nut is fixed, the screw is subjected to pulling force, the total axial force is Fb, and the axial force at the threaded connection screw x isF(x). If the position of the bottom end face of the nut is the origin 0 and, at the x position, the axial force isF(x), the screw elongation amount εb and the nut compression εn can be obtained from the following:(45)εbx=FxAbxEb(46)εnx=FxAnxEnFigure 1
Thread force.
(a) (b)Figure 2
Schematic diagram of the article.Figure 3
Thread deformation caused by various reasons.
(a) (b) (c) (d) (e)Figure 4
The comprehensive deformation of the thread.
(a) (b)Figure 5
The force on the thread.Figure 6
The lead angle of the thread.whereAb(x) and An(x) are the vertical cross-sectional areas of screws and nuts at the x position. Eb and En are, respectively, Young’s modulus of the screw body and Young’s modulus of the nut body. Find the displacement gradient for the expression, which is, respectively, expressed as(47)∂zbx∂x=1kubxx·∂2F∂x2(48)∂znx∂x=1kunxx·∂2F∂x2Here,kubxx=1/δb1·sin⁡β, and kunxx=1/δn1·sin⁡β.As shown in Figure1(a), the screw is subjected to the tensile force Fb, with the bottom of the nut as the coordinate origin, and the force at the x position is Fx, and then the elongation of the screw at x is ∫xLεbxdx=wb, and the compressed shortening amount of the nut at x is ∫xLεnxdx=wn. The relationship between wb, wn, zb, and zn is ∫xLεbxdx+∫xLεnxdx=zbx+znxx=L-zbx+znxx=x (see Figures 7 and 1(a),), and the partial derivative of this relation can be obtained by the following formula:(49)εbx+εnx=∂zbx∂x+∂znx∂xFigure 7
Illustration of the elastic deformation of the screwed portion of the threaded connection.Substituting (45), (46), (47), and (48) into (49) and simplifying it(50)AnEn+AbEbAbEbAnEn·kubxkunxkubx+kunxFx=∂2Fx∂x2Let(51)n=AnEn+AbEbAbEbAnEn·kubxkunxkubx+kunxThen(52)n2Fx=∂2Fx∂x2From mathematical knowledge, the equation is a differential equation. The general solution of the equation can be expressed as(53)Fx=C1sinh⁡nx+C2cosh⁡nxAs can be seen from Figure1, the axial force at the first thread at the connection surface of the nut and the screw is Fb, and the axial force at the last thread at the lower end of the thread joint surface of the nut and screw is 0; that is, the boundary condition isF(x=0)=Fb andF(x=L)=0. Taking these boundary conditions into the equation will give C1=-Fb(cosh⁡(nL))/sinh⁡(nL) and C2=Fb, so we get the expression of the threaded connection axial load as(54)Fx=Fbcosh⁡nx-cosh⁡nLsinh⁡nLsinh⁡nxTherefore, the axial force distribution density of the thread connection along thex direction can be expressed as(55)fx=FxFb=cosh⁡nx-cosh⁡nLsinh⁡nLsinh⁡nx


## 2.1.1. Bending Deformation

In the threaded connection, under the action of the load, the contact surface friction coefficient isμ, when the sliding force along the inclined plane is greater than the friction force along the inclined plane, the relative sliding occurs between the two inclined planes, and the axial unit width force (shown in Figure 3) is wz; if the influence of the lead angle is ignored, the force per unit width perpendicular to the thread surface can be expressed as(1)w=wzμsin⁡α+cos⁡αThe force per unit width perpendicular to the thread surface can be decomposed into thex-direction component force and they-direction component force, respectively(2)wcos⁡α=wzcos⁡αμsin⁡α+cos⁡αand(3)wsin⁡α=wzsin⁡αμsin⁡α+cos⁡αThe friction generated along the slope iswμ; i.e.,(4)wμ=wzμμsin⁡α+cos⁡αThe forcewμ is also decomposed intox-direction force andy-direction force, which are μsin⁡α and μcos⁡α, respectively.(5)wμsin⁡α=wzμsin⁡αμsin⁡α+cos⁡α(6)wμcos⁡α=wzμcos⁡αμsin⁡α+cos⁡αIn the unit width, the thread is regarded as a rectangular variable-section cantilever beam. Under the action of the above-mentioned force, the thread undergoes bending deformation, and the virtual work done by the bending momentE- on the beam dy section is(7)E-dθ=E-MwEbIydyAccording to the principle of virtual work, the deflectionδ1 (see Figure 3(a)) of the beam subjected to the load is(8)δ1=∫0cE-MwEbIydywhereE- is the bending moment of the unit load beam. Mw is the bending moment of the beam under the actual load.I(y) is the area moment of inertia of the beam at y. Eb is Young’s Modulus of the material.c is the length of the beam. Here, the forces are assumed as acting on the mean diameter of the thread.As shown in Figure5, the height h(y) of the beam section per unit width and the area moment of inertia I(y) of the section can be expressed by using the function interpolation.(9)hy=1+β1-1c-ych,0≤y≤cIy=112bh31+β1-1c-yc3,0≤y≤cβ1=ahwhereh is the beam end section height;b is the beam section width; β1 is the beam root section height and the beam end section height ratio; see Figure 5.From Figure5, the bending moment of the beam is related to they-axis component of w andwμ, and these components cause the beam to bend; therefore, the analytical solution shows that the bending moment of the unit width beam subjected to the friction force and the vertical load of the thread surface is(10)Mw=wzμsin⁡α+cos⁡α·cos⁡α·c+μsin⁡α·c-sin⁡αa2-ctan⁡α+μcos⁡αa2-ctan⁡αSubstituting (10) and (9) to (8) and integrating to obtain the analytical expression of the deflection δ1 (shown in Figure 3(a)) of the cantilever beam with variable cross-section under load one has(11)δ1=12Mwc2Ebbh3·121β12+1-1β1·1β1-12


## 2.1.2. Shear Deformation

Assume that the distribution of shear stress on any section is distributed according to the parabola [31] and the deformation δ2 (see Figure 3(b)) caused by the shear force within the width of unit 1 is(12)δ2=logeab·61+vcos⁡α+μsin⁡αcot⁡α5Eb·wzμsin⁡α+cos⁡α


## 2.1.3. Inclination Deformation of the Thread Root

Under the action of the load, the thread surface is subjected to a bending moment, and the root of the thread is tilted, as shown in Figure3(c). Due to the inclination of the thread, axial displacement occurs at the point of action of the thread surface force, and the axial displacement can be expressed as [31](13)δ3=wzμsin⁡α+cos⁡α·12c1-v2πEba2·cos⁡α·c+μsin⁡α·c-sin⁡αa2-ctan⁡α+μcos⁡αa2-ctan⁡α


## 2.1.4. Deformation due to Radial Expansion and Radial Shrinkage

According to the static analysis, the thread is subjected to radial forcewsin⁡α-wμcos⁡α (shown in Figure 4), and it is known from the literature [31] that the internal and external thread radial deformation (shown in Figure 3(d)) are(14)δ4b=1-νbdp2EbPtan⁡α·wzsin⁡α-μcos⁡αμsin⁡α+cos⁡αand(15)δ4n=D02+dp2D02-dp2+νndp2EnP·tan⁡α·wzsin⁡α-μcos⁡αμsin⁡α+cos⁡α


## 2.1.5. Shear Deformation of the Root

Assuming that the shear stress of the root section is evenly distributed, the displacement of theO point in the x direction caused by the shear deformation (shown in Figure 3(e)) is the same as the displacement of the thread in the x direction; this displacement can be expressed as [31](16)δ5=wzμsin⁡α+cos⁡α·21-v2·cos⁡α+μsin⁡απEb·PalogeP+a/2P-a/2+12loge4P2a2-1For ISO internal threads, the relationship between a, b, c, and pitchP is(17)a=0.833Pb=0.5Pc=0.289PSubstituting (17) into (9), (10), (11), (12), (13), (14), and (16) one gets the relation(18)δ1b=3.468MwbEbP2βb3·121βb2+1-1βb·1βb-12,βb=1.6684(19)Mwb=wzμsin⁡α+cos⁡α0.289P·cos⁡α+μ·sin⁡α-0.4165P-0.289P·tan⁡αsin⁡α-μ·cos⁡α(20)δ2b=0.51·61+vcos⁡α+μsin⁡αcot⁡α5Eb·wzμsin⁡α+cos⁡α(21)δ3b=wzμsin⁡α+cos⁡α·12c1-v2πEba2·0.289P·cos⁡α+μ·sin⁡α-0.4165P-0.289P·tan⁡αsin⁡α-μ·cos⁡α(22)δ4b=1-νbdp2EbPsin⁡α-μcos⁡αtan⁡α·wzμsin⁡α+cos⁡α(23)δ5b=1.8449wzμsin⁡α+cos⁡α·21-v2·cos⁡α+μsin⁡απEbFor ISO internal threads, the relationship betweena,b,c, and pitch P is(24)a=0.875Pb=0.5Pc=0.325PSubstituting (24) into (9), (10), (11), (12), (13), (15), and (16) type one gets the relation(25)δ1n=3.784MwnEnP2βn2·121βn2+1-1βn·1βn-12,βn=1.751(26)Mwn=wzμsin⁡α+cos⁡α0.325P·cos⁡α+μ·sin⁡α-0.4375P-0.325P·tan⁡αsin⁡α-μ·cos⁡α(27)δ2n=0.55962·61+vcos⁡α+μsin⁡αcot⁡α5Eb·wzμsin⁡α+cos⁡α(28)δ3n=wzμsin⁡α+cos⁡α·12c1-v2πEba2·0.325P·cos⁡α+μ·sin⁡α-0.4375P-0.325P·tan⁡αsin⁡α-μ·cos⁡α(29)δ4n=D02+dp2D02-dp2+νndp2EnP·sin⁡α-μcos⁡αtan⁡α·wzμsin⁡α+cos⁡α(30)δ5n=1.7928wzμsin⁡α+cos⁡α·21-v2·cos⁡α+μsin⁡απEbBy adding these deformations separately, the total deformation (shown in Figure4) of screw thread and nut thread can be obtained under the action of force wz.(31)δb=δ1b+δ2b+δ3b+δ4b+δ5b(32)δn=δ1n+δ2n+δ3n+δ4n+δ5nThe unit force per unit width of the axial direction can be expressed as(33)fΔ=wzwz=1Under the action of unit force of axial unit width, the total deformation of external thread and internal thread is(34)δb1=δ1b+δ2b+δ3b+δ4b+δ5bwzand(35)δn1=δ1n+δ2n+δ3n+δ4n+δ5nwzFor threaded connections, at thex-axis of the loadF, the axial deformation of screws and nuts can be expressed as(36)zbr=δb1·∂F∂r(37)znr=δn1·∂F∂rHerer is the length along the helical direction, and the relation between the axial height x and the length along the helix direction can be represented by the following formula according to the geometric relation shown in Figure 6.(38)r=xsin⁡βHere,β is the lead angle of the thread shown in Figure 6, and then(39)znx=δn1·∂F∂r=δn1·∂F∂x·∂x∂r=δn1·sin⁡β·∂F∂x(40)zbx=δb1·∂F∂r=δb1·∂F∂x·∂x∂r=δb1·sin⁡β·∂F∂xAssume(41)∂F∂x·zbx=1δb1·sin⁡β=kubxx(42)∂F∂x·znx=1δn1·sin⁡β=kunxxHere,kbx(x) and knx(x) represent the stiffness of the unit axial length of the nut and the screw, respectively, for the unit force.The axial total deformation of the threaded connection atx is denoted as(43)zxx=zbx+znxThe stiffness of the unit axial length of the threaded connection is expressed as(44)kuxx=∂F∂x·zxx=∂F∂x÷δb1·sin⁡β·∂F∂x+δn1·sin⁡β·∂F∂x=1δn1+δb1·sin⁡βAs shown in Figure1(a), the threaded connection structure includes a nut body and a screw body. The nut is fixed, the screw is subjected to pulling force, the total axial force is Fb, and the axial force at the threaded connection screw x isF(x). If the position of the bottom end face of the nut is the origin 0 and, at the x position, the axial force isF(x), the screw elongation amount εb and the nut compression εn can be obtained from the following:(45)εbx=FxAbxEb(46)εnx=FxAnxEnFigure 1
Thread force.
(a) (b)Figure 2
Schematic diagram of the article.Figure 3
Thread deformation caused by various reasons.
(a) (b) (c) (d) (e)Figure 4
The comprehensive deformation of the thread.
(a) (b)Figure 5
The force on the thread.Figure 6
The lead angle of the thread.whereAb(x) and An(x) are the vertical cross-sectional areas of screws and nuts at the x position. Eb and En are, respectively, Young’s modulus of the screw body and Young’s modulus of the nut body. Find the displacement gradient for the expression, which is, respectively, expressed as(47)∂zbx∂x=1kubxx·∂2F∂x2(48)∂znx∂x=1kunxx·∂2F∂x2Here,kubxx=1/δb1·sin⁡β, and kunxx=1/δn1·sin⁡β.As shown in Figure1(a), the screw is subjected to the tensile force Fb, with the bottom of the nut as the coordinate origin, and the force at the x position is Fx, and then the elongation of the screw at x is ∫xLεbxdx=wb, and the compressed shortening amount of the nut at x is ∫xLεnxdx=wn. The relationship between wb, wn, zb, and zn is ∫xLεbxdx+∫xLεnxdx=zbx+znxx=L-zbx+znxx=x (see Figures 7 and 1(a),), and the partial derivative of this relation can be obtained by the following formula:(49)εbx+εnx=∂zbx∂x+∂znx∂xFigure 7
Illustration of the elastic deformation of the screwed portion of the threaded connection.Substituting (45), (46), (47), and (48) into (49) and simplifying it(50)AnEn+AbEbAbEbAnEn·kubxkunxkubx+kunxFx=∂2Fx∂x2Let(51)n=AnEn+AbEbAbEbAnEn·kubxkunxkubx+kunxThen(52)n2Fx=∂2Fx∂x2From mathematical knowledge, the equation is a differential equation. The general solution of the equation can be expressed as(53)Fx=C1sinh⁡nx+C2cosh⁡nxAs can be seen from Figure1, the axial force at the first thread at the connection surface of the nut and the screw is Fb, and the axial force at the last thread at the lower end of the thread joint surface of the nut and screw is 0; that is, the boundary condition isF(x=0)=Fb andF(x=L)=0. Taking these boundary conditions into the equation will give C1=-Fb(cosh⁡(nL))/sinh⁡(nL) and C2=Fb, so we get the expression of the threaded connection axial load as(54)Fx=Fbcosh⁡nx-cosh⁡nLsinh⁡nLsinh⁡nxTherefore, the axial force distribution density of the thread connection along thex direction can be expressed as(55)fx=FxFb=cosh⁡nx-cosh⁡nLsinh⁡nLsinh⁡nx


## 2.2. Thread Connection Stiffness

The stiffness in the axial directionx of the bolted connection is equal to the axial force distribution of the threaded connection multiplied by the unit stiffness; i.e.,(56)kxx=fxkuxxThe overall stiffness of the bolt connection can be expressed as(57)Kc=∫0Lkxx=∫0LfxkuxxdxSubstituting (44) and (55) into (57), the stiffness of the bolt connection is expressed as(58)Kc=∫0Lkxx=∫0Lfxkuxxdx=1nδb1+δn1·sin⁡β·cosh⁡nL-1sinh⁡nL


## 3. FEA Model

A 3D finite element model (shown in Figure10) was established, and FEA was performed to analyze the influence of various parameters of the thread on the thread stiffness. These parameters include material, thread length, pitch, etc.The FEA software ANSYS 14.0 was used for analysis. During the analysis, the end face of the nut was fixed (shown in Figure9), the initial state of the model is shown in Figure 8, and an axial displacement Δx was forced to the end face of the screw. Then, the axial force Fx of the screw end face was extracted. The axial stiffness of the threaded connection was calculated by the FEM. The friction coefficient μ of the thread contact surface is set first. In FEA, the contact algorithm used is Augmented Lagrange. Figures 11(a)–11(c) are the force convergence curves for FEA of threaded connections. Figures 12(a)–12(c) are the effect of the reciprocal of the mesh size on the axial force obtained by FEA. We can see from Figures 12(a)–12(c) that as the mesh size decreases, the resulting axial force gradually decreases, but when the mesh size is small to a certain extent, the resulting axial force will hardly decrease. The axial force at this time is the axial force required by the author. With known displacements and axial force, the stiffness of the threaded connection can be calculated using the formula (59)Kc=FxΔxFigure 8
Initial state.Figure 9
Threaded connection deformed by axial forces.Figure 10
Finite element meshing of thread-bonded 3D finite element models.
(a)
Nut (b)
ScrewFigure 11
FEA convergence curve.
(a)
M10×1.5×9,μ=0.08, Δx=0.001mm (b)
M8×1.25×6.5,μ=0.08, Δx=0.001mm (c)
M6×1.0×6.1,μ=0.08, Δx=0.001mmFigure 12
Influence of the reciprocal of finite element mesh size on axial force.
(a)
M10×1.5×9,μ=0.08, Δx=0.001mm (b)
M8×1.25×6.5,μ=0.08, Δx=0.001mm (c)
M6×1.0×6.1, u=0.08,Δx=0.001mm


## 4. Tensile Test of Threaded Connections [1]

In order to verify the effectiveness of this paper method, the experimental data of the experimental device in [1] are used. In [1], the electronic universal testing machine is used to measure the load-defection data of samples, and the test sample is made of brass. The tension value Fxt can be read from the test machine. The axial deflection of thread connection can be represented by the displacement variation δL between two lines as shown in Figure 14, which can be measured by a video gauge [1]. In [1], in order to obtain the most accurate data possible, each size of the thread is in a small range of deformation during the tensile test, and each size of the thread tensile test is performed 10 times, and the average value is calculated as the final calculated data. Some samples in the experiment are shown in Figure 13.Figure 13
Part of experimental threaded connection samples and ISO internal thread and ISO external thread.Figure 14
Tensile test [1].The stiffness calculation formula is(60)Kc=FxtδL.The materials used to make nuts and screws are brass. Young’s modulus of brass is 107GPa, and Poisson’s ratio is 0.32 [1].


## 5. Results and Discussion

### 5.1. Stiffness of Threaded Connections

Croccolo, D. [12], Nassar SA [19], and Zou Q [32] studied the coefficient of friction of the thread. According to the study by Zou Q and Nassar SA, in the case of lubricating oil on the thread surface, the friction coefficient of the steel-steel thread connection thread is 0.08, and the friction coefficient of aluminum-aluminum thread connection thread is 0.1.In order to verify the correctness of the calculation results of the theory presented in this paper, a variety of threaded connections were used to calculate an experimental test.In the finite element analysis and theoretical calculations of this paper, Young’s modulus of steel isEb= En =200Gpa, and Poisson’s ratio of steel is 0.3, and the friction coefficient [12, 19, 32] is set to 0.08.The experimental data, FEA data, and Yamamoto method data in Tables1 and 5 are from literature [1]. As can be seen from Tables 1 and 5, the calculated values obtained in this paper are all higher than the experimental results. Perhaps the error is caused by the presence of a small amount of impurities on the surface of the thread and partial deformation of the thread inevitably and there is a slip between the threaded contact surfaces. The theoretical calculation results and FEA results in this paper have a small error.Table 1
Stiffness of threaded connections with different engaged lengths [1] (kN/mm).
No. Size code of threads Material Exp. Theory FEA This study Yamamoto method  1 M36×4×32 Brass 3627.6 4201.9 4282.1 3630.9 2 M36×4×20 2664.3 3152.9 3946.1 2761.6 3 M36×4×12 1801.3 2074.2 3129.3 1816.8In Table2, the effect of thread length on stiffness is presented. It can be seen that when the same nominal diameter M10, the same pitchP=1.5, and the same material steel are taken, when the thread engaged length is taken as 14 mm, 9 mm, and 6 mm, respectively, it is found that the longer the thread engaged length, the greater the stiffness and the smaller the length of the bond, the smaller the stiffness.Table 2
Stiffness of threaded connections with different engaged lengths (kN/mm).
No. Size code of threads Material Theory FEA This study Yamamoto method  1 M10×1.5×14 Steel 2095.26 2076.6 1965.21 2 M10×1.5×9 1759.23 2006.2 1770.02 3 M10×1.5×6 1354.70 1804.6 1551.27In the FEA, the method of this paper and the Yamamoto method, Young’s Modulus of aluminum alloy isE=68.9GPa; Poisson’s ratio of the aluminum alloy is 0.34. The friction coefficient of the steel- steel threaded connection [12, 21, 22, 32] is set to 0.08, and the friction coefficient of aluminum-aluminum threaded connection [12, 19, 32] is set to 0.1. Table 3 shows the effect of different materials on the stiffness of threaded connections. The two types of threaded connections are made of two different materials, the steel and aluminum alloys. It can be seen from the table that, under the condition of the same pitch, the same nominal diameter, and the same engaged length, the stiffness of the steel thread connection is larger than that when the material is aluminum.Table 3
Stiffness of threaded connections with different material (kN/mm).
No. Size code of thread Material Theory FEA This study Yamamoto method  1 M10×1.5×9 Steel 1759.23 2006.2 1770.02 2 M10×1.5×9 Aluminum alloy 607.51 682.11 597.99In Table4, it also shows the influence of different pitches on the stiffness of the thread connection. It can be seen that with the same engaged length, the same material, and the same nominal diameter, the pitch P is 1.5, 1.25, and 1, respectively, and we find that the smaller the pitch, the greater the stiffness.Table 4
Stiffness of threaded connections with different pitch (kN/mm).
No. Size code of threads Material Theory This study Yamamoto method  1 M10×1.5×9 Steel 1759.23 2006.2 2 M10×1.25×9 2083.17 2177.1 3 M10×1×9 2324.30 2375.2Table 5
Stiffness of threaded connections with different engaged lengths [1] (kN/mm).
No. Size code of threads Material Exp. Theory This study Yamamoto method  1 M36×3×12 Brass 2085.3 2593.2 3525.1 2 M36×2×12 Brass 2396.2 3418.3 4008.9When using FEM to analyze the influence of friction factors on the stiffness of threaded connections, the thread specification is M6×1×6.1 and the friction coefficients are 0.01, 0.05, 0.1, 0.2, 0.25, and 0.3. (as shown in Figures15–20)Figure 15
Axial displacement for screws with a friction coefficient ofμ=0.01. Total force reaction=1086.00 N.Figure 16
Axial displacement for screws with a friction coefficient ofμ=0.05. Total force reaction=1092.00 N.Figure 17
Axial displacement for screws with a friction coefficient ofμ=0.1. Total force reaction=1098.00 N.Figure 18
Axial displacement for screws with a friction coefficient ofμ=0.2. Total force reaction=1108.40 N.Figure 19
Axial displacement for screws with a friction coefficient ofμ=0.25. Total force reaction=1110.6 N.Figure 20
Axial displacement for screws with a friction coefficient ofμ=0.3. Total force reaction=1114.2 N.Figures21 and 22 show the results of stiffness calculations. The thread size is M10×1.5×9 and M6×1×6.1, the material is steel, Poisson’s ratio of the material is 0.3, Young’s Modulus of the material is 200 GPa, and the thread surface friction coefficient is taken as 0.01, 0.05, 0.1, 0.2, 0.25, and 0.3. Calculated using the theories of this paper, FEA and Yamamoto, respectively, and from Figures 21 and 22, we can see that the results of FEA are very similar to the results of the theoretical calculations of this paper, the variation trend of stiffness with friction coefficient is the same, and it increases with the increase of friction coefficient, and the results of the FEA are in good agreement with those of the FEA; however, Yamamoto theory does not consider the influence of the friction coefficient on stiffness, and this is obviously unreasonable.Figure 21
Effect of friction coefficient on stiffness. M10×1.5×9.Figure 22
Effect of friction coefficient on stiffness. M6×1×6.1.


### 5.2. Effect of Friction Coefficient on Axial Force Distribution

Take the thread size as M6×0.75×6.1, the axial loadFb is taken as 100N, 350N, and 550N, respectively, and take the friction coefficients 0, 0.3, 0.6, and 1, respectively, to calculate the axial force distribution of the thread. As can be seen from Figure 23, when the friction coefficient is 1, the curve bending degree is the greatest, when the friction coefficient is 0, the curve bending degree is the lightest, the curve bending degree is greater, indicating that the more uneven the distribution of axial force, the smaller the curve bending degree, indicating that the more uniform the distribution of axial force. We can see that the friction coefficient of thread surface has an effect on the distribution of axial force.Figure 23
Effect of Friction Coefficient on Axial Force Distribution.
(a)
F b=100N (b)
F b=350N (c)
F b=550N


## 5.1. Stiffness of Threaded Connections

Croccolo, D. [12], Nassar SA [19], and Zou Q [32] studied the coefficient of friction of the thread. According to the study by Zou Q and Nassar SA, in the case of lubricating oil on the thread surface, the friction coefficient of the steel-steel thread connection thread is 0.08, and the friction coefficient of aluminum-aluminum thread connection thread is 0.1.In order to verify the correctness of the calculation results of the theory presented in this paper, a variety of threaded connections were used to calculate an experimental test.In the finite element analysis and theoretical calculations of this paper, Young’s modulus of steel isEb= En =200Gpa, and Poisson’s ratio of steel is 0.3, and the friction coefficient [12, 19, 32] is set to 0.08.The experimental data, FEA data, and Yamamoto method data in Tables1 and 5 are from literature [1]. As can be seen from Tables 1 and 5, the calculated values obtained in this paper are all higher than the experimental results. Perhaps the error is caused by the presence of a small amount of impurities on the surface of the thread and partial deformation of the thread inevitably and there is a slip between the threaded contact surfaces. The theoretical calculation results and FEA results in this paper have a small error.Table 1
Stiffness of threaded connections with different engaged lengths [1] (kN/mm).
No. Size code of threads Material Exp. Theory FEA This study Yamamoto method  1 M36×4×32 Brass 3627.6 4201.9 4282.1 3630.9 2 M36×4×20 2664.3 3152.9 3946.1 2761.6 3 M36×4×12 1801.3 2074.2 3129.3 1816.8In Table2, the effect of thread length on stiffness is presented. It can be seen that when the same nominal diameter M10, the same pitchP=1.5, and the same material steel are taken, when the thread engaged length is taken as 14 mm, 9 mm, and 6 mm, respectively, it is found that the longer the thread engaged length, the greater the stiffness and the smaller the length of the bond, the smaller the stiffness.Table 2
Stiffness of threaded connections with different engaged lengths (kN/mm).
No. Size code of threads Material Theory FEA This study Yamamoto method  1 M10×1.5×14 Steel 2095.26 2076.6 1965.21 2 M10×1.5×9 1759.23 2006.2 1770.02 3 M10×1.5×6 1354.70 1804.6 1551.27In the FEA, the method of this paper and the Yamamoto method, Young’s Modulus of aluminum alloy isE=68.9GPa; Poisson’s ratio of the aluminum alloy is 0.34. The friction coefficient of the steel- steel threaded connection [12, 21, 22, 32] is set to 0.08, and the friction coefficient of aluminum-aluminum threaded connection [12, 19, 32] is set to 0.1. Table 3 shows the effect of different materials on the stiffness of threaded connections. The two types of threaded connections are made of two different materials, the steel and aluminum alloys. It can be seen from the table that, under the condition of the same pitch, the same nominal diameter, and the same engaged length, the stiffness of the steel thread connection is larger than that when the material is aluminum.Table 3
Stiffness of threaded connections with different material (kN/mm).
No. Size code of thread Material Theory FEA This study Yamamoto method  1 M10×1.5×9 Steel 1759.23 2006.2 1770.02 2 M10×1.5×9 Aluminum alloy 607.51 682.11 597.99In Table4, it also shows the influence of different pitches on the stiffness of the thread connection. It can be seen that with the same engaged length, the same material, and the same nominal diameter, the pitch P is 1.5, 1.25, and 1, respectively, and we find that the smaller the pitch, the greater the stiffness.Table 4
Stiffness of threaded connections with different pitch (kN/mm).
No. Size code of threads Material Theory This study Yamamoto method  1 M10×1.5×9 Steel 1759.23 2006.2 2 M10×1.25×9 2083.17 2177.1 3 M10×1×9 2324.30 2375.2Table 5
Stiffness of threaded connections with different engaged lengths [1] (kN/mm).
No. Size code of threads Material Exp. Theory This study Yamamoto method  1 M36×3×12 Brass 2085.3 2593.2 3525.1 2 M36×2×12 Brass 2396.2 3418.3 4008.9When using FEM to analyze the influence of friction factors on the stiffness of threaded connections, the thread specification is M6×1×6.1 and the friction coefficients are 0.01, 0.05, 0.1, 0.2, 0.25, and 0.3. (as shown in Figures15–20)Figure 15
Axial displacement for screws with a friction coefficient ofμ=0.01. Total force reaction=1086.00 N.Figure 16
Axial displacement for screws with a friction coefficient ofμ=0.05. Total force reaction=1092.00 N.Figure 17
Axial displacement for screws with a friction coefficient ofμ=0.1. Total force reaction=1098.00 N.Figure 18
Axial displacement for screws with a friction coefficient ofμ=0.2. Total force reaction=1108.40 N.Figure 19
Axial displacement for screws with a friction coefficient ofμ=0.25. Total force reaction=1110.6 N.Figure 20
Axial displacement for screws with a friction coefficient ofμ=0.3. Total force reaction=1114.2 N.Figures21 and 22 show the results of stiffness calculations. The thread size is M10×1.5×9 and M6×1×6.1, the material is steel, Poisson’s ratio of the material is 0.3, Young’s Modulus of the material is 200 GPa, and the thread surface friction coefficient is taken as 0.01, 0.05, 0.1, 0.2, 0.25, and 0.3. Calculated using the theories of this paper, FEA and Yamamoto, respectively, and from Figures 21 and 22, we can see that the results of FEA are very similar to the results of the theoretical calculations of this paper, the variation trend of stiffness with friction coefficient is the same, and it increases with the increase of friction coefficient, and the results of the FEA are in good agreement with those of the FEA; however, Yamamoto theory does not consider the influence of the friction coefficient on stiffness, and this is obviously unreasonable.Figure 21
Effect of friction coefficient on stiffness. M10×1.5×9.Figure 22
Effect of friction coefficient on stiffness. M6×1×6.1.


## 5.2. Effect of Friction Coefficient on Axial Force Distribution

Take the thread size as M6×0.75×6.1, the axial loadFb is taken as 100N, 350N, and 550N, respectively, and take the friction coefficients 0, 0.3, 0.6, and 1, respectively, to calculate the axial force distribution of the thread. As can be seen from Figure 23, when the friction coefficient is 1, the curve bending degree is the greatest, when the friction coefficient is 0, the curve bending degree is the lightest, the curve bending degree is greater, indicating that the more uneven the distribution of axial force, the smaller the curve bending degree, indicating that the more uniform the distribution of axial force. We can see that the friction coefficient of thread surface has an effect on the distribution of axial force.Figure 23
Effect of Friction Coefficient on Axial Force Distribution.
(a)
F b=100N (b)
F b=350N (c)
F b=550N


## 6. Conclusion

This study provides a new method of calculating the thread stiffness considering the friction coefficient and analyzes the influence of the thread geometry and material parameters on the thread stiffness and also analyzes the influence of the friction coefficient on the thread stiffness and axial force distribution.(1)
The results of the calculation of the thread stiffness calculated by the theoretical calculation method of this study are basically consistent with the results of the FEA. The results obtained by the test are smaller than the calculated results. This is due to the influence of the thread manufacturing on the experimental results.(2)
Thread-stiffness is closely related to material properties, pitch, and thread length. We can obtain higher stiffness by increasing Young’s modulus of the material, increasing the length of the thread, and reducing the pitch.(3)
We can also increase the friction coefficient of the thread joint surface to increase the stiffness of the thread connection, but we have found that using this method to increase the thread stiffness is limited.(4)
In order to make the axial load distribution of the thread uniform, we can reduce the friction coefficient of the thread surface, but we found that the use of this method to improve the distribution of the axial force of the thread has limited effectiveness.



---

*Source: 8424283-2019-01-23.xml*