1、 收稿日期 基金项目陕西省自然科学基础研究计划面上项目“近场爆炸下 增强碳纤维加固含缺陷埋地油气管道损伤评估及抗爆优化研究”（）；陕西省教育厅重点科研计划项目（新型智库项目）“天然气集输管道风险定量评估及管理机制研究”（）。第一作者樊玉光（），男，硕士，教授，现主要从事动力载荷及响应方面的科研和教学工作，。通信作者方军（），男，硕士生，现主要从事动力载荷及响应方面的学习和研究工作，。樊玉光，方军，崔莹，等 爆炸荷载下埋地管道动力响应试验及数值模拟研究 长江大学学报（自然科学版），（）：，（），（）：爆炸荷载下埋地管道动力响应试验及数值模拟研究樊玉光，方军，崔莹，党文波，王嘉乐西安石油大学机械工
2、程学院，陕西 西安 摘要基于设计开展油气管道浅埋爆炸试验与数值模拟，对油气管道在爆炸荷载下动力响应及 （）对防护效果进行研究。通过比对油气管道浅埋爆炸试验数据，验证了采用 （）方法建立管土耦合模型的合理性和有效性。研究结果表明，在爆炸荷载下，管道迎爆面发生凹陷变形，高应力区出现在正对爆心单元并向管道端部逐步延伸衰减。外部包裹 可以有效减缓管道变形，且一定程度上增加 层数可以有效降低管道等效应变，层 包裹加固的管道等效应变较层 包裹加固管道的等效应变降低 。同时，在炸药量小于 时，浅埋爆炸荷载条件下埋地油气管道的变形不明显，爆心距大于 时，管道最大有效应力衰减较快。关键词浅埋爆炸；数值模拟；油气
3、管道 中图分类号 文献标志码 文章编号 （），：，（），（），：；长江大学学报（自然科学版）年 第 卷 第期 （），DOI:10.16772/ki.1673-1409.2023.01.003近年来各种意外爆炸事故频发，加之战争与恐怖袭击的威胁，使得爆炸荷载已成为影响埋地油气管道安全的重要第三方荷载。为了提高管道的抗爆性能，合理包裹吸能材料是一个有效的思路。碳纤维复合材料（，）由于其重量轻、耐腐蚀性、耐冲击性等优点很适合用于油气管道加固。因此，研究埋地管道在地面爆炸作用下的力学响应以及防护很有意义。许多学者对油气管道在爆炸荷载下的力学响应及防护问题进行了大量研究，夏宇磬等结合爆炸应力波传播理论，
4、推导了可表征爆炸应力波穿过土岩界面后施加至埋地管道的爆破振动荷载衰减公式。王亚光等 针对饱和土中钢顶管在地面爆炸荷载作用下的动力响应问题展开研究，探究了炸药量、起爆位置、钢顶管埋置深度以及顶管半径等参数对于钢顶管位移响应之间的关系。龚相超等 设计和实施了系列爆炸波作用下管道的动态响应实验，分析管道在不同的爆心距和爆源埋深条件下，不同位置测点峰值应变的衰减规律。等 针对爆炸荷载下加压天然气管道地面管道相互作用以及管道的响应展开研究，提出了评估爆破安全距离的实用计算方法。等 通过数值模拟对 加固埋地管道展开研究，研究表明 包覆可以显着提高钢管在爆炸载荷下的性能。等 通过数值模拟对 加固埋地管道展开
5、研究，结果表明，可以显著减少爆炸引起的管道变形，使管道的损伤区长度、截面变形和最大纵向应变分别降低、和。通过上述国内外相关研究可知，虽然验证了 加固的有效性，但目前针对油气管道与土耦合响应数值模拟方法，以及 防护优化的研究相对较少。本文基于流固耦合的方法表述管道与土的相互作用，通过试验研究和数值模拟，探讨了埋地管道在爆炸荷载作用下的力学响应以及 的防护作用，并进一步讨论了 层数、炸药量、爆心距等因素对管道响应的影响。浅埋爆炸试验 试验材料试验选择与实际工程中所采用的 管线钢屈服强度一致的 油管作为试验对象，炸药量为 ，如图所示。管道试样的具体参数如表所示。将预拌和好的碳纤维浸渍胶涂抹在 管道上
6、，选择 一级 布进行粘贴加固，具体参数如表所示。图普通管道与 管道 表 管道试件基本参数表 长度外径壁厚屈服强度 极限强度 弹性模量 剪切模量 泊松比 长江大学学报（自然科学版）年月表 基本参数表 宽幅密度（）弹性模量 伸长率 试验布置如图所示，在试验场预先开挖 （长宽深）的埋设用坑，将两组试样及反力架一同放入坑内，土体回填包括两个阶段，第一次回填界限依据爆心距为 的标准确定以埋设炸药；第二次回填界限为自然地面，通过插入钢棒方便量测。试件实际就位如图所示。图试验剖面图图试件就位图 试验结果分析对浅埋爆炸荷载作用后的普通管道和 管道变形进行测量，如图所示。由图可知，由于爆炸超压的作用，普通管道与
7、 管道的迎爆面均发生了较大的凹陷变形，布未发生破损。对两组管道迎爆面的凹陷变形分别进行了量测，其中，普通管道的变形量为 ，管道变形量为 。图普通管道和 管道变形量 数值模拟 本构模型采用 有限元软件进行数值模拟，有限元模型涉及五种材料：土体、炸药、管道、空气。第 卷 第期樊玉光 等：爆炸荷载下埋地管道动力响应试验及数值模拟研究 土体材料模型试验场地土质为黏性土，密度 ，剪切模量 ，体积模量 ，选用 的 材料模型对浅埋爆炸荷载下土体的动态力学响应进行数值模拟，其理想塑性屈服函数为：（）（）式中：，为偏应力分量，；为平均应力，；常数、为无量纲曲线二次拟合曲线常数项。炸药材料模型与状态方程采用 高能
8、炸药燃烧模型和 状态方程来模拟 炸药爆轰。炸药密度为 ，爆速 。状态方程对于爆炸过程中爆轰产物的压力、体积和能量特性可以达到精准控制，可表示为：（）（）（）式中：为爆轰产物的压力，；为爆轰产物的相对体积；为爆轰产物单位体积初始内能，；、为材料常数。管道材料模型选用 材料库中的 材料模型来模拟钢管道在爆炸荷载作用下的动力特性。材料模型 属于单向受拉材料且在破坏前可以视为线弹性材料，因此选用 材料库中的 材料模型来模拟可以忽略应变率效应的 材料。空气材料模型和状态方程选取 提供的 材料模型和来模拟空气，空气密度设置为 。空气状态方程选用 ，可表示为：（）（）（）图有限元模型（尺寸单位：）（：）式中
9、：是单位初始体积的内能，；是空气质量密度，；是 参 考 质 量 密 度，；、和为状态方程参数。管土耦合有限元模型炸药在土中爆炸对埋地管道的响应问题是一个非线性问题，利用解析计算方法研究埋地管道在浅埋爆炸作用下的破坏是非常困难的。因此，数值模拟更适合解决这个问题。为了防止大变形和非线性结构分析中的单元变形，本文采用了任意拉格朗日欧拉（，）算法建立了管土三维有限元模型，有限元模型的示意图如图所示。为了减少计算量，考虑到问题对称性，建立二分之一几何结构和荷载的分析模型。其中，管道外径为 ，管道壁厚为 。整个模型在、方向的尺寸分别为 （长宽高）。有限元网格大小为 ，空气、土壤、炸药、管道使用 单元，共
10、划分 单元 个，布使用 单元，共划分 单元 个。约束垂直于对称长江大学学报（自然科学版）年月平面（和 平面）的节点的过渡位移，土体单元四周及底部对比实际试验条件为透射边界，顶部为自由边界，以描述半无限土体。数值模拟结果验证数值模拟所获取的管道超压分布及变形如图所示，由图可知，两类管道迎爆面均发生了凹陷变形。高应力区出现在迎爆面正对爆心单元，然后向管道端部延伸衰减。由图（）至图（）可知，不同时刻两组管道受到爆炸超压后的变形情况有所不同，炸药爆炸 两组管道开始受到爆炸超压影响，至 高应力区出现在迎爆面正对爆心单元，然后向管道端部延伸衰减，迎爆面发生较为明显塑性变形，至 迎爆面塑性变形进一步加大，背
11、爆面开始变形，管道中部开始呈现微小弯曲变形，最后在 时管道背爆面发生明显塑性变形，管道中部弯曲进一步加大。图是普通管道与 管道迎爆面凹陷时程曲线图，由图可知，首先两组管道迎爆面变形起始均未在零时刻，表明管道周边土体使得爆炸超压作用在管道上的时间得到了延缓，同时，管道的变形开始时间略晚于普通管道，表明 的存在 一定程 度上 延缓了管道变形的发 生。普 通 管 道 迎 爆 面 凹 陷 变 形 数 值 模 拟 结 果 为 ，试验结果为 ，两者误差为 ，管道迎爆面凹陷变形数值模拟结果为 ，试验结果为 ，两者误差为 。两组管道凹陷变形的数值模拟结果实际误差均未超出。综上所述，数值模拟结果是合理的。图不同
12、时刻浅埋爆炸荷载下两类管道超压分布及变形 图为土体变形图，由图可知，看出在爆炸荷载作用下土体发生变形，并逐渐向两边抛掷，形成爆坑。为了研究爆炸冲击波超压在土中规律，选取沿深度及距离爆心位置折合距离为 、的四个测点，对土中超压进行观测。绘制土中爆炸产生的冲击波压力时间曲线如图所示，由图可知，各个测点的爆炸波压力随着时间的增进而快速衰减，最终趋近于零，且不同测点的爆炸波峰值压力也存在着一定差异性，埋置深度和折合距离的增加，冲击波的峰值压强相应迅速降低。第 卷 第期樊玉光 等：爆炸荷载下埋地管道动力响应试验及数值模拟研究结果与讨论 层数对管道的影响图 普通管道与 管道迎爆面凹陷时程曲线图 为了在证实
13、 布提升管道抗爆性能的基础上，进一步探索不同 层数对管道抗爆性能的影响，建立三种不同的 层数的管道在同等浅埋爆炸荷载下的数值模型，并获取管道超压分布及变形如图 所示，由图 可知，层数的增多有效地提高了防爆性能。基于 数值模拟方法，获取三种不同层数的 包裹条件下管道的最大等效应变变化情况，结果如图 所示。由图 可知，随着 层数的增多，管道受爆 炸 荷 载 后，变 形 量 减 小。同 时，用 层 包裹加固的管道的等效应变要比层 包裹加固的管道的等效应变要小 。用 层 包裹加固的管道的等效应变要比层 包裹加固的管道的等效应变要小 ，可见在一定程度上增加 层数可以有效提升管道抗爆性能。图土体变形图图土
14、中不同测点压力时间曲线 炸药量的影响进一步讨论不同炸药量条件（、）下，不同 包裹层数的管道变形情况如图 所示。由图 可知，随着 量级的增加，不同 包裹的管道变形量均相应增加。当炸药量小于 时，缠绕层 布的管道斜率较大，而缠绕层和 层 布的管道斜率变化较小。分析其原因，主要是因为土体吸收了部分爆炸能量，传递到管道上的超压较小，同时包裹的 也起到了一定的缓冲作用。爆心距的影响爆心距是影响管道变形的一个重要因素。基于数值模拟结果，获取不同 包裹管道在不同爆心距条件下的迎爆面最大等效应力变化情况如图 所示。由图 可知，随着爆心距的增加，最大等效长江大学学报（自然科学版）年月应力逐渐下降。同时，而在同一
15、爆心距条件下，层数越多，相应的最大等效应力越小。在最大等效应力衰减方面，爆心距在 时，最大等效应力衰减较缓慢，大于 时衰减较快。因此，小于 的爆心距条件可以视为浅埋爆炸荷载下影响管道变形的重要因素之一。图 浅埋爆炸荷载下不同层数 管道超压分布情况 图 不同 层数管道最大等效应变变化 结论通过设计开展浅埋爆炸试验和数值模拟，研究了爆炸荷载下埋地管道的动力响应特征，并探讨了不同 层数对管道抗爆性能的影响，所得结论如下：）基于浅埋爆炸试验结果，验证了 方法在模拟土中浅埋爆炸荷载与管道相互作用的有效性和可行性，该方法可以很好地描述爆炸波在土中的传播规律，以及管土间的耦合作用。）包裹可以有效提升油气管道
16、的抗爆性能。从同等浅埋爆炸荷载下，不同包裹层数的 管道最大等效应变的变化可知，在一定范围内增加 层数可以有效提升埋地管道的抗爆性能。）随着炸药量的提升，管道迎爆面的凹陷变形量相应提升。同时，炸药量对于超过 层 包裹的管道迎爆面变形量降低有限，说明在小于 药量条件下，如果对管道抗爆无特殊要求，可以适当减少 包裹的层数。）随着爆心距的增加，管道迎爆面最大等效应力逐渐下降。爆心距大于 时，管道最大有效应力衰减较快。因此，小于 的爆心距条件会对管道安全形成较大的威胁，需要尽量规避。第 卷 第期樊玉光 等：爆炸荷载下埋地管道动力响应试验及数值模拟研究图 不同 层数炸药量与管道凹形图 不同 层数爆心距与迎爆面最大变形量的关系等效应力的关系 参考文献 ，（）：，：丁宇奇，戴子威，芦烨，等基于多相耦合的含缺陷埋地管道外爆破坏影响因素研究 压力容器，（）：，（）：，：伍颖，田中旭，李林涯 管道凹痕缺陷失效准则对比分析 西南石油大学学报（自然科学版），（）：，（），（）：崔莹，赵均海，屈展，等 加固埋地油气管道在爆炸荷载下的损伤判定研究 振动与冲击，（）：，（）：，：，（）：夏宇磬，蒋楠，周传波，等爆破荷
