Numerical Simulations of Dynamic Pipeline-Vessel Response on a Deepwater SLaying Vessel

The dynamic action induced on offshore pipelines by deepwater S-laying is significant, and directly determines how the pipeline structures are designed and installed. Existing research has not fully investigated the benefits of coupling models of pipeline and pipelaying vessel motions. Therefore, this paper presents a coupled time-domain numerical model for examining the effect of coupled dynamic reactions. The coupled model takes into account the motion of the pipelaying vessel, surface waves, ocean currents, wind forces, pipeline dynamics, and contact between the rollers and the pipeline. A proportional, integral, derivative (PID) controller was used for simulating the control of the pipelaying vessel. The hydrodynamic forces that the pipeline experiences were modeled using the Morison equation. The model was solved using Newmark’s method and verified using OrcaFlex software. The model was then used to analyze practical operations: the laying of a 22” gas export pipeline on the seabed by the pipelaying vessel HYSY201 in the Pingbei-Huangyan gas fields in the East China Sea. The effects of coupled factors on pipelaying vessel motions and pipeline dynamics were approximated. These effects included configurations, axial tensions, and bending moments. The results show a significant connection between the dynamic responses of the pipelines and pipelaying vessel motions.


Introduction
Submarine pipelines are considered to be the most efficient and feasible means of transporting large-scale gas and oil from offshore to onshore.Deepwater pipelaying is a standardized industrial process for offshore pipeline installation and maintenance [1,2].Installations of deepwater pipelines are implemented using the J-laying or S-laying techniques.During J-laying operations, the pipeline is lifted off the stinger in a nearly vertical position, which gives the pipeline a large bending radius from the sea surface to the seabed.Therefore, the touchdown point is relatively close to the pipelaying vessel, which makes it easy to monitor and position.The tension in the pipeline is small.However, the speed and efficiency of J-laying is lower than that of S-laying.Because all the J-laying operations are performed vertically, the vessel is not as stable as in S-laying [3].The S-laying technique still occupies a dominant position in the current pipelaying market because of its specific advantages.S-laying has a higher adaptability and workability in different sea states and at various water depths, and it has a higher pipelaying efficiency and lower costs.In 1967, Plunkett showed that a deepwater pipeline can be assembled in an angular position on a lay barge and maintained at a predetermined high-tensile force while lowering it to the ocean bottom.This process is called S-laying.The S-laying method was developed in the 1980s and became the most frequently-used installation method [4].In the S-laying method, the pipeline is supported using a stinger and passes over a regular roller sequence forming an S-shaped trajectory before landing on the seabed, as illustrated in Figure 1.The upper curved S-shaped parts that rest on stingers are called overbends and are followed by inflection points where the pipeline curvature is zero.Before approaching the seafloor, the pipeline is reversed into a sagbend.Horizontal tension in the pipeline is maintained by mounted tensioner devices on the top of the pipeline.Consequently, S-laying operations depend on roller configuration and stinger radius.Large plastic deformation occurs when pipes pass over stingers.The plastic deformation is caused by the combined effects of the pipelaying vessel motion, roller reaction forces, axial tensions, and bending moments [5,6].
the S-laying method, the pipeline is supported using a stinger and passes over a regular roller sequence forming an S-shaped trajectory before landing on the seabed, as illustrated in Figure 1.The upper curved S-shaped parts that rest on stingers are called overbends and are followed by inflection points where the pipeline curvature is zero.Before approaching the seafloor, the pipeline is reversed into a sagbend.Horizontal tension in the pipeline is maintained by mounted tensioner devices on the top of the pipeline.Consequently, S-laying operations depend on roller configuration and stinger radius.Large plastic deformation occurs when pipes pass over stingers.The plastic deformation is caused by the combined effects of the pipelaying vessel motion, roller reaction forces, axial tensions, and bending moments [5,6].There are several methods for analyzing pipelaying issues.These methods include static analysis, quasi-static analysis, and dynamic analysis.The natural catenary technique involves initially modeling a pipeline with asymptotic expansion, and is ideal for large, nonlinear deflections when tension has more effect on significant length than bending stiffness, and was first used in 1967 by Plunkett [7].Dixon used the stiffened catenary technique to solve the pipeline laying problem, and computed the required pipeline inclination and tension for the pipelaying vessel [8].Brewer examined the impact of pipelaying barge motions on suspended pipelines using quasi-static procedures and an extended stiffened catenary concept [9].In recent years, some researchers have analyzed the S-laying pipeline configuration using the stiffened catenary technique and proposed a co-rotational finite element formulation integrated with Bernoulli's beam theory for assessing offshore pipeline configurations and stresses in S-laying operations [10].
The abovementioned studies on modeling pipelaying issues focused on static and quasi-static analysis.Because offshore production has shifted to deepwater, it is necessary to address operationinduced dynamic pipeline laying factors using modeling.Zan et al. developed a real-time numerical model for dynamically analyzing offshore pipelaying using simulations; they included vessel motion, stratified current flow, nonuniform currents, and other dynamic factors.Their results revealed that dynamic forces were important factors that acted upon the pipeline [11].Gong et al. created an extensive finite element model for S-laying systems that considered the contact between stinger rollers and the pipeline, pipelaying vessel motion, ocean currents, and surface waves using OrcaFlex [12].OrcaFlex is a 3D, nonlinear, time-domain, finite element program developed by Orcina for static and dynamic analysis of marine risers and pipelaying operation systems.However, the calculations in OrcaFlex are based on the response amplitude operator (RAO) motions of the vessel.This implies that the vessel will only make slight movements at the balanced position, and be unaffected by the pipeline.Furthermore, the dynamic positioning system control model is not considered in OrcaFlex There are several methods for analyzing pipelaying issues.These methods include static analysis, quasi-static analysis, and dynamic analysis.The natural catenary technique involves initially modeling a pipeline with asymptotic expansion, and is ideal for large, nonlinear deflections when tension has more effect on significant length than bending stiffness, and was first used in 1967 by Plunkett [7].Dixon used the stiffened catenary technique to solve the pipeline laying problem, and computed the required pipeline inclination and tension for the pipelaying vessel [8].Brewer examined the impact of pipelaying barge motions on suspended pipelines using quasi-static procedures and an extended stiffened catenary concept [9].In recent years, some researchers have analyzed the S-laying pipeline configuration using the stiffened catenary technique and proposed a co-rotational finite element formulation integrated with Bernoulli's beam theory for assessing offshore pipeline configurations and stresses in S-laying operations [10].
The abovementioned studies on modeling pipelaying issues focused on static and quasi-static analysis.Because offshore production has shifted to deepwater, it is necessary to address operationinduced dynamic pipeline laying factors using modeling.Zan et al. developed a real-time numerical model for dynamically analyzing offshore pipelaying using simulations; they included vessel motion, stratified current flow, nonuniform currents, and other dynamic factors.Their results revealed that dynamic forces were important factors that acted upon the pipeline [11].Gong et al. created an extensive finite element model for S-laying systems that considered the contact between stinger rollers and the pipeline, pipelaying vessel motion, ocean currents, and surface waves using OrcaFlex [12].OrcaFlex is a 3D, nonlinear, time-domain, finite element program developed by Orcina for static and dynamic analysis of marine risers and pipelaying operation systems.However, the calculations in OrcaFlex are based on the response amplitude operator (RAO) motions of the vessel.This implies that the vessel will only make slight movements at the balanced position, and be unaffected by the pipeline.Furthermore, the dynamic positioning system control model is not considered in OrcaFlex [13].Liang et al. created a modified finite element model using the software of Abaqus to consider the complex surface contact properties within the overbend section.The contact state properties of various rollers in a single box of rollers were investigated and the emerging support torques (forces) from all roller boxes were computed and contrasted using commercial software [14].Xie determined a pipe's dynamic loading history for S-laying operations using a verified finite element model, and then computed residual plastic deformations for pipe cross-sections after the pipe had reached the seafloor [15].Jensen established a model using nonlinear partial differential equations for the dynamics of pipeline strings suspended from a pipelaying vessel to the seafloor during pipelaying operations.Jensen developed his model from an existing three-dimensional beam model that described tension forces, bending moments, and shearing in sea applications; Jensen added the effects of seabed interactions, hydrodynamic drag, and restoration forces [16].Wittbrodt employed an improved rigid finite element method to enhance computational effectiveness when a dynamic response simulation of pipelines is carried out [17].
However, none of the existing research has fully investigated, or evaluated the benefits of, coupled vessel motion and pipeline dynamics models using dynamic positioning systems.Therefore, the present paper addresses this knowledge gap by investigating a coupled vessel motion-pipeline dynamics time-domain numerical model for analyzing a pipeline's mechanical properties in S-laying operations.The coupled numerical model incorporates, among others, pipelaying vessel motions and the contacts between stinger rollers and the pipeline.A dynamic positioning system was used to control the pipelaying vessel's position within the coupled model.The hydrodynamic forces acting on the pipeline were considered using the Morison equation.Finally, the effects of coupled factors on the dynamic properties of offshore pipelines were extensively investigated using a practical example: the laying of a 22" gas export pipeline on the seafloor by the pipelaying vessel HYSY201 in the Pingbei-Huangyan oil and gas fields in the East China Sea.

S-Laying Vessel Description
The results described are for the pipelaying vessel HYSY201-a deepwater S-laying vessel owned by China's Offshore Oil Engineering Company Ltd. (Beijing, China) (see Figure 2).The main parameters for the pipelaying vessel with a 10% pipe load requirement are presented in Table 1.The pipelaying capacity for double-jointing systems with 6-60" pipes and a 2 × 200-ton tensioning system, supported by a 400-ton Abandonment & Recovery capacity, is provided.The weight of the pipes that can be stored on the upper deck is 9000 tons.Pipe transferring and handling is undertaken using conveyors, rollers, and two mobile gantry deck cranes [18].A permanent truss-type stinger and stinger adjustment system is located at the vessel aft.A dynamic positioning (DP) system has been designed that provides high redundancy and meets DP2 requirements for pipelaying operations in deepwater [19].
Processes 2018, 6, x FOR PEER REVIEW 3 of 22 [13].Liang et al. created a modified finite element model using the software of Abaqus to consider the complex surface contact properties within the overbend section.The contact state properties of various rollers in a single box of rollers were investigated and the emerging support torques (forces) from all roller boxes were computed and contrasted using commercial software [14].Xie determined a pipe's dynamic loading history for S-laying operations using a verified finite element model, and then computed residual plastic deformations for pipe cross-sections after the pipe had reached the seafloor [15].Jensen established a model using nonlinear partial differential equations for the dynamics of pipeline strings suspended from a pipelaying vessel to the seafloor during pipelaying operations.Jensen developed his model from an existing three-dimensional beam model that described tension forces, bending moments, and shearing in sea applications; Jensen added the effects of seabed interactions, hydrodynamic drag, and restoration forces [16].Wittbrodt employed an improved rigid finite element method to enhance computational effectiveness when a dynamic response simulation of pipelines is carried out [17].However, none of the existing research has fully investigated, or evaluated the benefits of, coupled vessel motion and pipeline dynamics models using dynamic positioning systems.Therefore, the present paper addresses this knowledge gap by investigating a coupled vessel motion-pipeline dynamics time-domain numerical model for analyzing a pipeline's mechanical properties in S-laying operations.The coupled numerical model incorporates, among others, pipelaying vessel motions and the contacts between stinger rollers and the pipeline.A dynamic positioning system was used to control the pipelaying vessel's position within the coupled model.The hydrodynamic forces acting on the pipeline were considered using the Morison equation.Finally, the effects of coupled factors on the dynamic properties of offshore pipelines were extensively investigated using a practical example: the laying of a 22″ gas export pipeline on the seafloor by the pipelaying vessel HYSY201 in the Pingbei-Huangyan oil and gas fields in the East China Sea.

S-Laying Vessel Description
The results described are for the pipelaying vessel HYSY201-a deepwater S-laying vessel owned by China's Offshore Oil Engineering Company Ltd. (Beijing, China) (see Figure 2).The main parameters for the pipelaying vessel with a 10% pipe load requirement are presented in Table 1.The pipelaying capacity for double-jointing systems with 6-60″ pipes and a 2 × 200-ton tensioning system, supported by a 400-ton Abandonment & Recovery capacity, is provided.The weight of the pipes that can be stored on the upper deck is 9000 tons.Pipe transferring and handling is undertaken using conveyors, rollers, and two mobile gantry deck cranes [18].A permanent truss-type stinger and stinger adjustment system is located at the vessel aft.A dynamic positioning (DP) system has been designed that provides high redundancy and meets DP2 requirements for pipelaying operations in deepwater [19].

Equations for Coupled Motion
Using Newton's second law, the time-domain equations of motion for the six degrees of freedom for a vessel were linearized to reach an equilibrium point: where M denotes the inertia matrix, x(t) represents the displacement vector, and t represents time.f rad (t) represents the radiation forces arising from the change in fluid momentum caused by the vessel's motion, f hs (t) represents the hydrostatic restoration forces caused by buoyancy and gravity, and f wave (t) represents the wave forces acting on the vessel.In addition, f exc (t) represents other loads: where f c (t) and f w (t) represent the current and wind forces, respectively, acting on the vessel.f dp (t) represents the thruster forces regulated by the DP system, f r (t) represents the reaction forces for the roller-stinger contact model, and f p (t) represents the tensional forces at the top of the pipeline that act on the vessel.
For the approximation of a linear wave, Cummins [20] used potential theory when studying hydrodynamic radiation within the time-domain of an ideal fluid and obtained the following representation.
where A(∞) represents a constant positive-definite matrix known as infinite-frequency added mass, which is linked to the displaced fluid because of the vessel's motion.The second term represents fluid-memory, which captures the transfer of energy from the vessel's motion into the liquid in free-surface radiated waves.The kernel of the convolution term, K(t), is the matrix of retardation or memory functions (impulse responses).The convolution term is known as a fluid-memory model.It represents the fluid memory effects that incorporate energy dissipation caused by the radiated waves generated by the motion of the vessel.This term is not an efficient term to compute numerically because it requires information from the previous time steps, and, in theory, from the start of the body motion.Therefore, most of the existing codes using this formulation truncate the integral in Equation (4): The accurate nature of the convolution term depends on the stored memory time quantity (t 0 ) and the modeling of the vessel's impulse reaction function quality.A pipelaying vessel, which radiates few waves, does not require long-term memory values.
When a vessel's motions are taken into account within the frequency dimension at the expense of other loads, it assumes the form [21] − where A(ω) and B(ω) represent the frequency-dependent added mass alongside the damping matrices.
The term X(jω) represents the motion amplitude and phase caused by wave excitations and F wave (jω) represents the linear forces caused by the waves.Ogilvie determined A(ω) and B(ω) directly by applying a Fourier transform in a sinusoidal regime [22]: Through the application of a Fourier transform, K(t) can be computed from the information for the damping matrices and added mass: In linear theory, when the motion of the pipelaying vessel and waves are assumed to be small, the hydrostatic restoring forces are linear and proportional to the respective displacements of the body, and are simply represented by where C is the hydrostatic restoring matrix.

Wave Forces
Wave forces possess first-order wave forces and second-order slow wave drifting forces.The first-order wave forces arise from incoming waves, thus accounting for the Froude-Krylov (pressure forces caused by undisturbed wave fields) forces as well as the diffraction forces (caused by modified wave fields as a result of their structure).Average low frequency and second-order wave forces play a significant role in floating body dynamics.The horizontal parts are referred to as second-order slow wave drift forces because, under their influence, floating unrestrained vessels would drift away, and undergo a slow drifting motion along the direction of the wave.In this paper, the total first-order and second-order wave forces and moments are computed using the frequency transfer function derived through diffraction analysis and a defined wave spectral density function.
To examine the impact caused by irregular waves, the Joint North Sea Wave Project (Jonswap) spectra are used [23].Jonswap wave spectral density refers to the function of two user-defined coefficients A and B, peak period (T p ), and significant wave heights (H s ).The density of the spectra is expressed as and the exponent q is calculated by Processes 2018, 6, 261 where S(ω) represents the spectral density, A is the generalized Phillips' constant, ω is the circular frequency of the wave component, and ω p is the circular wave frequency at the spectral peak.The Phillips' constant A can be expressed by The parameter t is presented as The Preakness parameter B is defined as To determine the surface wave effects on the dynamic behavior of the pipeline in deepwater S-laying operations, the corresponding Jonswap spectra for all specific wave conditions, four sea states containing varying significant wave heights, and peak periods were considered (Figure 3).represents the spectral density, A is the generalized Phillips' constant, ω is the circular frequency of the wave component, and p ω is the circular wave frequency at the spectral peak.
The Phillips' constant A can be expressed by ( ) The parameter t is presented as The Preakness parameter B is defined as To determine the surface wave effects on the dynamic behavior of the pipeline in deepwater Slaying operations, the corresponding Jonswap spectra for all specific wave conditions, four sea states containing varying significant wave heights, and peak periods were considered (Figure 3).The superposition principle is employed for building irregular sea states.The law of superposition states that problems can be broken down into different subproblems, with each subproblem carrying its solution.The solution for the entire issue is then considered as the total of the subproblem solutions.Therefore, when several wave components with a frequency i  The superposition principle is employed for building irregular sea states.The law of superposition states that problems can be broken down into different subproblems, with each subproblem carrying its solution.The solution for the entire issue is then considered as the total of the subproblem solutions.Therefore, when several wave components with a frequency ω i and wave amplitude ζ ai undergo superposition, the outcome is irregular waves.The waves' overall motion featuring n wave components can be explained by the velocity potential below.
where h represents the depth of water, k i represents the wave number, β i represents the main wave propagation direction, and ε i represents the wave's phase angle.
The amplitude of the wave can be obtained from the spectra S(ω i ) and spreading function D(β i ) of the wave: The first-order wave forces are computed within the time domain as where H (1) (ω) represents the first-order transfer function between wave elevation and excitation force, and was computed using WAMIT [24] software.WAMIT was developed by Lee and Newman, and uses the 3D numerical panel method to solve the linearized hydrodynamic radiation and diffraction problems for the interaction of surface waves with a vessel.The quadratic transfer function (QTF) is used for computing linear wave drift forces: where H (2) (ω) is the desired QTF, and is also computed using WAMIT software.Thus, the resulting wave forces are

Current Forces
The current forces acting upon the vessel were computed during every time step, based on definitions provided by OCIMF (Oil Companies International Marine Forum) [25].The formulae below for the surge, sway, and yaw directions were employed.
where ρ w represents the density of water; L represents the vessel's overall length; T represents the draught of the vessel; V c is the relative velocity between the vessel and the surface current of the sea; and C cx , C cy , and C cmz represent nondimensional current coefficients in a longitudinal direction, lateral direction, and yaw moment direction, and can be obtained through basin modeling or wind tunnel tests.In this study, the current coefficients for HYSY201 are based on data derived from wind tunnel tests (Figure 4).The current forces in the surge, sway, and yaw directions can be calculated.The DP system is designed only to control the vessel motion on the horizontal plane, which includes the surge, sway, and yaw.During pipelaying operations, the vessel motion in the other three directions has little influence on the operations.Thus, the forces in the other three directions can be ignored, and the current forces can be expressed as

Wind Forces
Wind forces acting on bodies were also computed for every time step.A varying or constant wind direction and velocity can be simulated.Different kinds of spectra exist within the model.Alternatively, user-defined time traces can be used.The OCIMF [25] formulae below were employed for the calculation of wind forces in the directions of surge, sway, and yaw.
where w V represents the relative velocity between the vessel and the wind, and wx C , wy C , and wmz C are the desired nondimensional wind coefficients in the longitudinal, lateral, and yaw moment directions.Wind tunnel tests are the most accurate means of estimating wind forces on offshore structures and vessels.The coefficients for HYSY201 are based on data derived from wind tunnel tests, as illustrated in Figure 5.
As in the calculation for current forces, the wind forces in the heave, roll, and pitch directions are ignored.Thus, the wind forces can be expressed as Figure 4.The current coefficients for HYSY201.

Wind Forces
Wind forces acting on bodies were also computed for every time step.A varying or constant wind direction and velocity can be simulated.Different kinds of spectra exist within the model.Alternatively, user-defined time traces can be used.The OCIMF [25] formulae below were employed for the calculation of wind forces in the directions of surge, sway, and yaw.
where V w represents the relative velocity between the vessel and the wind, and C wx , C wy , and C wmz are the desired nondimensional wind coefficients in the longitudinal, lateral, and yaw moment directions.
Wind tunnel tests are the most accurate means of estimating wind forces on offshore structures and vessels.The coefficients for HYSY201 are based on data derived from wind tunnel tests, as illustrated in Figure 5.
As in the calculation for current forces, the wind forces in the heave, roll, and pitch directions are ignored.Thus, the wind forces can be expressed as

Dynamic Positioning System Control Model
Dynamic positioning (DP) is a computer-guided system for automatically maintaining the

Dynamic Positioning System Control Model
Dynamic positioning (DP) is a computer-guided system for automatically maintaining the position and direction of a vessel using the vessel's thrusters and propellers.Position reference sensors coupled with gyro compasses, motion sensors, and wind sensors supply information regarding the direction of environmental forces influencing the vessel's position, magnitude, and position.This paper presents a simulation of the DP time domain to dynamically determine whether pipelaying vessels can maintain their positions.The simulation for the time domain was realized for DP vessel motion using a mathematical model [26].
A PID (proportional-integral-derivative) controller was used for simulating the pipelaying vessel's position.It has been extensively used in engineering and is a stable and simple control system.The controller corrects velocity and position errors, which are then eliminated by the thrusters.The PID's nonlinear horizontal-plane positioning feedback controller is given by [27] where f dp represents the desired control force from the thrusters, ε(t) represents the position error, K P represents the position feedback gains, K I is the desired integral feedback gain, K D represents the velocity feedback gains, f w represents the wind feedforward forces, x(t) is the desired filtered position, and x 0 (t) is the target position.The arrangement and position of the HYSY201 thrusters are illustrated in Figure 6 and summarized in Table 2.

Dynamic Positioning System Control Model
Dynamic positioning (DP) is a computer-guided system for automatically maintaining the position and direction of a vessel using the vessel's thrusters and propellers.Position reference sensors coupled with gyro compasses, motion sensors, and wind sensors supply information regarding the direction of environmental forces influencing the vessel's position, magnitude, and position.This paper presents a simulation of the DP time domain to dynamically determine whether pipelaying vessels can maintain their positions.The simulation for the time domain was realized for DP vessel motion using a mathematical model [26].
A PID (proportional-integral-derivative) controller was used for simulating the pipelaying vessel's position.It has been extensively used in engineering and is a stable and simple control system.The controller corrects velocity and position errors, which are then eliminated by the thrusters.The PID's nonlinear horizontal-plane positioning feedback controller is given by [27]    The vessel's positioning model was separated into low frequency (LF) components and high frequency (HF) components.Thrusters control the LF motions.The DP system is sensitive to high-frequency noise in the velocity and position signals because the noise undergoes amplification and is transported to the thrusters.To address this, the positions and velocities were modified using a Kalman filter.First, motion constructs were inserted using position reference systems and compared with the required position.Previously computed position predictions were then modified using a Kalman gain matrix.The Kalman gain matrix is related to cutoff frequencies, which were selected on the basis of natural durations of motions and on-wave frequency motions [28].

Pipeline Model
The model uses the finite element method and is based on continuum mechanics principles.Lagrangian descriptions were used to describe the pipeline motions.As shown in Figure 7, the particle motion can be represented as where X indicates the particle's position vector and η represents the particle's position at time t.Definition of the displacement vector u is achieved through For the Lagrangian formulation, measurement of strains is undertaken using the Green strain tensor E. The strain tensor is expressed as where dl 0 and dl n represent the line segment PQ's length after, and prior to, deformation, respectively.Additionally, an ideal stress measure is required during the analysis.In most instances, the symmetric Piola-Kirchhoff stress tensor l is used with the Green strain.
frequency noise in the velocity and position signals because the noise undergoes amplification and is transported to the thrusters.To address this, the positions and velocities were modified using a Kalman filter.First, motion constructs were inserted using position reference systems and compared with the required position.Previously computed position predictions were then modified using a Kalman gain matrix.The Kalman gain matrix is related to cutoff frequencies, which were selected on the basis of natural durations of motions and on-wave frequency motions [28].

Pipeline Model
The model uses the finite element method and is based on continuum mechanics principles.Lagrangian descriptions were used to describe the pipeline motions.As shown in Figure 7, the particle motion can be represented as ( ) , t = η η X (26) where X indicates the particle's position vector and η represents the particle's position at time t.Definition of the displacement vector u is achieved through = + η X u (27) For the Lagrangian formulation, measurement of strains is undertaken using the Green strain tensor E. The strain tensor is expressed as (28) where 0 dl and n dl represent the line segment PQ's length after, and prior to, deformation, respectively.Additionally, an ideal stress measure is required during the analysis.In most instances, the symmetric Piola-Kirchhoff stress tensor l is used with the Green strain.Malvern [29] is referred to for the detailed theory of the pipeline model.The virtual work equation can be used for expressing a finite body's equilibrium.Using Green strains and Piola-Kirchhoff stresses, the equation can be expressed as where the colon (:) notation indicates an inner product, δ represents virtual quantities, and V 0 and A 0 represent the volume and surface of the first reference configuration, respectively.Body forces f 0 and surface traction t 0 represent the unit volume and unit surface within the initial state of reference.Consequently, the equation for the dynamic equilibrium in terms of virtual work can be expressed as where ρ 0 represents mass density and c indicates a function of viscous damping density.Damping forces are proportional to velocity.The forces of inertia are proportional to the structure's acceleration and mass.The finite element nodal points may have up to six degrees of freedom (three translations and three rotations).Finite rigid bodies can be introduced at the nodes.Thus, the tensional forces at the top of the pipeline f p (t) can be obtained [30].

Hydrodynamic Forces Acting on the Pipeline
Hydrodynamic in-line forces from currents acting on the pipeline apply to an accelerated fluid setting in which the pipeline is kept vertical and stationary: Equation ( 31) is known as the Morison equation [31].In this equation, ρ w represents the density of water, D is the outer diameter of the pipeline, A represents the cross-sectional area of the pipeline, u represents the velocity of the fluid, C d represents the drag coefficient, and C m represents the hydrodynamic mass.The drag and mass coefficients depend on several parameters, including the ratio of surface roughness, a relative current value, the Keulegan-Carpenter number, and the Reynolds number.The drag and mass coefficients should be determined empirically [32].
When the pipeline is in motion and oceans currents are present during pipelaying operations, Equation ( 31) can be simplified for hydrodynamic force per unit length as where v r =v −v c represents relative velocity, v c represents the velocity of the ocean current, and v represents the pipeline's velocity.It should be noted that the Froude-Krylov force is independent of v because the force is linked to the fluid's absolute motion.It is assumed that ocean currents vary gradually; this implies that .v c ≈ 0 and the Froude-Krylov term can be discarded.

Roller: Simplifying Assumptions
On pipelaying vessels, the pipeline is usually supplied by multiple discrete rollers.For example, the stinger of HYSY201 has three sections: a tail section, an intermediate section, and a hitch section, with lengths of 19.6 m, 22.9 m, and 33.2 m, respectively.Ten boxes of regularly spaced rollers are on the stingers, as illustrated in Figure 8a; it can be seen that the pipeline is in contact with the stinger rollers as shown in Figure 8b.
In practice, the vertical heights of all roller boxes could be modified slightly to meet the pipelaying conditions prior to offshore pipeline installation.The roller boxes support the pipeline on the stinger, and pipelaying barges prevent the pipelines' downward lateral and vertical displacement.All rollers within the roller boxes are analyzed to ascertain if the pipelines are in contact with the rollers, and then support requirements are computed and applied to the rollers and pipes.In this paper, it is assumed that the rollers are cylindrical, and that the support acts in directions that are mutually perpendicular towards the pipelines' longitudinal axes and all supporting rollers.The contact forces of the roller acting on the pipeline can be calculated by linear stiffness k s [33]: where D represents the pipe diameter and δ represents the distance between two straight lines on different surfaces.In practice, the vertical heights of all roller boxes could be modified slightly to meet the pipelaying conditions prior to offshore pipeline installation.The roller boxes support the pipeline on the stinger, and pipelaying barges prevent the pipelines' downward lateral and vertical displacement.All rollers within the roller boxes are analyzed to ascertain if the pipelines are in contact with the rollers, and then support requirements are computed and applied to the rollers and pipes.In this paper, it is assumed that the rollers are cylindrical, and that the support acts in directions that are mutually perpendicular towards the pipelines' longitudinal axes and all supporting rollers.The contact forces of the roller acting on the pipeline can be calculated by linear stiffness ks [33]: where D represents the pipe diameter and δ represents the distance between two straight lines on different surfaces.

Model Solution
Newmark's integral method was used to solve Equation ( 1), in which the velocity and displacement of the vessel-pipeline system during the instant t + ∆t(i + 1) can be expressed as where S is the matrix of the pipelaying vessel displacement for the six degree freedom in surge, sway, heave, roll, pitch, and yaw directions; ∆t represents the time step length; and β and δ represent parameters associated with stability and accuracy, respectively.Using Equation ( 34), the acceleration and velocity at the instant i + 1 is approximated using the displacements at the instant i + 1 and the instant i: Substitution of Equation (35) into Equation ( 1) and simplification produces the equation ..

S i (36)
f i+1 represents all the forces input into the system, including environmental forces (f wave , f w , and f c ), radius forces (f rad ), hydrostatic forces (f hs ), tensional forces (f p ) at the top of the pipeline that act on the vessel, and roller-stinger contact forces (f r ).The responses of the system, including accelerations, velocities, and displacements at the instant i + 1 were derived by substituting Equation (36) into Equation ( 35) and simplifying.Then, these results were input into the DP system to be compared with the target position to obtain the deviation.Based on the deviation, the thrusters produce specific forces and feed them back to the Newmark method solver to run the calculation in the next step.Thus, the vessel will be maneuvered, step by step, to the target position.The solution procedure for the coupled model is shown in Figure 9.
where S is the matrix of the pipelaying vessel displacement for the six degree freedom in surge, sway, heave, roll, pitch, and yaw directions; Δt represents the time step length; and β and δ represent parameters associated with stability and accuracy, respectively.Using Equation ( 34), the acceleration and velocity at the instant i + 1 is approximated using the displacements at the instant i + 1 and the instant i:  1) and simplification produces the equation represents all the forces input into the system, including environmental forces ( wave f , w f , and c f ), radius forces ( rad f ), hydrostatic forces ( hs f ), tensional forces ( p f ) at the top of the pipeline that act on the vessel, and roller-stinger contact forces ( r f ).The responses of the system, including accelerations, velocities, and displacements at the instant i + 1 were derived by substituting Equation (36) into Equation (35) and simplifying.Then, these results were input into the DP system to be compared with the target position to obtain the deviation.Based on the deviation, the thrusters produce specific forces and feed them back to the Newmark method solver to run the calculation in the next step.Thus, the vessel will be maneuvered, step by step, to the target position.The solution procedure for the coupled model is shown in Figure 9.

Verification and Comparison
The model data originate from a practical engineering program for the Pingbei-Huangyan oil and gas fields, located in the East China Sea, the People's Republic of China.The oil and gas fields are approximately 425 km southeast of Shanghai.The pipeline parameters for the 22" gas pipeline are shown in Table 3.A fixed stinger was placed at the stern of the HYSY201 for launching the pipe into the water with a suitable curvature (radius of 200 m).A water depth of 110 m and a soft seabed with a normal soil stiffness of 10 5 kN/m 2 were assumed for the purpose of the analysis.
The overall configurations, axial tensions, and bending moments of the pipeline from HYSY201 through the stingers to the seafloor were computed using OrcaFlex [13] and the current model.The dynamic positioning system was not included while computing the motions of the vessel.It can be seen from Figure 10 that the results of the pipeline calculated using the OrcaFlex software and the present model are in good agreement.The offshore pipeline configurations are in good agreement.Similarly, the axial tensions on the pipeline are nearly identical.The axial tension at the top of the pipeline reaches 1117.5 kN.Additionally, the bending moments match independently throughout the whole pipeline, especially for the bending moments on the stinger and the pipelaying vessel; the variations in bending moments caused by the roller supports are highly coincident.The overall configurations, axial tensions, and bending moments of the pipeline from HYSY201 through the stingers to the seafloor were computed using OrcaFlex [13] and the current model.The dynamic positioning system was not included while computing the motions of the vessel.It can be seen from Figure 10 that the results of the pipeline calculated using the OrcaFlex software and the present model are in good agreement.The offshore pipeline configurations are in good agreement.Similarly, the axial tensions on the pipeline are nearly identical.The axial tension at the top of the pipeline reaches 1117.5 kN.Additionally, the bending moments match independently throughout the whole pipeline, especially for the bending moments on the stinger and the pipelaying vessel; the variations in bending moments caused by the roller supports are highly coincident.

Coupled Dynamic Pipelaying Analysis Results
The velocity values of the currents within the oil and gas fields of Pingbei-Huangyan were obtained from field measurements.Figure 11 illustrates the differences in current velocity for various annual return periods and water depths.The current velocity data for 1-year return periods at varying water depths were used in the dynamic analysis.The current direction was assumed to be

Coupled Dynamic Pipelaying Analysis Results
The velocity values of the currents within the oil and gas fields of Pingbei-Huangyan were obtained from field measurements.Figure 11 illustrates the differences in current velocity for various annual return periods and water depths.The current velocity data for 1-year return periods at varying water depths were used in the dynamic analysis.The current direction was assumed to be coincident with the heading of the pipelaying vessel.
obtained from field measurements.Figure 11 illustrates the differences in current velocity for various annual return periods and water depths.The current velocity data for 1-year return periods at varying water depths were used in the dynamic analysis.The current direction was assumed to be coincident with the heading of the pipelaying vessel.
Table 4 shows the environmental parameters for the hydrodynamic basin test for significant wave height (Hs), peak period peak (Tp), enhancement factor ( γ ), mean wind speed (Vw), and wave direction ( α ).The API (American Petroleum Institute) spectrum model was used to calculate the wind speed acting on the vessel.

Hs(m) Tp(s) γ
Vw(m/s) α (°)  Table 4 shows the environmental parameters for the hydrodynamic test for significant wave height (H s ), peak period peak (T p ), enhancement factor (γ), mean wind speed (V w ), and wave direction (α).The API (American Petroleum Institute) spectrum model was used to calculate the wind speed acting on the vessel.The hydrodynamic parameters used in the analysis, the drag coefficient, and the added mass coefficient are presented in Table 5 [34].
Table 5. Drag and added mass coefficients for the pipeline dynamic analysis.

Drag Coefficient 1.2 Added Mass Coefficient 1
The initial vessel position was determined (97.3 m, −1.52 m) and the heading angle was 1.22 • .The PID controller was switched on and propelled the vessel to the desired surface reference position (98.5 m, 0 m) with a specified heading angle of 0 • .The simulation time was 1000 s and the simulation step size was 0.1 s. Figure 12 shows the random wave elevation time history curve in 1000 s, which contains the maximum wave height (2.668 m) at the global time of 497.8 s.
were performed with identical initial conditions.Six degrees of freedom for the vessel motion are presented in Figure 13.
As shown in Figure 13a for surge, the vessel quickly reaches the target position of 98.5 m from the initial position 97.3 m in the uncoupled result.However, in the coupled result, the process takes longer, and the vessel goes through a larger deviation to 99.91 m before it approaches the target position of 98.5 m; this is because of the effect of pipeline forces.A comparison between coupling and uncoupling for the sway direction is shown in Figure 13b.It can be seen that in both uncoupled and coupled results, the vessel approaches the target position 0 m in a short time, and that sway is slightly less than 1 m around the target position.Furthermore, the dynamic motion that takes coupling into account is slightly smaller than the dynamic motion that does not consider coupling.
Figure 13c,d show that the time histories of coupled results in the heave and roll directions closely follow the uncoupled results.However, the coupled result in the pitch direction is smaller than the uncoupled result, as shown in Figure 13e.This is due to the influence of the pipeline on the vessel stern.
Figure 13f shows that for the yaw direction, the vessel gets to the target position of 0° from its initial position 1.22° in a short time in the uncoupled result, while it takes longer in the coupled result to get to the same position.In addition, the deviation in the coupled result reaches −1.44 m as the vessel approaches the target position.This is much larger than for the uncoupled result.To investigate the S-laying pipeline and DP coupling effect on the results of thruster force allocation, the coupled and decoupled simulations of the dynamic response of the pipelaying system were performed with identical initial conditions.Six degrees of freedom for the vessel motion are presented in Figure 13.
As shown in Figure 13a for surge, the vessel quickly reaches the target position of 98.5 m from the initial position 97.3 m in the uncoupled result.However, in the coupled result, the process takes longer, and the vessel goes through a larger deviation to 99.91 m before it approaches the target position of 98.5 m; this is because of the effect of pipeline forces.
A comparison between coupling and uncoupling for the sway direction is shown in Figure 13b.It can be seen that in both uncoupled and coupled results, the vessel approaches the target position 0 m in a short time, and that sway is slightly less than 1 m around the target position.Furthermore, the dynamic motion that takes coupling into account is slightly smaller than the dynamic motion that does not consider coupling.
Figure 13c,d show that the time histories of coupled results in the heave and roll directions closely follow the uncoupled results.However, the coupled result in the pitch direction is smaller than the uncoupled result, as shown in Figure 13e.This is due to the influence of the pipeline on the vessel stern.
Figure 13f shows that for the yaw direction, the vessel gets to the target position of 0 • from its initial position 1.22 • in a short time in the uncoupled result, while it takes longer in the coupled result to get to the same position.In addition, the deviation in the coupled result reaches −1.44 m as the vessel approaches the target position.This is much larger than for the uncoupled result.
As shown in Figure 14, the thrust force required under coupling conditions is much greater than that under uncoupling conditions.This is especially true for the surge force shown in Figure 14a, where only a small thrust force is required by the DP system to maintain the thruster's position when the coupling effect is not considered.When the coupling effect is considered, the total thrust force required increases to 1590 kN.It can be seen from Figure 14b that the difference in thrust force for coupling and without coupling is small in the sway direction.This small difference is the result of the pipeline being in the same direction as the pipelaying vessel.Thus, coupling has less impact on the sway direction force than on the surge direction force.To maintain the pipelaying vessel's heading and position, the torsional moment in the heading direction of the pipelaying vessel when coupling is considered is 1.3 times the torsional moment when coupling is not considered, as illustrated in Figure 14c.
The time histories for the pipeline top tension force are presented in Figure 15.The dynamic maximum pipeline tension is 1556 kN at the global time of 128 s, which is 39.24% higher than the static tension force (1117.5 kN).The dynamic minimum pipe tension is 880.5 kN at the global time of 488 s, which is about 29.94% smaller than the static tension force (1117.5 kN).Therefore, the designed tensioner capability has to be significantly higher than the static results, and the pipeline dynamic tension should be considered when planning pipeline installations.As shown in Figure 14, the thrust force required under coupling conditions is much greater than that under uncoupling conditions.This is especially true for the surge force shown in Figure 14a, where only a small thrust force is required by the DP system to maintain the thruster's position when the coupling effect is not considered.When the coupling effect is considered, the total thrust force required increases to 1590 kN.It can be seen from Figure 14b that the difference in thrust force for coupling and without coupling is small in the sway direction.This small difference is the result of the pipeline being in the same direction as the pipelaying vessel.Thus, coupling has less impact on the sway direction force than on the surge direction force.To maintain the pipelaying vessel's heading and position, the torsional moment in the heading direction of the pipelaying vessel when coupling is considered is 1.3 times the torsional moment when coupling is not considered, as illustrated in Figure 14c.As shown in Figure 14, the thrust force required under coupling conditions is much greater than that under uncoupling conditions.This is especially true for the surge force shown in Figure 14a, where only a small thrust force is required by the DP system to maintain the thruster's position when the coupling effect is not considered.When the coupling effect is considered, the total thrust force required increases to 1590 kN.It can be seen from Figure 14b that the difference in thrust force for coupling and without coupling is small in the sway direction.This small difference is the result of the pipeline being in the same direction as the pipelaying vessel.Thus, coupling has less impact on the sway direction force than on the surge direction force.To maintain the pipelaying vessel's heading and position, the torsional moment in the heading direction of the pipelaying vessel when coupling is considered is 1.3 times the torsional moment when coupling is not considered, as illustrated in Figure 14c.The time histories for the pipeline top tension force are presented in Figure 15.The dynamic maximum pipeline tension is 1556 kN at the global time of 128 s, which is 39.24% higher than the static tension force (1117.5 kN).The dynamic minimum pipe tension is 880.5 kN at the global time of 488 s, which is about 29.94% smaller than the static tension force (1117.5 kN).Therefore, the designed tensioner capability has to be significantly higher than the static results, and the pipeline dynamic tension should be considered when planning pipeline installations.area, where the arc length of the pipeline is approximately 427-530 m, exhibits almost no change.The bending moment at the upper curved part, where the arc length of the pipeline is approximately 390-427 m, has a small change when compared with the static bending moment.Then, the bending moment increases gradually at the lower curved part, where the arc length of the pipeline is approximately 50-390 m.These results indicate that the movement of the pipelaying vessel has an influence on the bending moment of the pipeline and a strong influence on the axial tension force.Figures 16 and 17 show the minimum and maximum values of the tension force and bending moment of the pipeline when the pipelaying vessel behaves stably for the period of 700 to 1000 s.As shown in Figure 16, the static axial tension force differs significantly from the maximum and minimum tension forces of the pipeline.The maximum axial tension force of the pipeline is 33.76% larger than the static axial tension force, and the minimum axial tension force of the pipeline is 23.88% smaller than the static tension force.The bending moment of the pipeline (Figure 17) in the stinger area, where the arc length of the pipeline is approximately 427-530 m, exhibits almost no change.The bending moment at the upper curved part, where the arc length of the pipeline is approximately 390-427 m, has a small change when compared with the static bending moment.Then, the bending moment increases gradually at the lower curved part, where the arc length of the pipeline is approximately 50-390 m.These results indicate that the movement of the pipelaying vessel has an influence on the bending moment of the pipeline and a strong influence on the axial tension force.

Conclusions
This paper presents a coupled dynamic model of a pipelaying vessel and the pipeline.The model takes into account the combined effects of wave, wind, current, and thruster forces to compute the dynamic motion of the vessel.In addition, current forces and roller contact forces with the pipeline were incorporated to calculate the pipeline's dynamic configurations, tension forces, and bending moments.The results from the model were verified using results from OrcaFlex software.A demonstration of the model was provided by assessing the placement of 22″ gas export pipelines on the seafloor by the HYSY201 pipelaying vessel in the Pingbei-Huangyan oil and gas fields in the East China Sea.The assessment revealed a significant relationship between the dynamic responses of the pipeline and the pipelaying vessel's motion.This shows that coupled disturbance models, which can compute the coupled motion of the pipeline and vessel are necessary.The dynamic positioning system should be included while computing the motions of the vessel.The most important influence of the pipeline feedback on the vessel is in the surge and the yaw directions.The total thrusting force increased significantly in the coupled model.The coupled model can be of value in planning offshore pipelaying activities and installation.

Conclusions
This paper presents a coupled dynamic model of a pipelaying vessel and the pipeline.The model takes into account the combined effects of wave, wind, current, and thruster forces to compute the dynamic motion of the vessel.In addition, current forces and roller contact forces with the pipeline were incorporated to calculate the pipeline's dynamic configurations, tension forces, and bending moments.The results from the model were verified using results from OrcaFlex software.A demonstration of the model was provided by assessing the placement of 22 gas export pipelines on the seafloor by the HYSY201 pipelaying vessel in the Pingbei-Huangyan oil and gas fields in the East China Sea.The assessment revealed a significant relationship between the dynamic responses of the pipeline and the pipelaying vessel's motion.This shows that coupled disturbance models, which can compute the coupled motion of the pipeline and vessel are necessary.The dynamic positioning system should be included while computing the motions of the vessel.The most important influence of the pipeline feedback on the vessel is in the surge and the yaw directions.The total thrusting force increased significantly in the coupled model.The coupled model can be of value in planning offshore pipelaying activities and installation.
ω and wave amplitude ai ζ undergo superposition, the outcome is irregular waves.The waves' overall motion featuring n wave components can be explained by the velocity potential below.S (m 2 s/rad)

Figure 6 .Table 2 .
Figure 6.The arrangement and position of the HYSY201 thrusters.

Figure 6 .
Figure 6.The arrangement and position of the HYSY201 thrusters.

Figure 7 .
Figure 7.The particle motion during deformation.Figure 7. The particle motion during deformation.

Figure 7 .
Figure 7.The particle motion during deformation.Figure 7. The particle motion during deformation.

Figure 8 .
Figure 8. Stinger for S-laying operations: (a) articulated stinger of the HYSY201 and (b) the pipeline in contact with the stinger roller.

Figure 8 .
Figure 8. Stinger for S-laying operations: (a) articulated stinger of the HYSY201 and (b) the pipeline in contact with the stinger roller.

Figure 9 .
Figure 9. Solution procedure for the coupled model.

Figure 9 .
Figure 9. Solution procedure for the coupled model.

Figure 10 .
Figure 10.Comparison between the present model results and OrcaFlex simulations: (a) configurations; (b) tension forces; and (c) bending moments.

Figure 10 .
Figure 10.Comparison between the present model results and OrcaFlex simulations: (a) configurations; (b) tension forces; and (c) bending moments.

Figure 11 .
Figure 11.The variation of the current velocity at different water depths.Figure 11.The variation of the current velocity at different water depths.

Figure 11 .
Figure 11.The variation of the current velocity at different water depths.Figure 11.The variation of the current velocity at different water depths.

Figure 12 .
Figure 12.The time history of the wave elevation.

Figure 12 .
Figure 12.The time history of the wave elevation.

Figure 14 .
Figure 14.Dynamic total force thrust for uncoupled and coupled results: (a) Surge direction force; (b) Sway direction force; and (c) Yaw direction moment.

Figure 14 .
Figure 14.Dynamic total force thrust for uncoupled and coupled results: (a) Surge direction force; (b) Sway direction force; and (c) Yaw direction moment.

Figure 15 .Figure 15 .
Figure 15.Time histories for the pipeline top tension force.

Table 2 .
The position of the thrusters for HYSY201.

Table 3 .
Parameters for the 22 gas pipeline.

Table 4 .
Environmental parameters for the hydrodynamic basin test.

Table 4 .
Environmental parameters for the hydrodynamic basin test.