管道受坠物冲击载荷作用的数值模拟基本方法

管道受坠物冲击载荷作用的数值模拟基本方法

【摘要】管道受冲击载荷作用是复杂的非线性接触问题，本文介绍了可以进行非线性分析的有限元软件、离散元软件的概况和基本原理，以及运用它们对冲击载荷作用引起的管道动力响应过程进行数值模拟的基本方法。

【关键词】管道冲击非线性数值模拟

随着我国油气管道建设的进一步深入，石油天然气管道穿越复杂地质条件的工程实践越来越多，这些管道多沿山体坡脚埋设，可能要经过滑坡、泥石流等自然灾害高发地段，由于自然灾害所产生的高速坠落的石块容易冲击管道导致管道失效。同时，在日益发展的海洋石油开采中，海底石油管道也容易在其安装与油气输送过程中，与锚泊作业以及货物运输等人类活动造成坠的落物体发生碰撞，造成管道损伤。因此，对管道受坠物冲击作用引起管线变形的规律和破坏机理进行深入研究具有重要意义。本文将以有限元方法为基础，介绍管道受坠物冲击载荷作用的数值模拟基本方法。

管道受坠物冲击载荷作用是管道-土体组成的体系在冲击荷载下的整体动力响应。无论是从静力学还是动力学的角度来分析结构的受力状态，管道与土体的相互作用都是不可忽略的，只有把管道与地基作为相互作用又相互制约的整体分析，才能得到比较符合实际的计算结果。随着数值非线性分析成为解决岩土工程问题的重要手段，有限元、离散元等方法在管土相互作用分析中也发挥着越来越大的作用，基于这些理论的数值模拟软件也得到了极大的发展。

1 非线性数值模拟软件

1.1 ANSYS/ABAQUS

ANSYS是一种大型通用有限元分析软件，融结构、流体、电场、磁场、声场分析于一体，由世界上最大的有限元分析软件公司之一的美国ANSYS开发，它能与多数CAD软件接口，实现数据的共享和交换，是现代产品设计中的高级CAE工具之一。有限元法（FEA，Finite Element Analysis）的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定合适的形函数，然后附上求解这个域总的满足条件，如结构的平衡条件、边界条件等，从而得到问题的解。这个解不是准确解，而是近似解，随着形函数精度的提高，有限元方法可以得到相对很高的计算精度高，而且能适应各种复杂形状，这样实际问题被较简单的问题所代替，有限元成为行之有效的工程分析手段。ANSYS主要分析类型包括：结构静力分析，结构动力学分析，结构非线性分析等。

ABAQUS同样也是一款功能强大的，以有限元理论为基础的工程模拟软件。与ANSYS相比，ABAQUS软件在求解非线性问题时具有非常明显的优势，其非线性涵盖材料非线性、几何非线性和状态非线性等多个方面，而且采用了人机

交互界面，其强大的非线性分析功能已经在设计和研究的用户中得到了广泛的认可。

1.2 3DEC

3DEC一款基于离散单元法基本理论来描述离散介质力学行为的数值分析软件，是英文3 Dimension Distinct Element Code的缩写。离散单元法是Cundall 于1971年提出来的，离散单元法也像有限单元法那样，将区域划分成单元。但是单元因受节理等不连续面控制，在以后的运动过程中，单元节点可以分离，即一个单元与其邻近单元可以接触，也可以分开。单元之间相互作用的力可以根据力和位移的关系求出，而个别单元的运动则完全根据单元所受的不平衡力和不平衡力矩的大小按牛顿运动定律确定。离散单元法是继有限元法之后，用于分析物质系统动力学问题的又一种强有力的数值计算方法，3DEC为解决众多涉及颗粒、结构、流体与电磁及其耦合等综合问题提供了一个平台，已成为过程分析、设计优化和产品研发的一种强有力的工具。它拥有极其丰富的岩土体或结构材料模型库，如常用的：弹性模型、各向异性模型、莫尔库仑模型、Drucker-Prager 模型、双线性塑性模型、应变软化模型、蠕变模型和用户自定义模型。

2 有限元模型建立

图1？裸置管道、埋地管道模型示意图

2.1 非线性接触求解过程

管道受冲击载荷作用是复杂的非线性接触问题，必须要考虑管道与地基的相互作用。接触问题是一类非线性问题，属于边界条件非线性问题，有别于材料非线性和几何非线性。接触问题中的边界条件不是在计算开始就给出，而是在计算过程中产生。两接触体之间接触面的面积与压力分布随外载荷变化而变化，而且与接触体刚性有关。一般采用罚函数法解决接触问题，罚函数法[2]是一种近似方法，实质是将接触非线性问题转化为材料非线性问题。由于计算简单，与显式算法完全兼容，所以使用比较广泛。

根据管道与土之间的相互作用特点，将接触面约束条件引入惩罚势能泛函，接触问题就等价于无约束优化问题，则相应泛函的极值条件为：

其中：M、C、K、pK、F 分别为系统的质量矩阵、阻尼矩阵、整体切向刚度矩阵、接触刚度矩阵、节点总体外载荷向量。

非线性分析软件ABAQUS认为，当两个表面直接间隙为零时，两个表面接触，并在相应节点上施加接触约束，对两个接触表面直接能够传递的接触压力的大小没有任何限制。而当接触面之间的接触压力变为零或负值时，两个接触面分离开来，同时解除相应节点上的接触约束，两表面可以发生任意的相对滑动。对于裸置管道，如图1中裸置管道模型，可以采用Explicit设置两个接触对：使用General contact定义坠物表面与管道上表面的接触，用Surface-to -surface contact

定义管道与土体的接触。而对于埋地管道，如图2中埋地管道模型，同样可以采用Explicit设置两个接触对：使用General contact定义坠物表面与土体表面的接触，用Surface-to-surface contact定义管道与土体的接触。

2.2 土体本构模型

f 为动剪应力的函数。