Thermal–Structural Coupling Analysis of Subsea Connector Sealing Contact

: Taking a subsea collet connector as an example, the contact characteristics of the sealing structure of the subsea connector under thermal–structural coupling were studied. Considering the heat transfer problem of the subsea connector in deep water, the heat transfer model of seawater layer between sealing structures was established


Introduction
The connection and sealing technology between the equipment of subsea production systems in the deep-water oil and gas field is an internationally recognized technical problem [1].The sealing contact performance of the subsea connector determines the success or failure of the connection.In particular, the temperature characteristics of the sealing structure under the action of internal high-temperature oil and gas and external low-temperature seawater are among the fundamental factors for the long-term reliable operation of the connector.The sealing structure is subjected to the combined action of external force and temperature load during underwater work, and it is very easy to produce excessive deformation of the contact surface, which affects the sealing performance [2].Therefore, it is necessary to conduct thermal-structural coupling analysis of the contact surface of the sealing structure of the subsea connector.
Scholars have done a lot of research on the contact under a structural load or temperature load.In 2007, Abid M et al. [3] analyzed the sealing performance of flange connectors under variable steady-state thermal load and thermal transient load and studied the relationship between the contact pressure of the sealing hub, stress relaxation of the bolt, etc., and the temperature load.Zhou Xianjun et al. [4] studied the heat transfer model of a bolted flange connection based on the nonlinearity of gaskets, obtained the distribution of the transient temperature field of the hub, and then analyzed the influence of transient temperature field changes on the bolt load and gasket stress.Al-Turki LI et al. [5] found that the bending strength of composite materials under bending loads mainly depends on the medium and the aging degree, and the decrease in bending strength under contact load is mainly due to the increase in the number of contact cycles.Avanzani A and Donzella G [6] combined the equivalent strain range calculated at the critical point of the seal with the life curve of the material to calculate the fatigue damage and working life of the sealing material in ultra-high pressure applications.In 2009, Omiya Yuya et al. [7] analyzed the thermal stress of bolted flanges under internal pressure and heat conduction conditions and proved that the sealing performance of the bolted flange structure increased with the increase in temperature based on numerical simulation and experiments.In 2015, Sawa T et al. [8] conducted many experiments on the sealing performance of the bolt flange at high temperature, obtained the compression-resilience performance and thermal expansion coefficient of sealing gasket at different temperatures, and carried out experimental verification.In 2016, Luo Yan-Yan et al. [9] studied the influence of heat differences and the temperature alteration ratio on the degeneration of aviation electrical connector performance in thermal fatigue experiments based on the theory of accelerated life experiments and analyzed the failure mechanism of thermal fatigue.In 2018, Wang Lu et al. [10] analyzed the axial bolt force, maximum bolt stress, and gasket contact stress during steady and transient heat processes and found that the maximum bolt stress was proportional to the temperature, while the average gasket contact stress was inversely proportional to the temperature.In addition, the heating rate significantly affected the maximum stress of bolts and gasket contact stress.Abdullah Oday Ibraheem et al. [11] used an axisymmetric model for numerical analysis to simulate the engagement of the single-disc friction clutch system and used a sequential thermo-mechanical coupling method to analyze the thermal stress of the automobile clutch under dry conditions.The research results showed that contact pressure had a significant influence on the thermal stress in both sliding and heating stages.In 2019, Zhang Lanzhu et al. [12] analyzed the stress distribution of metal-to-metal contact flange connections under different pressures, temperatures, and bending moments based on the finite element method and showed that under different thermal loads, when the initial bolt stress was sufficient to reach the metalto-metal contact, the maximum bending moment that the connection could withstand was determined by strength criteria.In 2019, Tang Liping et al. [13] found that the sealing performance of a metal polymer interface had more complex characteristics.Temperature has an important impact on the thermal deformation of the sealing structure in the temperature field, and an appropriate high temperature can improve the sealing ability.In 2020, Chen Jinlin et al. [14] considered the elastic deformation and thermal stress deformation of the micro-convex body based on fractal theory and analyzed the normal contact stiffness of the double friction interface of dry gas seals.These authors obtained the influence of the fractal dimension and characteristic dimension on the thermal and elastic normal contact stiffness, that is, the former had a positive ratio, and the latter was inversely proportional.Wang Qiang et al. [15] established a fluid dynamics and heat transfer model of porous media through the thermal stress module of the ANSYS Workbench, corrected the friction heating model, and calculated the thermal deformation of fingers according to the pressure and temperature results.
As can be seen from the above literature review, previous researchers mainly studied the mechanical properties of the seal structure under thermal-structural coupling and the structural contact stress under external load or thermal stress.However, the contact state and mechanical characteristics of the sealing structure under the action of thermal stress and structural stress were not considered sufficiently.Therefore, this paper analyzes the heat transfer process when the subsea connector is subjected to the action of internal hightemperature oil and external low-temperature seawater and simplifies the heat transfer by equivalent thermal conductivity and the composite heat transfer coefficient.The structural stress of the subsea connector under a mechanical load and thermal stress produced by temperature loads is coupled to analyze the contact characteristics and overall stress change of the subsea connector lenticular sealing structure in the thermal-structural coupling field.Because the contact stress of the sealing contact surface cannot be analyzed by a sticking strain gauge, the correctness of this method was verified by comparing the finite element numerical simulation with the theoretical analysis.

The Structure of the Subsea Connector
The metal lenticular sealing gasket is usually used as the main seal for the subsea collect connector and clamp-type connector, which are commonly used in deep sea environments.This article takes the subsea collect connector as an analysis case, as shown in Figure 1.The core sealing structure of the subsea collect connector is composed of a lenticular sealing gasket, top hub, and bottom hub.The external installation tool drives the press ring to move downward to make the fingers clasp the top and bottom hubs and exert pre-tightening force on the sealing structure to complete the connection and sealing.In the whole structure, the parts directly exposed to the oil and gas medium include the top hub, bottom hub, and lenticular sealing gasket.Therefore, these three parts are greatly affected by the change in the temperature load.The applicable range of underwater temperature of the whole underwater connector system is 3 °C ~ 150 °C, which is also the temperature load change range of the subsea connector thermal structure coupling analysis in this paper.When the sealing structure of the subsea connector is subjected to external load and temperature load, the structural stress and thermal stress will overlap.If the coupling stress is too large, it will have an irreversible impact on the sealing structure, especially on the sealing contact surface.Therefore, the mechanical and sealing contact characteristics of the sealing structure of the subsea connector under the thermal-structural coupling effect should be analyzed in detail.

Research on the Heat Transfer Model of the Subsea Connector
The subsea connector works at a relatively constant temperature.If the temperature of oil and gas inside the connector tends to be constant, the temperature of the components from the inside to the outside of the connector will also be constant, which is regarded as the steady-state temperature distribution.According to the second law of thermodynamics, subsea connectors are generally subjected to the action of internal high-temperature oil and gas and external low-temperature seawater during operation, so there must be heat transfer problems.The basic heat transfer modes of the subsea connector are thermal conduction, thermal convection, and thermal radiation.
As shown in Figure 2, the heat transfer model of the subsea connector was used to analyze the temperature field distribution of the sealing structure.In the finite element analysis of the ANSYS Workbench, the calculation methods of the three heat transfer modes are different.Among them, the thermal conduction is close to linear and is relatively simple; the convective heat transfer coefficient cannot be directly obtained, which increases the difficulty of applying boundary conditions; and the thermal radiation is proportional to the fourth power of the object temperature, so its finite element analysis is highly nonlinear, and its convergence cannot be guaranteed by direct calculation.Therefore, it is necessary to establish the corresponding heat transfer model.When establishing a heat transfer model, structural characteristics and calculation accuracy should be considered at the same time, and heat transfer factors that have little influence should be ignored.The heat transfer mode should be simplified, and equivalent thermal conductivity and composite thermal conductivity should be introduced to facilitate the application of boundary conditions.The temperature of the medium inside the subsea connector: at low temperature, the medium temperature is seawater temperature, and at high temperature, the medium temperature is oil and gas temperature; at the same time, the heat transfer coefficients of the internal and external fluid media can be obtained, so the third boundary condition of heat transfer is satisfied [16,17].

Space between hub and collet
Space outside the press ring

Heat Transfer Model of the Seawater Layer between the Hub and Collet
Figure 3 shows a simplified cross-sectional view of the position relationship between the hub and the collet.The seawater layer between the hub and the collet is very thin, but not negligible, otherwise, the temperature difference between the hub and the collet will increase.Since the hub outer surface and the collet inner surface are not flat, the introduced shape factor represents the energy flux between the two surfaces.Firstly, the energy exchange between the area elements of the two surfaces is obtained, and then the integral of the two surfaces is determined.Therefore, the shape factor is mainly related to the geometric state of the two surfaces.In addition, the shape factor only represents the percentage of radiative transfer on one surface and has nothing to do with the energy absorption capacity of the other surface.The heat transfer model between the hub and the collet can be simplified as the heat transfer problem between two concentric cylinders.As shown in Figure 3, L1 is the height of the seawater layer.ε is the emissivity, 0 ≤ ε ≤ 1; σ is the Stephen-Boltzmann constant; A1 is the area of the smaller surface, set as surface 1; T1 is the temperature of surface 1; T2 is the larger surface (set as the temperature of surface 2); φ1-2 is the shape factor from surface 1 to surface 2. The radiation shape factor [18,19] of the hub and collet surface is where r1 is the outer diameter of the flange; r2 is the inner diameter of the collet; . Considering the influence of radiation in the system on heat transfer, the emissivity of the system can be obtained, εs where ε1 and ε2 are the emissivity of the hub and collet.Since the materials selected for each part of the subsea connector are corrosion-resistant alloy steel and heat-resistant alloy steel, during the long-term operation in deep water, they will be corroded by seawater to produce an oxide layer, so the emissivity of each part is taken as 0.8.

Bottom hub
Top hub finger Heat transfer of the seawater between the hub and collet.
The flow of seawater is very slow because the seawater layer is wrapped in the press ring and the wall tube.Therefore, convective heat transfer can be ignored.Accordingly, only two heat transfer methods, thermal conduction and thermal radiation, are considered between the hub and the collet and can be obtained from the energy conservation on the hub surface where Q1-2 is the total heat transferred from the seawater layer between the hub and the collet; λk2 is the thermal conductivity of seawater at temperature T2 [18]; T1 is the temperature of the hub outer surface (as shown in Figure 3); T2 is the temperature of the collet inner surface (as shown in Figure 3); C0 is the black-body radiation coefficient.In the finite element analysis, compared with convective heat transfer and thermal radiation, the calculation and regulation of thermal conduction are easier.Therefore, the equivalent thermal conductivity is defined to describe the entire process of thermal conduction and thermal radiation.The right side of Equation ( 3) is the total energy transferred under the action of thermal conduction and thermal radiation.Using the equivalent thermal conductivity λe1, the heat transfer can be expressed as Through consideration of Equations ( 3) and (4) comprehensively, λe1, can be expressed as The relationship between the equivalent thermal conductivity λe1 and temperature is shown in Figure 4.

Heat Transfer Model of the Seawater Layer between the Hub and Lenticular Sealing Gasket
Figure 5a shows a schematic diagram of heat transfer between the hub inner surface and the lenticular sealing gasket.The heat transfer between the hub inner surface and the lenticular sealing gasket includes the thermal conduction q1 between the lenticular sealing gasket and seawater, the convective heat transfer q2 of the seawater, and the thermal radiation q3 from the lenticular sealing gasket to the top and bottom hub.The gap is narrow and the flow of seawater is extremely low because of the sealing state between the hub and the lenticular sealing gasket.Convective heat transfer on the surface has less influence than thermal conduction and radiative heat transfer, so it can be ignored.As shown in Figure 5b, the model of heat transfer is simplified as a composite heat transfer combining thermal conduction and thermal radiation.The heat transfer quantity of the seawater layer between the lenticular sealing gasket and the hub, Q3-4, can be expressed as: where λk4 is the thermal conductivity of seawater at temperature T4 [18]; T3 is the temperature of sealing gasket outer surface; T4 is the temperature of the top hub inner surface; r3 is the outer diameter of the sealing gasket; r4 is the inner surface of the hub; C0 is the blackbody radiation coefficient; φ3-4 is the shape factor from the sealing ring outer surface 3 to the flange inner surface 4, and the calculation method is the same as that in Equation ( 1).
The equivalent thermal conductivity of the heat transfer model, λe2, can be expressed as: The relationship between equivalent thermal conductivity λe2 and temperature is shown in Figure 6.

The Heat Transfer Model of the Outer Surface of the Finger and the Press Ring
Since the subsea connector has 12 fingers, the seawater between the fingers and the press ring will flow through the gaps between the fingers.This part of the seawater has strong fluidity, so it can be incorporated into the heat transfer model of the outer surface of the press ring for analysis.The heat transfer model on the outer surface of the press ring is shown in Figure 7.As shown in Figure 7, the outer surface of the press ring is in direct contact with the flowing seawater, and the main methods of heat transfer are convective heat transfer q2 generated by the seawater flowing through the outer surface and radiative heat transfer q3 with the outside world.The convective heat transfer loss, Qp, can be expressed as: where hp is the natural heat transfer coefficient of seawater [20]; Sp is the outer surface area of the press ring; T5 is the outer surface temperature of the press ring (as shown in Figure 7); T6 is the temperature of seawater (as shown in Figure 7).In addition, the heat loss caused by thermal radiation, QR, can be expressed as: where QR is the heat loss caused by thermal radiation from the outer surface of the press ring; C0 is the black-body radiation coefficient; ε is the emissivity of the press ring.Substituting the composite thermal conductivity of convection and thermal radiation surface, he, into Equation ( 9), Equation ( 9) can be written as:

T T h S T T h S T T C εS (10)
By simplifying Equation ( 10), the expression of the composite thermal conductivity, he, can be expressed as: 100 100 The composite thermal conductivity, he, which describes the heat transfer of the outer surface of the finger and the outer surface of the press ring, has nothing to do with the shape of the finger, and the press ring as can be seen in Equation (11).Since the outer surface of the finger is in full contact with seawater, the temperature of the finger is similar to the temperature of the press ring.Therefore, the two can be simulated and analyzed together.After calculation, the relationship between the composite thermal conductivity he and the temperature is shown in Figure 8.

Thermal-Structural Coupling Mathematical Model of the Subsea Connector
The subsea connector will produce structural stress under the action of external load.When subjected to the temperature load, the components of the connector will expand with the increase in temperature and contract with the decrease in temperature.Due to the constraints of each component, thermal stress will be generated, so it is necessary to study the coupling of thermal stress and structural stress [21].The stress on the sealing gasket of the subsea connector is divided into three parts: one is the thermal stress generated by the structure due to the combined action of internal and external temperature; the second is the structural stress caused by the internal oil and gas pressure acting on the lenticular sealing gasket and hub; the third is the structural stress generated by the axial preload in the contact area of the sealing gasket.

Three-Dimensional Stress Caused by the Steady-State Temperature Field
According to the characteristics of the sealing structure of the subsea connector, the sealing structure can be simplified to a thick-walled cylinder when the temperature field achieves a steady-state distribution, and the stress state at any point inside can be determined by three-dimensional stress, which includes the circumferential stress σθ, axial stress σz, and radial stress σr in the cylindrical coordinate system [22].In order to use the generalized Hooke law and the thermal stress function method to calculate thermal stress, the following assumptions need to be made for the structure.Suppose that (1) the material selected for the sealing structure of the connector is isotropic; (2) the connector is only subjected to the internal oil and gas pressure, and there is no other external load; (3) there is steady-state thermal conduction among the components of the connector sealing structure; (4) the length of the top hub and the bottom hub is infinite and the constraints on the hub, such as a jumper, are ignored; (5) the connector sealing structure satisfies the generalized Hooke law and small-deformation theory.
Usually, when the temperature of an object rises, it expands.Assuming that the change in temperature is τ (τ = T-T0:T0 is the initial temperature; T is the final temperature.), the expansion of any segment of the micro-body is not limited, and the strain of the isotropic body under free expansion is [23] where α * is the thermal expansion coefficient.However, each element is not free.There must be mutual restraint between them, so thermal stress must be generated.The strain of each element is caused by the combined action of the change of temperature and stress.
For the plane strain problem, the relationship between strain and stress under the action of temperature can be expressed by Hooke's law [23] as where ν is Poisson's ratio; E is the elastic modulus.For a plane problem, if the coordinates are x and y, which are independent variables and do not consider the volume force, then the balance equation can be written as: Then, the compatibility equation of the strain can be expressed as: In the plane strain state, it can be obtained by zz 0 ε  in Equation ( 13) By pluggin Equation ( 16) into Equation ( 13), the relationship between strain and stress in the plane strain state can be expressed as: By introducing Airy's thermal stress function χ, Equation ( 14) can satisfy Because the stress of the connector on the z-axis cannot be ignored, the plane strain model can be selected as its analysis model.Consideration of Equations ( 15), (17), and (18) in the plane strain state: where According to the simplified sealing structure, Δ can be converted into the expression in a cylindrical coordinate system by Equation ( 20) Then, the thermal stress component is: Since the sealing structure is a multi-connected domain, according to the heating condition of the thick-walled hollow cylinder with an inner diameter ri and an outer diameter ro, the boundary conditions are as follows: where a1, a2, a3 are unknown constants.In order to determine the three unknown constants, the Michell expression [23] is added and expressed in polar coordinates.When r = ri: Then, the Michell expression can be described as: The general temperature field distribution, τ, is calculated by Fourier series where Fj(r) and Gj(r) are the Fourier coefficients of τ; j = 0, 1, 2….If τ is the harm onic function in the plane, then Substituting Equation (26) into Equation ( 27), Fj and Gj must satisfy the following relations: Airy's thermal stress function χ is written in the form of Fourier series , cos sin where fj(r) and gj(r) are Fourier coefficients of Airy's thermal stress function.Equation ( 29) must satisfy Equations ( 19)- (25).Substituting fj(r) of Equation ( 29) into Equation (19), Equation ( 19) can be expressed as Substituting fj(r) of Equation ( 29) into Equation (23), Equation ( 23) can be expressed as: By comparing the corresponding coefficients, the following equation can be obtained: ( 2) , , Therefore, for Fj in Equation ( 26), only the temperature distribution of j = 0 and j = 1 affects the stress function.Therefore, when j = 1, applying the Michell expression to Fj and fj, we obtain: When j = 0, it can be obtained from Equation ( 25): (1) For gj, when j ≠ 1, gj = 0.When j = 1, the basic relations are: According to Michell's expression, we can obtain: To sum up, the thermal stress is calculated for the following two cases: First of all, for j = 0, the temperature distribution is only a function of r, so it will not cause shear stress in the micro-body.The basic relations of the stress function can be obtained from Equation (30): The solution of this differential equation is: Since τ (0) in Equation (37) must satisfy Equation ( 28), we can obtain: where K0, K1 are integral constants.Transformation of the above equation results in the following: where .
By substituting Equation (40) into Equation ( 22), the stress component of temperature variation along the radial direction is: According to the two boundary conditions where r = ri and r = ro, σ0 rr = 0 and Equation (34), the above constants are determined as: whereupon: If it is set at r = ri,  i τ τ and  0 τ at r = ro, we can obtain from Equation (43): Substituting Equation (47) into Equation (46), we can obtain: According to the axial force formula in cylindrical coordinate system [23], we can obtain: The subsea connector belongs to steady-state axisymmetric heat transfer [23], and its equation is: When r = ri,  i τ τ ; when r = ro,  o τ τ .Therefore, the solution of this equation is: Substituting Equation (51) into Equation (49), we can obtain: Substituting Equation (52) into Equation (49), we can obtain: The three-dimensional thermal stress

Three-Dimensional Stress Caused by Internal Pressure Load
In the working state, the subsea connector is subjected to multiple loads at the same time.Since the sealing structure is regarded as a thick-walled cylinder, the three-dimensional stress caused by the internal pressure load can be obtained from the Lame formula [24]: The main sealing structure studied in this paper is the metal lenticular sealing gasket.The initial contact between the spherical surface of the lenticular sealing gasket and the conical surface of the hub is a line contact.After the connector is clamped, the contact line is further compressed and deformed into a contact belt [25,26].However, because the contact position of the spherical surface and the conical surface is certain, the width of the contact belt is much smaller than the size of the entire sealing surface, and the width of the contact belt is also called the sealing width.In order to correctly describe the sealing performance of the lenticular sealing structure, the sealing width, contact stress mean value, and the contact stress distribution are essential key parameters.Therefore, the sealing mechanics of the lenticular sealing structure should be studied.
The contact surface between the spherical surface of the lenticular sealing gasket and the conical surface of the hub is expanded around the axis, which can be regarded as the contact between a cylinder and a plane, and it is expressed as a non-conforming surface [27] contacting the O0 point.As shown in Figure 9, the semi-sealing width is a, the radius of the sealing sphere is r, and the elastic modulus and Poisson's ratio of the hub and sealing gasket are E1, E2, μ1 and μ2 respectively.Therefore, the equivalent elastic modulus, E * , can be expressed as: According to Hertz contact [28], for this type of contact stress, it can be expressed as: where P is the line load, * . Then, the maximum contact stress, qmax, is: The mean value of contact stress, q , is: Substituting Equation (59) into Equation (58), it can be seen that: Meanwhile, the semi-sealing width, a, is: In the sealing width, any segment of the micro-unit is marked as ds.For the distributed force of normal contact subjected to q(s), the stress at any point (x0,z0) inside the sealing gasket is: For the contact stress distribution shown in Equation (57), on the contact interface, σx0 = σz0 = q(x0); outside the contact area, the stress components on the surfaces of the two objects are all zero.Substituting

( ( )
) q P a s πa s   into Equation (62) and integrating along the z0-axis, the principal stress inside the sealing gasket in the contact area can be obtained: ( ) McEwen uses the variables m0 and n0 to express the stress at the general point [28], m0 and n0 can be expressed as: where m0 and x0 have the same sign; n0 and z0 have the same sign.Thus, it can be obtained that: The contact stress of any point on the contact surface of lenticular sealing gasket and the structural stress state of any point in the contact area can be obtained in the O0-x0y0z0 coordinate system through this equation.

Conversion of Contact Stress
For the three-dimensional stress obtained under the action of contact stress, the main solution is the three-dimensional stress in the contact area.The structural stress of the lenticular sealing gasket is transformed into the coordinate system, and the structural stress and thermal stress caused by the internal pressure are unified under the same coordinate system.As shown in Figure 10a, the same as in Section 4.3.1, the Cartesian coordinate system O0-x0y0z0 is established with any point O0 on the contact circle as the origin.In addition, the Cartesian coordinate system O-xyz is established with the center of the contact circle as the origin O and the axis of the sealing gasket as z-axis;d 0 OO  is the positive direction of the x-axis.Under the assumption of two-dimensional stress, the micro-body shown in Figure 10b is taken at point O0, and the stress components σx0, σz0 and τx0z0 have been calculated in Section 4.3.1 of Chapter 4. Figure 10c is the orthographic projection of the micro-body.Figure 10d is the orthographic projection obtained by intercepting the micro-body with a plane parallel to Oxy, and the balance equation can be obtained: (66) can be written as: Figure 10e is the orthographic projection obtained by intercepting the micro-body with a plane parallel to Oyz.The same can be obtained: Substituting Equation (65) into Equation (67) and Equation ( 68), we can obtain: ) 0 sin( ) Then, the transformation relationship is: where Considering Equations ( 70)-(72) comprehensively, we can obtain: The cylindrical coordinate system O-rθz, as shown in Figure 11b, is established at the origin of the coordinate system O-xyz, where the r-axis coincides with the x-axis.Then, according to Equation (73), it can be obtained that: In the Cartesian coordinate system O-xyz and the cylindrical coordinate system Orθz, the O0 point has the same plane stress state, namely: From Equations ( 69) and ( 74)-(76), it can be obtained:

Thermal-Structural Coupling Stress
In order to obtain the stress caused by the combined action of multiple loads, it is necessary to couple various stresses.There are two conditions for stress superposition: one is that the system can be described as a linear second-order differential equation; the other is that the effects of various factors on the system cannot cause nonlinear phenomena [22].The structural stress and thermal stress of the lenticular sealing gasket are both within the elastic range and can be expressed by linear second-order differential equations.With the condition for using the superposition principle, the coupling stress caused by various loads, i.e., radial stress σr, circumferential stress σθ, and axial stress σz, can be obtained by applying the superposition principle:

Numerical Simulation of Thermal-Structural Coupling of the Subsea Connector
The simulation analysis of the coupling stress of the subsea connector requires that loads such as oil and gas temperature, pre-tightening force, and oil pressure be applied to the connector, which increases the complexity of the numerical simulation.The simulation analysis model of the subsea connector is limited by constraints, with a degree of freedom of 1 and set as a 2D axisymmetric structure.The plane has 912 elements and 3165 nodes, of which the element is mainly "quad 8".This article first carried out a numerical simulation of the temperature field of the subsea connector and then the node temperature of each component was used as a static load and was applied to the connector with the rest of the load.The thermal-structural coupling stress state of the connector under the action of the steady-state temperature field was obtained by statics.In the simulation, F22 material was used for hubs and fingers, and Incoloy 825 nickel-chromium iron corrosion resistant alloy material was selected for the sealing gasket.

Steady-State Temperature Field Analysis of the Subsea Connector
Under rated working conditions, the subsea connector has high-temperature oil and gas flowing inside and is exposed to low-temperature seawater outside.Set the surface temperature of hub and lenticular sealing gasket in contact with oil and gas to be 150 ℃ and seawater temperature to be 3 ℃.For the surface exposed to the seawater environment, the equivalent thermal conductivity and composite heat transfer coefficient obtained from Section 3 are used according to the third boundary condition.

The Overall Temperature Field Distribution of the Connector Sealing Structure
Figure 12 shows the overall steady-state temperature field of the subsea connector sealing structure.It can be seen that the overall temperature distribution of the connector gradually decreases from the inside to the outside.The temperature of the top and bottom hub and the inner wall of the lenticular sealing gasket is relatively high, which is above 143.39 ℃.The maximum temperature of the system is 150 ℃, which is on the innermost side of the lenticular sealing gasket and hub.The temperatures of the finger and the press ring are low, the temperature of the main body of the finger is below 61.506 ℃, the temperature of the main body of the press ring is below 14.566 ℃, and the lowest temperature of the system is 3.0697 ℃ at the end of the finger.
The lenticular sealing gasket is in direct contact with high-temperature oil and gas and is located inside of the connector leading to low heat emission efficiency and concentrated heat, so the temperature is the highest.Although the hub is in direct contact with high-temperature oil and gas, the outer wall of the hub is directly exposed to low-temperature seawater, so the heat emission efficiency of the hub is higher than that of the lenticular sealing gasket.Most of the surroundings of the finger are in low-temperature seawater.Since the finger temperature is relatively high, the tail is far from the main body and is not in contact with the hub and lenticular sealing gasket of higher temperature, so the temperature is the lowest.In order to show the temperature gradient inside the connector more intuitively, as shown in Figure 13a, a path from the midpoint Pi of the lenticular sealing gasket to the point Po on the outer surface of the press ring is established.The temperature change of each point on the path is shown in Figure 13b.Among the components of the subsea connector, the temperature of the lenticular sealing gasket drops by 22.32 ℃, and the temperature gradient is −0.79 ℃/mm; the temperature of the hub drops by 60.01 ℃, and the temperature gradient is −1.05 ℃/mm; the temperature of the finger drops by 19.54 °C, and the temperature gradient is −0.38°C/mm; the temperature of the press ring drops by 2.64 °C, and the temperature gradient is −0.03 °C/mm.In the seawater layer, the temperature of the seawater layer 1 between the sealing gasket and the hub drops by 9.32 ℃, and the temperature gradient is −1.33 ℃/mm; the temperature of the seawater layer 2 between the flange and the finger drops by 4.63 ℃, and the temperature gradient is −1.16 ℃/mm; the temperature of the seawater layer 3 between the finger and the press ring drops by 5.07 ℃, and the temperature gradient is −0.60 ℃/mm.It can be seen that on the path PiPo, the falling temperature and temperature gradient of the hub are the maximum values of all parts, and the values of the sealing gasket, finger, and press ring decrease in order, which is consistent with the previous analysis.

Temperature Field Distribution of the Lenticular Sealing Gasket
Since the lenticular sealing gasket is the easiest to use to conduct a comparative analysis between theory and simulation in the subsequent thermal-structural coupling process, its steady-state temperature field is discussed separately here.The temperature distribution of the lenticular sealing gasket is shown in Figure 14a.The temperature is uniformly distributed along the circumferential direction and gradually decreases from the inside to the outside and the temperature gradient is small.The highest temperature is 150 ℃, which is generated on the inner wall and at the end of the lenticular sealing gasket; the lowest temperature is 126.46 ℃, which is generated in contact with the seawater layer.As shown in Figure 14b, a path is established from the upper end face to the lower end face.The temperature on the extraction path is shown in Figure 14c, which clearly shows that the temperatures of the upper and lower parts of the lenticular sealing gasket are approximately symmetrically distributed.On this path, from 0~27.5 mm, the temperature continuously decreases from 150 ℃ to 142.68 ℃; from 27.5~65 mm, the temperature slowly rises to 145.32 ℃; from 65~92.5 mm, the temperature slowly drops again to 144.22 ℃; from 92.5~118 mm, the temperature quickly rises to 150 ℃.The overall temperature of the lenticular sealing gasket exceeds 126 ℃, so there is no need to worry about large internal stress changes due to an excessive temperature gradient.At the same time, Incoloy 825 has good high-temperature mechanical properties, and the temperature between 126~150 ℃ has little effect on its mechanical properties.It can be seen that the lenticular sealing gasket can maintain good sealing performance under the steady-state temperature distribution.

Analysis of Coupling Stress Examples under a Steady-State Temperature Field
The stress analysis of the lenticular sealing gasket of the subsea connector under the combined action of preload, internal pressure, and temperature load is carried out by using the composite stress expression of steady-state thermal structural coupling analyzed in Section 4.3.3,and the results are compared with the finite element simulation results.
The coupling stress analysis model of the lenticular sealing gasket is shown in Figure 15a, and this analysis was carried out by using the core sealing structure of the 1/12 subsea connector.Table 1 shows the material parameters of each part.The inner side of the structure is subjected to an internal pressure of 34.5 MPa.Frictionless support is used on both sides of the radial direction.The bottom of the bottom hub is set as the fixed end, and 1/12 of the axial pre-tightening force [29], which is 856.8 kN, is applied to the top of the top hub.At the same time, the steady-state temperature field analysis results of 150 ℃ inside and 3 ℃ outside the connector are introduced to simulate the thermal-structural coupling stress of the connector, and the sealing characteristics of the lenticular sealing gasket under the thermal-structural coupling action are obtained.The finite element simulation analysis in the previous section shows that the temperature difference between the inside and outside of the lenticular sealing gasket is τ = 10.17 ℃, and the initial maximum contact stress is qmax = 458.03MPa.These and other parameters are recorded in Table 1.As shown in Figure 15b, in order to study the sealing performance of the lenticular sealing gasket under the action of thermal-structural coupling, the contact position of the lenticular sealing gasket on any side can be used for analysis.From the analysis in Section 4.3.1, it can be seen that on the contact interface of the seal gasket, the principal stress is the contact stress, which cannot be superimposed with the structure thermal stress in Section 4.1.In the contact area of the sealing gasket, after the stress coordinate transformation in Section 4.3.2, the stress at any point has been given by Equation (77), which is the same as the calculation method of structure thermal stress, so thermal-structural coupling can be carried out.As shown in Figure 15b, in order to facilitate the finite element simulation analysis, the connecting line between the two endpoints on both sides of the sealing surface is taken as the analysis position, and the distance between the connecting line and the contact surface is 0.27 mm.
Substituting the known parameters into Equation (78), the theoretical calculation results of the three-dimensional stress in the cylindrical coordinate system are obtained, and the comparison with the results obtained by the simulation analysis is shown in Figure 16.The change trend of the three-dimensional stress in the sealing width of the theoretical analysis is basically the same as the simulation result.Since the lenticular sealing gasket is under pressure, the three-dimensional stress is negative.As shown in Figure 16a, the theoretical maximum value of radial stress σr along the negative direction of r-axis is 477.94MPa, and the simulation maximum value is 453.95MPa, both of which occur at the contact center.The theoretical analysis is 5.28% higher than the simulation analysis.As shown in Figure 16b, the maximum value of circumferential stress σθ occurs at the contact center.Along the negative direction of the θ-axis, the theoretical maximum value is 91.17MPa, and the simulation maximum value is 83.78 MPa.The theoretical analysis is 8.82% larger than the simulation analysis.As shown in Figure 16c, the maximum value of the axial stress σz also occurs at the contact center.Along the negative direction of the z-axis, the theoretical maximum value is 344.08 MPa, and the simulation maximum value is 326.71MPa.The theoretical analysis is 5.32% larger than the simulation analysis.At the same time, Figure 16 shows that at z0 = 0.27 mm, the sealing width of theoretical analysis is about 5.52 mm, and the sealing width of simulation analysis is about 6.28 mm.The theoretical analysis is 12.10% lower than the simulation analysis.The absolute values of the stress of σθ and σz on the negative semi-axis of the x-axis are less than those on the positive semi-axis of the x-axis.After analysis, it can be seen that the theoretical coupling stress analysis first performs the two-dimensional stress analysis and then performs the three-dimensional expansion, while the simulation directly performs three-dimensional finite element analysis.There is a certain algorithm error between the two.Since the two contact surfaces are assumed to be smooth in the theoretical analysis, the edge position of the sealing width is greatly affected.At the same time, due to the influence of friction, the stress at the bottom end of the contact position of the sealing gasket is greater than that at the top end, and σθ and σz are the most obvious stresses affected by this influence.

Coupling Stress Analysis in the Pressure Shock Mode
The pressure shock mainly simulates that the connector is subjected to a shock higher than the working pressure in a short time when the temperature is stable.As shown in Figure 17, the shock pressure used is 1.5 times of the rated working pressure (51.75 MPa). Figure 18 shows the change of the maximum contact stress of the lenticular sealing gasket and the maximum equivalent stress of each component.As shown in Table 2, the parameters under this working condition are extracted.After analysis, it can be seen that the change trend of the maximum contact stress of the lenticular sealing gasket with the pressure shock is similar to that under the condition of rapid pressure rise and pressure reduction.Under the steady-state temperature distribution, when the oil gas pressure load reaches 1.5 times of the rated working pressure, the maximum contact stress of the lenticular sealing gasket increases from 458.76 MPa to 492.59 MPa, which meets the contact stress conditions of the oil and gas sealing.Under the impact load of oil and gas pressure, the maximum equivalent stress of the finger is 113.00MPa, and the maximum equivalent stress of the hub is 195.21MPa, both of which are below the theoretical yield limit (310 MPa) of the 12Cr2Mo1 material.At the same time, the change of the maximum equivalent stress of the finger is 1.21 MPa, and the change of the maximum equivalent stress of the hub is 18.57MPa, which shows that the finger is affected by a smaller impact.The maximum equivalent stress at the spherical surface of the lenticular sealing gasket is 275.22 MPa, which is 23.37 MPa higher than the rated working pressure, and this stress further aggravates the elastoplastic and plastic deformation of the sealing spherical surface.It shows that the lenticular sealing gasket is more sensitive to high pressure under high temperature conditions and is prone to produce large plastic deformation.During the service period of the connector, it is necessary to avoid opening and closing the wellhead several times in a short time, so as to prevent the connector from pressure impact and fatigue wear on the sealing surface under high temperature.

Conclusions
In this article, the heat transfer model of the subsea connector is established and solved, and the steady-state temperature distribution of the connector is obtained.Considering the stress and deformation of the subsea connector under the thermal load, the thermal-structural coupling analysis model of steady-state temperature field is established.The theoretical analysis of thermal stress of the key sealing structures of the connector is compared with the numerical simulation, and the pressure shock mode under the steady-state temperature field is also analyzed.It is proved that the coupling mathematical model proposed in this article can be applied to the thermal-structural coupling theoretical analysis of the similar subsea oil and gas equipment.The main conclusions are as follows: (1) Considering the heat transfer problem of the subsea connector in deep water, the equivalent heat transfer models of seawater layer between the lenticular sealing gasket and hubs, between hubs and fingers, and outside the outer surface of the fingers and the press ring are established, and the relationship between equivalent thermal conductivity, composite heat transfer coefficient, and temperature is solved.(2) The mathematical model of steady-state thermal structural coupling of the subsea connector is established and verified by the simulation analysis.The theoretical analysis of the radial stress, circumferential stress, and axial stress is 5.28%, 8.82%, and 5.32% larger than the simulation analysis, and the sealing width is 12.10% smaller.(3) The steady-state temperature distribution of the subsea connector under rated working condition is simulated.The finite element simulation shows that in the center path of the connector's lenticular sealing gasket, the hub's falling temperature is 60.01 ℃ and the temperature gradient is −1.05 ℃/mm, which are the maximum values of all parts.The values of the sealing gasket, finger, and press ring decrease in turn.The temperature field of the lenticular sealing gasket is symmetrically distributed, and the temperature is always in the high temperature stage (126 ~ 150 ℃), which has little influence on the mechanical properties of the material, so as to ensure its stable sealing performance.(4) The numerical simulation of pressure shock mode under steady temperature field shows that the maximum equivalent stress of the sealing gasket exceeds the theoretical yield limit of the material under the combined action of temperature and pressure, resulting in plastic deformation, which indicates that the lenticular sealing gasket is sensitive to high pressure shock under high temperature.

Figure 1 .
Figure 1.The subsea connector.(a) Installation tool and top connector; (b) the structure of the subsea connector.

Figure 2 .
Figure 2. Heat transfer model of the subsea collet connector.
temperature in T 1 and T 2 (K)

Figure 4 .
Figure 4. Relationship between λe1 and the temperature.

Figure 5 .
Figure 5. Heat transfer between the hub and the lenticular sealing gasket.(a) Schematic diagram of heat transfer; (b) Simple diagram of heat transfer.

Figure 7 .
Figure 7. Heat transfer of the collet outer surface and the press ring outer surface.

Figure 8 .
Figure 8. Relationship between the composite thermal conductivity he and the temperature.

Figure 9 .
Figure 9. Contact principle diagram of the seal gasket.

Figure 10 .
Figure 10.Stress state analysis near the contact position.(a) The contact position; (b) The microbody; (c) Orthographic projection of the micro-body; (d) The orthographic projection obtained by intercepting the micro-body with a plane parallel to Oxy; (e) The orthographic projection obtained by intercepting the micro-body with a plane parallel to Oyz.
Figure11a, the coordinate system O0-x0y0z0 is transformed by the rota-

Figure 12 .
Figure 12.Steady-state temperature field of the whole subsea collet connector.

Figure 13 .
Figure 13.Steady-state temperature distribution of the subsea collet connector on the center path.(a) Temperature path PiPo; (b) The temperature change at each point along the path PiPo.

Figure 14 .
Figure 14.Steady-state temperature field of the lens seal gasket.(a) Temperature field nephogram; (b) Path of the temperature field; (c) Temperature distribution along the path.

Figure 15 .
Figure 15.Simulation and calculation model of the sealing structure.(a) Simulation model; (b) Location of stress calculation.

Figure 18 .
Figure 18.Maximum contact stress and equivalent stress changing with the oil and gas pressure in the pressure shock mode.(a) Maximum contact stress; (b) Maximum equivalent stress. )

Table 2 .
Contact stress and equivalent stress of the main sealing structures in every stage in the pressure shock mode.