A Semi-Analytical Method for the Sealing Performance Prediction of Subsea Pipeline Compression Connector

: To predict the sealing performance of the subsea pipeline compression connector, a semi-analytical method is proposed and veriﬁed. The leakage condition is obtained as a function of the minimum radial deﬂection. The semi-analytical method consists of three parts: a macroscopic analytical model for the interference process deduced based on the membrane theory, a mesoscopic two-dimensional ﬁnite element analysis of the internally turned sealing surface according to scanning electron microscope observation, and a formulation of the leakage ratio according to the Abbott– Firestone curve. The inﬂuences of the geometry parameters are analyzed. The prediction and identiﬁcation of the leakage condition near the minimum preloaded deﬂection are investigated. Experimental tests are also carried out to verify the proposed method. The results show that, within ﬁve times the pipeline threshold of the thin wall thickness, compression connectors designed by the proposed method can create a reliable seal with a rough internally turned surface, saving 57% of the design time while reducing the machining time and costs.


Introduction
Non-welded connections are increasingly used for the maintenance of subsea pipelines in marine industries [1]. The compression connector is one of the most efficient connectors that belongs to the mechanically attached connection (ASTM F1387 [2], ISO 8434-1 [3]). This technology has been developed by many well-known companies, including the Phastite ® connector from Parker Hannifin (Columbus, OH, USA), the Tube-Mac ® connector from PYPLOK ® Corp. (Stoney Creek, ON, Canada), the compression connector from HAELOK ® (Schlieren, Switzerland), etc. The connector is directly attached to the pipe by the mechanical deformation and creates a seal and a permanent joint by elastic-plastic deformation. The sealing principle is a metal static seal that uses the plastic deformation of the metal to block the leakage channel.
In studies on the leak ratio on randomly rough surfaces, the topography formed by the turning texture shows unique sealing properties [4]. This can be traced back to groundbreaking studies on the surface waviness theory [5] and further developments on surface roughness [6] and sealing design factors [7]. In the aspect of the gasket, Nakamura and Funabashi [8] observed the leakage on a spiral-like topography. Then, by establishing a simplified semi-empirical leakage model, Polycarpou and Etsion [9] demonstrated the importance of micro-profile parameters for the sealing performance of the metal gasket seal. To further understand the sealing performance of the turning topography, Geoffroy and Prat [10] established the seal leakage model of the face-turned metal gasket and found the transition from the radial leakage to the spiral leakage by studying the radial and circumferential leakage paths. Liao et al. [11] first presented the simplification of the radial leakage path and the radial leakage model to clarify the effects of the micro-profile parameters on the sealing behaviors. More recently, the influence of wavy deformations and the formation of synthetic "turned-like" surfaces have been extensively studied. Scholars have also conducted much research on sealing topography related to the studies of asperity [12] and the probabilistic models [13] of rough surfaces. By establishing a 2D stochastic model, Marie et al. [14] developed a practical model and applied the multiscale calculation method based on the leakage ratio of the sealing structure. Liao et al. [15] proposed a method to decompose 3D surfaces into multi-scale surfaces by using wavelets to predict surface functions and detect machining errors. Deltombe et al. [16] used multi-scale analysis to characterize sealing performance. By comparing the differences between machining and superfinishing, it was shown that micro-roughness plays a major role in the leakage. In addition, Bricaud et al. [17] studied the leakage between firmly pressed rigid surfaces by measuring mass flow under different contact geometries using leakage tests.
Scholars have also conducted much research on compression connections. Weddeling et al. [18] derived an analytical model to determine the charging energy to a clamped connection with a mandrel inside the pipe. In another study, Weddeling et al. [19] developed a formula for the transferable load of connection with different groove shapes using membrane theory. Wei et al. [20] deduced the formula of the sealing pressure of a compression connector based on the superposition theorem of elasticity, and obtained the sealing condition using the finite element method (FE). Henriksen et al. [21] used experimental data and FE simulations to describe and understand the forces and deformations during the connection process and identified the increase radial stiffness as the best indicator of a clamped pipe flange connection while also predicting the load capacity of the connection. Wang et al. [22] proposed an optimization method for a compression connector based on a static metal sealing mechanism and optimized the structure using an FE model with the zero-order method. Onyegiri et al. [23] investigated a compression connector used on sandwich pipes combining 2D axisymmetric and 3D FE models. Quispe et al. [24] proposed a 2D FE model for threaded compression connection with square teeth and a metal static seal to join subsea sandwich pipes. Yan and Fan [25] employed a multiscale model to study the effect of the fluid pressure on the sealing performance of a pipeline connection. Agrawal et al. [26] developed a new method for a compression connector and carried out the extrusion test and compression test, concluding that the connection strength is considerable with welding. Liu et al. [27] optimized the anchoring structure by using a nonlinear FE analysis method. Wu et al. [28] introduced an axial load-bearing capacity calculation method for steel pipe-to-sleeve grouted connections, which also showed the significance of the compression connection.
For the compression connection, there is no simple analytical method to determine the achievable joint strength based on the process and joint characteristics [19,29]. In practice, the FE method is still the only option in many cases, resulting in a long design cycle and high test costs. Better knowledge and understanding of the elasticity and plastic deformation of the compression connection is required. In addition, most of the research on metal seals are about face-turned seals. Few studies can be found on the internally turned seal. The purpose of this study is to propose an effective semi-analytical method to predict the leakage ratio on the internally turned sealing surface of the compression connector, and, by analyzing the sealing behavior on the mesoscopic scale, to predict the leakage ratio for all standard pipe series.
The paper is organized as follows: In Section 2, the structure and fitting process of the compression connector are explained. In Section 3, an analytical model of the radial interference process is derived based on the membrane theory and verified by a macroscopic FE model. In Section 4, the SEM observation of the internally turned sealing surface is performed, and a mesoscopic FE model for the seal is developed. Then, the change law of the sealing surface during the compression is analyzed. In Section 5, by transferring the deflection into the mesoscopic FE model, and using the Abbott-Firestone curve, the leakage ratio formulation is derived. In Section 6, a parameter analysis and experimental tests are presented.

Connector Fitting Process
The pipeline compression connector, the fitting tools, and the fitting process are shown in Figure 1. The connector is designed to maintain and/or repair typical single-wall subsea pipes with a diameter between 2 and 12 inches. deflection into the mesoscopic FE model, and using the Abbott-Firestone curve, the leakage ratio formulation is derived. In Section 6, a parameter analysis and experimental tests are presented.

Connector Fitting Process
The pipeline compression connector, the fitting tools, and the fitting process are shown in Figure 1. The connector is designed to maintain and/or repair typical single-wall subsea pipes with a diameter between 2 and 12 inches. As shown in Figure 1, the connector consists of a base and a pair of press rings, which are made of metal and are symmetrically arranged along the pipe. Inside the base, there are one or more ring-shaped ridges for sealing and anchors (optional configuration) for gripping. A hydraulic fitting tool is used to drive the press rings to create a static compression metal seal [30]

Analytical Model of the Fitting Process
During the fitting process, external loads are transferred through the press ring, the shaped ridges are compressed on the sealing surface of the pipe and partially press into the pipe surface, creating a form-fit [31] sealing condition. The geometries can be idealized as 2D axisymmetric thin-walled cylinders [32] and the membrane theory of symmetrically loaded cylinders based on Kirchhoff's hypothesis [33] is applicable. The internal forces can be solved by the physical relations between stress and strain, the differential equations of the deflections, and the relations between strain and displacement.

Formulation
A shell element of the thin-walled long cylinder with mid-surface radius and thickness is shown in Figure 2. The mid-surface of the shell is located at the cylindrical coordinates system with angle , axial position x, and radius . The corresponding displacements on the mid-surface are , , and .
is the shearing force, and are the membrane forces, and and are the bending moments. and are constant in the circumference. As shown in Figure 1, the connector consists of a base and a pair of press rings, which are made of metal and are symmetrically arranged along the pipe. Inside the base, there are one or more ring-shaped ridges for sealing and anchors (optional configuration) for gripping. A hydraulic fitting tool is used to drive the press rings to create a static compression metal seal [30].

Analytical Model of the Fitting Process
During the fitting process, external loads are transferred through the press ring, the shaped ridges are compressed on the sealing surface of the pipe and partially press into the pipe surface, creating a form-fit [31] sealing condition. The geometries can be idealized as 2D axisymmetric thin-walled cylinders [32] and the membrane theory of symmetrically loaded cylinders based on Kirchhoff's hypothesis [33] is applicable. The internal forces can be solved by the physical relations between stress and strain, the differential equations of the deflections, and the relations between strain and displacement.

Formulation
A shell element of the thin-walled long cylinder with mid-surface radius R and thickness δ is shown in Figure 2. The mid-surface of the shell is located at the cylindrical coordinates system with angle ϕ, axial position x, and radius ρ(= R + z). The corresponding displacements on the mid-surface are u, v, and w. P x is the shearing force, N ϕ and N x are the membrane forces, and M ϕ and M x are the bending moments. N ϕ and M ϕ are constant in the circumference. Assuming that external forces only contribute to the normal pressure on the surface, the following geometrical and equilibrium equations apply: The expressions for the strains can be written as and, when applying Hook's law and yield, With the load of applied, the basic differential equation of the symmetrically loaded cylinder can be written as where 4 , 12 1 .
The general solution of the constant coefficients equation is cos sin cos sin , in which is a particular solution of Equation (4) and to are the constants of integration that need to be determined in each special case by the conditions at the cylinder ends.
The forces and moments for are Both the base and the pipe can be considered long cylinders and subject to the radial external load, , which is symmetrically distributed about the cylinder central line. Fixing the origin in the cross-section where is applied, the constants of Equation (6) are Figure 2. A thin-walled cylindrical shell in the cylindrical coordinates.
Assuming that external forces only contribute to the normal pressure on the surface, the following geometrical and equilibrium equations apply: The expressions for the strains can be written as and, when applying Hook's law and yield, With the load of P x applied, the basic differential equation of the symmetrically loaded cylinder can be written as where The general solution of the constant coefficients equation is in which f (x) is a particular solution of Equation (4) and C 1 to C 4 are the constants of integration that need to be determined in each special case by the conditions at the cylinder ends. The forces and moments for w are Both the base and the pipe can be considered long cylinders and subject to the radial external load, P x , which is symmetrically distributed about the cylinder central line. Fixing the origin in the cross-section where P x is applied, the constants of Equation (6) are and the radial deflection w can be derived as Then, the deflection, slope, shearing load, and bending moment can be rewritten as According to Equation (9), the influenced axial range of the cylinder can be predicted as The half of the 2D axisymmetric connector-pipe assembly with the radial equivalent load, P z , of the press ring is shown in Figure 3 with regard to the symmetry about the central plane of the connector. The axial distance of the ridge and the anchor is much greater than the half influence range, L 1 χ p+ ; hence, the coupling effect between the ridge and the anchor can be neglected in the following analysis. In the fitting process, the applied equivalent load, P z , can be divided into two components: one is for eliminating the maximum radial clearance, h, and the other is for compressing the ridge into the pipe. Meanwhile, the ridge is subjected to the reaction load, P r , of the pipe. Due to the wedge surfaces inside the press ring, there is a short axial distance, L 0 , between P z and P r .
According to Equation (9), the influenced axial range of the cylinder can be predicted as The half of the 2D axisymmetric connector-pipe assembly with the radial equivalent load, , of the press ring is shown in Figure 3 with regard to the symmetry about the central plane of the connector. The axial distance of the ridge and the anchor is much greater than the half influence range, ≫ ; hence, the coupling effect between the ridge and the anchor can be neglected in the following analysis. In the fitting process, the applied equivalent load, , can be divided into two components: one is for eliminating the maximum radial clearance, ℎ, and the other is for compressing the ridge into the pipe. Meanwhile, the ridge is subjected to the reaction load, , of the pipe. Due to the wedge surfaces inside the press ring, there is a short axial distance, , between and . When cos sin , the deflection of the mid-surface of the cylindrical base, , in the radial direction can be expressed as When = , can be rewritten as: Figure 3. Free-body diagram of the assembly subjected to the radial equivalent load, P z .
When ψ(β b x) = e −β b x (cos β b x + sin β b x), the deflection of the mid-surface of the cylindrical base, w b1 , in the radial direction can be expressed as When χ 0 = β b L 0 , P b1 can be rewritten as: Once the ridge contacts the pipe surface, a sealing pressure occurs on the pipe surface. When the pressure reaches the pipe's elastic limit, σ sp , the reaction load, P r , can be given as In this case, according to the Equation (10), the pipe deflection is The deflection of the base is Substituting P r from Equation (15) into Equation (16) yields The total equivalent load on the base is Finally, the deflection can be rewritten as In particular, the radial deflection at the point of P z is The bending moment and slope can also be derived with the same method. Multiple ridges can be designed according to the superposition principle.

Numerical Modeling and Validation
The finite element method (FE) can be used in the design of the compression connector. In order to validate the proposed macroscopic analytical solutions, a macroscopic 2D axisymmetric FE model is created for the simulation using ABAQUS ® /Standard [34].
The base and pipe are idealized as cylinders with a given thickness. An isotropic bilinear hardening model (BISO) with von Mises yield criterion [35] is adopted in which the yield plateau is ignored. The constitutive relationship of the plastic stage takes the form of where σ T is the true stress, σ s is the yield strength, and ε T is the accumulated plastic strain. ξ is the hardening modulus, which can be determined by the tangent modulus, E t , and Young's modulus, E, with the relation of For general applicability, E t is averaged and rounded to 5.0 × 10 2 Mpa by referring to the properties of wildly used pipe steels: mild steels, Q235, X60, X70, and X80. The geometrical parameters and material properties used in the paper are listed in Table 1.
The FE model is shown in Figure 4. To consider the large mesh distortion without twist, the four-node bilinear axisymmetric quadrilateral element with reduced integration and hourglass control (CAX4R) is used. The mesh sensitivity is examined. The right corners of the base and the pipe are fixed. For general applicability, is averaged and rounded to 5.0 × 10 2 MPa by referring to the properties of wildly used pipe steels: mild steels, Q235, X60, X70, and X80. The geometrical parameters and material properties used in the paper are listed in Table 1.
The FE model is shown in Figure 4. To consider the large mesh distortion without twist, the four-node bilinear axisymmetric quadrilateral element with reduced integration and hourglass control (CAX4R) is used. The mesh sensitivity is examined. The right corners of the base and the pipe are fixed. Between the analytical solution and the FE model, a realistic assumption should be made. Because the analytical solution satisfies the Kirchhoff's hypotheses, one of the hypotheses is that the thickness of the cylinder, , is much smaller than the radius of the mid-surface, (the threshold is / 0.05). In practice, the / of some pipes is higher than 0.05, e.g., for six-inch pipes, the range is / ∈ [0.04, 0.54] (ASME B36 10M [36] and 19M [37]). To satisfy the maximum mean error of less than 10% between the analytical model and the FE model, in this study, it is assumed that the pipe thickness in the FE simulation is within five times the thin wall threshold; hence, 20 mm and / 0.25. During the macroscopic FE analysis, an equivalent radial external load, , is applied on the appropriate width of the external surface of the base part and is gradually loaded from one to two times the theoretical threshold. The mid-surface deflections of both the base, , and the pipe, , are observed, as shown in Figure 5. Between the analytical solution and the FE model, a realistic assumption should be made. Because the analytical solution satisfies the Kirchhoff's hypotheses, one of the hypotheses is that the thickness of the cylinder, δ, is much smaller than the radius of the mid-surface, R (the threshold is δ/R ≤ 0.05). In practice, the δ p /R p of some pipes is higher than 0.05, e.g., for six-inch pipes, the range is δ p /R p ∈ [0.04, 0.54] (ASME B36 10M [36] and 19M [37]). To satisfy the maximum mean error of less than 10% between the analytical model and the FE model, in this study, it is assumed that the pipe thickness in the FE simulation is within five times the thin wall threshold; hence, δ p = 20 mm and δ p /R p ≈ 0.25. During the macroscopic FE analysis, an equivalent radial external load, P z , is applied on the appropriate width of the external surface of the base part and is gradually loaded from one to two times the theoretical threshold. The mid-surface deflections of both the base, w x , and the pipe, w p , are observed, as shown in Figure 5. The deflections of the analytical method are consistent with the FE results. In the axial direction, the mean error of with the analytical method is 8.5% compared to that of the FE method, and the mean error of is 5.7%. The error should usually be within 10% [38]. In particular, within an axial distance of 20 mm from the seal (at x = 0), the mean error of is only 3.7%. However, the corresponding error of is much higher due to the local plastic deformation on the sealing surface under mesoscopic interference. With the help of a mesoscopic FE model, the sealing process can be analyzed, as discussed in Section 4.

FE model of the Sealing Surface
In this section, a mesoscale FE model is created to examine the plastic deformation of the internally turned surface during the fitting process.

Characterization of the Sealing Surface
The turning process is widely used in the manufacturing of flange surfaces, and these surfaces usually reveal spiral morphology [11,39]. Similarly, due to the difficulty of polishing a thin-walled cylindrical shell [40], the topography produced by the rough turning process is usually left on the ridge surface. The surface was observed using scanning electron microscopy (SEM), as shown in Figure 6. The internally turned surface has a distinctive anisotropic topography, which is characterized by regular helical and continuous streaks with slight defects and irregularities. Therefore, there are two potential leakage paths: the first one (Direction 1) is the longest helical path along the streak, and the second (Direction 2) is the shortest leak path in the axial direction. The deflections of the analytical method are consistent with the FE results. In the axial direction, the mean error of w x with the analytical method is 8.5% compared to that of the FE method, and the mean error of w p is 5.7%. The error should usually be within 10% [38]. In particular, within an axial distance of 20 mm from the seal (at x = 0), the mean error of w x is only 3.7%. However, the corresponding error of w p is much higher due to the local plastic deformation on the sealing surface under mesoscopic interference. With the help of a mesoscopic FE model, the sealing process can be analyzed, as discussed in Section 4.

FE Model of the Sealing Surface
In this section, a mesoscale FE model is created to examine the plastic deformation of the internally turned surface during the fitting process.

Characterization of the Sealing Surface
The turning process is widely used in the manufacturing of flange surfaces, and these surfaces usually reveal spiral morphology [11,39]. Similarly, due to the difficulty of polishing a thin-walled cylindrical shell [40], the topography produced by the rough turning process is usually left on the ridge surface. The surface was observed using scanning electron microscopy (SEM), as shown in Figure 6. The deflections of the analytical method are consistent with the FE results. In the axial direction, the mean error of with the analytical method is 8.5% compared to that of the FE method, and the mean error of is 5.7%. The error should usually be within 10% [38]. In particular, within an axial distance of 20 mm from the seal (at x = 0), the mean error of is only 3.7%. However, the corresponding error of is much higher due to the local plastic deformation on the sealing surface under mesoscopic interference. With the help of a mesoscopic FE model, the sealing process can be analyzed, as discussed in Section 4.

FE model of the Sealing Surface
In this section, a mesoscale FE model is created to examine the plastic deformation of the internally turned surface during the fitting process.

Characterization of the Sealing Surface
The turning process is widely used in the manufacturing of flange surfaces, and these surfaces usually reveal spiral morphology [11,39]. Similarly, due to the difficulty of polishing a thin-walled cylindrical shell [40], the topography produced by the rough turning process is usually left on the ridge surface. The surface was observed using scanning electron microscopy (SEM), as shown in Figure 6. The internally turned surface has a distinctive anisotropic topography, which is characterized by regular helical and continuous streaks with slight defects and irregularities. Therefore, there are two potential leakage paths: the first one (Direction 1) is the longest helical path along the streak, and the second (Direction 2) is the shortest leak path in the axial direction. The internally turned surface has a distinctive anisotropic topography, which is characterized by regular helical and continuous streaks with slight defects and irregularities. Therefore, there are two potential leakage paths: the first one (Direction 1) is the longest helical path along the streak, and the second (Direction 2) is the shortest leak path in the axial direction.

FE Model of the Sealing Surface
According to the turning topography [41,42] and ISO recommendations (ISO 16610-41 [43] and 16610-85 [44]), the surface profile can be considered a periodic connection from the nose profile of the turning tool. Consistently with the SEM results, the mesoscale sealing surface profile is created, as shown in Figure 7.

FE model of the Sealing Surface
According to the turning topography [41,42] and ISO recommendations (ISO 16610-41 [43] and 16610-85 [44]), the surface profile can be considered a periodic connection from the nose profile of the turning tool. Consistently with the SEM results, the mesoscale sealing surface profile is created, as shown in Figure 7. The turning process produces a helical profile on the inner surface in 3D. Chen and Shih [45] have demonstrated that the 2D axisymmetric model can be used to reduce the computational complexity of the helical 3D model while meeting the accuracy requirement. Therefore, a 2D axisymmetric FE model is used. According to Saint-Venant's principle [46], the computation domain is limited to the localized contact region with a portion of the base and is set as 0 mm. On the sealing surface, since leakage is more of a concern than stress concentration, a groove is positioned at the center and streaks are arranged at both ends. The topography layout is set to an acceptable minimum arrangement of six streaks with five equally spaced grooves. The cross-section of the groove is assumed to be part of a constant circle according to the ISO recommendation. For consistency, additional FE parameters are inherited from the previous macroscopic FE model (in Section 3.2). Since the pipe surface is preprocessed before connection, the surface roughness is 1-2 orders of magnitude smaller than the ridge, and it is simplified to an even surface (no roughing) in the analysis. To show subtle deformation and reduce computation, the mesh is optimized. The sensitivity of the mesh is examined. The following parameters for the study case are provided: the groove radius, 0.3 mm; the pitch of the streaks, 0.5 mm; the sealing width, ζ 3 mm.

Change Laws of the Leakage Channel
The interference can be obtained by applying an equivalent deflection/displacement to the ridge in the simulation. According to the FE simulation, the stress contour plots for the corresponding deformation phases are shown in Figure 8a. The deflections of the pipe surface, , the pipe mid-surface, , the streaks, , and the embedded depth of the streaks, ℎ , are plotted in Figure 8b. The contact pressure on the sealing surface is shown in Figure 8c. The turning process produces a helical profile on the inner surface in 3D. Chen and Shih [45] have demonstrated that the 2D axisymmetric model can be used to reduce the computational complexity of the helical 3D model while meeting the accuracy requirement. Therefore, a 2D axisymmetric FE model is used. According to Saint-Venant's principle [46], the computation domain is limited to the localized contact region with a portion of the base and is set as L 0 = 0 mm. On the sealing surface, since leakage is more of a concern than stress concentration, a groove is positioned at the center and streaks are arranged at both ends. The topography layout is set to an acceptable minimum arrangement of six streaks with five equally spaced grooves. The cross-section of the groove is assumed to be part of a constant circle according to the ISO recommendation. For consistency, additional FE parameters are inherited from the previous macroscopic FE model (in Section 3.2). Since the pipe surface is preprocessed before connection, the surface roughness is 1-2 orders of magnitude smaller than the ridge, and it is simplified to an even surface (no roughing) in the analysis. To show subtle deformation and reduce computation, the mesh is optimized. The sensitivity of the mesh is examined. The following parameters for the study case are provided: the groove radius, r c = 0.3 mm; the pitch of the streaks, ζ λ = 0.5 mm; the sealing width, ζ = 3 mm.

Change Laws of the Leakage Channel
The interference can be obtained by applying an equivalent deflection/displacement to the ridge in the simulation. According to the FE simulation, the stress contour plots for the corresponding deformation phases are shown in Figure 8a. The deflections of the pipe surface, w s , the pipe mid-surface, w p , the streaks, w k , and the embedded depth of the streaks, h b , are plotted in Figure 8b. The contact pressure on the sealing surface is shown in Figure 8c.
As shown in Figure 8a, both of the streaks of the ridge and the pipe surface are deformed during interference. The local upheaval of the pipe surface further fills the grooves and leads to the change in the cross-section of the leakage channels.
As shown in Figure 8b, the deflection of the pipe surface, w s , which is distinct from other deflections, has a plateau phase during the embedding, and correspondingly, the embedded depth, h b , levels off just at the end of the plateau. Hence, the plateau indicates the interference process. The process can be divided into three phases: synchronous (in this case, w b2 ∈ (0, 0.06)), interference (w b2 ∈ (0.6, 0.16)), and form-fit (w b2 > 0.16). As shown in Figure 8a, both of the streaks of the ridge and the pipe surface are deformed during interference. The local upheaval of the pipe surface further fills the grooves and leads to the change in the cross-section of the leakage channels.
As shown in Figure 8b, the deflection of the pipe surface, , which is distinct from other deflections, has a plateau phase during the embedding, and correspondingly, the embedded depth, ℎ , levels off just at the end of the plateau. Hence, the plateau indicates the interference process. The process can be divided into three phases: synchronous (in this case, 2 ∈ (0, 0.06)), interference ( 2 ∈ (0.6, 0.16)), and form-fit ( 2 > 0. 16). As shown in Figure 8c, the contact pressure is distributed axially with stress concentration areas (corresponding to streaks) and untouched areas (corresponding to the helical leakage channel).
The following change laws can be obtained: • In the synchronize phase, the ridge and pipe deflect synchronously, the streaks are barely embedded in the pipe surface, and the stress concentration occurs at the contact points. Materials are not yielded; • In the interference phase, streaks continue to be embedded in the pipe surface, the helical leakage channel is narrowed by the plastic deformations of the pipe, and the equispaced contact areas with high-contact pressure extend in the axial direction; • In the form-fit phase, embedding is mainly performed with constant ℎ ; however, the plastic deformation of the seal area (the difference between and ) continually increases as 2 increases, which makes the seal tighter. The untouched areas narrow considerably and will eventually disappear with the increase in As shown in Figure 8c, the contact pressure is distributed axially with stress concentration areas (corresponding to streaks) and untouched areas (corresponding to the helical leakage channel).
The following change laws can be obtained: • In the synchronize phase, the ridge and pipe deflect synchronously, the streaks are barely embedded in the pipe surface, and the stress concentration occurs at the contact points. Materials are not yielded; • In the interference phase, streaks continue to be embedded in the pipe surface, the helical leakage channel is narrowed by the plastic deformations of the pipe, and the equispaced contact areas with high-contact pressure extend in the axial direction; • In the form-fit phase, embedding is mainly performed with constant h b ; however, the plastic deformation of the seal area (the difference between w s and w p ) continually increases as w b2 increases, which makes the seal tighter. The untouched areas narrow considerably and will eventually disappear with the increase in w b2 .
To verify the availability of the mesoscale FE model, the analytical values and the mesoscale FE results need to be compared. For w x , the analytical preload value is 1.35 × 10 −2 mm, while the FE range in Figure 8a is from 1.20 × 10 −2 to 1.40 × 10 −2 mm; for w p , the analytical value is 1.54 × 10 −2 mm, while the FE range is from 0.26 × 10 −2 to 3.66 × 10 −2 mm. Therefore, the analytical variables of both the base and the pipe are included in the mesoscale FE model. The availability of the FE model is verified.

Semi-Analytical Model of the Leakage
In this section, the leakage ratio is derived from the macroscopic analytical model and the mesoscopic FE model using the Abbott-Firestone curve.
According to the FE results in Section 4.3, the leakage channel and the relevant dimensions are shown in the Figure 9. The height of the leakage paths in the radial direction is h a , and the embedded depth of the streaks is h b . Since the axial leakage vanishes during the interference, the direction of the leakage (in Section 4.1) is different before and after the interference [11,47]. Therefore, for the compression connector, the leakage ratio of the helical leak dominates and the boundary condition is limited to radial clearance h(w x ) < h a . 10 mm, while the FE range in Figure 8a is from 1.20 × 10 to 1.40 × 10 mm; for , the analytical value is 1.54 × 10 −2 mm, while the FE range is from 0.26 × 10 −2 to 3.66 × 10 −2 mm. Therefore, the analytical variables of both the base and the pipe are included in the mesoscale FE model. The availability of the FE model is verified.

Semi-Analytical Model of the Leakage
In this section, the leakage ratio is derived from the macroscopic analytical model and the mesoscopic FE model using the Abbott-Firestone curve.
According to the FE results in Section 4.3, the leakage channel and the relevant dimensions are shown in the Figure 9. The height of the leakage paths in the radial direction is ℎ , and the embedded depth of the streaks is ℎ . Since the axial leakage vanishes during the interference, the direction of the leakage (in Section 4.1) is different before and after the interference [11,47]. Therefore, for the compression connector, the leakage ratio of the helical leak dominates and the boundary condition is limited to radial clearance ℎ ℎ . To attribute the mesoscale profile parameters for calculating the seal leakage ratio, the Abbott-Firestone curve [48] is used. The effectiveness of the curve has been proved in describing the turning surface. The curve is shown in Figure 10. To attribute the mesoscale profile parameters for calculating the seal leakage ratio, the Abbott-Firestone curve [48] is used. The effectiveness of the curve has been proved in describing the turning surface. The curve is shown in Figure 10. The mathematical formula of the Abbott-Firestone curve can be represented by z(x), which can be mathematically acquired by the fitting method of a polynomial. However, according to the geometric characteristics of the FE model and the parameters in Table 1, z(x) can be given as A hypothetical 3D rectangular channel can be used as the leakage path in the mathematical model to calculate the leakage ratio, as shown in Figure 11. The mathematical formula of the Abbott-Firestone curve can be represented by z(x), which can be mathematically acquired by the fitting method of a polynomial. However, according to the geometric characteristics of the FE model and the parameters in Table 1, z(x) can be given as A hypothetical 3D rectangular channel can be used as the leakage path in the mathematical model to calculate the leakage ratio, as shown in Figure 11. The equivalent dimensions are width ζ , height ℎ and length , of which ζ and ℎ are regarded as definite values along .
represents the cross-sectional area of the hypothetical rectangular channel. Given the uncertainty of the pressure, the magnitude of the inlet pressure, , and the outlet pressure, , is assumed as . The equivalent dimensions are width ζ c , height h c and length L c , of which ζ c and h c are regarded as definite values along L c . A represents the cross-sectional area of the hypothetical rectangular channel. Given the uncertainty of the pressure, the magnitude of the inlet pressure, p i , and the outlet pressure, p o , is assumed as p i ≥ p o .
The integral value of the cross-sectional area of the streaks A o is Subjected to h(w x ), the true cross-section of the leakage channel can be calculated by the following relations: The equivalent height of the leakage channel can be expressed as The radius of the leakage path is given by The length of the leakage channel can be defined as The governing equation of the flow field based on the laminar flow can be defined using the Reynolds equations [49]: where q v is the volume flow ratio per unit width, η is the viscosity of the fluid, and p is the fluid pressure. The circumferential leakage ratio, Q c , can be represented in the following form: The leakage ratio, Q c , as a function of the radial deflection, w x , is given as: The following equation introduces the notation: The leakage ratio can be represented in the simplified form:

Parameter Analysis
The proposed semi-analytical method can reveal the relationship between the deflection and the leakage condition. Therefore, the key parameters are those that have an impact on the relationship as they influence deflection and are easy to implement in design. Via multi-scale parameter analysis, key design and performance rules can be identified, and the minimum preloaded deflection can be obtained.
On the macroscopic level, the effects of the radius and clearance on deflections are first investigated. Because R b is determined by R p , δ, and h, by taking R p as an independent variable while others are constant, the influences of R p on w x and w p at x = 0 are plotted with logarithmic coordinates using Equations (15) and (20), as shown in Figure 12a. However, according to Equations (12) and (13), h is positively correlated with the deflection component, w b1 , and its influence on the total deflection, w x , varies. As a reference for the design, with the increase in h, the corresponding change in h/w x and its gradient is calculated, as shown in Figure 13b. Note that w p is not shown in Figure 12b since it is irrelevant to h from Equation (15).
The leakage ratio can be represented in the simplified form:

Parameter Analysis
The proposed semi-analytical method can reveal the relationship between the deflection and the leakage condition. Therefore, the key parameters are those that have an impact on the relationship as they influence deflection and are easy to implement in design. Via multi-scale parameter analysis, key design and performance rules can be identified, and the minimum preloaded deflection can be obtained.
On the macroscopic level, the effects of the radius and clearance on deflections are first investigated. Because is determined by , , and ℎ , by taking as an independent variable while others are constant, the influences of on and at x = 0 are plotted with logarithmic coordinates using Equations (15) and (20), as shown in Figure 12a. However, according to Equations (12) and (13), ℎ is positively correlated with the deflection component, 1 , and its influence on the total deflection, , varies. As a reference for the design, with the increase in ℎ, the corresponding change in ℎ/ and its gradient is calculated, as shown in Figure 13b. Note that is not shown in Figure 12b since it is irrelevant to ℎ from Equation (15). The increment of is gradually related to ℎ, and, in this instance, when ℎ > 0.2 mm, the gradient of ℎ/ approaches zero; • There is a threshold of ℎ (when ℎ = 0.2) at which the growth of is directly proportional to the growth of ℎ.
Similarly to the radius, , the thickness, , also shows a prominent influence on the deflections. The thickness of the base, , and the pipe, , are analyzed under three  The following conclusions are drawn from Figure 13: • As the increases, all deflections decrease with a decreased change rate; • has a minimum influence on (the curve in Case @), and is irrelevant to , so the change in has a negligible impact on the interference. Therefore, consideration may be given to appropriately increasing or decreasing the thickness, , as needed.
In addition to the above parameters, deflections are also affected by the sealing • For all standard diameters in the ratio range of δ p /R p < 0.25, within four-inch diameters (R p ≤ 50 mm), changes of w x and w p can be considered independent of R p . The following conclusions are drawn from Figure 12b: • The increment of w x is gradually related to h, and, in this instance, when h > 0.2 mm, the gradient of h/w x approaches zero; • There is a threshold of h (when h = 0.2) at which the growth of w x is directly proportional to the growth of h.
Similarly to the radius, R, the thickness, δ, also shows a prominent influence on the deflections. The thickness of the base, δ b , and the pipe, δ p , are analyzed under three conditions: in Case 1, δ b is constant and δ p is variable; Case 2 has the opposite conditions to Case 1; in Case 3, δ b = δ p is considered. Note that, according to Equation (15), δ b is irrelevant to w p . The influences of δ on the deflections are shown in Figure 13. The curve in Case 1 (w x ) is the influence of δ p on w x while δ b = 20 mm, and the w p in Case 1 is the influence on w p . The w x in Case 2 is the influence of δ b on w x while δ p = 20 mm. The w x in Case 3 is the influence of δ p on w x while δ b = δ p , and the w p in Case 3 is the influence on w p .
The following conclusions are drawn from Figure 13: • As the δ p increases, all deflections decrease with a decreased change rate; • δ b has a minimum influence on w x (the curve w x in Case @), and is irrelevant to w p , so the change in δ b has a negligible impact on the interference. Therefore, consideration may be given to appropriately increasing or decreasing the thickness, δ b , as needed.
In addition to the above parameters, deflections are also affected by the sealing width. For generality, a ratio of the pipe thickness to the sealing width, δ p /ζ, is used. The deflection as a function of δ p /ζ is shown in Figure 14 with logarithmic coordinates. Moreover, since the ratio also affects flexural rigidity, the flexural rigidity of the pipe, D p , is also shown in the figure. For deflections, they are influenced by a threshold of the ratio: when / 1, deflections decrease as the ratio increases; when / 1, the deflection change is practically unrelated to the ratio. This shows that when the thickness exceeds the sealing width, / hardly affects deflection.
Next, the effects of and on deflections are also investigated as depicted in the logarithmic coordinates in Figure 15. Since is irrelevant to , its curve is not shown. • For the flexural rigidity, a linear relationship exists between the ratio and D p in logarithmic coordinates; • For deflections, they are influenced by a threshold of the ratio: when δ p /ζ ≤ 1, deflections decrease as the ratio increases; when δ p /ζ > 1, the deflection change is practically unrelated to the ratio. This shows that when the thickness exceeds the sealing width, δ p /ζ hardly affects deflection.
Next, the effects of D b and D p on deflections are also investigated as depicted in the logarithmic coordinates in Figure 15. Since w p is irrelevant to D b , its curve is not shown. As shown in Figure 15, exhibits distinct effects on and :  exhibits an almost linear relationship with (the curve wp-Dp) in logarithmic coordinates but a decreasing effect on (wx-Dp);  When 10 3 , all effects on deflections are approximately linear, but when 10 3 , the effect on and is significantly different;  When 10 4 , the only effect of the flexural rigidity is that of on ;  The effect on of the flexural rigidity of (wx-Db) is comparable to that of the pipe (wx-Dp).
On the mesoscopic level, the leakage condition can be analyzed. The influence of on the nonlinear change in the cross-sectional area of the leakage channel is investigated first. Based on the mesoscopic FE model (in Section 4), the true width, , and the equivalent height, ℎ , of are shown in Figure 16 with increments, where the corresponding radius of the leakage path, , is also depicted to indicate the relationship between macroscopic and mesoscopic changes. As shown in Figure 15, D p exhibits distinct effects on w p and w x : • D p exhibits an almost linear relationship with w p (the curve w p -D p ) in logarithmic coordinates but a decreasing effect on w x (w x -D p ); • When D p < 10 3 , all effects on deflections are approximately linear, but when D p ≥ 10 3 , the effect on w x and w p is significantly different; • When D p ≥ 10 4 , the only effect of the flexural rigidity is that of D p on w p ; • The effect on w x of the flexural rigidity of On the mesoscopic level, the leakage condition can be analyzed. The influence of w x on the nonlinear change in the cross-sectional area of the leakage channel A is investigated first. Based on the mesoscopic FE model (in Section 4), the true width, ζ c , and the equivalent height, h c , of A are shown in Figure 16 with w x increments, where the corresponding radius of the leakage path, R c , is also depicted to indicate the relationship between macroscopic and mesoscopic changes.
deflections decrease as the ratio increases; when / > 1, the deflection change is practically unrelated to the ratio. This shows that when the thickness exceeds the sealing width, / hardly affects deflection.
Next, the effects of and on deflections are also investigated as depicted in the logarithmic coordinates in Figure 15. Since is irrelevant to , its curve is not shown. As shown in Figure 15, exhibits distinct effects on and : • exhibits an almost linear relationship with (the curve wp-Dp) in logarithmic coordinates but a decreasing effect on (wx-Dp); • When < 10 3 , all effects on deflections are approximately linear, but when ≥ 10 3 , the effect on and is significantly different; • When ≥ 10 4 , the only effect of the flexural rigidity is that of on ; • The effect on of the flexural rigidity of (wx-Db) is comparable to that of the pipe (wx-Dp).
On the mesoscopic level, the leakage condition can be analyzed. The influence of on the nonlinear change in the cross-sectional area of the leakage channel is investigated first. Based on the mesoscopic FE model (in Section 4), the true width, , and the equivalent height, ℎ , of are shown in Figure 16 with increments, where the corresponding radius of the leakage path, , is also depicted to indicate the relationship between macroscopic and mesoscopic changes.  When w x > 0.45 mm, ζ c and h c are virtually unchanged, indicating that the leakage ratio is stable.
(η = 1.77 × 10 5 Pa·s), air (η = 1.01 × 10 3 Pa·s), and oil (η = 1.50 Pa·s), with inlet pressure, p i = 1 MPa, and outlet pressure, p o = 101 KPa, can be calculated as a function of w x , as shown in Figure 17 in the logarithmic coordinates. The corresponding R c is also shown. Other unlisted parameters are consistent with the above analyses.
• Coincidentally with the macroscopic results in Section 4.3, the true height of the channel (ℎ − ℎ ) and ℎ also have a notable reduction during the interference within the range ∈ (0.34, 0.45), while ℎ shows a clear descent stage; • When > 0.45 mm, and ℎ are virtually unchanged, indicating that the leakage ratio is stable.
Based on the above findings, the minimum preloaded deflection can be predicted. Since varies with different sealed mediums, the typical leakage ratios of water ( = 1.77 × 10 5 Pa· s), air ( = 1.01 × 10 3 Pa· s), and oil ( = 1.50 Pa· s), with inlet pressure, = 1 MPa, and outlet pressure, = 101 KPa, can be calculated as a function of , as shown in Figure 17 in the logarithmic coordinates. The corresponding is also shown. Other unlisted parameters are consistent with the above analyses. As shown in Figure 17, the air leakage ratio is distinct from other mediums, which makes a great difference to the sealing performance of connectors. Accordingly, the following conclusions are drawn: • For air medium, according to the highest Class VI of ASME B16.104 [50], the connector can be considered leak-tight below, = 0.24 L/h in the study case, while the corresponding minimum is in the range from 0.37 to 0.38 mm (0.25 to 0.10 L/h, shadowed area in the figure); • Similar to the equivalent height curve, ℎ , in Figure 16, there is also an accelerated descent range of the leakage ratio and all leakage ratios are reduced to a minimum with > 0.45 mm; • The macroscopic minimum, , calculated from the proposed analytical model is 0.42 mm (the dash-double-dotted line in the figure), which is covered by the ASME criteria; • The ASME criteria and the theoretical minimum, , are included within the leakage decrease range (from 0.37 to 0.45). Therefore, the results of the semi-analytical method are self-consistent, and the case connector is theoretically leak-tight to all typical mediums. The theoretical feasibility of the proposed semi-analytical method to predict the minimum preloaded deflection is verified.
At the end of the parameter analysis, considering the negligible effect of the thickness, the leakage ratio as a function of the rate, / , is explored in Figure 18 with logarithmic coordinates. The corresponding and ASME criteria of the air leakage ratio (from 0.53 to 0.20 L/h, shadowed area in the figure) are also shown. 84 Leak-tight As shown in Figure 17, the air leakage ratio is distinct from other mediums, which makes a great difference to the sealing performance of connectors. Accordingly, the following conclusions are drawn:

•
For air medium, according to the highest Class VI of ASME B16.104 [50], the connector can be considered leak-tight below, Q c = 0.24 L/h in the study case, while the corresponding minimum w x is in the range from 0.37 to 0.38 mm (0.25 to 0.10 L/h, shadowed area in the figure); • Similar to the equivalent height curve, h c , in Figure 16, there is also an accelerated descent range of the leakage ratio and all leakage ratios are reduced to a minimum with w x > 0.45 mm; • The macroscopic minimum, w x , calculated from the proposed analytical model is 0.42 mm (the dash-double-dotted line in the figure), which is covered by the ASME criteria; • The ASME criteria and the theoretical minimum, w x , are included within the leakage decrease range (from 0.37 to 0.45). Therefore, the results of the semi-analytical method are self-consistent, and the case connector is theoretically leak-tight to all typical mediums. The theoretical feasibility of the proposed semi-analytical method to predict the minimum preloaded deflection is verified.
At the end of the parameter analysis, considering the negligible effect of the thickness, the leakage ratio as a function of the rate, δ p /ζ, is explored in Figure 18 with logarithmic coordinates. The corresponding R c and ASME criteria of the air leakage ratio (from 0.53 to 0.20 L/h, shadowed area in the figure) are also shown.
The following conclusions can be drawn from Figure 18: • When δ p /ζ ∈ (1,3], w x and Q c begin descending with the increase in δ p /ζ; • When δ p /ζ ∈ (3,5], the leakage is noticeably decreasing; • When δ p /ζ > 5, the effects on w x , Q c , and r c are moderated, the leakage is stable, and the sealing state needs to be further verified by prototype tests; • When δ p /ζ ≤ 3, the leakage may occur on the internally turned topography, the recommended design value is δ p /ζ ≥ 5. The following conclusions can be drawn from Figure 18:  When / ∈ (1, 3], and begin descending with the increase in / ;  When / ∈ (3,5], the leakage is noticeably decreasing;  When / 5, the effects on , , and are moderated, the leakage is stable, and the sealing state needs to be further verified by prototype tests;  When / 3, the leakage may occur on the internally turned topography, the recommended design value is / 5.

Experimental Test
In order to verify the feasibility of the proposed semi-analytical method on the sealing performance prediction, according to ISO 21329:2004 [51], a working pressure test and a burst pressure test are carried out. All tests are conducted under the temperature of 25 °C and the maximum gauge pressure is 60 MPa with an accuracy rating of 0.4 (0.24 MPa). Furthermore, a customized 80 MPa manual hydraulic pump and an air compressor are used. The schematic diagram of the experimental setup is shown in Figure 19.
Leak-tight Figure 18. Leakage ratios as a function of rate δ p /ζ.

Experimental Test
In order to verify the feasibility of the proposed semi-analytical method on the sealing performance prediction, according to ISO 21329:2004 [51], a working pressure test and a burst pressure test are carried out. All tests are conducted under the temperature of 25 • C and the maximum gauge pressure is 60 MPa with an accuracy rating of 0.4 (0.24 MPa). Furthermore, a customized 80 MPa manual hydraulic pump and an air compressor are used. The schematic diagram of the experimental setup is shown in Figure 19. Prototypes are developed based on the proposed semi-analytical approach using seamless wrought steel pipes (ASME B36 10M [36]) of 3 to 8 inches. Geometrical parameters are listed in Table 2. A few redesigned prototypes are shown in Figure 20.  Prototypes are developed based on the proposed semi-analytical approach using seamless wrought steel pipes (ASME B36 10M [36]) of 3 to 8 inches. Geometrical parameters are listed in Table 2. A few redesigned prototypes are shown in Figure 20.

Working Pressure Test
The purpose of the test is to verify whether the designed prototype can reach the designed working pressure. The mediums are clean air and water. For the air pressure test, the prototype is pressurized with air at the required depth of an experimental tank. During the working pressure test, the internal pressure gradually increases to the required working pressure, rising below a maximum rate of 100 MPa/min. When the working pressure is reached, to stabilize the system and periodically check the assembly for the leakage it is necessary to maintain constant pressure for a sufficient time (10 min for air pressure and 24 h for hydraulic pressure in the study). Neither structural failure nor visual indications of the leakage are permitted. The setups of the air pressure test and the hydraulic pressure test are shown in Figure 21. A sectional view of the tested prototype is shown in Figure 22.

Working Pressure Test
The purpose of the test is to verify whether the designed prototype can reach the designed working pressure. The mediums are clean air and water. For the air pressure test, the prototype is pressurized with air at the required depth of an experimental tank. During the working pressure test, the internal pressure gradually increases to the required working pressure, rising below a maximum rate of 100 MPa/min. When the working pressure is reached, to stabilize the system and periodically check the assembly for the leakage it is necessary to maintain constant pressure for a sufficient time (10 min for air pressure and 24 h for hydraulic pressure in the study). Neither structural failure nor visual indications of the leakage are permitted. The setups of the air pressure test and the hydraulic pressure test are shown in Figure 21. A sectional view of the tested prototype is shown in Figure 22.  From Figure 22, the internal deflection of the prototypes conforms to the design expectation, there are no visual indications of leakage during sampling inspection with different sealing widths under variable structure configurations, and specimens are in good contact with seals. The sealing surface condition is not detailed here since it has already From Figure 22, the internal deflection of the prototypes conforms to the design expectation, there are no visual indications of leakage during sampling inspection with different sealing widths under variable structure configurations, and specimens are in good contact with seals. The sealing surface condition is not detailed here since it has already been discussed in our previous studies [20,22].
Typical test results for the prototypes are shown in Table 3, indicating that all prototypes are leak tight. Therefore, the prototypes developed by the proposed semi-analytical method meet the working pressure requirement. The purpose of the test is to find the failure pressure of the designed prototypes. To avoid the risk of the compression air, only the hydraulic pressure was applied. Following the leak test, the pressure is continuously raised below the maximum pressurization rate, and the pressure is maintained for 5 min per 1 MPa until there is evidence of the structural failure or leakage. The mean failure pressures for the prototypes are listed in Table 4. For comparison, the results of the same prototype types with polished sealing surfaces from our previous design by FE methods [20,22] are also listed.  Table 4, there is no significant difference between the analytical design and the previous FE design, indicating that the semi-analytical solution does not have a considerable effect on the previous solution. In addition, the experiments did not show any leakage. All failures are structural, and most are manifested as the axial pipe movement. In the failure, the internal pressure no longer increased and decreased slightly in some tests. This type of failures involves axial anchoring failure and pull-out deformation, which will be discussed in a follow-up study.
To summarize, the results show that even with a rough radial sealing surface, the compression connector designed by the proposed method can quickly create a reliable radial seal. On the other hand, the computational time has been reduced, saving more than 57% design time in practice. Meanwhile, the internal finish turning machining is no longer necessary, which further reduces the time and cost of the manufacturing process.

Conclusions
A semi-analytical method for the sealing performance prediction of the subsea pipeline compression connector is proposed in this paper. Prototypes are fabricated on the basis of the design results with this method and the failure pressures are tested in lab. The following conclusions can be drawn: An analytical model of the fitting process of the compression connector is derived, based on the membrane theory, in which the assembly is idealized as pairs of thin-walled cylinders under uniformly distributed radial load. The expressions of the deflections with the contact pressure are obtained;

1.
A macroscopic 2D FE model is created to verify the analytical model. An agreement is found between the analytical and FE results with a mean error of less than 8.5%. The analytical model is proven to be good enough with the thickness-radius ratio being no more than 0.25 (δ/R ≤ 0.25); 2.
The geometric parameter effects on the interference process are investigated. The influence and thresholds of the geometric parameters on sealing performance are presented, which provide the design guidance for the subsea pipeline compression connector; 3.
The radial sealing surface of the connector is observed by SEM. The surface is an anisotropic topography characterized by regular helical and continuous streaks with slight defects and irregularities. Potential leak paths in two orthogonal directions are identified as axial and circumferential; 4.
A mesoscopic 2D FE model is created based on the SEM result and the macroscopic FE model. The effects of the leakage channel in the sealing process are analyzed. Then, the analytical minimum deflections and those of FE are compared and shown to be consistent with each other. Since the macroscopic analytical minimum deflections of both the base and the pipe are included in the mesoscale results, it is theoretically verified that the semi-analytical method can identify and predict the leakage condition of the connector; 5.
The Abbott-Firestone curve of the sealing surface is introduced based on SEM observation and the FE model. A semi-analytical leakage model was obtained. Since the true height of the leakage channel shows a clear descent stage, the minimum preloaded deflection of the leak-tight range is identified. On the rough internally turned sealing surface, leakage may occur when the ratio of the pipe thickness to the sealing width is δ p /ζ < 3, and the recommended leak-tight design value is δ p /ζ > 5; 6.
The results with the proposed semi-analytical method correlate well with the working pressure test and the burst pressure test. The designed connectors can quickly create a reliable radial metal static seal with the rough turning topography and save more than 57% of design time compared to the common FE method, while reducing machining time and costs. 7.
The proposed semi-analytical method can be used to optimize the overall structural performance of compression connectors with a low computational cost. It allows the prediction and identification of the leakage condition of all standard pipe series with δ p /R p ≤ 0.25. The model also provides a good understanding of the radial sealing process. These results are helpful to highlight the different sealing performances with other type of seals and guide the further improvements.
In order to perform further quantitative analyses of the leakage ratio, it is necessary to explore the tangential contact mechanics and the deflection behaviors of the connector in subsequent studies.

Data Availability Statement:
The data used to support the findings of this study are included within the article.

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

A
Cross-sectional area of the hypothetical rectangle channel (mm 2 ) A o Integral value of the total Abbott-Firestone curve (mm 2 ) D, D b , D p Flexural rigidity of cylindrical shell, base, and pipe, respectively (N·mm) E Young's modulus (MPa) E t Tangent modulus (MPa) h Radial clearance between sealing surfaces (mm) h a Height of streaks (mm) h b Depth of embedded streaks (mm) h c Equivalent height of the hypothetical rectangular channel (mm) L 0 Axial distance between radial load and ridge (mm) L 1 Axial distance between anchors and ridge (mm) L c Length of the leakage channel (mm) M ϕ , M x Bending moments of axial section and a section perpendicular to the axis of a cylindrical shell, respectively (N·mm) N ϕ , N x Membrane forces of axial section and a section perpendicular to the axis of a cylindrical shell, respectively (N) Inlet the outlet pressures of the leakage channel (MPa) P b1 , P b2 Components of radial load P z (N/mm) P r Reaction load (N/mm) P x Shearing force of intensity parallel to z axis of a section perpendicular to the axis of a cylindrical shell (N) P z Total equivalent radial load on the base (N/mm) Q c Helical leakage ratio (L/h) r c Corner radius of the tool nose (mm) R, R b , R p Radius of mid-surface of cylinder, base, and pipe, respectively (mm) R c Radius of leakage path (mm) u, v, w Displacements of the mid-surface in the cylindrical coordinate (mm) w b1 , w b2 Deflection components of P z (mm) w x , w p Radial deflection of the base and pipe, respectively (mm) w k Deflection of the streak peaks (mm) w r Deflection component of P r (mm) w s Deflection of pipe surface at the seal (mm) x, y, z Global rectangular coordinates of 2D axisymmetric FE model ϕ, x, ρ(=R+z) Cylindrical coordinates