海洋柔性结构涡激振动的流固耦合机理和响应

对近几十年来国内外在涡激振动的基础研究包括机理认识和动响应分析等方面的进展进行了论述,尤其针对海洋油气平台中的立管、隔水管等细长柔性结构的涡激振动.描述了涡激振动这种典型的非线性流固耦合现象所具有的特征,包括自激、自限制、展向相关、尾迹水动力与结构动力的流固耦合等及其主要影响参数.介绍了目前常用的结构响应预测方法和相关实验.通过讨论当前理论研究和实际工程中的热点问题,诸如多模态宽带振动、浮体运动与水下立管的耦合、响应抑制措施、双向振动、高雷诺数下的大尺度物理实验等,对今后该领域的研究方向进行了力所能及的展望.     

Fluid-solid coupling and dynamic response of vortex-induced vibration of slender ocean cylinders

对近几十年来国内外在涡激振动的基础研究包括机理认识和动响应分析等方面的进展进行了论述,尤其针对海洋油气平台中的立管、隔水管等细长柔性结构的涡激振动.描述了涡激振动这种典型的非线性流固耦合现象所具有的特征,包括自激、自限制、展向相关、尾迹水动力与结构动力的流固耦合等及其主要影响参数.介绍了目前常用的结构响应预测方法和相关实验.通过讨论当前理论研究和实际工程中的热点问题,诸如多模态宽带振动、浮体运动与水下立管的耦合、响应抑制措施、双向振动、高雷诺数下的大尺度物理实验等,对今后该领域的研究方向进行了力所能及的展望.     

剪切流场中含内流立管横向涡激振动特性

立管是海洋工程中输送油气或其他矿产资源的必备结构, 外部洋流引起的立管涡激振动影响着立管的疲劳寿命, 危害深海资源开发. 本文基于欧拉?伯努利梁方程, 结合半经验时域水动力模型, 建立剪切流与内流耦合作用下海洋立管涡激振动预报模型, 运用有限元方法和Newmark-β逐步积分法求解方程, 首先将数值模拟结果与实验数据进行对比, 验证模型正确性. 然后, 运用此模型, 对剪切流作用下含内流的顶张立管在不同内流速度和密度下的横向涡激振动响应特性进行研究, 主要分析了立管的横向振动模态、振动频率以及均方根位移等涡激振动参数随内流速度和密度等参数的变化规律. 结果表明, 在剪切流场中, 含内流海洋立管在横向上表现出多模态多频率的涡激振动;立管横向振动的最大均方根位移随内流速度和密度的增大而增大, 特别是当内流速度较大时, 横向最大均方根位移增大明显;立管横向振动的主导频率随内流速度和密度的增大而减小, 并且内流密度的增大同样会引起模态转换和频率转换.      

基于能量传递规律的海洋立管涡激振动抑制研究

涡激振动是造成海洋立管疲劳损伤的重要因素, 抑制振动能够保障结构安全, 延长使用寿命. 多数涡激振动抑制方法基于干扰流场的方式, 但在复杂环境条件下, 仅通过干扰流场对振动的抑制效果有限. 因此, 从结构层面考虑开展了海洋立管涡激振动抑制研究. 基于能量传递的理论, 阐述了立管涡激振动过程中的能量传递规律. 振动能量以行波形式由能量输入区传播至能量耗散区, 主要在能量耗散区被消耗. 通过局部增大能量耗散区的阻尼, 增加振动能量在传播过程中的消耗, 实现涡激振动抑制. 为了求解立管涡激振动响应, 构建了尾流振子预报模型, 并根据实验结果验证了理论模型的可靠性. 基于理论计算得到的能量系数, 判定立管涡激振动的能量输入区和能量耗散区. 通过对比立管增大阻尼前后的响应, 分析了涡激振动抑制效果. 研究结果表明: 在能量输入区增大阻尼对涡激振动的抑制效果并不显著; 在能量耗散区增大阻尼使能量衰减系数达到临界值之后, 能够显著降低立管上部和底部的涡激振动位移; 当能量衰减系数超过临界值后, 继续增大耗散区阻尼对涡激振动抑制效果的提升不明显.      

海洋热塑性增强管(RTP)涡激振动数值计算

基于Van der Pol尾流振子模型和多体系统传递矩阵法(transfer matrix method for multibody systems, MSTMM), 建立了可以快速预测海洋热塑性增强管(reinforced thermoplastic pipe, RTP)振动特性和涡激振动响应的动力学模型. 仿真结果与ANSYS软件仿真结果以及文献实验数据对比, 验证了本文模型的准确性. 研究了考虑RTP立管刚性接头, 不同顶张力, 不同来流分布等情况对RTP立管涡激振动响应的影响. 计算结果表明: 流速越大, 立管涡激振动激发出的模态越高; 立管涡激振动主要受低阶模态控制; 立管的刚性接头对立管的湿模态影响较小, 但是对较高阶模态为主所激发出的涡激振动振幅分布影响较大; 剪切流对沿立管轴向的涡激振动振幅分布影响较大, 低流速能量小所引起的涡激振动幅值较小, 但是当剪切流流速达到能激发出较高阶模态时, 相比同等流速的均匀流所引起的涡激振动振幅要大.      

多用途输液立管涡激振动实验研究与疲劳分析

为研究输运不同流体的海洋立管在海流作用下的振动规律,在大型波浪流水槽中进行涡激振动模型实验。实验分别将四种不同质量比的立管模型竖直固定于支架上,立管外部承受不同速度的流体作用,上端施加顶张力。立管模型上均匀布置六个测点,根据每个测点布置的两个应变计,分别测得来流向和横向两个方向振动响应。通过小波变换对实验数据进行去噪处理,利用振型分解法求解立管各点涡激振动位移。考察输运不同流体对立管自振频率以及涡激振动响应的影响,并利用雨流计数法对模型进行疲劳分析。实验结果表明,随质量比增加立管涡激振动频率降低;低质量比的立管更容易产生大位移。     

倾斜圆柱结构涡激振动研究进展

圆柱结构涡激振动现象在生活中十分常见, 如海洋工程中的管道、土木工程中的高耸建筑、桥梁斜拉索, 核工程中的热交换器等频繁受到涡激振动影响, 诱发结构的疲劳损伤, 甚至破坏失效. 现阶段, 人们对垂直来流作用下圆柱结构涡激振动机理已有较为全面的认识. 然而, 当圆柱倾斜置于流场中, 结构后缘的尾流形态与垂直放置差异显著, 结构与流体的耦合作用机理更为复杂. 为简化倾斜圆柱涡激振动问题, 提出了不相关原则, 来流速度被分解为垂直圆柱结构轴向和平行圆柱结构轴向的两个速度分量, 仅考虑垂直结构轴向速度分量的影响, 忽略平行结构轴向速度分量的影响. 近年来, 针对倾斜圆柱涡激振动及不相关原则的适用性, 出现了大量实验和数值模拟研究成果. 为了深化对倾斜圆柱结构涡激振动相关机理的认知, 本文全面阐述了倾斜圆柱结构涡激振动响应规律、尾迹流场模式和流体力特性等方面的研究进展, 分析了不相关原则的适用范围, 探讨了倾斜圆柱结构涡激振动抑制措施, 并对今后该领域的研究进行了力所能及的展望.      

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法. 鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究. 先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解. 对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析. 分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成. 当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导. 当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导. 随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.      

基于非线性吸能机理的涡激振动减振理论与实验研究

流致振动现象广泛存在于机械、航空、土木和石油等重要工程领域, 为防止工程结构因流致振动行为而造成疲劳破坏, 有必要对稳定性、动力学响应及其振动控制做深入研究. 本文提出了一种由弹簧和质量块构成的非线性吸能器(nonlinear targeted energy transfer, NTET), 研究了该非线性吸能器对弹性支承圆柱体涡激振动的被动控制影响机制. 基于能量法推导了圆柱体涡激振动非线性被动控制的耦合动力学方程, 通过设计非线性弹簧?质量块构型的NTET, 进一步开展了涡激振动控制的实验研究, 并与理论预测结果进行了较好的对比, 获得提升涡激振动控制效果的最佳参数值. 研究发现, NTET的质量、弹簧刚度以及弹簧预应力等参数会对涡激振动控制效果产生显著的影响. 本文研究结果表明, 该耦合系统中圆柱体和NTET均表现出周期性的稳态振动响应, NTET质量的改变会显著影响系统的耦合频率. 在无预应力状态下, NTET质量越大、刚度越小时, 有更好的减振效果. 当弹簧预应力逐渐增大时, NTET的非线性刚度逐渐变弱, 会降低涡激振动控制性能. 参数分析表明: 随着涡激振动控制性能的提升, 圆柱体的振幅逐渐较小, NTET的振幅逐渐增大, 能量传递效率逐渐提高. 研究结果可为工程中涡激振动控制策略的高效设计提供有用的理论支撑和实验数据.      

深海细长柔性立管涡激振动数值分析方法研究进展

深海细长柔性立管的涡激振动与刚性圆柱的涡激振动问题相比,既有联系又有区别.当立管的长细比(L/D)较大时,会展现出一些刚性圆柱涡激振动没有的现象.本文介绍了近几十年来深海开发的细长柔性立管涡激振动的研究进展,主要侧重于数值计算分析研究领域,对目前常用的数值分析方法进行了归纳和总结.在此基础之上,详细介绍了目前理论研究和实际工程中的热点问题,诸如流场剖面的影响、立管顶端浮体运动的耦合、涡激振动抑制措施、振动多模态特性等.     

Numerical scheme for riser motion calculation during 3-D VIV simulation

This paper presents a numerical scheme for riser motion calculation and its application to riser VIV simulations. The discretisation of the governing differential equation is studied first. The top tensioned risers are simplified as tensioned beams. A centered space and forward time finite difference scheme is derived from the governing equations of motion. Then an implicit method is adopted for better numerical stability. The method meets von Neumann criteria and is shown to be unconditionally stable. The discretized linear algebraic equations are solved using a LU decomposition method. This approach is then applied to a series of benchmark cases with known solutions. The comparisons show good agreement. Finally the method is applied to practical riser VIV simulations. The studied cases cover a wide range of riser VIV problems, i.e. different riser outer diameter, length, tensioning conditions, and current profiles. Reasonable agreement is obtained between the numerical simulations and experimental data on riser motions and cross-flow VIV a/D. These validations and comparisons confirm that the present numerical scheme for riser motion calculation is valid and effective for long riser VIV simulation.     

Numerical study of the suppression mechanism of vortex-induced vibration by symmetric Lorentz forces

In this paper, the electro-magnetic control of vortex-induced vibration (VIV) of a circular cylinder is investigated numerically based on the stream function–vorticity equations in the exponential–polar coordinates attached on the moving cylinder for Re=150. The effects of the instantaneous wake geometries and the corresponding cylinder motion on the hydrodynamic forces for one entire period of vortex shedding are discussed using a drag–lift phase diagram. The drag–lift diagram is composed of the upper and lower closed curves due to the contributions of the vortex shedding but is magnified, translated and turned under the action of the cylinder motion. The Lorentz force for controlling the vibration cylinder is classified into the field Lorentz force and the wall Lorentz force. The symmetric field Lorentz force will symmetrize the flow passing over the cylinder and decreases the lift oscillation, which, in turn, suppresses the VIV, whereas the wall Lorentz force has no effect on the lift. The cylinder vibration increases as the work performed by the lift dominates the energy transfer. Otherwise, the cylinder vibration decreases. If the net transferred energy per motion is equal to zero, the cylinder will vibrate steadily or be fixed.     

Numerical study on the suppression of the vortex-induced vibration of an elastically mounted cylinder by a traveling wave wall

In the present paper, the commercial CFD code “Fluent” was employed to perform 2-D simulations of an entire process that included the flow around a fixed circular cylinder, the oscillating cylinder (vortex-induced vibration, VIV) and the oscillating cylinder subjected to shape control by a traveling wave wall (TWW) method. The study mainly focused on using the TWW control method to suppress the VIV of an elastically supported circular cylinder with two degrees of freedom at a low Reynolds number of 200. The cross flow (CF) and the inline flow (IL) displacements, the centroid motion trajectories and the lift and drag forces of the cylinder that changed with the frequency ratios were analyzed in detail. The results indicate that a series of small-scale vortices will be formed in the troughs of the traveling wave located on the rear part of the circular cylinder; these vortices can effectively control the flow separation from the cylinder surface, eliminate the oscillating wake and suppress the VIV of the cylinder. A TWW starting at the initial time or at some time halfway through the time interval can significantly suppress the CF and IL vibrations of the cylinder and can remarkably decrease the fluctuations of the lift coefficients and the average values of the drag coefficients; however, it will simultaneously dramatically increase the fluctuations of the drag coefficients.     

Hydrodynamics of a flexible cylinder under modulated vortex-induced vibrations

Both amplitude modulation and frequency modulation of Vortex-induced Vibration (VIV) are observed in a recent model test of a flexible cylinder under oscillatory flow, but its hydrodynamics has not yet been broached in detail. This paper employs the Forgetting Factor Least Squares (FF-LS) method for identification of time-varying hydrodynamics of a flexible cylinder under modulated VIV. The FF-LS method’s applicability to accurately identify time-varying hydrodynamic coefficients is demonstrated through an elastically mounted rigid cylinder under flow with a given modulated motion. Furthermore, we propose a framework to predict instantaneous amplitude (envelope) and frequency using time-varying hydrodynamic coefficients to establish their analytical relationship. This prediction method is further extended to a highly tensioned flexible cylinder through Fourier series expansion in the spatial domain. By performing the identification procedure for all sampled data of a flexible cylinder undergoing oscillatory flow, we obtain the corresponding time-varying hydrodynamics in the cross-flow direction considering the amplitude and frequency modulation. The results show that, under modulated VIV, hydrodynamic coefficients of the flexible cylinder also show time-varying characteristics. We further investigate differences between identified hydrodynamic coefficients and those obtained from the database of a cylinder with modulated motion under flow. Prediction results using these identified time-varying coefficients reveal that the time-varying excitation coefficients mainly influence the amplitude modulation, and the time-varying added-mass coefficients contain the major information of frequency modulation. These results further suggest including the temporal derivative of the instantaneous amplitude as one determining parameter in building databases to improve the prediction of modulated VIV.     

大柔性圆柱体两自由度涡激振动试验研究

基于模型试验研究了柔性圆柱体两自由度涡激振动问题, 研究了顺流向涡激振动和横向涡激振动的频率与振幅关系, 提出了考虑流固耦合的两自由度涡激振动非线性分析模型. 研究表明, 在不同的流速(雷诺数)范围, 柔性圆柱体顺流向涡激振动与横向涡激振动的频率比和幅值比是不同的; 在非锁定区, 圆柱体的顺流向振动频率与横向振动频率相同, 在锁定区, 圆柱体的顺流向振动频率是横向振动频率的两倍; 在非锁定区, 顺流向振幅与横向振幅比约为1, 而在锁定区, 顺流向振幅与横向振幅比约为1/3~2/3.      

Application of wake oscillators to two-dimensional vortex-induced vibrations of circular cylinders in oscillatory flows

A nonlinear time-domain simulation model for predicting two-dimensional vortex-induced vibration (VIV) of a flexibly mounted circular cylinder in planar and oscillatory flow is presented. This model is based on the utilization of van der Pol wake oscillators, being unconventional since wake oscillators have typically been applied to steady flow VIV predictions. The time-varying relative flow–cylinder velocities and accelerations are accounted for in deriving the coupled hydrodynamic lift, drag and inertia forces leading to the cylinder cross-flow and in-line oscillations. The system fluid–structure interaction equations explicitly contain the time-dependent and hybrid trigonometric terms. Depending on the Keulegan–Carpenter number (KC) incorporating the flow maximum velocity and excitation frequency, the model calibration is performed, entailing a set of empirical coefficients and expressions as a function of KC and mass ratio. Parametric investigations in cases of varying KC, reduced flow velocity, cylinder-to-flow frequency ratio and mass ratio are carried out, capturing some qualitative features of oscillatory flow VIV and exploring the effects of system parameters on response prediction characteristics. The model dependence of hydrodynamic coefficients on the Reynolds number is studied. Discrepancies and limitations versus advantages of the present model with different feasible solution scenarios are illuminated to inform the implementation of wake oscillators as a computationally efficient prediction model for VIV in oscillatory flows.     

A robust high-resolution method for the time–frequency analysis of vortex-induced-vibration signals

The dynamic interaction between ocean current and marine riser is complex in nature, and the riser׳s vortex-induced vibration (VIV) due to the current often strongly exhibits a non-stationary phenomenon. For investigating the time-varying dominant frequencies of the VIV motion, a joint time–frequency analysis is necessary. Traditional methods for time–frequency analysis include the Short Time Fourier Transform (STFT) and Wavelet Transform (WT) methods, though both methods have significant drawbacks. Specifically, the STFT method suffers frequency resolution and leakage problems, while the WT method is highly sensitive to its basic wavelet selection. This paper newly introduces a robust high-resolution method, named the STPT-SS method, which is the Short Time Prony Transform (STPT) using a State-Space (SS) model. In particular, the STPT algorithm contributes to the high-resolution feature of the proposed method, and the SS model to the robustness. Using test VIV data that include a synthesized signal and measurements from laboratory and field experiments, the STPT-SS method is found to significantly outperform the STFT and WT methods in the time–frequency analysis.     

Effects of tension on vortex-induced vibration (VIV) responses of a long tensioned cylinder in uniform flows

The effects of tension on vortex-induced vibra-tion (VIV) responses for a tension-dominated long cylinder with an aspect ratio of 550 in uniform flows are experimen-tally investigated in this paper. The results show that elevated tension suppresses fluctuations of maximum displacement with respect to flow velocity and makes chaotic VIV more likely to appear. With respect to periodic VIV, if elevated tension is applied, the dominant vibration frequency in the in-line (IL) direction will switch from a fundamental vibra-tion frequency to twice the value of the fundamental vibration frequency, which results in a ratio of the dominant vibration frequency in the IL direction to that in the cross-flow direc-tion of 2.0. The suppression of the elevated tension in the fluctuation of the maximum displacement causes the axial tension to become an active control parameter for the VIV maximum displacement of a tension-dominated long riser or tether of an engineering structure in deep oceans. However, the axial tension must be optimized before being used since the high dominant vibration frequency due to the elevated tension may unfavorably affect the fatigue life of the riser or tether.     

Three-dimensional numerical simulation of vortex-induced vibration of an elastically mounted rigid circular cylinder in steady current

Vortex-induced vibration (VIV) of an elastically mounted rigid circular cylinder in steady current is investigated by solving the three-dimensional Navier–Stokes equations. The cylinder is allowed to vibrate only in the cross-flow direction. The aim of this study is to investigate the variation of the vortex shedding flow in the axial direction of the cylinder and to study the transition of the flow from two-dimensional (2D) to three-dimensional (3D) for VIV of a cylinder. Simulations are carried out for a constant mass ratio of 2, the Reynolds numbers ranging from 150 to 1000 and the reduced velocities ranging from 2 to 12. The three-dimensionality of the flow is found to be the strongest in the upper branch of the VIV response and weakest in the initial branch. The 2S and 2P vortex shedding modes are found to coexist along the cylinder span in the upper branch, leading to strong variations of the lift coefficient in the axial direction of the cylinder. The difference between the flow transition from 2D to 3D in the VIV lock-in regime and that in the wake of a stationary cylinder is identified. The transition mode B found in the wake of a stationary cylinder is also found in the wake of a vibrating cylinder. The critical Reynolds number for flow transition from 2D to 3D of a cylinder undergoing cross-flow VIV at a reduced velocity of 6 is found to be greater than that for a stationary cylinder. For a constant reduced velocity of 6, the wake flow changes from 2D to 3D as the Reynolds number is increased from 250 to 300. Some 2D numerical simulations are performed and it is found that the 2D Navier–Stokes (NS) equations are not able to predict the VIV in the turbulent flow regime, while the 2D Reynolds-averaged Navier–Stokes (RANS) equations improve the results.