Numerical simulation on the dynamic response of submarine pipelines under blast loading
-
摘要: 为研究爆炸载荷对海底管道破坏作用机理,考虑爆心距、管中所充流体等因素,利用有限元软件ANSYS/LS-DYNA对爆炸载荷作用下的裸露充流体、埋地充流体、空管裸露、空管埋地4种工况下的管道分别进行动力学响应的数值模拟。研究表明:海底管道在裸露时受到的爆炸载荷作用力明显比埋地时大;在相同条件下,充流体管道比空管更能抵抗爆炸冲击;当受爆炸载荷冲击时,管道迎爆面和背爆面均可能发生破坏失效;爆心距越小,海管正对爆心垂向加速度越大,加速度峰值产生的时间越短,对管道破坏力越大。研究结果可为海底管道优化设计、安全施工及运行提供参考。Abstract: To study the damage mechanism of the blast loading on submarine pipelines, numerical simulation was conducted on the dynamic response of pipelines under blast loading using the finite element software ANSYS/LS-DYNA. Considering the distance from the blast center, the fluid in pipelines and other factors, the pipelines in the following operational modes were simulated respectively, i.e., exposed pipeline with fluid, buried pipeline with fluid, exposed pipeline without fluid and buried pipeline without fluid. The following simulation results were obtained. First, the acting force of blast loading on exposed submarine pipelines is much higher than that on buried submarine pipelines. Second, under the same conditions, pipelines with fluid are more resistant to blast shock than those without fluid. Third, both front side and back side of a pipeline may suffer damage when it is shocked by blast loading. Fourth, the shorter the distance from the blast center is, the higher the vertical acceleration of the submarine pipeline right facing the blast center is, the shorter time it takes to generate peak acceleration and the greater the destructive power on pipelines is. The research results can provide reference for the optimization design, safety construction and operation of submarine pipelines.
-
Keywords:
- submarine pipelines /
- blast loading /
- dynamic response /
- numerical simulation
-
近年来, 随着海洋石油的大力开发, 海底管道的铺设量迅速增加, 海底管道安全问题受到国内外广泛关注。在爆炸载荷对管道破坏作用方面, 国内外众多学者进行了广泛研究。王飞等[1]利用有限元软件ANSYS/LS-DYNA模拟了爆炸载荷对天然气管道破坏作用, 分析了爆炸载荷作用下管道的变形、破坏及内部应变随时间的发展过程; 刘建民等[2]模拟了埋地管道在炸药地面爆炸和侵入浅层爆炸两种情况, 对比了二者的动力响应, 以及炸药用量、爆心距、管径和设置方式等因素对埋地管道的破坏作用。李又绿等[3]根据准静定理论计算爆腔半径, 在爆腔内壁施加压力时程, 建立了岩土介质-爆腔-输气管道动力相互作用模型, 通过改变装药量、爆心距, 对比了埋地输气管道动力响应过程中管道质点速度和动应力峰值的变化情况; 纪冲等[4]对钢管在外接触爆炸载荷作用下的非线性动态响应过程进行了三维数值模拟, 得到管壁迎爆面的变形破坏及对面管壁在爆炸破片碰撞下的后效作用过程; 高进东等[5]对炸药在土中爆炸情况下埋地管道的动力响应问题, 分析了不同装药量和爆心距对埋地管道的损伤作用。Kouretzi等[6]采用解析方法, 将管道简化为薄壁弹性柱状壳, 推导出爆炸荷载作用下管道应力分布函数表达式; Jama等[7-8]采用数值模拟与试验相结合的方法, 分析了钢管的整体变形和局部变形情况。
上述研究大多集中于炸药和管道在空气或土壤中发生爆炸的影响, 而针对海底(水中)环境下爆炸载荷对管道影响的研究较少, 且极少对比分析爆炸对空管和充流体管道的影响。基于此, 运用有限元软件ANSYS/LS-DYNA[9-14]对海底管道裸露铺设、掩埋铺设以及管道是否充入流体等4种工况进行数值模拟, 进而探究海底管道在爆炸载荷下的动力响应和抗爆性能的影响因素。
1. 模型建立
1.1 有限元模型
对海底管道管内充流体裸露、埋地充流体、空管裸露、空管埋地4种工况分别建立简化的物理模型。炸药尺寸规格为10 cm×10 cm×10 cm的TNT, 将其设置于管道正上方, 炸药与管道的距离为H; 海管内所充流体采用淡水, 空管时内充空气, 初始速度均为0;海底土壤为黏土, 管道铺设于海底土壤表面或埋在土中(埋深10 cm); 管材选用X70钢管, 管道外径、内径分别为20 cm、17 cm。
1.2 网格划分
利用ANSYS/LS-DYNA建立有限元模型(图 1), 所有网格划分采用Solid164六面体实体单元[15], 采用任意拉格朗日-欧拉算法(Arbitrary Lagrangian-Eulerian Method, ALE)和流固耦合方法分析爆炸作用下海底管道的动力响应。该计算模型x方向长度为30 cm, z方向长度为200 cm, 顶部设置为自由边界, 四周设置为无反射边界, 底部施加y方向约束, 取1/4模型进行计算。炸药和海水采用Euler单元描述, 在两种材料交界面处采用共用节点的方式, 管道设置为Lagrange单元, 海水采用ALE网格, 炸药结构定义为流体, 管内流体与管道之间的相互作用运用有限元软件ANSYS/LS-DYNA中集成的流固耦合算法进行数值模拟。
1.3 材料模型及参数
1.3.1 炸药模型及参数
在爆炸载荷作用下海底管道动力响应流固耦合计算模型中, 炸药爆炸压力采用JWL状态方程以及Mat_High_Explosive_Burn材料模型[15-16]进行描述:
(1) 式中: p1为炸药爆炸压力, Pa; V为当前相对体积; E1为单位体积炸药的内能, J/m3; A、B、R1、R2、ω分别为JWL状态方程参数。
在该计算模型中, 炸药材料的各项参数分别为: 密度ρ1=1.63×103 kg/m3, 爆速D=6.93×103 m/s, A=3.712×1011 Pa, B=7.43×109 Pa, E1=7.0×109 J/m3, R1=4.15, R2=0.90, ω =0.30。
1.3.2 海水模型及参数
海水采用Mat_Null材料模型与多项式状态方程[15-16]进行描述, 其形式根据压缩状态而定。
当海水为压缩状态和膨胀状态时, 其状态方程分别为:
(2) (3) (4) (5) 式中: p2为海水压强, Pa; E2为海水的比内能, J/kg; ρ2为当前海水密度, kg/m3; ρ 0为海水初始密度, kg/m3;A1、A2、A3、B0、B1、T1、T2为AUTODYN材料库赋值的常数; h为海水深度, m; g为重力加速度, 取9.8 m/s2;p0为大气压强, Pa。
海水材料的各项参数分别为: ρ0=1.0×103 kg/m3, A1=2.2×109 Pa, A2=9.54×109 Pa, A3=1.457×1010 Pa, B0=B1=0.28, T1=2.2×109 Pa, T2=0。
1.3.3 空气模型及参数
空气采用Mat_Null材料模型以及线性多项式状态方程Eos_Linear_Polynomal[15-17]进行描述:
(6) (7) 式中: p3为空气压力, Pa; ρ3为空气当前密度, kg/m3; ρk0为空气初始密度, kg/m3; E3为材料当前单位体积内能, J/m3; C0~C6为状态方程参数。
在该计算模型中, 空气的材料的各项参数分别为: ρk0=1.29 kg/m3; E3=0.25×109 J/m3; C1=C2=C3=C6=0;C4=C5=0.4;C0=-1.0×10-6。
1.3.4 钢管模型及参数
钢管材料本构方程采用各向同性随动及混合硬化弹塑性模型Mat_Plastic_Kinematic[17-19]进行分析, 其本构关系为:
(8) (9) 式中: σy为钢管当前屈服应力, Pa; σ0为钢管初始屈服应力, Pa; ε为应变率, s-1; εeff为有效塑性应变; β为硬化参数; C、p均为Cowper-Symonds应变率参数, 其取值分别为C =40.4 s-1、p =50;Ep为塑性硬化模量, 取1;E4、Etan分别为弹性模量和切线模量, Pa。
在该计算模型中, 钢管选用X70钢管, 其相关参数分别为: 密度ρ4=1.63×103 kg/m3, 泊松比υ=0.3, E4=2.10×1011 Pa, Etan=1.35×1010 Pa, 屈服应力σs=4.80×108 Pa, 抗拉强度σb=5.40×108 Pa, 失效应变为σf=0.2。
1.3.5 海底土壤模型及参数
海底土壤采用Mat_Siol_and_Form材料模型来描述[20], 介质应力屈服函数为:
(10) 式中: f为土壤屈服应力, Pa; Sij为Cauchy偏应力张量, Pa; δij为Kronecker符号; a0、a1、a2均为土壤弹性屈服函数系数, 其取值分别为6.20×10-3、1.67×10-7、0.112 6。
对于海底土壤, vcr表示冲击条件下材料是否发生大变形, 当vcr=0表明发生大变形, 当vcr=1表明未发生大变形; ref为压力初始化时是否采用几何状态, 当ref=0表明未采用, 当ref=1表明采用。
在该计算模型中, 土壤的各项材料参数分别为: 密度ρ5=1.88×103 kg/m3; 杨氏模量E5=5.14×107 Pa; vcr、ref均取0, vcr=0时的体积模量K=2.40×1010; 土壤材料拉伸断裂后的压降p5=6.90×103 Pa。
2. 计算结果
2.1 抗爆性能
在4种工况下, 对不同爆心距(0.4~1.9 m)所对应的爆炸载荷作用在海底管道上的最大等效应力等差取值(间隔0.25 m), 结果表明: 在爆心距相同的情况下, 裸露管道受到的最大等效应力明显大于埋地管道; 当爆心距较小时, 空管受到的最大等效应力明显大于充流体管道, 表明土壤和管内流体对海底管道具有一定的保护作用(图 2)。
2.2 爆炸应力
爆心距为0.4 m的裸露空管在不同爆炸时刻的应力分布情况(图 3): 最大等效应力产生在管道正对爆心迎爆面上, 且先受到爆轰波的作用产生塑性变形, 并进一步导致背爆面拉应力逐渐增大, 同时逐渐向管壁周围扩展, 共同加剧了管道的破坏程度。
2.3 垂向质点加速度
质点加速度是衡量爆炸振动强度的有效参数。根据不同工况下正对爆心垂向加速度随时间变化情况可以得到: 爆心距越小, 海底管道正对爆心垂向加速度则越大, 相应地对管道破坏作用也越强, 且爆炸后加速度峰值产生的时间也越短; 尤其是当海底管道裸露时, 正对爆心垂向加速度大于埋入海底土壤管道的正对爆心垂向加速度, 进一步证明了掩埋铺设的海底管道比裸露海底管道更安全。
3. 工程算例
以埕岛油田某段海底管道为例, 其管道规格为Φ 304.8×12.7, 管道输送压力为4.0 MPa, 材料为碳钢X60, 管道屈服极限、弹性极限分别为289 MPa、206 MPa。
按照爆炸载荷作用下海底管道动力响应流固耦合计算模型(图 1)炸药尺寸与位置, x、y、z方向长度分别为150 cm、360 cm、150 cm。以裸露海底管道为例, 选择爆心距R分别为1 m、1.25 m、1.5 m、1.75 m及1.9 m的不同工况建立模型(图 4)。
采用最大等效应力对海底管道进行强度评估, 由不同爆心距下裸露海底管道在运行状态下受到爆炸冲击作用最大等效应力的模拟效果(图 5)可见: 正对爆心的区域受到的爆炸冲击作用最大。
采用该方法, 可以得出4种工况下海底管道最大等效应力的分布情况(表 1), 从而判定海底管道失效标准, 即管道的最大等效应力是否达到屈服强度极限。
表 1 4种工况下海底管道最大等效应力分布情况对表 1中4种工况进行数值模拟(图 6), 并运用最小二乘法得出裸露运行管道、埋地运行管道、裸露空管、埋地空管的最大等效应力与爆心距的关系式分别为σc=1.655 0×109 R-3.277 5、σc=3.348 7× 109 R-3.459 4、σc=4.563 7×108 R-2.949 1、σc=3.649 4×108 R-2.952 4。
在不同工况下, 随着爆心距增大, 海底管道的最大等效应力呈明显下降趋势。最大等效应力曲线在爆心距小于1.5 m范围内下降较快, 随后趋于平稳。取海底管道的屈服极限为289 MPa, 分别代入裸露运行管道、埋地运行管道、裸露空管、埋地空管4种工况下的最大等效应力与爆心距的关系式, 反算得出4种工况下爆心距的安全距离分别为1.15 m、1.1 m、1.26 m、1.16 m。
处于运行中的海底管道内存在介质, 吸收了部分爆炸冲击波能量, 因此, 运行的海底管道的最大等效应力在相同爆心距下明显小于空管; 同样, 埋地海底管道由于上覆土吸收了部分爆炸冲击波能量, 因此, 裸露管道的最大等效应力在相同爆心距下也明显大于埋地管道, 尤其当爆心距为1~1.5 m时数值差距较大, 随着爆心距增大, 二者差距逐渐减小。可见, 运行管道抗爆能力强于空管, 埋地管道抗爆能力优于裸露管道。
4. 结论
(1) 根据管道不同位置处所受应力、应变、加速度的计算结果, 得到了管道在爆炸载荷作用下的失效位置, 主要包括管道正对爆心处、管道正对爆心对面处、管道正对爆心截面外壁左右两侧处。
(2) 模拟得出不同工况下管道受到的最大等效应力, 运行管道抗爆能力强于空管、埋地管道抗爆能力优于裸露管道, 因此海底管道铺设时需要掩埋; 一旦出现海流冲刷致使管道裸露时, 应当尽可能施工掩埋, 同时管道停运时, 也要适当维持管道内压, 提高安全性。
(3) 根据有限元模拟结果, 拟合得出最大等效应力函数表达式, 并对管道进行强度校核, 反算出裸露运行管道、埋地运行管道、裸露空管、埋地空管爆心距的安全距离, 从而为海底管道的工程设计和施工提供参考。
-
表 1 4种工况下海底管道最大等效应力分布情况
-
[1] 王飞, 王连来, 刘广初. 爆炸荷载对天然气管道(空管)的破坏作用研究[J]. 爆破, 2006, 23(4): 20-24. doi: 10.3963/j.issn.1001-487X.2006.04.005 WANG F, WANG L L, LIU G C. Numerical simulationon damage of pipeline(empty) for natural gas under blastloading[J]. Blasting, 2006, 23(4): 20-24. doi: 10.3963/j.issn.1001-487X.2006.04.005
[2] 刘建民, 陈文涛. 爆炸荷载下埋地管道动力响应分析研究[J]. 工程爆破, 2008, 14(2): 20-25. doi: 10.3969/j.issn.1006-7051.2008.02.005 LIU J M, CHEN W T. Dynamic response study of buried pipelinesubjected to blast loads[J]. Engineering Blasting, 2008, 14(2): 20-25. doi: 10.3969/j.issn.1006-7051.2008.02.005
[3] 李又绿, 姚安林, 赵师平, 等. 爆炸载荷对埋地输气管道的动力响应和极限荷载分析[J]. 焊管, 2009, 32(11): 63-69. doi: 10.3969/j.issn.1001-3938.2009.11.014 LI Y L, YAO A L, ZHAO S P, et al. Analysis on dynamicresponse and extreme loads of buried gas pipeline under blastloading[J]. Welded Pipe and Tube, 2009, 32(11): 63-69. doi: 10.3969/j.issn.1001-3938.2009.11.014
[4] 高进东, 韦春艳, 周轶. 爆炸载荷下埋地管线动力学分析数值模拟[J]. 科技导报, 2013, 31(26): 30-34. doi: 10.3981/j.issn.1000-7857.2013.26.003 GAO J D, WEI C Y, ZHOU Y. Dynamic analysis of buriedpipelines under blast loads with numerical simulation method[J]. Science & Technology Review, 2013, 31(26): 30-34. doi: 10.3981/j.issn.1000-7857.2013.26.003
[5] 纪冲, 龙源, 唐献述, 等. 爆炸载荷下X70钢管道的局部破坏效应[J]. 高压物理学报, 2013, 27(4): 567-574. https://www.cnki.com.cn/Article/CJFDTOTAL-GYWL201304017.htm JI C, LONG Y, TANG X S, et al. Local damage effects of X70steel pipe subjected to constant explosion loading[J]. ChineseJournal of High Pressure Physics, 2013, 27(4): 567-574. https://www.cnki.com.cn/Article/CJFDTOTAL-GYWL201304017.htm
[6] KOURETZIS G P, BOUCKOVALAS G D, GANTES CJ. Analytical calculation of blast-induced strains to buried pipelines[J]. International Journal of Impact Engineering, 2007, 34(10): 1683-1704. doi: 10.1016/j.ijimpeng.2006.08.008
[7] JAMA H H, BAMBACH M R, NURICK G N, et al. Numericalmodelling of square tubular steel beams subjected to transverseblast loads[J]. Thin-Walled Structures, 2009, 47(12): 1523-1534. doi: 10.1016/j.tws.2009.06.004
[8] JAMA H H, NURICK G N, BAMBACH M R, et al. Steel squarehollow sections subjected to transverse blast loads[J]. Thin-Walled Structures, 2012, 53: 109-122. doi: 10.1016/j.tws.2012.01.007
[9] 王万鹏. 混凝土加固钢壳结构内爆响应的数值模拟研究[D]. 西安: 西北工业大学, 2005: 125-131. WANG W P. Study on numerical simulation of inside explosiveof steel shell structures reinforced by concrete[D]. Xi'an: Northwestern Polytechnical University, 2005: 125-131.
[10] 时党勇, 李裕春, 张胜民, 等. 基于ANSYS/LS-DYNA 8.1进行显示动力分析[M]. 北京: 清华大学出版社, 2004: 184-245. SHI D Y, LI Y C, ZHANG S M, et al. Explicit dynamicanalysis based on ANSYS/LS-DYNA8.1[M]. Beijing: Tsinghua University Press, 2004: 184-245.
[11] WANG G H, WANG Y X, LU W B, et al. On the determinationof the mesh size for numerical simulations of shock wavepropagation in near field underwater explosion[J]. AppliedOcean Research, 2016, 59: 1-9.
[12] 余晓菲, 刘土光, 张涛, 等. 水下爆炸载荷作用下加筋圆柱壳的响应分析[J]. 振动与冲击, 2006, 25(5): 106-115. doi: 10.3969/j.issn.1000-3835.2006.05.027 YU X F, LIU T G, ZHANG T, et al. Response analysis ofcylindrical shell subjected to underwater explosion shockwaves[J]. Journal of Vibration and Shock, 2006, 25(5): 106-115. doi: 10.3969/j.issn.1000-3835.2006.05.027
[13] O'DANIEL L J, KRAUTHAMMER T. Assessment ofnumerical simulation capabilities for medium-structure interaction systems under explosive loads[J]. Computers & Structures, 1997, 63(5): 875-887.
[14] 路胜卓, 张博一, 王伟, 等. 爆炸作用下薄壁柱壳结构动力响应实验研究[J]. 南京理工大学学报, 2011, 35(5): 621-626. doi: 10.3969/j.issn.1005-9830.2011.05.008 LU S Z, ZHANG B Y, WANG W, et al. Experimental researchon dynamic response mechanism of thin cylindrical shell underblast loading[J]. Journal of Nanjing University of Science andTechnology, 2011, 35(5): 621-626. doi: 10.3969/j.issn.1005-9830.2011.05.008
[15] 张社荣, 李宏璧, 王高辉, 等. 空中和水下爆炸冲击波数值模拟的网格尺寸效应对比分析[J]. 水利学报, 2015, 46(3): 298-306. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201503008.htm ZHANG S R, LI H B, WANG G H, et al. Comparative analysisof mesh size effects on numerical simulation of shock wave inair blast and underwater explosion[J]. Journal of Hydraulic Engineering, 2015, 46(3): 298-306. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201503008.htm
[16] 蔡珣, 张亚栋. 爆炸荷载下大口径管道破坏作用的数值模拟[J]. 工业建筑, 2014(增刊1): 236-240. https://www.cnki.com.cn/Article/CJFDTOTAL-GYJZ2014S1058.htm CAI X, ZHANG Y D. Numerical stimulation on damageof large diameter pipe under blast loading[J]. Industrial Construction, 2014(S1): 236-240. https://www.cnki.com.cn/Article/CJFDTOTAL-GYJZ2014S1058.htm
[17] 贾普荣, 李年, 朱维斗. 反向加载时X70钢的力学性能[J]. 机械强度, 2005, 27(1): 82-84. doi: 10.3321/j.issn:1001-9669.2005.01.016 JIA P R, LI N, ZHU W D. Mechanical properties of X70pipeline steel at reverse loading[J]. Journal of MechanicalStrength, 2005, 27(1): 82-84. doi: 10.3321/j.issn:1001-9669.2005.01.016
[18] 杨政, 郭万林, 董蕙茹, 等. X70管线钢静态和动态韧性试验研究[J]. 西安交通大学学报, 2003, 37(5): 488-491. doi: 10.3321/j.issn:0253-987X.2003.05.013 YANG Z, GUO W L, DONG H R, et al. Static and dynamictests of X70 pipeline steel[J]. Journal of Xi'an JiaotongUniversity, 2003, 37(5): 488-491. doi: 10.3321/j.issn:0253-987X.2003.05.013
[19] 阎石, 齐宝欣, 辛志强, 等. 高温与爆炸作用下轻钢柱动力响应与破坏模式数值分析[J]. 土木工程学报, 2010, 43(增刊1): 484-489. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2010S1084.htm YAN S, QI B X, XIN Z Q, et al. Nunmerical analysis ondynamic response and failure mode under elevated temperatureand blast loading[J]. China Civil Engineering Journal, 2010, 43(S1): 484-489. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2010S1084.htm
[20] 郭毅之, 金先龙, 丁峻宏, 等. 沉管隧道地震响应分析中的三维接触模型与算法研究[J]. 应用力学学报, 2006, 23(1): 48-52. https://www.cnki.com.cn/Article/CJFDTOTAL-YYLX200601010.htm GUO Y Z, JIN X L, DING J H, et al. Three-dimensionalcontact model and algorithm for earthquake response ofimmersed tunnel[J]. Chinese Journal of Applied Mechanics, 2006, 23(1): 48-52. https://www.cnki.com.cn/Article/CJFDTOTAL-YYLX200601010.htm