变张力细长柔性圆柱体涡激振动响应预测方法技术

本发明专利技术公开了一种变张力细长柔性圆柱体涡激振动响应预测方法，主要包括以下步骤：1)建立结构和流场互为耦合的振动控制方程；2)基于有限差分法对耦合的振动控制方程进行求解；3)基于分析数据，对实例进行计算分析。本发明专利技术针对变张力细长柔性圆柱体涡激振动响应预测进行了全面研究。建立一个完整的流体‑结构耦合模型，用于分析顶张式立管、钢悬链式立管等变张力柔性圆柱体涡激振动响应的预测问题。该模型实时考虑了结构与流场之间互为耦合的影响，可很好地预测变张力柔性圆柱体的振动响应特性。

Prediction of Vortex-Induced Vibration Response of Slender Flexible Cylinder with Variable Tension

【技术实现步骤摘要】
变张力细长柔性圆柱体涡激振动响应预测方法
本专利技术涉及细长柔性圆柱体涡激振动响应领域，具体地说，特别涉及一种变张力细长柔性圆柱体涡激振动响应预测方法。
技术介绍
近些年，伴随着海洋油气工程逐步向深海发展，出现了大量的超细长柔性结构，如：用于油气开发的各种形式的海洋立管(顶张式立管、钢悬链式立管等)以及用于浮式海洋平台(半潜式平台、立柱式平台等)定位的细长锚链。这种细长柔性结构在一定的洋流作用下，会在结构尾部形成交替脱落的漩涡，周期性的漩涡脱落会在结构上产生周期性的水动力载荷，从而诱发结构发生振动，称为“涡激振动”(Gao,Y,Fu,SX,Xiong,YM,Zhao,Y,andLiu,LM,Experimentalstudyonresponseperformanceofvortex-inducedvibrationofaflexiblecylinder,ShipsandOffshoreStructures,2017,12(1):116-134)。涡激振动响应是导致结构产生疲劳的主要因素之一，因此有必要对其加以研究。细长柔性圆柱体涡激振动响应的研究方法主要可分为实验方法以及数值方法。与数值方法相比，实验方法获取的数据更为可靠、得到的现象更为直观，但实验方法同样存在一些缺陷，比如：实验研究成本昂贵、实验结构尺寸受实验场地空间大小限制以及复杂流剖面难以模拟等。数值方法通常包括计算流体动力学(Computationalfluiddynamics,CFD)方法以及经验模型方法(高云，邹丽，宗智，基于有限差分法刚性圆柱体涡激振动响应特性数值研究，计算物理，2019，36(1)：53-59)。近些年，随着计算机计算速度的快速发展，CFD方法被广泛应用于预测圆柱体涡激振动响应中。由于对整个轴线方向的全尺度细长柔性圆柱体涡激振动响应进行数值模拟非常困难，目前绝大多数CFD研究还仅局限于对刚性圆柱体的涡激振动响应展开数值模拟。在解决细长柔性圆柱体涡激振动响应问题时，由于实验方法的高实验研究成本特性以及CFD方法的高数值研究成本特性，因此有必要建立一种快速预报细长柔性圆柱体涡激振动响应的方法，而尾流振子模型方法则正是近些年广泛应用于细长柔性圆柱体涡激振动响应预报的一种半经验模型方法。目前针对细长柔性圆柱体涡激振动响应的研究主要集中在恒张力细长柔性圆柱体，即：假定圆柱体的轴向张力沿轴线方向恒定不变。而实际海洋工程中的立管，如张紧式立管以及钢悬链式立管，其张力分布沿结构轴线方向呈变化趋势：从立管上端到立管下端其张力逐渐减小。对于这种变张力分布的细长柔性圆柱体，其涡激振动响应还有待进一步展开研究。
技术实现思路
为了解决现有技术的问题，本专利技术实施例提供了一种变张力细长柔性圆柱体涡激振动响应预测方法。所述技术方案如下：本专利技术提供了一种变张力细长柔性圆柱体涡激振动响应预测方法，主要包括以下步骤：1)建立结构和流场互为耦合的振动控制方程；2)基于有限差分法对耦合的振动控制方程进行求解；3)基于分析数据，对实例进行计算分析。进一步地，所述步骤1)具体包括：a、建立均匀流中变张力细长柔性圆柱体模型b、在所述均匀流中变张力细长柔性圆柱体模型中，考虑一变张力细长柔性圆柱体在横流方向的涡激振动响应问题；假设结构长度为L，直径为D，结构两端采用铰接边界条件，流场为速度大小为U的均匀来流；坐标系X方向取顺流方向，Z方向为铅直方向，Y方向为横流向振动方向；柔性圆柱体上受到的张力沿着轴线Z方向发生变化，记为Θ(Z)，最大张力以及最小张力分别位于圆柱体上端以及下端处，分别记为：Θmax以及Θmin；柔性圆柱体的弯曲刚度为EI。进一步地，将所述的细长柔性圆柱体看作变张力梁模型，建立如下振动方程：式(1)中T为时间，p(Z,T)为升力，表示为：p(Z,T)＝ρU2DCL0q(Z,T)/4，其中q(Z,T)表示尾部流场的运动，CL0为静止圆柱体横向升力系数；式(1)中系统单位长度质量m包括结构质量ms以及附加流体质量mf，表示为：m＝ms+mf,mf＝CMρD2π/4，其中ρ表示流体密度，CM表示流体附加质量系数；式(1)中阻尼R包括结构阻尼Rs和流体阻尼Rf，表示为：R＝Rs+Rf,Rf＝γΩfρD2，其中γ为黏滞力系数，与流体拖曳力系数CD的关系为：γ＝CD/(4πSt)，St数为斯脱哈尔数；对于海洋工程中的低质量阻尼比介质(如水)，通常结构阻尼与流体阻尼比起是小量，因此本专利技术忽略了结构阻尼的影响，只考虑了流体阻尼；Ωf为依据斯脱哈尔关系式得到的斯脱哈尔漩涡泄放频率，写作：Ωf＝2πStU/D；采用改进的Vanderpol方程来满足尾流振子的非线性特性，表达式如下：式(2)中ε和A为实验测试得到的经验小参数；将方程(1)和(2)转换成无量纲形式，令:y＝Y/D,z＝Z/D,t＝T·Ωf，并将其代入到方程(1)和(2)中，整理得到结构与尾部流场耦合振动方程的无量纲形式如下：式(3)中，μ为质量比，表示为μ＝m/(ρD2)，无量纲参数M表示为：M＝CL0/(16π2St2μ)；式(3)中b为无量纲弯曲刚度，表示为：c(z)为无量纲张力，表示为本专利技术假设无量纲张力c(z)沿着轴线方向呈线性变化，且假定最大无量纲张力cmax和最小无量纲张力cmin之间的关系为cmax＝2cmin，这样仅需要通过平均张力cmean的大小便可确定整个轴线上的张力分布情况；由于平均张力cmean＝(cmax+cmin)/2，因此cmin＝2cmean/3，cmax＝4cmean/3；那么任意一点z处对应的无量纲张力表示为：c(z)＝cmin+(cmax-cmin)×z/(L/D)。进一步地，所述步骤2)具体如下：采用具备二阶精度的有限差分法对方程(3)进行求解，将结构无量纲总长度L/D划分为M段；将无量纲总时间ttotal划分为N段；因此，数值计算时空间步长Δz＝L/(D×M)；时间步长Δt＝ttotal/N；被划分后的M+1个空间点记为：z＝zi(i＝0,1,2,…,M)；被划分后的N+1时间点记为：t＝tj(j＝0,1,2,..,N)；假设tn时刻zm位置处参数y和q表示为以及因此式(3)中各偏导数项的二阶精度差分格式表示为：将式(4)代入式(3)，并整理得到：y的初始条件设为在整个轴线上圆柱体振动位移以及速度均为0，即：y＝0以及q的初始条件设为：q为一小幅扰动以及圆柱体两端采用铰接边界条件，即结构在两个端点处位移为0且弯矩为0，表示如下：由式(5)可以看出：当2≤m≤M-2时y直接求出；当m＝0以及m＝M时，则需结合位移边界条件对其加以求解，由两端处位移为0得到：当m＝1以及m＝M-1时，则需结合弯矩为0边界条件，由两端弯矩为0得到：联合式(8)以及式(5)得到m＝1以及m＝M-1时y的表达式如下：当tn+1时刻zm处的y已知，便依据式(9)的第二个表达式求得tn+1时刻zm处的q，依此类推对式(5)进行迭代求解。本专利技术实施例提供的技术方案带来的有益效果是：本专利技术针对变张力细长柔性圆柱体涡激振动响应预测进行了全面研究。建立一个完整的流体-结构耦合模型，用于分析顶张式立管、钢悬链式立管等变张力柔性圆柱体涡激振动响应的预测问题。该模型实时考虑了结构与流场之间互为耦合的影响，可本文档来自技高网...

【技术保护点】
1.一种变张力细长柔性圆柱体涡激振动响应预测方法，其特征在于，主要包括以下步骤：1)建立结构和流场互为耦合的振动控制方程；2)基于有限差分法对耦合的振动控制方程进行求解；3)基于分析数据，对实例进行计算分析。

【技术特征摘要】
1.一种变张力细长柔性圆柱体涡激振动响应预测方法，其特征在于，主要包括以下步骤：1)建立结构和流场互为耦合的振动控制方程；2)基于有限差分法对耦合的振动控制方程进行求解；3)基于分析数据，对实例进行计算分析。2.如权利要求1所述的方法，其特征在于，所述步骤1)具体包括：a、建立均匀流中变张力细长柔性圆柱体模型b、在所述均匀流中变张力细长柔性圆柱体模型中，考虑一变张力细长柔性圆柱体在横流方向的涡激振动响应问题；假设结构长度为L，直径为D，结构两端采用铰接边界条件，流场为速度大小为U的均匀来流；坐标系X方向取顺流方向，Z方向为铅直方向，Y方向为横流向振动方向；柔性圆柱体上受到的张力沿着轴线Z方向发生变化，记为Θ(Z)，最大张力以及最小张力分别位于圆柱体上端以及下端处，分别记为：Θmax以及Θmin；柔性圆柱体的弯曲刚度为EI。3.如权利要求2所述的方法，其特征在于，将所述的细长柔性圆柱体看作变张力梁模型，建立如下振动方程：式(1)中T为时间，p(Z,T)为升力，表示为：p(Z,T)＝ρU2DCL0q(Z,T)/4，其中q(Z,T)表示尾部流场的运动，CL0为静止圆柱体横向升力系数；式(1)中系统单位长度质量m包括结构质量ms以及附加流体质量mf，表示为：m＝ms+mf,mf＝CMρD2π/4，其中ρ表示流体密度，CM表示流体附加质量系数；式(1)中阻尼R包括结构阻尼Rs和流体阻尼Rf，表示为：R＝Rs+Rf,Rf＝γΩfρD2，其中γ为黏滞力系数，与流体拖曳力系数CD的关系为：γ＝CD/(4πSt)，St数为斯脱哈尔数；对于海洋工程中的低质量阻尼比介质(如水)，通常结构阻尼与流体阻尼比起是小量，因此本发明忽略了结构阻尼的影响，只考虑了流体阻尼；Ωf为依据斯脱哈尔关系式得到的斯脱哈尔漩涡泄放频率，写作：Ωf＝2πStU/D；采用改进的Vanderpol方程来满足尾流振子的非线性特性，表达式如下：式(2)中ε和A为实验测试得到的经验小参数；将方程(1)和(2)转换成无量纲形式，令:y＝Y/D,z＝Z/D,t＝T·Ωf，并将其代入到方...

【专利技术属性】
技术研发人员：高云，张壮壮，刘黎明，杨斌，邹丽，宗智，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：四川,51