Numerical Analysis of Mechanical Behaviors of Composite Tensile Armored Flexible Risers in Deep-Sea Oil and Gas

: As oil and natural gas production continue to go deeper into the ocean, the ﬂexible riser, as a connection to the surface of the marine oil and gas channel, will confront greater problems in its practical application. Composite materials are being considered to replace steel in the unbonded ﬂexible pipe in order to successfully meet the lightweight and high-strength criteria of ultra-deep-water oil and gas production. The carbon-ﬁber-reinforced material substitutes the steel of the tensile armor layer with a greater strength-to-weight ratio. However, its performance in deep-water environments is less researched. To investigate the mechanical response of a carbon ﬁber composite ﬂexible riser in the deep sea, this study establishes the ABAQUS quasi-static analysis model to predict the performance of the pipe. Considering the special constitutive relations of composite materials, the tensile stiffness of steel pipe and carbon ﬁber-reinforced composite ﬂexible pipe are predicted. The results show that the replacement of steel strips with carbon ﬁber can provide 85.06% tensile stiffness while reducing the weight by 77.7%. Moreover, carbon-ﬁber-reinforced strips have a lower radial modulus, which may not be sufﬁcient to cause buckling under axial compression, so the instability of the carbon ﬁber composite armor layer under axial compression is further studied in this paper; furthermore, the characteristics of axial stiffness are analyzed, and the effects of the friction coefﬁcient and hydrostatic pressure are discussed


Introduction
The offshore flexible riser is a pipeline that connects the subsea oil and gas wellhead to the gathering station (ship) and is crucial to the development of offshore oil and gas resources [1][2][3].It has gradually replaced traditional steel pipe due to its high corrosion resistance, good terrain adaptability, long continuous length, ease of installation, and other advantages [4].The unbonded flexible pipe is a composite structure formed by a metal winding structure and polymer extrusion; they work together but are independent of each other.Flexible pipes have been frequently utilized in offshore oil and gas production throughout the past few decades.During installation and operation, risers may be subjected to axial compression and bending loads, and complex constructions composed of numerous layers are extremely deformable when subjected to axisymmetric loads.The failure modes of tensile armor wire include fracture, radial, and lateral buckling.The unbonded flexible pipes typically consist of a carcass layer, a pressure armor layer, tensile armor layers, inner and outer sheaths, anti-wear layers, and an anti-buckling layer.During installation and operation, the riser may be subjected to axial compression and bending loads.Fracture, radial buckling, and lateral buckling are the mechanisms of failure for tensile armor layers.Theoretical analysis, numerical simulation, and experimental verification are the three main methods in engineering applications of flexible risers.The theoretical analysis is to obtain the mechanical property parameters of the structure by solving Equations [5][6][7][8].Typical unbonded flexible (Figure 1) pipes have complex structural layers, making it difficult to perform a structural mechanical study of the entire pipe.Therefore, a large number of scholars use numerical simulation to analyze flexible risers [9][10][11][12][13][14][15][16][17][18][19].Because the working environment of the flexible riser is located in the deep sea, it is difficult to conduct field experimental research, so there are only a few experimental studies [20][21][22][23].
ci. Eng.2023, 11, x FOR PEER REVIEW 2 of 20 verification are the three main methods in engineering applications of flexible risers.The theoretical analysis is to obtain the mechanical property parameters of the structure by solving Equations [5][6][7][8].Typical unbonded flexible (Figure 1) pipes have complex structural layers, making it difficult to perform a structural mechanical study of the entire pipe.Therefore, a large number of scholars use numerical simulation to analyze flexible risers [9][10][11][12][13][14][15][16][17][18][19].Because the working environment of the flexible riser is located in the deep sea, it is difficult to conduct field experimental research, so there are only a few experimental studies [20][21][22][23].In recent years, the growth of the deep-sea oil and gas industry has highlighted the importance of considering the weight of flexible pipes in their design.The weight of the pipe has a significant impact on the size of the platform and the magnitude of the axial tensile load.To address this challenge, researchers and industry practitioners have been exploring the potential of new materials, such as fiber-reinforced composites.These materials offer a higher strength-to-mass ratio than traditional metals, meaning that less material is required to achieve the same strength.Additionally, fiber-reinforced composites exhibit a lower modulus in the non-fiber direction, resulting in greater flexibility and enhanced resistance to corrosion and fatigue compared to metals.As such, fiber-reinforced composites hold great promise for improving the performance and efficiency of flexible pipe systems in deep-sea applications [24].The application of carbon fiber-reinforced plastic composites in offshore platforms is growing rapidly, especially in flexible risers and pressure vessels [25,26].
The application of composite materials in flexible risers began in the eighties of the last century.In the 1980s, the French Institute for Petroleum Research (IFP) and the French Institute for Aerospace Research (Aerospatiale) made the initial effort to design and assess composite production risers [27].Belye et al. [28] conducted a comparative study on steel and composite drilling risers and developed a design method.Lotveit and Ward [29] first proposed the use of aramid fibers as tensile sheath lines in unbonded flexible risers but did not consider the poor corrosion resistance of aramid materials.In 1993, Coflexip Stena Offshore manufactured the first flexible pipe with a stretched composite sheath line using a fiberglass-reinforced epoxy composite.In the mid-1990s, under the supervision of the National Institute of Standards and Technology Advanced Technology program, a joint industrial project was started to develop a composite flexible riser [30][31][32] and a composite drilling riser [33][34][35].Kalman et al. [36] proposed a lightweight design that uses composite tape In recent years, the growth of the deep-sea oil and gas industry has highlighted the importance of considering the weight of flexible pipes in their design.The weight of the pipe has a significant impact on the size of the platform and the magnitude of the axial tensile load.To address this challenge, researchers and industry practitioners have been exploring the potential of new materials, such as fiber-reinforced composites.These materials offer a higher strength-to-mass ratio than traditional metals, meaning that less material is required to achieve the same strength.Additionally, fiber-reinforced composites exhibit a lower modulus in the non-fiber direction, resulting in greater flexibility and enhanced resistance to corrosion and fatigue compared to metals.As such, fiber-reinforced composites hold great promise for improving the performance and efficiency of flexible pipe systems in deep-sea applications [24].The application of carbon fiber-reinforced plastic composites in offshore platforms is growing rapidly, especially in flexible risers and pressure vessels [25,26].
The application of composite materials in flexible risers began in the eighties of the last century.In the 1980s, the French Institute for Petroleum Research (IFP) and the French Institute for Aerospace Research (Aerospatiale) made the initial effort to design and assess composite production risers [27].Belye et al. [28] conducted a comparative study on steel and composite drilling risers and developed a design method.Lotveit and Ward [29] first proposed the use of aramid fibers as tensile sheath lines in unbonded flexible risers but did not consider the poor corrosion resistance of aramid materials.In 1993, Coflexip Stena Offshore manufactured the first flexible pipe with a stretched composite sheath line using a fiberglass-reinforced epoxy composite.In the mid-1990s, under the supervision of the National Institute of Standards and Technology Advanced Technology program, a joint industrial project was started to develop a composite flexible riser [30][31][32] and a composite drilling riser [33][34][35].Kalman et al. [36] proposed a lightweight design that uses composite tape instead of steel tensile armored wire.In 2001, composite risers were first used on the Heidrun TLP platform in the Norwegian Sea and their successful deployment accelerated the commercial pace of fiber-reinforced materials for deep-sea applications.The 2002 OTC conference went into the deep ocean and presented many ideas, conceptual designs, and analyses for composite applications in the ocean, indicating that composites will gradually replace traditional steel in deep-sea applications in the future.Ochoa and Salama [37,38] proposed several difficulties in the use of composite materials to replace metal materials, two of which are the in-depth excavation of the damage mechanism of composite materials and the effective application of nondestructive monitoring technology (NDE) on composite materials.In 2003, the Norwegian Register of Shipping gave a design specification for a composite riser [39], and Ochoa gave a design guideline for a composite rigid riser in a report for the US Federal Minerals Management Service (MMS) [40].Ramirez et al. conducted a collapse pressure test on the composite flexible riser for offshore drilling; it was found that the carbon fiber tube can be stratified under pressure but there is only one experimental sample [41].Almeida et al. [42,43] evaluated the deterioration and failure of carbon fiber composite armored pipes and studied the failure of composite pipes under static pressure through numerical and experimental analysis.The damage degradation of the composite riser was evaluated by a genetic algorithm for the first time.Toh et al. [44] analyzed the composite flexible pipe and believed that the weight of the flexible composite pipe could be reduced and the bearing pressure of the offshore platform could be reduced through appropriate design on the premise of maintaining its mechanical properties.The results show that the composite riser has a higher safety margin than the traditional steel pipe under the same sea conditions.Amaechi et al. [25] used the ANSYS Composite PrepPost module to study the influence of composite and put forward the optimization design method.Zhang et al. [45] proposed a cylindrical composite tensile armor layer.The optimum fiber helix angle was obtained by analysis and calculation.
Although several theoretical and experimental studies have been conducted on the application of new composite materials, previous research has primarily focused on the mechanical response of fiber-reinforced risers under axial tensile loads.Most of these studies have only considered the composite cylinder, with little attention paid to the fine tensile armor structure.The link between a riser's global response and its local behavior has not been fully explored and there is almost no research on the buckling performance of composite risers.
Since fiber-reinforced materials are prone to buckling failure due to their small radial stiffness resulting from their orthogonal anisotropic material properties, studying the compression response of composite risers is of great significance.Fiber-reinforced materials are prone to buckling failure because they have small radial stiffness due to their orthogonal anisotropic material properties, so it is of great significance to study the compression response of composite risers.Firstly, this paper summarizes the vertical axial theoretical analysis method and uses ABAQUS software to establish an unbonded flexible riser.Because the experiment of the riser is affected by the actual length and production environment, it is difficult to experiment, and the accuracy of the model is verified by the experimental data.Then, the tensile stiffness and ultimate tensile strength of the carbon fiber riser were analyzed based on ABAQUS.Finally, the mechanical response of composite flexible risers under compression load was considered, and the buckling resistance of the risers under different external pressures was studied.

Theoretical Formulations
Non-bonded flexible risers are constructed with multiple layers that are allowed to slip relative to one another.Because each layer has unique geometric forms and material properties, they play different roles in the stress process, resulting in a complex analysis process.Currently, there are two main methods for analyzing the section characteristics of flexible risers: analytical and numerical methods.These methods can be used to study both axisymmetric and bending responses.The mechanical properties analysis of these two parts are closely related; for example, the interlayer contact pressure generated by the axisymmetric response can have a significant impact on the bending response.Previous investigations have shown that the axisymmetric response of flexible risers is linearly variable, whereas the bending response exhibits clear hysteresis behavior.
This chapter is based on the analysis of some basic assumptions: (1) The flexible riser has large axial and radial stiffness and small deformation.(2) Assume that the friction between each layer is zero, and the layers can slide freely without friction energy consumption.(3) Ignore residual stresses inherited from the manufacturing process.(4) The axial force and torque are uniformly distributed on each layer and the pressure distribution inside and outside the pipeline does not change as the pipe length increases.
According to the different characteristics of the pipeline, the pipeline is divided into three types: (1) Isotropic cylindrical shell-inner sheath, outer sheath, and anti-wear layer.

Theoretical Calculation of Cylindrical Shells
In the case of flexible risers, the outer layer, inner layer, and anti-wear layer, which primarily consist of polymer materials, serve the functions of preventing leakage, reducing friction, and providing protection but contribute very little to bearing external loads.For the purposes of this study, these layers are uniformly classified as polymer layers.Polymer materials can be considered isotropic and linearly elastic within a small deformation range.Moreover, due to the relatively low elastic modulus of polymer materials, the thick-walled cylinder theory is utilized to accurately simulate the radial displacement of these layers.The definition of coordinates and symbols for the bounded flexible pipe as shown in Figure 2.
two parts are closely related; for example, the interlayer contact pressure generated by axisymmetric response can have a significant impact on the bending response.Previo investigations have shown that the axisymmetric response of flexible risers is linearly v iable, whereas the bending response exhibits clear hysteresis behavior.
This chapter is based on the analysis of some basic assumptions: (1) The flexible ri has large axial and radial stiffness and small deformation.(2) Assume that the frict between each layer is zero, and the layers can slide freely without friction energy c sumption.(3) Ignore residual stresses inherited from the manufacturing process.(4) T axial force and torque are uniformly distributed on each layer and the pressure distri tion inside and outside the pipeline does not change as the pipe length increases.
According to the different characteristics of the pipeline, the pipeline is divided i three types: (1) Isotropic cylindrical shell-inner sheath, outer sheath, and anti-wear lay (2) Small-angle spiral winding layer-tensile armor layer.(3) Large-angle spiral laye carcass layer and pressure armor layer.

Theoretical Calculation of Cylindrical Shells
In the case of flexible risers, the outer layer, inner layer, and anti-wear layer, wh primarily consist of polymer materials, serve the functions of preventing leakage, red ing friction, and providing protection but contribute very little to bearing external loa For the purposes of this study, these layers are uniformly classified as polymer laye Polymer materials can be considered isotropic and linearly elastic within a small def mation range.Moreover, due to the relatively low elastic modulus of polymer materi the thick-walled cylinder theory is utilized to accurately simulate the radial displacem of these layers.The definition of coordinates and symbols for the bounded flexible p as shown in Figure 2. The stiffness of the inner sheath, the anti-wear layer, and the outer layer can be c culated by the analysis of the isotropic cylinder.This section gives an expression for isotropic layer.Assuming that the strain is the same along the length of the pipe, the re tionship between strain and displacement can be written as: where R is the radius of the middle plane and L is the length of the pipe.The strain ene can be derived by: The stiffness of the inner sheath, the anti-wear layer, and the outer layer can be calculated by the analysis of the isotropic cylinder.This section gives an expression for the isotropic layer.Assuming that the strain is the same along the length of the pipe, the relationship between strain and displacement can be written as: where R is the radius of the middle plane and L is the length of the pipe.The strain energy can be derived by: The work performed by the external load on the specified displacement increment can be expressed as: where N, T, M x , and M y are the axial force, axial torque, and bending moment for two radial directions, respectively.∆p radial pressure difference in the tube layer, ∆p = p in − p out .

Theoretical Calculation of Large-Angle Spiral Layer
According to the generalized Hooke's law, the stress-strain relation of the thin-walled cylinder in the orthogonal coordinate system can be expressed as: According to the plane stress assumption, τ 12 , τ 13 , τ 23 = 0, so the strain energy can be expressed as: where σ 1 is axial stress, σ 2 is annular stress, and σ 12 is shear stress.

Theoretical Calculation of Tensile Armor Layer
It is assumed that each spiral armor flat steel is linearly elastic and only its axial stress is considered, other stresses are ignored, so the elastic energy of this layer is the sum of all the elastic energy of each spiral armor flat.
Here, n is the number of flat steels in the helical armor layer, A i is the cross-sectional area of its i th stripe, and α is the angle of its spiral winding.

Elastic Constants of Composite Material
At the macroscopic mechanical scale, composite monolayers can be regarded as an orthotropic material.The constitutive relation is expressed in the local coordinate system: where the strength and stiffness of the fiber direction are stronger, called longitudinal, and represented by "1" or "R"; the direction perpendicular to the fiber is called transverse and represented by "2" or "T"; and the thickness direction is represented by "3" or "Z".The same coordinate system is needed to establish the mechanical equilibrium equation for the whole structure.The corresponding coordinate directions of the composite strips are shown in Figure 3.The micromechanics method was used to figure out the properties of composite tensile armor strips.The subscripts F and M denote the fiber and matrix.Based on the prediction formula of the transverse modulus of composite material proposed by Halpin and Tsai [46], the elastic modulus of composite strips in the local coordinate system can be written as follows: where E 1 , E 2 , E 3 are the elastic modulus in all three directions; G 12 , G 13 , G 23 are the shear modulus in all three directions; and µ 12 , µ 13 , µ 23 are Poisson's ratio in all three directions.The subscripts f and m denote the fiber and matrix.The generalized Hooke's law which is the stress-strain relation can be written in the local coordinate system.

Overall Stiffness of a Typical Unbonded Flexible Riser
According to the principle of virtual work, the relationship between work performed by external forces and strain energy can be obtained, as shown in Equation (6).
The equilibrium equation of the whole flexible pipe can be obtained by synthesizing the balance equation of each layer.To be consistent, the same flexible pipe length, axial displacement, axial torsion, and bending were used for each layer when deriving the equilibrium equations for each layer.Moreover, to consider the separation between layers, we allow different radial displacements for each layer.By superimposing the balance equations obtained in the previous sections, we can obtain the global stiffness matrix of the eight-layer flexible riser. [ where [K] is the superimposed overall stiffness matrix; [∆] is the displacement matrix; and [F] is the load matrix.To solve Equation ( 12), using the interlayer pressures ∆p 1 , ∆p 2 , . . . . . .∆p 8 , the interlayer pressure of this equation in the flexible tube can be solved by an iterative method, and the external pressure directly acts on the outermost layer of the flexible tube, namely the outer sheath.The internal pressure acts on the inner sheath.

Compatibility Equation
The radial displacement and deformation of the flexible pipe layer are shown in Figure 4. Assuming that there is no gap between the pipe layers, the pipe should meet the continuum equation of radial displacement.
Where ∆r out and ∆r in are the changes in R in and ∆ t thickness of each layer defined in Figure 4, respectively.Due to the large rigidity of the spiral armor layer (5, 7), pressure armor layer (3), and carcass layer (1), each layer is assumed to have a constant thickness and move along radial stiffness.We can obtain: The other layers (2, 4, 6, 8) can be represented as: an iterative method, and the external pressure directly acts on the outermost layer of the flexible tube, namely the outer sheath.The internal pressure acts on the inner sheath.

Compatibility Equation
The Where Δr out and Δr in are the changes in R and ∆ thickness of each layer defined in Figure 4, respectively.Due to the large rigidity of the spiral armor layer (5, 7), pressure armor layer (3), and carcass layer (1), each layer is assumed to have a constant thickness and move along radial stiffness.We can obtain: The other layers (2, 4, 6, 8) can be represented as: The equilibrium of pressure in the flexible riser can be described by the compatibility equation, which states that the sum of the contact pressures between the layers is equal to the pressure difference between the outer layers of the riser: The pressure difference equation of spiral armor flat steel layer deduced by Lanteigne is: Considering the relationship between the radius and displacement of each layer of flexible pipe, the following equation can be derived: Substituting the deformation coordination relations derived above into the equation, seven deformation coordination equations can be obtained: The equilibrium of pressure in the flexible riser can be described by the compatibility equation, which states that the sum of the contact pressures between the layers is equal to the pressure difference between the outer layers of the riser: The pressure difference equation of spiral armor flat steel layer deduced by Lanteigne is: Considering the relationship between the radius and displacement of each layer of flexible pipe, the following equation can be derived: Substituting the deformation coordination relations derived above into the equation, seven deformation coordination equations can be obtained:

Geometric and Material Parameters
A flexible riser has a complex structure and high cost, and its performance test usually depends on special experimental equipment.It is not practical to use the composite riser since it is situated in a deep-sea environment that is unsuitable for effective experimentation.Therefore, this paper adopts a numerical simulation research method to simulate the composite riser and solve the buckling critical of CFRP risers under the action of axial compression load.In this paper, a 3D finite element model of the flexible riser is established by using ABAQUS, and the length of the pipe is 1000 mm.The geometric and material parameters of each layer are given in Table 1.Some geometric characteristics of each layer refer to Witz's research data [13].
To better simulate the mechanical response of the actual riser under axial compression load during installation and operation, the anti-wear layer is added to the riser studied in this paper to more accurately simulate the actual contact between layers.The internal carcass layer and pressure armor layer are simplified as orthotropic cylindrical shells, the anti-wear layer, inner sheath, and outer sheath are modeled as an isotropic cylinder, and the tensile armor layer is modeled as spiral strip winding.For an accurate simulation of the inter-layer contact stress, the previous research is optimized, and the models of all layers are established as C3D8I elements.The key to simplifying the carcass layer into an orthotropic cylindrical shell is the accurate definition of orthotropic material parameters, here we have to solve for nine parameters: E 1 , E 2 , and E 3 are the elastic modulus of the three axes; G 12 , G 23 , and G 13 are the shear modulus; v 12 , v 13 , and v 23 are Poisson's ratio in three directions.The above parameters in ABAQUS correspond to the engineering constants, which correspond to the defined material properties.The FE model of the flexible pipe adopts carbon-fiber-reinforced strips as orthogonal anisotropic materials defined by engineering constants.The material parameters of the remaining layers are shown in Table 2.The material performance is characterized by nine parameters, including the elastic modulus in three directions, shear modulus in three directions, and Poisson's ratio in three directions, as shown in Table 3.The parameters about T300/5208 graphite-epoxy come from Tawfik's [47] research data.All layers, except for the carcass layer and the pressure armor layer, are isotropic cylindrical shells.In addition, the ABAQUS includes two solvers: implicit solver ABAQUS/Standard and explicit solver ABAQUS/Explicit.The standard implicit solver can be easily applied to solve general problems.However, when solving nonlinear problems with a large number of contact behaviors, the iterative analysis will incur a huge computational cost and the analysis results may not converge.The explicit solver solves nonlinear problems without iteration and the result converges more easily.

Boundary Conditions and Settings
In the finite element model, the normal direction of the interlayer contact type of the riser is selected as hard contact, and the axis is selected as a penalty function.The friction coefficient is selected as 0.1 according to Saevik and Berge [22].ABAQUS/Explicit is used to avoid nonconvergent results at the center of the cross-section at both ends, adopt kinematic coupling at the end of the layer, and apply all loads and boundary conditions to the reference point, as shown in Figure 5.To reduce the instability in the numerical simulation, smooth step loading is used to avoid the fluctuation of the results caused by the discontinuity of loading.Moreover, to ensure the effectiveness of the simulation results, the ratio of kinetic energy (ALLKE) to internal energy (ALLIE) is controlled below 5%.
the analysis results may not converge.The explicit solver solves nonlinear problems out iteration and the result converges more easily.

Boundary Conditions and Settings
In the finite element model, the normal direction of the interlayer contact type riser is selected as hard contact, and the axis is selected as a penalty function.The fr coefficient is selected as 0.1 according to Saevik and Berge [22].ABAQUS/Explicit is to avoid nonconvergent results at the center of the cross-section at both ends, adopt matic coupling at the end of the layer, and apply all loads and boundary conditions reference point, as shown in Figure 5.To reduce the instability in the numerical si tion, smooth step loading is used to avoid the fluctuation of the results caused b discontinuity of loading.Moreover, to ensure the effectiveness of the simulation re the ratio of kinetic energy (ALLKE) to internal energy (ALLIE) is controlled below 5

Axial Tensile Response
In previous research, numerous studies have examined the tensile response of tional steel flexible risers.To verify the proposed FE model, we referred to the experim data of Witz [13] on a 1000 mm flexible riser in the initial stage of the study.Witz's st a seminal work on the experimental investigation of flexible hoses.It provides detail formation on the dimensions and material parameters of each layer of the one-inch fl riser and the test results for its axial symmetric response and bending response.

Model Verification
Table 4 shows some representative results.In this paper, the 3D finite element m is established by ABAQUS, and the numerical simulation results are compared with results and experimental data to verify the accuracy of the numerical model.The t   4 shows some representative results.In this paper, the 3D finite element model is established by ABAQUS, and the numerical simulation results are compared with these results and experimental data to verify the accuracy of the numerical model.The tensile stiffness mentioned in the paper refers to the equivalent tensile stiffness, which can be expressed as: K tensile = F/ε.Under the condition that the end is allowed to twist, the tensile strain curve of the 2.5-inch flexible riser measured experimentally is shown in Figure 6, including three loading periods in total.It can be seen that the curve presents weak nonlinearity with certain hysteresis behavior, and the axial stiffness of the hose (i.e., the slope of the curve) is relatively small when the strain is low, and larger when the strain is higher.The test curve of the subsequent loading cycle is different from that of the first cycle, which is common in the test of the flexible riser and is mainly attributed to the winding process of making the spiral strip.
Figure 6 shows that both theoretical analysis and numerical simulation results exhibit a linear trend.However, the axial stiffness values obtained from these methods are generally higher than those obtained from experimental measurements.Numerical simulations often result in stiffness values that are greater than those obtained in experiments.There are several reasons for this.Firstly, numerical simulations typically involve simplifications and assumptions to make the calculations computationally feasible.These simplifications may not fully capture the complexities and non-linearities of the physical system, leading to an overestimation of stiffness.Secondly, the boundary conditions used in numerical simulations may not perfectly match the experimental conditions, leading to differences in stiffness.Thirdly, the material properties used in numerical simulations may not accurately reflect the properties of the actual materials, leading to differences in stiffness between the simulated and actual systems.Finally, there may be errors or uncertainties in the experimental measurements of stiffness that contribute to the difference between the numerical and experimental stiffness values.The measured quantity in the experiment is relatively small compared to the length of the pipe, and a minor deviation may significantly affect the final measurement outcome.Under the condition that the end is allowed to twist, the tensile strain 2.5-inch flexible riser measured experimentally is shown in Figure 6, includin ing periods in total.It can be seen that the curve presents weak nonlinearity hysteresis behavior, and the axial stiffness of the hose (i.e., the slope of the c tively small when the strain is low, and larger when the strain is higher.The the subsequent loading cycle is different from that of the first cycle, which i the test of the flexible riser and is mainly attributed to the winding process o spiral strip.Figure 6 shows that both theoretical analysis and numerical simulation r a linear trend.However, the axial stiffness values obtained from these meth erally higher than those obtained from experimental measurements.Nume tions often result in stiffness values that are greater than those obtained in There are several reasons for this.Firstly, numerical simulations typically in fications and assumptions to make the calculations computationally feasibl plifications may not fully capture the complexities and non-linearities of the tem, leading to an overestimation of stiffness.Secondly, the boundary cond numerical simulations may not perfectly match the experimental condition differences in stiffness.Thirdly, the material properties used in numerica may not accurately reflect the properties of the actual materials, leading to stiffness between the simulated and actual systems.Finally, there may be err tainties in the experimental measurements of stiffness that contribute to t between the numerical and experimental stiffness values.The measured qu experiment is relatively small compared to the length of the pipe, and a mi may significantly affect the final measurement outcome.Figure 7 illustrates the axial stress distribution of the entire pipe and the tensile armor layers in different flexible pipe models under tensile conditions.The results show that the two tensile armor layers primarily bear the tensile load of the pipeline, while the other layers contribute very little to bearing the load.Additionally, a comparison of the stress distribution in the two tensile armor layers reveals that the axial stress in the inner layer is slightly higher than that in the outer layer.This is primarily due to the fact that the number of armored steel wires in the inner layer is less than that in the outer layer, resulting in higher axial stress in the inner armored layer.
Figure 7 illustrates the axial stress distribution of the entire pipe and the tensile armor layers in different flexible pipe models under tensile conditions.The results show that the two tensile armor layers primarily bear the tensile load of the pipeline, while the other layers contribute very little to bearing the load.Additionally, a comparison of the stress distribution in the two tensile armor layers reveals that the axial stress in the inner layer is slightly higher than that in the outer layer.This is primarily due to the fact that the number of armored steel wires in the inner layer is less than that in the outer layer, resulting in higher axial stress in the inner armored layer.

Carbon Fiber Composite Materials Armored Flexible Pipe
The mechanical properties of composite materials are primarily determined by the fibers used.For comparison with steel armored flexible pipes, the high-modulus carbon fiber composite armored strip was selected.Carbon-fiber-reinforced composite materials are only 20% as dense as conventional steel, making them a highly effective weight-saving alternative with a reduction of up to 80%.To predict the axial stiffness of the carbon-fiberreinforced material armored flexible pipe (CFRP), the numerical model validated in Section 4.1.1 is used in this section.The results of the numerical simulations are compared

Carbon Fiber Composite Materials Armored Flexible Pipe
The mechanical properties of composite materials are primarily determined by the fibers used.For comparison with steel armored flexible pipes, the high-modulus carbon fiber composite armored strip was selected.Carbon-fiber-reinforced composite materials are only 20% as dense as conventional steel, making them a highly effective weight-saving alternative with a reduction of up to 80%.To predict the axial stiffness of the carbonfiber-reinforced material armored flexible pipe (CFRP), the numerical model validated in Section 4.1.1 is used in this section.The results of the numerical simulations are compared with those of the traditional steel armored flexible pipe (SP).
As can be seen from Figure 8, the tensile stiffness of the CFRP riser is very close to that of the steel riser, and the axial stiffness only decreases by 14.94% when the weight loss is close to 80%.Marine flexible risers need strong tensile performance.With the increase in the depth of ocean exploration, the vertical load increases due to the self-weight of the riser, which becomes the main limiting factor.The application of carbon fiber in the deep sea can greatly reduce the axial tensile load of the riser without reducing the axial height of the riser, which increases its application prospects in the deep sea.
fibers used.For comparison with steel armored flexible pipes, the hig fiber composite armored strip was selected.Carbon-fiber-reinforced co are only 20% as dense as conventional steel, making them a highly effec alternative with a reduction of up to 80%.To predict the axial stiffness o reinforced material armored flexible pipe (CFRP), the numerical mode tion 4.1.1 is used in this section.The results of the numerical simulat with those of the traditional steel armored flexible pipe (SP).
As can be seen from Figure 8, the tensile stiffness of the CFRP riser of the steel riser, and the axial stiffness only decreases by 14.94% when close to 80%.Marine flexible risers need strong tensile performance.With depth of ocean exploration, the vertical load increases due to the self-w which becomes the main limiting factor.The application of carbon fiber greatly reduce the axial tensile load of the riser without reducing the axia which increases its application prospects in the deep sea.

Ultimate Tensile Strength
Axial ultimate tensile load is also an important index to evaluate the performance of flexibility vertically.When studying the axial stiffness of the flexible riser, the excessive strain of the flexible riser will cause a self-locking phenomenon of the carcass layer, leading to a significant increase in axial stiffness.The results show that the interlocking of carcasses occurs at about 4% longitudinal elongation [48].
Figure 9 shows the failure behavior of two risers.The numerical simulation results show that at point P1 of SP, the inner steel tensile armor wire enters the plastic deformation stage from the elastic stage, and the axial stiffness of the riser decreases.At point P2, the outer tensile armor wires also enter the plastic stage.The steel armored flexible riser had reached its tensile strength limit at the F1 spot at which point the standpipe axial strain was 1.72% and the riser failed.Compared to SP, the failure of the CFRP showed a completely linear behavior, which was attributed to the high brittleness of fiber reinforcement material, leading to no plastic behavior during deformation.When the axial strain of the riser reached 2.157%, the CFRP flexible riser failed.The CFRP has slightly less axial tensile stiffness than the SP at the same section size.The ultimate tensile strength of CFRP during simulations was 2083 kN, which is 1.29 times the ultimate tensile load of SP (1620 kN).Therefore, CFRP has a better tensile capacity and can work in a higher tension load environment.
1.72% and the riser failed.Compared to SP, the failure of the CFRP showed a completely linear behavior, which was attributed to the high brittleness of fiber reinforcement material, leading to no plastic behavior during deformation.When the axial strain of the riser reached 2.157%, the CFRP flexible riser failed.The CFRP has slightly less axial tensile stiffness than the SP at the same section size.The ultimate tensile strength of CFRP during simulations was 2083 kN, which is 1.29 times the ultimate tensile load of SP (1620 kN).Therefore, CFRP has a better tensile capacity and can work in a higher tension load environment.

Effect of Coefficient Friction
As can be seen from Figure 10, the friction coefficient between layers does not affect the axial stiffness of the riser basically when the axial tensile response is studied.This is because the axial stiffness of the riser is mainly provided by the tensile armor layer, and the other layers rarely bear the axial tensile load.Therefore, the influence of the friction coefficient can be neglected in the simulation of axial tensile response.

Effect of Coefficient Friction
As can be seen from Figure 10, the friction coefficient between layers does not affect the axial stiffness of the riser basically when the axial tensile response is studied.This is because the axial stiffness of the riser is mainly provided by the tensile armor layer, and the other layers rarely bear the axial tensile load.Therefore, the influence of the friction coefficient can be neglected in the simulation of axial tensile response.linear behavior, which was attributed to the high brittleness of fiber reinforcement material, leading to no plastic behavior during deformation.When the axial strain of the riser reached 2.157%, the CFRP flexible riser failed.The CFRP has slightly less axial tensile stiffness than the SP at the same section size.The ultimate tensile strength of CFRP during simulations was 2083 kN, which is 1.29 times the ultimate tensile load of SP (1620 kN).Therefore, CFRP has a better tensile capacity and can work in a higher tension load environment.

Effect of Coefficient Friction
As can be seen from Figure 10, the friction coefficient between layers does not affect the axial stiffness of the riser basically when the axial tensile response is studied.This is because the axial stiffness of the riser is mainly provided by the tensile armor layer, and the other layers rarely bear the axial tensile load.Therefore, the influence of the friction coefficient can be neglected in the simulation of axial tensile response.

Axial Compression Response
For deep-water flexible riser systems, the riser is bearing axial load, bending, and pressure during installation and operation, which may lead to local buckling.As this is a complex mechanical response with inter-layer contact between layers of the flexible riser, previous studies adjust for the buckling of steel tensile armor layers.With the gradual application of carbon-fiber-reinforced materials in the deep sea, using CFRP instead of steel tensile belts has practical applications.It is necessary to study the buckling of carbon fiber tensile armor layers.This part aims to use the finite element simulation software ABAQUS to analyze the flexion resistance of the carbon fiber tensile armor layer, and evaluate the flexion resistance of CFRP by comparing them with SP.
To confirm the accuracy and dependability of the theoretical model described in this work, a typical 4-inch unbonded flexible pipe is used as an illustration.First, de Sousa's experimental data [18] were used to calibrate the FE modeling processes.The geometrical and material parameters are shown in Table 5, and carbon fiber still uses the Table 3 model, T300 carbon fiber.ABAQUS/Explicit is used to avoid nonconvergent results.Create two reference points at the center of the cross-section at both ends, adopt kinematic coupling at the end of the layer, and apply all loads and boundary conditions to the reference point, as shown in Figure 11.There is no unified specification for selecting the contact friction coefficient of flexible pipe, and we set it as 0.1.For normal contact behavior, we choose the default 'hard contact' in the software.coupling at the end of the layer, and apply all loads and boundary conditions to the erence point, as shown in Figure 11.There is no unified specification for selecting the tact friction coefficient of flexible pipe, and we set it as 0.1.For normal contact beha we choose the default 'hard contact' in the software.Radial buckling under axial compression load is actually a local instability effect.manifested in the sudden decrease in axial stiffness, sudden change in the torsion ang the stable armored wire, and sudden increase in radial displacement.A comparison o axial compression-shortening curves from the FEA and the test is shown in Figure 12 good correlation is achieved.Initially, for very low loads, flexible pipes have similar stiff Radial buckling under axial compression load is actually a local instability effect.It is manifested in the sudden decrease in axial stiffness, sudden change in the torsion angle of the stable armored wire, and sudden increase in radial displacement.A comparison of the axial compression-shortening curves from the FEA and the test is shown in Figure 12a.A good correlation is achieved.Initially, for very low loads, flexible pipes have similar stiffness as when they are subjected to axial tension, so the deformation rate produced by the load is very small.As the load increases, the stiffness decreases when the interlayer gap begins to form.When the axial expression load of the FEA is up to 273.04 kN, radial buckling of the tensile armor layer occurs, which corresponds to point A in Figure 12a.
In addition, when the steel tension armor wire is unstable, its torsion direction will change.In this paper, the expression load applied end node is extracted to establish the relationship between the axial compression load and torsion angle per unit length.At point B in Figure 12b, the tensile armor layers buckled, and the corresponding compression load was 273.04 kN.Radial buckling is a local effect, and the reference point at the end selected may have some difference from the actual local measurement date.The axial compression-radial expansion curves from FEA are also consistent with the experimental results, as shown in Figure 12c.
Compared with the experimental data, the finite element model in this paper has the same trend with only a small difference.The error of radial buckling load between FEA and experimental data is only 3.82%, as shown in Table 6.The comparison shows that this model is effective for the prediction of riser buckling.Figure 12 and Table 6 demonstrate that the finite element model aligns well with the experimental test in terms of axial force and displacement response.The model is still considered to be consistent with the experiment.The discrepancies observed between the two can be attributed to variations in the material properties, interlayer friction, and contact interactions of the carcass and pressure armor layers.Furthermore, data uncertainty may also arise from the limited information regarding the test procedures and measurements provided by de Sousa, as well as some simplifications made during the numerical simulation process.Nevertheless, the numerical model's accuracy is still deemed satisfactory for the purposes of this study.
as when they are subjected to axial tension, so the deformation rate produced by the load is very small.As the load increases, the stiffness decreases when the interlayer gap begins to form.When the axial expression load of the FEA is up to 273.04 kN, radial buckling of the tensile armor layer occurs, which corresponds to point A in Figure 12a.In addition, when the steel tension armor wire is unstable, its torsion direction will change.In this paper, the expression load applied end node is extracted to establish the relationship between the axial compression load and torsion angle per unit length.At point B in Figure 12b, the tensile armor layers buckled, and the corresponding compression load was 273.04 kN.Radial buckling is a local effect, and the reference point at the end selected may have some difference from the actual local measurement date.The axial compression-radial expansion curves from FEA are also consistent with the experimental results, as shown in Figure 12c.
Compared with the experimental data, the finite element model in this paper has the same trend with only a small difference.The error of radial buckling load between FEA and experimental data is only 3.82%, as shown in Table 6.The comparison shows that this model is effective for the prediction of riser buckling.Figure 12 and Table 6 demonstrate that the finite element model aligns well with the experimental test in terms of axial force and displacement response.The model is still considered to be consistent with the experiment.The discrepancies observed between the two can be attributed to variations in the material properties, interlayer friction, and contact interactions of the carcass and pressure armor layers.Furthermore, data uncertainty may also arise from the limited information regarding the test procedures and measurements provided by de Sousa, as well as some simplifications made during the numerical simulation process.Nevertheless, the numerical model's accuracy is still deemed satisfactory for the purposes of this study.13a presents the axial shortening of different material armored risers under axial compression force.At the beginning of load application, the steel armored pipe (SP) has greater axial stiffness compared to carbon fiber composite materials armored pipe (CFRP) because steel has a larger modulus.As the load gradually increases, the transverse movement of steel tensile armored strips occurs because the interlayer friction is insufficient to support the load.Carbon fiber materials armored strips have a low lateral modulus, and the armor wires move laterally when the load is applied, so the stiffness of CFRP is linear, showing a different trend compared with SP.SP and CFRP begin to buckle radially at points A1 and A2.The axial strain increases rapidly after buckling, and the buckling loads are 273.04kN and 237.26 kN, respectively.Compared with CFRP, SP has better buckling resistance.Figure 13b presents the change in twist angle per unit length at the end of the under compressive load.The twist angle of SP changes linearly before radial buck while the twist angle of CFRP changes nonlinearly due to the anisotropic material pa eters.Radial buckling of SP and CFRP occurs at points B1 and B2, respectively.As id fied in Figure 13b, at the early stage of loading, SP and CFRP have similar r Figure 13b presents the change in twist angle per unit length at the end of the pipe under compressive load.The twist angle of SP changes linearly before radial buckling, while the twist angle of CFRP changes nonlinearly due to the anisotropic material parameters.Radial buckling of SP and CFRP occurs at points B1 and B2, respectively.As identified in Figure 13b, at the early stage of loading, SP and CFRP have similar radial displacements.With the increase in load, the radial stiffness of SP is higher than that of CFRP, and the increase rate in the radial displacement of SP is smaller than that in CFRP.CFRP is subjected to radial buckling failure under smaller compression loads, which correspond to points C1 and C2 in Figure 13c.

Effect of Hydrostatic Pressure
Flexible risers operate in deep water, where the outer side of the pipe is subjected to the hydrostatic pressure of seawater.The effect of hydrostatic pressure on the outer sheath prevents the radial movement of the flexible riser and increases the contact pressure between layers to increase the friction between layers.In this section, the flexible riser is subjected to a combined load of external pressure and axial compression to evaluate the anti-buckling ability of the riser under external hydrostatic pressure.
Figure 14 presents the axial compression response under different external hydrostatic pressures.The hydrostatic force is applied to the outer sheath, which restrains the radial deformation of the flexible risers, CFRP and SP.Therefore, the axial stiffness of CFRP and SP also increases with the increase in external pressure.When the external pressure increases from 0 MPa to 2 MPa, the axial stiffness of the pipeline increases significantly at the initial stage of load loading, and CFRP and SP have a larger buckling load at 2 MPa.However, as the external pressure increases further, the stiffness of the flexible riser changes very little, with only a small increase.With the change in external pressure, CFRP and SP have the same trend.CFRP has a small axial stiffness, and its stiffness increases more obviously under hydrostatic pressure.
Figure 14 presents the axial compression response under different external hydrostatic pressures.The hydrostatic force is applied to the outer sheath, which restrains the radial deformation of the flexible risers, CFRP and SP.Therefore, the axial stiffness of CFRP and SP also increases with the increase in external pressure.When the external pressure increases from 0 MPa to 2 MPa, the axial stiffness of the pipeline increases significantly at the initial stage of load loading, and CFRP and SP have a larger buckling load at 2 MPa.However, as the external pressure increases further, the stiffness of the flexible riser changes very little, with only a small increase.With the change in external pressure, CFRP and SP have the same trend.CFRP has a small axial stiffness, and its stiffness increases more obviously under hydrostatic pressure.

Effect of Coefficient Friction
In the previous study, it was found that the interlayer friction coefficient had no effect on the axial tensile stiffness of the flexible riser.However, in this section, we investigate the effect of the interlayer friction coefficient on the compression response of the riser.As shown in Figure 15, as the interlayer friction coefficient increases, the transverse deformation of the steel tensile armor wires is inhibited, and the occurrence of radial buckling nodes is delayed, which results in an increase in the critical buckling load of the riser.The critical buckling loads of both the CFRP and SP show a linear increase.

Effect of Coefficient Friction
In the previous study, it was found that the interlayer friction coefficient had no effect on the axial tensile stiffness of the flexible riser.However, in this section, we investigate the effect of the interlayer friction coefficient on the compression response of the riser.As shown in Figure 15, as the interlayer friction coefficient increases, the transverse deformation of the steel tensile armor wires is inhibited, and the occurrence of radial buckling nodes is delayed, which results in an increase in the critical buckling load of the riser.The critical buckling loads of both the CFRP and SP show a linear increase.

Conclusions
This study aims to investigate the mechanical behavior of unbonded flexible risers with composite armor layers under axial tension and compression.To simulate the instability of unbonded flexible pipe armor wires under axial compression, a non-linear finite element model is presented in this work.Furthermore, the effect of external hydrostatic pressure on the buckling mechanism of the flexible riser in a deep-water environment is also explored.The benefits of replacing the steel tensile armor in unbonded flexible pipes with carbon fiber composite materials are discussed.The following conclusions are finally obtained in this paper: (1) The model proposed in this article demonstrates its effectiveness in predicting the tensile response of the riser.When utilizing carbon fiber composite material to replace steel in the tensile armor of the flexible riser, the tensile stiffness of the pipe experiences a decrease of only 15% while maintaining the same section size.Moreover, the ultimate tensile stiffness of the composite flexible pipe is significantly increased.The substitution of carbon fiber material results in a 77.7% reduction in the weight of the tensile armor layer; (2) The axial tensile rigidity of the flexible riser remains unaffected by the coefficient of friction between layers.However, the coefficient of friction significantly impacts the com-

Conclusions
This study aims to investigate the mechanical behavior of unbonded flexible risers with composite armor layers under axial tension and compression.To simulate the instability of unbonded flexible pipe armor wires under axial compression, a non-linear finite element model is presented in this work.Furthermore, the effect of external hydrostatic pressure on the buckling mechanism of the flexible riser in a deep-water environment is also explored.The benefits of replacing the steel tensile armor in unbonded flexible pipes with carbon fiber composite materials are discussed.The following conclusions are finally obtained in this paper: (1) The model proposed in this article demonstrates its effectiveness in predicting the tensile response of the riser.When utilizing carbon fiber composite material to replace steel in the tensile armor of the flexible riser, the tensile stiffness of the pipe experiences a decrease of only 15% while maintaining the same section size.Moreover, the ultimate tensile stiffness of the composite flexible pipe is significantly increased.The substitution of carbon fiber material results in a 77.7% reduction in the weight of the tensile armor layer; (2) The axial tensile rigidity of the flexible riser remains unaffected by the coefficient of friction between layers.However, the coefficient of friction significantly impacts the compression response of the pipe.An increase in the friction coefficient delays the lateral movement of the tensile armor wires, leading to a gradual increase in the critical buckling load of the flexible riser; (3) The finite element analysis (FEA) conducted in this study revealed that the failure of the flexible pipe is caused by excessive radial expansion of the tensile armor wires.This type of failure is identified by the rupture of both the high-strength tape and outer sheath layers of the pipe.The consequence is an abrupt reduction in axial stiffness and a rapid increase in radial expansion; (4) The radial expansion of the tensile armor is constrained as the hydrostatic pressure increases.Due to the lower radial stiffness of the carbon fiber tensile armor layer, the axial stiffness of the composite flexible riser increases more significantly under external pressure.The critical load of the carbon fiber armor layer can reach up to 86.8% of that of the steel armor layer.
The utilization of composite materials in deep-sea flexible risers has revealed promising prospects for composite flexible risers.The further objectives of this study involve conducting experimental studies on composite flexible risers and considering the progressive damage characteristics of composite materials.

Figure 1 .
Figure 1.The structure of typical unbonded flexible riser.

Figure 1 .
Figure 1.The structure of typical unbonded flexible riser.

Figure 2 .
Figure 2. The definition of coordinates and symbols for the bounded flexible pipe: (a) global; local.

Figure 2 .
Figure 2. The definition of coordinates and symbols for the bounded flexible pipe: (a) global; (b) local.

Figure 3 .
Figure 3. Coordinate system of an orthotropic composite strip.The micromechanics method was used to figure out the properties of composite tensile armor strips.The subscripts F and M denote the fiber and matrix.Based on the prediction formula of the transverse modulus of composite material proposed byHalpin and

Figure 3 .
Figure 3. Coordinate system of an orthotropic composite strip.

Figure 4 .
Figure 4. Radial displacement and deformation of pipe layer: (a) vertical view; (b) profile view.

Figure 4 .
Figure 4. Radial displacement and deformation of pipe layer: (a) vertical view; (b) profile view.

Figure 5 .
Figure 5. Load application and boundary conditions.

Figure 5 .
Figure 5. Load application and boundary conditions.

4. Case Study 4 . 1 .
Axial Tensile ResponseIn previous research, numerous studies have examined the tensile response of traditional steel flexible risers.To verify the proposed FE model, we referred to the experimental data of Witz[13] on a 1000 mm flexible riser in the initial stage of the study.Witz's study is a seminal work on the experimental investigation of flexible hoses.It provides detailed information on the dimensions and material parameters of each layer of the one-inch flexible riser and the test results for its axial symmetric response and bending response.4.1.1.Model Verification Table

Figure 7 .
Figure 7. Axial stress nephogram of finite element model: (a) the overall of the riser; (b) tensile armor layers.

Figure 7 .
Figure 7. Axial stress nephogram of finite element model: (a) the overall of the riser; (b) tensile armor layers.

Figure 8 .
Figure 8. Tensile response of CFRP and SP.

Figure 8 .
Figure 8. Tensile response of CFRP and SP.

Figure 9 .
Figure 9. Tensile failure response of CFRP and SP.

Figure 10 .
Figure 10.Axial tension response with different friction coefficients.

Figure 9 .
Figure 9. Tensile failure response of CFRP and SP.

Figure 9 .
Figure 9. Tensile failure response of CFRP and SP.

Figure 10 .
Figure 10.Axial tension response with different friction coefficients.Figure 10.Axial tension response with different friction coefficients.

Figure 10 .
Figure 10.Axial tension response with different friction coefficients.Figure 10.Axial tension response with different friction coefficients.

Figure 12 .
Figure 12.Model verification: (a) axial compression-shortening curves of the flexible riser; (b) axial compression-twist curves of the flexible riser; (c) axial compression-radial expansion curves of the flexible riser.

Figure 12 .
Figure 12.Model verification: (a) axial compression-shortening curves of the flexible riser; (b) axial compression-twist curves of the flexible riser; (c) axial compression-radial expansion curves of the flexible riser.

Figure 13 .
Figure 13.Axial compression response: (a) axial compression-shortening curves of the CFR SP; (b) axial compression-twist curves of the CFRP and SP; (c) axial compression-radial expa curves of the CFRP and SP.

Figure 13 .
Figure 13.Axial compression response: (a) axial compression-shortening curves of the CFRP and SP; (b) axial compression-twist curves of the CFRP and SP; (c) axial compression-radial expansion curves of the CFRP and SP.

Table 1 .
Parameters of the unbonded flexible riser with a composite reinforced layer.

Table 2 .
Material properties of traditional unbonded flexible risers.

Table 4 .
Axial tensile stiffness obtained by different scholars.

Table 6 .
Difference between finite element model and experiment.