内孤立波与Spar型浮式风机平台耦合作用研究_康信龙.pdf

第 卷第期 年月中 国 海 洋 大 学 学 报 （）：，内孤立波与 型浮式风机平台耦合作用研究*康信龙，王俊荣，*，谢波涛，何镇阳（中国海洋大学工程学院，山东 青岛 ；中国海洋大学山东省海洋工程重点实验室，山东 青岛 ；中海石油（中国）有限公司北京研究中心，北京 ）摘要：针对内孤立波与系泊浮式平台的耦合作用问题，本文以单立柱 型浮式风机平台为研究对象，研究了平台在内孤立波作用下的耦合运动特性。利用重叠网格技术及自编耦合动力程序，以计算流体动力学求解 流体控制方程，建立了 型浮式风机平台与内孤立波（，）流固耦合作用数值模型。通过系列数值模拟，研究了 型浮式风机平台在内孤立波作用下的水平动态载荷和纵荡、横荡规律。研究表明，内孤立波对 平台会产生较显著的水平作用力，并产生大幅的水平位移；在水深一定时，同一流体分层比下，随着内孤立波振幅提高，纵荡运动幅值增大；当振幅一定时，随着流体分层比的增大，纵荡幅值减小；通过调整系泊系统，发现非线性系泊系统下平台纵荡幅值要大于线性系泊系统下的情况，即线性系泊系统会低估纵荡运动幅值。关键词：平台；内孤立波；重叠网格；流固耦合作用；浮式风机平台中图法分类号：文献标志码：文章编号：（）：引用格式：康信龙，王俊荣，谢波涛，等内孤立波与 型浮式风机平台耦合作用研究中国海洋大学学报（自然科学版），（）：，（）：*基金项目：国家自然科学基金项目（）资助 （）收稿日期：；修订日期：作者简介：康信龙（），男，硕士生，研究方向：数值仿真。：*通讯作者：海上风能资源丰富，近年来海上风电开发已成为新能源开发的热点，对于深远海风电的开发，往往优先考虑浮式平台基础，如现已大规模应用在 海域的 单点柱型（）浮式风机平台。南海作为中国未来海上风电开发的主战场，同时也是内孤立波频发的海域。内孤立波具有振幅大、流速快等特点，这种强非线性波浪会对平台的安全运营产生威胁。目前针对浮式结构物与内孤立波作用的荷载及运动响应特性研究，仍是海洋工程界面临的挑战性难题。以往的研究主要是将某一类内孤立波传播理论模型与 公式结合来研究浮式结构物的荷载动力响应，如 等将 公 式 与 （）理论结合研究了 平台的运动响应问题。尤云翔等基于 理论采用改进的 公式研究了两层流体中张力腿平台和半潜式平台的动力响应特性。马孟达和黄文昊等在考虑 、（）、（）三类内孤立波理论的适用性条件后，分别建立了张力腿（）平台、平台的荷载理论模型，并在此基础上研究了结构物的运动响应特性。等通过开展系列模型实验，研究了一浮体在内波作用下的垂荡、纵荡和纵摇特点，同时指出系泊张力的幅值同内波振幅正相关。在现有的文献中，计算结构物的内孤立波荷载主要是采用 公式，但其中采用的拖曳系数和惯性力系数值大部分是参考表面波来选取，其中会产生一定偏差。而计算流体力学（）方法凭借在处理非线性波浪问题上具有的显著优势，开始逐渐应用在浮式结构与内波相互作用的研究中。文献 中采用 方法研究了 平台、平台、半潜式平台的内波荷载特性，研究表明三类平台所受的水平力和垂向力主要由波浪力、黏性压差力和摩擦力组成，其中摩擦力与其他两种力相比可以忽略不计。等通过建立三维数值内波水槽，研究了不同浪向角下一小比尺半潜平台的整体受力情况及结构周围的复杂流场特性。需要指出的是，在这些文献中，往往只研究固定浮式结构物的荷载，并未考虑结构物的运动状态，无法准确得出结构物所受动态载荷。同时，对于实际工程中往往需要准确预报平台所受的内波动态荷载才能更好地评估平台运动响应特性，此时需要考虑到结构物对内波流场的影响，即考虑两者的耦合作用，并基于此开展内波作用下浮式平台的运动特性研究。本文以 数值模拟方法，采用重叠网格技术以期康信龙，等：内孤立波与 型浮式风机平台耦合作用研究六自由度运动方式（）研究了一 型浮式风机平台在内孤立波作用下的运动特性。依据文献 给出的 、内孤立波理论的适用性条件，建立振幅及波形可控的内孤立波数值水槽。在此基础上，开展 平台与不同振幅、分层比下内孤立波的耦合作用机理研究，分析平台所受的水平动态载荷及纵、横荡特点；同时通过调整系泊的非线性刚度，模拟分析非线性系泊对平台纵荡的影响。数值方法本文采用两层流体模式来研究内孤立波，建立如图所示的两层流体系统。设稳定分层的两层流体的总深度为，上、下层深度分别为和，上、下层密度分别为和。建立直角坐标系 ，其中平面 位于未扰动的静止的分层界面上，轴与平台中心轴重合且垂直向上为正。内孤立波沿 轴正向传播，表示内孤立波界面位移。采用速度入口边界法生成内孤立波，并经过传播区形成稳定内孤立波后，研究 平台所受的内孤立波动态载荷及耦合运动特性。图内孤立波数值水槽示意图 控制方程与湍流模型计算中假设流体为不可压黏性流，则结构物周围流场控制方程包括连续性方程：。（）动量方程：。（）式中：为哈密顿算子；为流体速度矢量；为物质导数；为时间；为流体元单位质量的体积力；为流体密度；为压力；为流体运动黏度。设内孤立波界面位移为，相速度为，则上、下层流体层平均速度需要满足：（）。（）式中：为分层流体水平流速；为分层流体深度；为流体层数序号。速度入口处上下层速度分别取和，采用体积分数（，）模型追踪内孤立波的界面变化，消波方式采用阻尼消波和开渠边界相结合 的方式。流场湍流模型采用的是标准-模型，壁面函数采用标准壁面函数，模型中湍动能和湍流耗散率的方程为：（）（）（）。（）（）（）（）（）（）。（）式中：为平均速度梯度产生的湍流动能；为模型修正系数；、和为耗散方程的经验参数，一般由试验确定或经验公式计算得到，式中各 封 闭常 数为 .，.，（）（）；.；（）；.，.；为有效黏度系数；和分别为湍动能及其耗散率的湍流普朗特数；为平均应变率。流固耦合原理计算物体耦合运动时需要使用到动网格技术，对于大幅度偏移运动问题，使用重叠网格技术具有一定的优势，其在计算过程中能够保持网格质量不变，使流固耦合运动计算更稳定。流固耦合的实现过程是：每一时间步内，将 求解得到的外部激励及系泊约束代入到平台运动微分方程，求解得到平台下一步位移、速度及加速度，再将数据传递回 ，更新平台位置及流场信息，重复迭代计算（见图）。图耦合运动求解流程图 中国海洋大学学报 年计算域与网格布置几何模型本文选取典型 风机平台 作为研究对象，同时考虑到计算成本及计算效率采用 的缩尺比，表给出了该 平台的主尺度。表 平台模型计算参数 计算参数 原型尺度 模型尺度 吃水 下部直径 上部直径 质量 系泊系统平台模型采用等效水平系泊提供回复力，系泊系统由根线性弹簧组成，坐标轴建立在平台的重心位置，每根弹簧与坐标轴呈 角，其布置方式如图所示。原型尺度下，系泊刚度 ，预张力 ，系泊系统阻尼比为 。图系泊布置方式示意图 计算域网格和边界条件图给出了计算域整体网格划分。计算域分为背景网格和前景网格两部分，背景网格区域尺寸为 ；前景网格区域尺寸为 。考虑多自由度运动时，尽量减少网格数量，同时保证计算精度，背景网格采用结构化网格，网格单元数量为 万；前景网格采用非结构化网格，对平台周围网格加密处理，网格单元数量为 万，两者组装后网格单元总数量为 万。（）计算域和边界条件示意图 ；（）浮式风机平台基础 ；（）计算域整体网格划分 ；（）平台周围网格 ）图计算域示意图及网格划分 计算域左侧边界设为速度入口，通过 调整入口速度；右侧为压力出口边界，给定出口静压力为；上、下表面设为对称边界；在平台立柱表面，采用无滑移壁面；前景区域和背景区域之间采用重叠网格边界，期康信龙，等：内孤立波与 型浮式风机平台耦合作用研究个计算域之间的数据传递主要采用界面耦合技术，通过插值方法实现重叠界面的数据交换。数值模拟工况设置本文 参 照 南 海 北 部 流 花 海 域 内 孤 立 波 实 测 数据，得出该海域水深在 ，内孤立波振幅在 ，主要以 为主。本文以此为依据，设计作业水深为 ，典型下凹型内孤立波振幅为 。上层流体密度为 ，下层流体密度为 ，模拟不同内孤立波振幅下平台的耦合运动，分析不同工况下平台的运动特点，具体计算工况见表。根据具体工况，需要采用适合的内孤立波理论模型，即根据文献 中确定的类内孤立波理论 、的适用性条件，计算得到式（）中的界面位移。表数值模拟工况 工况 分层比 振幅水深比 、数值模拟结果与分析内孤立波数值模拟结果为了生成稳定传播的内孤立波，首先需要校验本文所建内孤立波水槽的合理性，本文主要从内孤立波的生成传播特性、波形、最大波致流水平流速三方面进行说明。图给出了在，时，即 工况下内孤立波的数值模拟结果，可以看出在内孤立波的传播过程中，内孤立波波形保持稳定不变，波高衰减很小。因此，本文所产生的内孤立波是稳定的。图内孤立波的传播特性 图分别给出了分层比时种工况下内孤立 波 数 值 模 拟 结 果，工 况 的 振 幅 水 深 比 ，此时为弱非线性和弱色散性，采用 理论；工况 的振幅水深比 ，此时为中等非线性和弱色散性，采用 理论；工况 的振幅水深比 ，此时为强非线性和弱色散性，采用 理论；从波形对比结果，可以看出 数值模拟结果分别与理论解拟合较好。因此，从波形结果也可以看出本文所采用的的数值造波方法是准确可靠的。图同时给出了在工况、和 下的最大波致流水平流速垂向分布，即波谷处的水平流速垂向分布，可以看出 数值模拟的上、下层流体水平瞬时速度与理论计算拟合较好，主要区别体现在两层流体界面与波谷间的过渡区域，与理论解在过渡区域速度突变不同，很好地模拟出过渡区域速度渐变的过程，也为内波与结构物的耦合作用模拟提供合理的内波环境。图内孤立波波形对比 中国海洋大学学报 年图最大波致流水平流速对比 荷载与动力响应时程变化特性为方便陈述，本文分别定义?、?、?为平台无因次内孤立波水平力、平台纵荡、平台横荡。其中为平台迎流面积，为平台吃水深度。图给出了，时平台所受无因次水平力的时程变化曲线。结果表明，平台所受水平力在内孤立波波谷经过平台轴线的附近某个时刻达到峰值，随后逐渐减小，并且在内孤立波经过平 台 后，伴 随 着 尾 波 列 的 作 用 会 产 生 小 幅 的振荡。图给出了，时平台纵荡、横荡的时程曲线图。由图可知，平台会产生较大幅度的纵荡，其运动趋势与平台所受水平力基本一致，纵荡位移随着内孤立波的传播开始增大，在达到峰值后开始减小，当内孤立波经过平台后由于受尾波列的作用还会产生小幅振荡。同时在横向方向，平台还会产生小幅度的振荡，其运动峰值总要滞后于纵荡运动峰值，主要原因是平台横向运动的能量来源于纵荡方向的能量传递效应，且能量主要还是表现在纵荡方向。图，时平台所受水平力时程曲线图 ，.图，时平台纵荡（）和横荡（）时程曲线图 （）（），.期康信龙，等：内孤立波与 型浮式风机平台耦合作用研究荷载与动力响应幅值特性本节进一步探究了内孤立波作用下平台所受荷载及运动响应幅值的规律。图 给出了水深比分别为、和时，平台所受的水平力峰值随波幅变化的情况。可以看出同一分层比下，随着内孤立波振幅水深比的增大，水平力的峰值增大。同时从图中还可以发现，分层比对平台动态荷载幅值会产生影响，随着分层比的增大，即上层流体厚度变厚，同一振幅下水平载荷开始减小。图 则给出了种分层比下平台纵荡和横荡幅值变化特性，平台的运动响应与其受力有关，同一水深比条件下，平台纵荡幅值随波幅的增大而增大，且增长幅度因分层比的不同而不同，同时随着分层比的增大，即流体分界面下移，下层流体厚度变小，此时根据其所受水平力特性，纵荡幅值 减小；而横 荡运 动 较 为 复 杂，在 分层比为或时，横荡峰值与波幅呈正相关，而在分层比为的条件下会出现拐点及局部最小值。图 不同分层比下平台水平力幅值变化特性 图 不同分层比下平台纵荡（）与横荡（）幅值变化特性 （）（）系泊系统调整系泊系统对平台运动往往起着重要的作用，特别是对纵荡，为研究系泊刚度对平台纵荡运动影响，设定三维数值水槽模型参数为：上层水深 ；下层水深 。表调整系泊计算工况 工况 振幅水深比 为了研究系泊系统刚度改变对平台纵荡的影响，对系泊系统刚度做出调整，选取了种类型的系泊刚度，分别为线性刚度、强非线性刚度、弱非线性刚度，其位移恢复力曲线见图 并采用表中的工况对其进行计算分析。图 平台系泊位移恢复力曲线 中国海洋大学学报 年图 给出了在分层比时采用种不同系泊刚度的平台纵荡历时曲线图，从图中可以看出在同一波幅下，平台的纵荡趋势一致，但非线性刚度对平台的纵荡运动峰值影响显著，平均峰值比线性刚度大 左右，这主要是因为在平台运动的前期，线性刚度系泊给予结构的拉力要大于非线性系统所给予的拉力，因此系泊系统允许平台的运动位移量要小；同时还发现在非线性刚度下，弱非线性系统运动幅值要小于强非线性系统运动幅值，这是由于受流场及各自系泊系统的影响，两者对平台的运动影响程度不同。因此，由以上分析可知，系泊系统采用等效线性刚度弹簧时，往往会低估平台的纵荡幅值。（），；（），；（），）图 调整系泊系统平台纵荡时历曲线 结论（）内孤立波作用下，平台会受到很大的突发性水平冲击荷载并表现出较大的水平偏移。当分层比，振幅水深比 时，得到的最大波致流水平速度为 ，与文献 中该海域观测到的最大波致流速相近，此时平台水平 偏移 量 达 到 的水深。（）随内孤立波振幅的增大，平台所受的水平力、纵荡、横荡峰值都呈增大的趋势。其中平台纵荡运动幅度较为明显，伴随小幅的横荡运动，且横荡峰值较纵荡滞后出现；随着上层流体厚度变厚，平台所受的水平力、纵荡幅值均减小。（）调整系泊系统分析可知，平台纵荡运动趋势是一致，主要区别体现在纵荡幅值；相比于非线性系统，采用线性系泊系统会低估纵荡峰值平台约；同时在非线性系统下，弱非线性系统对平台的约束能力要优于强非线性系统。研究表明，在中国南海北部海域进行海上风电开发时，需要充分考虑内孤立波对浮式风电平台的影响，本文为内孤立波与浮式结构物耦合运动分析提供了一种新的手段。参考文献：，（）：尤云祥，李巍，时忠民，等 海洋内孤立波中张力腿平台的水动力特性 上海交通大学学报，（）：，（）：马孟达，尤云祥，张新曙海洋内孤立波作用下张力腿平台动力响应特性水动力学研究与进展辑，（）：期康信龙，等：内孤立波与 型浮式风机平台耦合作用研究 ，（）：黄文昊，尤云祥，王旭，等圆柱型结构内孤立波载荷实验及其理论模型 力学学报，（）：，（）：，：，：王旭，林忠义，尤云祥 内孤立波与直立圆柱体相互作用特性数值模拟哈尔滨工程大学学报，（）：，（）：，（）：，（）：，：，（）：，（）：，（）：石新刚，刘耀华，兰志刚，等南海北部流花海域内孤立波特征研究 热带海洋学报，（）：，（）：，（）：，（，；，；，）：（）（），；，；，；，：；责任编辑高蓓