爆炸地震波作用下埋地管道流固耦合响应研究.docxVIP

爆炸地震波作用下埋地管道流固耦合响应研究 0 管道内爆振动的定义 埋地管道作为一种高效的气田运输方式，具有交通便利、用地少等优点，在世界范围内得到了广泛应用。 埋地管道在外爆振动载荷下的破坏，主要与爆炸所产生的振动，即爆炸地震波、管土相互作用和管内输送介质有关 综上所述，现有研究主要是将管内流体对管道的作用效果等效为内压载荷，而实际外爆振动对管道的破坏，与管道内流体流动状态密切相关。为此，本文通过建立土体-管道-管内流体的流固耦合动力学模型，对埋地管道在外爆地震波作用下的响应进行分析，通过分析管内流体不同输送压力与流体流速对管道振动的影响，探讨不同影响因素下埋地管道外爆振动响应。 1 埋地管道外爆炸流固耦合动态模型的构建 1.1 地震波传播区域 埋地炸药爆炸之后气体向外膨胀冲击土体，在土体中传播冲击波并使爆心附近(一般认为10～15倍装药半径范围内)的土体压碎、破裂，同时其波速、压力、能量随着距离的增加而很快衰减，此时冲击波传播区域为“爆炸近区”;在离开爆心10～15倍装药半径后，岩体中冲击波的能量衰减到某一临界值时，冲击波开始转换为没有陡峭波阵面的应力波(弹塑性波)，该区域称为“爆炸中区”(小于400倍装药半径);而后非弹性过程逐渐终止，显示出弹性效应，并以地震波的形式向外传播，此区域称为“爆炸远区”(400～500倍装药半径) 管道采用壳单元，土体、炸药、流体采用实体单元进行离散，管道与土体采用拉格朗日网格描述，炸药和流体采用欧拉网格描述;土体除上表面为自由边界外，其他面均施加无反射边界条件，其目的是为了消除爆炸波在边界处的反射;爆炸地震波传播方向指向管道，在管道中间截面取迎爆点A和背爆点B两点作为管道径向位移与应力响应的监测点，模型整体施加重力，有限元模型如图2所示。 外爆炸药采用TNT高能炸药材料模型 式中P———TNT爆炸压力，Pa; e———自然对数; E 管道材料本构方程采用随动硬化双线性弹塑模型，遵循von Mises屈服准则，其表达式为: 式中σ———管道应力，Pa; 管道材料选用X70钢的参数。土体材料采用Drucker-Prager材料模型，该模型可有效模拟土壤，它考虑了材料的剪切膨胀特性，允许材料各向同时硬化或软化，并可以描述土壤的长期塑性变形。表1列出土体材料参数。 管内流体采用多状态线性方程描述: 式中P———流体压力，Pa; C ρ，ρ 1.2 管土界面接触模拟 管土之间采用接触描述，由于爆炸源距管道较远，爆炸地震波的传播方向可认为与管道轴向垂直。管土接触可认为是土体表面与管道外表面的法向接触，而忽略轴向的位移;因此，将管道与土体轴向摩擦系数视为无穷大，即土体表面与管道外表面之间没有切向滑动位移，同时管土界面法向接触采用罚函数的方式进行模拟，管土接触数据传递如图3所示。管内流体与管道内壁的相互作用采用流固耦合算法，耦合界面载荷传递如图4所示。 管道与土体作用过程中，管道接触面穿透土体目标面时产生一段穿透距离x 当弹簧罚刚度k 管内壁与流体之间的耦合界面上需要满足运动学条件: 式中D D 流体的力通过下式传递到管道结构上: τ 管道结构接收到流体的力会产生变形，并通过下式反馈给流体: 式中d H d 2 爆炸波向外扩散的特性 与自然地震波不同，爆炸地震波持续时间短、衰减快、幅值大;自然地震波持续时间长而且不会在短时间内衰减，因此不能将两种载荷等同对待 由图5可以看出，炸药爆炸使土体瞬间受到挤压，在爆心处形成一个空腔，左侧土体向外形成鼓包;爆炸波向外扩散，由于左侧和下侧土体边界设定为无反射条件，图5(d)中爆炸波直接向外扩散而没有反射回来，而上侧土体为自由边界，所以出现了波反射的现象。整个过程中爆炸波加速度随时间的推移而减小，从50μs开始时的1.9158×10 由图6可以看出，爆炸地震波在管道附近持续了6 s左右，最大加速度在第2.9 s出现，其值为96 m/s 3 结构方程模型中压力下的表现 管道模型建模所根据的真实输油管道的设计压力为8 MPa，因此为了研究流动流体对管道的影响，设定管内流体的入口压力为8 MPa、流速为1 m/s。由于加速度载荷在第2.9 s时达到峰值，则取此时刻管内流体的压力和流速分布如图7所示。 从图7可以看出，流体壁面压力为8 MPa，壁面流速由于不滑移条件的存在而显示为0，内部流速为1 m/s，所以可以认为流体传递给管道壁面的压力为8 MPa。 图8示出第2.9 s时空管和充液管道的变形和应力分布。从图8(a),(b)可以看出，管道各任意截面的径向变形分布一致，说明管道在爆炸地震波作用下各个截面变形基本一致;从图8(c),(d)可以看出，管道应力分布并不均匀，这是由于管道在爆炸地震波作用下受到土体不规则的挤压力所致，而图8(d)所示应力值大于图8(c)的应力值，是由于