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

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

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

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

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

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

DYNAMIC RESPONSE OF THE CATENARY RISER IN DEEP-SEA MINING SYSTEMS 1)

基于有限元模型,模拟、分析深海采矿系统悬链线立管在海流和水面船运动约束下的动力响应。结果表明：悬链线立管的最大等效应力和最大位移随时间呈周期性变化,且存在半个周期的相位差;当水面船运动到最高点时,悬链线立管位移达到最大值。悬链线立管最大 等效应力和最大位移随水面船运动位移幅值的增加而增加,随运动周期的增加而减小。进一步对比发现,水面船运动位移和周期对立管等效应力的影响大于对其位移的影响。     

非对称槽道中涡旋波的特性研究

利用PIV流场显示技术，对振荡流体在非对称槽道中涡旋波的产生、发展和消失的规律进 行了实验研究和分析，测得了涡旋波流场的速度矢量图，阐明了涡旋波流场周期性变化的特 点. 结合涡动力学方程，深入分析并揭示了做周期性运动的流体能在槽道中产生波的特性这 一规律，从中发现：流体周期变化的非定常性和不对称的槽道结构是形成涡旋波流动的主要 因素. 本文对涡旋波流场中各个旋涡的速度分布和涡量进行了测量和计算，分析了涡旋波 强化传质的机理，并研究了Re数对涡旋波流动的影响     

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

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

内孤立波作用下潜深对潜体运动响应和载荷特性的影响研究

内孤立波常发生于海洋密度跃层, 因其峰高谷深、携带能量巨大, 在传播过程中会导致跃层上下的海水流动呈现剪切状态, 并引起突发性的强流. 潜体在水下悬停时极有可能会遭遇内孤立波, 由于内孤立波的流场特性, 置于跃层上下的悬浮潜体所产生运动响应和水动力载荷变化不尽相同, 甚者会出现掉深现象. 为探究潜深对波体耦合作用的影响, 基于不可压缩N-S方程和mKdV理论, 采用速度入口造波, 结合重叠网格技术和流固耦合方法, 建立了分层流中内孤立波耦合水下潜体多自由度运动的数值模型, 通过该模型分析了不同潜深下悬浮潜体的运动响应和载荷特性. 结果表明: 在内孤立波作用下, 位于密度跃层上方和跃层中的潜体顺着波的前进方向运动, 先下沉后抬升, 位于跃层下方的潜体则会逆流持续下沉; 潜体与波面的垂向距离越小, 对其纵荡、垂荡和速度的影响越显著, 而位于密度跃层中的潜体在分界面处沿着波形运动, 其运动响应和载荷变化受影响较小; 潜体在跃层上、下流体中所受水平力的方向相反, 水平力峰值小于垂向力峰值, 且位于跃层下方的潜体一直受到低头力矩, 最终导致掉深.      

柔性圆柱涡激振动流体力系数识别及其特性

涡激振动是诱发海洋立管、浮式平台系泊缆和海底悬跨管道等柔性圆柱结构疲劳损伤的重要因素.目前,海洋工程中用于柔性圆柱涡激振动预报的流体力系数主要来源刚性圆柱横流向受迫振动的实验数据,存在一定缺陷和误差.本文综合考虑横流向与顺流向振动耦合作用,建立了柔性圆柱涡激振动流体力模型,运用有限元法和最小二乘法确定升力系数、脉动阻力系数和附加质量系数.为了准确识别柔性圆柱涡激振动流体力系数,设计并开展了拖曳水池模型实验,实验用柔性圆柱模型的质量比为1.82,长径比为195.5.通过与刚性圆柱流体力系数对比,深入分析了柔性圆柱流体力系数的特性.结果表明:柔性圆柱在一阶模态控制区,流体力系数随约化速度变化趋势与刚性圆柱大致相似;二阶模态控制区,升力系数和脉动阻力系数显著增大;附加质量系数在响应频率较低时与振动位移的相关性增强;当响应频率较低时,振动位移较大区域为能量耗散区,当响应频率较高时,振动位移较大区域为能量输入区.     

分层流体中运动源生成的内波研究进展

针对两类密度分布模型------连续分层流体和间断分层流体, 综述了在运动潜体生成的Kelvin型和非Kelvin型内尾迹研究方面的现状, 内容侧重于运动源生成内波的解析理论和分层拖曳水槽中内尾迹实验方面的研究成果. 介绍了在连续分层流体中运动源生成的Kelvin型非线性内波的一般方程和在间断分层流体中Kelvin型内波的势流分析的一般方法; 概述了运动源诱生的先锋内孤立子、代数孤立子和平孤立波3类特殊非线性内波的研究进展, 其中运动潜体生成的平孤立内波被作者实验证实是一类极限孤立波, 并首次建立了共轭流动模型予以描述; 综合分析了在密度线性分布流体中潜体运动生成内波的动力学过程多样性特征, 其中包括内尾迹近场和远场的时空结构、不稳定结构、涡旋与湍流耦合结构以及湍流与内波相互作用结构等.      

?????????????????е???????

研究了埋置于弹性地基内充液压力管道中非线性波的传播. 假设管壁是线弹 性的,地基反力采用Winkler线性地基模型,管中流体为不可压缩理想流体. 假定系统初始 处于内压为$P_0$的静力平衡状态,动态的位移场及内压和流速的变化是叠加在静 力平衡状态上的扰动. 基于质量守恒和动量定理,建立了管壁和流体耦合作用的非 线性运动方程组; 进而用约化摄动法, 在长波近似情况下得到了KdV方程,表征 着系统有孤立波解.     

Dynamic response of a deepwater riser subjected to combined axial and transverse excitation by the nonlinear coupled model

In offshore engineering long slender risers are simultaneously subjected to both axial and transverse excitations. The axial load is the fluctuating top tension which is induced by the floater’s heave motion, while the transverse excitation comes from environmental loads such as waves. As the time-varying axial load may trigger classical parametric resonance, dynamic analysis of a deepwater riser with combined axial and transverse excitations becomes more complex. In this study, to fully capture the coupling effect between the planar axial and transverse vibrations, the nonlinear coupled equations of a riser’s dynamic motion are formulated and then solved by the central difference method in the time domain. For comparison, numerical simulations are carried out for both linear and nonlinear models. The results show that the transverse displacements predicted by both models are similar to each other when only the random transverse excitation is applied. However, when the combined axial dynamic tension and transverse wave forces are both considered, the linear model underestimates the response because it ignores the coupling effect. Thus the coupled model is more appropriate for deep water. It is also found that the axial excitation can significantly increase the riser’s transverse response and hence the bending stress, especially for cases when the time-varying tension is located at the classical parametric resonance region. Such time-varying effects should be taken into account in fatigue safety assessment.     

Experimental study on dynamic response of smooth and straked risers in tandem arrangement coupling interference effect

An experiment was conducted in a combined wave–current water flume on two tandem risers subjected to uniform flow. The riser model has an effective length of 2.0 m. The aspect ratio is 111.11. The upstream riser is smooth and the downstream riser fitted with three-strand helical strakes with pitch 17.5D and height 0.25D. By varying the external flow velocities and spacing ratios, through comparisons with the dynamic response of isolated smooth and isolated straked riser, the paper observe how interference effect impacts the dynamic characteristics and dynamic response of two risers in tandem arrangement, reveal how the suppression efficiency of the three-strand helical strakes responds to spacing ratio and external flow velocity, and explore the wake excitation effect of inter-riser fluids on the downstream riser and their dynamic feedback to the upstream riser. The results show that the dominant frequency of the upstream smooth riser is sensitivity to the change of the spacing ratio is low, and the displacement response is offset or enhanced in different degree due to the difference of the interference efficiency. The downstream straked riser dominates frequency and displacement higher than the isolated straked riser. The wake vortex of the upstream smooth riser acts on the downstream riser, occupying a dominant position in the vibration of the downstream riser. It degrades the vibration suppression efficiency of the three-strand helical strakes: the suppression efficiency is the highest at spacing ratio of 8, being a merely 70.57%. With the increase of the spacing ratio, the CF displacement of the riser gradually decreases, and the IL displacement gradually increases. The interference efficiency partition and suppression efficiency at different spacing ratios reflects the dynamic feedback of the upstream smooth riser is much smaller than the interference effect of the downstream suppression riser.     

Monitoring VIV fatigue damage on marine risers

Long flexible cylinders (e.g., risers, tendons and mooring lines) exposed to the marine environment encounter ocean currents leading to vortex-induced vibration (VIV). These oscillations, often driven at high frequencies over extended periods of time, may result in structural failure of the member due to fatigue damage accumulation. Recent developments in instrumentation and installation of data acquisition systems on board marine risers have made accurate measurement of riser responses possible. This paper aims at using the data from these data acquisition devices (typically strain gages and accelerometers) in order to understand the evolution of the riser VIV, with the final aim of estimating the fatigue damage. For this purpose we employ systematic techniques to reconstruct riser VIV response using the data from the available sensors. The reconstructed riser response allows estimation of the dynamic axial stresses due to bending and consequently the estimates of the fatigue damage along the entire riser. The above methods can take into account the fatigue damage arising from complicated riser motions involving the presence of traveling waves even with the use of very few sensors. An alternate approach using a Van der Pol wake oscillator model is also explored to obtain fatigue life estimates caused by riser VIV.     

A systematic approach to riser VIV response reconstruction

Vortex-induced vibration (VIV) of long flexible cylindrical structures (e.g. risers, pipelines, tendons, mooring lines) enduring ocean currents is ubiquitous in the offshore industry. Though significant effort has gone into understanding this complicated fluid–structure interaction problem, major challenges remain in modeling and predicting the response of such structures (for example a riser). The work presented in this paper provides a systematic approach to estimate and analyze the vortex-induced motions of a marine riser. A systematic framework is developed, which allows reconstruction of the riser motion from a limited number of sensors placed along its length. A full reconstruction criterion is developed, which classifies when the measurements from the sensors contain all information pertinent to riser VIV response, and when they do not, in which case additional, analytical methods must be employed. Reconstruction methods for both scenarios are developed and applied to experimental data. Finally, a systematic study on the error during the reconstruction is also undertaken. The methods developed in this paper can be applied to: improve understanding of the vortex shedding mechanisms, including the presence of traveling waves and higher-harmonic forces; develop tools for in-situ estimation of fatigue damage on marine risers; and estimate the vortex-induced forces on marine risers.     

Finite memory quadratic Volterra model for the response prediction of a slender marine structure under a Morison load

A quadratic Volterra model with a finite nonlinear memory effect was introduced and applied to the time series prediction of a slender marine structure exposed to the Morison load. First, the unknown nonlinear single-input–single-output dynamic system was identified using the nonlinear autoregressive with exogenous input (NARX) technique based on the prepared datasets of the wave elevation and system response, which was obtained by running nonlinear time domain analysis for a certain short term sea state. The structure of NARX was designed in such a way that the linear part had infinite memory, whereas the nonlinear part had finite memory of a certain length. Second, the frequency domain Volterra kernels, both linear and quadratic, were derived analytically by applying the harmonic probing method to the identified system. To derive the frequency response functions, the sigmoidal function used in NARX to realize the nonlinear relationship between the input and output was expanded to polynomials based on the Taylor series expansion, so that the harmonics of same frequencies were easily matched between the input and output. Finally, the time series of the system response under arbitrarily given short term sea states were predicted using the quadratic Volterra series. The proposed methodology was used to predict the nonlinear dynamic response of a 2-dimentional free standing catenary riser exposed to a random ocean wave load, and the comparison between the prediction and simulation results was made on the probability distribution of the maximum excursion of riser top. The results show that the proposed methodology can successfully capture the nonlinear effects of the dynamic response of a slender marine structure induced by the quadratic term of the Morison formula.     

被均匀流缓慢调制的有限水深毛细重力波

非线性的存在会产生高次谐波,这些谐波又反作用于原来的低次谐波,使波幅发生缓慢变化,从而产生缓慢调制现象．这里从考虑均匀流作用下的毛细重力水波基本方程出发,在不可压缩、无旋、无黏条件假设下,使用多重尺度分析方法推导出了在均匀流影响下有限深水毛细重力波振幅所满足的非线性Schr?dinger方程(NLSE)．分析了NLSE解的调制不稳定性．给出了毛细重力波调制不稳定的条件和钟型孤立波产生的条件．分析了无量纲最大不稳定增长率随无量纲水深和表面张力的变化趋势．同时给出了无量纲不稳定增长率随无量纲微扰动波数变化的曲线,呈现出了先增大后减小的趋势．最后指出均匀顺流减小了无量纲不稳定增长率及最大增长率,逆流增大了它们．由表面张力作用产生的毛细波及重力与表面张力共同作用产生的毛细重力波,与流的相互作用可以改变海表粗糙度和海洋上层流场结构,进而影响海气界面动量、热量及水汽的交换．了解海表这些短波动力机制,对卫星遥感的精确测量、海气相互作用的研究及海气耦合模式的改进等有重要意义．      

INTERNAL SOLITARY WAVE LOADS EXPERIMENTS AND ITS THEORETICAL MODEL FOR A CYLINDRICAL STRUCTURE

在大型重力式密度分层水槽中, 对内孤立波与圆柱型结构的相互作用特性开展了系列实验. 基于两层流体中 内孤立波的KdV,eKdV和MCC理论, 建立了圆柱型结构内孤立波载荷的理论预报模型, 给出了该载荷理论预报模型中3类内孤立波理论的适用性条件.研究表明, 圆柱型结构内孤立波水平载荷包括水平Froude-Krylov力、附加质量力和拖曳力3个部分, 可以由Morison公式计算, 而内孤立波垂向载荷主要为垂向Froude-Krylov力, 可以由内孤立波诱导动压力计算.系列实验结果表明, 附加质量系数可以取为常数1.0, 拖曳力系数与内孤立波诱导速度场的雷诺数之间为指数函数关系, 而且基于理论预报模型的数值结果与系列实验结果吻合.