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

摘要
申请专利号：	CN201310264437.1	申请日：	2013.06.27
公开号：	CN103353911A	公开日：	2013.10.16
当前法律状态：	授权	有效性：	有权
法律详情：	授权\|\|\|实质审查的生效IPC(主分类):G06F 17/50申请日:20130627\|\|\|公开
IPC分类号：	G06F17/50	主分类号：	G06F17/50
申请人：	天津大学
发明人：	刘润; 熊昊; 郭绍曾; 李彪
地址：	300072 天津市南开区卫津路92号
优先权：
专利代理机构：	天津市北洋有限责任专利代理事务所 12201	代理人：	杜文茹
PDF下载：	PDF下载

内容摘要

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

权利要求书

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，  μ₁

2，  μ₂

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计算模拟海底管线的动态屈曲整体过程。

说明书

模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法
技术领域
本发明涉及一种海底管线整体屈曲二维动力模拟方法。特别是涉及一种基于模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法。
背景技术
为了满足生产工艺的要求，海底油气海底管线内部通常具有较高的输送压力和温度。而海底管线通常为无缝的钢质管线，在高压和高温的作用下管壁内部出现应力的累积，由于海底管线受到地基土的约束作用而无法完全自由变形释放应力，当累积的应力达到某一临界状态时，海底管线将发生整体屈曲。这种随机且不可控的整体屈曲是关系到海底管线运营安全的重大技术问题，也是海底管线设计中必须考虑的关键节点。
通常海底管线在制造和铺设过程中，会具有初始缺陷，而具有初始缺陷的海底管线在高温高压联合作用下，更易发生整体屈曲。较大的整体屈曲一方面可能导致海底管线中的弯曲应力增长，接近或达到钢材的屈服强度，对海底管线的安全运营造成威胁；另一方面发生变形的海底管线容易受到渔业活动和船只航行的影响，增加安全隐患；此外，屈曲变形可能导致海底管线配重层、保温层结构遭到破坏甚至进水，影响海底管线的正常使用。
开展海底管线在温压联合作用下的全尺寸室内试验和现场试验均有较大的难度，因此研发合适的数值模拟方法具有重要的意义。这一问题的难度体现在以下几个方面：一是研究的对象具有明显尺度特征即海底管线沿轴线方向的尺度远大于其截面尺度；二是合理的初始缺陷的引入使建立的几何模型具有反映海底管线真实情况的特征；三是采用可以捕捉海底管线在温压联合作用下变形动态发展的分析方法。
发明内容
本发明所要解决的技术问题是，提供一种能够实现对具有初始缺陷海底管线在高温高压下发生整体屈曲动态过程数值模拟的模态法引入初始缺陷的海底管线整体屈曲二维动力模拟方法。
本发明所采用的技术方案是：一种模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，包括如下阶段：
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，μ₁
2，μ₂
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，μ₁
2，μ₂
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为不同温度荷载情况下海底管线沿轴向各点的轴向应力，由图可以看出海底管线各点轴向应力随着温度升高的整个变化过程。

资源描述

《模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法.pdf》由会员分享，可在线阅读，更多相关《模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法.pdf（11页珍藏版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 103353911 A(43)申请公布日 2013.10.16CN103353911A*CN103353911A*(21)申请号 201310264437.1(22)申请日 2013.06.27G06F 17/50(2006.01)(71)申请人天津大学地址 300072 天津市南开区卫津路92号(72)发明人刘润熊昊郭绍曾李彪(74)专利代理机构天津市北洋有限责任专利代理事务所 12201代理人杜文茹(54) 发明名称模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法(57) 摘要一种模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法：根据海底管线的几何和物理。

2、力学参数应用有限元软件ABAQUS，计算海底管线特征屈曲模态；根据海底管线与海底土体的几何参数和物理力学应用有限元软件ABAQUS，重新建立用于动态屈曲分析的模型，并引入海底管线的初始缺陷；显示动力法分析海底管线的动态整体屈曲过程。本发明可以有效的模拟细长的海底管线结构；并引入完全光滑的初始缺陷，避免应力集中，同时通过叠加引入海底管线的诸多模态可达到模拟真实高阶模态初始缺陷的目标；而且可以模拟海底管线在高温高压下的动态变化过程，对海底管线整体屈曲计算具有很好的收敛性，使得计算结果更能反映真实情况。(51)Int.Cl.权利要求书2页说明书5页附图3页(19)中华人民共和国国家知识产权局(1。

3、2)发明专利申请权利要求书2页说明书5页附图3页(10)申请公布号 CN 103353911 ACN 103353911 A1/2页21.一种模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，包括如下阶段：1）根据海底管线的几何和物理力学参数应用有限元软件ABAQUS，计算海底管线特征屈曲模态；2）根据海底管线与海底土体的几何参数和物理力学应用有限元软件ABAQUS，重新建立用于动态屈曲分析的模型，并引入海底管线的初始缺陷；3）显示动力法分析海底管线的动态整体屈曲过程，包括如下步骤：（1）在ABAQUS软件的INTERACTION模块中建立海底管线与海底土体表面的接触关系；（。

4、2）在ABAQUS软件的LOAD模块中施加温度荷载与内压荷载；（3）采用显示动力分析步Dynamic-Explicit模拟海底管线的动态屈曲整体过程。2.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段1）所述的计算海底管线特征屈曲模态，包括如下步骤：（1）根据海底管线的几何参数应用有限元软件ABAQUS的PART模块建立海底管线的二维梁单元模型；（2）根据海底管线的物理力学参数应用有限元软件ABAQUS的MATERIAL模块赋予海底管线的材料属性；（3）应用有限元软件ABAQUS的MESH模块划分步骤1）得到的二维梁单元模型的网格；（4）应用有限元软件。

5、ABAQUS的MATERIAL模块建立步骤1）得到的二维梁单元模型的边界条件；（5）应用有限元软件ABAQUS的LOAD模块对步骤1）得到的二维梁单元模型施加温度荷载；（6）设置有限元软件ABAQUS的STEP模块中的分析步类型为Buckle，通过Buckle计算海底管线1-20阶特征屈曲模态；（7）应用有限元软件ABAQUS的关键词编辑器在*Restart语句后面添加*Nodefile语句，计算完成后，检查输出文件是否有海底管线模型的输出结果文件名.Fil文件。3.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段2）包括如下步骤：（1）应用有限元软件A。

6、BAQUS的PART模块，根据海底管线的几何参数建立海底管线的二维梁单元模型，根据海底土体的几何参数建立海底土体的二维平面应变单元模型；（2）应用有限元软件ABAQUS的MATERIAL模块，根据海底管线的物理力学参数赋予海底管线的材料属性，根据海底土体的物理力学参数赋予海底土体的材料属性；（3）应用有限元软件ABAQUS的MESH模块划分步骤1）得到的海底管线的二维梁单元模型和海底土体的二维平面应变单元模型的网格，新建的海底管线的二维梁单元模型网格必须与阶段1）中的二维梁单元模型的网格完全一致；（4）应用有限元软件ABAQUS的MATERIAL模块建立步骤1）得到的海底管线的二维梁单元模型和。

7、海底土体的二维平面应变单元模型的边界条件；（5）应用有限元软件ABAQUS的关键词编辑器在*Step语句中添加*Imperfection语句。4.根据权利要求3所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟权利要求书CN 103353911 A2/2页3方法，其特征在于，所述的在有限元软件ABAQUS的关键词编辑器中*Step语句中添加*Imperfection语句的格式如下：*Imperfection，file=（Fil文件名），step=（Buckle分析步名）1，12，23，3n，n其中n是引入模态的阶数，n是引入模态阶数对应的比例因子，n的取值在0-1之间，n的取值建议在。

8、110之内。5.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段3）所述的在ABAQUS软件的INTERACTION模块中建立海底管线与海底土体表面的接触关系包括：法向接触行为选择硬接触，切向接触行为选择罚函数，输入海底管线与土体的摩擦系数为0.11.5。6.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段3）所述的温度荷载选择20300。7.根据权利要求1所述的模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，其特征在于，阶段3）所述的内压荷载选择0MPa200MPa。8.根据权利要求1所述的模态法引入初始缺陷。

9、海底管线整体屈曲二维动力模拟方法，其特征在于，阶段3）所述的模拟海底管线的动态屈曲整体过程具体是：设置有限元软件ABAQUS的STEP模块中的分析步类型为Dynamic-Explicit，通过Dynamic-Explicit计算模拟海底管线的动态屈曲整体过程。权利要求书CN 103353911 A1/5页4模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法技术领域0001 本发明涉及一种海底管线整体屈曲二维动力模拟方法。特别是涉及一种基于模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法。背景技术0002 为了满足生产工艺的要求，海底油气海底管线内部通常具有较高的输送压力和温度。而海底。

10、管线通常为无缝的钢质管线，在高压和高温的作用下管壁内部出现应力的累积，由于海底管线受到地基土的约束作用而无法完全自由变形释放应力，当累积的应力达到某一临界状态时，海底管线将发生整体屈曲。这种随机且不可控的整体屈曲是关系到海底管线运营安全的重大技术问题，也是海底管线设计中必须考虑的关键节点。0003 通常海底管线在制造和铺设过程中，会具有初始缺陷，而具有初始缺陷的海底管线在高温高压联合作用下，更易发生整体屈曲。较大的整体屈曲一方面可能导致海底管线中的弯曲应力增长，接近或达到钢材的屈服强度，对海底管线的安全运营造成威胁；另一方面发生变形的海底管线容易受到渔业活动和船只航行的影响，增加安全隐患；此外。

11、，屈曲变形可能导致海底管线配重层、保温层结构遭到破坏甚至进水，影响海底管线的正常使用。0004 开展海底管线在温压联合作用下的全尺寸室内试验和现场试验均有较大的难度，因此研发合适的数值模拟方法具有重要的意义。这一问题的难度体现在以下几个方面：一是研究的对象具有明显尺度特征即海底管线沿轴线方向的尺度远大于其截面尺度；二是合理的初始缺陷的引入使建立的几何模型具有反映海底管线真实情况的特征；三是采用可以捕捉海底管线在温压联合作用下变形动态发展的分析方法。发明内容0005 本发明所要解决的技术问题是，提供一种能够实现对具有初始缺陷海底管线在高温高压下发生整体屈曲动态过程数值模拟的模态法引入初始缺陷的海。

12、底管线整体屈曲二维动力模拟方法。0006 本发明所采用的技术方案是：一种模态法引入初始缺陷海底管线整体屈曲二维动力模拟方法，包括如下阶段：0007 1）根据海底管线的几何和物理力学参数应用有限元软件ABAQUS，计算海底管线特征屈曲模态；0008 2）根据海底管线与海底土体的几何参数和物理力学应用有限元软件ABAQUS，重新建立用于动态屈曲分析的模型，并引入海底管线的初始缺陷；0009 3）显示动力法分析海底管线的动态整体屈曲过程，包括如下步骤：0010 （1）在ABAQUS软件的INTERACTION模块中建立海底管线与海底土体表面的接触关系；0011 （2）在ABAQUS软件的LOAD模块。

13、中施加温度荷载与内压荷载；0012 （3）采用显示动力分析步Dynamic-Explicit模拟海底管线的动态屈曲整体过程。说明书CN 103353911 A2/5页50013 阶段1）所述的计算海底管线特征屈曲模态，包括如下步骤：0014 （1）根据海底管线的几何参数应用有限元软件ABAQUS的PART模块建立海底管线的二维梁单元模型；0015 （2）根据海底管线的物理力学参数应用有限元软件ABAQUS的MATERIAL模块赋予海底管线的材料属性；0016 （3）应用有限元软件ABAQUS的MESH模块划分步骤1）得到的二维梁单元模型的网格；0017 （4）应用有限元软件ABAQUS的M。

14、ATERIAL模块建立步骤1）得到的二维梁单元模型的边界条件；0018 （5）应用有限元软件ABAQUS的LOAD模块对步骤1）得到的二维梁单元模型施加温度荷载；0019 （6）设置有限元软件ABAQUS的STEP模块中的分析步类型为Buckle，通过Buckle计算海底管线1-20阶特征屈曲模态；0020 （7）应用有限元软件ABAQUS的关键词编辑器在*Restart语句后面添加*Nodefile语句，计算完成后，检查输出文件是否有海底管线模型的输出结果文件名.Fil文件；0021 阶段2）包括如下步骤：0022 （1）应用有限元软件ABAQUS的PART模块，根据海底管线的几何参数建立海。

15、底管线的二维梁单元模型，根据海底土体的几何参数建立海底土体的二维平面应变单元模型；0023 （2）应用有限元软件ABAQUS的MATERIAL模块，根据海底管线的物理力学参数赋予海底管线的材料属性，根据海底土体的物理力学参数赋予海底土体的材料属性；0024 （3）应用有限元软件ABAQUS的MESH模块划分步骤1）得到的海底管线的二维梁单元模型和海底土体的二维平面应变单元模型的网格，新建的海底管线的二维梁单元模型网格必须与阶段1）中的二维梁单元模型的网格完全一致；0025 （4）应用有限元软件ABAQUS的MATERIAL模块建立步骤1）得到的海底管线的二维梁单元模型和海底土体的二维平面应变单。

16、元模型的边界条件；0026 （5）应用有限元软件ABAQUS的关键词编辑器在*Step语句中添加*Imperfection语句。0027 所述的在有限元软件ABAQUS的关键词编辑器中*Step语句中添加*Imperfection语句的格式如下：0028 *Imperfection，file=（Fil文件名），step=（Buckle分析步名）0029 1，10030 2，20031 3，30032 0033 n，n0034 其中n是引入模态的阶数，n是引入模态阶数对应的比例因子，n的取值在0-1之间，n的取值建议在110之内。0035 阶段3）所述的在ABAQUS软件的INTERACTION。

17、模块中建立海底管线与海底土体表面的接触关系包括：法向接触行为选择硬接触，切向接触行为选择罚函数，输入海底管线与说明书CN 103353911 A3/5页6土体的摩擦系数为0.11.5。0036 阶段3）所述的温度荷载选择20300。0037 阶段3）所述的内压荷载选择0MPa200MPa。0038 阶段3）所述的模拟海底管线的动态屈曲整体过程具体是：设置有限元软件ABAQUS的STEP模块中的分析步类型为Dynamic-Explicit，通过Dynamic-Explicit计算模拟海底管线的动态屈曲整体过程。0039 本发明的模态法引入初始缺陷的海底管线整体屈曲二维动力模拟方法，是基于概率。

18、理论的模态方法引入海底管线的初始缺陷，进而开展二维模型的显示动力数值分析，该方法可以有效的模拟细长的海底管线结构；并引入完全光滑的初始缺陷，避免应力集中，同时通过叠加引入海底管线的诸多模态可达到模拟真实高阶模态初始缺陷的目标；而且可以模拟海底管线在高温高压下的动态变化过程，对海底管线整体屈曲计算具有很好的收敛性，使得计算结果更能反映真实情况。附图说明0040 图1是本发明方法的流程图；0041 图2是海底管线屈曲模态分析模型图；0042 图3是海底管线整体屈曲分析模型图；0043 图4是不同温度荷载情况下海底管线沿轴向各点屈曲幅值变化；0044 图5是不同温度荷载情况下海底管线沿轴向各点轴向应。

19、力。0045 图中0046 a：二维梁单元海底管线部件 b：二维平面应变单元海底土体部件具体实施方式0047 下面结合实施例和附图对本发明的模态法引入初始缺陷的海底管线整体屈曲二维动力模拟方法做出详细说明。0048 本发明的模态法引入初始缺陷的海底管线整体屈曲二维动力模拟方法，提出了一种基于概率理论的海底管线初始缺陷的模态引入方法，以此为基础建立了二维梁单元显示动力模拟温压联合作用下的海底管线整体屈曲数值方法。本发明可以有效模拟海底管线的细长结构特点，可以较真实的还原海底管线初始缺陷形状，并可以实现对海底管线整体屈曲变形的动态模拟。0049 基于概率理论的模态分析法引入海底管线的初始缺陷，利用。

20、大型通用有限元软件ABAQUS的二维显示动力方法分析模拟海底管线在高温高压下的动态整体屈曲变形过程，本发明的模态法引入初始缺陷的海底管线整体屈曲二维动力模拟方法包括如下三个阶段。0050 1）根据海底管线的几何和物理力学参数应用有限元软件ABAQUS，计算海底管线特征屈曲模态，包括如下步骤：0051 （1）根据海底管线的几何参数应用有限元软件ABAQUS的PART模块建立海底管线的二维梁单元模型；0052 （2）根据海底管线的物理力学参数应用有限元软件ABAQUS的MATERIAL模块赋予海底管线的材料属性；说明书CN 103353911 A4/5页70053 （3）应用有限元软件ABAQ。

21、US的MESH模块划分步骤1）得到的二维梁单元模型的网格；0054 （4）应用有限元软件ABAQUS的MATERIAL模块建立步骤1）得到的二维梁单元模型的边界条件；0055 （5）应用有限元软件ABAQUS的LOAD模块对步骤1）得到的二维梁单元模型施加温度荷载；0056 （6）设置有限元软件ABAQUS的STEP模块中的分析步类型为Buckle（屈曲特征值求解器），通过Buckle计算海底管线1-20阶特征屈曲模态；0057 （7）应用有限元软件ABAQUS的关键词编辑器在*Restart语句后面添加*Nodefile语句，计算完成后，检查输出文件是否有海底管线模型的输出结果文件名.Fil。

22、文件；0058 2）根据海底管线与海底土体的几何参数和物理力学应用有限元软件ABAQUS，重新建立用于动态屈曲分析的模型，并引入海底管线的初始缺陷，包括如下步骤：0059 （1）应用有限元软件ABAQUS的PART模块，根据海底管线的几何参数建立海底管线的二维梁单元模型，根据海底土体的几何参数建立海底土体的二维平面应变单元模型；0060 （2）应用有限元软件ABAQUS的MATERIAL模块，根据海底管线的物理力学参数赋予海底管线的材料属性，根据海底土体的物理力学参数赋予海底土体的材料属性；0061 （3）应用有限元软件ABAQUS的MESH模块划分步骤1）得到的海底管线的二维梁单元模型和海底。

23、土体的二维平面应变单元模型的网格，新建的海底管线的二维梁单元模型网格必须与阶段1）中的二维梁单元模型的网格完全一致；0062 （4）应用有限元软件ABAQUS的MATERIAL模块建立步骤1）得到的海底管线的二维梁单元模型和海底土体的二维平面应变单元模型的边界条件；0063 （5）应用有限元软件ABAQUS的关键词编辑器在*Step语句中添加*Imperfection语句，0064 格式如下：0065 *Imperfection，file=（Fil文件名），step=（Buckle分析步名）0066 1，10067 2，20068 3，30069 0070 n，n0071 其中n是引入模态的阶。

24、数，n是引入模态阶数对应的比例因子，n的取值在0-1之间，n的取值建议在1-10之内；0072 3）显示动力法分析海底管线的动态整体屈曲过程，包括如下步骤：0073 （1）在ABAQUS软件的INTERACTION模块中建立海底管线与海底土体表面的接触关系，其中法向接触行为建议选择硬接触，切向接触行为选择罚函数，输入海底管线与海底土体的摩擦系数为0.11.5；0074 （2）在ABAQUS软件的LOAD模块中施加温度荷载与内压荷载，温度荷载选择20300，内压荷载选择0MPa-200MPa；0075 （3）采用显示动力分析步Dynamic-Explicit模拟海底管线的动态屈曲整体过程说明。

25、书CN 103353911 A5/5页80076 设置有限元软件ABAQUS的STEP模块中的分析步类型为Dynamic-Explicit（显示动力求解器），通过Dynamic-Explicit求解器计算模拟海底管线的动态屈曲整体过程。0077 下面给出一个最佳实施例0078 某工程位于我国渤海，采用钢管运输油气，温度荷载为100、内压为80MPa海底管线与海底土体工程资料如表1与表2：0079 表1海底土体几何和物理力学参数0080 0081 表2海底管线几何和物理力学参数0082 0083 （1）进行海底管线屈曲模态分析0084 首先建立模型Model-1，创建二维梁单元海底管线部件Par。

26、t-1，如图1所示，总长250m，截面尺寸与材料属性见表2；建立Buckle分析步，选择子空间求解器；设置温度荷载100，内压荷载80MPa；在*Restart语句后添加*Nodefile语句；提交计算，最后检查工作目录下是否有Job-1.fil文件生成。0085 （2）显示动力法分析海底管线整体屈曲0086 新建模型Model-2，创建二维梁单元海底管线部件Part-1与二维平面应变单元海床土体部件Part-2，海底管线部件必须与Model-1一致，如图2所示。海床部件取600m3m，截面尺寸与材料属性见表1与表2；创建Part-1与Part-2的接触关系，摩擦系数取0.4；创建Step-1。

27、用于动态整体屈曲分析，选取Dynamic-Explicit求解器，设置海底管线两端为固定约束并施加100的温差与80MPa的内压；在*Step语句前添加*Imperfection语句：0087 *Imperfection，file=Job-1，step=10088 1，0.150089 最后提交Model-2进行计算分析。0090 （3）结果分析0091 图3为不同温度荷载情况下海底管线沿轴向各点屈曲幅值变化情况，由图可得海底管线随温度升高的整个屈曲破坏过程；图4为不同温度荷载情况下海底管线沿轴向各点的轴向应力，由图可以看出海底管线各点轴向应力随着温度升高的整个变化过程。说明书CN 103353911 A1/3页9图1说明书附图CN 103353911 A2/3页10图2图3图4说明书附图CN 103353911 A10。

展开阅读全文