Finite-Element Analysis on Energy Dissipation and Sealability of Premium Connections under Dynamic Loads

: In the process of high ﬂ ow rate fracture and high gas production, the sealing performance of the premium connection decreases due to the dynamic load and vibration of downhole tubing strings, which may cause accidents. Existing static analysis methods cannot e ﬀ ectively explain this phenomenon. The main objective of this paper is to propose a novel analytical method for evaluating the sealing performance of a premium connection. In this paper, a dynamic model of sealing surfaces of the premium connection is established based on the vibration equation of elastic rod, and the hysteresis characteristics and energy dissipation mechanism of sealing surfaces are analyzed. Considering the in ﬂ uence of spherical radius, internal pressure, axial cyclic load amplitude, and modal vibration, a spherical-conical premium connection ﬁ nite element model is established to analyze the in ﬂ uence laws of the connection’s energy dissipation and sealing performance. The results show that the sealing performance of the premium connection under dynamic load can be e ﬀ ectively analyzed by using energy dissipation theory compared with traditional static contact analysis. Compared with the vibration of the tubing string, the dynamic loads caused by the change of ﬂ uid pressure and ﬂ ow rate in the tubing string have a signi ﬁ cant in ﬂ uence on the connection’s sealing performance. When the internal pressure and axial cyclic loads are 80 MPa, 400 kN, or 60 MPa and 500 kN respectively, serious plastic deformation occurs in the thread and sealing surfaces, and the energy dissipation of the sealing surfaces increases signi ﬁ cantly, which could lead to sealing failure.


Introduction
The premium connection serves as the connection structure for downhole tubing, connecting tubing to thousands of meters of the tubing string. As shown in Figure 1, the premium connection consists mainly of sealing surfaces, torque shoulders, and threads. The sealing surface is the main sealing structure, the torque shoulder has the function of resisting torque overload and auxiliary sealing, and the threads are used to connect the tubing. Connections are the most dangerous part of the tubing string. Due to complex downhole loads, connections leak frequently and even lead to accidents [1,2]. According to field statistics, connection failure accounts for 85-95% [3,4] of all tubing string failures. Earlier research on the premium connection was mainly based on static load, and the sealing performance was studied by analyzing the equivalent stress, contact pressure, and contact length of the connection. Li et al. [5] used a combination of experimental and theoretical research to study the friction characteristics of bolted connection interfaces. Xu et al. [6][7][8] used a combination of finite element simulation and theoretical models to calculate the contact pressure and contact area of the sealing surfaces, taking into account the micro profile of the sealing surfaces. The relationship between the average contact pressure, circumferential leakage width, and sealing surfaces was obtained. Nasraoui et al. [9] constructed a simplified numerical model of bolt bonding surfaces and established a contact stiffness matrix. Xu et al. [10] established contact pressure distribution and sealing surfaces length of premium connection under different buckling torques by analytical method and further evaluated the sealing condition of the connection by a comprehensive consideration of various factors. Zhiqian et al. [11] analyzed the sealing performance of the connection by combining the elastic-plastic mechanical contact theory and calculated geometric morphology parameters of the rough surfaces. Chen et al. [12,13] analyzed the effect of preload on the threaded seal and the influence of contact stress at threads, sealing surfaces and shoulders on the sealing performance. Affected by loads such as the weight of tubing string, internal and external fluid pressure, viscous damping when fluid flows in tubing string and Coulomb friction force, premium connections are subjected to loads such as axial force, internal and external pressure, bending moment, etc. During high flow rate fracturing and high production gas operation, the dynamic load change of tubing caused by the change of fluid pressure and flow rate in tubing causes vibration of tubing [14][15][16]. Under dynamic load, the sealing surfaces will slip cyclically, resulting in connection loosening or poor sealing performance. However, the traditional static contact analysis method cannot effectively explain this problem. As shown in Figure 1, the sealing surfaces can be considered as two contact surfaces, which are mainly subjected to normal and tangential loads. Normal loads are generated by the interference fit of sealing surfaces and internal and external pressure, which affect contact morphology and pressure of sealing surfaces. Tangential loads are generated by the axial force of the tubing strings and will cause the horizontal slip of the sealing surfaces. Under cyclic loading, the sealing surfaces of premium connections will undergo microslip or gross slip, which will lead to energy dissipation in the structural system [17,18]. As shown in Figure 2, under the action of external incremental load, microslip (stick and slip coexist) occurs on the sealing surfaces, and the slope of the forcedisplacement curve decreases gradually, resulting in stiffness degradation. As the load increases further, gross slip occurs on the sealing surfaces. At this time, the force-displacement curve approaches the level and the force remains constant, and the gross slip increases gradually [19]. As shown in Figure 3, under cyclic loading, the microslip and gross slip are converted into each other during loading and unloading, resulting in force-displacement hysteresis [17,20,21]. The area of the force-displacement hysteresis curve is the energy dissipation generated during a cycle of loading [22]. Loosening of the connection can be explained by the phenomenon of force-displacement hysteresis of contact surfaces. In the early stage, Iwan [23] built an energy dissipation model of bolted connections based on spring-damped vibration theories. Wang et al. [24] optimized the Iwan model to apply it to the study of energy dissipation in more complex situations. Zhang et al. [25,26] studied the normal stiffness, normal damping, and energy dissipation of the bolt bonding surfaces under different external loads.  The premium connection has more complex geometry and load conditions than the bolted joint. At present, there are fewer reports about the sealing property of premium connections under dynamic load. Yang et al. [10] established a contact stress model of sealing surfaces and analyzed the energy dissipation law under different sealing conicities. Yu et al. [27] established a shear layer microslip model of the sealing surfaces at an early stage and analyzed the influence of the shear layer coefficient on the microslip state and energy dissipation of the sealing surfaces. Yu et al. [28]. studied the sealing performance of special threaded joints from a micro level, and based on fractal contact theory, analyzed the effects of fractal dimension and amplitude coefficient on the contact behavior of rough surfaces. The above researches provide some new ideas for the sealing analysis of premium connection, but the results are limited by simplifying models and assuming conditions, such as ignoring the threads and torque shoulders.
In summary, to address the limitations of existing sealing performance evaluation methods, this paper considers the influence of threads and torque shoulders on the contact pressure of the sealing surface of premium connections. On the basis of the traditional evaluation of sealing performance, energy dissipation theory is added to evaluate the sealing performance of premium connections under dynamic loads. We analyze the impact of different parameters on sealing performance under different working conditions based on energy dissipation values. Firstly, based on the vibration equation of the elastic rod, a dynamic model of the connection sealing surfaces is established to analyze the energy dissipation mechanism. Then, based on friction contact analysis and vibration mode analysis, a finite element model (FEM) of spherical-conical premium connection is established by ABQUS numerical simulation software. Finally, the influence of different spherical radii, internal pressure, axial cyclic load amplitudes, and modal shapes on energy dissipation and sealing performance of premium connection under dynamic load is analyzed. The sealing performance evaluation method established based on the energy dissipation theory can provide technical guidance for the design of spherical-conical premium connections and serve as a reference basis for selecting suitable premium connections on oilfield construction sites.

Materials and Methods
The energy dissipation of the sealing surfaces is reflected by the force-displacement hysteresis curve. In order to obtain the force-displacement hysteresis curve, a non-linear vibration model of the sealing surfaces must be established first. As shown in Figure 4, the two sealing surfaces of the Pin and the Box are simplified to an elastic rod and a rigid body, respectively, and the contact slip behavior between them is represented by a shear layer without thickness.

Dynamic Model Establishment
Based on the vibration balance equation of the rod, the governing balance equations for the stick and slip regions can be obtained as [18,27]: Stick regions: Microslip regions: Boundary conditions: The continuity conditions satisfied at the stick and slip boundary positions, x = ln, is: In equation: E-is the Elastic modulus of elastic rod, GPa; A-is the section area of the elastic rod, mm 2 ; θ-is the taper of the sealing surface, °; k-is the shear layer stiffness in the stick region, MPa; L-is the horizontal length of the shear layer, mm; ln-is the length of the stick zone, mm; u-is the displacement of a point on the shear layer, mm; p(x)-is the non-uniform normal pressure distribution, N/mm; F-is the tangential load on the elastic rod, kN; µ-is the coefficient of friction The normal pressure decreases linearly in the direction of the elastic rod as follows: In formula: P0-is the maximum normal pressure, N/mm; kP-is the slope of normal pressure.
From Equations (1) and (2), the boundary condition Equation (3), and the continuity condition (4), the displacements in the stick and slipping regions of the elastic rod are obtained as follows [29]: F that can be obtained by simultaneous Equations (6) and (7): In Equation (8), when ln = L, the minimum load at which the microslip occurs can be obtained: When ln = 0, the maximum tangential force for gross slip can be obtained: Based on Masing's hypothesis of steady-state cyclic hysteresis response, the forcedisplacement curves during unloading and reloading can be obtained [30]. Masing's hypothesis holds that the unloading and reloading of the system's steady-state hysteresis response are geometrically similar to the initial force curve, but only amplified twice.
Elastic rod unloading process [31]: In the formula, F0 and FA are the starting force at unloading and the force during unloading, N; u(L)0 and u(L)A are the initial displacement at unloading and the displacement during unloading, mm.
Reloading process [31]: (7) and (8), the force-displacement hysteresis curve of the sealing surfaces under cyclic load can be obtained. The area included in the hysteresis curve is the energy dissipation generated by each cycle of cyclic loading [18].
where, h E  is the microslip energy dissipation of per cycle, J; un(x) is the critical displacement when the stick region is converted into slip region, mm.

Verification of Dynamic Model
To verify the above model, referring to the sealing surfaces' structure and constraints shown in Figure 4, a FEM of friction contact of the sealing surfaces as shown in Figure 5 is established by ABAQUS. The FEM selects the C3D8I element suitable for the contact analysis. Mesh refinement is applied to particular areas, such as sealing surfaces and shoulder regions, with a mesh size of 0.10. Areas not deemed critical are assigned a relatively rough mesh in order to reduce computational time. Abaqus/Standard Solver is used to analyze finite element models. The interaction attribute is set using the Coulomb friction model to set the contact attribute, the tangential behavior friction formula uses a penalty function, and the normal behavior is set to hard contact [10,32]. Set the Box (rigid body) with higher hardness as the main contact surface and set the Pin (elastic body) as the secondary contact surface. Select a penalty friction formula to simulate a small amount of relative movement during the sticking of the contact surfaces. The left end of the rigid body is set to be fully constrained and a cyclic tangential load is applied at the right end of the elastic body, and a normal pressure is applied at the contact surface. According to the dimensions and material parameters of premium connections and literature [18], the parameters in Table 1 are respectively substituted into the formula (11) and (12) and the FEM and the force-displacement hysteresis curves of the two models are calculated (as shown in Figure 6). From the figure, it can be seen that the force-displacement hysteresis curves of both models approximate olive-shaped curves, and the degree of overlap between the two curves is relatively high. Under monotonic loading, the force-displacement curves consist of linear and nonlinear segments. As the load increases, initially, the displacement of both models shows linear increments with nearly the same magnitude. Then, with further increases in load, the displacement increases nonlinearly, with the FEM exhibiting slightly greater increments than the shear layer model, indicating that the FEM dissipates more energy than the shear layer model. The energy dissipation values of the two models are 0.0399 J and 0.0414 J, respectively, with a difference of 3.5%, thus validating the effectiveness of the FEM.

FEM of Premium Connection
From the analysis in the previous section, it is known that the energy dissipation of the sealing surfaces is characterized by a force-displacement hysteresis curve, and the effectiveness of using numerical simulation methods to analyze the energy dissipation of the sealing surfaces is verified. In order to study the influence of sealing structure and downhole loads on energy dissipation and sealing performance of premium connections, a friction contact FEM is established in Section 2-a φ 88.9 × 6.45 P110 full-size sphericalconical premium connection. As shown in Figure 7, the sealing surface of the coupling part is a 1:16 conical surface, the sealing surface of the tubing end is a spherical surface with a radius of 10 mm, the torque shoulder is 25°, the connection thread bearing angle is −3°, the guide surface is 10° and the thread taper is 1:16.

FEM for Friction Contact
Due to the elastic and plastic deformation of the sealing surfaces and thread under the action of make-up and external loads, it is necessary to enable the geometric nonlinearity and material nonlinearity settings in ABAQUS software. The relevant material parameters for premium connections are shown in Table 2. The C3D8I element is selected for premium connection, and mesh refinement is performed on the threads, sealing surfaces, and torque shoulders, as shown in Figure 8. To accelerate the solution speed of nonlinear analysis and ensure the accuracy of calculation, ABAQUS/Standard solver is selected for solution analysis [33,34]. After make-up, the premium connection's contact includes the threads, sealing surfaces, and torque shoulder. In nonlinear contact analysis, surface-to-surface contact elements are selected and contact pairs are established according to the method in Section 3. As shown in Figure 9, the coupling's face is set with fixed full constraints, a reference point is arranged at the axis position of the tubing and coupled with the face of the tubing to apply an axial cyclic load to the reference point, and a constant internal pressure is applied to the inner wall of the tubing and the coupling.

FEM for Vibration
A FEM of tubing with external threads (each tubing is 9.8 m in length) is established and fully constrained at both ends of the tubing to form a vibration model of the tubing string containing the connection. The material parameters, meshing, and contact settings of the FEM refer to Section 2.1. The finite element meshed model is shown in Figure 10.

Mesh Independent Verification
In order to minimize the errors in the calculation results caused by the quality and quantity of the mesh during the analysis process, the paper adopts a mesh-independent verification method. At the sealing surface of the premium connection, this article uses C3D8I elements for mesh division, with mesh sizes of 0.05, 0.1, 0.15, 0.2, and 0.25 mm, respectively [35], the maximum Von Mises stress distribution at the sealing surface of the premium connections under the optimal tightening torque is shown in Figure 11. When the mesh size decreases from 0.25 mm to 0.15 mm, the magnitude of the maximum stress increases monotonically from 0 to 2.12%. When the mesh size is between 0.05 and 0.15 mm, the maximum stress change amplitude at the sealing surface is between 4.0 and 4.4%. The results show that when the mesh size is less than 0.15 mm, the change in mesh size has a small impact on the maximum contact stress value at the sealing surface. As in finite element simulation, the smaller the mesh size, the more accurate the convergence of the obtained results. Taking into account the calculation time and efficiency, this paper adopts a mesh size of 0.10 for analysis.

Result
Structural parameters and load conditions are the main factors affecting the sealing performance of premium connections. For spherical-conical premium connections, the radius determines the initial contact conditions such as the contact area, length, and contact pressure of the sealing surfaces. In addition, the dynamic load and vibration of the downhole string will change the contact state of the sealing surfaces, threatening the sealing performance of the connection. This section considers the influence of the radius of the sealing surfaces, the pressure inside the tubing, the axial force amplitude on the tubing string (As shown in Table 3), and the modal vibration mode under the axial cyclic load as shown in Figure 12. The finite element method is used to analyze the von Mises stress, contact pressure, and force-displacement hysteresis curve, obtaining the energy dissipation and sealing performance of the premium connection under various influencing factors.   Figure 13 shows the von Mises stress nephogram of premium connection with different spherical radii under 60 MPa internal pressure and 400 kN axial cyclic load. As can be seen from Figure 13, under axial compression load, there is significant stress concentration at the thread portion of the connection. When the spherical radius is larger than 15 mm, the maximum von Mises stress of the connection is 834.8 MPa, which is near the strength limit of the material (835 MPa). However, under axial tension load, the overall von Mises stress distribution of the connection is uniform and does not exceed the maximum yield limit of the material (828 MPa). This is because there is a preload at the thread in the axial direction during the make-up of the premium connection. Due to the same axial compression load and preload direction, the von Mises stress in the connection increases, and the stress concentration at the thread is significant.

Axial compression
Axial tension (a) The spherical radius is 10 mm

Axial compression
Axial tension (b) The spherical radius is 15 mm

Axial compression
Axial tension (c) The spherical radius is 20 mm  Figure 14 shows the contact pressure curves of sealing surfaces with different spherical radii during axial cyclic loading. In the figure, 0-1 s is the internal pressure loading process, and 1 s-2 s is the axial cyclic load loading process. The wave peak is the maximum axial compression load, and the wave trough is the maximum axial tension load. As can be seen from Figure 14, during initial internal pressure loading, the larger the spherical radius, the lower the contact pressure on the sealing surfaces. During the entire axial load cycle, the contact pressure hardly changes with the increase of the spherical radius, and only decreases by 1.8% during axial tension. Extracting the friction force and node displacement at each analysis step on the sealing surfaces within a cycle, the force-displacement hysteresis curves with different spherical radius is obtained as shown in Figure 15, and the gross slip and energy dissipation curves as shown in Figure 16. From Figure 15, it can be seen that the force-displacement hysteresis curve of the spherical-conical seal structure approximates a parallelogram, and its gross slip is much greater than the microslip, indicating that the energy dissipation of the seal structure is mainly caused by gross slip [36].
As can be seen from Figure 16, energy dissipation increases as the spherical radius increases. When the spherical radius is increased from 10 mm to 20 mm, the gross slip increases by 5.2% and 5.4%, and the energy dissipation increases by 2.9% and 3.8%, respectively (as shown in Table 4). This is because the larger the spherical radius of the surface, the lower the contact pressure, resulting in a decrease in the friction required for gross slip of the sealing surfaces, and an increase in the amount of gross slip at the connection sealing surfaces. At the same time, since gross slip is the main factor causing energy dissipation, the larger the spherical radius, the greater the energy dissipation value. This will lead to increased wear on the sealing surfaces, which will reduce the sealing performance of the connection.
From the above analysis, it can be seen that under the same axial cyclic load, it is difficult to analyze the impact of the spherical radius on the sealing performance of the connection by the contact pressure of the sealing surfaces; however, the energy dissipation theory can effectively solve this problem.    Figure 17 shows the interaction of 400 kN axial cyclic load with 60 MPa, 70 MPa, and 80 MPa internal pressure. As can be seen from Figure 17, with the increase of internal pressure, the von Mises stress in the connection increases. When the internal pressure is 80 MPa, the maximum von Mises stress of the connection under axial compression load is 844.4 MPa, which exceeds the strength limit of the material. The stress concentration of the thread is a potential danger point.   Figure 18 shows the contact pressure curve of sealing surfaces under the combined action of axial cyclic load and different internal pressures. Table 5 shows the maximum contact pressure of sealing surfaces under different internal pressures. When the internal pressure increases from 60 MPa to 80 MPa, the contact pressure decreases by 0.5% and 4.6%. During the axial compression stage, the contact pressure decreased by 1.5% and 12% respectively (as shown in Table 5). Existing research has shown that the greater the internal pressure [37], the greater the contact pressure on the sealing surfaces. This is because previous studies were mostly based on statics analysis, ignoring the impact of dynamic loads [38]. From Figure 16, it can be seen that serious plastic deformation has occurred between the thread and the sealing surfaces when the internal pressure is 80 MPa. It can be predicted that under the long-term axial cyclic load, the connection will become loose due to the gradual accumulation of plastic deformation [39], which will lead to a decrease in the contact pressure of the sealing surfaces.  Figure 19 is the force-displacement hysteresis curve of the sealing surfaces under the combined action of 400 kN axial cyclic load and different internal pressures. Figure 20 is the corresponding gross slip and energy dissipation curve. It can be seen from Figure 18 that with the increase of internal pressure, the gross slip of sealing surfaces increases, and the energy dissipation increases. The internal pressure increases from 60 MPa to 80 MPa, the gross slip increases by 2.9% and 13.9%, and the energy dissipation increases by 7.8% and 18.2%, respectively (as shown in Table 6). Therefore, when 400 kN axial cyclic load is combined with 80 MPa internal pressure, the energy dissipation of the connection increases significantly, which easily leads to connection loosening and the risk of seal failure exists. In addition, the contact pressure decreases with the loosening of the connection, which explains the decrease of the contact pressure with the increase of internal pressure in the previous section.    Figure 21 shows the von Mises stress nephogram of the premium connection under 60 MPa internal pressure combined with 400 kN, 500 kN, and 600 kN axial cyclic load, respectively. It can be seen from Figure 18 that the von Mises stress of connection increases with the increase of axial compressive load amplitude. However, the magnitude of the axial tension load has little effect on the von Mises stress of the connection. When the axial cyclic load is 500 kN, the maximum von Mises stresses of the connection under the axial compressive load are 864.1 MPa, which exceeds the strength limit of the material. Therefore, when the axial cyclic load is greater than 500 kN, the thread stress concentration is a potential risk point.

Axial compression
Axial tension (a) Axial cyclic load 400 kN

Axial compression
Axial tension (b) Axial cyclic load 500 kN

Axial compression
Axial tension (c) Axial cyclic load 600 kN  Figure 22 shows the contact pressure curve of the sealing surfaces under the combined action of internal pressure of 60 MPa and different axial cyclic loads. Table 6 shows the maximum contact pressure of the sealing surfaces under different axial cyclic loads. It can be seen from Figure 22 that the contact pressure on the sealing surfaces decreases with the increase of the axial tension load amplitude. Compression load has little effect on the contact pressure of sealing surfaces. Table 6 shows that the amplitude of the axial tension load increases from 400 kN to 600 kN and the contact pressure decreases by 6% and 6.5%, respectively (As shown in Table 7). This is because, under the action of axial tension, the threaded part of the connection deforms, which causes the displacement of the sealing surfaces and the reduction of interference between the spherical and conical surfaces, finally resulting in the reduction of contact pressure.   Figure 23 shows the force-displacement hysteresis curve of the sealing surfaces under the combined action of internal pressure of 60 MPa and different axial cyclic loads, and Figure 24 shows the corresponding gross slip and energy dissipation curves. It can be seen from Figure 23 that with the increase of the amplitude of the axial cyclic load, the gross slip of the sealing surfaces increases, and the energy dissipation value increases. The axial tension load increases from 400 kN to 600 kN, the gross slip increases by 39.5% and 27.6%, and the energy dissipation increases by 50% and 25%, respectively (as shown in Table 8). Therefore, when the axial cyclic load is 500 kN, the energy dissipation of the connection increases significantly, which is easy to cause connection sealing failure.
Literature [34] assessed the load of sealing failure by using contact pressure. The result shows that there is a risk of sealing failure (up to 95% of the sealing limit) when the load is 60 MPa internal pressure combined with a 600 kN tension load. Compared with the results in this section, the axial load causing sealing failure is larger. This is because the influence of gross slip of sealing surfaces is not considered in the literature [40].

The Influence of Modal Shapes on Energy Dissipation
First, the modal analysis of the tubing string with premium connection is carried out, and then the modal analysis results are applied as displacement fields to the FEM of the premium connection to obtain the contact pressure and force-displacement hysteresis curves of the sealing surfaces in two vibration cycles. Figure 25 shows the first-order and second-order modes of the premium connection. The modal data of the premium connection are extracted and loaded on the connection in the form of cyclic displacement by editing the ABAQUS input file. Cyclic displacement has two cycles and is divided into 100 incremental steps.
(a) First-order mode (b) Second-order mode  Figure 26 shows the contact pressure nephogram at the sealing surfaces of the connection during two displacement cycles, in which every 0.5 s is a cycle, between 0 s and 0.25 s is the axial compression, and between 0.25 s and 0.5 s is the axial tension. As can be seen from Figure 26, the maximum contact pressure under the first-order mode decreased by 8.7 MPa. The maximum contact pressure under the second-order mode decreased by 8.4 MPa. Due to the fact that the second-order modal shape is smaller than the first-order modal shape, the contact pressure at the sealing surfaces decreases with the increase of the vibration mode order. Figure 27 shows the average contact pressure curve of the sealing surfaces. It can be seen from the figure that the average contact pressure gradually decreases with the increase of the cycle, and the variation amplitude of the average contact pressure under the first-order mode is greater than that under the second-order mode. This phenomenon is also caused by the fact that the second-order modal shape is smaller than the first-order modal shape.
(a) First-order modal shape (b) Second-order modal shape  Extracting the friction force and node displacement of each analysis step within a cycle to obtain the force-displacement hysteresis curve shown in Figure 28. As can be seen from Figure 29, under the vibration modal shape, the force-displacement curve at the sealing surfaces of the spherical-conical premium connection presents a needle-like shape, and there is almost no gross slip of the sealing surfaces. According to Mindlin's [30] theory, the contact edge of the sealing surfaces undergoes microslip and the center of the sealing surfaces is still in the stick state. Due to the fact that the amplitude of the contact pressure change on the sealing surfaces under the first-order mode is greater than that under the second-order mode (Figure 27), and according to Formula (13), it is known that the normal pressure p(x) is positively correlated with the energy dissipation, resulting in the energy dissipation under the first-order mode being greater than that under the second-order mode.
Compared with the force-displacement curve of the connection sealing surfaces under axial cyclic load (Figures 15, 19, and 23), only stick-slip occurs on the sealing surfaces under the first-order and the second-order modes and the energy dissipation value is min-imal (the difference between these is 10 2 orders of magnitude). This indicates that the dynamic load caused by changes in fluid pressure and flow rate in the tubing has a greater impact on the sealing performance of the connection than the impact of string vibration.

Discussion
Through the aforementioned research, it can be found that the method of analyzing the contact pressure of the sealing surface of premium connections under dynamic loading does not effectively explain the decline of connection sealing performance. Firstly, the change in spherical radius under dynamic loading has an insignificant impact on the contact pressure of the sealing surface. Secondly, although internal pressure and contact pressure are positively correlated, the conclusion drawn from traditional analysis methods is that higher internal pressure leads to better sealing performance, contrary to practical situations. It is obvious that all materials have limited strength, exceeding which failure occurs and leads to seal degradation under excessive loading. The method of analyzing the sealing performance of connections using energy dissipation in this article can explain this issue. For spherical-conical premium connections, the gross slip of the sealing surfaces under dynamic load and the microslip of the sealing surfaces under modal vibration mode are analyzed. Both types of slip mentioned above will generate energy dissipation, leading to a decrease in sealing performance.
However, there are still numerous unresolved issues regarding the energy dissipation and sealing performance of premium connections, in particular, the lack of experimental verification. The sealing surfaces of premium connections are typically located in a completely enclosed structure, making it difficult to directly measure the contact pressure and sliding length of the sealing surfaces. This makes it difficult to obtain the forcedisplacement hysteresis curve of the sealing surfaces. When the tubing string is working downhole, there is liquid or gas in the tubing string, and there is repeated momentum and energy exchange at the interface between the fluid and the tubing string, resulting in fluid-structure interaction excitation, which increases the difficulty of the modal analysis of connections. In addition, the relative slip between the sealing surfaces not only generates energy dissipation but also causes wear, further impacting the sealing performance of the connections. In future studies, the accuracy of FEM analysis will be further validated through force-displacement experiments on connections.
The main purpose of this paper is to propose a novel analytical method for evaluating the sealing performance of a premium connection. This method takes into account the dynamic loads and modal shapes and utilizes finite element analysis to obtain the forcedisplacement hysteresis curve of the sealing surface. Through analyzing the slip status and energy dissipation of the connection, the impact of dynamic loads on the sealing performance is revealed. The accuracy of the FEM is preliminarily validated by the vibration equation of the rod. The research results can guide the design, manufacturing, inspection, quality control, and production parameter optimization of the premium connections to ensure their sealing performance. In addition, the results can be applied to the production of oil with a high wax content and production processes in difficult conditions [41]. In future research, finite element simulation will be combined with physical experiments to conduct a more in-depth analysis.

Conclusions
Based on the dynamic model of sealing surfaces, a spherical-conical premium connection FEM is established to analyze the influence laws of connection energy dissipation and sealing performance under different spherical radii, internal pressure, axial cyclic load amplitude, and mode vibration. The results show that the energy dissipation theory can effectively explain the sealing failure of premium connections under dynamic load. The main conclusions are as follows: (1) The force-displacement curve of the spherical-conical premium connection under dynamic load is similar to a parallelogram, and its energy dissipation is mainly caused by gross slip. The force-displacement curve under the modal vibration mode assumes a needle-like shape, and its energy dissipation is mainly caused by microslip. Compared with dynamic load, the energy dissipation value of the sealing surfaces under the first-order, the second-order modes is minimal. Therefore, the dynamic load impact on the sealing performance of the connection is greater than string vibration. (2) During the cyclic load, the spherical radius is positively correlated with energy dissipation but has little impact on the contact pressure of the sealing surfaces. Under the same axial cyclic load, it is impossible to analyze the influence of spherical radius on the sealing performance of connections by contact pressure, but the energy dissipation theory can effectively solve this problem. (3) Under different cyclic loads, the change amplitude of contact pressure is much smaller than that of energy dissipation. When analyzing the influence of cyclic load amplitude on the sealing performance of connections, the results tend to be conservative in terms of the contact pressure of sealing surfaces. (4) Under the combined action of internal pressure and axial cyclic load, it will cause plastic deformation accumulation on the thread and sealing surfaces. When the internal pressure and axial cyclic loads are 80 MPa, 400 kN, or 60 MPa, 500 kN, respectively, the thread with stress concentration is a potential danger point. At the same time, the energy dissipation significantly increases, leading to a decline in the sealing performance of the connection which has a risk of sealing failure. Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.