CN102353506A - 一种深水顶张式立管竖向振动分析方法 - Google Patents
一种深水顶张式立管竖向振动分析方法 Download PDFInfo
- Publication number
- CN102353506A CN102353506A CN2011101610896A CN201110161089A CN102353506A CN 102353506 A CN102353506 A CN 102353506A CN 2011101610896 A CN2011101610896 A CN 2011101610896A CN 201110161089 A CN201110161089 A CN 201110161089A CN 102353506 A CN102353506 A CN 102353506A
- Authority
- CN
- China
- Prior art keywords
- standpipe
- partiald
- displacement
- centerdot
- axial
- 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
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明涉及海洋深水立管的研究方法,具体涉及一种深水顶张式立管竖向振动分析方法。该方法考虑了横向弯曲振动大位移引起的轴向位移对立管轴向振动的影响,将横向弯曲挠度引起的轴向位移转化为附加张力引入轴向刚度;同时,考虑了剪切变形对立管轴向位移的影响,将弯曲刚度与剪切刚度的比值引入轴向振动数学模型,使深水顶张式立管的轴向振动分析更加符合立管的实际受力和变形状态。
Description
技术领域
本发明涉及海洋深水立管的研究方法,具体涉及一种深水顶张式立管竖向振动分析方法。
背景技术
深水顶张式立管是深水海洋立管的一种主要类型,立管的下端与海底井口连接,管壁张力直接作用在井口上,而管内的流体源源不断地从海底穿过立管流向浮式平台。竖向振动是浮式平台垂荡运动引起的深水顶张式立管轴向振动,由于现有分析方法采用的是仅受张力而不受横向力作用的等截面直杆的轴向振动理论,其数学模型中不包括显式的张力,仅包括轴向位移,因此,该模型不包括横向振动的影响,且立管的轴向位移响应不是通过轴向荷载的强迫振动计算得到的,而是通过端部的轴向位移计算求解的,故而,也被称为参数(位移本身是对荷载响应,因此,称为参数)振动。但是,深水顶张式立管在海洋环境中不仅受张力作用,同时还受到海洋环境荷载的作用,因此,轴向振动和横向弯曲振动是同时发生的。在小变形条件下,横向弯曲引起的轴向位移较小,可以忽略。但是,对于深水立管的大位移横向弯曲振动,由于横向弯曲引起的轴向位移已经超出了小变形的范围,轴向振动分析时应考虑横向弯曲引起的轴向位移影响。因此,现有的轴向振动(参数振动)分析方法与深水顶张式立管的实际受力(同时受张力和横向荷载作用)和变形(轴向变形和横向变形耦合)状态有较大的区别。
现有的深水顶张式立管轴向振动分析方法是基于等截面直杆的轴向振动理论,该理论建立在小变形假设和没有横向弯曲变形的基础上。基于该理论的深水顶张式立管轴向振动方程可表示为:
式中:u为立管轴向位移;
x为立管轴向坐标;
t为时间;
为立管单位长度质量;
c为立管材料的阻尼系数;
EA为立管截面拉压刚度。
公式(1)是没有横向弯曲变形且符合小变形假设的深水顶张式立管轴向振动方程,没有考虑横向弯曲大位移对立管轴向振动的影响,即横向弯曲大位移引起的轴向位移。
发明内容
本发明的目的在于针对现有技术的缺陷,提供一种深水顶张式立管竖向振动分析方法,使深水顶张式立管的竖向振动分析更加符合立管的实际受力和变形状态。
本发明的技术方案如下:一种深水顶张式立管竖向振动分析方法,该方法采用考虑横向弯曲振动大位移引起轴向位移的深水顶张式立管轴向振动分析模型,其模型方程如下:
式中:u为立管轴向位移;
x为立管轴向坐标;
t为时间;
c为立管材料的阻尼系数;
EA为立管截面拉压刚度;
k为剪切不均匀系数;
T为立管张力;
y为立管横向位移;
EI为立管截面拉弯刚度;
利用上述方程,计算立管轴向振动的加速度、速度、位移、应力和应变随时间的变化。
进一步,如上所述的深水顶张式立管竖向振动分析方法,该方法计算立管轴向振动的加速度、速度、位移、应力和应变随时间变化的具体过程如下:
(1)将立管划分为若干个单元;
(2)将单元的位移函数表示为插值函数的形式:
u=[S]{q},
y=[N]{a};
其中,[S]为轴向位移插值函数矩阵,{q}为单元节点轴向位移向量;[N]为横向位移插值函数矩阵,{a}为单元节点横向位移向量;
(3)采用伽辽金方法,将深水顶张式立管轴向振动模型的方程转换为矩阵方程如下:
式中:
阻尼矩阵[C]=α[M]+β[K],α,β为瑞雷阻尼系数;
为竖向速度矢量;
{q}为竖向位移矢量;
{a}为横向位移矢量;
l为单元长度;
n为单元数量;
{P}为横向弯曲位移引起的轴向扰动;
(4)基于立管的横向弯曲振动方程求出每个时间增量的横向位移{a},并将其代入步骤(3)中{P}的表达式,计算出{P};
(5)由步骤(3)中质量矩阵[M]、刚度矩阵[K]和阻尼矩阵[C]的表达式分别计算出立管的质量矩阵、刚度矩阵和阻尼矩阵;
此处,{ΔPi}为与时间增量Δt对应的横向弯曲位移增量引起的荷载增量;
(7)由加速度增量、速度增量和位移增量按下式计算第一个时间增量后的加速度、速度和位移,此时i=0:
qi+1=qi+Δqi
(8)根据步骤(7)中得到的位移计算出立管的应力和张力;
(9)重复步骤(4)~(8),直至时间t达到要求的时长,即可计算出立管轴向振动的加速度、速度、位移、应力和轴向应变随时间的变化。
更进一步,如上所述的深水顶张式立管竖向振动分析方法,步骤(4)中所述的立管的横向弯曲振动方程如下:
式中:y为立管横向弯曲位移;
x为立管的轴向坐标;
t为时间;
EI为立管横截面抗弯刚度;
T为立管张力,是时间和立管轴向坐标的函数,即:T=T(x,t);
c为阻尼系数;
k为剪切不均匀系数;
εT为张力引起的立管轴向应变;
q(x,t)为作用在立管上的流体荷载。
本发明的有益效果如下:本发明在深水顶张式立管的轴向振动分析中,考虑了横向弯曲振动大位移引起的轴向位移对立管轴向振动的影响,将横向弯曲挠度引起的轴向位移转化为附加张力引入轴向刚度;同时,考虑了剪切变形对立管轴向位移的影响,将弯曲刚度与剪切刚度的比值引入轴向振动数学模型,使深水顶张式立管的轴向振动分析更加符合立管的实际受力和变形状态。
附图说明
图1为本发明的深水顶张式立管的竖向振动分析方法流程图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细描述。
本发明所提供的深水顶张式立管竖向振动分析方法采用考虑横向弯曲振动大位移引起轴向位移的深水顶张式立管轴向振动分析模型,其模型方程如下:
式中:u为立管轴向位移;
x为立管轴向坐标;
t为时间;
c为立管材料的阻尼系数;
EA为立管截面拉压刚度;
k为剪切不均匀系数;
T为立管张力;
y为立管横向位移;
EI为立管截面拉弯刚度;
利用上述方程,计算立管轴向振动的加速度、速度、位移、应力和应变随时间的变化。
深水顶张式立管竖向振动分析方法的具体步骤如下:
(S1)将立管划分为若干个单元;
(S2)将单元的位移函数表示为插值函数的形式:
u=[S]{q} (3)
y=[N][a} (4)
公式(3)、(4)中,[S]为轴向位移插值函数矩阵,{q}为单元节点轴向位移向量;[N]为横向位移插值函数矩阵,{a}为单元节点横向位移向量;
(S3)采用伽辽金方法(此为本领域的公知常识),
将深水顶张式立管轴向振动分析模型的方程(公式(2))转换为矩阵方程如下:
式中:
阻尼矩阵[C]=α[M]+β[K],α,β为瑞雷阻尼系数 (9)
{q}为竖向位移矢量;
{a}为横向位移矢量;
l为单元长度;
n为单元数量;
{P}为横向弯曲位移引起的轴向扰动,相当于轴向强迫振动的外力;
(S4)基于立管的横向弯曲振动方程求出每个时间增量的横向位移{a},并将其代入步骤(S3)中{P}的表达式(公式(10)),计算出{P};
此处,横向位移{a}是与每个时间增量的竖向位移对应的横向位移,即每个时间增量先计算横向位移,然后再依据下述各步骤计算竖向位移;
(S5)由步骤(S3)中质量矩阵[M]、刚度矩阵[K]和阻尼矩阵[C]的表达式(公式(7、8、9))分别计算出立管的质量矩阵、刚度矩阵和阻尼矩阵;
公式(11)中,{ΔPi}为与时间增量Δt对应的横向弯曲位移增量引起的荷载增量;
(S7)由加速度增量、速度增量和位移增量按下式计算第一个时间增量后的加速度、速度和位移,此时i=0:
qi+1=qi+Δqi
(S8)根据步骤(S7)中得到的位移计算出立管的应力和张力;
(S9)重复步骤(S4)~(S8),直至时间t达到要求的时长,即可计算出立管轴向振动的加速度、速度、位移、应力和轴向应变随时间的变化。
上述方法中步骤(S4)所述的立管的横向弯曲振动方程如下:
式中:y为立管横向弯曲位移;
x为立管的轴向坐标;
t为时间;
EI为立管横截面抗弯刚度;
T为立管张力,是时间和立管轴向坐标的函数,即:T=T(x,t);
c为阻尼系数;
k为剪切不均匀系数;
εT为张力引起的立管轴向应变;
q(x,t)为作用在立管上的流体荷载。
采用伽辽金方法(此为本领域的公知常识),
上式中,L是立管的总长度,δy是弯曲位移的变分,在公式(13)中的作用是权函数,此为公知常识;
将横向弯曲振动模型的方程(公式(12))转换为矩阵方程如下:
式中:
阻尼矩阵[C]=α[M]+β[K],α,β为瑞雷阻尼系数; (17)
为速度矢量;
{a}为位移矢量;
l为单元长度;
n为单元数量;
将立管的质量矩阵、阻尼矩阵、刚度矩阵和荷载向量代入公式(14),采用逐步积分法按下式求出立管第一个时间增量Δt位移增量{Δai},此时i=0:
此处,{ΔFi}是与时间增量Δt对应的荷载增量;
由位移增量按下式计算第一个时间增量后的位移,此时i=0:
ai+1=ai+Δai
依据上述方法,计算出时间t达到要求的时长内每个时间增量的立管横向位移。
上述立管横向位移具体的计算方法可以参照申请人同期进行的专利申请“一种深水顶张式立管弯曲振动分析方法”。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (3)
2.如权利要求1所述的深水顶张式立管竖向振动分析方法,其特征在于:该方法计算立管轴向振动的加速度、速度、位移、应力和应变随时间变化的具体过程如下:
(1)将立管划分为若干个单元;
(2)将单元的位移函数表示为插值函数的形式:
u=[S]{q},
y=[N]{a};
其中,[S]为轴向位移插值函数矩阵,{q}为单元节点轴向位移向量;[N]为横向位移插值函数矩阵,{a}为单元节点横向位移向量;
(3)采用伽辽金方法,将深水顶张式立管轴向振动分析模型的方程转换为矩阵方程如下:
式中:
阻尼矩阵[C]=α[M]+β[K],α,β为瑞雷阻尼系数;
{q}为竖向位移矢量;
{a}为横向位移矢量;
l为单元长度;
n为单元数量;
{P}为横向弯曲位移引起的轴向扰动;
(4)基于立管的横向弯曲振动方程求出每个时间增量的横向位移{a},并将其代入步骤(3)中{P}的表达式,计算出{P};
(5)由步骤(3)中质量矩阵[M]、刚度矩阵[K]和阻尼矩阵[C]的表达式分别计算出立管的质量矩阵、刚度矩阵和阻尼矩阵;
此处,{ΔPi}为与时间增量Δt对应的横向弯曲位移增量引起的荷载增量;
(7)由加速度增量、速度增量和位移增量按下式计算第一个时间增量后的加速度、速度和位移,此时i=0:
qi+1=qi+Δqi
(8)根据步骤(7)中得到的位移计算出立管的应力和张力;
(9)重复步骤(4)~(8),直至时间t达到要求的时长,即可计算出立管轴向振动的加速度、速度、位移、应力和轴向应变随时间的变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110161089.6A CN102353506B (zh) | 2011-06-16 | 2011-06-16 | 一种深水顶张式立管竖向振动分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110161089.6A CN102353506B (zh) | 2011-06-16 | 2011-06-16 | 一种深水顶张式立管竖向振动分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102353506A true CN102353506A (zh) | 2012-02-15 |
CN102353506B CN102353506B (zh) | 2014-07-02 |
Family
ID=45577117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110161089.6A Expired - Fee Related CN102353506B (zh) | 2011-06-16 | 2011-06-16 | 一种深水顶张式立管竖向振动分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102353506B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102607787A (zh) * | 2012-03-23 | 2012-07-25 | 中国海洋大学 | 一种测试内流对海洋立管动力特性影响的试验方法 |
CN108491615A (zh) * | 2018-03-17 | 2018-09-04 | 中国海洋大学 | 一种三层顶张式立管动力响应的有限元分析方法 |
CN109388907A (zh) * | 2018-10-31 | 2019-02-26 | 中船动力研究院有限公司 | 一种具有预设纵向振动动柔度的轴系的设计方法 |
CN110378060A (zh) * | 2019-07-26 | 2019-10-25 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101539477B (zh) * | 2009-05-08 | 2010-11-03 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
-
2011
- 2011-06-16 CN CN201110161089.6A patent/CN102353506B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101539477B (zh) * | 2009-05-08 | 2010-11-03 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
Non-Patent Citations (2)
Title |
---|
N. SRINIL ET AL.: "Reduced-order modelling of vortex-induced vibration of catenary riser", 《OCEAN ENGINEERING》 * |
唐世振等: "深水立管两向自由度涡激振动的数值分析", 《振动与冲击》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102607787A (zh) * | 2012-03-23 | 2012-07-25 | 中国海洋大学 | 一种测试内流对海洋立管动力特性影响的试验方法 |
CN102607787B (zh) * | 2012-03-23 | 2017-05-17 | 中国海洋大学 | 一种测试内流对海洋立管动力特性影响的试验方法 |
CN108491615A (zh) * | 2018-03-17 | 2018-09-04 | 中国海洋大学 | 一种三层顶张式立管动力响应的有限元分析方法 |
CN108491615B (zh) * | 2018-03-17 | 2022-09-02 | 中国海洋大学 | 一种三层顶张式立管动力响应的有限元分析方法 |
CN109388907A (zh) * | 2018-10-31 | 2019-02-26 | 中船动力研究院有限公司 | 一种具有预设纵向振动动柔度的轴系的设计方法 |
CN109388907B (zh) * | 2018-10-31 | 2022-10-14 | 中船动力研究院有限公司 | 一种具有预设纵向振动动柔度的轴系的设计方法 |
CN110378060A (zh) * | 2019-07-26 | 2019-10-25 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
CN110378060B (zh) * | 2019-07-26 | 2021-02-09 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102353506B (zh) | 2014-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Static analysis of deep-water marine riser subjected to both axial and lateral forces in its installation | |
CN101539477B (zh) | 一种深水顶张式立管涡激振动与疲劳分析的方法 | |
Wang et al. | A nonlinear model for deepwater steel lazy-wave riser configuration with ocean current and internal flow | |
CN102353506A (zh) | 一种深水顶张式立管竖向振动分析方法 | |
CN102222168A (zh) | 一种深水钻井立管参激横向振动分析方法 | |
CN102252897B (zh) | 一种深水顶张式立管弯曲振动分析方法 | |
CN102787677B (zh) | 工程结构抗震烈度预测方法 | |
Esfahanian et al. | Fluid-Structure Interaction in microchannel using Lattice Boltzmann method and size-dependent beam element | |
CN101713707A (zh) | 深海钢悬链线立管波浪响应分析方法 | |
Ng et al. | Fluid/structure interaction study on the variation of radial gate’s gap height in dam | |
CN103902754A (zh) | 一种深海钢悬链线立管触地点动力响应分析方法 | |
Liu et al. | Nonlinear vibration model and response characteristic of drilling risers in deep-sea under soft suspension evacuation condition | |
CN102445318A (zh) | 一种顶张式立管顺流向振动分析方法 | |
CN102353516A (zh) | 一种深水立管的时域脉动拖曳力确定方法 | |
CN102507082B (zh) | 一种深水立管的时域涡激升力确定方法 | |
CN102289536B (zh) | 一种深水钻井立管参数横力弯曲分析方法 | |
CN116911026A (zh) | 一种岩石渗流-应力耦合损伤本构模型的构建方法 | |
CN105184102A (zh) | 一种管土动态耦合作用分析方法 | |
CN102507115B (zh) | 一种深水顶张式立管大变形弯曲振动的分析方法 | |
Pavlou | Flow-riser interaction in deep-sea mining: an analytic approach for multi-layered FRP risers | |
CN105808505A (zh) | 一种计算弧形地连墙侧位移的方法 | |
Li et al. | Hydroelastic simulation of breaking wave impact on a flexible coastal seawall | |
Muthukkumaran et al. | Three dimensional analysis of piles on sloping ground subjected to passive load induced by surcharge | |
CN106294927B (zh) | 一种载流管道动力特性的计算方法 | |
CN102507083A (zh) | 一种尾流立管的脉动拖曳力确定方法 |
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: 20140702 Termination date: 20150616 |
|
EXPY | Termination of patent right or utility model |