CN103353911B - 模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法 - Google Patents

模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法 Download PDF

Info

Publication number
CN103353911B
CN103353911B CN201310264437.1A CN201310264437A CN103353911B CN 103353911 B CN103353911 B CN 103353911B CN 201310264437 A CN201310264437 A CN 201310264437A CN 103353911 B CN103353911 B CN 103353911B
Authority
CN
China
Prior art keywords
subsea pipeline
dynamic
buckling
finite element
pipeline
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
Application number
CN201310264437.1A
Other languages
English (en)
Other versions
CN103353911A (zh
Inventor
刘润
熊昊
郭绍曾
李彪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201310264437.1A priority Critical patent/CN103353911B/zh
Publication of CN103353911A publication Critical patent/CN103353911A/zh
Application granted granted Critical
Publication of CN103353911B publication Critical patent/CN103353911B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Stored Programmes (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法：根据海底管线的几何和物理力学参数应用有限元软件ABAQUS，计算海底管线特征屈曲模态；根据海底管线与海底土体的几何参数和物理力学应用有限元软件ABAQUS，重新建立用于动态屈曲分析的模型，并引入海底管线的初始缺陷；显示动力法分析海底管线的动态整体屈曲过程。本发明可以有效的模拟细长的海底管线结构；并引入完全光滑的初始缺陷，避免应力集中，同时通过叠加引入海底管线的诸多模态可达到模拟真实高阶模态初始缺陷的目标；而且可以模拟海底管线在高温高压下的动态变化过程，对海底管线整体屈曲计算具有很好的收敛性，使得计算结果更能反映真实情况。

Description

模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法
技术领域
本发明涉及一种海底管线整体屈曲二维动力模拟方法。特别是涉及一种基于模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法。
背景技术
为了满足生产工艺的要求，海底油气海底管线内部通常具有较高的输送压力和温度。而海底管线通常为无缝的钢质管线，在高压和高温的作用下管壁内部出现应力的累积，由于海底管线受到地基土的约束作用而无法完全自由变形释放应力，当累积的应力达到某一临界状态时，海底管线将发生整体屈曲。这种随机且不可控的整体屈曲是关系到海底管线运营安全的重大技术问题，也是海底管线设计中必须考虑的关键节点。
通常海底管线在制造和铺设过程中，会具有初始缺陷，而具有初始缺陷的海底管线在高温高压联合作用下，更易发生整体屈曲。较大的整体屈曲一方面可能导致海底管线中的弯曲应力增长，接近或达到钢材的屈服强度，对海底管线的安全运营造成威胁；另一方面发生变形的海底管线容易受到渔业活动和船只航行的影响，增加安全隐患；此外，屈曲变形可能导致海底管线配重层、保温层结构遭到破坏甚至进水，影响海底管线的正常使用。
开展海底管线在温压联合作用下的全尺寸室内试验和现场试验均有较大的难度，因此研发合适的数值模拟方法具有重要的意义。这一问题的难度体现在以下几个方面：一是研究的对象具有明显尺度特征即海底管线沿轴线方向的尺度远大于其截面尺度；二是合理的初始缺陷的引入使建立的几何模型具有反映海底管线真实情况的特征；三是采用可以捕捉海底管线在温压联合作用下变形动态发展的分析方法。
发明内容
本发明所要解决的技术问题是，提供一种能够实现对具有初始缺陷海底管线在高温高压下发生整体屈曲动态过程数值模拟的模态法引入初始缺陷的海底管线整体屈曲二维动力模拟方法。
本发明所采用的技术方案是：一种模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，包括如下阶段：
1）根据海底管线的几何和物理力学参数应用有限元软件ABAQUS，计算海底管线特征屈曲模态；
2）根据海底管线与海底土体的几何参数和物理力学应用有限元软件ABAQUS，重新建立用于动态屈曲分析的模型，并引入海底管线的初始缺陷；
3）显示动力法分析海底管线的动态整体屈曲过程，包括如下步骤：
（1）在ABAQUS软件的INTERACTION模块中建立海底管线与海底土体表面的接触关系；
（2）在ABAQUS软件的LOAD模块中施加温度荷载与内压荷载；
（3）采用显示动力分析步Dynamic-Explicit模拟海底管线的动态屈曲整体过程。
阶段1）所述的计算海底管线特征屈曲模态，包括如下步骤：
（1）根据海底管线的几何参数应用有限元软件ABAQUS的PART模块建立海底管线的二维梁单元模型；
（2）根据海底管线的物理力学参数应用有限元软件ABAQUS的MATERIAL模块赋予海底管线的材料属性；
（3）应用有限元软件ABAQUS的MESH模块划分步骤1）得到的二维梁单元模型的网格；
（4）应用有限元软件ABAQUS的MATERIAL模块建立步骤1）得到的二维梁单元模型的边界条件；
（5）应用有限元软件ABAQUS的LOAD模块对步骤1）得到的二维梁单元模型施加温度荷载；
（6）设置有限元软件ABAQUS的STEP模块中的分析步类型为Buckle，通过Buckle计算海底管线1-20阶特征屈曲模态；
（7）应用有限元软件ABAQUS的关键词编辑器在*Restart语句后面添加*Nodefile语句，计算完成后，检查输出文件是否有海底管线模型的[输出结果文件名].Fil文件；
阶段2）包括如下步骤：
（1）应用有限元软件ABAQUS的PART模块，根据海底管线的几何参数建立海底管线的二维梁单元模型，根据海底土体的几何参数建立海底土体的二维平面应变单元模型；
（2）应用有限元软件ABAQUS的MATERIAL模块，根据海底管线的物理力学参数赋予海底管线的材料属性，根据海底土体的物理力学参数赋予海底土体的材料属性；
（3）应用有限元软件ABAQUS的MESH模块划分步骤1）得到的海底管线的二维梁单元模型和海底土体的二维平面应变单元模型的网格，新建的海底管线的二维梁单元模型网格必须与阶段1）中的二维梁单元模型的网格完全一致；
（4）应用有限元软件ABAQUS的MATERIAL模块建立步骤1）得到的海底管线的二维梁单元模型和海底土体的二维平面应变单元模型的边界条件；
（5）应用有限元软件ABAQUS的关键词编辑器在*Step语句中添加*Imperfection语句。
所述的在有限元软件ABAQUS的关键词编辑器中*Step语句中添加*Imperfection语句的格式如下：
*Imperfection，file=（Fil文件名），step=（Buckle分析步名）
1，μ1
2，μ2
3，μ3
……
n，μn
其中n是引入模态的阶数，μn是引入模态阶数对应的比例因子，μn的取值在0-1之间，n的取值建议在1～10之内。
阶段3）所述的在ABAQUS软件的INTERACTION模块中建立海底管线与海底土体表面的接触关系包括：法向接触行为选择硬接触，切向接触行为选择罚函数，输入海底管线与土体的摩擦系数为0.1～1.5。
阶段3）所述的温度荷载选择20℃～300℃。
阶段3）所述的内压荷载选择0MPa～200MPa。
阶段3）所述的模拟海底管线的动态屈曲整体过程具体是：设置有限元软件ABAQUS的STEP模块中的分析步类型为Dynamic-Explicit，通过Dynamic-Explicit计算模拟海底管线的动态屈曲整体过程。
本发明的模态法引入初始缺陷的海底管线整体屈曲二维动力模拟方法，是基于概率理论的模态方法引入海底管线的初始缺陷，进而开展二维模型的显示动力数值分析，该方法可以有效的模拟细长的海底管线结构；并引入完全光滑的初始缺陷，避免应力集中，同时通过叠加引入海底管线的诸多模态可达到模拟真实高阶模态初始缺陷的目标；而且可以模拟海底管线在高温高压下的动态变化过程，对海底管线整体屈曲计算具有很好的收敛性，使得计算结果更能反映真实情况。
附图说明
图1是本发明方法的流程图；
图2是海底管线屈曲模态分析模型图；
图3是海底管线整体屈曲分析模型图；
图4是不同温度荷载情况下海底管线沿轴向各点屈曲幅值变化；
图5是不同温度荷载情况下海底管线沿轴向各点轴向应力。
图中
a：二维梁单元海底管线部件b：二维平面应变单元海底土体部件
具体实施方式
下面结合实施例和附图对本发明的模态法引入初始缺陷的海底管线整体屈曲二维动力模拟方法做出详细说明。
本发明的模态法引入初始缺陷的海底管线整体屈曲二维动力模拟方法，提出了一种基于概率理论的海底管线初始缺陷的模态引入方法，以此为基础建立了二维梁单元显示动力模拟温压联合作用下的海底管线整体屈曲数值方法。本发明可以有效模拟海底管线的细长结构特点，可以较真实的还原海底管线初始缺陷形状，并可以实现对海底管线整体屈曲变形的动态模拟。
基于概率理论的模态分析法引入海底管线的初始缺陷，利用大型通用有限元软件ABAQUS的二维显示动力方法分析模拟海底管线在高温高压下的动态整体屈曲变形过程，本发明的模态法引入初始缺陷的海底管线整体屈曲二维动力模拟方法包括如下三个阶段。
1）根据海底管线的几何和物理力学参数应用有限元软件ABAQUS，计算海底管线特征屈曲模态，包括如下步骤：
（1）根据海底管线的几何参数应用有限元软件ABAQUS的PART模块建立海底管线的二维梁单元模型；
（2）根据海底管线的物理力学参数应用有限元软件ABAQUS的MATERIAL模块赋予海底管线的材料属性；
（3）应用有限元软件ABAQUS的MESH模块划分步骤1）得到的二维梁单元模型的网格；
（4）应用有限元软件ABAQUS的MATERIAL模块建立步骤1）得到的二维梁单元模型的边界条件；
（5）应用有限元软件ABAQUS的LOAD模块对步骤1）得到的二维梁单元模型施加温度荷载；
（6）设置有限元软件ABAQUS的STEP模块中的分析步类型为Buckle（屈曲特征值求解器），通过Buckle计算海底管线1-20阶特征屈曲模态；
（7）应用有限元软件ABAQUS的关键词编辑器在*Restart语句后面添加*Nodefile语句，计算完成后，检查输出文件是否有海底管线模型的[输出结果文件名].Fil文件；
2）根据海底管线与海底土体的几何参数和物理力学应用有限元软件ABAQUS，重新建立用于动态屈曲分析的模型，并引入海底管线的初始缺陷，包括如下步骤：
（1）应用有限元软件ABAQUS的PART模块，根据海底管线的几何参数建立海底管线的二维梁单元模型，根据海底土体的几何参数建立海底土体的二维平面应变单元模型；
（2）应用有限元软件ABAQUS的MATERIAL模块，根据海底管线的物理力学参数赋予海底管线的材料属性，根据海底土体的物理力学参数赋予海底土体的材料属性；
（3）应用有限元软件ABAQUS的MESH模块划分步骤1）得到的海底管线的二维梁单元模型和海底土体的二维平面应变单元模型的网格，新建的海底管线的二维梁单元模型网格必须与阶段1）中的二维梁单元模型的网格完全一致；
（4）应用有限元软件ABAQUS的MATERIAL模块建立步骤1）得到的海底管线的二维梁单元模型和海底土体的二维平面应变单元模型的边界条件；
（5）应用有限元软件ABAQUS的关键词编辑器在*Step语句中添加*Imperfection语句，
格式如下：
*Imperfection，file=（Fil文件名），step=（Buckle分析步名）
1，μ1
2，μ2
3，μ3
……
n，μn
其中n是引入模态的阶数，μn是引入模态阶数对应的比例因子，μn的取值在0-1之间，n的取值建议在1-10之内；
3）显示动力法分析海底管线的动态整体屈曲过程，包括如下步骤：
（1）在ABAQUS软件的INTERACTION模块中建立海底管线与海底土体表面的接触关系，其中法向接触行为建议选择硬接触，切向接触行为选择罚函数，输入海底管线与海底土体的摩擦系数为0.1～1.5；
（2）在ABAQUS软件的LOAD模块中施加温度荷载与内压荷载，温度荷载选择20℃～300℃，内压荷载选择0MPa-200MPa；
（3）采用显示动力分析步Dynamic-Explicit模拟海底管线的动态屈曲整体过程
设置有限元软件ABAQUS的STEP模块中的分析步类型为Dynamic-Explicit（显示动力求解器），通过Dynamic-Explicit求解器计算模拟海底管线的动态屈曲整体过程。
下面给出一个最佳实施例
某工程位于我国渤海，采用钢管运输油气，温度荷载为100℃、内压为80MPa海底管线与海底土体工程资料如表1与表2：
表1海底土体几何和物理力学参数
表2海底管线几何和物理力学参数
（1）进行海底管线屈曲模态分析
首先建立模型Model-1，创建二维梁单元海底管线部件Part-1，如图1所示，总长250m，截面尺寸与材料属性见表2；建立Buckle分析步，选择子空间求解器；设置温度荷载100℃，内压荷载80MPa；在*Restart语句后添加*Nodefile语句；提交计算，最后检查工作目录下是否有Job-1.fil文件生成。
（2）显示动力法分析海底管线整体屈曲
新建模型Model-2，创建二维梁单元海底管线部件Part-1与二维平面应变单元海床土体部件Part-2，海底管线部件必须与Model-1一致，如图2所示。海床部件取600m×3m，截面尺寸与材料属性见表1与表2；创建Part-1与Part-2的接触关系，摩擦系数取0.4；创建Step-1用于动态整体屈曲分析，选取Dynamic-Explicit求解器，设置海底管线两端为固定约束并施加100℃的温差与80MPa的内压；在*Step语句前添加*Imperfection语句：
*Imperfection，file=Job-1，step=1
1，0.15
最后提交Model-2进行计算分析。
（3）结果分析
图3为不同温度荷载情况下海底管线沿轴向各点屈曲幅值变化情况，由图可得海底管线随温度升高的整个屈曲破坏过程；图4为不同温度荷载情况下海底管线沿轴向各点的轴向应力，由图可以看出海底管线各点轴向应力随着温度升高的整个变化过程。

Claims (8)

1.一种模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，包括如下阶段：
1)根据海底管线的几何和物理力学参数应用有限元软件ABAQUS，计算海底管线特征屈曲模态；
2)根据海底管线与海底土体的几何参数和物理力学应用有限元软件ABAQUS，重新建立用于动态屈曲分析的模型，并引入海底管线的初始缺陷；
3)显示动力法分析海底管线的动态整体屈曲过程，包括如下步骤：
(1)在ABAQUS软件的INTERACTION模块中建立海底管线与海底土体表面的接触关系；
(2)在ABAQUS软件的LOAD模块中施加温度荷载与内压荷载；
(3)采用显示动力分析步Dynamic-Explicit模拟海底管线的动态屈曲整体过程。
2.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段1)所述的计算海底管线特征屈曲模态，包括如下步骤：
(1)根据海底管线的几何参数应用有限元软件ABAQUS的PART模块建立海底管线的二维梁单元模型；
(2)根据海底管线的物理力学参数应用有限元软件ABAQUS的MATERIAL模块赋予海底管线的材料属性；
(3)应用有限元软件ABAQUS的MESH模块划分步骤1)得到的二维梁单元模型的网格；
(4)应用有限元软件ABAQUS的MATERIAL模块建立步骤1)得到的二维梁单元模型的边界条件；
(5)应用有限元软件ABAQUS的LOAD模块对步骤1)得到的二维梁单元模型施加温度荷载；
(6)设置有限元软件ABAQUS的STEP模块中的分析步类型为Buckle，通过Buckle计算海底管线1-20阶特征屈曲模态；
(7)应用有限元软件ABAQUS的关键词编辑器在*Restart语句后面添加*Nodefile语句，计算完成后，检查输出文件是否有海底管线模型的[输出结果文件名].Fil文件。
3.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段2)包括如下步骤：
(1)应用有限元软件ABAQUS的PART模块，根据海底管线的几何参数建立海底管线的二维梁单元模型，根据海底土体的几何参数建立海底土体的二维平面应变单元模型；
(2)应用有限元软件ABAQUS的MATERIAL模块，根据海底管线的物理力学参数赋予海底管线的材料属性，根据海底土体的物理力学参数赋予海底土体的材料属性；
(3)应用有限元软件ABAQUS的MESH模块划分步骤1)得到的海底管线的二维梁单元模型和海底土体的二维平面应变单元模型的网格，新建的海底管线的二维梁单元模型网格必须与阶段1)中的二维梁单元模型的网格完全一致；
(4)应用有限元软件ABAQUS的MATERIAL模块建立步骤1)得到的海底管线的二维梁单元模型和海底土体的二维平面应变单元模型的边界条件；
(5)应用有限元软件ABAQUS的关键词编辑器在*Step语句中添加*Imperfection语句。
4.根据权利要求3所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，所述的在有限元软件ABAQUS的关键词编辑器中*Step语句中添加*Imperfection语句的格式如下：
*Imperfection，file＝(Fil文件名)，step＝(Buckle分析步名)
1，μ1
2，μ2
3，μ3
……
n，μn
其中n是引入模态的阶数，μn是引入模态阶数对应的比例因子，μn的取值在0-1之间，n的取值在1～10之内。
5.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段3)所述的在ABAQUS软件的INTERACTION模块中建立海底管线与海底土体表面的接触关系包括：法向接触行为选择硬接触，切向接触行为选择罚函数，输入海底管线与土体的摩擦系数为0.1～1.5。
6.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段3)所述的温度荷载选择20℃～300℃。
7.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段3)所述的内压荷载选择0MPa～200MPa。
8.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段3)所述的模拟海底管线的动态屈曲整体过程具体是：设置有限元软件ABAQUS的STEP模块中的分析步类型为Dynamic-Explicit，通过Dynamic-Explicit计算模拟海底管线的动态屈曲整体过程。
CN201310264437.1A 2013-06-27 2013-06-27 模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法 Expired - Fee Related CN103353911B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310264437.1A CN103353911B (zh) 2013-06-27 2013-06-27 模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310264437.1A CN103353911B (zh) 2013-06-27 2013-06-27 模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法

Publications (2)

Publication Number Publication Date
CN103353911A CN103353911A (zh) 2013-10-16
CN103353911B true CN103353911B (zh) 2016-01-13

Family

ID=49310282

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310264437.1A Expired - Fee Related CN103353911B (zh) 2013-06-27 2013-06-27 模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法

Country Status (1)

Country Link
CN (1) CN103353911B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103714206A (zh) * 2013-12-17 2014-04-09 天津大学 基于模态引入初始缺陷的海底管线二维静态数值模拟方法
CN103714207A (zh) * 2013-12-17 2014-04-09 天津大学 基于模态引入初始缺陷的海底管线三维动力数值模拟方法
CN103729504A (zh) * 2013-12-17 2014-04-16 天津大学 基于模态引入初始缺陷的海底管线三维静态数值模拟方法
CN106934160B (zh) * 2017-03-14 2020-04-21 天津大学 海底管道整体屈曲数值模拟中管道长度的确定方法
CN107291982A (zh) * 2017-05-18 2017-10-24 天津大学 PSI与改进Riks法耦合的管线水平向整体屈曲模拟方法
CN107941538B (zh) * 2017-10-27 2019-11-01 天津大学 一种海底管线铺设安装模拟测试方法
CN110765505B (zh) * 2019-09-17 2023-06-30 中国石油天然气集团有限公司 一种含表面划伤复合凹陷油气管道极限内压的预测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101561839A (zh) * 2008-05-13 2009-10-21 赵汝江 海床上石油管道悬跨安全性的分析测算方法
CN101702186A (zh) * 2009-11-17 2010-05-05 中冶南方工程技术有限公司 高炉炉壳结构设计方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101561839A (zh) * 2008-05-13 2009-10-21 赵汝江 海床上石油管道悬跨安全性的分析测算方法
CN101702186A (zh) * 2009-11-17 2010-05-05 中冶南方工程技术有限公司 高炉炉壳结构设计方法

Also Published As

Publication number Publication date
CN103353911A (zh) 2013-10-16

Similar Documents

Publication Publication Date Title
CN103353911B (zh) 模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法
Lopatin et al. Buckling of the composite sandwich cylindrical shell with clamped ends under uniform external pressure
Hong et al. A lateral global buckling failure envelope for a high temperature and high pressure (HT/HP) submarine pipeline
Liu et al. Numerical studies on global buckling of subsea pipelines
Jiang et al. Probability-interval hybrid reliability analysis for cracked structures existing epistemic uncertainty
Shi et al. Buckling resistance of grid-stiffened carbon-fiber thin-shell structures
Viola et al. GDQFEM numerical simulations of continuous media with cracks and discontinuities
Mirzaei et al. Experimental and numerical analysis of dynamic rupture of steel pipes under internal high-speed moving pressures
Shariyat et al. Three-dimensional non-linear elasticity-based 3D cubic B-spline finite element shear buckling analysis of rectangular orthotropic FGM plates surrounded by elastic foundations
Xia et al. Modified interval perturbation finite element method for a structural-acoustic system with interval parameters
Liang et al. Constructing simplified models for dynamic analysis of monopile-supported offshore wind turbines
Cuamatzi-Melendez et al. Finite element and theoretical analyses of bisymmetric collapses in flexible risers for deepwaters developments
Gu et al. OpenSees-SNOPT framework for finite-element-based optimization of structural and geotechnical systems
Huang et al. Numerical simulation of deformations and forces of a floating fish cage collar in waves
CN103714206A (zh) 基于模态引入初始缺陷的海底管线二维静态数值模拟方法
DK201170431A (en) Finite element modelling of borehole seismic responses in layered anisotropic formations and its use for elastic inversion
Petrolo Flutter analysis of composite lifting surfaces by the 1D Carrera Unified Formulation and the doublet lattice method
CN103714207A (zh) 基于模态引入初始缺陷的海底管线三维动力数值模拟方法
Yang et al. Multi-objective shape optimization design for liquefied natural gas cryogenic helical corrugated steel pipe
Han et al. Transient response analysis by model order reduction of a Mokpo-Jeju submerged floating tunnel under seismic excitations
Ebrahimi et al. Finite element investigation on the tensile armor wire response of flexible pipe for axisymmetric loading conditions using an implicit solver
Vasilikis et al. Buckling design of confined steel cylinders under external pressure
CN103729504A (zh) 基于模态引入初始缺陷的海底管线三维静态数值模拟方法
Li et al. A strain energy-based equivalent layer method for the prediction of critical collapse pressure of flexible risers
Zhang et al. Analysis of torsional instability and reliability of marine flexible pipelines

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160113

Termination date: 20200627