CN102252897A - 一种深水顶张式立管弯曲振动分析方法 - Google Patents
一种深水顶张式立管弯曲振动分析方法 Download PDFInfo
- Publication number
- CN102252897A CN102252897A CN2011101612213A CN201110161221A CN102252897A CN 102252897 A CN102252897 A CN 102252897A CN 2011101612213 A CN2011101612213 A CN 2011101612213A CN 201110161221 A CN201110161221 A CN 201110161221A CN 102252897 A CN102252897 A CN 102252897A
- Authority
- CN
- China
- Prior art keywords
- standpipe
- centerdot
- partiald
- matrix
- displacement
- 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)
Abstract
本发明涉及海洋深水立管的研究方法,具体涉及一种深水顶张式立管弯曲振动分析方法。该方法考虑了大位移引起的截面转动对立管弯曲刚度的影响,将截面转动引起的张力作用方向变化转化为立管弯曲刚度的一部分,比现有方法的几何刚度更完整,且更符合深水立管的实际受力和变形状态;同时,该方法还考虑了剪切变形对立管弯曲刚度的影响,将张力与剪切刚度的比值引入弯曲刚度,修正了现有方法的不合理假定,使深水顶张式立管的弯曲振动分析更加符合实际。
Description
技术领域
本发明涉及海洋深水立管的研究方法,具体涉及一种深水顶张式立管弯曲振动分析方法。
背景技术
深水顶张式立管是深水海洋立管的一种主要类型,立管的下端与海底井口连接,管壁张力直接作用在井口上,而管内的流体源源不断地从海底穿过立管流向浮式平台。弯曲振动是深水顶张式立管在海洋环境荷载作用下的主要振动形式,也是引起立管疲劳损伤的主要因素。因此,弯曲振动响应分析是深水顶张式立管设计的重要内容,也是疲劳分析的基础。目前,深水立管的弯曲振动分析方法是基于Eular-Bernoulli梁的纯弯曲理论,该理论建立在小变形假设的基础上,即梁的挠度远远小于其横截面尺寸。
基于该理论的深水顶张式立管弯曲振动方程可表示为:
式中:y为立管横向弯曲位移;
x为立管的轴向坐标;
t为时间;
EI为立管横截面抗弯刚度;
T为立管张力,是时间和立管轴向坐标的函数,即:T=T(x,t);
c为阻尼系数;
q(x,t)为作用在立管上的流体荷载。
式(1)是小变形条件下的深水顶张式立管弯曲振动方程,没有考虑大位移引起的立管截面转动对横向弯曲振动的影响,即张力作用方向改变引起的弹性回复力,而仅仅考虑张力作为几何刚度对弯曲刚度的贡献:
式中:[KG]——立管的几何刚度矩阵;
[N]——立管单元的插值函数矩阵。
现有的深水顶张式立管弯曲振动分析方法的主要缺陷如下:
1.没有考虑深水顶张式立管弯曲振动大位移引起的横截面转动对弯曲刚度的影响。
现有技术是建立在纯弯曲理论和小变形假设基础上的,因此,仅考虑张力的几何刚度部分,而没有考虑由于截面转动引起的张力作用方向变化对立管弯曲刚度的影响。
2.没有考虑剪切变形对深水立管弯曲振动的影响。
现有技术基于纯弯曲理论,因此,没有考虑深水顶张式立管剪切变形。如果仅仅从深水顶张式立管的外观几何形状来看,它属于细长梁,其长细比较大,可以采用纯弯曲模型。但立管是管状结构,其截面的剪切刚度较小,因此,其剪切变形的影响较大,特别是张力与截面剪切刚度的比值较大,从而对弯曲刚度影响较大。
发明内容
本发明的目的在于针对现有技术的缺陷,提供一种深水顶张式立管弯曲振动分析方法,使深水顶张式立管的弯曲振动分析更加符合实际。
本发明的技术方案如下:一种深水顶张式立管弯曲振动分析方法,该方法同时考虑大位移引起的截面转动和剪切变形,提供的深水顶张式立管弯曲振动分析模型的方程如下:
式中:y为立管横向弯曲位移;
x为立管的轴向坐标;
t为时间;
EI为立管横截面抗弯刚度;
T为立管张力,是时间和立管轴向坐标的函数,即:T=T(x,t);
c为阻尼系数;
k为剪切不均匀系数;
εT为张力引起的立管轴向应变;
q(x,t)为作用在立管上的流体荷载;
利用上述方程,计算立管弯曲振动的加速度、速度、位移、应力和应变随时间的变化。
进一步,如上所述的深水顶张式立管弯曲振动分析方法,该方法计算立管弯曲振动的加速度、速度、位移、应力和应变随时间变化的具体过程如下:
(1)将立管划分为若干个单元;
(2)将单元的位移函数表示为插值函数的形式:y=[N]{a},
其中,[N]是单元插值函数矩阵,{a}是单元节点位移向量;
(3)采用伽辽金方法,将深水顶张式立管弯曲振动分析模型的方程转换为矩阵方程如下:
式中:
阻尼矩阵[C]=α[M]+β[K],α,β为瑞雷阻尼系数;
为加速度矢量;
{a}为位移矢量;
l为单元长度;
n为单元数量;
(4)基于立管的初始张力求出立管的初始轴向应变εT,并将其带入步骤(3)中刚度矩阵[K]的表达式,计算出立管的初始刚度矩阵;
(5)由步骤(3)中质量矩阵[M]、阻尼矩阵[C]和荷载向量{F}的表达式分别计算出立管的质量矩阵、阻尼矩阵和荷载向量;
{ΔFi}是与时间增量Δt对应的荷载增量;
(7)由加速度增量、速度增量和位移增量按下式计算第一个时间增量后的加速度、速度和位移,此时i=0:
ai+1=ai+Δai
(8)根据步骤(7)中得到的位移计算出立管的应力和轴向应变;
(9)将步骤(8)中得到的轴向应变代入步骤(3)中刚度矩阵[K]的表达式,计算出第一个时间增量后的立管刚度矩阵;
(10)重复步骤(5)~(9),直至时间t达到要求的时长,即可计算出立管弯曲振动的加速度、速度、位移、应力和轴向应变随时间的变化。
本发明的有益效果如下:本发明在深水顶张式立管的弯曲振动分析中,考虑了大位移引起的截面转动对立管弯曲刚度的影响,将截面转动引起的张力作用方向变化转化为立管弯曲刚度的一部分,比现有方法的几何刚度更完整,且更符合深水顶张式立管的实际受力和变形状态;同时,本发明还考虑了剪切变形对立管弯曲刚度的影响,将张力与剪切刚度的比值引入弯曲刚度,修正了现有方法的不合理假定,使深水顶张式立管的弯曲振动分析更加符合实际。
附图说明
图1为本发明的深水顶张式立管的弯曲振动分析方法流程图;
图2为立管小变形的微元体模型示意图;
图3为立管大变形的微元体模型示意图;
图4为纯弯曲梁的微元体模型示意图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细描述。
本发明采用同时考虑大位移引起的截面转动和剪切变形的深水顶张式立管弯曲振动分析模型,方程如下:
式中:y为立管横向弯曲位移;
x为立管的轴向坐标;
t为时间;
EI为立管横截面抗弯刚度;
T为立管张力,是时间和立管轴向坐标的函数,即:T=T(x,t);
c为阻尼系数;
k为剪切不均匀系数;
εT为张力引起的立管轴向应变;
q(x,t)为作用在立管上的流体荷载。
上式中,2.6kεT即为张力和剪切变形对立管弯曲刚度的影响。
深水顶张式立管弯曲振动分析方法的具体步骤如下:
(S1)将立管划分为若干个单元;
(S2)将单元的位移函数表示为插值函数的形式:
y=[N]{a} (4)
其中,[N]是单元插值函数矩阵,{a}是单元节点位移向量;
(S3)采用伽辽金方法(此为本领域的公知常识),
上式中,L是立管的总长度,δy是弯曲位移的变分,在公式(5)中的作用是权函数,此为公知常识;
将深水顶张式立管弯曲振动分析模型的方程(公式(3))转换为矩阵方程如下:
式中:
阻尼矩阵[C]=α[M]+β[K],α,β为瑞雷阻尼系数; (9)
荷载向量
{a}为位移矢量;
l为单元长度;
n为单元数量;
(S4)基于立管的初始张力求出立管的初始轴向应变εT,并将其带入步骤(S3)中刚度矩阵[K]的表达式(公式(8)),计算出立管的初始刚度矩阵;
(S5)由步骤(S3)中质量矩阵[M]、阻尼矩阵[C]和荷载向量{F}的表达式(公式(7)、(9)、(10))分别计算出立管的质量矩阵、阻尼矩阵和荷载向量;
(S6)将立管的质量矩阵、阻尼矩阵、刚度矩阵和荷载向量代入步骤(S3)中的矩阵方程(公式(6)),采用逐步积分法按下式求出立管第一个时间增量Δt(此时t=Δt)后的加速度增量速度增量和位移增量{Δai},此时i=0:
此处,{ΔFi}是与时间增量Δt对应的荷载增量;
(S7)由加速度增量、速度增量和位移增量按下式计算第一个时间增量后的加速度、速度和位移,此时i=0:
ai+1=ai+Δai
(S8)根据步骤(S7)中得到的位移计算出立管的应力和轴向应变;
(S9)将步骤(S8)中得到的轴向应变代入步骤(S3)中刚度矩阵[K]的表达式(公式(8)),计算出第一个时间增量后的立管刚度矩阵;
(S10)重复步骤(S5)~(S9),直至时间t达到要求的时长,即可计算出立管弯曲振动的加速度、速度、位移、应力和轴向应变随时间的变化(即时程)。
本发明考虑了深水顶张式立管弯曲振动的横向大位移引起的横截面转动对立管弯曲刚度的影响,比现有方法的几何刚度更完整,且更符合深水顶张式立管的实际受力和变形状态。
现有方法建立在小变形假设的基础上,即假定立管变形后的微元段两截面仍保持平行,如图2所示。但深水立管的横向弯曲振动属于大变形(大位移),其变形后的微元段两横截面有相对转动,如图3所示。因此,变形后的张力不再保持水平方向,其横向分量将提供部分回复力,从而增大了弯曲刚度。本发明将张力对回复力的贡献从仅仅考虑微元体两截面因剪切变形发生相对错动(见图2)的几何刚度发展成同时考虑截面相对错动和相对转动(见图3)耦合的几何刚度。
另外,本发明考虑了深水立管剪切变形对弯曲刚度的影响,修正了现有方法的不合理假定。
现有方法采用图2微元体模型,这与纯弯曲梁(不考虑剪切变形)不同,两横截面之间有上下的错动。而纯弯曲梁的微元体模型如图4所示。
如果按照图4的纯弯曲梁模型,则张力对弯曲刚度没有影响,即几何刚度等于零[KG]=0。现有技术为了能够考虑张力梁的几何刚度而采用了图2的模型,这意味着考虑了剪切变形,但却采用了纯弯曲的数学模型(此为公知技术):
式中:θ为弯曲引起的梁截面转角,即挠度的导数等于转角。这表明,现有技术的数学模型和立学模型不一致。本发明克服了现有技术的这一不足,采用了考虑剪切变形的剪弯梁数学模型:
式中:γ为梁截面的剪应变,即挠度的导数等于截面转角与剪应变之和。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (2)
1.一种深水顶张式立管弯曲振动分析方法,其特征在于:该方法同时考虑大位移引起的截面转动和剪切变形,提供的深水顶张式立管弯曲振动分析模型的方程如下:
式中:y为立管横向弯曲位移;
x为立管的轴向坐标;
t为时间;
EI为立管横截面抗弯刚度;
T为立管张力,是时间和立管轴向坐标的函数,即:T=T(x,t);
为立管单位长度的质量;
c为阻尼系数;
k为剪切不均匀系数;
εT为张力引起的立管轴向应变;
q(x,t)为作用在立管上的流体荷载;
利用上述方程,计算立管弯曲振动的加速度、速度、位移、应力和应变随时间的变化。
2.如权利要求1所述的深水顶张式立管弯曲振动分析方法,其特征在于:该方法计算立管弯曲振动的加速度、速度、位移、应力和应变随时间变化的具体过程如下:
(1)将立管划分为若干个单元;
(2)将单元的位移函数表示为插值函数的形式:y=[N]{a},
其中,[N]是单元插值函数矩阵,{a}是单元节点位移向量;
(3)采用伽辽金方法,将深水顶张式立管弯曲振动分析模型的方程转换为矩阵方程如下:
式中:
阻尼矩阵[C]=α[M]+β[K],α,β为瑞雷阻尼系数;
{a}为位移矢量;
l为单元长度;
n为单元数量;
(4)基于立管的初始张力求出立管的初始轴向应变εT,并将其带入步骤(3)中刚度矩阵[K]的表达式,计算出立管的初始刚度矩阵;
(5)由步骤(3)中质量矩阵[M]、阻尼矩阵[C]和荷载向量{F}的表达式分别计算出立管的质量矩阵、阻尼矩阵和荷载向量;
此处,{ΔFi}是与时间增量Δt对应的荷载增量;
(7)由加速度增量、速度增量和位移增量按下式计算第一个时间增量后的加速度、速度和位移,此时i=0:
ai+1=ai+Δai
(8)根据步骤(7)中得到的位移计算出立管的应力和轴向应变;
(9)将步骤(8)中得到的轴向应变代入步骤(3)中刚度矩阵[K]的表达式,计算出第一个时间增量后的立管刚度矩阵;
(10)重复步骤(5)~(9),直至时间t达到要求的时长,即可计算出立管弯曲振动的加速度、速度、位移、应力和轴向应变随时间的变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110161221 CN102252897B (zh) | 2011-06-16 | 2011-06-16 | 一种深水顶张式立管弯曲振动分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110161221 CN102252897B (zh) | 2011-06-16 | 2011-06-16 | 一种深水顶张式立管弯曲振动分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102252897A true CN102252897A (zh) | 2011-11-23 |
CN102252897B CN102252897B (zh) | 2013-04-03 |
Family
ID=44980290
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110161221 Expired - Fee Related CN102252897B (zh) | 2011-06-16 | 2011-06-16 | 一种深水顶张式立管弯曲振动分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102252897B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103292774A (zh) * | 2013-05-31 | 2013-09-11 | 广州大学 | 一种桥梁动态挠度测量方法 |
CN106909710A (zh) * | 2017-01-11 | 2017-06-30 | 中国海洋大学 | 深水顶张式立管全耦合动力分析方法 |
CN107505210A (zh) * | 2017-07-04 | 2017-12-22 | 国家***第二海洋研究所 | 深水顶张式立管的复杂弯曲试验方法 |
CN110378060A (zh) * | 2019-07-26 | 2019-10-25 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
CN111307614A (zh) * | 2020-03-31 | 2020-06-19 | 广西交科集团有限公司 | 连续梁分段抗弯及剪切刚度测定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11294678A (ja) * | 1998-04-10 | 1999-10-29 | Denso Corp | 熱交換器の配管構造 |
CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
-
2011
- 2011-06-16 CN CN 201110161221 patent/CN102252897B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11294678A (ja) * | 1998-04-10 | 1999-10-29 | Denso Corp | 熱交換器の配管構造 |
CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
Non-Patent Citations (1)
Title |
---|
唐世振等: "《深水立管两向自由度涡激振动的数值分析》", 《振动与冲击》, vol. 29, no. 9, 31 December 2010 (2010-12-31), pages 206 - 211 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103292774A (zh) * | 2013-05-31 | 2013-09-11 | 广州大学 | 一种桥梁动态挠度测量方法 |
CN103292774B (zh) * | 2013-05-31 | 2016-04-13 | 广州大学 | 一种桥梁动态挠度测量方法 |
CN106909710A (zh) * | 2017-01-11 | 2017-06-30 | 中国海洋大学 | 深水顶张式立管全耦合动力分析方法 |
CN106909710B (zh) * | 2017-01-11 | 2020-03-17 | 中国海洋大学 | 深水顶张式立管全耦合动力分析方法 |
CN107505210A (zh) * | 2017-07-04 | 2017-12-22 | 国家***第二海洋研究所 | 深水顶张式立管的复杂弯曲试验方法 |
CN110378060A (zh) * | 2019-07-26 | 2019-10-25 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
CN110378060B (zh) * | 2019-07-26 | 2021-02-09 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
CN111307614A (zh) * | 2020-03-31 | 2020-06-19 | 广西交科集团有限公司 | 连续梁分段抗弯及剪切刚度测定方法 |
CN111307614B (zh) * | 2020-03-31 | 2022-06-10 | 广西交科集团有限公司 | 连续梁分段抗弯及剪切刚度测定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102252897B (zh) | 2013-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102252897B (zh) | 一种深水顶张式立管弯曲振动分析方法 | |
CN101539477B (zh) | 一种深水顶张式立管涡激振动与疲劳分析的方法 | |
WO2020113898A1 (zh) | 一种模拟深水无隔水管钻井钻柱纵-横-扭耦合非线性振动试验装置及其方法 | |
CN102222168A (zh) | 一种深水钻井立管参激横向振动分析方法 | |
CN107423486B (zh) | 一种在进船阶段浮托安装作业模型建模仿真方法 | |
CN102072805B (zh) | 阶梯剪切流下斜置立管的涡激振动旋转测试装置 | |
CN106909710A (zh) | 深水顶张式立管全耦合动力分析方法 | |
Galloway et al. | Experimental and numerical results of rotor power and thrust of a tidal turbine operating at yaw and in waves | |
CN102012306A (zh) | 双向剪切流下斜置立管的涡激振动旋转测试装置 | |
CN105335569A (zh) | 特高压输电线路多联v型绝缘子串力学特性仿真模拟方法 | |
CN101713707B (zh) | 深海钢悬链线立管波浪响应分析方法 | |
CN102410920B (zh) | 阶梯剪切流下顶部可运动斜置立管涡激振动旋转测试装置 | |
CN102445318B (zh) | 一种顶张式立管顺流向振动分析方法 | |
CN102353506B (zh) | 一种深水顶张式立管竖向振动分析方法 | |
CN103902754A (zh) | 一种深海钢悬链线立管触地点动力响应分析方法 | |
CN102288376A (zh) | 均匀流和阶梯流下顶部可运动竖置立管涡激振动旋转装置 | |
Yue et al. | Analysis of the wave load and dynamic response of a new semi-submersible wave-energy-powered aquaculture platform | |
CN102507082B (zh) | 一种深水立管的时域涡激升力确定方法 | |
CN102354333A (zh) | 一种钢悬链式立管出平面运动的分析方法 | |
CN102289536B (zh) | 一种深水钻井立管参数横力弯曲分析方法 | |
Samadi et al. | Hydrodynamic response simulation of Catenary mooring in the spar truss floating platform under Caspian Sea conditions | |
CN102507115B (zh) | 一种深水顶张式立管大变形弯曲振动的分析方法 | |
CN102117363A (zh) | 浮式Spar平台的扶正分析方法 | |
CN105675255A (zh) | 一种平台隔水管耦合水池模拟实验*** | |
Damblans et al. | Investigations on galloping of non-circular cross sections |
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: 20130403 Termination date: 20140616 |
|
EXPY | Termination of patent right or utility model |