Time and frequency domain coupled analysis of deepwater floating production systems

The dynamic analysis of a deepwater floating structure is complicated by the fact that there can be significant coupling between the dynamics of the floating vessel and the attached risers and mooring lines. Furthermore, there are significant nonlinear effects, such as geometric nonlinearities, drag forces, and second order (slow drift) forces on the vessel, and for this reason the governing equations of motion are normally solved in the time domain. This approach is computationally intensive, and the aim of the present work is to develop and validate a more efficient linearized frequency domain approach. To this end, both time and frequency domain models of a coupled vessel/riser/mooring system are developed, which each incorporate both first and second order motions. It is shown that the frequency domain approach yields very good predictions of the system response when benchmarked against the time domain analysis, and the reasons for this are discussed. It is found that the linearization scheme employed for the drag forces on the risers and mooring lines yields a very good estimate of the resulting contribution to slow drift damping.     

Hull/mooring/riser coupled dynamic analysis and sensitivity study of a tanker-based FPSO

A computer program is developed for hull/mooring/riser coupled dynamic analysis of a tanker-based turret-moored FPSO (Floating Production Storage and Offloading) in waves, winds, and currents. In this computer program, the floating body is modeled as a rigid body with six degrees of freedom. The first- and second-order wave forces, added mass, and radiation damping at various yaw angles are calculated from the second-order diffraction/radiation panel program WAMIT. The wind and current forces for various yaw angles of FPSO are modeled following the empirical method suggested by OCIMF (Oil Company International Marine Forum).

The mooring/riser dynamics are modeled using a rod theory and finite element method (FEM), with the governing equations described in a generalized coordinate system. The dynamics of hull, mooring lines, and risers are solved simultaneously at each time step in a combined matrix for the specified connection condition. For illustration, semi-taut chain-steel wire-chain mooring lines and steel catenary risers are employed and their effects on global FPSO hull motions are investigated. To better understand the physics related to the motion characteristics of a turret-moored FPSO, the role of various hydrodynamic contributions is analyzed and assessed including the effects of hull and mooring/riser viscous damping, second-order difference-frequency wave-force quadratic transfer functions, and yaw-angle dependent wave forces and hydrodynamic coefficients. To see the effects of hull and mooring/riser coupling and mooring/riser damping more clearly, the case with no drag forces on those slender members is also investigated. The numerical results are compared with MARIN's wave basin experiments.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

A hybrid time/frequency domain approach for efficient coupled analysis of vessel/mooring/riser dynamics

The dynamic analysis of a deepwater floating structure is complex due to dynamic coupling between the platform and the moorings/risers. Furthermore, the system response at the incident wave frequency and at the resonant low frequency is coupled due to geometric and hydrodynamic nonlinearities. As such, it is generally held that a fully coupled time-domain analysis should be used for an accurate prediction of the dynamic response. However, in a recent work, it is found that for an ultra-deepwater floating system, a fully coupled frequency-domain analysis can provide highly accurate response predictions. One reason is the accuracy of the drag linearization procedure over the motions at two time scales, another is the minimal geometric nonlinearity of the moorings/risers in deepwater. In this paper, the frequency-domain approach is investigated for intermediate water depths, and it is found that the accuracy reduces substantially as geometric nonlinearity becomes important. Therefore, a novel hybrid approach is developed, in which the low-frequency motion is simulated in the time domain while the wave frequency motion is solved in the frequency domain at regular intervals. Coupling between the two analyses is effected by the fact that (i) the low-frequency motion affects the line geometry for the wave frequency motion, and (ii) the wave frequency motion affects the modeling of the drag forces, which damp the low-frequency motion. The method is found to be nearly as accurate as fully coupled time domain analysis even for a system with a preponderance of nonlinear and coupling effects, but requiring only one-tenth of the computational effort.     

Dynamic Analysis of Semi-Submersible Production Platform Under the Failure of Mooring Lines

This paper quantitatively studies the transient dynamic response of a semi-submersible production platform with the loss of one or several positioning mooring lines.A semi-submersible platform,production risers,and positioning mooring lines are all included in the numerical simulation.Increased motion of the semi-submersible platform,tension variation of the remaining mooring lines/risers and the risk of mooring line or riser clashing are all investigated through fully coupled time-domain analysis.Combined environmental loads are selected from irregular waves and the steady current varying from very rough to extreme sea conditions.Three dimension radiation/diffraction theories and Morison’s equation are applied to calculate first-order wave force and second-order mean drift force of floating semi-submersible platform.Nonlinear time-domain finite element methods are employed to analyze the behavior of mooring lines and risers.Results show that the failure of mooring lines seriously reduce the platform’s stability performance.The tension of the rest lines is also increased accordingly.Remaining lines which are closer to the failed lines will have larger tension increase to compensate.Line-Line distance provides practical information for the risk of clashing investigation.     

Time-domain simulation of berthing problem between FPSO and shuttle tanker in waves

A time-domain simulation method based on potential flow model has been developed to investigate the berthing problem between two floating bodies in wave. The boundary value problem is formulated with respect to an earth-fixed coordinate system because the relative positions of the two vessels continuously change during the berthing operation. The classical finite element method is used to solve the Laplace equation in the fluid domain with moving boundary. The linearized free-surface boundary conditions are integrated in time by applying 4th-order Adams–Bashforth–Moulton method. A simple re-mesh algorithm with local and global mesh systems is introduced to update mesh by considering large horizontal movement of the berthing vessel. The developed numerical method is used to investigate the berthing problem between a FPSO and shuttle tanker in waves. The focus is on the wave-induced motion response during the berthing process. The characteristics of the motion responses in berthing operation are examined with various wave frequencies, berthing speeds and wave headings.     

A finite difference approximation for dynamic calculation of vertical free hanging slender risers in re-entry application

The dynamic calculations of slender marine risers, such as Finite Element Method (FEM) or Modal Expansion Solution Method (MESM), are mainly for the slender structures with their both ends hinged to the surface and bottom. However, for the re-entry operation, risers held by vessels are in vertical free hanging state, so the displacement and velocity of lower joint would not be zero. For the model of free hanging flexible marine risers, the paper proposed a Finite Difference Approximation (FDA) method for its dynamic calculation. The riser is divided into a reasonable number of rigid discrete segments. And the dynamic model is established based on simple Euler-Bernoulli Beam Theory concerning tension, shear forces and bending moments at each node along the cylindrical structures, which is extendible for different boundary conditions. The governing equations with specific boundary conditions for riser’s free hanging state are simplified by Keller-box method and solved with Newton iteration algorithm for a stable dynamic solution. The calculation starts when the riser is vertical and still in calm water, and its behavior is obtained along time responding to the lateral forward motion at the top. The dynamic behavior in response to the lateral parametric excitation at the top is also proposed and discussed in this paper.     

The effect of nonlinear multi-contact coupling with gap between risers and guide frames on global spar motion analysis

Nonlinear multi-contact coupling between vertical risers and guide frames inside the spar moon-pool is studied. The existing numerical model for hull/mooring/riser coupled dynamics analysis treats riser as an elastic rod truncated at the keel (truncated riser model), and the risers are free to slide in vertical direction with constant tension, while restricted in horizontal direction. The truncated riser model neglects the portion of the riser inside the moon-pool and tends to overestimate the spar pitch motion. In the newly developed model, the risers are extended through the moon-pool with realistic boundary conditions at multiple guide frames, and thus additional contact forces and moments on the spar hull are considered. The gap effects between the riser buoyancy-cans and riser guide frames are also modeled using three different types of gap springs. Their different dynamic characteristics are extensively studied. The new riser model also considers the Coulomb damping between buoyancy-cans and riser guide frames and also allows the impact-like contact force calculation on risers for ensuring fatigue analysis.     

Nonlinear Dynamic Analysis of Deepwater Drilling Risers Subjected to Random Loads

Excited by ocean currents, random wave and vessel motion, deepwater drilling risers exhibit significant dynamic response. In time domain, a method is proposed to calculate the nonlinear dynamic response of deepwater drilling risers subjected to random wave and dynamic large displacement vessel motion boundary condition. Structural and functional loads, external and internal pressure, free surface effect of irregular wave, hydrodynamic forces induced by current and wave, as well as wave and low frequency (drift) motion of the drilling vessel are all accounted for. An example is presented which illustrates the application of the proposed method. The study shows that long term drift motion of the vessel has profound effect on the envelopes of bending stress and lateral displacement, as well as the range of lower flex joint angle of the deepwater riser. It can also be concluded that vessel motion is the principal dynamic loading of nonlinear dynamic response for the deepwater risers rather than wave force.     

三维物体上二阶波浪绕射的实时模拟

采用二阶时域理论对非线性波浪在任意三维物体周围的绕射问题进行了研究，对自由表面边界条件进行Taylor级数展开，应用摄动展开可以建立相应的边值问题，而且此边值问题的计算域不随时间变化，运用基于B－样条的边界元方法求解每一时刻的波浪场，二阶自由表面边界条件在时间上进行数值积分，在自由表面加了一个人工阻尼层以避免波浪的反射，速度势分解为已知的入射势和未知的散射势，初始条件采用二阶Stokes波浪场，通过加入物体表面边界条件，得到散射势在时间和空间上的发展，本文对圆柱所受规则波的二阶波浪力和波浪爬高进行了计算，数值结果表明此理论计算准确，效率高，数值稳定。     

张力腿式浮式风机耦合动力响应理论模型与数值实现

针对张力腿结构动力响应特点,提出一套张力腿式浮式风机动力耦合理论模型,编制计算程序,并通过与实验结果的对比,验证理论模型以及程序的适用性。张力腿式系泊系统稳定性好,但张力载荷随平台位移的变化显著;张力腿结构的动态响应与支撑平台和风机系统动态响应存在显著的耦合作用。针对张力腿浮式风机动力耦合特性,采用谱方法对张力腿结构建立耦合动力模型;考虑细长体在波浪中受到的粘性力,通过哈密顿原理推导出张力腿运动控制方程,利用自主开发的计算程序DARwind和数值计算方法求解系统动态响应。最后,将数值计算结果同UMAINE的试验数据进行对比分析,以验证程序的合理性。     

Time domain simulation of side-by-side floating bodies using a 3D numerical wave tank approach

The coupled system of two side-by-side fixed and/or floating bodies interacting with a large amplitude nonlinear wave is studied using a direct time domain solution method. The numerical model is based on a three-dimensional mixed Eulerian–Lagrangian (MEL) method under certain simplifying approximations permitting Rankine panel scheme to be implemented over a time-invariant boundary surface to solve the boundary value problem for the unknown velocity potentials. A 4th order Adams–Bashforth–Moulton scheme is used for time marching of rigid-body motion histories of the individual bodies and evolution of the free-surface including the gap region in which large resonant fluid motions occur. A systematic study has been carried out to evaluate the performance of the developed time domain method in simulating the forces and motions as well as the fluid motion in the gap region for the two body system under various arrangements and in different wave-headings. At first, the computed numerical results have been validated and verified with computational and experimental results available in literature for standard geometries such as vertical truncated cylinders and rectangular boxes. Secondly, effectiveness of the damping lid model which is introduced to suppress wave resonance in the gap region is investigated including its influence on maximum sway forces on fixed and floating rectangular barges in side-by-side configurations. Thirdly, comparative studies on absolute and relative motion response for two cases (two rectangular barges, and a FLNG-FPSO + shuttle tanker) in side-by-side arrangement are detailed to bring out the importance of nonlinearities arising due to steep nonlinear incident waves. Finally, coupled motions of the two-body system of an FPSO and a shuttle tanker floating in side-by-side configuration in a steep nonlinear wave field are studied in which the two bodies are connected through hawsers, and also the FPSO is moored to the ground. Additionally there is a fender between the two bodies.     

Nonplanar multi-modal vibrations of fluid-conveying risers under shear cross flows

Risers/pipes conveying fluid are a typical kind of slender structures commonly used in marine engineering. It is of great academic significance and application value for us to evaluate and understand the vibration characteristics and nonlinear responses of these risers under the combined action of internal and external fluid flows. In this paper, the nonplanar vibrations and multi-modal responses of pinned-pinned risers in shear cross flow are numerically studied. With this objective in mind, the van der Pol wake oscillators are used to simulate the dynamical behavior of the vortex shedding in the wake. Two nonlinear equations of motion of the riser are proposed to govern the lateral responses of the riser structure. The nonplanar nonlinear equations for the riser and wake are then discretized by employing Galerkin's method and solved by using a fourth-order Runge–Kutta integration algorithm. Theoretical results show that the coupled frequencies for cross-flow (CF) and in-line (IL) motions and the corresponding coupled damping ratio could be influenced by the external and/or internal fluid velocities. Based on extensive calculations, the dynamical behavior of the riser with various internal and external flow velocities are presented in the form of bifurcation diagrams, time traces, phase portraits, oscillation trajectories and response spectrum curves. It is shown that some interesting dynamical phenomena, such as ‘lock-in’ state, ‘figure-of-eight’ trajectory and quasi-periodic oscillation, could occur in such a fluid-structure interaction system. Our results also demonstrate that the shear parameter can significantly affect the dynamic responses of the riser. When the shear parameter of the cross flow is large, multi-modal quasi-periodic responses of the riser can be excited, showing some new features undetected in the system of fluid-conveying risers in uniform cross flow.     

An Efficient Hydrodynamic Model for Surface Waves

In the present study,a semi-implicit finite difference model for non-bydrostatic,free-surface flows is analyzed and discussed.The governing equations are the three-dimensional free-surface Reynolds-averaged Navier-Stokes equations defined on a general,irregular domain of arbitrary scale.At outflow,a combination of a sponge layer technique and a radiation boundary condition is applied to minimize wave reflection.The equations are solved with the fractional step method where the hydrostatic pressure component is determined first,while the non-hydrostatic component of the pressure is computed from the pressure Poisson equation in which the coefficient matrix is positive definite and symmetric.The advectiou and horizontal viscosity terms are discretized by use of a semi-Lagrangian approach.The resulting model is computationally efficient and unrestricted to the CFL condition.The developed model is verified against analytical solutions and experimental data,with excellent agreement.     

Numerical Prediction of Vortex Induced Vibrations on Top Tensioned Riser in Consideration of Internal Flow

In consideration of the effect of the internal flowing fluid and the external marine environmental condition on the vortex-induced vibration （VIV） of top tensioned riser （Till）, the differential equation is derived based on work-energy principles and the riser near wake dynamics is modeled by Facchinetti＇ s wake oscillator model. Then Galerkin＇ s finite element approximation is implemented to derive the nonlinear matrix equation of the coupled equations and file corresponding numerical programs are compiled which solve the coupled equations directly in the time domain. The comparison of the predicted results with the recent experimental results and the prediction of SHEAR7 is performed. The results show the validity of the proposed method on the prediction of VIV of deep water risers. The effect of internal flow on the dynamic characteristics and dynmnic response of the riser is analyzed and several valuable conelusions are drawn.     

Nonlinear Random Motion Analysis of A Tension Leg Platform Considering the Set-down Motion of A Floating Body

The dynamic analysis of a Tension Leg Platform (TLP) in random wave is investigated by considering the set-down of a floating body. The nonlinear restoring stiffness is derived with the set-down motion of a floating body and the coupled motion of the tension leg and platform and the differential equations of the motion are established. The study focuses on the influence of the set-down motion on the nonlinear response of the platform. By considering different significant wave heights and currents, motion responses of the platform are calculated and compared. The analysis shows that the set-down motion significantly increases the heave motion with low frequency and the equilibrium position of the heave motion with the set-down motion is much lower than that without set-down motion. The results in this paper indicate that the set-down motion has a major impact on the safety of the platform inproduction operation, and it is also a threat to the strength of tension legs and risers.     

Earthquake response of sea-based storage tanks

An effective method for the linear analysis of dynamic response of submerged underwater oil storage tanks resting on a horizontal seabed to horizontal earthquake excitations is presented. The tank is axisymmetric in shape, has a flexible wall/roof, and is filled with oil and water. A general hybrid-finite element solution procedure has been formulated, wherein the tank structure, the interior fluids, as well as the near-field of the exterior water region are discretized into a toroidal finite-element network. The tank displacement is calculated as a superposition of the first few modes of the structure's free vibration. Contribution from the hydrodynamic interaction to the coupled motion is obtained by solving the Laplace equation with the appropriate boundary conditions, which includes a matching to the exterior far-field pressure (analytic) representation to simplify the computation process. The effects of fluids surrounding and inside the tank are studied. It is demonstrated that these effects have, in general, a significant impact on the tank earthquake response analysis.A comprehensive and predictive computer program for use in such tank response analysis has been developed for engineering applications.     

Dynamic Response Analysis of Risers Subjected to Combined Wave and Current Loading

A dynamic response analysis in the frequency domain is presented for risers subjected to combined wave and current loading. Considering the effects of current, a modified wave spectrum is adopted to compute the linearized drag force. An additional drag force convolution term is added to the linearized drag force spectrum, therefore the error is reduced which arises from the truncation of higher order terms in the drag force auto-correlation function. An expression of linearized drag force spectrum is given taking the relative velocity into account. It is found that the additional term is a fold convolution integral. In this paper dynamic responses of risers are investigated, while the influence of floater motion on risers is considered. The results demonstrate that the accuracy of the present method reaches the degree required in time domain analysis.     

预修正快速傅里叶变换高阶边界元方法在波浪对物体作用问题中的应用

边界元方法被广泛应用于波浪对海上婕筑物作用领域,但由于传统边界元方法的存储量和计算量均为未知量的平方量级,很难满足大范围多未知量问题的计算需要.本文基于高阶边界元方法,应用预修正快速傅里叶变换方法,使计算量与存储量分别降低至O(NlogN)和O(N)量级,并得到一个连续的压强分布以适应结构设计的要求,同时可以通过使用满...     

顶部浮体激励下钢悬链线立管的动力响应分析

考虑钢悬链线立管大变形的特性以及内流的影响,利用Kane方法和拉格朗日应变理论建立钢悬链线立管的二维动力学模型。将顶部浮体的激励简化为一个沿顺流向的周期作用力施加于立管顶部,采用平均加速度法求解立管的动力响应。探讨顶部浮体不同激励频率下立管的非线性动力响应以及内流流速的大小对管道共振响应幅值的影响。