CN102353516A - 一种深水立管的时域脉动拖曳力确定方法 - Google Patents
一种深水立管的时域脉动拖曳力确定方法 Download PDFInfo
- Publication number
- CN102353516A CN102353516A CN2011102908065A CN201110290806A CN102353516A CN 102353516 A CN102353516 A CN 102353516A CN 2011102908065 A CN2011102908065 A CN 2011102908065A CN 201110290806 A CN201110290806 A CN 201110290806A CN 102353516 A CN102353516 A CN 102353516A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- drag
- deep water
- rho
- water standpipe
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及海洋深水立管的研究方法,具体涉及一种深水立管的时域脉动拖曳力确定方法。该方法将由深水立管的振动速度和加速度产生的粘性阻力和附加质量力引入了脉动拖曳力的计算,并采用迭代方法计算时域的脉动拖曳力,从而建立了一个完整的流固耦合时域脉动拖曳力计算方法。由于本发明增加了深水立管振动速度和加速度产生的粘性阻力和附加质量力,比现有方法更符合圆柱体顺流向涡激振动的受力状态。
Description
技术领域
本发明涉及海洋深水立管的研究方法,具体涉及一种深水立管的时域脉动拖曳力确定方法。
背景技术
脉动拖曳力是圆柱体受到的流体流动方向上的一种交变流体荷载,由于流体的流速是恒定的,因此,脉动拖曳力是圆柱体尾流处的涡旋泄放引起的。对于弹性圆柱体而言,脉动拖曳力将使圆柱体产生沿着流体流动方向的振动。而对于刚性圆柱体,脉动拖曳力并不引起圆柱体振动。这就是说,不论圆柱体是否振动,脉动拖曳力始终是存在的。
脉动拖曳力的交变频率与涡旋泄放的频率有关,而涡旋泄放的频率取决于流体流经圆柱体的速度大小,脉动拖曳力的大小取决于流体流经圆柱体的速度。因此,对于刚性圆柱体来说,脉动拖曳力的大小和频率仅仅与流体的流动速度有关。而对于弹性圆柱体,由于圆柱体本身的振动,不仅使得流体流经圆柱体的速度发生变化,而且使得流体在圆柱体的扰动下产生粘性阻力和附加质量力。因此,对于弹性圆柱体的脉动拖曳力计算不能仅仅考虑流场的流速一个因素,而应考虑圆柱体振动的速度和加速度。
深水立管或海底管线属于弹性圆柱体,其涡激振动为弯曲振动,因此,沿圆柱体轴向的不同位置,其振动速度和加速度是不同的。如果采用相同的脉动拖曳力计算,则误差是可想而知的。
申请人在前一项发明专利“一种深水顶张式立管我即振动与疲劳分析的方法”,ZL200910136583.X中提出了一个涡旋泄放引起的脉动拖曳力模型,而除此之外,传统尾流振子模型并不包括脉动拖曳力的计算,也不涉及顺流向涡激振动。总体来说,现有技术存在如下缺陷:
(1)没有考虑深水立管(弹性圆柱体)顺流向振动的加速度效应,忽略了深水立管周围流场压力梯度产生的附加质量力。因此,计算得到的深水立管振动频率偏低;
(2)没有考虑深水立管顺流向振动的速度效应,忽略了深水立管与流体相对速度变化引起的粘性阻力,因此,计算结果偏小;
(3)不能计算脉动拖曳力沿深水立管轴向的变化,即不能考虑深水立管轴向不同位置的流固耦合引起的脉动拖曳力变化;
(4)没有考虑脉动拖曳力在锁定区和非锁定区的频率变化,仅仅采用锁定区的频率。
发明内容
本发明的目的在于针对现有技术的缺陷,提供一种考虑深水立管顺流向振动加速度引起的附加质量力和振动速度引起的粘性阻力的深水立管的时域脉动拖曳力确定方法,解决弹性圆柱体轴向不同位置的脉动拖曳力计算,实现深水立管顺流向涡激振动的时域非线性分析。
本发明的技术方案如下:一种深水立管的时域脉动拖曳力确定方法,分为非涡旋泄放锁定区和涡旋泄放锁定区,建立的脉动拖曳力时域模型如下:
a)非涡旋泄放锁定区
b)涡旋泄放锁定区
式中,FD--脉动拖曳力;
CD--拖曳力系数;
ρ--流体密度;
D--深水立管直径;
U--流体流速;
--深水立管顺流向振动速度;
t--时间;
ü--深水立管顺流向振动加速度;
针对深水立管顺流向振动速度和加速度,采用迭代方法对以上公式进行计算,分别得到非涡旋泄放锁定区和涡旋泄放锁定区的脉动拖曳力。
进一步,如上所述的深水立管的时域脉动拖曳力确定方法,其中,所述的采用迭代方法进行计算的具体步骤如下:
1)给定深水立管顺流向振动速度和加速度及计算时间的初值:
tj=t0=0,
式中:j--时间步数,计算开始时j=0;
i--迭代次数,每个时间步开始时i=0;
2)计算给定流速下的约化速度:
式中:Vr--约化速度;
U--流体流速;
fn--深水立管的固有频率;
D--深水立管直径;
式中:m--深水立管的质量;
c--深水立管的阻尼系数;
k--深水立管的弯曲刚度;
然后,重复步骤3)~5)的计算;
tj=tj+1=tj+Δt,
n为第j时间步的总迭代次数;
重复步骤3)~5)的计算,直至计算时长满足需要。
本发明的有益效果如下:本发明将由深水立管的振动速度和加速度产生的粘性阻力和附加质量力引入了脉动拖曳力的计算,并采用迭代方法计算时域的脉动拖曳力,从而建立了一个完整的流固耦合时域脉动拖曳力计算方法。与现有的脉动拖曳力计算方法相比,本发明具有下述优点:
(1)增加了深水立管振动速度和加速度产生的粘性阻力和附加质量力,比现有方法更符合圆柱体顺流向涡激振动的受力状态;
(2)可以计算时域的脉动拖曳力,从而可以进行非线性的顺流向涡激振动分析;
(3)可以计算深水立管轴向不同位置的脉动拖曳力时程,因而可以得到脉动拖曳力沿深水立管(弹性圆柱体)轴向的变化规律。
附图说明
图1为本发明的方法流程图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细描述。
传统的涡激振动理论认为,深水立管的顺流向涡激振动非常小,可以忽略不计。因此,深水立管的涡激振动就是垂直流体流动方向的振动被人们所接受。“涡激振动”仅仅意味着深水立管的垂直流向振动,这可以从目前的多款商业软件的计算方法得到证明。Orcaflex、VIVA和Shear7是目前海洋工程界公认的深水立管涡激振动分析软件,被广泛应用于深水立管的涡激振动与疲劳分析。这些软件只提供横向涡激振动分析结果。当然现有的其它分析方法中也没有顺流向涡激振动的分析方法,唯一的计算方法是本发明申请人已获得的前一项发明专利中提出的顺流向涡激振动分析方法,其中包括脉动拖曳力模型。
本发明同时考虑圆柱体顺流向振动速度和加速度效应、并考虑了涡旋泄放锁定区和非锁定区的脉动拖曳力频率变化,提出了基于下述模型的深水立管时域脉动拖曳力计算方法:
a)非涡旋泄放锁定区
b)涡旋泄放锁定区
式中,FD--脉动拖曳力;
CD--拖曳力系数;
ρ--流体密度;
D--深水立管直径;
U--流体流速;
t--时间;
ü--深水立管顺流向振动加速度。
上述公式中,第一项为涡旋泄放引起的脉动拖曳力,第二项为圆柱体振动速度引起的粘性阻力,第三项为圆柱体振动加速度引起的附加质量力。上述两个公式一个用于非涡旋泄放锁定区,另一个用于涡旋泄放锁定区,两个公式的第一项涡旋泄放频率不同,这个现象与现有的涡激振动理论不同,是本发明的主要研究结果。
由于上述两个公式中包含圆柱体的顺流向速度和加速度,因此,必须采用迭代方法计算,如图1所示,具体计算步骤如下:
1)给定深水立管顺流向振动速度和加速度及计算时间的初值(公知技术):
tj=t0=0,
式中:j--时间步数,计算开始时j=0;
i--迭代次数,每个时间步开始时i=0;
2)计算给定流速(流速是设计或分析给定的条件,为已知值)下的约化速度(公知技术):
式中:Vr--约化速度;
U--流体流速;
fn--深水立管的固有频率;
D--深水立管直径;
上述两个公式分别对应非涡旋泄放锁定区和涡旋泄放锁定区的脉动拖曳力计算模型,各参数的含义与上面公式(1)、(2)中对应参数的含义相同。
式中:m--深水立管的质量;
c--深水立管的阻尼系数;
k--深水立管的弯曲刚度;
然后,重复步骤3)~5)的计算;
tj=tj+1=tj+Δt,
n为第j时间步的总迭代次数;Δt一般取0.02秒;
重复步骤3)~5)的计算,直至计算时长满足需要(时长根据需要设定,如100秒)。
本发明相对于现有技术具有如下特点:
1、考虑了深水立管振动的惯性耦合效应
涡激振动是深水立管在流体定常流动作用下的一种振动形式,包括顺流向振动和横流向振动,其中的顺流向振动是由脉动拖曳力引起的。由于深水立管振动加速度对流场产生的扰动作用(流体的加速度为零),深水立管周围的流场将产生压力梯度,该压力梯度对圆柱体的作用被称为附加质量力。因此,深水立管(弹性圆柱体)的脉动拖曳力应在刚性圆柱体脉动拖曳力的基础上增加由圆柱体的振动加速度引起的附加质量力。
2、考虑了深水立管振动的粘性耦合效应
深水立管的顺流向振动速度对流场的扰动作用使深水立管与流体之间的相对速度发生周期性的变化,从而导致深水立管受到交变的流体粘性阻力作用。因此,深水立管(弹性圆柱体)的脉动拖曳力应在刚性圆柱体脉动拖曳力的基础上增加由圆柱体与流体之间周期性变化的相对速度引起的交变粘性阻力。
3、解决了深水立管弯曲振动的展向流固耦合效应
对于弹性圆柱体,如深水立管和海底管线,其涡激振动为弯曲振动,因此,沿圆柱体轴向的不同位置,其振动速度和加速度是不同的。如果采用相同的脉动拖曳力计算,则误差是可想而知的。而本发明提出的脉动拖曳力计算方法可以解决上述问题,使计算结果更符合实际自然现象。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。
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次迭代的脉动拖曳力:
式中:m--深水立管的质量;
c--深水立管的阻尼系数;
k--深水立管的弯曲刚度;
然后,重复步骤3)~5)的计算;
tj=tj+1=tj+Δt,
n为第j时间步的总迭代次数;
重复步骤(3~5)的计算,直至计算时长满足需要。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110290806 CN102353516B (zh) | 2011-09-29 | 2011-09-29 | 一种深水立管的时域脉动拖曳力确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110290806 CN102353516B (zh) | 2011-09-29 | 2011-09-29 | 一种深水立管的时域脉动拖曳力确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102353516A true CN102353516A (zh) | 2012-02-15 |
CN102353516B CN102353516B (zh) | 2013-12-25 |
Family
ID=45577127
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110290806 Expired - Fee Related CN102353516B (zh) | 2011-09-29 | 2011-09-29 | 一种深水立管的时域脉动拖曳力确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102353516B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109163878A (zh) * | 2018-07-17 | 2019-01-08 | 中山大学 | 植被拖拽力系数测量方法、装置、计算机设备和存储介质 |
CN109932368A (zh) * | 2019-04-02 | 2019-06-25 | 北京林业大学 | 一种柔性植被群拖曳力系数的测量方法 |
CN110598337A (zh) * | 2019-09-17 | 2019-12-20 | 中国海洋大学 | 一种圆柱体涡激振动的流固耦合时域分析方法 |
CN111175018A (zh) * | 2020-01-09 | 2020-05-19 | 哈尔滨电机厂有限责任公司 | 一种判定高比速水力机械模型肘管下游侧刚度不足的方法 |
CN115859748A (zh) * | 2023-02-15 | 2023-03-28 | 山东科技大学 | 一种拖曳式温盐深测量仪柔性缆绳涡激振动分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1268498A (en) * | 1968-07-05 | 1972-03-29 | Continental Elektroindustrie Aktiengesellschaft | Apparatus for detecting a fluid pressure differential |
DE2317811A1 (de) * | 1973-04-09 | 1974-10-17 | Werner O Dipl Ing Sommer | Vorrichtung zur messung von druecken und kraeften, die von stroemenden medien verursacht werden |
CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
-
2011
- 2011-09-29 CN CN 201110290806 patent/CN102353516B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1268498A (en) * | 1968-07-05 | 1972-03-29 | Continental Elektroindustrie Aktiengesellschaft | Apparatus for detecting a fluid pressure differential |
DE2317811A1 (de) * | 1973-04-09 | 1974-10-17 | Werner O Dipl Ing Sommer | Vorrichtung zur messung von druecken und kraeften, die von stroemenden medien verursacht werden |
CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
Non-Patent Citations (1)
Title |
---|
唐世振: "考虑顺流向振动的深水顶张力立管涡激振动分析", 《中国海洋大学博士论文》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109163878A (zh) * | 2018-07-17 | 2019-01-08 | 中山大学 | 植被拖拽力系数测量方法、装置、计算机设备和存储介质 |
CN109163878B (zh) * | 2018-07-17 | 2020-01-17 | 中山大学 | 植被拖拽力系数测量方法、装置、计算机设备和存储介质 |
CN109932368A (zh) * | 2019-04-02 | 2019-06-25 | 北京林业大学 | 一种柔性植被群拖曳力系数的测量方法 |
CN109932368B (zh) * | 2019-04-02 | 2021-06-18 | 北京林业大学 | 一种柔性植被群拖曳力系数的测量方法 |
CN110598337A (zh) * | 2019-09-17 | 2019-12-20 | 中国海洋大学 | 一种圆柱体涡激振动的流固耦合时域分析方法 |
CN111175018A (zh) * | 2020-01-09 | 2020-05-19 | 哈尔滨电机厂有限责任公司 | 一种判定高比速水力机械模型肘管下游侧刚度不足的方法 |
CN111175018B (zh) * | 2020-01-09 | 2021-08-20 | 哈尔滨电机厂有限责任公司 | 一种判定高比速水力机械模型肘管下游侧刚度不足的方法 |
CN115859748A (zh) * | 2023-02-15 | 2023-03-28 | 山东科技大学 | 一种拖曳式温盐深测量仪柔性缆绳涡激振动分析方法 |
CN115859748B (zh) * | 2023-02-15 | 2023-05-05 | 山东科技大学 | 一种拖曳式温盐深测量仪柔性缆绳涡激振动分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102353516B (zh) | 2013-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101539477B (zh) | 一种深水顶张式立管涡激振动与疲劳分析的方法 | |
CN102353516B (zh) | 一种深水立管的时域脉动拖曳力确定方法 | |
Lian et al. | Analysis on flow induced motion of cylinders with different cross sections and the potential capacity of energy transference from the flow | |
CN109827734A (zh) | 一种评估内外流作用下深海立管涡激振动的方法 | |
CN102507084B (zh) | 一种尾流立管的时域涡激升力确定方法 | |
CN103353382A (zh) | 一种深海顶张力立管参激-涡激联合振动疲劳分析方法 | |
CN102507082B (zh) | 一种深水立管的时域涡激升力确定方法 | |
CN105205280A (zh) | 一种单跨管道横流向涡激振动的预报方法 | |
CN102445318B (zh) | 一种顶张式立管顺流向振动分析方法 | |
CN105279373A (zh) | 一种海底自由悬跨管道纯顺流向涡激振动预报方法 | |
CN106934138A (zh) | 内输气液两相段塞流的海底管道悬跨段振动特性分析方法 | |
CN102222168A (zh) | 一种深水钻井立管参激横向振动分析方法 | |
Thorsen | Time domain analysis of vortex-induced vibrations | |
CN111177926B (zh) | 铺管过程中的涡激振动预报方法 | |
Gupta et al. | Identification of vortex-induced-response parameters in time domain | |
CN102507083B (zh) | 一种尾流立管的脉动拖曳力确定方法 | |
Li et al. | Two-phase flow induced vibration of subsea span pipeline | |
Van Der Heijden et al. | Fatigue analysis of subsea jumpers due to slug flow | |
Major et al. | Effect of drill string rotation on the dynamic response of drilling risers | |
Lefebure et al. | Analytical and CFD study of the influence of control parameters on the maximum efficiency of a hydro-power conversion system based on vortex-induced vibrations | |
Yin et al. | Experimental and numerical analysis of forced motion of a circular cylinder | |
Slingsby | Dynamic interaction of subsea pipeline spans due to vortex-induced vibrations | |
Lei et al. | Numerical analysis of the effects of parametric excitation on riser VIV | |
CN102507122A (zh) | 一种深水浮筒平台的流致振荡分析方法 | |
Thorsen et al. | Time domain simulation of vortex-induced vibrations based on phase-coupled oscillator synchronization |
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: 20131225 Termination date: 20140929 |
|
EXPY | Termination of patent right or utility model |