CN108491615A - 一种三层顶张式立管动力响应的有限元分析方法 - Google Patents
一种三层顶张式立管动力响应的有限元分析方法 Download PDFInfo
- Publication number
- CN108491615A CN108491615A CN201810221139.7A CN201810221139A CN108491615A CN 108491615 A CN108491615 A CN 108491615A CN 201810221139 A CN201810221139 A CN 201810221139A CN 108491615 A CN108491615 A CN 108491615A
- Authority
- CN
- China
- Prior art keywords
- matrix
- standpipe
- outer tube
- formula
- pipe
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Earth Drilling (AREA)
Abstract
本发明公开了一种基于有限元分析的三层顶张紧式海洋立管的全耦合动力响应计算方法,属于深水油气开发工程管道技术分析领域。本发明采用有限元分析方法,分别建立三层管的有限元模型,将单元模型简化为欧拉梁单元;分别计算三个不同单元的质量矩阵与刚度矩阵,分别组合装配三层立管的质量矩阵与刚度矩阵;根据立管的质量矩阵与刚度矩阵,分别分析计算立管的频率与振型,进而选取扶正器的布置形式;根据立管所在海域环境,计算作用在外管上的水平拖曳力与水平惯性力;根据计算得出的油管与内管的横向位移响应,计算油管与内管各节点的应力、弯矩以及扶正器的横向作用力。本发明能正确分析出多层顶张紧式立管的横向受力与横向位移响应,具有概念明确、计算精度高等特点。
Description
技术领域
本发明公开了一种基于有限元分析的三层顶张紧式海洋立管的全耦合动力响应计算方法,属于深水油气开发工程管道技术分析领域。
背景技术
深水顶张式立管是联系海底油井与水面设施的纽带与桥梁,立管为海底油井与水面设施之间提供油、气、水的输送,是深水油气开发的特有装备,尽管浅水油气开发装备中也有立管,但与深水立管不可同日而语,其结构形式与水动力性能几乎发生了质的变化。
深水顶张式立管由油管和套管组成,油管与套管最常用的方式为管中管形式,最外层管为外管,外管的作用为保护内管与油管免受海流载荷、隔绝海水腐蚀、支撑油管的作用,在深水采油的过程中,伴随着天然气的产生,内管与油管之间形成的通道用来输送天然气,油管的作用为连接水面采油树和海底井口、输送石油。
目前的顶张式立管设计的分析方法都是按照截面刚度等效的方法将油管、内管和套管等效为外径与外套管相同(保证水动力荷载相等)的单层管来计算的,但是基于拉压刚度(EA)等效与弯曲刚度(EI)等效得到的单层管截面参数是不同的,因此,采用等效管替代三层管耦合的计算结果与立管各功能管的实际受力状态是有差异的,采用等效管的方法不能模拟出各层管的实际受力状态、实际横向位移、各层管之间的实际耦合作用。
因此,提供一种正确解决各层立管之间耦合作用的方法对于深水油气开发具有重要的意义。
发明内容
本发明的主要目的是提供一种正确解决各层立管之间全耦合动力响应分析的方法,以解决现有技术不能很好模拟各层管之间全耦合作用的问题。
本发明的三层顶张式立管动力响应的有限元分析方法依次按以下步骤进行:
(1)首先根据下列公式1建立立管的动力特性方程表达式:
式中:n=1,2,3分别表示油管、内管、外管;
[Mn]、[Kn]分别表示第n层管的质量矩阵与刚度矩阵;
ωn为第n层管的固有频率;
为第n层管的振型向量;
(2)根据下列公式2建立立管的动力响应方程表达式:
式中:n=1,2,3分别表示油管、内管、外管;
[Mn]、[Cn]、[Kn]分别表示第n层管的质量矩阵、阻尼矩阵与刚度矩阵;
{F}t为n层立管所受外力的大小;
(3)建立立管的有限元模型,根据下列公式求解每根立管的单元质量矩阵与刚度矩阵,进而装配形成每根立管的质量矩阵[Mn]与刚度矩阵[Kn]:
其中,N1N2N3N4的表达式如下:
N1(x)=1-3γ2+2γ3 N2(x)=(γ-2γ2+γ3)l
N3(x)=3γ2-2γ3 N4(x)=(-γ2+γ3
其中,γ的表达式如下:
式中:mnkn为单元的质量矩阵与刚度矩阵;
ln为第n层管划分的单元长度;
ρn、An分别为单元的密度与截面积,En、In分别为单元的弹性模量与截面惯性矩,当单元内的质量密度ρnAn和刚度EnIn为常量时,单元质量矩阵与刚度矩阵可以直接表达出来;
Ni(x)为Hermite形函数,该函数可保证相邻单元在边界上横向位移与转角连续;
xi为单元节点的起点位置,i表示单元编号;
(4)根据下式计算每层立管的阻尼矩阵:
[Cn]=a0[Mn]+a1[Kn]
式中:δ为阻尼比;ωni、ωnj分别为第n层管的i阶、j阶频率;an0和an1为第n层管的瑞雷阻尼系数;
(5)根据立管实际工作所处海域的海况,由以下公式计算出外层立管所受的海流作用力、顺流向涡激力的频率:
式中:CD′为顺流向脉动拖曳力系数,取CD′=0.05~0.1;CD为拖曳力系数,取CD=0.6~2.0;CM为惯性力系数,取CM=1.0~2.0;ρ为外管密度;D为外管直径;u为海水流速;为外管顺流向振动速度;为外管顺流向振动加速度;
ωs为涡泄圆频率;St为斯托哈尔数,取St=0.18~0.2;
(6)根据步骤(1)计算得出的各层管的固有频率与步骤(5)计算得出的顺流向涡激力的频率,找出各层管与顺流向涡激力的频率相近的固有频率与阵型。
(7)外层管与内层管之间,在外层管阵型峰值点处设置扶正器,内层管与油管之间,在内层管峰值处设置扶正器;
(8)在进行立管耦合动力分析前,应确定立管两端的链接情况,由此决定剔除立管质量模型矩阵与刚度模型矩阵中的哪些元素,当为铰接时,应剔除相应节点处位移矩阵元素,当为固定链接时,应剔除相应节点处位移与转角矩阵元素;
(9)三层管之间布置有扶正器,在扶正器处,当管与管发生碰撞时,扶正器处相邻管节点的位移、速度、加速度都是相同的,扶正器只传递横向作用力,通过变换立管质量矩阵、阻尼矩阵、刚度矩阵与外力矩阵,就可以进行耦合动力分析,以下公式给出各层管所受横向等效力的大小:
式中:{F3i}、{F2i}、{F1i}分别表示油管、内管与外管在单元节点i处所受等效横向作用力的大小;{fd(2i)}表示在节点i处,油管与内管之间扶正器横向作用力的大小;{fc(2i)}表示在节点i处,外管与内管之间扶正器横向作用力的大小;{fx(2i)}表示在节点i处,外管所受海流横向作用力的大小;
(10)通过联立求解步骤(9)中的三个公式,可以得出最外层管的动力响应表达式:
(11)由于在扶正器处,相邻两层管的横向响应同步,由步骤(10)可知,可通过将油管与内管扶正器处单元节点i处的质量矩阵、阻尼矩阵与刚度矩阵的2i行与2i列加到外管上,再将油管与内管质量矩阵、阻尼矩阵、刚度矩阵的2i行与2i列剔除,形成耦合矩阵进行求解;
(12)根据以下公式进行三层顶张紧式立管的耦合动力分析:
式中:M1、M2、M3分别为外管、内管与油管的质量矩阵;C1、C2、C3分别为外管、内管与油管的阻尼矩阵;K1、K2、K3分别为外管、内管与油管的刚度矩阵;{F1}、{F2}、{F3}分别为外管、内管与油管的受力;分别为外管、内管、油管的横向响应加速度;为横向响应速度;x1、x2、x3为横向响应位移;
(13)通过步骤(12)可以求得外管的节点位移,用插值法求解单元在外分布载荷作用下的最大弯曲变形,以此来判断外层立管与内层立管是否会发生碰撞;
(14)通过步骤(12)可以求得三根立管随时间的动力响应,由此来确定各层管各个单元节点的应力与弯矩,进行疲劳分析。
(15)公式3中,在计算外管的单元质量矩阵时,应考虑附加质量的影响:
式中:Ca为附加惯性力系数;ρw为海水密度;
(16)公式4中,在计算外管的单元刚度矩阵时,应考虑顶张紧力的影响:
式中:[KG]e为单元几何刚度矩阵;T为单元所受顶张紧力的大小。
本发明通过建立各层管的有限元模型,进而通过组合与变换各层管的总质量矩阵、阻尼矩阵、刚度矩阵和外载荷矩阵实现各层管之间的全耦合作用,正确的解决了现有技术存在的不足与缺陷,可以简单、准确、高效得计算出各层管的全耦合动力响应。
附图说明
图1本发明三层立管纵向结构示意图。
图2本发明三层立管力的传递示意图。
图3通过本发明方法计算得出的耦合动力响应结果图。
具体实施方式
下面通过附图和实施例对本发明作进一步详细说明,但本发明的保护范围不局限于所述内容。
实施例1:一种三层顶张式立管动力响应的有限元分析方法
具体内容如下:
海域水深300m,波高3m,波浪周期8s,海流顶部流速0.8m/s,立管底部海流流速0.2m/s。
表1立管设计参数
表2立管设计物理参数
(1)首先根据下列公式1建立立管的动力特性方程表达式:
式中:n=1,2,3分别表示油管、内管、外管;
[Mn]、[Kn]分别表示第n层管的质量矩阵与刚度矩阵;
ωn为第n层管的固有频率;
为第n层管的振型向量;
(2)根据下列公式2建立立管的动力响应方程表达式:
式中:n=1,2,3分别表示油管、内管、外管;
[Mn]、[Cn]、[Kn]分别表示第n层管的质量矩阵、阻尼矩阵与刚度矩阵;
{F}t为n层立管所受外力的大小;
(3)建立立管的有限元模型,根据下列公式求解每根立管的单元质量矩阵与刚度矩阵,进而装配形成每根立管的质量矩阵[Mn]与刚度矩阵[Kn]:
其中,N1N2N3N4的表达式如下:
N1(x)=1-3γ2+2γ3 N2(x)=(γ-2γ2+γ3)l
N3(x)=3γ2-2γ3 N4(x)=(-γ2+γ3)l
其中,γ的表达式如下:
式中:mnkn为单元的质量矩阵与刚度矩阵;
ln为第n层管划分的单元长度;
ρn、An分别为单元的密度与截面积,En、In分别为单元的弹性模量与截面惯性矩,当单元内的质量密度ρnAn和刚度EnIn为常量时,单元质量矩阵与刚度矩阵可以直接表达出来;
Ni(x)为Hermite形函数,该函数可保证相邻单元在边界上横向位移与转角连续;
xi为单元节点的起点位置,i表示单元编号;
(4)根据下式计算每层立管的阻尼矩阵:
[Cn]=a0[Mn]+a1[Kn]
式中:δ为阻尼比;ωni、ωnj分别为第n层管的i阶、j阶频率;an0和an1为第n层管的瑞雷阻尼系数;
(5)根据立管实际工作所处海域的海况,由以下公式计算出外层立管所受的海流作用力、顺流向涡激力的频率:
式中:CD′为顺流向脉动拖曳力系数,取CD′=0.05~0.1;CD为拖曳力系数,取CD=0.6~2.0;CM为惯性力系数,取CM=1.0~2.0;ρ为外管密度;D为外管直径;u为海水流速;为外管顺流向振动速度;为外管顺流向振动加速度;
ωs为涡泄圆频率;St为斯托哈尔数,取St=0.18~0.2;
(6)根据步骤(1)计算得出的各层管的固有频率与步骤(5)计算得出的顺流向涡激力的频率,找出各层管与顺流向涡激力的频率相近的固有频率与阵型。
(7)外层管与内层管之间,在外层管阵型峰值点处设置扶正器,内层管与油管之间,在内层管峰值处设置扶正器;
(8)在进行立管耦合动力分析前,应确定立管两端的链接情况,由此决定剔除立管质量模型矩阵与刚度模型矩阵中的哪些元素,当为铰接时,应剔除相应节点处位移矩阵元素,当为固定链接时,应剔除相应节点处位移与转角矩阵元素;
(9)如附图2所示,三层管之间布置有扶正器,在扶正器处,当管与管发生碰撞时,扶正器处相邻管节点的位移、速度、加速度都是相同的,扶正器只传递横向作用力,通过变换立管质量矩阵、阻尼矩阵、刚度矩阵与外力矩阵,就可以进行耦合动力分析,以下公式给出各层管所受横向等效力的大小:
式中:{F3i}、{F2i}、{F1i}分别表示油管、内管与外管在单元节点i处所受等效横向作用力的大小;{fd(2i)}表示在节点i处,油管与内管之间扶正器横向作用力的大小;{fc(2i)}表示在节点i处,外管与内管之间扶正器横向作用力的大小;{fx(2i)}表示在节点i处,外管所受海流横向作用力的大小;
(10)通过联立求解步骤(9)中的三个公式,可以得出最外层管的动力响应表达式:
(11)由于在扶正器处,相邻两层管的横向响应同步,由步骤(10)可知,可通过将油管与内管扶正器处单元节点i处的质量矩阵、阻尼矩阵与刚度矩阵的2i行与2i列加到外管上,再将油管与内管质量矩阵、阻尼矩阵、刚度矩阵的2i行与2i列剔除,形成耦合矩阵进行求解;
(12)根据以下公式进行三层顶张紧式立管的耦合动力分析:
式中:M1、M2、M3分别为外管、内管与油管的质量矩阵;C1、C2、C3分别为外管、内管与油管的阻尼矩阵;K1、K2、K3分别为外管、内管与油管的刚度矩阵;{F1}、{F2}、{F3}分别为外管、内管与油管的受力;分别为外管、内管、油管的横向响应加速度;为横向响应速度;x1、x2、x3为横向响应位移;
由于上式耦合求解过程过于繁琐,所以我们应当借助数学软件求解,Matlab软件具有解复杂方程且能绘制复杂函数曲线的功能,推荐使用其对模型求解与图像绘制,求解程序如下:
dsp=zeros(sdof,nt/dt);
vel=zeros(sdof,nt/dt);
acc=zeros(sdof,nt/dt);
dsp(:,1)=dsp0;
vel(:,1)=vel0;
acc(:,1)=acc0;
ekk=mm*4/(dt^2)+cc*2/dt+kk;
for in=1:nt/dt
efd=fd(:,in)+(mm*4/(dt^2)+cc*2/dt)*dsp(:,in)+(mm*4/dt+cc)*vel(:,in)+mm*acc(:,in);
dsp(:,in+1)=ekk\efd;
vel(:,in+1)=2*dsp(:,in+1)/dt-2*dsp(:,in)/dt-vel(:,in);
acc(:,in+1)=4*dsp(:,in+1)/(dt^2)-4*dsp(:,in)/(dt^2)-4*vel(:,in)/dt-acc(:,in);
end
(13)通过步骤(12)可以求得外管的节点位移,用插值法求解单元在外分布载荷作用下的最大弯曲变形,以此来判断外层立管与内层立管是否会发生碰撞;
(14)如附图3所示,通过步骤(12)可以求得三根立管随时间的动力响应,由此来确定各层管各个单元节点的应力与弯矩,进行疲劳分析。
(15)公式3中,在计算外管的单元质量矩阵时,应考虑附加质量的影响:
式中:Ca为附加惯性力系数;ρw为海水密度;
(16)公式4中,在计算外管的单元刚度矩阵时,应考虑顶张紧力的影响:
式中:[KG]e为单元几何刚度矩阵;T为单元所受顶张紧力的大小。
由附图3所示,该发明方法能够准确分析出三层立管的耦合动力响应,得出三层管之间合理的扶正器布置方式与每层管的横向位移响应曲线,进而分析出立管之间是否会发生碰撞以及疲劳分析。
Claims (2)
1.三层顶张式立管动力响应的有限元分析方法,其特征在于,包括以下步骤:
(1)首先根据下列公式1建立立管的动力特性方程表达式:
式中:n=1,2,3分别表示油管、内管、外管;
[Mn]、[Kn]分别表示第n层管的质量矩阵与刚度矩阵;
ωn为第n层管的固有频率;
为第n层管的振型向量;
(2)根据下列公式2建立立管的动力响应方程表达式:
式中:n=1,2,3分别表示油管、内管、外管;
[Mn]、[Cn]、[Kn]分别表示第n层管的质量矩阵、阻尼矩阵与刚度矩阵;
{F}t为n层立管所受外力的大小;
(3)建立立管的有限元模型,根据下列公式求解每根立管的单元质量矩阵与刚度矩阵,进而装配形成每根立管的总质量矩阵[Mn]与总刚度矩阵[Kn]:
其中,N1、N2、N3、N4的表达式如下:
N1(x)=1-3γ2+2γ3 N2(x)=(γ-2γ2+γ3)l
N3(x)=3γ2-2γ3 N4(x)=(-γ2+γ3)l
其中,γ的表达式如下:
式中:mn、kn为单元的质量矩阵与刚度矩阵;
ln为第n层管划分的单元长度;
ρn、An分别为单元的密度与截面积,En、In分别为单元的弹性模量与截面惯性矩,当单元内的质量密度ρnAn和刚度EnIn为常量时,单元质量矩阵与刚度矩阵可以直接表达出来;
Ni(x)为Hermite形函数,该函数可保证相邻单元在边界上横向位移与转角连续;
xi为单元节点的起点位置,i表示单元编号;
(4)根据下式计算每层立管的阻尼矩阵:
[Cn]=a0[Mn]+a1[Kn]
式中:δ为阻尼比;ωni、ωnj分别为第n层管的i阶、j阶频率;an0和an1为第n层管的瑞雷阻尼系数;
(5)根据立管实际工作所处海域的海况,由以下公式计算出外层立管所受的海流作用力、顺流向涡激力的频率:
式中:CD′为顺流向脉动拖曳力系数,取CD′=0.05~0.1;CD为拖曳力系数,取CD=0.6~2.0;CM为惯性力系数,取CM=1.0~2.0;ρ为外管密度;D为外管直径;u为海水流速;为外管顺流向振动速度;为外管顺流向振动加速度;
ωs为涡泄圆频率;St为斯托哈尔数,取St=0.18~0.2;
(6)根据步骤(1)计算得出的各层管的固有频率与步骤(5)计算得出的顺流向涡激力的频率,找出各层管与顺流向涡激力的频率相近的固有频率与阵型。
(7)外层管与内层管之间,在外层管阵型峰值点处设置扶正器,内层管与油管之间,在内层管峰值处设置扶正器;
(8)在进行立管耦合动力分析前,应确定立管两端的链接情况,由此决定剔除立管质量模型矩阵与刚度模型矩阵中的哪些元素,当为铰接时,应剔除相应节点处位移矩阵元素,当为固定链接时,应剔除相应节点处位移与转角矩阵元素;
(9)三层管之间布置有扶正器,在扶正器处,当管与管发生碰撞时,扶正器处相邻管节点的位移、速度、加速度都是相同的,扶正器只传递横向作用力,通过变换立管质量矩阵、阻尼矩阵、刚度矩阵与外力矩阵,就可以进行耦合动力分析,以下公式给出各层管所受等效横向力的大小:
式中:{F3i}、{F2i}、{F1i}分别表示油管、内管与外管在单元节点i处所受等效横向作用力的大小;{fd(2i)}表示在节点i处,油管与内管之间扶正器横向作用力的大小;{fc(2i)}表示在节点i处,外管与内管之间扶正器横向作用力的大小;{fx(2i)}表示在节点i处,外管所受海流横向作用力的大小;
(10)通过联立求解步骤(9)中的三个公式,可以得出最外层管的动力响应表达式:
(11)由于在扶正器处,相邻两层管的横向响应同步,由步骤(10)可知,可通过将油管与内管扶正器处单元节点i处的质量矩阵、阻尼矩阵与刚度矩阵的2i行与2i列加到外管上,再将油管与内管质量矩阵、阻尼矩阵、刚度矩阵的2i行与2i列剔除,形成耦合矩阵进行求解;
(12)根据以下公式进行三层顶张紧式立管的耦合动力分析:
式中:M1、M2、M3分别为外管、内管与油管的质量矩阵;C1、C2、C3分别为外管、内管与油管的阻尼矩阵;K1、K2、K3分别为外管、内管与油管的刚度矩阵;{F1}、{F2}、{F3}分别为外管、内管与油管的受力;分别为外管、内管、油管的横向响应加速度;为横向响应速度;x1、x2、x3为横向响应位移;
(13)通过步骤(12)可以求得外管的节点位移,用插值法求解单元在外分布载荷作用下的最大弯曲变形,以此来判断外层立管与内层立管是否会发生碰撞;
(14)通过步骤(12)可以求得三根立管随时间的动力响应,由此来确定各层管各个单元节点的应力与弯矩,进行疲劳分析。
2.根据权利要求1所述的三层顶张式立管动力响应的有限元分析方法,其特征在于:
(1)公式3中,在计算外管的单元质量矩阵时,应考虑附加质量的影响:
式中:Ca为附加惯性力系数;ρw为海水密度;
(2)公式4中,在计算外管的单元刚度矩阵时,应考虑顶张紧力的影响:
式中:[KG]e为单元几何刚度矩阵;T为单元所受顶张紧力的大小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810221139.7A CN108491615B (zh) | 2018-03-17 | 2018-03-17 | 一种三层顶张式立管动力响应的有限元分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810221139.7A CN108491615B (zh) | 2018-03-17 | 2018-03-17 | 一种三层顶张式立管动力响应的有限元分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108491615A true CN108491615A (zh) | 2018-09-04 |
CN108491615B CN108491615B (zh) | 2022-09-02 |
Family
ID=63339775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810221139.7A Active CN108491615B (zh) | 2018-03-17 | 2018-03-17 | 一种三层顶张式立管动力响应的有限元分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108491615B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110378060A (zh) * | 2019-07-26 | 2019-10-25 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
CN110727977A (zh) * | 2019-05-30 | 2020-01-24 | 浙江海洋大学 | 一种顶张式立管耦合运动的数值模拟方法 |
CN113673126A (zh) * | 2021-07-28 | 2021-11-19 | 中国石油大学(北京) | 一种钻井用多层套管环空带压有限元计算方法和装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040031614A1 (en) * | 2002-04-26 | 2004-02-19 | Kleinhans John W. | Marine bottom tensioned riser and method |
US20060080069A1 (en) * | 2004-10-12 | 2006-04-13 | William Fujimoto | System for analyzing fastener loads |
CN101475048A (zh) * | 2009-01-13 | 2009-07-08 | 中国海洋大学 | 一种新型海洋深水浮筒平台 |
CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
CN201580542U (zh) * | 2009-12-02 | 2010-09-15 | 中国海洋大学 | 一种三角形桁架单柱平台 |
CN102353506A (zh) * | 2011-06-16 | 2012-02-15 | 中国海洋大学 | 一种深水顶张式立管竖向振动分析方法 |
CN102445318A (zh) * | 2011-09-30 | 2012-05-09 | 中国海洋大学 | 一种顶张式立管顺流向振动分析方法 |
CN106909710A (zh) * | 2017-01-11 | 2017-06-30 | 中国海洋大学 | 深水顶张式立管全耦合动力分析方法 |
-
2018
- 2018-03-17 CN CN201810221139.7A patent/CN108491615B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040031614A1 (en) * | 2002-04-26 | 2004-02-19 | Kleinhans John W. | Marine bottom tensioned riser and method |
US20060080069A1 (en) * | 2004-10-12 | 2006-04-13 | William Fujimoto | System for analyzing fastener loads |
CN101475048A (zh) * | 2009-01-13 | 2009-07-08 | 中国海洋大学 | 一种新型海洋深水浮筒平台 |
CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
CN201580542U (zh) * | 2009-12-02 | 2010-09-15 | 中国海洋大学 | 一种三角形桁架单柱平台 |
CN102353506A (zh) * | 2011-06-16 | 2012-02-15 | 中国海洋大学 | 一种深水顶张式立管竖向振动分析方法 |
CN102445318A (zh) * | 2011-09-30 | 2012-05-09 | 中国海洋大学 | 一种顶张式立管顺流向振动分析方法 |
CN106909710A (zh) * | 2017-01-11 | 2017-06-30 | 中国海洋大学 | 深水顶张式立管全耦合动力分析方法 |
Non-Patent Citations (4)
Title |
---|
HUANG WEIPING .ETAL: ""Vortex-Induced Vibration of Two Parallel Risers Experimental Test and Numerical Simulation"", 《中国海洋大学学报:英文版》 * |
吴学敏 等: ""深水顶张式立管参数振动与涡激振动耦合振动分析方法研究"", 《中国海上油气》 * |
孙传栋 等: ""深水顶张式生产立管动力响应的有限元计算"", 《海洋石油》 * |
毛国辉 等: ""车-桥耦合振动的联合方程分析方法"", 《重庆交通大学学报(自然科学版)》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110727977A (zh) * | 2019-05-30 | 2020-01-24 | 浙江海洋大学 | 一种顶张式立管耦合运动的数值模拟方法 |
CN110727977B (zh) * | 2019-05-30 | 2023-06-23 | 浙江海洋大学 | 一种顶张式立管耦合运动的数值模拟方法 |
CN110378060A (zh) * | 2019-07-26 | 2019-10-25 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
CN110378060B (zh) * | 2019-07-26 | 2021-02-09 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
CN113673126A (zh) * | 2021-07-28 | 2021-11-19 | 中国石油大学(北京) | 一种钻井用多层套管环空带压有限元计算方法和装置 |
CN113673126B (zh) * | 2021-07-28 | 2024-02-13 | 中国石油大学(北京) | 一种钻井用多层套管环空带压有限元计算方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108491615B (zh) | 2022-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fan et al. | Dynamic analysis of a hang-off drilling riser considering internal solitary wave and vessel motion | |
CN108491615A (zh) | 一种三层顶张式立管动力响应的有限元分析方法 | |
Zou et al. | Response of a submerged floating tunnel subject to flow-induced vibration | |
Lee et al. | Dynamic response analysis of submerged floating tunnels by wave and seismic excitations | |
CN106909710A (zh) | 深水顶张式立管全耦合动力分析方法 | |
Jin et al. | The effect of key design parameters on the global performance of submerged floating tunnel under target wave and earthquake excitations | |
CN106934138A (zh) | 内输气液两相段塞流的海底管道悬跨段振动特性分析方法 | |
Fu et al. | Experimental investigation on VIV of the flexible model under full scale Re number | |
Lakshmynarayanana et al. | Hydroelastic analysis of a flexible barge in regular waves using coupled CFD-FEM modelling | |
CN106202774A (zh) | 一种波浪作用过程中深水网箱浮架变形的数值模拟方法 | |
Jing et al. | Coupled dynamic modeling of a moored floating platform with risers | |
CN102507082B (zh) | 一种深水立管的时域涡激升力确定方法 | |
CN113283138B (zh) | 一种基于深度学习的深海养殖平台动力响应分析方法 | |
Jin et al. | Coupled time-domain hydro-elastic simulation for submerged floating tunnel under wave excitations | |
Adaalah et al. | Analysis of mooring line tension in floating collar net cage | |
CN102507115B (zh) | 一种深水顶张式立管大变形弯曲振动的分析方法 | |
Jin et al. | Digital Twin Method for Global Motion and Stress Monitoring of a Steel Lazy Wave Riser | |
Spraul et al. | Parameter calibration in dynamic simulations of power cables in shallow water to improve fatigue damage estimation | |
Kamble et al. | 3D VIV fatigue analysis using CFD simulation for long marine risers | |
Wang et al. | Analysis of marine risers subjected to shoal/deep water in the installation process | |
Grządziela | Ship impact modeling of underwater explosion | |
CN110727977B (zh) | 一种顶张式立管耦合运动的数值模拟方法 | |
Chatjigeorgiou et al. | Numerical and experimental investigation of the wave loading on a three-legged offshore wind turbine jacket platform | |
Johnson et al. | Simplified modeling for the nonlinear response of tension leg platforms in deep water | |
Gudmestad et al. | The influence of internal solitary waves on steel catenary riser and steel lazy wave risers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |