CN102507084A - 一种尾流立管的时域涡激升力确定方法 - Google Patents
一种尾流立管的时域涡激升力确定方法 Download PDFInfo
- Publication number
- CN102507084A CN102507084A CN2011102913735A CN201110291373A CN102507084A CN 102507084 A CN102507084 A CN 102507084A CN 2011102913735 A CN2011102913735 A CN 2011102913735A CN 201110291373 A CN201110291373 A CN 201110291373A CN 102507084 A CN102507084 A CN 102507084A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- wake flow
- standpipe
- time
- vibration
- 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.)
- Granted
Links
Images
Landscapes
- Measuring Volume Flow (AREA)
Abstract
本发明涉及海洋深水立管的研究方法,具体涉及一种尾流立管的时域涡激升力确定方法。本发明考虑了上游立管涡街对尾流立管涡激升力幅值和尾流立管涡旋脱落频率的影响,并采用迭代方法计算时域的涡激升力,从而建立了一个尾流立管的流固耦合时域涡激升力计算方法,解决了尾流立管的时域涡激升力的计算问题,为尾流立管的涡激振动与疲劳设计分析提供了时域分析的途径和依据。
Description
技术领域
本发明涉及海洋深水立管的研究方法,具体涉及一种尾流立管的时域涡激升力确定方法。
背景技术
涡激升力是圆柱体受到的垂直于流体流动方向的一种交变流体荷载,是由圆柱体尾流处的涡旋泄放引起的。对于弹性圆柱体而言,涡激升力将使圆柱体产生垂直于流体流动方向的振动。而对于刚性圆柱体,涡激升力并不引起圆柱体振动。这就是说,不论圆柱体是否振动,涡激升力始终是存在的。
涡激升力的交变频率等于涡旋泄放的频率,而涡旋泄放的频率取决于流体流经圆柱体的速度大小,涡激升力的大小取决于流体流经圆柱体的速度。因此,对于刚性圆柱体来说,涡激升力的大小和频率仅仅与流体的流动速度有关。而对于弹性圆柱体,由于圆柱体本身的振动,不仅使得流体流经圆柱体的速度发生变化,而且使得流体在圆柱体的扰动下产生粘性阻力和附加质量力。因此,对于弹性圆柱体的涡激升力计算不能仅仅考虑流场的流速一个因素(传统涡激升力计算方法的缺点),而应考虑圆柱体振动的速度和加速度。
对于顺流向排列的两根立管,从上游立管脱落的涡旋在尾流处形成涡街,从而对尾流立管产生干扰。同时,尾流立管自身的涡旋泄放也受到上游立管涡街的影响而改变其频率和强度。因此,尾流立管的横流向涡激振动强度远远大于孤立的单根立管。这个现象已被国内外的专家学者注意到,并开展了相应的研究。图1是本发明发明人的研究结果,其中幅度大的曲线为尾流立管的涡激升力时程,虚线是单根立管的涡激升力时程,幅度小的曲线是上游立管的涡激升力时程。从图中可以明显地看出,尾流立管的涡激升力远远大于单个孤独立管和上游立管。这将使尾流立管产生大幅度的横流向振动,从而引起较大的疲劳损伤。但是,现有技术尚没有尾流立管涡激升力的计算方法。现行的立管涡激振动与疲劳设计均采用单根立管的涡激升力计算方法,因此,计算结果偏于不安全。
发明内容
本发明的目的在于针对现有技术的缺陷,提供一种尾流立管的时域涡激升力确定方法,为尾流立管的涡激振动与疲劳设计分析提供了时域分析的途径和依据。
本发明的技术方案如下:一种尾流立管的时域涡激升力确定方法,建立的涡激升力时域模型如下:
式中:FL--涡激升力;
CL--升力系数;
CD--拖曳力系数;
ρ--流体密度;
λ--上下游立管轴线之间的距离L与直径之比L/D;
D--尾流立管直径;
U--流速;
t--时间;
q,b--由试验确定的系数;
针对尾流立管顺流向振动速度、横流向振动速度和加速度,采用迭代方法对上式进行计算,得到尾流立管的时域涡激升力。
进一步,如上所述的尾流立管的时域涡激升力确定方法,其中,所述的采用迭代方法进行计算的具体步骤如下:
1)给定尾流立管顺流向振动速度、横流向振动速度和加速度及计算时间的初值:
tj=t0=0,
式中:j--时间步数,计算开始时j=0;
i--迭代次数,每个时间步开始时i=0;
2)计算给定流速下的约化速度:
式中:Vr--约化速度;
U--流体流速;
fn--尾流立管的固有频率;
D--尾流立管直径;
当5≤Vr≤7时,将tj,代入下式计算第j时间步内第i次迭代的脉动拖曳力:
在以上两个公式中,
FD--脉动拖曳力;
CD--拖曳力系数;
ρ--流体密度;
λ--上下游立管轴线之间的距离L与直径之比L/D;
D--尾流立管直径;
U--流体流速;
t--时间;
ü--深水立管顺流向振动加速度;
p,b--由试验确定的系数;
式中:m-尾流立管的质量;
c-尾流立管的阻尼系数;
k-尾流立管的弯曲刚度;
然后,重复步骤3)~5)的计算;
如果 则令:
n为第j时间步内计算脉动拖曳力的总迭代次数;
式中:m-尾流立管的质量;
c-尾流立管的阻尼系数;
k-尾流立管的弯曲刚度;
然后,重复步骤6)~8)的计算;
tj=tj+1=tj+Δt,
m为第j时间步内计算涡激升力的总迭代次数;
重复步骤3)~8)的计算,直至计算时长满足需要。
本发明的有益效果如下:本发明考虑了上游立管涡街对尾流立管涡激升力幅值和尾流立管涡旋脱落频率的影响,并采用迭代方法计算时域的涡激升力,从而建立了一个尾流立管的流固耦合时域涡激升力计算方法,解决了尾流立管的时域涡激升力的计算问题,为尾流立管的涡激振动与疲劳设计分析提供了时域分析的途径和依据。
附图说明
图1为顺流向排列两根立管涡激升力时程示意图;
图2为本发明的方法流程图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细描述。
顺流向排列的两根立管在定常流动流体的作用下将产生涡激振动,下游立管由于受到上游立管尾流的作用,其横流向振动的频率和幅值与孤立的单根立管有较大的区别。而现有技术只有单根立管的涡激升力计算方法,没有尾流立管的涡激升力计算方法。因此,本发明提出了尾流立管的涡激升力计算方法。
本发明考虑了上游立管的涡街对尾流立管涡旋泄放的影响,考虑了上下游立管之间的距离,提出了基于下述模型的尾流立管时域涡激升力计算方法:
式中:FL--涡激升力;
CL--升力系数;
CD--拖曳力系数;
ρ--流体密度;
λ--上下游立管轴线之间的距离L与直径之比L/D;
D--尾流立管直径;
U--流速;
--尾流立管顺流向振动速度;
t--时间;
q,b--由试验确定的系数。
公式(1)中的第一项为涡旋泄放引起的升力,其中的(1-λ-q)-1是考虑上游立管涡街对涡激升力幅值的影响,这是本发明的关键,它取决于两立管的间距;第二项为圆柱体振动速度引起的粘性阻力,第三项为圆柱体振动加速度引起的附加质量力。
尾流立管涡旋泄放频率的计算公式中(1-λ-b)-1是考虑上游立管涡街对尾流立管涡旋泄放频率的影响。
由于公式(1)包含尾流立管的顺流向振动速度、横流向振动速度和加速度,因此,必须采用迭代方法计算,如图2所示,具体计算步骤如下:
1)给定尾流立管顺流向振动速度、横流向振动速度和加速度及计算时间的初值(公知技术):
tj=t0=0,
式中:j--时间步数,计算开始时j=0;
i--迭代次数,每个时间步开始时i=0;
2)计算给定流速(流速是设计或分析给定的条件,为已知值)下的约化速度(公知技术):
式中:Vr--约化速度;
U--流体流速,单位m/s;
fn--深水立管的固有频率,单位Hz;
D--深水立管直径,单位m;
在以上两个公式中,
FD--脉动拖曳力;
CD--拖曳力系数;
ρ--流体密度;
λ--上下游立管轴线之间的距离L与直径之比L/D;
D--尾流立管直径;
U--流体流速;
t--时间;
ü--深水立管顺流向振动加速度;
p,b--由试验确定的系数;采用水池或水槽试验,也可采用计算流体力学方法计算得到,相关试验为公知技术;
式中:m-尾流立管的质量;
c-尾流立管的阻尼系数;
k-尾流立管的弯曲刚度;
--第j时间步内第i+1次迭代的顺流向振动加速度;
5)如果ε为预先设定的计算精度(根据需要确定,如ε=1×10-5),则继续进行迭代计算,令:
然后,重复步骤3)~5)的计算;
n为第j时间步内计算脉动拖曳力的总迭代次数;
上式中各参数的含义与公式(1)中对应参数的含义相同。
式中:m-尾流立管的质量;
c-尾流立管的阻尼系数;
k-尾流立管的弯曲刚度;
然后,重复步骤6)~8)的计算;
tj=tj+1=tj+Δt,
m为第j时间步内计算涡激升力的总迭代次数;Δt一般取0.02秒;
重复步骤3)~8)的计算,直至计算时长满足需要(时长根据需要设定,如100秒)。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (2)
2.如权利要求1所述的尾流立管的时域涡激升力确定方法,其特征在于:所述的采用迭代方法进行计算的具体步骤如下:
1)给定尾流立管顺流向振动速度、横流向振动速度和加速度及计算时间的初值:
tj=t0=0,
式中:j--时间步数,计算开始时j=0;
i--迭代次数,每个时间步开始时i=0;
2)计算给定流速下的约化速度:
式中:Vr--约化速度;
U--流体流速;
fn--尾流立管的固有频率;
D--尾流立管直径;
当5≤Vr≤7时,将tj,代入下式计算第j时间步内第i次迭代的脉动拖曳力:
在以上两个公式中,
FD--脉动拖曳力;
CD--拖曳力系数;
ρ--流体密度;
λ--上下游立管轴线之间的距离L与直径之比L/D;
D--尾流立管直径;
U--流体流速;
t--时间;
ü--深水立管顺流向振动加速度;
p,b--由试验确定的系数;
式中:m-尾流立管的质量;
c-尾流立管的阻尼系数;
k-尾流立管的弯曲刚度;
然后,重复步骤3)~5)的计算;
如果 则令:
n为第j时间步内计算脉动拖曳力的总迭代次数;
式中:m-尾流立管的质量;
c-尾流立管的阻尼系数;
k-尾流立管的弯曲刚度;
然后,重复步骤6)~8)的计算;
tj=tj+1=tj+Δt,
m为第j时间步内计算涡激升力的总迭代次数;
重复步骤3)~8)的计算,直至计算时长满足需要。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110291373.5A CN102507084B (zh) | 2011-09-30 | 2011-09-30 | 一种尾流立管的时域涡激升力确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110291373.5A CN102507084B (zh) | 2011-09-30 | 2011-09-30 | 一种尾流立管的时域涡激升力确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102507084A true CN102507084A (zh) | 2012-06-20 |
CN102507084B CN102507084B (zh) | 2014-02-26 |
Family
ID=46219191
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110291373.5A Expired - Fee Related CN102507084B (zh) | 2011-09-30 | 2011-09-30 | 一种尾流立管的时域涡激升力确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102507084B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102507082A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种深水立管的时域涡激升力确定方法 |
CN103512717A (zh) * | 2013-09-26 | 2014-01-15 | 中国石油集团工程设计有限责任公司 | 大型低温蒸发器在两相流作用下的管束振动预测方法 |
CN110110408A (zh) * | 2019-04-25 | 2019-08-09 | 西南石油大学 | 刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法 |
CN110243571A (zh) * | 2019-07-11 | 2019-09-17 | 中国海洋大学 | 一种尾流立管涡激升力频率确定方法 |
CN111985138A (zh) * | 2020-08-21 | 2020-11-24 | 哈尔滨工业大学(威海) | 一种柔性结构横流与顺流方向涡激振动耦合响应预测方法 |
CN112834164A (zh) * | 2020-12-31 | 2021-05-25 | 中国海洋大学 | 一种考虑约化速度和间距的尾流立管涡激升力确定方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59142431A (ja) * | 1983-02-03 | 1984-08-15 | Nishi Nippon Riyuutai Giken:Kk | 流力試験用回流水槽 |
US5303882A (en) * | 1993-02-22 | 1994-04-19 | The United States Of America As Represented By The Secretary Of The Navy | Corner vortex suppressor |
CN101089577A (zh) * | 2007-07-12 | 2007-12-19 | 上海交通大学 | 竖直于拖曳水池中柔性管件模型的涡激振动试验装置 |
CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
-
2011
- 2011-09-30 CN CN201110291373.5A patent/CN102507084B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59142431A (ja) * | 1983-02-03 | 1984-08-15 | Nishi Nippon Riyuutai Giken:Kk | 流力試験用回流水槽 |
US5303882A (en) * | 1993-02-22 | 1994-04-19 | The United States Of America As Represented By The Secretary Of The Navy | Corner vortex suppressor |
CN101089577A (zh) * | 2007-07-12 | 2007-12-19 | 上海交通大学 | 竖直于拖曳水池中柔性管件模型的涡激振动试验装置 |
CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
Non-Patent Citations (1)
Title |
---|
王增波等: "立管扰流装置Helical Strakes的数值模拟分析", 《石油化工高等学校学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102507082A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种深水立管的时域涡激升力确定方法 |
CN102507082B (zh) * | 2011-09-30 | 2014-02-26 | 中国海洋大学 | 一种深水立管的时域涡激升力确定方法 |
CN103512717A (zh) * | 2013-09-26 | 2014-01-15 | 中国石油集团工程设计有限责任公司 | 大型低温蒸发器在两相流作用下的管束振动预测方法 |
CN110110408A (zh) * | 2019-04-25 | 2019-08-09 | 西南石油大学 | 刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法 |
CN110110408B (zh) * | 2019-04-25 | 2022-06-07 | 西南石油大学 | 刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法 |
CN110243571A (zh) * | 2019-07-11 | 2019-09-17 | 中国海洋大学 | 一种尾流立管涡激升力频率确定方法 |
CN111985138A (zh) * | 2020-08-21 | 2020-11-24 | 哈尔滨工业大学(威海) | 一种柔性结构横流与顺流方向涡激振动耦合响应预测方法 |
CN111985138B (zh) * | 2020-08-21 | 2023-12-19 | 哈尔滨工业大学(威海) | 一种柔性结构横流与顺流方向涡激振动耦合响应预测方法 |
CN112834164A (zh) * | 2020-12-31 | 2021-05-25 | 中国海洋大学 | 一种考虑约化速度和间距的尾流立管涡激升力确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102507084B (zh) | 2014-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102507084B (zh) | 一种尾流立管的时域涡激升力确定方法 | |
CN101539477B (zh) | 一种深水顶张式立管涡激振动与疲劳分析的方法 | |
Chaplin et al. | Ringing of a vertical cylinder in waves | |
CN105243275A (zh) | 一种海底多跨管道涡激振动的预报方法 | |
CN105205280A (zh) | 一种单跨管道横流向涡激振动的预报方法 | |
CN102353516B (zh) | 一种深水立管的时域脉动拖曳力确定方法 | |
CN103353382A (zh) | 一种深海顶张力立管参激-涡激联合振动疲劳分析方法 | |
Denys et al. | Fluid structure interaction of piano key weirs | |
Erpicum et al. | Energy dissipation on a stepped spillway downstream of a Piano Key Weir–Experimental study | |
CN102445318B (zh) | 一种顶张式立管顺流向振动分析方法 | |
CN105279373A (zh) | 一种海底自由悬跨管道纯顺流向涡激振动预报方法 | |
CN102507082B (zh) | 一种深水立管的时域涡激升力确定方法 | |
CN203287174U (zh) | 一种深水钻井工况下隔水管振动特性模拟试验装置 | |
Xu et al. | Determination of hydrodynamic coefficients of the floating cage collar with forced oscillation experiments | |
CN102507083B (zh) | 一种尾流立管的脉动拖曳力确定方法 | |
Yang et al. | Investigation on scour scale of piggyback pipeline under wave conditions | |
CN111177926B (zh) | 铺管过程中的涡激振动预报方法 | |
Repalle et al. | CFD simulation of wave run-up on a spar cylinder | |
Van Der Heijden et al. | Fatigue analysis of subsea jumpers due to slug flow | |
CN102507122B (zh) | 一种深水浮筒平台的流致振荡分析方法 | |
CN102509024B (zh) | 一种深水浮筒平台的自激振荡分析方法 | |
Guo et al. | Experimental study on vortex-induced vibration of risers transporting fluid | |
Peng et al. | Research on the Vortex-induced Vibration of A Cylinder Under Oscillatory Flow | |
Mao et al. | The response property of inverted siphon in long distance water distribution system | |
Münch et al. | Hydro elastic behavior of vibrating blades |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140226 Termination date: 20140930 |
|
EXPY | Termination of patent right or utility model |