CN106909710B - 深水顶张式立管全耦合动力分析方法 - Google Patents
深水顶张式立管全耦合动力分析方法 Download PDFInfo
- Publication number
- CN106909710B CN106909710B CN201710016913.6A CN201710016913A CN106909710B CN 106909710 B CN106909710 B CN 106909710B CN 201710016913 A CN201710016913 A CN 201710016913A CN 106909710 B CN106909710 B CN 106909710B
- Authority
- CN
- China
- Prior art keywords
- pipe
- unit
- models
- riser
- nodes
- 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.)
- Expired - Fee Related
Links
Images
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
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Earth Drilling (AREA)
Abstract
本发明涉及油气开发工程管道分析技术领域,尤其涉及深水顶张式立管全耦合动力分析方法。其分别建立油管、内层立管和外层立管的有限元模型;分别计算每根管模型的质量矩阵和刚度矩阵元素,形成每根管模型的质量矩阵阵和刚度矩阵,计算每根管的固有频率和振型,根据实际工程海域的海洋环境条件下计算涡激升力和脉动拖曳力的频率。本发明正确地描述了顶张式立管横向振动的实际受力和变形状态,提出了正确的分析方法。
Description
技术领域
本发明涉及油气开发工程管道分析技术领域,尤其涉及深水顶张式立管全耦合动力分析方法。
背景技术
深水顶张式立管是一种由三根不同直径的钢管套装在一起的同心管中管结构,最内层的是油管、中间的和最外层的为立管。
目前常用的顶张式立管其油管为4.5英寸,中间立管为8~9英寸,外层立管为12~13英寸。水深越大,中间和外层立管的直径越大。对于这样的管中管结构,波浪和海流是直径作用在外层立管表面的,而内层立管和油管则通过分别设置在两层管之间的导向环而受到外部荷载的作用。
深水顶张式立管在深水油气开发中的作用是采油,其中的油管用于将井下的原油输送至平台上,因此,在承受其他外部荷载(与内、外层立管承受的荷载类型相同)的同时还承受原油的压力作用。为了原油的正常输送,油管的直径与井下的井管相同,而井管的直径受抽油泵的限制,尺寸较小,因此,油管的尺寸较小。
采油的同时,伴有天然气的产生,因此,内层立管是作为气举线使用的,即它与油管之间形成的环形空间输送井下的天然气。
现阶段,顶张式立管的设计分析采用等效管模型,采用弯曲刚度等效的原则将三根管等效为一根管。而这种设计分析方式存在着拉压刚度不等效、压力不等效、动力响应不等效的分析缺陷,严重影响立管使用的安全性。
发明内容
本发明要解决的技术问题是如何克服现有技术的不足,提供一种深水顶张式立管全耦合动力分析方法
本发明为实现上述目的采用的技术方案是:深水顶张式立管全耦合动力分析方法,包括以下步骤:
(1)分别建立油管、内层立管和外层立管的有限元模型;由下列公式1分别计算每根管模型的质量矩阵和刚度矩阵元素,形成每根管模型的质量矩阵阵[M]和刚度矩阵[K]:
式中:M,K分别为单元质量矩阵和刚度矩阵的元素;
EI为截面抗弯刚度;N为单元插值函数;
i,j为结点自由度,当i,j取值为1时,表示轴向变形自由度;当i,j取值为2、3时,表示截面两个主轴方向的变形自由度,所述变形自由度为弯曲变形引起的挠度;当i,j取值分别为4、5、6时,表示绕x、y、z轴的转动自由度,所述转动自由度为转角变形;
n=1,2,3分别表示油管、内层立管和外层立管;
(2)由下列公式2分别计算每根管的固有频率ω和振型{φ}:
([K]-ω2[M]){φ}={0}
(3)根据实际工程海域的海洋环境条件下由下列公式3计算涡激升力和脉动拖曳力的频率:
式中:fL,fD分别为涡激升力和脉动拖曳力的频率;
St为斯托哈尔数,范围为0.18~0.2;Uw为海水流速;Douter为外层立管直径;
(4)分别找出三根管与波浪力频率和涡激力频率最相近的频率及振型,所述涡激力包括升力和拖曳力;
(5)分别找出外层立管的振型节点、油管和内层立管的峰值点,并将这些点设置为有限元模型的单元节点,若这些结点中的某些结点之间的距离小于外层管直径,则删除距离小于外层管直径的模型结点,以确保单元长度大于管模型直径;
(6)步骤(5)中得到的单元结点作为立管中导向环位置,赋以同一个导向环位置的三根管模型结点相同的结点编号,使模型两端的结点也具有相同的结点编号;
(7)根据步骤(6)中结点之间的最小距离对三根管模型进行单元划分,由此得到的单元结点具有相同的结点坐标,但不能具有相同的结点编号;
(8)依据上述步骤得到的关于三根管模型的几何和物理参数,通过计算机程序计算弯曲刚度和拉压刚度;
(9)继续对三根管模型施加约束及载荷,将管模型的下端的所有自由度进行约束,上端施加竖直方向的弹性约束,波浪和海流荷载仅施加在外层立管模型上;
(10)对上述施加约束和载荷后的三个管模型进行有限元分析,分别得到三根管的结点位移、单元应力应变,并分别进行强度和疲劳校核;
(11)利用单元插值函数计算每个单元的最大弯曲变形,分别校核三根管模型之间是否发生碰撞,完成对深水顶张式立管的动力分析。
进一步,外层立管的固有频率计算应考虑海水的附加质量,即在计算质量矩阵元素时,在公式1的第一式中增加附加质量系数:
本发明正确地描述了顶张式立管横向振动的实际受力和变形状态,提出了正确的分析方法。
附图说明
图1本发明立管径向截面结构示意图。
图2本发明立管局部剖视图。
图3立管振型曲线和节点位置示意图。
图中,1、油管;2、内层立管;3、外层立管;4、导向环。
具体实施方式
下面结合附图,对本发明进一步说明。
如图1所示,表示深水顶张向截面图,由图可以看出深水顶张式立管是一种有三根不同直径的钢管套装在一起的同心管中管结构,最内层的是油管1,中间的和最为层的为立管。油管1用于采油,内层立管2是气举线,外层立管3为隔水管,是承受海洋环境荷载的主要结构。为了防止管与管之间发生碰撞并传递载荷,油管与内层立管之间以及内层立管与外层立管之间设置了导向环4。
目前常用的顶张式立管其油管为4.5英寸,中间立管为8~9英寸,外层立管为12~13英寸。水深越大,中间和外层立管的直径越大。对于这样的管中管结构,波浪和海流是直接作用在外层立管表面的,而内层立管和油管则通过分别设置在两层管之间的导向环而受到外部载荷的作用。
本发明采用全耦合的分析方法来计算结构的动力响应,是解决目前顶张式立管设计分析中存在的计算模型与实际结构的受力和变形状态不等效的问题,从而更准确的设计导向环的数量和位置,达到减小管与管碰撞摩擦的风险,改善油管和立管的受力与变形。
基于上述问题,本发明提出的技术方案是:深水顶张式立管全耦合动力分析方法,包括以下步骤:
(1)分别建立油管、内层立管和外层立管的有限元模型;由下列公式1分别计算每根管模型的质量矩阵和刚度矩阵元素,形成每根管模型的质量矩阵阵[M]和刚度矩阵[K]:
式中:M,K分别为单元质量矩阵和刚度矩阵的元素;
EI为截面抗弯刚度;N为单元插值函数;
i,j为结点自由度,当i,j取值为1时,表示轴向变形自由度;当i,j取值为2、3时,表示截面两个主轴方向的变形自由度,所述变形自由度为弯曲变形引起的挠度;当i,j取值分别为4、5、6时,表示绕x、y、z轴的转动自由度,所述转动自由度为转角变形;
n=1,2,3分别表示油管、内层立管和外层立管;
(2)由下列公式2分别计算每根管的固有频率ω和振型{φ}:
([K]-ω2[M]){φ}={0}
(3)根据实际工程海域的海洋环境条件下由下列公式3计算涡激升力和脉动拖曳力的频率:
式中:fL,fD分别为涡激升力和脉动拖曳力的频率;
St为斯托哈尔数,范围为0.18~0.2;Uw为海水流速;Douter为外层立管直径;
(4)分别找出三根管与波浪力频率和涡激力频率最相近的频率及振型,所述涡激力包括升力和拖曳力;
(5)分别找出外层立管的振型节点、油管和内层立管的峰值点,如图3所示,并将这些点设置为有限元模型的单元节点,若这些结点中的某些结点之间的距离小于外层管直径,则删除距离小于外层管直径的模型结点,以确保单元长度大于管模型直径;
(6)步骤(5)中得到的单元结点作为立管中导向环位置,赋以同一个导向环位置的三根管模型结点相同的结点编号,使模型两端的结点也具有相同的结点编号;
(7)根据步骤(6)中结点之间的最小距离对三根管模型进行单元划分,由此得到的单元结点具有相同的结点坐标,但不能具有相同的结点编号;
(8)依据上述步骤得到的关于三根管模型的几何和物理参数,通过计算机程序计算弯曲刚度和拉压刚度;
(9)继续对三根管模型施加约束及载荷,将管模型的下端的所有自由度进行约束,上端施加竖直方向的弹性约束,波浪和海流荷载仅施加在外层立管模型上;
(10)对上述施加约束和载荷后的三个管模型进行有限元分析,分别得到三根管的结点位移、单元应力应变,并分别进行强度和疲劳校核;
(11)利用单元插值函数计算每个单元的最大弯曲变形,分别校核三根管模型之间是否发生碰撞,完成对深水顶张式立管的动力分析。
需要提出的是,外层立管的固有频率计算应考虑海水的附加质量,即在计算质量矩阵元素时,在公式1的第一式中增加附加质量系数:
基于上述方案,本发明能够解决传统的分析方法中存在的拉压刚度不等效、压力不等效、动力响应不等效的问题,而且能够为导向环的设置提供最优方案。
上述实施例只是为了说明本发明的技术构思及特点,其目的是在于让本领域内的普通技术人员能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡是根据本发明内容的实质所作出的等效的变化或修饰,都应涵盖在本发明的保护范围内。
Claims (2)
1.深水顶张式立管全耦合动力分析方法,其特征在于,包括以下步骤:
(1)分别建立油管、内层立管和外层立管的有限元模型;由下列公式1分别计算每根管模型的质量矩阵和刚度矩阵元素,形成每根管模型的质量矩阵阵[M]和刚度矩阵[K]:
式中:M,K分别为单元质量矩阵和刚度矩阵的元素;
EI为截面抗弯刚度;N为单元插值函数;
i,j为结点自由度,当i,j取值为1时,表示轴向变形自由度;当i,j取值为2、3时,表示截面两个主轴方向的变形自由度,所述变形自由度为弯曲变形引起的挠度;当i,j取值分别为4、5、6时,表示绕x、y、z轴的转动自由度,所述转动自由度为转角变形;
n=1,2,3分别表示油管、内层立管和外层立管;
(2)由下列公式2分别计算每根管的固有频率ω和振型{φ}:
([K]-ω2[M]){φ}={0}
(3)根据实际工程海域的海洋环境条件下由下列公式3计算涡激升力和脉动拖曳力的频率:
式中:fL,fD分别为涡激升力和脉动拖曳力的频率;
St为斯托哈尔数,范围为0.18~0.2;Uw为海水流速;Douter为外层立管直径;
(4)分别找出三根管与波浪力频率和涡激力频率最相近的频率及振型,所述涡激力包括升力和拖曳力;
(5)分别找出外层立管的振型节点、油管和内层立管的峰值点,并将这些点设置为有限元模型的单元节点,若这些结点中的某些结点之间的距离小于外层管直径,则删除距离小于外层管直径的模型结点,以确保单元长度大于管模型直径;
(6)步骤(5)中得到的单元结点作为立管中导向环位置,赋以同一个导向环位置的三根管模型结点相同的结点编号,使模型两端的结点也具有相同的结点编号;
(7)根据步骤(6)中结点之间的最小距离对三根管模型进行单元划分,由此得到的单元结点具有相同的结点坐标,但不能具有相同的结点编号;
(8)依据上述步骤得到的关于三根管模型的几何和物理参数,通过计算机程序计算弯曲刚度和拉压刚度;
(9)继续对三根管模型施加约束及载荷,将管模型的下端的所有自由度进行约束,上端施加竖直方向的弹性约束,波浪和海流荷载仅施加在外层立管模型上;
(10)对上述施加约束和载荷后的三个管模型进行有限元分析,分别得到三根管的结点位移、单元应力应变,并分别进行强度和疲劳校核;
(11)利用单元插值函数计算每个单元的最大弯曲变形,分别校核三根管模型之间是否发生碰撞,完成对深水顶张式立管的动力分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710016913.6A CN106909710B (zh) | 2017-01-11 | 2017-01-11 | 深水顶张式立管全耦合动力分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710016913.6A CN106909710B (zh) | 2017-01-11 | 2017-01-11 | 深水顶张式立管全耦合动力分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106909710A CN106909710A (zh) | 2017-06-30 |
CN106909710B true CN106909710B (zh) | 2020-03-17 |
Family
ID=59207334
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710016913.6A Expired - Fee Related CN106909710B (zh) | 2017-01-11 | 2017-01-11 | 深水顶张式立管全耦合动力分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106909710B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108229043B (zh) * | 2018-01-12 | 2021-05-11 | 中国海洋大学 | 考虑涡激效应的深海spar型浮式风机疲劳损伤分析方法 |
CN108491615B (zh) * | 2018-03-17 | 2022-09-02 | 中国海洋大学 | 一种三层顶张式立管动力响应的有限元分析方法 |
CN109165447B (zh) * | 2018-08-27 | 2022-12-13 | 武汉理工大学 | 一种空间钢管结构壁厚损失的评价方法及系统 |
CN110727977B (zh) * | 2019-05-30 | 2023-06-23 | 浙江海洋大学 | 一种顶张式立管耦合运动的数值模拟方法 |
CN110378060B (zh) * | 2019-07-26 | 2021-02-09 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
CN114290644A (zh) * | 2021-12-17 | 2022-04-08 | 烟台市佶楷矿山设备科技有限公司 | 一种非金属长托辊的制造方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
CN102252897A (zh) * | 2011-06-16 | 2011-11-23 | 中国海洋大学 | 一种深水顶张式立管弯曲振动分析方法 |
CN102445318A (zh) * | 2011-09-30 | 2012-05-09 | 中国海洋大学 | 一种顶张式立管顺流向振动分析方法 |
CN102507115A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种深水顶张式立管大变形弯曲振动的分析方法 |
CN102507082A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种深水立管的时域涡激升力确定方法 |
US8764346B1 (en) * | 2010-06-07 | 2014-07-01 | Nagan Srinivasan | Tension-based tension leg platform |
-
2017
- 2017-01-11 CN CN201710016913.6A patent/CN106909710B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
US8764346B1 (en) * | 2010-06-07 | 2014-07-01 | Nagan Srinivasan | Tension-based tension leg platform |
CN102252897A (zh) * | 2011-06-16 | 2011-11-23 | 中国海洋大学 | 一种深水顶张式立管弯曲振动分析方法 |
CN102445318A (zh) * | 2011-09-30 | 2012-05-09 | 中国海洋大学 | 一种顶张式立管顺流向振动分析方法 |
CN102507115A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种深水顶张式立管大变形弯曲振动的分析方法 |
CN102507082A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种深水立管的时域涡激升力确定方法 |
Non-Patent Citations (2)
Title |
---|
深水顶张式生产立管动力响应的有限元计算;孙传栋 等;《海洋石油》;20090922;正文第85-90页 * |
深水顶张式立管参数振动与涡激振动祸合振动;吴学敏 等;《中国海上油气》;20140830;正文第100-105页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106909710A (zh) | 2017-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106909710B (zh) | 深水顶张式立管全耦合动力分析方法 | |
Wang et al. | A nonlinear model for deepwater steel lazy-wave riser configuration with ocean current and internal flow | |
Amaechi et al. | Numerical studies on CALM buoy motion responses and the effect of buoy geometry cum skirt dimensions with its hydrodynamic waves-current interactions | |
CN106934138B (zh) | 内输气液两相段塞流的海底管道悬跨段振动特性分析方法 | |
Sun et al. | Global analysis of a flexible riser | |
Wang et al. | Mathematical model of steel lazy-wave riser abandonment and recovery in deepwater | |
US9777539B2 (en) | Composite component deployment configurations | |
CN108491615B (zh) | 一种三层顶张式立管动力响应的有限元分析方法 | |
Wang et al. | A nonlinear mechanical model for deepwater steel lazy-wave riser transfer process during installation | |
CN104422424B (zh) | 一种隔水管偏移及偏移方向监测方法 | |
Lou et al. | Experimental and numerical study of dynamic performance of CVAR subjected to regular wave and platform motion | |
Greene et al. | The influence of drilling rig and riser system selection on wellhead fatigue loading | |
CN102507115B (zh) | 一种深水顶张式立管大变形弯曲振动的分析方法 | |
Yue et al. | Tension and expansion analysis of pipe-in-pipe risers: part B finite element modeling | |
CN105952389B (zh) | 一种布置有表面织构的浮力块及其安全控制方法 | |
Tang et al. | Time-domain nonlinear wave-current interaction with a steep wave riser considering internal flow effect | |
Bomfimsilva et al. | On the feasibility of a novel concept for a free standing riser | |
Buberg | Design and analysis of steel catenary riser systems for deep waters | |
Yu et al. | Nonlinear dynamic response of touchdown zone for steel catenary riser under multiple internal solitary waves | |
Chen et al. | Numerical analysis of configuration for steel lazy-wave riser in deepwater | |
Yan et al. | Reliability optimization design of the steel tube umbilical cross section based on particle swarm algorithm | |
Lei et al. | Calculation method of residual fatigue life of in-service unbonded flexible pipe | |
Lobo et al. | Case Study: Use of Digital Twins to Mitigate Risks in Pull-Out Operations of Damaged Risers | |
WO2013167710A2 (en) | Riser assembly and method | |
CN110727977B (zh) | 一种顶张式立管耦合运动的数值模拟方法 |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Huang Weiping Inventor after: Song Hong Inventor after: Chang Shuang Inventor after: Fu Xuepeng Inventor before: Huang Weiping Inventor before: Song Hong Inventor before: Fu Xuepeng |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200317 Termination date: 20210111 |
|
CF01 | Termination of patent right due to non-payment of annual fee |