Loading and contact stress distribution on the thread teeth in tubing and casing premium threaded connections are of great importance for design optimization, pretightening force control, and thread failure prevention. This paper proposes an analytical method based on the elastic mechanics. This is quite different from other papers, which mainly rely on finite element analysis. The differential equation of load distribution on the thread teeth was established according to equal pitch of the engaged thread after deformation and solved by finite difference method. Furthermore, the relation between load acting on each engaged thread and mean contact stress on its load flank is set up based on the geometric description of thread surface. By comparison, this new analytical method with the finite element analysis for a modified API 177.8 mm premium threaded connection is approved. Comparison of the contact stress on the last engaged thread between analytical model and FEM shows that the accuracy of analytical model will decline with the increase of pretightening force after the material enters into plastic deformation. However, the analytical method can meet the needs of engineering to some extent because its relative error is about 6.2%~18.1% for the in-service level of pretightening force.
1. Introduction
Premium threaded connections have been widely used for maintaining structural and sealing integrity of tubing and casing strings in HPHT wells. The overwhelming majority of premium threaded connections are characterized by a tapered metal-to-metal seal structure for leakage resistance and an interior rotary torque shoulder to control final make-up torques, providing additional sealing performance (shown in Figure 1). These characteristics are different from API buttress thread connections and extreme-line thread connections [1–3]. In general, the thread teeth of premium threaded connection are used to bear axial tensile load, due to the load flank angle α design (shown in Figure 3) of usually -4°~3°. We also noticed that the taper angle of metal-to-metal seal is usually 1.79°~2.86°, so the axial component of sealing contact stress could be neglected and the action force of sealing and shoulder replaced approximatively by an axial pretightening force Ft acting on the A-B cross section plane (shown in Figure 1). Consequently, it is significant to investigate the load and contact stress on thread teeth caused by this pretightening force Ft for design optimization, pretightening torque control, and thread failure prevention.
The schematic of a premium threaded connection.
The load distribution and contact stress on thread teeth were investigated in many previous references in the last century, whose methods can mainly be divided into three categories. The first category is developed on the basis of analytical method. Maduschka [4] and Sopwith [5] first proposed the models of the load distribution in cylindrical screw threads, whose results clearly indicate that the first three or four engaged threads carry more than half of the preload induced by the make-up torque. Taking axial and tangential forces, as well as bending moments, into consideration, Yazawa and Hongo [6] developed a model to investigate the load distribution of a bolt-nut connection. Wang and Marshek [7] also established a spring model to predict the load distribution of bolt-nut connection. Furthermore, Bruschelli and Latorrata [8] developed an analytical model to study the load distribution for conical threads connections. Chen et al. [9] presented a practical and convenient model for cylindrical pipe threaded connection, which can be applied to calculate the load and deformation on each thread tooth just with tightening torque and thread numbers. The vast majority of the analytical models mentioned above are established from screw threaded connection or cylindrical pipe threaded connection and therefore are not applicable to premium connection: its tapered thread connection type contains special metal-to metal seal, interior rotary torque shoulder, and the specialized tooth form of buttress.
The second category method is developed based on finite element method (FEM). Capable of considering material nonlinear behavior and complex boundary conditions easily, FEM makes studying the load distribution, contact stress, and stress concentration on the thread tooth possible. O’Hara [10] established a 2D axisymmetric model of thread connection based on thread characteristic lengths and the numerical results were generally validated by Heywood’s equations. Bretl and Cook [11] proposed a unique FEM technique for threaded connections and their simulated results agreed with those from analytical and experimental methods both for conventional and for tapered threads. Zhao [12, 13] developed an elastic-plastic FEM and investigated the load distribution and SCFs (stress concentration factors) of the bolt-nut connection. Macdonald and Deans [14] acquired the SCFs and induced stresses on the drilling string thread and they found out that the high local stresses are located at the root of the first engaged tooth in the pin and the last engaged tooth in the box. Bahai et al. [15] also calculated the SCFs of API drill string threaded connector subjected to preloading, axial loading, and bending loads. Chen and Shih [16] developed a 3D FEM for bolted joints and analyzed the helical and friction effect on the load distribution of each thread. Baragetti [17] proposed a numerical finite element procedure for rotary shouldered connections (RSCs) and quantified the effects of taper on the loads carried by each engaged thread and contact stress on the thread flanks. Fukuoka and Nomura [18] established the FE model with accurate thread geometry by deriving a mathematical expression of the helical thread geometry of a single-thread screw. The numerical results indicate that axial load along engaged threads shows a different distribution pattern from those obtained by axisymmetric FE analysis and elastic theory. Although FE models have incomparable advantages for threaded connections analysis, most of them are based on 2D axisymmetric assumptions and the helix angle is usually ignored.
The third category method is developed based on experimental testing. Stoeckly and Macke [19] developed a testing apparatus to measure the axial displacement of the threads and compare their experimental data by modifying Sopwith equations. They found out that the frictional coefficient between the lubricated bolt and nut has little effect on the thread load distribution. Yuan et al. [20] tested the make-up and break-out operations in oilfields by placing the foil gauges on the inner surface of the male thread and the outer surface of the female thread along the axial and tangential direction and concluded that those operations have a strong impact on the service life of oil tubing threaded connection. Plácido et al. [21] carried out some full and reduced scale experimental tests on aluminum drill pipes to investigate their fatigue mechanism under cyclic bending and constant tensile loads. Slack et al. [22] proposed the ultrasonic reflection techniques to directly measure the contact stress between thread and sealing surface. It turned out to be an accurate and effective method to evaluate the sealing performance and to detect sealing damage in premium connections [23]. Though experimental testing can be an effective method for validating analytical models and FEM, those tests focused on limited aspects of threaded connection problems and are difficult for application in practice.
Based on elastic mechanics theory and some simplifications, this paper proposed an analytical method to calculate the load distribution and contact stress on each engaged thread tooth for premium threaded connections. The 2D FEM validated the analytical model, and the comparison of the contact stress on the last engaged thread between two methods investigated the accuracy of the analytical model and its practicality in engineering.
2. Model Development and Solution2.1. Loading Analysis on the Premium Threaded Connection
During tightening, the engaged threads of a premium connection will transfer the pretightening force Ft caused mainly by shoulder contact, thus leading to the compressed pin and the stretched coupling. To facilitate the analytical model, we assumed all thread teeth were intact and let the axis of premium threaded connection as x direction to establish a one-dimensional coordinate system OX with itsorigin being at the big end of the pin thread. Meanwhile, the engaged thread teeth were numbered 1~Nt consecutively from origin to the small end of the pin thread (shown in Figure 2).
The schematic of load distribution for premium connection.
The schematic of simplified cantilever beam model of an isosceles trapezoid for buttress type thread tooth.
The contact stresses between load flanks of engaged thread interact with each other as action and reaction, so the axial load distribution F(x) along x direction is equal within pin thread and coupling thread. Let f(x) denote the load distribution intensity; then
(1)fx=dFxdx.
The load distribution intensity f(x) equals the axial load acting on the engaged thread within unit axial length, which reflects the distribution characteristics of pretightening force along the engaged thread. It is obvious that axial load distribution F(x) can be calculated by the integral formula as follows:
(2)Fx=∫0xfxdx.
The axial load acting on each engaged thread can be further obtained as
(3)FiP=FiP-Fi-1P=∫i-1PiPfxdx,i=1,2,…,Nt.
2.2. Deformation of Thread Tooth in Premium Threaded Connection
The thread tooth profile of most premium connections typically adopts modified API buttress type thread, whose deformation can be calculated by simplifying itself into a cantilever beam model of an isosceles trapezoid [24] (shown in Figure 3).
The necessary parameters for calculating the elastic deformation of simplified thread tooth can be obtained from this basic geometry parameter of buttress type thread:
(4)a=P-Frscosγ,b=a-ctanα-γ+tanβ+γ,c=hBcosγ2.
According to [24], the deformations of simplified thread tooth are mainly caused by five reasons (shown in Figure 4). Now, let angle α-γ be equal to ϑ in Figure 4, and the thread tooth deformation δa caused by tooth bending can be expressed as
(5)δax=Nxcosα-γE31-ν24·1-2-ba2+2lnabctan3α-γhhhhh-4ca2tanα-γ.
The elastic deformations of simplified thread tooth (a) deformation caused by tooth bending or shear force; (b) deformation caused by tooth root incline; (c) deformation caused by shear deformation of tooth root; (d) deformation caused by radial effect.
The thread tooth deformation δb caused by shear force can be expressed as
(6)δbx=Nxcosα-γE61+ν5ctanα-γlnab.
The thread tooth deformation δc caused by tooth root incline can be expressed as
(7)δcx=Nxcosα-γE12c1-ν2πa2c-b2tanα-γ.
The thread tooth deformation δd caused by shear deformation of tooth root can be expressed as
(8)δdx=Nxcosα-γE21-ν2π·PalnP+0.5aP-0.5a+12ln4P2a2-1.
The deformations δrp, δrc caused by radial effect for pin thread tooth and coupling thread tooth can be, respectively, expressed as
(9)δrpx=Nxcosα-γEtan2α-γ2P·d22x+4r02d22x-4r02-νd2x,(10)δrcx=Nxcosα-γEtan2α-γ2P·4R02+d22x4R02-d22x-νd2x.
Consequently, the total elastic deformations δtp, δtc for pin thread tooth and coupling thread tooth along its cone direction (Z direction) are, respectively, as follows:
(11)δtpx=δa+δb+δc+δd+δrp=k+krpxNxcosα-γE,(12)δtcx=δa+δb+δc+δd+δrc=k+krcxNxcosα-γE.
In (11) and (12), k=ka+kb+kc+kd. The deformation coefficients ka, kb, kc, kd, krp, and krc can be determined from (5) to (10).
In (5)~(12), N(x) denotes the normal load acting on the load flank of thread tooth within unit axial length dx. Because the load flank angle and thread taper are both very small in premium threaded connection, we ignored the axial components of frictional force on load flank and that of the contact stress between root and crest. The relation between N(x) and f(x) can be briefly expressed as
(13)Nxcosαtanψxdx=fxdx.
In (13), the helix angle ψ(x) is
(14)ψx=arctanPπd2x,d2x=d2s-2xtanγ.
By submitting (1) and (14) into (13), we can get
(15)Nx=Pπd2xcosαdFxdx.
Then, by submitting (15) into (11) and (12), the total elastic deformations for pin thread tooth and coupling thread tooth are, respectively, as follows:
(16)δtpx=λpxdFxdx=k+krpxEcosα-γcosαPπd2xdFxdx,(17)δtcx=λcxdFxdx=k+krcxEcosα-γcosαPπd2xdFxdx.
In (16) and (17), the coefficients λp and λc depend on the geometry parameters of connection and thread tooth as well as the used material properties, whose derivatives for x are, respectively, as follows:
(18)λp′=dλpdx=cosα-γEcosαdk+krpP/πd2dx=cosα-γEcosαdk+krpdxPπd2+k+krpdP/πd2dx=cosα-γEcosαd22+4r02d22-4r02-16d22r02d22-4r022-νtan2α-γtanγPhhhhhhhhhhhhh·d22+4r02d22-4r02-16d22r02d22-4r022-νhhhhhhhhhhhhh·Pπd2-k+krp2Ptanγπd22d22+4r02d22-4r02-16d22r02d22-4r022-ν=tanγcosα-γπEd2xcosαd24x-16r04-16r02d22xd22x-4r022-νtan2α-γ2Pd2xhhhhhhhhhhhhhhhhh·d24x-16r04-16r02d22xd22x-4r022-νhhhhhhhhhhhhhhhhh-k+krp2Pd2xd24x-16r04-16r02d22xd22x-4r022-ν,(19)λc′=dλcdx=cosα-γEcosαdk+krcP/πd2dx=cosα-γEcosαdk+krcdxPπd2+k+krcdP/πd2dx=cosα-γEcosα4R02+d224R02-d22+16R02d224R02-d222-νtan2α-γtanγPhhhhhhhhhhhh·4R02+d224R02-d22+16R02d224R02-d222-νhhhhhhhhhhhh·Pπd2-k+krc2Ptanγπd224R02+d224R02-d22+16R02d224R02-d222-ν=tanγcosα-γπEd2xcosα16R04-d24x+16R02d22x4R02-d22x2-νtan2α-γ16R04-d24x+16R02d22x4R02-d22x2hhhhhhhhhhhhhhhhh·16R04-d24x+16R02d22x4R02-d22x2-νhhhhhhhhhhhhhhhhh-16R04-d24x+16R02d22x4R02-d22x2k+krc2Pd2x16R04-d24x+16R02d22x4R02-d22x2-ν.
2.3. Differential Equation of the Axial Load F(x)
Under the action of pretightening force, the pin thread and coupling thread are always engaged to each other, indicating that their thread pitches are equal after deformation. However, the compressed pin will make its own thread pitch decrease while the stretched coupling will make its own thread pitch increase. The deformation of pin thread tooth should increase its own pitch while the deformation of coupling thread tooth should reduce its own pitch, thus guaranteeing the equal pitch and the perfectly engaged state for pin thread and coupling thread.
According to Hooke’s law, the compressed and stretched elastic deformations of pin and coupling body along the thread cone direction are, respectively, as follows:
(20)χp=4πEcosγ∫0xFxd2x-hB2-4r02dx,(21)χc=4πEcosγ∫0xFx4R02-d2x+hB2dx.
Moreover, the deformation compatibility equation can be obtained from the equal pitch after deformation for pin thread and coupling thread:
(22)χpx+χcx=δtpx+δtcx-δtp0+δtc0.
Then, by submitting (16), (17), (20), and (21) into (22), we get
(23)4πEcosγ∫0x1d2x-hB2-4r024πEcosγhhhhh+14R02-d2x+hB2Fxdx=λp+λcdFxdx.
Taking the derivative with respect to x for (23) and considering that λp+λc≠0, the differential equation of axial load distribution can be obtained as
(24)d2Fxdx2+mxdFxdx+nxFx=0,mx=λp′+λc′λp+λc,nx=-4πEλp+λccosγhhhhhh·1d2x-hB2-4r02+14R02-d2x+hB2.
The boundary conditions for (24) are
(25)F0=0,FL=Ft.
2.4. Solution of Differential Equation
The axial load distribution equation (24) is a nonlinear differential equation and its analytical solution is very difficult to get, so we solved its numerical solution with finite difference method. First, we divided the engaged thread according to unit length of thread pitch P and numbered the node 0~Nt successively from the big end of the pin thread, so the coordinate for node i was xi=iP. Then we took Fi to denote the node force F(xi) and considered the boundary conditions equation (25). Finally, the difference equation of (24) can be obtained as
(26)AF1F2⋯FN-1⊤=B.
In (26), the coefficient matrix is(27)A=-2P2+nx11P2+mx12P01P2-mx22P-2P2+nx21P2+mx22P⋯1P2-mxi2P-2P2+nxi1P2+mxi2P⋯⋯01P2-mxN-12P-2P2+nxN-1.
And the constant term column matrix is
(28)B=00⋯-1P2+mxN-12PFt⊤.
The difference equation (26) was solved based on Gaussian main elimination method.
2.5. Geometric Description for Load Flank of Thread Tooth
For calculating the mean contact stress on the load flank of thread tooth, the geometric description of thread surface is necessary. Setting the premium thread connection as the x-axis and its intersection with the plane of thread load flank as the origin, we established global coordinate system O-xyz, whose unit vectors are i, j, and k, respectively, shown in Figure 5.
Geometric description schematic of thread load flank.
For determining the unit normal vector n and unit tangent vector τ of thread surface, the cylindrical coordinate system O-rϑx was also established and its vectors are
(29)u=-j,v=-k,w=i.
The plane OB2C2D2, which stands for thread load flank, can be obtained by rotating the quadrangle OBCD in yoz plane α around the y-axis and then ψ around the z-axis successively. If quadrangle OBCD is a 1 × 1 square, its vertices after 2 rotations are O(0, 0, 0), B2(sinα, 0, cosα), C2(sinα+sinψ, cosψ, cosα), and D2(sinψ, cosψ, 0). Thus, the equation of load flank plane can be obtained as
(30)x-ytanψ-ztanα=0.
The unit normal vector n and unit tangent vector τ of thread load flank are
(31)n=i-tanψj-tanαk1+tan2ψ+tan2α=tanψu+tanαv+w1+tan2ψ+tan2α,τ=-sinψi-cosψj=cosψu-sinψw.
And the infinitesimal area of thread surface dΩ can be also expressed with the infinitesimal area dA in the axial cross section as
(32)dΩ=1+tan2ψ+tan2αdA.
2.6. Mean Contact Stress on Load Flank of Thread Tooth
Supposing that the contact stress on load flank was uniform and that the helix angle was constant in each engaged thread, the axial component of contact stress integral on the load flank within each engaged thread should be equal to the axial load acting on itself; that is,
(33)∬Ωipini+μtpiτidΩ·i=FiP.
By submitting (31) and (32) into (33), the mean contact stress on load flank of ith tooth is
(34)pi=FiP1-μt1+tan2ψi+tan2αsinψi∬AidA.
In (34),
(35)∬AidA=∫02π∫d2i-1P-ϑptanγ/π-hB/2d2i-1P-ϑptanγ/π+hB/2rdrdϑ=πhBd2s-2i+1Ptanγ.
At last, the mean contact stress on the load flank of ith tooth is
(36)pi=FiPπhB1-μt1+tan2ψi+tan2αsinψihhhhhhhh·d2s-2i+1Ptanγ1-μt1+tan2ψi+tan2αsinψi-1.
3. Calculation Example3.1. Analytical Method
For simplicity, the geometrical parameters of a premium connection are mainly taken from API 177.8 mm, P110 grade, buttress thread casing in API Spec.5B. Its geometry parameters and material properties are listed in Table 1.
Geometry parameters and material properties of premium connection.
Symbol
Value
Unit
Symbol
Value
Unit
R0
98
mm
Frs
2.39
mm
r0
78.54
mm
hB
1.58
mm
P
5.08
mm
Nt
15
Integer
d2s
176.63
mm
E
206
GPa
γ
1.79
Degree
ν
0.3
Dimensionless
α
3
Degree
μt
0.08
Dimensionless
β
10
Degree
fymn
758.42
MPa
Suppose the pretightening force Ft from sealing surface and shoulder contact is 1000 KN; the axial load acting on each engaged thread FiP and the mean contact stress on load flank pi can be calculated according to the analytical model proposed in this paper; the results are shown in Figures 6 and 7, respectively.
Load acting on each engaged thread by the proposed analytical model.
Mean contact stress on load flank of thread tooth in premium connection obtained by analytical method and finite element method.
Figure 6 shows that the load acting on engaged thread tooth increases with the increase of thread number and that the last 4~6 engaged threads bear more than 90% of the pretightening force. Particularly for the last engaged thread, it carries 44.54% of the pretightening force. The load distribution is similar to the practical situation in threaded connection, indicating that the proposed analytical method is reasonable. Meanwhile, contact stress distribution on load flank of thread teeth shows a similar pattern to the load distribution, whose maximum value is 522.74 MPa on the last engaged thread, shown in Figure 7. It is obvious that thread surface would sustain a galling risk if the pretightening force is increased to some extent.
3.2. Finite Element Method
In order to validate the proposed analytical model in this paper, we analyzed the same sample with the finite element software ANSYS 14.5. Considering the characteristics of axisymmetric load and boundary conditions, a 2D finite element model of modified API 177.8 mm premium thread connection has been established (shown in Figure 8).
Finite element model of modified API 177.8 mm premium connection.
To make the simulated results more approachable to practical situations, the materials of pin and coupling were both defined as bilinear isotropic hardening model in the FEM. The frictional contact type between pin thread and coupling thread has been selected, indicating that the contact between thread teeth is taken as an elastic-plastic and frictional contact problem. The eight-node quadrilateral element Quad 8 node-183, which can be used as an axisymmetric element to simulate plastic deformation, was selected for this analysis. The element size was controlled and then local mesh refinement was applied to the thread tooth edges. Figures 9 and 10 show the mesh on the whole model and the teeth edges. Because the contact between engaged thread teeth is a material and boundary nonlinear problem, static analysis has been selected for solving it. The boundary conditions set for this analysis are an axial displacement confined at the big end of coupling thread and application of 313.93 MPa compressed stress, which is equivalent to 1000 KN of pretightening force, on the small end of pin thread. By using ANSYS, we can acquire the stress distribution on the thread teeth (shown in Figure 11).
Mesh on the whole finite element model.
Refined mesh at teeth edges: (a) left region and (b) right region.
Von-Mises stress distribution on the thread teeth in ANSYS.
Figure 11 shows the Von-Mises stress distribution on each contact thread tooth. In Figure 11, we can find that the last 4~6 engaged threads bear almost all the pretightening force, which is identical to the load distribution calculated by the proposed analytical model in this paper.
3.3. Comparison and Discussion
The mean contact stress on the load flank of thread teeth calculated from the proposed analytical model and that obtained by FEM are compared in Figure 7. One can find that both curves are identical in trend and that the minimum and maximum values obtained by the proposed analytical method are also close to those of FEM. Meanwhile, we can also find that there exists some discrepancy in the middle region of the two curves, which may be attributed to the simplifications and assumptions of the analytical model. The helix angle in the analytical model may be a reason because the 2D FEM cannot contain this factor. However, the analytical model can still predict a relatively meaningful stress value for the last engaged threads, which is the easiest one to fail.
Figure 12 shows the comparison of mean contact stresses on the last engaged threads under different pretightening forces. It can be seen from Figure 12 that mean contact stress on the last engaged thread obtained by analytical model keeps increasing linearly with the increase of pretightening force, while that calculated from FEM tends to be stable after a section of linear increase. This discrepancy can mainly because that elastic deformation has been supposed for the analytical model but elastic-plastic deformation considered in FEM. Because the allowable in-service pretightening force is about 1350 KN~1750 KN, the error of the analytical model for the stress on the last engaged thread is about 6.2%~18.1%. Consequently, the proposed analytical method has certain accuracy and can, to some extent, meet the needs of engineering.
Comparison of mean contact stresses on the last engaged thread tooth under different pretightening forces.
4. Conclusions
This paper proposes an analytical model for investigating loading and contact stress on the thread teeth in tubing and casing premium threaded connection on the basis of elastic mechanics.
The proposed analytical model is validated by the application of the FEM to the same sample. The pretightening force mainly applies to the last 4~6 engaged threads, about 50% of that to the last engaged threads.
Comparing the contact stress of the last engaged thread between the analytical model and FEM shows that the accuracy of the analytical model will decline with the increase of pretightening force after the material enters into plastic deformation. In practice, the relative error is about 6.2%~18.1%. This indicates that the analytical method can, to some extent, meet the needs of engineering.
NomenclatureFt:
Total force transferred by engaged threads
Nt:
Total number of engaged thread teeth
N(x):
Normal load acting on the load flank of thread tooth within unit axial length dx
f(x):
Axial load distribution intensity
F(x):
Axial load distribution
FiP:
Axial load acted on each engaged thread
P:
Pitch of thread
α:
Load flank angle in the axial cross section of thread
β:
Leading flank angle in the axial cross section of thread
ψ(x):
Helix angle
d2(x):
Pitch diameter
d2s:
Pitch diameter at the big end of pin thread
γ:
Half angle of thread cone
a:
Root width of isosceles trapezoid tooth
b:
The width of isosceles trapezoid tooth at pitch diameter
c:
The distance from pitch diameter of isosceles trapezoid tooth to its root
Frs:
Root width of buttress type tooth
hB:
Completed thread height
δtp:
Total elastic compressed deformation of pin thread tooth
δtc:
Total elastic stretched deformation of coupling thread tooth
δa, δb, δc, δd:
Elastic deformation of thread tooth caused by tooth bending, shear force, root incline, and shear deformation of tooth root, respectively
δrp, δrc:
Elastic deformation of thread tooth caused by radial effect for pin thread tooth and coupling thread tooth, respectively
ka, kb, kc, kd, krp, krc, k, λp, λc:
Relevant deformation coefficients of thread tooth
E:
Elastic modulus of material
ν:
Poisson ratio of material
fymn:
Specified minimum yield strength of material
r0:
Internal radius of pipe
R0:
External radius of coupling
χp, χc:
Compressed and stretched elastic deformation along thread cone direction for pipe and coupling body, respectively
m,n:
Coefficients of differential equation
Fi:
Axial load on ith node
pi:
Mean contact stress on the load flank of ith tooth
μt:
Frictional coefficient on the thread surface
ψi:
Mean helix angle on ith tooth.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors are grateful for the support of the National Natural Science Foundation of China (Grant nos. 50904050 and 51244007) and the financial support of Tianjin Pipe (Group) Co., Ltd.
CarcagnoG.The design of tubing and casing premium connections for HTHP wellsProceedings of the SPE High Pressure/High Temperature Sour Well Design Applied Technology Workshop2005Society of Petroleum EngineersBradleyA.NagasakuS.VergerE.Premium connection design testing and installation for HPHT sour wellsSPE High Pressure/High Temperature Sour Well Design Applied Technology Workshop2005Society of Petroleum EngineersSuginoM.NakamuraK.YamaguchiS.DalyD.BriquetG.VergerE.Development of an innovative high-performance premium threaded connection for OCTGProceedings of the Offshore Technology Conference (OTC '10)May 2010188318912-s2.0-79955795705MaduschkaL.Beanspruchung von Schraubenverbindungen und zweckmäßige Gestaltung der Gewindeträger19367629930510.1007/BF025849122-s2.0-51649213923SopwithD.The distribution of load in screw threads19481591373383YazawaS.HongoK.Distribution of load in the screw thread of a bolt-nut connection subjected to tangential forces and bending moments1988312174180WangW.MarshekK. M.Determination of the load distribution in a threaded connector having dissimilar materials and varying thread stiffness199511711810.1115/1.28032732-s2.0-0029251567BruschelliL.LatorrataV.The influence of “shell behavior” on load distribution for thin-walled conical joints200067229830610.1115/1.13048412-s2.0-0034205949ChenS.-J.AnQ.ZhangY.GaoL.-X.LiQ.Loading analysis on the thread teeth in cylindrical pipe thread connection201013230312021031202810.1115/1.40007292-s2.0-77955181320O'HaraP.Finite element analysis of threaded connectionsProceedings of the Army Symposium on Solid Mechanics197499119DTIC DocumentBretlJ. L.CookR. D.Modelling the load transfer in threaded connections by the finite element method19791491359137710.1002/nme.16201409092-s2.0-0018730518ZhaoH.A numerical method for load distribution in threaded connections1996118227427910.1115/1.28268802-s2.0-0030165973ZhaoH.Stress concentration factors within bolt-nut connectors under elasto-plastic deformation199820965165910.1016/S0142-1123(98)00036-X2-s2.0-0032187553MacdonaldK. A.DeansW. F.Stress analysis of drillstring threaded connections using the finite element method19952113010.1016/1350-6307(95)00007-D2-s2.0-0029257045BahaiH.GlinkaG.EsatI. I.Numerical and experimental evaluation of SIF for threaded connectors199654683584510.1016/0013-7944(95)00251-02-s2.0-0030214641ChenJ.-J.ShihY.-S.Study of the helical effect on the thread connection by three dimensional finite element analysis1999191210911610.1016/S0029-5493(99)00134-X2-s2.0-0033155564BaragettiS.Effects of taper variation on conical threaded connections load distribution2002124232032910.1115/1.14564592-s2.0-0036623945FukuokaT.NomuraM.Proposition of helical thread modeling with accurate geometry and finite element analysis200813011601120410.1115/1.28264332-s2.0-47049084707StoecklyE.MackeH.Effect of taper on screw thread load distribution1952741103112YuanG.YaoZ.HanJ.WangQ.Stress distribution of oil tubing thread connection during make and break process200411453754510.1016/j.engfailanal.2003.10.0072-s2.0-2142753045PlácidoJ. C. R.Valadã de MirandaP. E.NettoT. A.PasqualinoI. P.MiscowG. F.de Carvalho PinheiroB.Fatigue analysis of aluminum drill pipes2005844094152-s2.0-33644537917SlackM.SalkinH.LangerF.Technique to assess directly make-up contact stress inside tubular connectionsSPE/IADC Drilling Conference1990Society of Petroleum EngineersHamiltonK. A.WaggB.RothT.Using ultrasonic techniques to accurately examine seal surface contact stress in premium connectionsProceedings of the SPE Annual Technical Conference and ExhibitionNovember 2007Society of Petroleum Engineers346134712-s2.0-40349085931YamamotoA.1980Tokyo, JapanYoukendo(Japanese)