CN116702648B - 一种三维钻柱的涡激振动计算方法及装置 - Google Patents
一种三维钻柱的涡激振动计算方法及装置 Download PDFInfo
- Publication number
- CN116702648B CN116702648B CN202310699786.XA CN202310699786A CN116702648B CN 116702648 B CN116702648 B CN 116702648B CN 202310699786 A CN202310699786 A CN 202310699786A CN 116702648 B CN116702648 B CN 116702648B
- Authority
- CN
- China
- Prior art keywords
- drill string
- dimensional
- preset
- parameters
- fluid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 29
- 239000012530 fluid Substances 0.000 claims abstract description 137
- 238000000034 method Methods 0.000 claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims description 103
- 238000005553 drilling Methods 0.000 claims description 35
- 238000013016 damping Methods 0.000 claims description 26
- 238000012545 processing Methods 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 12
- 230000003068 static effect Effects 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 4
- 230000004044 response Effects 0.000 abstract description 7
- 238000004088 simulation Methods 0.000 abstract description 7
- 238000004422 calculation algorithm Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 13
- 238000006073 displacement reaction Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 230000001133 acceleration Effects 0.000 description 7
- 238000004891 communication Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000000638 solvent extraction Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000005452 bending Methods 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Earth Drilling (AREA)
Abstract
本申请公开了一种三维钻柱的涡激振动计算方法及装置,方法包括基于预设的单元节点个数将三维钻柱切分为至少两个二维流场平面,并确定出对应的钻柱截面参数;基于每个钻柱截面参数建立二维数值模型,并对每个二维数值模型进行划分处理;从每个二维数值模型中确定出流体区域,并计算出每个流体区域的流场特性;基于流场特性确定出三维钻柱的受力参数以及三维钻柱的速度参数;根据速度参数以及预设的旋转角速度,得到三维钻柱的运动参数。通过将三维结构划分为多个二维平面,并结合预设的求解算法求解钻柱结构动力学方程,实现三维钻柱涡激振动数值模拟,以实现三维钻柱结构动力学响应的快速求解,不仅有效降低数值计算量,还可保障计算精度。
Description
技术领域
本申请属于海洋工程技术领域,特别的涉及一种三维钻柱的涡激振动计算方法及装置。
背景技术
50多年来,国际大洋科学钻探计划在全球各大洋已钻井4000余口,取回岩芯40多万米,开辟了探索地球深部的途径,推动了地球科学的革命性进步。近年来,随着我国与国际大洋发现计划(IODP)的深入合作,我国在大洋科学钻探领域已取得了长足进步。以往科学钻探主要通过立管钻探技术实现,随着IODP向深海发展,立管钻探工作时所需的附加设备以及安装回收引起的工时等都会增加成本,不仅如此,由于立管与海水直接接触,其承受载荷复杂、失效的风险较大,一旦发生失效,立管内的泥浆泄漏也将会对海洋环境造成巨大破坏。同时,作业水深增加必然导致立管使用数量增加,从而增大外载荷,对钻井平台具备的承载能力也是一个挑战。基于此,安全高效的无立管钻探技术成为当前的迫切需求。当进行无立管钻探工作时,海流流经钻柱,其表面发生周期性的漩涡脱落,引起涡激振动现象(Vortex Induced Vibration,简称VIV)。钻柱自身的旋转将进一步影响钻柱表面漩涡的生成和脱落,导致复杂的漩涡脱落和VIV响应。
在早期的研究中,学者们通过实验进行钻柱涡激振动研究,但实验结果很难实现流场可视化,且成本较高。二十一世纪以来计算机技术和计算流体力学快速发展,使数值模拟计算成为了研究钻柱涡激振动的另一种重要途径。然而,对于海洋钻柱这一类细长柔性旋转体,在涡激振动过程中会发生弯曲变形,采用一般的二维数值模拟无法正确体现其结构动力学响应机理;同时基于当前的计算机发展水平,采用三维数值模拟非常困难,不仅计算量非常大,而且精度无法保障。
发明内容
本申请为解决上述提到的对于海洋钻柱这一类细长柔性旋转体,在涡激振动过程中会发生弯曲变形,采用一般的二维数值模拟无法正确体现其结构动力学响应机理;同时基于当前的计算机发展水平,采用三维数值模拟非常困难,不仅计算量非常大,而且精度无法保障等技术缺陷,提出一种三维钻柱的涡激振动计算方法及装置,其技术方案如下:
第一方面,本申请实施例提供了一种三维钻柱的涡激振动计算方法,包括:
基于预设的单元节点个数将三维钻柱切分为至少两个二维流场平面,并确定出与每个二维流场平面对应的钻柱截面参数;其中,每个单元节点对应于一个二维流场平面;
基于每个钻柱截面参数建立相应的二维数值模型,并根据预设的距离区间,对每个二维数值模型进行划分处理;
从经过划分处理后的每个二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个流体区域的流场特性;
基于所有流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数;
根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。
在第一方面的一种可选方案中,根据预设的距离区间,对每个二维数值模型进行划分处理,包括:
在每个二维数值模型中确定出三维钻柱的截面中心;
将与三维钻柱的截面中心之间的距离处于预设的第一距离区间的区域划分为刚性区域;
将与三维钻柱的截面中心之间的距离处于预设的第二距离区间的区域划分为动网格区域;
将与三维钻柱的截面中心之间的距离处于预设的第三距离区间的区域划分为静网格区域;
对经过划分处理后的每个二维数值模型进行标记处理;
其中,预设的第一距离区间中的最大值等于预设的第二距离区间中的最小值,预设的第二距离区间中的最大值等于预设的第三距离区间的最小值。
在第一方面的又一种可选方案中,从经过划分处理后的每个二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个流体区域的流场特性,包括:
将经过划分处理后的每个二维数值模型中除去刚性区域的所有区域作为每个二维数值模型的初始流体区域;
基于预设的边界条件,对每个二维数值模型的初始流体区域进行裁剪处理,得到每个二维数值模型的目标流体区域;
根据每个二维数值模型的目标流体区域以及预设的流体参数,计算出每个流体区域的流场特性。
在第一方面的又一种可选方案中,基于所有流体区域的流场特性确定出三维钻柱的受力参数,包括:
从每个流体区域的流场特性中分别提取出三维钻柱的表面压力载荷;
对每个三维钻柱的表面压力载荷中处于顺流向的压力载荷进行汇总处理,得到三维钻柱的顺流向载荷矩阵;
对每个三维钻柱的表面压力载荷中处于横流向的压力载荷进行汇总处理,得到三维钻柱的横流向载荷矩阵。
在第一方面的又一种可选方案中,根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数,包括:
根据所有钻柱截面参数分别计算出三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵;
将三维钻柱的顺流向载荷矩阵、三维钻柱的横流向载荷矩阵、三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵代入至预设的钻柱动力方程,并基于预设的求解参数对钻柱动力方程进行计算,得到三维钻柱的速度参数。
在第一方面的又一种可选方案中,当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数,包括:
基于三维钻柱的运动参数绘制振幅曲线,并判断振幅曲线是否处于预设的振幅区间;
当检测到振幅曲线处于预设的振幅区间时,确定三维钻柱的运动参数满足预设条件,并输出三维钻柱的运动参数;或
基于三维钻柱的运动参数绘制振幅曲线,并判断振幅曲线中是否存在至少两个振幅值相同的点;
当检测到振幅曲线中存在至少两个振幅值相同的点时,确定三维钻柱的运动参数满足预设条件,并输出三维钻柱的运动参数。
在第一方面的又一种可选方案中,在根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数之后,还包括:
当检测到三维钻柱的运动参数不满足预设条件时,基于三维钻柱的运动参数对每个流体区域进行更新处理;
基于预设的流体参数,计算出更新处理后的每个流体区域的流场特性;
基于所有更新处理后的流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数;
根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。
第二方面,本申请实施例提供了一种三维钻柱的涡激振动计算装置,包括:
第一处理模块,用于基于预设的单元节点个数将三维钻柱切分为至少两个二维流场平面,并确定出与每个二维流场平面对应的钻柱截面参数;其中,每个单元节点对应于一个二维流场平面;
第二处理模块,用于基于每个钻柱截面参数建立相应的二维数值模型,并根据预设的距离区间,对每个二维数值模型进行划分处理;
第一计算模块,用于从经过划分处理后的每个二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个流体区域的流场特性;
第二计算模块,用于基于所有流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数;
参数输出模块,用于根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。
在第二方面的一种可选方案中,第二处理模块用于:
在每个二维数值模型中确定出三维钻柱的截面中心;
将与三维钻柱的截面中心之间的距离处于预设的第一距离区间的区域划分为刚性区域;
将与三维钻柱的截面中心之间的距离处于预设的第二距离区间的区域划分为动网格区域;
将与三维钻柱的截面中心之间的距离处于预设的第三距离区间的区域划分为静网格区域;
对经过划分处理后的每个二维数值模型进行标记处理;
其中,预设的第一距离区间中的最大值等于预设的第二距离区间中的最小值,预设的第二距离区间中的最大值等于预设的第三距离区间的最小值。
在第二方面的又一种可选方案中,第一计算模块用于:
将经过划分处理后的每个二维数值模型中除去刚性区域的所有区域作为每个二维数值模型的初始流体区域;
基于预设的边界条件,对每个二维数值模型的初始流体区域进行裁剪处理,得到每个二维数值模型的目标流体区域;
根据每个二维数值模型的目标流体区域以及预设的流体参数,计算出每个流体区域的流场特性。
在第二方面的又一种可选方案中,第二计算模块用于:
从每个流体区域的流场特性中分别提取出三维钻柱的表面压力载荷;
对每个三维钻柱的表面压力载荷中处于顺流向的压力载荷进行汇总处理,得到三维钻柱的顺流向载荷矩阵;
对每个三维钻柱的表面压力载荷中处于横流向的压力载荷进行汇总处理,得到三维钻柱的横流向载荷矩阵。
在第二方面的又一种可选方案中,第二计算模块还用于:
根据所有钻柱截面参数分别计算出三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵;
将三维钻柱的顺流向载荷矩阵、三维钻柱的横流向载荷矩阵、三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵代入至预设的钻柱动力方程,并基于预设的求解参数对钻柱动力方程进行计算,得到三维钻柱的速度参数。
在第二方面的又一种可选方案中,参数输出模块用于:
基于三维钻柱的运动参数绘制振幅曲线,并判断振幅曲线是否处于预设的振幅区间;
当检测到振幅曲线处于预设的振幅区间时,确定三维钻柱的运动参数满足预设条件,并输出三维钻柱的运动参数;或
基于三维钻柱的运动参数绘制振幅曲线,并判断振幅曲线中是否存在至少两个振幅值相同的点;
当检测到振幅曲线中存在至少两个振幅值相同的点时,确定三维钻柱的运动参数满足预设条件,并输出三维钻柱的运动参数。
在第二方面的又一种可选方案中,装置还包括:
在根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数之后,当检测到三维钻柱的运动参数不满足预设条件时,基于三维钻柱的运动参数对每个流体区域进行更新处理;
基于预设的流体参数,计算出更新处理后的每个流体区域的流场特性;
基于所有更新处理后的流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数;
根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。
第三方面,本申请实施例还提供了一种三维钻柱的涡激振动计算装置,包括处理器以及存储器;
处理器与存储器连接;
存储器,用于存储可执行程序代码;
处理器通过读取存储器中存储的可执行程序代码来运行与可执行程序代码对应的程序,以用于实现本申请实施例第一方面或第一方面的任意一种实现方式提供的三维钻柱的涡激振动计算方法。
第四方面,本申请实施例提供了一种计算机存储介质,计算机存储介质存储有计算机程序,计算机程序包括程序指令,程序指令当被处理器执行时,可实现本申请实施例第一方面或第一方面的任意一种实现方式提供的三维钻柱的涡激振动计算方法。
在本申请实施例中,可在计算三维钻柱的涡激振动数值时,基于预设的单元节点个数将三维钻柱切分为至少两个二维流场平面,并确定出与每个二维流场平面对应的钻柱截面参数;其中,每个单元节点对应于一个二维流场平面;基于每个钻柱截面参数建立相应的二维数值模型,并根据预设的距离区间,对每个二维数值模型进行划分处理;从经过划分处理后的每个二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个流体区域的流场特性;基于所有流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数;根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。通过将三维结构划分为多个二维平面,并结合预设的求解算法求解钻柱结构动力学方程,实现三维钻柱涡激振动数值模拟,以实现三维钻柱结构动力学响应的快速求解,不仅有效降低数值计算量,还可保障计算精度。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例中所需使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的一种三维钻柱的涡激振动计算方法的整体流程图;
图2为本申请实施例提供的一种二维数值模型的划分效果示意图;
图3为本申请实施例提供的又一种二维数值模型的划分效果示意图;
图4为本申请实施例提供的一种三维钻柱的涡激振动流场涡量等值线示意图;
图5为本申请实施例提供的一种三维钻柱的涡激振动振幅示意图;
图6为本申请实施例提供的一种三维钻柱的涡激振动计算装置的结构示意图;
图7为本申请实施例提供的又一种三维钻柱的涡激振动计算装置的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述。
在下述介绍中,术语“第一”、“第二”仅为用于描述的目的,而不能理解为指示或暗示相对重要性。下述介绍提供了本申请的多个实施例,不同实施例之间可以替换或者合并组合,因此本申请也可认为包含所记载的相同和/或不同实施例的所有可能组合。因而,如果一个实施例包含特征A、B、C,另一个实施例包含特征B、D,那么本申请也应视为包括含有A、B、C、D的一个或多个所有其他可能的组合的实施例,尽管该实施例可能并未在以下内容中有明确的文字记载。
下面的描述提供了示例,并且不对权利要求书中阐述的范围、适用性或示例进行限制。可以在不脱离本申请内容的范围的情况下,对描述的元素的功能和布置做出改变。各个示例可以适当省略、替代或添加各种过程或组件。例如所描述的方法可以以所描述的顺序不同的顺序来执行,并且可以添加、省略或组合各种步骤。此外,可以将关于一些示例描述的特征组合到其他示例中。
请参阅图1,图1示出了本申请实施例提供的一种三维钻柱的涡激振动计算方法的整体流程图。
如图1所示,该三维钻柱的涡激振动计算方法至少可以包括以下步骤:
步骤102、基于预设的单元节点个数将三维钻柱切分为至少两个二维流场平面,并确定出与每个二维流场平面对应的钻柱截面参数。
在本申请实施例中,三维钻柱的涡激振动计算方法可应用于安装有ICEM CFD软件以及Fluent软件的控制终端,该控制终端可先根据用户输入的单元节点个数或是程序预设的单元节点个数,将需要计算涡激振动数值的三维钻柱划分为至少两个二维流场平面,并可将每个二维流场平面所对应的钻柱截面参数输入至ICEM CFD软件,以由该ICEM CFD软件根据每个钻柱截面参数生成各自对应的二维数值模型,以及对每个二维数值模型进行划分处理。接着,控制终端可将每个经过划分处理后的二维数值模型所对应的文件导入至Fluent软件,以由该Fluent软件结合所有经过划分处理后的二维数值模型、预设的边界参数以及预设的求解参数,计算出三维钻柱的涡激振动数值。可以理解的是,控制终端可在该Fluent软件运行之前,获取用户输入的基于C语言所编译的UDF二次开发程序,并通过Fluent软件中的User Defined Functions窗口将该UDF程序嵌入该Fluent软件,以通过嵌入有UDF程序的Fluent软件,保障输出结果准确反应出旋转圆柱涡激振动响应情况。
具体地,在对三维钻柱计算涡激振动数值时,控制终端可以但不局限于结合预设的单元节点,将三维钻柱结构沿着钻柱轴向离散切分为至少两个二维流场平面,其中,每个单元节点可对应一个二维流场平面,每个二维流场平面中包括有该三维钻柱沿着钻柱轴向所划分出的一个截面(一般为圆形),且每个截面所对应的钻柱截面参数各不相同。此处,钻柱截面参数可以但不局限于包括截面面积、截面的直径、截面相对Y轴的惯性矩、截面相对Z轴的惯性矩以及截面相对X轴的极惯性矩,该截面面积以及截面的直径可根据各个单元节点所在三维钻柱上的位置进行确定,该截面相对Y轴的惯性矩、截面相对Z轴的惯性矩以及截面相对X轴的极惯性矩则可根据确定出的截面面积以及空间直角坐标系进行确定。需要说明的是,每个截面可以但不局限于以截面中心为原点建立空间直角坐标系,该空间直角坐标系可选择以截面所在平面作为XOY平面,且不限定于此。
需要注意的是,在本申请实施例中是将三维钻柱结构视为有一定壁厚的立管,为了便于三维流场和钻柱的涡激振动数值计算,可基于切片理论以及Euler-Bernoulli梁模型,将该三维钻柱结构离散为n个梁单元(也可理解为将一根简支梁离散成n段),其中,梁单元与单元节点的关系为单元节点个数等于梁单元个数加1,也即该三维钻柱设置有n+1个单元节点,且每个单元节点处均可沿着钻柱轴向将三维钻柱的三维流场空间离散为n+1个二维流场平面。
步骤104、基于每个钻柱截面参数建立相应的二维数值模型,并根据预设的距离区间,对每个二维数值模型进行划分处理。
具体地,在将三维钻柱划分为多个二维流场平面,并确定出与每个二维流场平面对应的钻柱截面参数之后,控制终端可从每个钻柱截面参数中筛选出截面的直径,接着将该每个截面的直径依次输入至ICEM CFD软件中,以由该ICEM CFD软件根据每个截面的直径建立相应的二维数值模型。其中,二维数值模型可理解为矩形区域(尺寸可以但不局限于为50D*70D,D为截面的直径,且阻塞率可以但不局限于设置为2%),在该矩形区域中设置有与截面的直径一致的截面(用直径D为该截面的直径进行表示)以及与该截面连接的弹簧阻尼系统。此处,每个截面均可模拟发生横向运动或纵向运动或横向运动以及纵向运动,该横向运动可理解为截面的中心沿着与矩形区域横边平行的轴线向左或是向右运动,且该截面的中心设置在与矩形区域横边平行的轴线上;该纵向运动可理解为截面的中心沿着与矩形区域竖边平行的轴线向上或是向下运动,且该截面的中心设置在与矩形区域竖边平行的轴线上。
进一步的,控制终端在ICEM CFD软件中建立出与每个截面的直径所对应的二维数值模型之后,可基于该ICEM CFD软件对每个二维数值模型分别进行划分处理,以在每个二维数值模型中确定出刚性区域、动网格区域以及静网格区域。
作为本申请实施例的一种可选,根据预设的距离区间,对每个二维数值模型进行划分处理,包括:
在每个二维数值模型中确定出三维钻柱的截面中心;
将与三维钻柱的截面中心之间的距离处于预设的第一距离区间的区域划分为刚性区域;
将与三维钻柱的截面中心之间的距离处于预设的第二距离区间的区域划分为动网格区域;
将与三维钻柱的截面中心之间的距离处于预设的第三距离区间的区域划分为静网格区域;
对经过划分处理后的每个二维数值模型进行标记处理。
具体地,控制终端可先在每个二维数值模型中确定出三维钻柱的截面中心,接着可按照预设的第一距离区间,将与该三维钻柱的截面中心之间的距离处于预设的第一距离区间的区域划分为刚性区域,该预设的第一距离区间可以但不局限于表示为[D,L1],D为截面的直径。接着按照预设的第二距离区间,将与该三维钻柱的截面中心之间的距离处于预设的第二距离区间的区域划分为动网格区域,该预设的第二距离区间可以但不局限于表示为[L1,L2],D为截面的直径,其中,截面可在该动网格区域内进行运动。接着,按照预设的第三距离区间,将与该三维钻柱的截面中心之间的距离处于预设的第三距离区间的区域划分为静网格区域,该预设的第三距离区间可以但不局限于表示为[L2,L3],D为截面的直径。此处,L1小于L2,L2小于L3。
可以理解的是,在对每个二维数值模型进行区域划分处理之后,为了实现动网格模型,还可对该经过区域划分处理后的每个二维数值模型进行网格划分处理,网格类型可以但不局限于为三角形网格,并可由控制终端将与该经过网格划分处理后的每个二维数值模型所对应的文件导入至Fluent软件。
此处还可参阅图2示出的本申请实施例提供的一种二维数值模型的划分效果示意图。如图2所示,在每个二维数值模型中包括截面、截面壁面以及边界层部分区域被划分为刚性区域,截面附近区域被划分为动网格区域,其他所有区域被划分为静网格区域,且每个区域全部被划分为多个三角形网格。
当然,在本申请实施例中ICEM CFD软件还可直接建立每个二维数值模型,对该每个二维数值模型进行划分处理,并可在控制终端将该经过划分处理后的每个二维数值模型型所对应的文件导入至Fluent软件之后,由该Fluent软件结合每个单元节点所对应的截面的直径对该圆柱绕流数值模型进行调整,此处不限定于此。
进一步的,为了便于区分各个经过划分处理后的二维数值模型,还可在对每个二维数值模型进行划分处理之后,控制终端基于ICEM CFD软件依次对每个经过划分处理后的二维数值模型进行标号,例如第一个经过划分处理后的二维数值模型标号为1,第二个经过划分处理后的二维数值模型标号为2,并以此类推。
步骤106、从经过划分处理后的每个二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个流体区域的流场特性。
具体地,在对每个二维数值模型进行划分处理,并将每个经过划分处理后的二维数值模型由ICEM CFD软件导入至Fluent软件之后,控制终端可在基于Fluent软件所导出的每个二维数值模型中将除去刚性区域的其他所有区域作为初始流体区域,并结合预先设定在Fluent软件中的边界条件,对该初始流体区域进行裁剪处理,以得到目标流体区域。其中,边界条件可以但不局限于包括左边界条件、右边界条件以及上下侧边界,该左边界条件为速度入口(velocity-inlet),该右边界条件为压力出口(pressure-out),该上下侧边界为对称边界(symmetry)。
此处还可参阅图3示出的本申请实施例提供的又一种二维数值模型的划分效果示意图。如图3所示,每个二维数值模型中包括有截面、横向设置的与该截面连接的弹簧阻尼系统以及竖向设置的与该截面连接的弹簧阻尼系统,该二维数值模型的左侧对应为速度入口,该二维数值模型的右侧对应为压力出口,该二维数值模型的上下侧对应为对称边界,且该对称边界到截面的中心距离保持一致(图式为25D,D为截面的直径)。在图3中,速度入口侧与截面的中心之间的距离为20D,自由出流侧与截面的中心之间的距离为50D。
进一步的,在得到每个二维数值模型的目标流体区域之后,控制终端可利用Fluent软件,结合每个目标流体区域以及预设的流体参数,计算出每个流体区域的流场特性。其中,预设的流体参数可以但不局限于由用户预先输入至该Fluent软件中,该预设的流体参数具体可包括流体的密度以及运动粘度,此处所模拟的流体对象一般可设置为水。可以理解的是,在基于Fluent软件计算每个流体区域的流场特性的过程中,具体可由Fluent软件基于有限体积法把每个目标流体区域的流场离散成有限个节点,接着可用RANS方法对流体对应的雷诺平均N-S方程进行求解,以得到相应的流场特性,且该流场特性可以但不局限于包括整个流场中流体的速度、压力等信息。需要说明的是,在本申请实施例中应用Fluent软件计算流体区域的流场特性可为本领域常见的技术手段,此处不过多赘述。
当然,在基于Fluent软件计算每个流体区域的流场特性之前,控制终端还可将预先设定的用于计算每个流场特性的求解参数输入至Fluent软件中,以便于该Fluent软件结合每个目标流体区域、预设的流体参数以及求解参数,快速准确的计算出流体区域的流场特性。其中,该求解参数可以但不局限于包括求解算法、离散格式以及求解精度中任意至少一种。
还可以理解的是,在计算出每个流体区域的流场特性之前,控制终端还可控制Fluent软件结合预设的流体参数以及每个目标流体区域,对每个目标流体区域进行初始化,以保障模拟出的每个目标流体区域的可靠性。其中,初始化操作可以但不局限于由控制终端根据自动控制程序对该Fluent软件进行操作,此处不限定于此。
步骤108、基于所有流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数。
具体地,在计算出每个流体区域的流场特性之后,控制终端可利用Fluent软件中所嵌入的UDF程序,通过调用宏Compute_Force_And_Moment,从该每个流体区域的流场特性中提取出三维钻柱的表面压力载荷(也即理解为流体施加在设置在三维钻柱的每个单元节点上的压力分布数据),接着可从该每个三维钻柱的表面压力载荷中按照不同流向分别划分出处于顺流向(也可理解为竖流向)的压力载荷以及处于横流向的压力载荷。
进一步的,控制终端可对所有的处于顺流向的压力载荷进行汇总处理,例如但不局限于将其代入至预设的顺流向载荷矩阵中,以得到三维钻柱的顺流向载荷矩阵;接着,还可对所有处于横流向的压力载荷进行汇总处理,例如但不局限于将其代入至预设的横流向载荷矩阵中,以得到三维钻柱的横流向载荷矩阵。
作为本申请实施例的又一种可选,控制终端在计算三维钻柱的速度参数的过程中,可利用Fluent软件中所嵌入的UDF程序,通过预设的三维钻柱的动力学参数以及所有的钻柱截面参数,分别计算出三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵,其中,该三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵可结合上述提到的Euler-Bernoulli梁模型,理解为与该Euler-Bernoulli梁模型中的三维钻柱所对应的质量矩阵、刚度矩阵以及阻尼矩阵,且该质量矩阵可以但不局限于表示为,该刚度矩阵可以但不局限于表示为/>,该阻尼矩阵可以但不局限于表示为/>。
此处,三维钻柱的刚度矩阵包括有弹性刚度矩阵/>以及几何刚度矩阵/>,该弹性刚度矩阵/>的单元弹性刚度矩阵/>可表示为:
该几何刚度矩阵的单元几何刚度矩阵/>可表示为:
三维钻柱的刚度矩阵的单元刚度矩阵可结合上述提到的单元弹性刚度矩阵以及单元几何刚度矩阵/>,表示为:
三维钻柱的质量矩阵的单元质量矩阵/>可表示为:/>
在上述表达式中,E可对应为三维钻柱的弹性模量,A可对应为每个钻柱截面参数中的截面面积,可对应为每个钻柱截面参数中的截面相对Y轴的惯性矩,/>可对应为每个钻柱截面参数中的截面相对Z轴的惯性矩,/>可对应为每个钻柱截面参数中的截面相对X轴的极惯性矩,/>可对应为每两个相邻的单元节点之间的长度,/>可对应为三维钻柱的材料密度。
需要说明的是,在本申请实施例中三维钻柱的阻尼矩阵采用瑞利阻尼形式,其与质量矩阵和刚度矩阵之间存在如下的线性关系(也即通过下式得到三维钻柱的阻尼矩阵):
参数a以及参数b可以但不局限于通过下式得到:
上式中,以及/>可对应为三维钻柱的振动频率,/>可对应为三维钻柱的阻尼比,且为了简化计算,该/>以及/>可分别为三维钻柱的一阶固有振动频率以及二阶固有振动频率。
还需要说明的是,预设的三维钻柱的动力学参数可以但不局限于包括如下所示的参数表中的所有参数:
三维钻柱的动力学参数表
/>
进一步的,在分别计算出三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵之后,控制终端还可利用Fluent软件中所嵌入的UDF程序,将三维钻柱的顺流向载荷矩阵、三维钻柱的横流向载荷矩阵、三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵代入至预设的钻柱动力方程,并基于预设的求解参数对钻柱动力方程进行计算,得到三维钻柱的速度参数。其中,预设的钻柱动力方程可以但不局限于表示如下:
上式中,、/>以及/>可分别对应为三维钻柱处于顺流向的加速度向量、速度向量以及位移向量,/>、/>以及/>可分别对应为三维钻柱处于横流向的加速度向量、速度向量以及位移向量,/>可对应为三维钻柱的顺流向载荷矩阵,/>可对应为三维钻柱的横流向载荷矩阵。
进一步的,控制终端还可利用Fluent软件中所嵌入的UDF程序,基于预设的法,以顺流向的速度参数为例,可建立由t时刻三维钻柱处于顺流向的加速度向量、速度向量以及位移向量到/>时刻三维钻柱处于顺流向的加速度向量、速度向量以及位移向量的递推关系,可以但不局限于表示如下:
接着,假设在时刻有:
其中,上式中,/>
联立可得时刻的位移为:
接着,可得到:
接着,时刻三维钻柱的加速度以及速度可分别表示为:
其中,,/>,/>,/>,/>,,/>,/>,/>对应为0.5,/>对应为0.25。
可以理解的是,预设的求解参数可以但不局限于包括湍流模型、求解算法、离散格式、求解精度以及时间步长,该三维钻柱的速度参数可以但不局限于包括分别处于顺流向的位移、速度、加速度以及处于横流向的位移、速度以及加速度。
步骤110、根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。
具体地,在计算出三维钻柱的速度参数之后,控制终端可利用Fluent软件中所嵌入的UDF程序,调用宏DEFINE_CG_MOTION,将该三维钻柱的速度传递给Fluent软件中的每个二维数值模型,并通过宏DEFINE_CG_MOTION赋予每个二维数值模型中预设的旋转角速度,以使该二维数值模型中的三维钻柱可按照速度以及旋转角速度进行模拟运动,进而实现三维钻柱边界的位移和旋转。
其中,控制终端可按照预设间隔时刻,依次计算出与每个间隔时刻对应的三维钻柱的运动参数,其中,三维钻柱的运动参数可包括升力系数、阻力系数、横流向振幅以及顺流向振幅等参数。
可以理解的是,控制终端还可在Fluent软件中预先选择三维钻柱的运动参数的具体类型,以便于该Fluent软件可结合每个二维数值模型中三维钻柱的运动情况计算出相应类型的结果,例如控制终端可在Fluent软件中预先设置提取升力系数()、阻力系数(/>)、横流向振幅(/>)和顺流向振幅(/>)等计算数据。
可能的,在得到与每个间隔时刻对应的三维钻柱的运动参数之后,控制终端可根据该与每个间隔时刻对应的三维钻柱的运动参数中的振幅参数绘制出振幅曲线,并判断该振幅曲线是否处于预设的振幅区间。当检测到该振幅曲线处于预设的振幅区间时,确定三维钻柱的运动参数满足预设条件,并输出该三维钻柱的运动参数。
或是,控制终端还可在振幅曲线中判断是否存在至少两个振幅值相同的点,例如可在振幅曲线中规划任意一条处于该振幅最小值至振幅最大值之间的直线,该直线中所有点的纵坐标均保持一致,并判断该直线与振幅曲线是否存在至少两个交点。当检测到振幅曲线中存在至少两个振幅值相同的点时,确定三维钻柱的运动参数满足预设条件,并输出三维钻柱的运动参数。
作为本申请实施例的又一种可选,当检测到三维钻柱的运动参数不满足预设条件时,表明该三维钻柱的运动参数仍存在比较大的波动,则可结合当前时刻的三维钻柱的运动参数对每个流体区域进行更新处理。其中,控制终端还可在Fluent软件的Dynamic Mesh窗口中预先选择弹簧光顺模型(Smoothing)和网格重构模型(Remeshing),以设置动网格模型,便于每个二维数值模型的更新。
进一步的,控制终端可重复上述步骤,以计算出新的三维钻柱的运动参数,并继续判断该三维钻柱的运动参数是否满足预设条件。其中,计算出新的三维钻柱的运动参数的过程具体为:
基于预设的流体参数,计算出更新处理后的每个流体区域的流场特性;
基于所有更新处理后的流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数;
根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。
可以理解的是,在使该圆柱绕流数值模型中的圆柱可按照速度以及旋转角速度进行模拟运动之后,控制终端还可以但不局限于通过Fluent软件导出相应的case文件和data文件,并可运用Tecplot等后处理软件分析绘制出涡量云图和速度云图,或是还可运用Origin等后处理软件分析绘制振幅、流体力系数趋势图,以便于后续为海洋无立管圆柱的涡激振动机理和响应特性分析提供基础,并为无立管海洋圆柱的研发设计和无立管海洋圆柱涡激振动抑制技术的发展提供了理论指导。
此处可参阅图4示出的本申请实施例提供的一种三维钻柱的涡激振动流场涡量等值线示意图。如图4所示,三维钻柱的涡激振动流场涡量等值线中,z可理解为单元节点在三维钻柱上的位置,L可理解为三维钻柱的高度,在两端处(也即z/L=0或1)为标准的2S模式(一个振动周期内三维钻柱表面脱落两次漩涡),且漩涡在横流方向上间距很小;在中段处(也即z/L=0.4或0.5)漩涡在横流方向的间距增大,并且规律性下降。
此处还可参阅图5示出的本申请实施例提供的一种三维钻柱的涡激振动振幅示意图。如图5所示,上半部分为三维钻柱涡激振动横流向振幅图,该图中横坐标对应为z/L,z可理解为单元节点在三维钻柱上的位置,L可理解为三维钻柱的高度,纵坐标对应为振幅值,可以明显看出,三维钻柱从初始位置到最大位移处的过程中,该三维钻柱的振动模态从一阶变成二阶再变为一阶。下半部分为三维钻柱涡激振动顺流向振幅图,该图中横坐标对应为z/L,z可理解为单元节点在三维钻柱上的位置,L可理解为三维钻柱的高度,纵坐标对应为振幅值,可以明显看出,由于水流产生的阻力,三维钻柱始终沿顺流方向弯曲,且该三维钻柱在处于顺流向的运动方式与处于横流向的运动方式存在明显差异。
请参阅图6,图6示出了本申请实施例提供的一种三维钻柱的涡激振动计算装置的结构示意图。
如图6所示,该三维钻柱的涡激振动计算装置至少可以包括第一处理模块601、第二处理模块602、第一计算模块603、第二计算模块604以及参数输出模块605,其中:
第一处理模块601,用于基于预设的单元节点个数将三维钻柱切分为至少两个二维流场平面,并确定出与每个二维流场平面对应的钻柱截面参数;其中,每个单元节点对应于一个二维流场平面;
第二处理模块602,用于基于每个钻柱截面参数建立相应的二维数值模型,并根据预设的距离区间,对每个二维数值模型进行划分处理;
第一计算模块603,用于从经过划分处理后的每个二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个流体区域的流场特性;
第二计算模块604,用于基于所有流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数;
参数输出模块605,用于根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。
在一些可能的实施例中,第二处理模块用于:
在每个二维数值模型中确定出三维钻柱的截面中心;
将与三维钻柱的截面中心之间的距离处于预设的第一距离区间的区域划分为刚性区域;
将与三维钻柱的截面中心之间的距离处于预设的第二距离区间的区域划分为动网格区域;
将与三维钻柱的截面中心之间的距离处于预设的第三距离区间的区域划分为静网格区域;
对经过划分处理后的每个二维数值模型进行标记处理;
其中,预设的第一距离区间中的最大值等于预设的第二距离区间中的最小值,预设的第二距离区间中的最大值等于预设的第三距离区间的最小值。
在一些可能的实施例中,第一计算模块用于:
将经过划分处理后的每个二维数值模型中除去刚性区域的所有区域作为每个二维数值模型的初始流体区域;
基于预设的边界条件,对每个二维数值模型的初始流体区域进行裁剪处理,得到每个二维数值模型的目标流体区域;
根据每个二维数值模型的目标流体区域以及预设的流体参数,计算出每个流体区域的流场特性。
在一些可能的实施例中,第二计算模块用于:
从每个流体区域的流场特性中分别提取出三维钻柱的表面压力载荷;
对每个三维钻柱的表面压力载荷中处于顺流向的压力载荷进行汇总处理,得到三维钻柱的顺流向载荷矩阵;
对每个三维钻柱的表面压力载荷中处于横流向的压力载荷进行汇总处理,得到三维钻柱的横流向载荷矩阵。
在一些可能的实施例中,第二计算模块还用于:
根据所有钻柱截面参数分别计算出三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵;
将三维钻柱的顺流向载荷矩阵、三维钻柱的横流向载荷矩阵、三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵代入至预设的钻柱动力方程,并基于预设的求解参数对钻柱动力方程进行计算,得到三维钻柱的速度参数。
在一些可能的实施例中,参数输出模块用于:
基于三维钻柱的运动参数绘制振幅曲线,并判断振幅曲线是否处于预设的振幅区间;
当检测到振幅曲线处于预设的振幅区间时,确定三维钻柱的运动参数满足预设条件,并输出三维钻柱的运动参数;或
基于三维钻柱的运动参数绘制振幅曲线,并判断振幅曲线中是否存在至少两个振幅值相同的点;
当检测到振幅曲线中存在至少两个振幅值相同的点时,确定三维钻柱的运动参数满足预设条件,并输出三维钻柱的运动参数。
在一些可能的实施例中,装置还包括:
在根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数之后,当检测到三维钻柱的运动参数不满足预设条件时,基于三维钻柱的运动参数对每个流体区域进行更新处理;
基于预设的流体参数,计算出更新处理后的每个流体区域的流场特性;
基于所有更新处理后的流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数;
根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。
本领域的技术人员可以清楚地了解到本申请实施例的技术方案可借助软件和/或硬件来实现。本说明书中的“单元”和“模块”是指能够独立完成或与其他部件配合完成特定功能的软件和/或硬件,其中硬件例如可以是现场可编程门阵列(Field-ProgrammableGate Array,FPGA)、集成电路(Integrated Circuit,IC)等。
请参阅图7,图7示出了本申请实施例提供的又一种三维钻柱的涡激振动计算装置的结构示意图。
如图7所示,该三维钻柱的涡激振动计算装置700可以包括至少一个处理器701、至少一个网络接口704、用户接口703、存储器705以及至少一个通信总线702。
其中,通信总线702可用于实现上述各个组件的连接通信。
其中,用户接口703可以包括按键,可选用户接口还可以包括标准的有线接口、无线接口。
其中,网络接口704可以但不局限于包括蓝牙模块、NFC模块、Wi-Fi模块等。
其中,处理器701可以包括一个或者多个处理核心。处理器701利用各种接口和线路连接三维钻柱的涡激振动计算装置700内的各个部分,通过运行或执行存储在存储器705内的指令、程序、代码集或指令集,以及调用存储在存储器705内的数据,执行路由三维钻柱的涡激振动计算装置700的各种功能和处理数据。可选的,处理器701可以采用DSP、FPGA、PLA中的至少一种硬件形式来实现。处理器701可集成CPU、GPU和调制解调器等中的一种或几种的组合。其中,CPU主要处理操作系统、用户界面和应用程序等;GPU用于负责显示屏所需要显示的内容的渲染和绘制;调制解调器用于处理无线通信。可以理解的是,上述调制解调器也可以不集成到处理器701中,单独通过一块芯片进行实现。
其中,存储器705可以包括RAM,也可以包括ROM。可选的,该存储器705包括非瞬时性计算机可读介质。存储器705可用于存储指令、程序、代码、代码集或指令集。存储器705可包括存储程序区和存储数据区,其中,存储程序区可存储用于实现操作系统的指令、用于至少一个功能的指令(比如触控功能、声音播放功能、图像播放功能等)、用于实现上述各个方法实施例的指令等;存储数据区可存储上面各个方法实施例中涉及到的数据等。存储器705可选的还可以是至少一个位于远离前述处理器701的存储装置。如图7所示,作为一种计算机存储介质的存储器705中可以包括操作系统、网络通信模块、用户接口模块以及三维钻柱的涡激振动计算应用程序。
具体地,处理器701可以用于调用存储器705中存储的三维钻柱的涡激振动计算应用程序,并具体执行以下操作:
基于预设的单元节点个数将三维钻柱切分为至少两个二维流场平面,并确定出与每个二维流场平面对应的钻柱截面参数;其中,每个单元节点对应于一个二维流场平面;
基于每个钻柱截面参数建立相应的二维数值模型,并根据预设的距离区间,对每个二维数值模型进行划分处理;
从经过划分处理后的每个二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个流体区域的流场特性;
基于所有流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数;
根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。
在一些可能的实施例中,根据预设的距离区间,对每个二维数值模型进行划分处理,包括:
在每个二维数值模型中确定出三维钻柱的截面中心;
将与三维钻柱的截面中心之间的距离处于预设的第一距离区间的区域划分为刚性区域;
将与三维钻柱的截面中心之间的距离处于预设的第二距离区间的区域划分为动网格区域;
将与三维钻柱的截面中心之间的距离处于预设的第三距离区间的区域划分为静网格区域;
对经过划分处理后的每个二维数值模型进行标记处理;
其中,预设的第一距离区间中的最大值等于预设的第二距离区间中的最小值,预设的第二距离区间中的最大值等于预设的第三距离区间的最小值。
在一些可能的实施例中,从经过划分处理后的每个二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个流体区域的流场特性,包括:
将经过划分处理后的每个二维数值模型中除去刚性区域的所有区域作为每个二维数值模型的初始流体区域;
基于预设的边界条件,对每个二维数值模型的初始流体区域进行裁剪处理,得到每个二维数值模型的目标流体区域;
根据每个二维数值模型的目标流体区域以及预设的流体参数,计算出每个流体区域的流场特性。
在一些可能的实施例中,基于所有流体区域的流场特性确定出三维钻柱的受力参数,包括:
从每个流体区域的流场特性中分别提取出三维钻柱的表面压力载荷;
对每个三维钻柱的表面压力载荷中处于顺流向的压力载荷进行汇总处理,得到三维钻柱的顺流向载荷矩阵;
对每个三维钻柱的表面压力载荷中处于横流向的压力载荷进行汇总处理,得到三维钻柱的横流向载荷矩阵。
在一些可能的实施例中,根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数,包括:
根据所有钻柱截面参数分别计算出三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵;
将三维钻柱的顺流向载荷矩阵、三维钻柱的横流向载荷矩阵、三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵代入至预设的钻柱动力方程,并基于预设的求解参数对钻柱动力方程进行计算,得到三维钻柱的速度参数。
在一些可能的实施例中,当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数,包括:
基于三维钻柱的运动参数绘制振幅曲线,并判断振幅曲线是否处于预设的振幅区间;
当检测到振幅曲线处于预设的振幅区间时,确定三维钻柱的运动参数满足预设条件,并输出三维钻柱的运动参数;或
基于三维钻柱的运动参数绘制振幅曲线,并判断振幅曲线中是否存在至少两个振幅值相同的点;
当检测到振幅曲线中存在至少两个振幅值相同的点时,确定三维钻柱的运动参数满足预设条件,并输出三维钻柱的运动参数。
在一些可能的实施例中,在根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数之后,还包括:
当检测到三维钻柱的运动参数不满足预设条件时,基于三维钻柱的运动参数对每个流体区域进行更新处理;
基于预设的流体参数,计算出更新处理后的每个流体区域的流场特性;
基于所有更新处理后的流体区域的流场特性确定出三维钻柱的受力参数,并根据三维钻柱的受力参数、所有钻柱截面参数以及预设的钻柱动力方程,计算出三维钻柱的速度参数;
根据三维钻柱的速度参数以及预设的旋转角速度,得到三维钻柱的运动参数,并当检测到三维钻柱的运动参数满足预设条件时,输出三维钻柱的运动参数。
本申请还提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述方法的步骤。其中,计算机可读存储介质可以包括但不限于任何类型的盘,包括软盘、光盘、DVD、CD-ROM、微型驱动器以及磁光盘、ROM、RAM、EPROM、EEPROM、DRAM、VRAM、闪速存储器设备、磁卡或光卡、纳米系统(包括分子存储器IC),或适合于存储指令和/或数据的任何类型的媒介或设备。
需要说明的是,对于前述的各方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本申请并不受所描述的动作顺序的限制,因为依据本申请,某些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和模块并不一定是本申请所必须的。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置,可通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些服务接口,装置或单元的间接耦合或通信连接,可以是电性或其它的形式。
作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储器中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储器中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本申请各个实施例方法的全部或部分步骤。而前述的存储器包括:U盘、只读存储器(Read-Only Memory, ROM)、随机存取存储器(RandomAccess Memory,RAM)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通进程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储器中,存储器可以包括:闪存盘、只读存储器(Read-Only Memory, ROM)、随机存取器(Random AccessMemory,RAM)、磁盘或光盘等。
以上者,仅为本公开的示例性实施例,不能以此限定本公开的范围。即但凡依本公开教导所作的等效变化与修饰,皆仍属本公开涵盖的范围内。本领域技术人员在考虑说明书及实践这里的公开后,将容易想到本公开的其实施方案。本申请旨在涵盖本公开的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未记载的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本公开的范围和精神由权利要求限定。
Claims (8)
1.一种三维钻柱的涡激振动计算方法,其特征在于,包括:
基于预设的单元节点个数将三维钻柱切分为至少两个二维流场平面,并确定出与每个所述二维流场平面对应的钻柱截面参数;其中,每个所述单元节点对应于一个所述二维流场平面;基于每个所述钻柱截面参数建立相应的二维数值模型,并根据预设的距离区间,对每个所述二维数值模型进行划分处理;
从经过划分处理后的每个所述二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个所述流体区域的流场特性;
基于所有所述流体区域的流场特性确定出所述三维钻柱的受力参数,并根据所述三维钻柱的受力参数、所有所述钻柱截面参数以及预设的钻柱动力方程,计算出所述三维钻柱的速度参数;
根据所述三维钻柱的速度参数以及预设的旋转角速度,得到所述三维钻柱的运动参数,并当检测到所述三维钻柱的运动参数满足预设条件时,输出所述三维钻柱的运动参数;
其中,所述基于所有所述流体区域的流场特性确定出所述三维钻柱的受力参数,包括:
从每个所述流体区域的流场特性中分别提取出所述三维钻柱的表面压力载荷;
对每个所述三维钻柱的表面压力载荷中处于顺流向的压力载荷进行汇总处理,得到所述三维钻柱的顺流向载荷矩阵;
对每个所述三维钻柱的表面压力载荷中处于横流向的压力载荷进行汇总处理,得到所述三维钻柱的横流向载荷矩阵;
其中,所述根据所述三维钻柱的受力参数、所有所述钻柱截面参数以及预设的钻柱动力方程,计算出所述三维钻柱的速度参数,包括:
根据所有所述钻柱截面参数分别计算出所述三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵;将所述三维钻柱的顺流向载荷矩阵、所述三维钻柱的横流向载荷矩阵、所述三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵代入至预设的钻柱动力方程,并基于预设的求解参数对所述钻柱动力方程进行计算,得到所述三维钻柱的速度参数。
2.根据权利要求1所述的方法,其特征在于,所述根据预设的距离区间,对每个所述二维数值模型进行划分处理,包括:
在每个所述二维数值模型中确定出所述三维钻柱的截面中心;
将与所述三维钻柱的截面中心之间的距离处于预设的第一距离区间的区域划分为刚性区域;将与所述三维钻柱的截面中心之间的距离处于预设的第二距离区间的区域划分为动网格区域;
将与所述三维钻柱的截面中心之间的距离处于预设的第三距离区间的区域划分为静网格区域;
对经过划分处理后的每个所述二维数值模型进行标记处理;
其中,所述预设的第一距离区间中的最大值等于所述预设的第二距离区间中的最小值,所述预设的第二距离区间中的最大值等于所述预设的第三距离区间的最小值。
3.根据权利要求2所述的方法,其特征在于,所述从经过划分处理后的每个所述二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个所述流体区域的流场特性,包括:
将经过划分处理后的每个所述二维数值模型中除去所述刚性区域的所有区域作为每个所述二维数值模型的初始流体区域;
基于预设的边界条件,对每个所述二维数值模型的初始流体区域进行裁剪处理,得到每个所述二维数值模型的目标流体区域;
根据每个所述二维数值模型的目标流体区域以及预设的流体参数,计算出每个所述流体区域的流场特性。
4.根据权利要求1所述的方法,其特征在于,所述当检测到所述三维钻柱的运动参数满足预设条件时,输出所述三维钻柱的运动参数,包括:
基于所述三维钻柱的运动参数绘制振幅曲线,并判断所述振幅曲线是否处于预设的振幅区间;
当检测到所述振幅曲线处于所述预设的振幅区间时,确定所述三维钻柱的运动参数满足预设条件,并输出所述三维钻柱的运动参数;或
基于所述三维钻柱的运动参数绘制振幅曲线,并判断所述振幅曲线中是否存在至少两个振幅值相同的点;
当检测到所述振幅曲线中存在至少两个振幅值相同的点时,确定所述三维钻柱的运动参数满足预设条件,并输出所述三维钻柱的运动参数。
5.根据权利要求1所述的方法,其特征在于,在所述根据所述三维钻柱的速度参数以及预设的旋转角速度,得到所述三维钻柱的运动参数之后,还包括:
当检测到所述三维钻柱的运动参数不满足所述预设条件时,基于所述三维钻柱的运动参数对每个所述流体区域进行更新处理;
基于预设的流体参数,计算出更新处理后的每个所述流体区域的流场特性;
基于所有更新处理后的所述流体区域的流场特性确定出所述三维钻柱的受力参数,并根据所述三维钻柱的受力参数、所有所述钻柱截面参数以及预设的钻柱动力方程,计算出所述三维钻柱的速度参数;
根据所述三维钻柱的速度参数以及预设的旋转角速度,得到所述三维钻柱的运动参数,并当检测到所述三维钻柱的运动参数满足预设条件时,输出所述三维钻柱的运动参数。
6.一种三维钻柱的涡激振动计算装置,其特征在于,包括:
第一处理模块,用于基于预设的单元节点个数将三维钻柱切分为至少两个二维流场平面,并确定出与每个所述二维流场平面对应的钻柱截面参数;其中,每个所述单元节点对应于一个所述二维流场平面;
第二处理模块,用于基于每个所述钻柱截面参数建立相应的二维数值模型,并根据预设的距离区间,对每个所述二维数值模型进行划分处理;
第一计算模块,用于从经过划分处理后的每个所述二维数值模型中确定出相应的流体区域,并基于预设的流体参数,计算出每个所述流体区域的流场特性;
第二计算模块,用于基于所有所述流体区域的流场特性确定出所述三维钻柱的受力参数,并根据所述三维钻柱的受力参数、所有所述钻柱截面参数以及预设的钻柱动力方程,计算出所述三维钻柱的速度参数;
参数输出模块,用于根据所述三维钻柱的速度参数以及预设的旋转角速度,得到所述三维钻柱的运动参数,并当检测到所述三维钻柱的运动参数满足预设条件时,输出所述三维钻柱的运动参数;
其中,所述基于所有所述流体区域的流场特性确定出所述三维钻柱的受力参数,包括:
从每个所述流体区域的流场特性中分别提取出所述三维钻柱的表面压力载荷;
对每个所述三维钻柱的表面压力载荷中处于顺流向的压力载荷进行汇总处理,得到所述三维钻柱的顺流向载荷矩阵;
对每个所述三维钻柱的表面压力载荷中处于横流向的压力载荷进行汇总处理,得到所述三维钻柱的横流向载荷矩阵;
其中,所述根据所述三维钻柱的受力参数、所有所述钻柱截面参数以及预设的钻柱动力方程,计算出所述三维钻柱的速度参数,包括:
根据所有所述钻柱截面参数分别计算出所述三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵;将所述三维钻柱的顺流向载荷矩阵、所述三维钻柱的横流向载荷矩阵、所述三维钻柱的质量矩阵、刚度矩阵以及阻尼矩阵代入至预设的钻柱动力方程,并基于预设的求解参数对所述钻柱动力方程进行计算,得到所述三维钻柱的速度参数。
7.一种三维钻柱的涡激振动计算装置,其特征在于,包括处理器以及存储器;
所述处理器与所述存储器连接;
所述存储器,用于存储可执行程序代码;
所述处理器通过读取所述存储器中存储的可执行程序代码来运行与所述可执行程序代码对应的程序,以用于执行如权利要求1-5任一项所述方法的步骤。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机可读存储介质中存储有指令,当所述指令在计算机或处理器上运行时,使得所述计算机或处理器执行如权利要求1-5任一项所述方法的步骤。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310056576 | 2023-01-19 | ||
CN2023100565769 | 2023-01-19 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116702648A CN116702648A (zh) | 2023-09-05 |
CN116702648B true CN116702648B (zh) | 2024-05-10 |
Family
ID=87830716
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310699786.XA Active CN116702648B (zh) | 2023-01-19 | 2023-06-14 | 一种三维钻柱的涡激振动计算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116702648B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117113777B (zh) * | 2023-09-08 | 2024-04-09 | 武汉理工大学 | 一种考虑内流的钻柱涡激振动计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110795827A (zh) * | 2019-10-11 | 2020-02-14 | 扬州大学 | 非线性能量阱作用下弹性支撑刚性柱体结构的涡激振动仿真方法 |
CN111008490A (zh) * | 2019-10-11 | 2020-04-14 | 扬州大学 | 采用非线性能量阱抑制双柱体结构流激振动的流固耦合预测方法 |
CN112084560A (zh) * | 2020-09-09 | 2020-12-15 | 长沙理工大学 | 一种多幅桥面桥梁涡激振动数值模拟计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107605666B (zh) * | 2017-09-11 | 2019-01-11 | 北京金风科创风电设备有限公司 | 具有抑制涡激振动功能的围护结构及抑制涡激振动的方法 |
-
2023
- 2023-06-14 CN CN202310699786.XA patent/CN116702648B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110795827A (zh) * | 2019-10-11 | 2020-02-14 | 扬州大学 | 非线性能量阱作用下弹性支撑刚性柱体结构的涡激振动仿真方法 |
CN111008490A (zh) * | 2019-10-11 | 2020-04-14 | 扬州大学 | 采用非线性能量阱抑制双柱体结构流激振动的流固耦合预测方法 |
CN112084560A (zh) * | 2020-09-09 | 2020-12-15 | 长沙理工大学 | 一种多幅桥面桥梁涡激振动数值模拟计算方法 |
Non-Patent Citations (3)
Title |
---|
Numerical investigation of wake and flow-induced vibrations of a rotating cylinder in flow;Yanxu Bao等;《Ocean Engineering》;20220916;全文 * |
深水钻井隔水管涡激振动特性的数值模拟研究;王嘉松;蒋世全;许亮斌;;石油钻采工艺;20150120(第01期);全文 * |
附属圆柱与强制旋转对圆柱流致振动的联合作用研究;包燕旭等;《第三十三届全国水动力学研讨会论文集》;20221028;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116702648A (zh) | 2023-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Otoguro et al. | Space–time VMS flow analysis of a turbocharger turbine with isogeometric discretization: computations with time-dependent and steady-inflow representations of the intake/exhaust cycle | |
CN116702648B (zh) | 一种三维钻柱的涡激振动计算方法及装置 | |
Saud Afzal et al. | Three-dimensional numerical modeling of pier scour under current and waves using level-set method | |
CN107060746A (zh) | 一种复杂裂缝性油藏流动模拟的方法 | |
Volz et al. | Numerical modelling of non-cohesive embankment breach with the dual-mesh approach | |
Lee et al. | Solving the shallow water equations using 2d sph particles for interactive applications | |
CN105069184A (zh) | 一种基于浸入边界法的搅拌反应釜模拟方法 | |
Zou et al. | Response of a submerged floating tunnel subject to flow-induced vibration | |
CN116702649B (zh) | 一种旋转圆柱的涡激振动计算方法及装置 | |
JP2014081900A (ja) | シミュレーションプログラム、シミュレーション方法及びシミュレーション装置 | |
Moreira et al. | Dam spillways and the SPH method: Two case studies in Portugal | |
Bai et al. | 3D mechanical analysis of subsea manifold installation by drill pipe in deep water | |
Zhang et al. | Numerical study of distributed hydrodynamic forces on a circular heave plate by large-eddy simulations with volume of fluid method | |
CN107832524A (zh) | 一种确定储层改造体积的方法、装置及系统 | |
Frère et al. | Discontinuous Galerkin methodology for Large-Eddy Simulations of wind turbine airfoils | |
Zhang et al. | ICCM2016: The implementation of two-dimensional multi-block lattice Boltzmann method on GPU | |
Nizamani et al. | Analysis of loads caused by waves on the deck near the free surface of the offshore platform using computational fluid dynamics | |
Furnes | Flow in Pelton turbines | |
Duan et al. | A high-order flux reconstruction method for 3D mixed overset meshes | |
Zhang et al. | CFD modeling of propeller tip vortex over large distances | |
CN116702647B (zh) | 一种附体-旋转组合圆柱结构的涡激振动计算方法及装置 | |
CN116842864A (zh) | 一种近壁状态下的涡激振动计算方法及装置 | |
CN109505591B (zh) | 确定缝洞型油藏未充填溶洞渗透率界限的方法及系统 | |
Wu et al. | Research progress on dam-break floods | |
Christakis et al. | A volume of fluid numerical model for wave impacts at coastal structures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |