CN106909710A - 深水顶张式立管全耦合动力分析方法 - Google Patents
深水顶张式立管全耦合动力分析方法 Download PDFInfo
- Publication number
- CN106909710A CN106909710A CN201710016913.6A CN201710016913A CN106909710A CN 106909710 A CN106909710 A CN 106909710A CN 201710016913 A CN201710016913 A CN 201710016913A CN 106909710 A CN106909710 A CN 106909710A
- Authority
- CN
- China
- Prior art keywords
- node
- pipe
- standpipe
- outer layer
- model
- 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
- 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)
- Earth Drilling (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (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的第一式中增加附加质量系数:
式中:为外层立管单元长度的附加质量,其中,Ca为惯性力系数,一般取1.0;ρw为海水密度,Aouter为外层立管的横截面积。
本发明正确地描述了顶张式立管横向振动的实际受力和变形状态,提出了正确的分析方法。
附图说明
图1本发明立管径向截面结构示意图。
图2本发明立管局部剖视图。
图3立管振型曲线和节点位置示意图。
具体实施方式
下面结合附图,对本发明进一步说明。
如图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的第一式中增加附加质量系数:
式中:为外层立管单元长度的附加质量,其中,Ca为惯性力系数,一般取1.0;ρw为海水密度,Aouter为外层立管的横截面积。
基于上述方案,本发明能够解决传统的分析方法中存在的拉压刚度不等效、压力不等效、动力响应不等效的问题,而且能够为导向环的设置提供最优方案。
上述实施例只是为了说明本发明的技术构思及特点,其目的是在于让本领域内的普通技术人员能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡是根据本发明内容的实质所作出的等效的变化或修饰,都应涵盖在本发明的保护范围内。
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计算涡激升力和脉动拖曳力的频率:
fD=2fL
式中:fL,fD分别为涡激升力和脉动拖曳力的频率;
St为斯托哈尔数,一般取0.18~0.2;Uw为海水流速;Douter为外层立管直径;
(4)分别找出三根管与波浪力频率和涡激力(升力和拖曳力)频率最相近的频率及振型;
(5)分别找出外层立管的振型节点、油管和内层立管的峰值点,并将这些点设置为有限元模型的单元节点,若这些结点中的某些结点之间的距离小于外层管直径,则删除距离较小的模型结点,以确保单元长度大于管模型直径;
(6)步骤(5)中得到的单元结点作为立管中导向环位置,赋以同一个导向环位置的三根管模型结点相同的结点编号,使模型两端的结点也具有相同的结点编号;
(7)根据步骤(6)中结点之间的最小距离对三根管模型进行单元划分,由此得到的单元结点可以具有相同的结点坐标,但不能具有相同的结点编号;
(8)依据上述步骤得到的关于三根管模型的几何和物理参数,通过计算机程序计算弯曲刚度和拉压刚度;
(9)继续对三根管模型施加约束及载荷,将管模型的下端的所有自由度进行约束,上端施加竖直方向的弹性约束,波浪和海流荷载仅施加在外层立管模型上;
(10)对上述施加约束和载荷后的三个管模型进行有限元分析,分别得到三根管的结点位移、单元应力应变,并分别进行强度和疲劳校核;
(11)利用单元插值函数计算每个单元的最大弯曲变形,分别校核三根管模型之间是否发生碰撞,完成对深水顶张式立管的动力分析。
2.根据权利要求1所述的深水顶张式立管全耦合动力分析方法,其特征在于:外层立管的固有频率计算应考虑海水的附加质量,即在计算质量矩阵元素时,在公式1的第一式中增加附加质量系数:
式中:为外层立管单元长度的附加质量,其中,Ca为惯性力系数,一般取1.0;ρw为海水密度,Aouter为外层立管的横截面积。
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 true CN106909710A (zh) | 2017-06-30 |
CN106909710B 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) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108229043A (zh) * | 2018-01-12 | 2018-06-29 | 中国海洋大学 | 考虑涡激效应的深海spar型浮式风机疲劳损伤分析方法 |
CN108491615A (zh) * | 2018-03-17 | 2018-09-04 | 中国海洋大学 | 一种三层顶张式立管动力响应的有限元分析方法 |
CN109165447A (zh) * | 2018-08-27 | 2019-01-08 | 武汉理工大学 | 一种空间钢管结构壁厚损失的评价方法及*** |
CN110378060A (zh) * | 2019-07-26 | 2019-10-25 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
CN110727977A (zh) * | 2019-05-30 | 2020-01-24 | 浙江海洋大学 | 一种顶张式立管耦合运动的数值模拟方法 |
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 | 中国海洋大学 | 一种顶张式立管顺流向振动分析方法 |
CN102507082A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种深水立管的时域涡激升力确定方法 |
CN102507115A (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 | 中国海洋大学 | 一种顶张式立管顺流向振动分析方法 |
CN102507082A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种深水立管的时域涡激升力确定方法 |
CN102507115A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种深水顶张式立管大变形弯曲振动的分析方法 |
Non-Patent Citations (2)
Title |
---|
吴学敏 等: "深水顶张式立管参数振动与涡激振动祸合振动", 《中国海上油气》 * |
孙传栋 等: "深水顶张式生产立管动力响应的有限元计算", 《海洋石油》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108229043A (zh) * | 2018-01-12 | 2018-06-29 | 中国海洋大学 | 考虑涡激效应的深海spar型浮式风机疲劳损伤分析方法 |
CN108229043B (zh) * | 2018-01-12 | 2021-05-11 | 中国海洋大学 | 考虑涡激效应的深海spar型浮式风机疲劳损伤分析方法 |
CN108491615A (zh) * | 2018-03-17 | 2018-09-04 | 中国海洋大学 | 一种三层顶张式立管动力响应的有限元分析方法 |
CN108491615B (zh) * | 2018-03-17 | 2022-09-02 | 中国海洋大学 | 一种三层顶张式立管动力响应的有限元分析方法 |
CN109165447A (zh) * | 2018-08-27 | 2019-01-08 | 武汉理工大学 | 一种空间钢管结构壁厚损失的评价方法及*** |
CN109165447B (zh) * | 2018-08-27 | 2022-12-13 | 武汉理工大学 | 一种空间钢管结构壁厚损失的评价方法及*** |
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 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
CN114290644A (zh) * | 2021-12-17 | 2022-04-08 | 烟台市佶楷矿山设备科技有限公司 | 一种非金属长托辊的制造方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106909710B (zh) | 2020-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106909710A (zh) | 深水顶张式立管全耦合动力分析方法 | |
CN101539477B (zh) | 一种深水顶张式立管涡激振动与疲劳分析的方法 | |
CN1118618A (zh) | 挠性/刚性上升管*** | |
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 | |
Aubault et al. | WindFloat: A floating foundation for offshore wind turbines—Part III: Structural analysis | |
Huang et al. | Numerical simulation of deformations and forces of a floating fish cage collar in waves | |
CN101713707B (zh) | 深海钢悬链线立管波浪响应分析方法 | |
CN105243275A (zh) | 一种海底多跨管道涡激振动的预报方法 | |
CN102445318B (zh) | 一种顶张式立管顺流向振动分析方法 | |
Liu et al. | Nonlinear vibration model and response characteristic of drilling risers in deep-sea under soft suspension evacuation condition | |
CN102252897B (zh) | 一种深水顶张式立管弯曲振动分析方法 | |
CN108491615B (zh) | 一种三层顶张式立管动力响应的有限元分析方法 | |
CN107944100B (zh) | 综合生产脐带缆的线型设计方法 | |
Huang et al. | Application and analysis of elastomeric supporting style for topside modules of a spread mooring FPSO | |
Chandrasekaran et al. | Analysis and design of offshore triceratops under ultra-deep waters | |
Song et al. | Design of the World's deepest hybrid riser system for the cascade & chinook development | |
CN102289536B (zh) | 一种深水钻井立管参数横力弯曲分析方法 | |
CN102507115B (zh) | 一种深水顶张式立管大变形弯曲振动的分析方法 | |
Song et al. | Engineering design of deepwater free standing hybrid riser | |
CN103115198A (zh) | 一种深水海底增强管线 | |
Shi et al. | Coupled analysis approach in OTEC system design | |
Buberg | Design and analysis of steel catenary riser systems for deep waters | |
Damiani | JacketSE: An Offshore Wind Turbine Jacket Sizing Tool; Theory Manual and Sample Usage with Preliminary Validation | |
Orimolade | Steel lazy wave risers from turret moored FPSO | |
Lei et al. | Calculation method of residual fatigue life of in-service unbonded flexible pipe |
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 |