Next Article in Journal
Smart-Grid-Aware Load Regulation of Multiple Datacenters towards the Variable Generation of Renewable Energy
Previous Article in Journal
The Effects of Textural Parameters of Zeolite and Silica Materials on the Protective and Functional Properties of Polymeric Nonwoven Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Study of Vortex-Induced Vibration of FRP Composite Risers with Large Length to Diameter Ratio Under Different Environmental Situations

1
School of Civil and Architectural Engineering, Shandong University of Technology, Zibo 255000, China
2
Department of Mechanical Engineering and Mechanics, Lehigh University, Bethlehem, PA 18015, USA
3
School of Civil and Environmental Engineering, University of Illinois Urbana-Champaign, Champaign, IL 61801-2352, USA
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2019, 9(3), 517; https://doi.org/10.3390/app9030517
Submission received: 12 December 2018 / Revised: 28 January 2019 / Accepted: 30 January 2019 / Published: 3 February 2019
(This article belongs to the Section Acoustics and Vibrations)

Abstract

:
The increasing attraction of crude resources on the ocean floor leads to the development of risers, which are an indispensable component of an offshore oil and gas exploitation structure. High-grade steel is normally utilized for risers, however, fiber reinforced polymer (FRP) composite risers have become a feasible and economical alternative due to the desirable properties of the material. A comparative study of vortex-induced vibration (VIV) of FRP composite risers with two designed geometries were conducted and an X80 steel riser was used as the benchmark. The length to diameter ratios of all three risers were set to 100 and their VIV responses under different environmental situations, including the natural frequencies, global displacements, global stresses and the stress distributions in every composite layers, were obtained using computational fluid dynamics (CFD) simulation with coupled fluid–structure interaction. The VIV characteristics of both FRP composite risers and their distinction compared with that of steel risers were analyzed and discussed. From this paper, the dynamic characteristics in different environmental situations of VIV for FRP composite risers can be further understood.

1. Introduction

The consumption of oil and gas is exponentially increasing throughout the world. The limitation of crude resources of oil and gas on land makes the resources on the ocean floor more and more attractive. In order to achieve the offshore oil and natural gas extraction, many different types of production systems including conventional fixed platforms, tension-leg platforms, compliant towers, semi-submersibles and floating production facilities, have been employed. While these constructions vary mainly in the ways in which their platforms are positioned over sub-sea wells, the risers, both drilling risers which are used to guide a drilling stem and conduct the drilling fluid upwards, and production risers which raise the extracted oil or natural gas to a floating platform [1,2], are indispensable.
In offshore engineering, the material usually employed for structural components is high-grade steel. Nowadays, the utilization of fiber reinforced polymer (FRP) composites is gradually increasing due to their superior mechanical properties, low density and other desirable properties, such as better thermal insulation, and corrosion and fatigue resistance [3,4]. Therefore, it has been widely recognized that using FRP composites instead of steel for risers will lead to the significant reduction of operational costs and maintenance costs of the whole offshore production platform systems [5,6,7].
In order to apply FRP composites in offshore risers successfully, their durability in sea water is crucial. Venkatesan et al. [6] found that none of the long-term properties of carbon fiber reinforced polymer (CFRP) composites experienced any significant variations after exposure to pure water and seawater at different temperatures. On the other hand, many researchers [3,4,8,9] found that the long-term tensile strength of a CFRP reduced to about 80–95% of its short-term values. For thermoplastic composite carbon/polyether ether ketone (PEEK), the transverse strength decreased obviously in boiling water [10] and therefore, the maximum service temperature of thermoplastic FRP composites has to be well below the glass transition temperatures of their polymer matrices. Moreover, in order to obtain reliable results using a numerical model, constitutive laws and damage models for FRP components have to be selected to optimize the accuracy of the numerical predictions [11].
Besides the research of mechanical properties of FRP composite materials in seawater, the responses and performances, including load distributions, fatigue, and vortex-induced vibrations (VIVs) of composite risers under environmental and functional loads have also been investigated [12,13,14,15,16,17,18,19,20,21,22,23,24,25]. The improvement of composite risers in performance over a pure metallic mandrel had been proved through a scaled-down test [23]. For the global response, the tension force decreases with an increasing water depth and the maximum bending moment occurs on the bottom, followed by the joints at the sea surface in general [12,13,14,16,20,22]. The axial tension and the bending moment throughout an entire composite riser are smaller due to its lower overall weight compared with steel risers [14,16]. The excellent fatigue resistance of FRP composites contributes to the better fatigue durability of composite risers, which was confirmed by previous studies [13,14,16]. More specifically, the structural composite body of the composite riser joint is likely to have an infinite fatigue life [13,14,16], but still may vary depending on the choice of constituent materials and the manufacturing process. In terms of the VIV of composite risers, it is found that the fundamental frequency of a composite riser with the same boundary condition and geometry was higher than that of a steel riser [18] since the composite riser was much lighter and the value was relatively small [13,14]. The structural damping and tension variations were also significant to the VIV of composite risers [13], i.e., higher tension and damping could reduce VIV-induced fatigue [13,17]. A. F. Omar et al. [18] found that the maximum VIV stresses of a composite riser were much smaller than those of a steel riser, indicating the better fatigue lives of composite risers. Similarly, William Toh et al. [25] and L.B. Tan et al. [24] conducted comprehensive studies on both local and global scales on composite and steel risers using a finite element method and a larger safety margin for composite risers were observed. Although a composite riser without VIV suppression had a moderate VIV-induced fatigue damage, adding strakes and buoyancy modules to a composite riser could lead to even better VIV responses [13,24].
In addition, the design of a composite riser (thick composite tube) can be optimized to specific requirements by modifying some of its design variables, such as fiber and matrix combinations, fiber orientations, thicknesses of the liner and composite lamina, and stacking sequence, to enhance the benefits offered by FRP composites risers using different methods, e.g., finite element analysis [22,26,27,28,29], simple-input displacement-based method [30,31] and optimization algorithms [21,32,33,34].
Besides all these research studies, joint industry projects emerged several decades ago, and Table 1 summarizes some of them.
From all these research studies and joint industry projects, the feasibility, advantages and effectiveness of using FRP composites for risers have been proven, although no offshore riser currently being employed is made entirely of composites, except the demonstration composite drilling riser joint (a tube segment) installed in the field on the Heidrun TLP [38].
In this paper, three different risers with large length to diameter ratio, the X80 steel riser (riser 1), AS4-epoxy composites with titanium liner with geometries of (liner/90/(0/90)10) (riser 2) and (liner/03/(+53, −53)5/904) (riser 3) are selected to study their VIV characteristics. The CFD simulation with coupled fluid–structure interaction was achieved using ANSYS Workbench 17.1. Based on the combination of different risers, tension forces, buoyancies, gravities and current conditions, nine study cases were considered. The VIV responses under different situations of all these three risers were obtained, including the natural frequencies, global displacements, global stresses and the stress distributions in every layers of both FRP composite risers.

2. Materials and Methods

2.1. Materials for the Risers

In this study, X80 steel was utilized for steel riser and for composite risers, titanium alloy (Ti-6Al-4V) was used for the inner liner and AS4-epoxy was employed for the composite layers. The properties of these materials can be found in Table 2 and Table 3 as well as in Figure 1 and Figure 2. For the X80 steel (for riser 1) and titanium (for liner in risers 2 and 3), bilinear kinematic hardening model was selected to consider the elastic–plastic stage. For steel riser 1, 67% of the yield stress of X80 was identified as the allowable stress according to American Bureau of Shipping (ABS) standard [45]. For the titanium liner in risers 2 and 3, 59.5% of the yield stress of titanium was regarded as the allowable stress [46]. In terms of the properties of unidirectional lamina of AS4-Epoxy FRP composites, an elastic model was considered, i.e., the plastic stage was disregarded. The maximum stress failure criterion [47] considered failures in the fiber direction, transverse direction and in-plane shear separately for all composite layers. In Table 3, subscripts 1, 2 and 3 represent the fiber direction, in-plane transverse direction and through-thickness direction of a composite lamina, respectively (Figure 2).

2.2. Riser Geometries

Figure 3a,b respectively show the cross-sections of the steel and composite riser joints in this study. The extracted fluids were transported by production tubings within the riser, which were protected by the riser joints and the riser joints carried all the structural loads. The internal diameters (IDs) of the steel and composite riser joints were fixed at 250 mm while, as the thickness of each joint was determined by the design to accommodate all the loads considered, their outer diameters (ODs) depended on the previous design results [20,21,22,27] (Table 4). From Table 4, the thicknesses for the steel riser (riser 1), AS4-epoxy composites and titanium liner with (liner/90/(0/90)10) (riser 2) and AS4-epoxy composites and titanium liner with (liner/03/(+53, −53)5/904) (riser 3) were 25 mm, 39.5 mm and 30.5 mm, respectively. In other words, riser 2 was a conventional design of FRP composite riser while riser 3 was an optimized design of FRP composite riser. Both FRP composite risers (59.6 kg/m and 45.7 kg/m) achieved significant weight saving compared to the steel riser (169 kg/m) and the optimized FRP composite riser 3 even led to a 23.3% weight reduction to the conventional FRP composite riser 2. For all the risers in this study the length to inner diameter ratios were 100. The fiber angle of 0° and 90° were the axial and hoop directions of FRP composite risers, respectively (Figure 4).

2.3. FE Modelling and Study Cases

ANSYS Workbench 17.1 was utilized for the comparative study of vortex-induced vibration (VIV) for the FRP composite risers and steel riser with length to diameter ratio of 100. Initially, the natural frequencies of all the risers were analyzed using pre-stress modal analysis. For FRP composite risers, ACP (Pre), Static Structural and Modal modules were involved while for steel risers, only Static Structural and Modal modules were required. In order to achieve the VIV simulation in ANSYS Workbench 17.1, the following modules were included: Geometry, ACP (Pre), Transient Structural, Fluid Flow, ACP (Post) and System Coupling [19]. Generally speaking, Geometry was utilized to create the geometries of risers and flow zones; ACP (Pre) and ACP (post) were employed to include the data of materials and lamination information for composite risers (risers 2 and 3) and to show the stress and deformation response of every composite laminae, respectively; Transient Structural was used to simulate the dynamic responses of risers; System Coupling was applied to achieve the data transfer between Transient Structural and Fluid Flow. As mentioned in Section 2.2, the geometries of composite risers 2 and 3 were obtained from the design results of Wang [20,21,27] and the finite element modelling method had been verified by Kim [14] and Amaechi [26]. The difference here was the load conditions which was current flow with different velocities.
In this study, layered-structure method (LSM) [19] was employed for FRP composite riser simulation and their geometries are listed in Table 4. For the FRP composite risers, the data of material properties, the thicknesses of liner and composite laminate, reinforced fiber angles, the layer numbers of each lamina with different fiber orientations and the laminate stacking sequences were inputted in the module of ACP (Pre) and the stress and deformation response of every composite laminae were obtained in ACP (Post). The possible bond-slip effect is a significant phenomenon of FRP layers, and this can be checked indirectly in post-processing stage by comparing the maximum interface forces and stresses with the reference pull-out resistances [48]. However, in this paper, the numerical model of the composite riser using LSM assumed a rigid connection between the components, i.e., perfect bonding was assumed and no bond-slip effect was considered to simplify the modelling.
Then, the geometry and material data of FRP composite riser in ACP (Pre) were transferred to the Transient Structural module. Transient Structural was utilized to simulate the dynamic response of risers. In this module, the gravity, top-tension force, buoyancy, end support condition (simple support at the top end, fixed support at the bottom end) and the fluid solid interface were employed. After the whole simulation, the displacement history, global stress distribution etc., could be obtained.
Fluid Flow (Fluent) was used to simulate the flow conditions. 10D (Y-direction) × 20D (X-direction) × 25 m (Z-direction) was set as the flow dimension, where D was the outside diameter of the risers. Here X, Y and Z directions were in-line flow direction, cross flow direction and depth direction, respectively. Large eddy simulation (LES) with bounded second order implicit transient formulation was selected as the mathematical model for turbulence in computational fluid dynamics (CFD), which reduced the computational cost with acceptable accuracy. The flow properties, such as flow density, kinematic viscosity coefficient and flow velocity listed in Table 5 were inputted in this module. The inlet, outlet, left, right, top and bottom sides were defined as the boundary conditions of velocity inlet, outflow, symmetry and wall, respectively.
To achieve the data transfer between Transient Structural and Fluid Flow, i.e., the two-way fluid structure interaction (FSI), module of System Coupling was utilized, and the simulation time (80 s) and step size (0.2 s) were set in this module.
The fluid grid of 5D × 5D area around the riser was refined and set as the dynamic mesh zone with smoothing and re-meshing methods. Based on the results of convergence study, there were 1.15 million elements in the refined zone and 0.52 million elements for the other flow area, i.e., 1.67 million elements in total were employed in the flow zone. For the risers, 300 elements in axial direction and 20 elements in hoop direction were employed. In the though-thickness direction, one element, 22 elements and 18 elements for risers 1, 2 and 3, respectively. Note that, the laminate data were transferred from ACP (Pre) to Transient Structural and the number of elements in though-thickness direction was set as the number of layers of the structure automatically. In total, there were 6000 elements, 132,000 elements and 108,000 elements for risers 1, 2 and 3 respectively.
In this study, a tension was applied to the riser’s top to keep its vertical position and the gravity, buoyancy and current load were also employed. The tension forces were set as 1.5 and 2 times the structural weight of steel and composite risers, respectively. The gravity and buoyancy loads were determined by the geometry of each riser. The current data of the Gulf of Mexico [2,14] was utilized according to the API Recommended Practice [49], which were 0.36 m/s for a 1-year winter storm, 1.22 m/s for a 100-year hurricane and 2.13 m/s for a 100-year loop current conditions. In the Gulf of Mexico, the sea water density is 1024 kg/m3 and the kinematic viscosity coefficient υ is 1.06 × 10−6. Table 5 shows the different cases in this study. Reynolds numbers for all the study cases were determined by Equation (1), where U is the flow velocity, D is the outside diameter of the riser, υ is kinematic viscosity coefficient. Turbulent intensity employed were calculated by Equation (2), where I is turbulent intensity and Re is Reynolds number.
Re = UD υ
I = 0.16 Re ( 1 / 8 )
From Table 5, it is clear that the top-tension forces required for composite risers were much smaller than that of the steel riser, although the tension force ratios of composite risers were larger (1.5 and 2 respectively). This illustrates the benefit of composite risers in the weight reduction and the benefit would be more significant with the increasing length of risers.

3. Results and Discussion

In this section, the VIV responses under different situations of the three different risers were obtained and discussed. The large length to diameter ratio, fluid–structure interaction, different current velocities and different design geometries of FRP composite risers would affect each riser’s VIV characteristics significantly.

3.1. Natural Frequencies

Natural frequency of the riser is a key characteristic in VIV analysis, because the “lock-in” would occur when the natural frequency of a riser is close to the vortex shedding frequency. The “lock-in” phenomenon would lead to fatigue failure of the riser because of the large and damaging vibration. The vortex shedding frequency can be calculated by Equation (3) which includes the Strouhal number St, the outside diameter of the riser D and the flow velocity U. According to the relationship between Re and St, the Strouhal number St is 0.2 for Re between 300 and 300,000 and is 0.2–0.3 for Re between 300,000 and 3,500,000 [50]. More specifically, the St are 0.2 for cases 1–3, 0.21 for cases 4–6 and 0.22 for cases 7–9.
f s = S t U D
In addition, the possibility of a riser’s “lock-in” using reduced velocity Ur was estimated using Equation (4). When Ur ranges 4–8 [50], “lock-in” may occur. The Ur of all the risers are listed in Table 6.
U r = U f n D
The natural frequencies and first three modal shapes for all the risers are presented in Figure 5 and Table 6. Note that all the risers have similar modal shapes and therefore only the first three modal shapes of FRP composite riser 3 are showed as an example.
From the results listed in Table 6, it can be seen the natural frequencies of cases 7 (riser 1) and case 9 (riser 3) are close to their vortex shedding frequencies, and similarly the reduced velocity Ur of these two cases are close to the range of 4–8. Hence, for these two cases, “lock-in” might occur. We have to notice that Equation (3) of vortex shedding frequency and Equation (4) of reduced velocity are calculated based on the rigid cylinder without deformation and fluid structure interaction.

3.2. Time History of Displacements for All Three Risers

The time history of displacements in in-line flow direction (X) and cross flow direction (Y) for all risers in different cases are plotted in Figure 6 and Figure 7, respectively. In order to show the comparison more clearly, the displacements in both directions were normalized by the O.D. of the risers. For all the risers, three typical locations, i.e., −6 m, −12.5 m and −19 m, were selected.
From Figure 6, it is clear that the in-line flow direction (X) displacements maintained at a relatively stable value and vibrate slightly after the initial growth for all three risers and at all locations. More specifically, for velocities equal to 0.36 m/s, the maximum displacements in flow direction/outside diameter (Ax/D) of risers 1, 2 and 3 were around 0.0017, 0.003 and 0.0065 at −6 m, 0.002, 0.0036 and 0.008 at −12.5 m and 0.0009, 0.0016 and 0.0035 at −19 m, respectively. For velocities equal to 1.22 m/s, the maximum displacements in flow direction/outside diameter (Ax/D) of risers 1, 2 and 3 were 0.02, 0.035 and 0.075 at −6 m, 0.024, 0.042 and 0.095 at −12.5 m and 0.011, 0.019 and 0.042 at −19 m, respectively. For velocities equal to 2.13 m/s, the maximum displacements in flow direction/outside diameter (Ax/D) of risers 1, 2 and 3 were 0.06, 0.10 and 0.23 at −6 m, 0.075, 0.128 and 0.274 at −12.5 m and 0.033, 0.057 and 0.12 at −19 m, respectively. To summarize: (1) In the cases with the same current velocity, the displacement of riser 3 was the largest, followed by riser 2, and riser 1 was the smallest for all locations of the three risers. (2) With the increase of current velocity, the Ax/D rose for all three risers. Take riser 3 at 12.5 m for example, the Ax/D(s) was about 11 times (1.22 m/s) and 32 times (2.13 m/s) of that at 0.36 m/s. (3) The end support conditions affected the displacement significantly, i.e., simple support lead to larger displacement. Using riser 3 at 2.13 m/s as an example, the Ax/D(s) was about 0.23 at −6 m (6 m to the simple support end) and 0.12 at −19 m (6 m to the fixed support end).
From Figure 7, it is clear that the cross-flow direction (Y) displacements of all risers vibrated obviously. The vibrations for the cases with 0.36 m/s were comparatively more stable than those for the cases with 1.22 m/s and 2.13 m/s. This was because for cases 1–3 (0.36 m/s), the Re (s) was 101,887, 111,736 and 105,623 and therefore the vortex streets for these cases were fully turbulent [50]. For cases 4–9, the Re (s) was between 300,000 and 3,500,000 and therefore the laminar boundary layer underwent turbulent transition and the wake was narrower and disorganized [50]. More specifically, for velocities equal to 0.36 m/s, the maximum range of displacements in cross-flow direction/outside diameter (Ay/D) of riser 1, 2 and 3 were around 0.0004, 0.001 and 0.002 at −6 m, 0.0005, 0.0018 and 0.003 at −12.5 m and 0.0002, 0.0006 and 0.0013 at −19 m, respectively. For velocities equal to 1.22 m/s, the maximum ranges of displacements in cross-flow direction/outside diameter (Ay/D) of risers 1, 2 and 3 were 0.005, 0.014 and 0.024 at −6 m, 0.006, 0.017 and 0.026 at −12.5 m and 0.003, 0.0075 and 0.012 at −19 m, respectively. For velocities equal to 2.13 m/s, maximum ranges of displacements in cross-flow direction/outside diameter (Ay/D) of risers 1, 2 and 3 were 0.018, 0.035 and 0.059 at −6 m, 0.022, 0.043 and 0.073 at −12.5 m and 0.009, 0.018 and 0.033 at −19 m, respectively. To summarize: (1) In the cases with the same current velocity, the vibration range of riser 3 was the largest, followed by riser 2, and riser 1 was the smallest for all locations of the three risers. (2) With the increase of current velocity, the Ay/D rose for all three risers. Take riser 3 at 12.5 m for an example, the Ay/D(s) was about 9 times (1.22 m/s) and 24 times (2.13 m/s) of that at 0.36m/s. (3) the end support conditions affected the displacement significantly, i.e., the simple support lead to larger vibration. Using riser 3 at 2.13 m/s as an example, the Ay/D(s) were about 0.059 at −6 m (6 m to the simple support end) and 0.033 at −19 m (6 m to the fixed support end).

3.3. Maximum Displacements and Von Mises Stresses of All Three Risers

Similarly, the time history of displacements in in-line flow direction (X), cross flow direction (Y) and von Mises stresses for all locations of risers can be acquired. According to these data, the maximum displacements and stresses at different locations of risers are obtained and presented in Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12.
Figure 8a shows the maximum Ax/D for a velocity of 0.36 m/s, Figure 8b presents the maximum Ax/D for velocity of a 1.22 m/s, and Figure 8c illustrates the maximum Ax/D for a velocity of 2.13 m/s.
From Figure 8, for all cases riser 3 had the maximum Ax/D and riser 1 had the minimum Ax/D. Also it could be found that the maximum Ax/D occurred above the middle point of each riser, which indicated the simple support lead to larger displacement compared with the fixed support. Based on Figure 8 and the procedure of creating it, the maximum Ax/D for risers 1, 2 and 3 for 0.36 m/s were 0.0023 at −11 m (60.6 s), 0.0039 at −11 m (74.6 s) and 0.0085 at −11 m (62.6 s), respectively. The maximum Ax/D for risers 1, 2 and 3 for 1.22 m/s were 0.0263 at −11 m (71.0 s), 0.0455 at −11 m (28.8 s) and 0.100 at −11 m (68.0 s), respectively. The maximum Ax/D for risers 1, 2 and 3 for 2.13 m/s were 0.0805 at −11 m (46.6 s), 0.140 at −11 m (39.8 s) and 0.290 at −11 m (54.4 s), respectively.
Figure 9a shows the maximum Ay/D for velocity of 0.36 m/s, Figure 9b presents the maximum Ay/D for velocity of 1.22 m/s, and Figure 9c illustrates the maximum Ay/D for velocity of 2.13 m/s.
From Figure 9, for all cases riser 3 had the maximum ranges of Ay/D and riser 1 had the minimum ranges of Ay/D. Also it could be found that the maximum ranges of Ay/D occurred above the middle point of each riser, which indicated the simple support lead to larger vibrations compared with the fixed support. Based on Figure 9 and the procedure of creating it, the maximum ranges of Ay/D for risers 1, 2 and 3 for 0.36 m/s were 0.00048 at −11 m, 0.0019 at −11 m and 0.003 at −11 m, respectively. The maximum ranges of Ay/D for risers 1, 2 and 3 for 1.22 m/s were 0.0063 at −11 m, 0.018 at −11 m and 0.0282 at −9 m, respectively. The maximum ranges of Ay/D for risers 1, 2 and 3 for 2.13 m/s were 0.0233 at −10 m, 0.044 at −11 m and 0.076 at −11 m, respectively. Here it is noted that the maximum Ay/D and the minimum Ay/D occurred at different times, and therefore no time was measured for the maximum range of Ay/D.
From Figure 8 and Figure 9, the in-line displacement was much larger than the cross flow displacement. Based on the time when the maximum Ax/D occur for all cases, Figure 10 presents the total displacement distributions of riser 3 at those specified times. Here we have to note, the displacement distributions of risers 1 and 2 were similar to those of riser 3, and the only difference was the specific maximum displacement value for each riser. More specifically, the maximum displacements for riser 1 were 0.00071 m at 60.6 s, 0.0079 m at 71 s and 0.024 m at 46.6 s for flow velocities of 0.36 m/s, 1.22 m/s and 2.13 m/s respectively, the maximum displacements for riser 2 were 0.0013 m at 74.6 s, 0.015 m at 28.8 s and 0.046 m at 39.8 s for flow velocities of 0.36 m/s, 1.22 m/s and 2.13 m/s respectively, and the maximum displacements for riser 3 were 0.0027 m at 62.6 s, 0.031 m at 68 s and 0.090 m at 54.4 s for flow velocities of 0.36 m/s, 1.22 m/s and 2.13 m/s respectively.
Figure 11 shows the maximum von Mises stresses of all three risers for the velocity of 0.36 m/s, 1.22 m/s, and 2.13 m/s, respectively.
From Figure 11, for all cases riser 3 had the maximum von Mises stresses and riser 2 had the minimum von Mises stresses. Here we have to note that, only the von Mises stresses for steel riser (riser 1) could be utilized to verify the failure, more specifically, the 67% of the yield stress of X80 steel (371 MPa) was identified as the allowable stress according to American Bureau of Shipping (ABS) standard [45]. Therefore, for all cases, steel risers would not incur stress failure. The minimum factors of safety, which are defined as allowable strength/stress, were 82.4, 31.6 and 12.1 for steel riser 1 with current velocities of 0.36 m/s, 1.22 m/s and 2.13 m/s, respectively. However, the von Mises stresses for FRP composite risers (risers 2 and 3) could not be utilized to verify the failure because the global stresses could not represent the stress distributions in each composite laminae. If the equivalent stresses of a composite riser is required, the effective modulus method (EMM) rather than the layered-structure method (LSM) should be utilized [19]. In this paper, the global stress distributions for FRP composite risers here were only used to find the locations and time of the maximum stresses occurred, and this information was used for failure verification of FRP composite risers in the next section.
For all cases, the maximum von Mises stresses occurred at the top/bottom end of all risers, followed by the middle parts (above the middle point) of the risers, while at about −19 m, the minimum von Mises stresses happened. The reasons for this trend were that the fixed support at the bottom lead to the maximum stresses and the middle parts had the largest deformation attributed to large stresses as well. For the small velocity situation, top-tension force lead to the maximum von Mises stress at the top end, especially for the steel riser whose tension force was the largest. For the location around −19 m, all risers had comparatively small deformation and it was relatively far from the fixed support, therefore, minimum von Mises stresses occurred.
Based on Figure 11 and the procedure of creating it, the maximum von Mises stresses for risers 1, 2 and 3 for 0.36 m/s were 4.5 MPa at 0 m (72.6 s), 3.0 MPa at −25 m (30.6 s) and 5.7 MPa at −7 m (24.6 s), respectively. The maximum von Mises stresses for risers 1, 2 and 3 for 1.22 m/s were 11.7 MPa at −25 m (71.0 s), 11.3 MPa at −25 m (52.6 s) and 22.7 MPa at −25 m (68.0 s), respectively. The maximum von Mises stresses for risers 1, 2 and 3 for 2.13 m/s were 30.7 MPa at −25 m (46.6 s), 29.4 MPa at −25 m (39.8 s) and 59.2 MPa at −25 m (54.4 s), respectively.
According to the time when the maximum von Mises stresses occurred for all cases, Figure 12 presents the stress distributions of all risers at those specified time.

3.4. Stress Distributions in Every Layer of the FRP Composite Risers

For each case, there were two significant times, when the maximum displacement occurred and the maximum stresses occurred. Comparing the stresses in the FRP composite layers at these two time points, the maximum stresses in every lamina of FRP composite risers (risers 2 and 3) were obtained and the factor of safety (FS) are presented in Figure 13, Figure 14 and Figure 15. Here, the factor of safety was defined as allowable strength/stress.
In composite riser design, the factor of safety was utilised to verify the stress failure. In this paper, the maximum stress failure criterion [47] considered failures in the fiber direction, transverse direction and in-plane shear separately for all composite layers. The minimum FS required was 1.53 for composite layers and 1.68 for titanium liner [46].
For cases with velocity of 0.36 m/s for riser 2, the maximum von Mises stress in the titanium liner is 3 MPa and FS = 293. Figure 13a,b, respectively, show the minimum FSs in the fiber and transverse directions for all the layers of riser 2. Here, it was noted that for the riser geometry with only orthotropic reinforcements (riser 2), the shear stresses in all layers were small enough to be neglected. The minimum FS in the fiber direction was 493 (layer 20 in Figure 13a) while that in the transverse direction is 168 (layers 21 in Figure 13b). It is evident that, for riser 2 with 0.36 m/s, the in-plane transverse stresses were the most critical stresses and they were far away from the failure stresses.
For cases with velocity of 0.36 m/s for riser 3, the maximum von Mises stress in the titanium liner was 5.2 MPa and FS = 169. Figure 13c–e, respectively, show the minimum FSs in the fiber, transverse and shear directions for all the layers of riser 3. The minimum FS in the fiber direction were 304 (layer 3 in Figure 13c), 73.7 (layer 17 in Figure 13d in the transverse direction and 273 in shear direction (layer 13 in Figure 13e). It was evident that, for riser 3 with 0.36 m/s, the in-plane transverse stresses were the most critical stresses and much smaller than those in riser 2. They were, however, far away from the failure stresses.
Comparing the FS (s) in riser 2 and riser 3, it was obvious that the tailor-designed composite riser 3 could take advantage of the fiber strength more efficiently.
For cases with velocity of 1.22 m/s for riser 2, the maximum von Mises stress in the titanium liner was 11.2 MPa and FS = 78. Figure 14a,b, respectively, show the minimum FSs in the fiber and transverse directions for all the layers of riser 2. Here, it was noted that for the riser geometry with only orthotropic reinforcements (riser 2), the shear stresses in all layers were small enough to be neglected. The minimum FS in the fiber direction was 120 (layer 20 in Figure 14a) while that in the transverse direction was 43.7 (layers 21 in Figure 14b). It was evident that, for riser 2 with 1.22 m/s, the in-plane transverse stresses were the most critical stresses and they were far away from the failure stresses.
For cases with velocity of 1.22 m/s for riser 3, the maximum von Mises stress in the titanium liner was 22.4 MPa and FS = 39. Figure 14c–e, respectively, show the minimum FSs in the fiber, transverse and shear directions for all the layers of riser 3. The minimum FS in the fiber direction were 71.4 (layer 3 in Figure 14c), 24.2 (layers 17 in Figure 14d) in the transverse direction and 60.3 in the shear direction (layer 13 in Figure 14e). It was evident that, for riser 3 with 1.22 m/s, the in-plane transverse stresses were the most critical stresses and much smaller than those in riser 2. But they were far away from the failure stresses.
Comparing the FS (s) in riser 2 and riser 3, it was obvious that tailor-designed composite riser 3 could take advantage of the fiber strength more efficiently.
For cases with velocity of 2.13 m/s for riser 2, the maximum von Mises stress in the titanium liner was 29.3 MPa and FS = 30. Figure 15a,b, respectively, show the minimum FSs in the fiber and transverse directions for all the layers of riser 2. Here, it was noted that for the riser geometry with only orthotropic reinforcements (riser 2), the shear stresses in all layers were small enough to be neglected. The minimum FS in the fiber direction was 35.3 (layer 20 in Figure 15a) while that in the transverse direction was 16.1 (layers 21 in Figure 15b). It was evident that, for riser 2 with 2.13 m/s, the in-plane transverse stresses were the most critical stresses and they were far away from the failure stresses.
For cases with velocity of 2.13 m/s for riser 3, the maximum von Mises stress in the titanium liner was 58.1 MPa and FS = 15. Figure 15c–e, respectively, show the minimum FSs in the fiber, transverse and shear directions for all the layers of riser 3. The minimum FS in the fiber direction were 22.5 (layer 3 in Figure 15c), 12.9 (layers 17 in Figure 15d in the transverse direction and 22.8 in shear direction (layer 13 in Figure 15e). It was evident that, for riser 3 with 2.13 m/s, the in-plane transverse stresses were the most critical stresses and much smaller than those in riser 2. But they were far away from the failure stresses.
Comparing the FS (s) in riser 2 and riser 3, it was obvious that tailor-designed composite riser 3 could take advantage of the fiber strength more efficiently.

4. Conclusions

This paper presents a comparative VIV study of FRP composite risers (AS4-epoxy with titanium liner) with two geometries ([liner/90/(0/90)10] and [liner/03/(+53, −53)5/904]) using a steel riser as the benchmark. The CFD simulation with coupled fluid–structure interaction was employed for nine different study cases considering a combination of different risers, tension forces, buoyancies, gravities and current conditions to investigate each riser’s VIV characteristics including their natural frequencies, global displacements, global stresses and the stress distributions in every layer of the FRP composite risers. The results showed that:
(1) VIV responses of displacements for all three risers at 2.13 m/s were about 30 times of those at 0.36 m/s and 3 times of those at 1.22 m/s; VIV responses of global stresses for all three risers at 2.13 m/s were about 10 times of those at 0.36 m/s and 3 times of those at 1.22 m/s; and only in the 100 year loop current condition (2.13 m/s), “lock-in” might occur for steel riser (riser 1) and optimized FRP composite riser (riser 3).
(2) Displacements in in-line flow direction maintained at a relatively stable value and vibrated slightly after the initial growth, while the vibration in cross flow direction was obvious. In terms of the total displacement, that of the optimized FRP composite riser was the largest and that of steel riser was the smallest, i.e., the maximum total displacement for riser 3 was about three times of that for riser 1 and two times of that for riser 2 at all flow velocities.
(3) End support conditions affected the displacement significantly, i.e., simple support lead to larger displacement compared to the fixed support.
(4) Maximum von Mises stresses occurred at the top/bottom end of all risers, followed by the middle parts (above the middle point) of the risers, while at about −19 m, the minimum von Mises stresses happened. When the maximum von Mises stresses occurred at the top end, the top tension force played a key role, and when the maximum von Mises stresses occurred at the bottom end, the bending force was the most significant one. For the middle part, both tension and bending forces were essential.
(5) The conventional FRP composite riser (riser 2) had the largest FS (168, 43.7 and 16.1 for 0.36 m/s, 1.22 m/s and 2.13 m/s, respectively) while the optimized FRP composite riser (73.7, 24.2 and 12.9 for 0.36 m/s, 1.22 m/s and 2.13 m/s, respectively) and steel riser (82.4, 31.6 and 12.1 for 0.36 m/s, 1.22 m/s and 2.13 m/s, respectively) had the similar FS, indicating the optimized composite riser could take advantage of the fiber strength more efficiently.

Author Contributions

Conceptualization, C.W.; methodology, C.W. and S.G.; validation, C.W. and S.G.; formal analysis, S.G.; investigation, M.S.; data curation, S.G.; writing—original draft preparation, C.W. and S.G.; writing—review and editing, C.W.; supervision, C.W. and B.H.; project administration, C.W. and Z.J.; funding acquisition, C.W.

Funding

This research was funded by the Shandong Provincial Natural Science Foundation (ZR2017LEE031), National Natural Science Foundation of China (51806130) and the Science and Technology Development Project of Zibo (2017KJ040017).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. American Petroleum Institute (API). Design of Risers for Floating Production Systems (FPSs) and Tension-Leg Platforms (TLPs); American Petroleum Institute: Washington, DC, USA, 1998. [Google Scholar]
  2. Chakrabarti, S.K. Handbook of Offshore Engineering, 1st ed.; Elsevier: London, UK, 2005; Volume 1. [Google Scholar]
  3. Balazs, G.L.; Borosnyoi, A. Long-term behavior of FRP. In Composites in Construction: A Reality; Edoardo, C., Gaetano, M., Antonio, N., Eds.; American Society of Civil Engineers: Reston, VA, USA, 2001; pp. 84–91. [Google Scholar]
  4. Jansons, J.O.; Glejbol, K.; Rytter, J.; Aniskevich, A.N.; Arnautov, A.K.; Kulakov, V.L. Effect of water absoption, elevated temperatures and fatigue on the mechanical properties of carbon-fiber-reinforced epoxy composites for flexible risers. Mech. Compos. Mater. 2002, 38, 299–310. [Google Scholar] [CrossRef]
  5. Salama, M.M.; Martinussen, E.; Spencer, B.; Hanna, S.; Hsu, T.M.; Stjern, G.; Franco, E.; Storhaug, T.; Echtermeyer, A. Composite risers are ready for field applications-status of technology, field demonstration and life cycle economics. In Proceedings of the 13th Annual Deep Offshore Technology Conference, Rio de Janerio, Brazil, 17–19 October 2001. [Google Scholar]
  6. Venkatesan, R.; Dwarakadasa, E.S.; Ravindran, M. Study on behavior of carbon fiber-reinforced composite for deep sea applications. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 6–9 May 2002. [Google Scholar]
  7. Ochoa, O.O.; Salama, M.M. Offshore composites: Transition barriers to an enabling technology. Compos. Sci. Technol. 2005, 65, 2588–2596. [Google Scholar] [CrossRef]
  8. Ochoa, O.O.; Ross, G.R. Hybrid composites: Models and tests for environmental aging. J. Reinf. Plast. Compos. 1998, 17, 787–799. [Google Scholar] [CrossRef]
  9. Ross, G.R.; Ochoa, O.O. Environmental effects on unsymmetric composite laminates. J. Thermoplast. Compos. Mater. 1991, 4, 266–284. [Google Scholar] [CrossRef]
  10. Bismarck, A.; Hofmeier, M.; Dörner, G. Effect of hot water immersion on the performance of carbon reinforced unidirectional poly(ether ether ketone) (peek) composites: Stress rupture under end-loaded bending. Compos. Part A 2007, 38, 407–426. [Google Scholar] [CrossRef]
  11. Bedon, C.; Louter, C. Numerical investigation on structural glass beams with GFRP-embedded rods, including effects of pre-stress. Compos. Struct. 2018, 184, 650–661. [Google Scholar] [CrossRef]
  12. Karayaka, M.; Wu, S.; Wang, S.; Lu, X.; Partha Ganguly, B.P. Composite production riser dynamics and its effects on tensioners, stress joints, and size of deep water tension leg platform. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 4–7 May 1998. [Google Scholar]
  13. Huang, K.Z. Composite TTR design for an ultradeepwater TLP. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 2–5 May 2005. [Google Scholar]
  14. Kim, W.K. Composite Production Riser Assessment. Ph.D. Thesis, Texas A&M University, College Station, TX, USA, 2007. [Google Scholar]
  15. Beyle, A.I.; Gustafson, C.G.; Kulakov, V.L.; Tarnopol’skii, Y.M. Composite risers for deep-water offshore technology: Problems and prospects. 1. Metal-composite riser. Mech. Compos. Mater. 1997, 33, 403–414. [Google Scholar] [CrossRef]
  16. Ward, E.G.; Ochoa, O.O.; Kim, W.; Gilbert, R.M.; Jain, A.; Miller, C.; Denison, E. A Comparative Risk Analysis of Composite and Steel Production Risers; Texas A&M University: College Station, TX, USA, 2007. [Google Scholar]
  17. Rakshit, T.; Atluri, S.; Dalton, C. VIV of a composite riser at moderate Reynolds number using CFD. J. Offshore Mech. Arct. Eng. 2008, 130, 011009. [Google Scholar] [CrossRef]
  18. Omar, A.F.; Karayka, M.; Murray, J.J. A comparative study of the performance of top-tensioned composite and steel risers under vertex-induced loading. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 3–6 May 1999. [Google Scholar]
  19. Wang, C.; Sun, M.; Shankar, K.; Xing, S.; Zhang, L. CFD simulation of vortex induced vibration for FRP composite riser with different modeling methods. Appl. Sci. 2018, 8, 684. [Google Scholar] [CrossRef]
  20. Wang, C.; Shankar, K.; Morozov, E.V. Global design and analysis of deep sea FRP composite risers under combined environmental loads. Adv. Compos. Mater. 2017, 26, 79–98. [Google Scholar] [CrossRef]
  21. Wang, C.; Shankar, K.; Ashraf, M.A.; Morozov, E.V.; Ray, T. Surrogate-assisted optimisation design of composite riser. Proc. Inst. Mech. Eng. Part L 2016, 230, 18–34. [Google Scholar] [CrossRef]
  22. Wang, C.; Shankar, K.; Morozov, E.V. Tailored design of top-tensioned composite risers for deep-water applications using three different approaches. Adv. Mech. Eng. 2017, 9, 1–18. [Google Scholar] [CrossRef]
  23. Chen, Y.; Seemann, R.; Krause, D.; Tay, T.-E.; Tan, V.B. Prototyping and testing of composite riser joints for deepwater application. J. Reinf. Plast. Compos. 2016, 35, 95–110. [Google Scholar] [CrossRef]
  24. Tan, L.B.; Chen, Y.; Jaiman, R.K.; Sun, X.; Tan, V.B.C.; Tay, T.E. Coupled fluid–structure simulations for evaluating a performance of full-scale deepwater composite riser. Ocean Eng. 2015, 94, 19–35. [Google Scholar] [CrossRef]
  25. Toh, W.; Tan, L.B.; Jaiman, R.K.; Tay, T.-E. A comprehensive study on composite risers: Material solution, local end fitting design and global response. Mar. Struct. 2018, 61, 155–169. [Google Scholar] [CrossRef]
  26. Amaechi, C.V.; Nathaniel, G.; Agbomerie, C.O.; Hou, X.; Ye, J. Composite risers for deep waters using a numerical modelling approach. Compos. Struct. 2018, 210, 486–499. [Google Scholar] [CrossRef]
  27. Wang, C.; Shankar, K.; Morozov, E.V. Tailored local design of deep sea FRP composite risers. Adv. Compos. Mater. 2015, 24, 375–397. [Google Scholar] [CrossRef]
  28. Harte, A.M.; McNamara, J.F.; Roddy, I.D. Application of optimisation methods to the design of high performance composite pipelines. J. Mater. Process. Technol. 2003, 142, 58–64. [Google Scholar] [CrossRef]
  29. Harte, A.M.; McNamara, J.F.; Roddy, I. Evaluation of optimisation techniques in the design of composite pipelines. J. Mater. Process. Technol. 2001, 118, 478–484. [Google Scholar] [CrossRef]
  30. Yazdani Sarvestani, H.; Hoa, S.V.; Hojjati, M. Effects of shear loading on stress distributions at sections in thick composite tubes. Compos. Struct. 2016, 140, 433–445. [Google Scholar] [CrossRef]
  31. Yazdani Sarvestani, H.; Hoa, S.V.; Hojjati, M. Three-dimensional stress analysis of orthotropic curved tubes-part 1: Single-layer solution. Eur. J. Mech. 2016, 60, 327–338. [Google Scholar] [CrossRef]
  32. Ghiasi, H.; Pasini, D.; Lessard, L. Optimum stacking sequence design of composite materials part I: Constant stiffness design. Compos. Struct. 2009, 90, 1–11. [Google Scholar] [CrossRef]
  33. Ghiasi, H.; Fayazbakhsh, K.; Pasini, D.; Lessard, L. Optimum stacking sequence design of composite materials part II: Variable stiffness design. Compos. Struct. 2010, 93, 1–13. [Google Scholar] [CrossRef]
  34. Da Silva, R.F.; Teófilo, F.A.; Parente, E., Jr.; de Melo, A.M.; de Holanda, Á.S. Optimization of composite catenary risers. Mar. Struct. 2013, 33, 1–20. [Google Scholar] [CrossRef]
  35. Ahlstone, A.G. Light Weight Marine Riser Pipe. U.S. Patent 3,768,842, 30 October 1973. [Google Scholar]
  36. Sparks, C.P.; Odru, P.; Bono, H.; Metivaud, G. Mechanical testing of high-performance composite tubes for tlp production risers. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 2–5 May 1988. [Google Scholar]
  37. Salama, M.M.; Johnson, D.B.; Long, J.R. Composite production riser-testing and qualification. SPE Prod. Facil. 1998, 13, 170–177. [Google Scholar] [CrossRef]
  38. Salama, M.M.; Stjern, G.; Storhaug, T.; Spencer, B.; Echtermeyer, A. The first offshore field installation for a composite riser joint. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 6–9 May 2002. [Google Scholar]
  39. Smith, K.L.; Leveque, M.E. Ultra-Deepwater Production Systems Technical-Progress Report; ConocoPhillips Company: Houston, TX, USA, 2003. [Google Scholar]
  40. Smith, K.L.; Leveque, M.E. Ultra-Deepwater Production Systems-Final Report; ConocoPhillips Company: Houston, TX, USA, 2005. [Google Scholar]
  41. Picard, D.; Hudson, W.; Bouquier, L.; Dupupet, G.; Zivanovic, I. Composite carbon thermoplastic tubes for deepwater application. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 30 April–3 May 2007. [Google Scholar]
  42. Alexander, C.; Vyvial, B.; Cederberg, C.; Baldwin, D. Evaluating the performance of a composite-reinforced steel drilling riser via full-scale testing for HPHT service. In Proceedings of the 6th International Offshore Pipeline Forum, Houston, TX, USA, 19–20 October 2011. [Google Scholar]
  43. Mintzas, A.; Hatton, S.; Simandjuntak, S.; Little, A.; Zhang, Z. An integrated approach to the design of high performance carbon fibre reinforced risers—From micro to macro—Scale. In Proceedings of the Deep Offshore Technology International Conference, Houston, TX, USA, 6–9 May 2013. [Google Scholar]
  44. Roberts, D.; Hatton, S.A. Development and qualification of end fittings for composite riser pipe. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 6–9 May 2013. [Google Scholar]
  45. American Bureau of Shipping (ABS). Guide for building and classing subsea riser systems. In Design Requirements and Loads; American Bureau of Shipping: Houston, TX, USA, 2008. [Google Scholar]
  46. Det Norske Veritas (DNV). Offshore Standard for Composite Components (DNV-OS-C501); Det Norske Veritas: Oslo, Norway, 2009. [Google Scholar]
  47. Kaw, A.K. Mechanics of Composite Materials, 2nd ed.; CSC Press: Wilmington, DE, USA, 2006. [Google Scholar]
  48. Bedon, C.; Louter, C. Structural glass beams with embedded GFRP, CFRP or steel reinforcement rods: Comparative experimental, analytical and numerical investigations. J. Build. Eng. 2019, 22, 227–241. [Google Scholar] [CrossRef]
  49. American Petroleum Institute (API). Recommended Practice for Planning, Designing and Constructing Fixed Offshore Platforms-Working Stress Design; American Petroleum Institute: Washington, DC, USA, 2000. [Google Scholar]
  50. Blevins, R.D. Flow-Induced Vibratin, 2nd ed.; Krieger Publishing Company: Malabar, FL, USA, 2001. [Google Scholar]
Figure 1. Bilinear kinematic hardening model of (a) X80 steel and (b) titanium.
Figure 1. Bilinear kinematic hardening model of (a) X80 steel and (b) titanium.
Applsci 09 00517 g001
Figure 2. Lamina co-ordinate system.
Figure 2. Lamina co-ordinate system.
Applsci 09 00517 g002
Figure 3. Cross-section for (a) steel riser joint and (b) composite riser joint.
Figure 3. Cross-section for (a) steel riser joint and (b) composite riser joint.
Applsci 09 00517 g003
Figure 4. FE models for (a) flow zone and (b) steel and FRP composite risers.
Figure 4. FE models for (a) flow zone and (b) steel and FRP composite risers.
Applsci 09 00517 g004aApplsci 09 00517 g004b
Figure 5. First three mode shapes of fiber reinforced polymer (FRP) composite riser 3.
Figure 5. First three mode shapes of fiber reinforced polymer (FRP) composite riser 3.
Applsci 09 00517 g005
Figure 6. Time history of displacements in in-line flow direction (X) at three typical locations.
Figure 6. Time history of displacements in in-line flow direction (X) at three typical locations.
Applsci 09 00517 g006aApplsci 09 00517 g006b
Figure 7. Time history of displacements in cross flow direction (Y) at three typical locations.
Figure 7. Time history of displacements in cross flow direction (Y) at three typical locations.
Applsci 09 00517 g007aApplsci 09 00517 g007b
Figure 8. Maximum displacements/outside diameter (Ax/D) in in-line flow direction (X) of all three risers.
Figure 8. Maximum displacements/outside diameter (Ax/D) in in-line flow direction (X) of all three risers.
Applsci 09 00517 g008
Figure 9. Maximum displacements cross flow direction (Y) of all three risers.
Figure 9. Maximum displacements cross flow direction (Y) of all three risers.
Applsci 09 00517 g009
Figure 10. Maximum displacement of the riser 3: (a) at 62.6 s with 0.36 m/s, (b) at 68.0 s with 1.22 m/s and (c) at 54.4 s with 2.13 m/s.
Figure 10. Maximum displacement of the riser 3: (a) at 62.6 s with 0.36 m/s, (b) at 68.0 s with 1.22 m/s and (c) at 54.4 s with 2.13 m/s.
Applsci 09 00517 g010
Figure 11. Maximum von Mises stresses of all three risers.
Figure 11. Maximum von Mises stresses of all three risers.
Applsci 09 00517 g011
Figure 12. Maximum von Mises stresses of the risers (a) riser 1 at 72.6 s with 0.36 m/s, (b) riser 2 at 30.6 s with 0.36 m/s, (c) riser 3 at 24.6 s with 0.36 m/s, (d) riser 1 at 71.0 s with 1.22 m/s, (e) riser 2 at 52.6 s with 1.22 m/s, (f) riser 3 at 68.0 s with 1.22 m/s, (g) riser 1 at 46.6 s with 2.13 m/s, (h) riser 2 at 39.8 s with 2.13 m/s, (i) riser 3 at 54.4 s with 2.13 m/s.
Figure 12. Maximum von Mises stresses of the risers (a) riser 1 at 72.6 s with 0.36 m/s, (b) riser 2 at 30.6 s with 0.36 m/s, (c) riser 3 at 24.6 s with 0.36 m/s, (d) riser 1 at 71.0 s with 1.22 m/s, (e) riser 2 at 52.6 s with 1.22 m/s, (f) riser 3 at 68.0 s with 1.22 m/s, (g) riser 1 at 46.6 s with 2.13 m/s, (h) riser 2 at 39.8 s with 2.13 m/s, (i) riser 3 at 54.4 s with 2.13 m/s.
Applsci 09 00517 g012aApplsci 09 00517 g012b
Figure 13. Factor of safety of composite layers for 0.36 m/s in (a) riser 2-fiber direction, (b) riser 2-transverse direction, (c) riser 3-fiber direction, (d) riser 3-transverse direction and (e) riser 3-shear direction.
Figure 13. Factor of safety of composite layers for 0.36 m/s in (a) riser 2-fiber direction, (b) riser 2-transverse direction, (c) riser 3-fiber direction, (d) riser 3-transverse direction and (e) riser 3-shear direction.
Applsci 09 00517 g013aApplsci 09 00517 g013b
Figure 14. Factor of safety of composite layers for 1.22 m/s in (a) riser 2-fiber direction, (b) riser 2-transverse direction, (c) riser 3-fiber direction, (d) riser 3-transverse direction and (e) riser 3-shear direction.
Figure 14. Factor of safety of composite layers for 1.22 m/s in (a) riser 2-fiber direction, (b) riser 2-transverse direction, (c) riser 3-fiber direction, (d) riser 3-transverse direction and (e) riser 3-shear direction.
Applsci 09 00517 g014aApplsci 09 00517 g014b
Figure 15. Factor of safety of composite layers for 2.13 m/s in (a) riser 2-fiber direction, (b) riser 2-transverse direction, (c) riser 3-fiber direction, (d) riser 3-transverse direction and (e) riser 3-shear direction.
Figure 15. Factor of safety of composite layers for 2.13 m/s in (a) riser 2-fiber direction, (b) riser 2-transverse direction, (c) riser 3-fiber direction, (d) riser 3-transverse direction and (e) riser 3-shear direction.
Applsci 09 00517 g015
Table 1. A summary of previous joint industry projects on composite risers.
Table 1. A summary of previous joint industry projects on composite risers.
YearProject/FunderRiser Type *Materials
1973Ahlstone [35]Dglass fiber with epoxy
1980sInstitut Francais du Petrole and Aerospatiale of France [36]Pglass and carbon fibers with epoxy
1990sNational Institute of Standards and Technology, Advanced Technology Programs [37]P and Dcarbon and E-glass fibers with epoxy
1995–2001Norske Conoco AS and Kvaerner Oilfield Products [38]Dglass and carbon fiber with epoxy
2003ConocoPhillips, Kvaerner Oilfield Products and ChevronTexaco [39,40]P and Dcarbon fiber with epoxy
2007Doris Engineering, Freyssinet, Total and Soficar [41]Pcarbon fiber with thermoplastic PA11
2008–2011Research Partnership to Secure Energy for America [42]Dcarbon fiber with epoxy
2011Magma [43,44]P and Dcarbon fiber with PEEK
* riser type: P—production riser, D—drilling riser.
Table 2. Properties of the materials for the steel riser and liner of composite riser.
Table 2. Properties of the materials for the steel riser and liner of composite riser.
MaterialDensity
(kg/m3)
Modulus (GPa)Tangent
Modulus (GPa)
νYield Stress (MPa)Ultimate Stress (MPa)Elongation at Break (%)
X80 steel78502071.250.35556255.868
Titanium4430113.80.530.34288095014
Table 3. Properties of unidirectional lamina of AS4-Epoxy for composite risers.
Table 3. Properties of unidirectional lamina of AS4-Epoxy for composite risers.
MaterialDensity
(kg/m3)
E1
(GPa)
E2 = E3 (GPa)G12 = G13
(GPa)
ν12 =
ν13
G23
(GPa)
ν23 σ 1 T
(MPa)
σ 1 C
(MPa)
σ 2 T
(MPa)
σ 2 C
(MPa)
τ12 (MPa)
AS4-Epoxy1530135.49.374.960.323.200.461732125649.4167.271.2
Subscripts 1) fiber direction, 2) in-plane transverse and 3) through-thickness direction. Superscript T = Tension, C = Compression.
Table 4. Geometries of risers.
Table 4. Geometries of risers.
RiserI.D.
(m)
O.D.
(m)
Length
(m)
Lay-upst_liner
(mm)
t_0
(mm)
t_±θ
(mm)
t_90
(mm)
±θ
(º)
10.250.325X80 steel
20.250.32925[liner/90/(0/90)10]21.385 2.15
30.250.31125[liner/03/(+53, −53)5/904]21.701.641.7553
Subscript numbers in Lay-ups represent the numbers of lamina with different fiber reinforcement angle.
Table 5. Study cases.
Table 5. Study cases.
Case No.RiserReFlow
Density
(kg/m3)
Velocity (m/s)kinematic Viscosity CoefficientTurbulent Intensity (%)Gravity
(N)
Tension Force (N)Buoyancy (N)
11101,88710240.361.06 × 10−63.7841,53962,30917,734
22111,73610240.361.06 × 10−63.7414,59229,18321,328
33105,62310240.361.06 × 10−63.7711,20022,40019,058
41345,28310241.221.06 × 10−63.2541,53962,30917,734
52378,66010241.221.06 × 10−63.2114,59229,18321,328
63357,94310241.221.06 × 10−63.2311,20022,40019,058
71602,83010242.131.06 × 10−63.0341,53962,30917,734
82661,10410242.131.06 × 10−63.0014,59229,18321,328
93624,93410242.131.06 × 10−63.0211,20022,40019,058
Table 6. Natural frequencies, vortex shedding frequencies and reduced velocities of the risers.
Table 6. Natural frequencies, vortex shedding frequencies and reduced velocities of the risers.
PropertyRiser 1Riser 2 Riser 3
Case 1Case 4Case 7Case 2Case 5Case 8Case 3Case 6Case 9
natural frequency f n /(Hz)2.012.451.89
vortex shedding frequency fs/(Hz)0.240.851.560.220.781.420.230.821.51
reduced velocity Ur0.602.023.530.451.512.640.612.063.61

Share and Cite

MDPI and ACS Style

Wang, C.; Ge, S.; Sun, M.; Jia, Z.; Han, B. Comparative Study of Vortex-Induced Vibration of FRP Composite Risers with Large Length to Diameter Ratio Under Different Environmental Situations. Appl. Sci. 2019, 9, 517. https://doi.org/10.3390/app9030517

AMA Style

Wang C, Ge S, Sun M, Jia Z, Han B. Comparative Study of Vortex-Induced Vibration of FRP Composite Risers with Large Length to Diameter Ratio Under Different Environmental Situations. Applied Sciences. 2019; 9(3):517. https://doi.org/10.3390/app9030517

Chicago/Turabian Style

Wang, Chunguang, Shiquan Ge, Mingyu Sun, Zhirong Jia, and Baomin Han. 2019. "Comparative Study of Vortex-Induced Vibration of FRP Composite Risers with Large Length to Diameter Ratio Under Different Environmental Situations" Applied Sciences 9, no. 3: 517. https://doi.org/10.3390/app9030517

APA Style

Wang, C., Ge, S., Sun, M., Jia, Z., & Han, B. (2019). Comparative Study of Vortex-Induced Vibration of FRP Composite Risers with Large Length to Diameter Ratio Under Different Environmental Situations. Applied Sciences, 9(3), 517. https://doi.org/10.3390/app9030517

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop