刚性柱附近浅水爆炸荷载特性研究

刘靖晗 唐廷 韦灼彬 于小存 李凌锋 张元豪

引用本文:
Citation:

刚性柱附近浅水爆炸荷载特性研究

    作者简介: 刘靖晗(1992-),男,博士研究生,主要从事港口工程、防护工程研究.E-mail:1226001717@qq.com;
    通讯作者: 唐廷, tangting1980@126.com
  • 中图分类号: O383

Pressure Characteristics of Shallow Water Explosion near the Rigid Column

    Corresponding author: TANG Ting, tangting1980@126.com ;
  • CLC number: O383

  • 摘要: 刚性柱附近浅水爆炸时冲击波传播、气泡射流受多种因素影响。考虑水面、水底、刚性柱与水下爆炸冲击波及气泡的耦合作用,基于LS-DYNA有限元软件,建立浅水爆炸全耦合模型,通过经验公式验证有限元模型的正确性。研究表明:采用炸药直径1/3~1/2中心渐变网格能够较好地保证数值模拟精度。在冲击波传播阶段,刚性柱迎爆区冲击波峰值上升并产生切断现象,冲击波下降段被“截断”,而背爆区冲击波峰值衰减约50%,同时正压作用时间增加;在气泡脉动阶段,气泡在收缩阶段产生指向刚性柱的气泡射流,当刚性柱与炸药之间的距离约为一个气泡半径时,刚性柱附近的脉冲荷载增幅最大,脉冲荷载最大测点水深较爆心上移。
  • 图 1  有限元计算模型

    Figure 1.  FEM calculation model

    图 2  不同网格尺寸下的冲击波峰值压力

    Figure 2.  Shock peak pressures simulated by different mesh sizes

    图 3  不同网格尺寸下冲击波峰值压力偏差

    Figure 3.  Error of shock peak pressures for different mesh sizes

    图 4  刚性柱有限元计算模型

    Figure 4.  FEM calculation model of rigid column

    图 5  刚性柱附近冲击波传播过程

    Figure 5.  Propagation of shock wave near a rigid column

    图 6  冲击波荷载时程曲线

    Figure 6.  Time history curve of shock wave load

    图 7  刚性柱附近冲击波峰值压力的变化

    Figure 7.  Peak pressure of shock wave near a rigid column

    图 8  刚性柱附近冲击波比冲量的变化

    Figure 8.  Specific impulse of shock wave near a rigid column

    图 9  刚性柱附近气泡射流压力等值线

    Figure 9.  Pressure peak contour of bubble impulse near a rigid column

    表 1  有限元计算模型材料参数

    Table 1.  Material parameters of FEM calculation model

    Materialρ/(kg·m–3)C0, C1, C2, C3C4C5C6E/(J·kg–1)
    Air1.2900.40.402.5×105
    Materialρ/(kg·m–3)CS1S2S3γ
    Water100014802.56–1.9860.22680.5
    Materialρ/(kg·m–3)A/GPaB/GPaωR1R2
    TNT16303747.330.34.150.95
    Materialρ/(kg·m–3)E/MPaG/MPa
    Soil180022.48
    下载: 导出CSV

    表 2  数值模拟与经验公式比较

    Table 2.  Comparison of the numerical and theoretical results

    S/mPmaxRmT
    Theoretical results/MPaNumerical results/MPaError/%Theoretical results/mNumerical results/mError/%Theoretical results/sNumerical results/sError/%
    2135.05137.361.715.005.122.400.560.535.40
    4 61.7162.391.10
    6 39.0338.780.63
    8 28.2027.213.50
    10 21.9120.805.07
    12 16.7817.213.49
    14 14.9814.026.42
    下载: 导出CSV

    表 3  刚性柱附近冲击波荷载比较

    Table 3.  Comparison of the pressure near a rigid column

    deMeasuring areaPmaxIb
    Near rigid column/MPaNo rigid column/MPaIncrease/%Near rigid column/(kN·s·m–2)No rigid column/(kN·s·m–2)Increase/%
    0.39In front of column278.68137.36102.88111.5383.7133.23
    Behind the column34.8280.88–56.9539.7757.01–30.24
    0.98In front of column87.6446.986.8743.1636.6517.76
    Behind the column19.9538.78–48.5626.3029.28–10.18
    1.95In front of column34.7421.0864.8018.7716.5913.14
    Behind the column11.9819.02–37.0113.3615.94–16.19
    下载: 导出CSV

    表 4  刚性柱迎爆区气泡脉冲荷载比较

    Table 4.  Peak pressure of bubble impulse near a rigid column

    dePmax/MPaDepth of the maximum
    impulse/m
    Increase/%
    0.399.66615.83
    0.986.33729.45
    1.951.58 6.04
    下载: 导出CSV
  • [1] 高勇军, 王伟策, 陈小波, 等. 浅层水中爆炸冲击波压力的测试与分析 [J]. 爆破, 1999(1): 9–13. doi: 10.3969/j.issn.1001-8352.1999.01.003
    GAO Y J, WANG W C, CHEN X B, et al. Testing and analysis on shock wave pressure generated by explosion under shallow water [J]. Blasting, 1999(1): 9–13. doi: 10.3969/j.issn.1001-8352.1999.01.003
    [2] 顾文彬, 叶序双, 张朋祥, 等. 浅层水中爆炸水底影响的试验研究 [J]. 解放军理工大学自然科学版, 2001, 2(2): 55–58. doi: 10.7666/j.issn.1009-3443.20010213
    GU W B, YE X S, ZHANG P X, et al. Experimental studies of bottom influence in shallow-layer water explosion [J]. Journal of PLA University of Science and Technology, 2001, 2(2): 55–58. doi: 10.7666/j.issn.1009-3443.20010213
    [3] 顾文彬, 孙百连, 阳天海, 等. 浅层水中沉底爆炸冲击波相互作用数值模拟 [J]. 解放军理工大学自然科学版, 2003, 4(6): 64–68. doi: 10.7666/j.issn.1009-3443.20030615
    GU W B, SUN B L, YANG T H, et al. Numerical simulation of explosive shockwave interaction in shallow-layer water [J]. Journal of PLA University of Science and Technology, 2003, 4(6): 64–68. doi: 10.7666/j.issn.1009-3443.20030615
    [4] 韦灼彬, 唐廷, 王立军. 港口水下爆炸荷载冲击特性研究 [J]. 振动与冲击, 2014, 33(6): 18–22.
    WEI Z B, TANG T, WANG L J. Shock characteristics of underwater explosion in port [J]. Journal of Vibration and Shock, 2014, 33(6): 18–22.
    [5] 方斌, 朱锡, 张振华. 垂直刚性面边界条件下水下爆炸气泡运动的理论研究 [J]. 海军工程大学学报, 2007, 19(2): 81–85. doi: 10.3969/j.issn.1009-3486.2007.02.018
    FANG B, ZHU X, ZHANG Z H. Theoretical study of the pulsation of an underwater explosion bubble with a vertical rigid plane [J]. Journal of Naval University of Engineering, 2007, 19(2): 81–85. doi: 10.3969/j.issn.1009-3486.2007.02.018
    [6] 方斌, 朱锡, 陈细弟, 等. 水平刚性面下方水下爆炸气泡垂向运动的理论研究 [J]. 爆炸与冲击, 2006, 26(4): 345–350. doi: 10.3321/j.issn:1001-1455.2006.04.010
    FANG B, ZHU X, CHEN X D, et al. Pulsation dynamics of an underwater explosion bubble vertical migrating to a horizontal rigid plane [J]. Explosion and Shock Waves, 2006, 26(4): 345–350. doi: 10.3321/j.issn:1001-1455.2006.04.010
    [7] 牟金磊, 朱石坚, 刁爱民, 等. 边界条件对水下爆炸气泡运动特性的影响分析 [J]. 振动与冲击, 2014, 33(13): 92–97.
    MU J L, ZHU S J, DIAO A M. Analysis on the characteristics of UNDEX bubbles under different boundary conditions [J]. Journal of Vibration and Shock, 2014, 33(13): 92–97.
    [8] 牟金磊, 朱锡, 黄晓明. 近壁面水下爆炸冲击波载荷参数研究 [J]. 海军工程大学学报, 2011, 23(1): 23–27. doi: 10.3969/j.issn.1009-3486.2011.01.005
    MU J L, ZHU X, HUANG X M. Parameters of shock waves from underwater explosion near structures [J]. Journal of Naval University of Engineering, 2011, 23(1): 23–27. doi: 10.3969/j.issn.1009-3486.2011.01.005
    [9] 朱锡, 李海涛, 牟金磊, 等. 水下近距爆炸作用下船体梁的动态响应特性 [J]. 高压物理学报, 2010, 24(5): 343–350. doi: 10.11858/gywlxb.2010.05.005
    ZHU X, LI H T, MU J L, et al. Dynamic response of characteristics of ship-like beam subjected underwater explosion in near field [J]. Chinese Journal of High Pressure Physics, 2010, 24(5): 343–350. doi: 10.11858/gywlxb.2010.05.005
    [10] 张阿漫, 姚熊亮. 复杂边界附近气泡的动态特性研究 [J]. 力学季刊, 2008, 29(1): 24–32.
    ZHANG A M, YAO X L. Dynamic of bubble near complex boundary [J]. Chinese Quarterly of Mechanics, 2008, 29(1): 24–32.
    [11] 张阿漫, 姚熊亮. 近自由面水下爆炸气泡的运动规律研究 [J]. 物理学报, 2008, 57(1): 339–353. doi: 10.3321/j.issn:1000-3290.2008.01.054
    ZHANG A M, YAO X L. The law of the underwater explosion bubble motion near free surface [J]. Acta Physica Sinica, 2008, 57(1): 339–353. doi: 10.3321/j.issn:1000-3290.2008.01.054
    [12] BENJAMIN T B, ELLIS A T. The collapse of cavitation bubbles and the pressures thereby produced against solid boundaries [J]. Philosophical Transactions of the Royal Society of London, 1966, 260(1110): 221–240. doi: 10.1098/rsta.1966.0046
    [13] RAJENDRAN R, NARASIMHAN K. Deformation and fracture behaviour of plate specimens subjected to underwater explosion–a review [J]. International Journal of Impact Engineering, 2006, 32(12): 1945–1963. doi: 10.1016/j.ijimpeng.2005.05.013
    [14] RAJENDRAN R, LEE J M. Blast loaded plates [J]. Marine Structures, 2009, 22(2): 99–127. doi: 10.1016/j.marstruc.2008.04.001
    [15] HUNG C F, HWANGFU J J. Experimental study of the behaviour of mini-charge underwater explosion bubbles near different boundaries [J]. Journal of Fluid Mechanics, 2010, 651: 55. doi: 10.1017/S0022112009993776
    [16] WARDLAW A B, LUTON J A. Fluid-structure interaction mechanisms for close-in explosions [J]. Shock and Vibration, 2015, 7(5): 265–275.
    [17] 陈永念. 舰船水下爆炸数值仿真及抗爆结构研究 [D]. 上海: 上海交通大学, 2008.
    CHEN Y N. Study on damage mechanism in ship underwater explosion and structure anti-shock [D]. Shanghai: Shanghai Jiao Tong University, 2008.
    [18] WANG G, WANG Y, LU W, et al. On the determination of the mesh size for numerical simulations of shock wave propagation in near field underwater explosion [J]. Applied Ocean Research, 2016, 38(59): 1–9.
    [19] 张社荣, 李宏璧, 王高辉, 等. 水下爆炸冲击波数值模拟的网格尺寸确定方法 [J]. 振动与冲击, 2015, 34(8): 93–100.
    ZHANG S R, LI H B, WANG G H, et al. A method to determine mesh size in numerical simulation of shock wave of underwater explosion [J]. Journal of Vibration and Shock, 2015, 34(8): 93–100.
  • [1] 孙立志李治源吕庆敖 . 四层剪切运动间隔靶板对射流干扰的数值模拟. 高压物理学报, 2013, 27(3): 423-430. doi: 10.11858/gywlxb.2013.03.017
    [2] 董琪韦灼彬唐廷张宁 . 爆炸深度对浅水爆炸气泡脉动的影响. 高压物理学报, 2018, 32(2): 024102-1-024102-9. doi: 10.11858/gywlxb.20170580
    [3] 孟祎张舵郑监 . 距离参数对浅水中气泡运动影响规律的实验和数值模拟. 高压物理学报, 2018, 32(6): 065103-1-065103-10. doi: 10.11858/gywlxb.20180583
    [4] 崔杰李世铭黄超姚熊亮张阿漫 . 射流冲击模型在水下爆炸实验中的应用研究. 高压物理学报, 2012, 26(5): 523-530. doi: 10.11858/gywlxb.2012.05.007
    [5] 言克斌黄正祥刘荣忠 . 射流侵彻陶瓷/橡胶/钢复合靶的数值仿真与实验研究. 高压物理学报, 2014, 28(4): 467-472. doi: 10.11858/gywlxb.2014.04.013
    [6] 李晓杰莫非闫鸿浩张程娇 . 爆炸焊接斜碰撞过程的数值模拟. 高压物理学报, 2011, 25(2): 173-176 . doi: 10.11858/gywlxb.2011.02.014
    [7] 黄超汪斌刘仓理张阿漫姚熊亮 . 非球形水下爆炸气泡坍塌机制. 高压物理学报, 2012, 26(5): 501-507. doi: 10.11858/gywlxb.2012.05.004
    [8] 李健荣吉利 . 近壁面气泡运动特性的数值计算. 高压物理学报, 2010, 24(3): 168-174 . doi: 10.11858/gywlxb.2010.03.002
    [9] 王新征张松林邹广平 . 内部短药柱爆炸作用下钢筒破裂特征的数值分析. 高压物理学报, 2010, 24(1): 61-66 . doi: 10.11858/gywlxb.2010.01.011
    [10] 缪广红李亮江向阳刘文震李雪交汪泉余勇沈兆武 . 双面爆炸焊接的数值模拟. 高压物理学报, 2018, 32(4): 045202-1-045202-8. doi: 10.11858/gywlxb.20180513
    [11] 缪广红王章文李亮江向阳刘文震程扬帆汪泉余勇马宏昊沈兆武 . 爆炸复合边界效应的数值模拟. 高压物理学报, 2017, 31(1): 93-96. doi: 10.11858/gywlxb.2017.01.014
    [12] 李彪彪王辉袁宝慧 . 水防护层对射流头部速度衰减作用研究. 高压物理学报, 2018, 32(3): 035105-1-035105-5. doi: 10.11858/gywlxb.20170691
    [13] 张钟文黄正祥祖旭东贾鑫 . 异型夹层结构复合靶板抗射流侵彻性能研究. 高压物理学报, 2013, 27(6): 928-935. doi: 10.11858/gywlxb.2013.06.021
    [14] 陈昊陶钢 . 铜射流侵彻穿孔处的温度及微观组织研究. 高压物理学报, 2011, 25(4): 344-350 . doi: 10.11858/gywlxb.2011.04.010
    [15] 吴鹏李如江阮光光周杰雷伟聂鹏松石军磊于金升李优赵海平 . 弹着点位置对V形反应装甲干扰射流的影响. 高压物理学报, 2018, 32(1): 015105-1-015105-7. doi: 10.11858/gywlxb.20170568
    [16] 曹玉忠卢泽生管怀安张幼平 . 抗爆容器内爆炸流场数值模拟. 高压物理学报, 2001, 15(2): 127-133 . doi: 10.11858/gywlxb.2001.02.009
    [17] 刘伟李火坤刘成梅刘玮琳 . 基于FLUENT的动态高压微射流内部孔道流场的数值模拟. 高压物理学报, 2012, 26(1): 113-120. doi: 10.11858/gywlxb.2012.01.017
    [18] 朱孙科陈二云马大为乐贵高 . 非定常欠膨胀射流的正格式数值模拟. 高压物理学报, 2011, 25(5): 429-434 . doi: 10.11858/gywlxb.2011.05.008
    [19] 陈伟马宏昊沈兆武薛冰 , . 起爆方式对环向聚能射流成型影响的数值模拟. 高压物理学报, 2015, 29(6): 419-424. doi: 10.11858/gywlxb.2015.06.003
    [20] 何远航李海军张庆明 . 高速水射流与凝聚炸药相互作用的数值模拟. 高压物理学报, 2005, 19(2): 169-173 . doi: 10.11858/gywlxb.2005.02.012
  • 加载中
图(9)表(4)
计量
  • 文章访问数:  322
  • 阅读全文浏览量:  298
  • PDF下载量:  8
出版历程
  • 收稿日期:  2018-12-20
  • 录用日期:  2019-02-25
  • 网络出版日期:  2019-09-18
  • 刊出日期:  2019-10-01

刚性柱附近浅水爆炸荷载特性研究

    作者简介:刘靖晗(1992-),男,博士研究生,主要从事港口工程、防护工程研究.E-mail:1226001717@qq.com
    通讯作者: 唐廷, tangting1980@126.com
  • 1. 海军工程大学,湖北 武汉 430033
  • 2. 海军勤务学院,天津 300450
  • 3. 联勤保障部队第四工程代建管理办公室,辽宁 沈阳 110005

摘要: 刚性柱附近浅水爆炸时冲击波传播、气泡射流受多种因素影响。考虑水面、水底、刚性柱与水下爆炸冲击波及气泡的耦合作用,基于LS-DYNA有限元软件,建立浅水爆炸全耦合模型,通过经验公式验证有限元模型的正确性。研究表明:采用炸药直径1/3~1/2中心渐变网格能够较好地保证数值模拟精度。在冲击波传播阶段,刚性柱迎爆区冲击波峰值上升并产生切断现象,冲击波下降段被“截断”,而背爆区冲击波峰值衰减约50%,同时正压作用时间增加;在气泡脉动阶段,气泡在收缩阶段产生指向刚性柱的气泡射流,当刚性柱与炸药之间的距离约为一个气泡半径时,刚性柱附近的脉冲荷载增幅最大,脉冲荷载最大测点水深较爆心上移。

English Abstract

  • 在水下爆炸过程中,边界条件对冲击波传播和气泡脉冲的影响很大。高勇军[1]、顾文彬[2-3]等先后通过浅水爆炸试验和数值仿真方法研究了浅水环境中水下冲击波形成和传播过程,系统分析了自由面和水底界面对气泡脉动和冲击波传播的影响,发现了水面和水底反射稀疏波的切断效应。韦灼彬等[4]采用数值仿真软件模拟港口水下爆炸冲击波的传播过程,在考虑水面和水底影响的情况下,分析冲击波和气泡脉冲沿港池水深的变化规律,发现浅水环境流场边界复杂,不仅包含水面、水底边界,还包含码头、桥梁等刚性边界。方斌[5-6]、牟金磊[7-8]、朱锡[9]等以气泡脉动为重点,分别研究水平、垂直刚性壁面以及船型梁附近气泡的脉动运动,得到了气泡脉动解析解和射流规律,张阿漫等[10-11]考虑自由面、刚性边界等约束对水下爆炸气泡脉动的影响,利用边界元法模拟近边界水下爆炸,建立了气泡、壁面和自由面三者之间的耦合动力学模型。国外学者Benjamin等[12]通过高速摄影技术,首次拍摄了刚性界面附近的气泡射流现象。Rajendran等[13-14]进行了一系列钢板附近水下爆炸试验,研究不同爆距条件下钢板的响应及破坏现象。Hung等[15]进行了复杂边界下小当量水下爆炸试验,观察到近水面、刚性壁、弹性壁附近气泡的运动特征。Wardlaw等[16]考虑弹性、刚性界面等边界因素,分析了冲击波反射、气穴效应、气泡脉动以及二次荷载等现象,发现刚性界面会加大气泡和冲击波荷载,从而促进气穴效应和毁伤效应。从上述研究中可以看出,刚性界面对冲击波传播和气泡脉冲均有较大影响,已有工作主要集中在刚性壁边界,对高桩码头、桥梁等透空式结构开展刚性柱附近浅水爆炸荷载特性研究具有一定意义。

    本研究考虑水面、水底、刚性柱边界影响,通过LS-DYNA有限元软件对浅水爆炸荷载传播问题进行分析,与经验公式对比验证数值仿真的准确性,得到冲击波传播过程和气泡射流现象,总结刚性柱附近冲击波及气泡脉冲荷载特性,为进一步开展复杂边界条件下水下爆炸毁伤效应研究打下基础。

    • 通过有限元软件LS-DYNA模拟浅水爆炸,设水下深度为20 m,空气域高度为5 m,水底黏土厚度为5 m,100 kg TNT球形装药在距水面10 m深处,模型尺寸为30 m×30 m×20 m,有限元模型如图1所示。其中空气、水、炸药采用欧拉单元,黏土采用拉格朗日单元,黏土区域重叠空白材料(Mat_Ale_Vacuum),流域边界采用环境单元(Eulerian Ambient),保证四周无边界约束,并初始化重力场和静水压力。

      图  1  有限元计算模型

      Figure 1.  FEM calculation model

      假设空气、水、炸药均为均匀连续,空气采用线性多项式状态方程,C0C6为状态方程参数,E为初始单位质量内能;水采用Grüneisen状态方程,CS1S3为状态方程参数,γ为Grüneisen常数;炸药采用标准JWL方程,ABωR1R2为状态方程参数;黏土采用线弹性模型,E为弹性模量,G为剪切模量,详细材料参数如表1所示。

      Materialρ/(kg·m–3)C0, C1, C2, C3C4C5C6E/(J·kg–1)
      Air1.2900.40.402.5×105
      Materialρ/(kg·m–3)CS1S2S3γ
      Water100014802.56–1.9860.22680.5
      Materialρ/(kg·m–3)A/GPaB/GPaωR1R2
      TNT16303747.330.34.150.95
      Materialρ/(kg·m–3)E/MPaG/MPa
      Soil180022.48

      表 1  有限元计算模型材料参数

      Table 1.  Material parameters of FEM calculation model

    • 考虑浅水爆炸水面、水底反射主要影响冲击波比冲量,冲击波峰值与无限水域冲击波荷载一致,基于大量试验得到水下爆炸冲击波峰值的经验公式[17]

      ${P_{\max }} = 52.16{\left(\frac{{{W^{1/3}}}}{S}\right)^{1.13}}$

      式中:Pmax为冲击波峰压,单位MPa;W为炸药装量,单位kg;S为测点距离,单位m。

      有限元模拟的准确性与单元网格尺寸密切相关[18-19],选取炸药网格尺寸分别为10、15、20和40 cm,研究网格尺寸对仿真结果的影响。4种网格尺寸数值仿真结果以及经验公式中冲击波峰值压力随测距变化规律见图2。随着网格尺寸增加,近场冲击波峰值压力减小明显,而远场冲击波仅有略微减小,冲击波衰减速度减缓,这是由于网格较大时冲击波高频段损失所导致的。图3显示了5种网格尺寸的数值模拟结果与经验公式的相对偏差。冲击波峰值压力与经验公式计算结果的相对偏差随网格尺寸增加呈先减小后增大的趋势,网格尺寸在15~20 cm范围内,相对偏差较小,但随测距变化而产生波动,因此本数值仿真采用15~20 cm网格尺寸,即1/3~1/2炸药直径的中心渐变网格,所得仿真结果与经验公式计算结果的相对偏差随测距的变化在允许的误差范围之内,从而确保数值仿真的准确性。

      图  2  不同网格尺寸下的冲击波峰值压力

      Figure 2.  Shock peak pressures simulated by different mesh sizes

      图  3  不同网格尺寸下冲击波峰值压力偏差

      Figure 3.  Error of shock peak pressures for different mesh sizes

    • 自由水域水下爆炸时,气泡膨胀最大半径Rm和气泡第一次脉动周期T的计算公式[13]

      ${R_{\rm m}} = 3.36{\left(\frac{W}{{H + 10}}\right)^{1/3}}$

      $T = 2.08\frac{{{W^{1/3}}}}{{{{(H + 10)}^{5/6}}}}$

      式中:H为水深,单位m;Rm为气泡最大半径,单位m;T为气泡第一次脉动周期,单位s。

      表2为100 kg TNT炸药距水面10 m爆炸时气泡最大半径、第一次脉动周期、冲击波峰值压力的数值模拟结果与理论结果比较。气泡脉动和冲击波峰值压力的数值模拟结果与经验公式计算结果较为吻合,在S=14 m时,相对偏差最大为6.42%,说明数值仿真能清晰描述气泡脉动和冲击波荷载规律。

      S/mPmaxRmT
      Theoretical results/MPaNumerical results/MPaError/%Theoretical results/mNumerical results/mError/%Theoretical results/sNumerical results/sError/%
      2135.05137.361.715.005.122.400.560.535.40
      4 61.7162.391.10
      6 39.0338.780.63
      8 28.2027.213.50
      10 21.9120.805.07
      12 16.7817.213.49
      14 14.9814.026.42

      表 2  数值模拟与经验公式比较

      Table 2.  Comparison of the numerical and theoretical results

    • 针对浅水环境分别在距离炸药2、5、10 m处设置刚性柱,建立基于任意拉格朗日欧拉法(ALE)的有限元模型,如图4所示。系统中设立尺寸为1 m×1 m×20 m的全约束钢性柱,炸药装药量、刚性柱与炸药距离、刚性柱截面尺寸均为影响水下爆炸荷载结果的重要因素。为了统一变量,定义无量纲比例距离为距离参数与气泡最大半径Rm的比值,3种工况下刚性柱与炸药的比例距离de分别为0.39、0.98和1.95,刚性柱截面的比例长度为0.2。空气、水、空白物质采用多物质ALE算法,土壤采用拉格朗日算法,通过关键字*Constrained_Lagrange_in_Solid定义流固耦合,模拟炸药、水、空气与刚性柱之间的流固耦合作用。

      图  4  刚性柱有限元计算模型

      Figure 4.  FEM calculation model of rigid column

    • 起爆后,炸药瞬间变为高温高压的气态物质(爆轰产物),急剧膨胀压缩水体产生初始冲击波向四周传播。冲击波传播过程中受到刚性柱影响会产生反射、透射、绕射现象。图5为距刚性柱de=0.39处水下爆炸冲击波传播过程:t=1.49 ms时,刚性柱附近的冲击波产生反射现象;t=3.29 ms时,反射稀疏波与入射波在刚性壁迎爆区叠加使峰值压力骤降并产生切断现象,同时冲击波发生绕射在柱后汇聚继续传播;t=7.20 ms时,绕射冲击波在传播一定距离后压力云图与左侧冲击波相近。

      图  5  刚性柱附近冲击波传播过程

      Figure 5.  Propagation of shock wave near a rigid column

      图6(a)图6(b)分别为刚性柱的迎爆区(S=2 m)和背爆区(S=3 m)冲击波压力时程曲线:迎爆区冲击波在达到峰值荷载后下降段出现明显“截断”,即冲击波切断现象,冲击波比冲量衰减;背爆区冲击波下降段变缓,因此冲击波绕射延长了正压作用时间。

      图  6  冲击波荷载时程曲线

      Figure 6.  Time history curve of shock wave load

      在距炸药比例距离de分别为0.39、0.98和1.95处设置刚性柱,冲击波峰值压力沿测距S的变化规律见图7。对比无刚性柱边界条件下冲击波峰值压力变化规律可以发现:刚性柱迎爆区冲击波峰值压力较自由场增加,背爆区较自由场减小。当冲击波传播到刚性柱附近时,刚性柱反射的稀疏波与入射波在迎爆区叠加导致压力上升;冲击波经过反射后绕射能量衰减,因此刚性柱背部压力下降。图8显示了刚性柱附近冲击波比冲量Ib的变化规律。刚性柱附近冲击波比冲量的变化规律与冲击波峰值压力基本一致。随着刚性柱与炸药比例距离增加,冲击波传播过程中受水面和水底反射影响,冲击波传播方向复杂,受刚性柱反射、绕射影响减弱,因此,刚性柱附近冲击波峰值荷载以及比冲量变化不大。

      图  7  刚性柱附近冲击波峰值压力的变化

      Figure 7.  Peak pressure of shock wave near a rigid column

      图  8  刚性柱附近冲击波比冲量的变化

      Figure 8.  Specific impulse of shock wave near a rigid column

      表3列出了刚性柱附近冲击波峰值和比冲量较无刚性柱增幅(衰减)。可见:与无刚性柱时相比,刚性柱附近冲击波峰值变化显著;由于迎爆区切断现象“截断”了冲击波荷载,背爆区绕射增加了荷载时长,冲击波比冲量变化幅度较小,并且随着刚性柱比例距离的增加,刚性柱影响减弱。

      deMeasuring areaPmaxIb
      Near rigid column/MPaNo rigid column/MPaIncrease/%Near rigid column/(kN·s·m–2)No rigid column/(kN·s·m–2)Increase/%
      0.39In front of column278.68137.36102.88111.5383.7133.23
      Behind the column34.8280.88–56.9539.7757.01–30.24
      0.98In front of column87.6446.986.8743.1636.6517.76
      Behind the column19.9538.78–48.5626.3029.28–10.18
      1.95In front of column34.7421.0864.8018.7716.5913.14
      Behind the column11.9819.02–37.0113.3615.94–16.19

      表 3  刚性柱附近冲击波荷载比较

      Table 3.  Comparison of the pressure near a rigid column

    • 水下爆炸产生冲击波时会伴随气泡产生。当气泡在刚性柱附近运动时,膨胀阶段的气泡受刚性柱阻碍,在刚性柱与气泡之间形成低压区;收缩阶段时,气泡由于压力差向刚性柱移动,远离刚性柱的气泡收缩速度较大并向气泡内部凹陷,产生指向刚性柱的射流,在射流方向流场压力骤增,对结构产生脉冲荷载。

      图9为刚性柱附近气泡射流时的流场压力等值线。气泡发生射流时形成指向刚性柱的脉冲荷载,随着刚性柱与炸药比例距离增加,气泡射流现象趋近无刚性柱条件。

      图  9  刚性柱附近气泡射流压力等值线

      Figure 9.  Pressure peak contour of bubble impulse near a rigid column

      选取5~10 m水深的测点,比较刚性柱迎爆区气泡脉冲荷载峰值,如表4所示。刚性柱附近水下爆炸气泡受浮力和刚性柱Bjerknes力的共同作用,气泡上浮同时发生射流。刚性柱距炸药较近时,刚性柱迎爆区6~7 m水深测点处气泡脉冲荷载最大;随着比例距离增加,气泡脉冲荷载最大的测点水深(危险深度)下移;当刚性柱与炸药比例距离大于1.98时,不同测深的气泡脉冲峰值基本一致。当刚性柱与炸药比例距离为0.98时,脉冲荷载较自由场增幅最大。这是由于当刚性柱与炸药的距离小于气泡最大半径时,气泡受刚性柱约束无法充分膨胀,气泡射流时刚性柱对脉冲荷载的增强效应有限;当刚性柱与炸药距离大于气泡最大半径时,随着比例距离增加,刚性柱对气泡脉动的影响减弱,脉冲荷载增幅减小。

      dePmax/MPaDepth of the maximum
      impulse/m
      Increase/%
      0.399.66615.83
      0.986.33729.45
      1.951.58 6.04

      表 4  刚性柱迎爆区气泡脉冲荷载比较

      Table 4.  Peak pressure of bubble impulse near a rigid column

    • 通过数值模拟方法研究刚性柱附近浅水爆炸冲击波和气泡射流的荷载特性,与经验公式对比验证得出下述结论。

      (1)网格尺寸对水下爆炸冲击波荷载的数值模拟精度有显著影响,随着网格尺寸增加,冲击波峰值呈现下降趋势,并且冲击波衰减速度减慢,采用炸药直径1/3~1/2的中心渐变式网格能够得到较为准确的数值模拟结果。

      (2)在水下爆炸冲击波阶段,冲击波在刚性柱附近发生反射、绕射现象。刚性柱迎爆区冲击波峰值压力显著增加,增幅约60%~100%;由于冲击波切断现象,比冲量增幅较小,约10%~30%。刚性柱背爆区压力降低约50%,同时由于正压作用时间增加,比冲量衰减约20%。

      (3)在水下爆炸气泡脉动阶段,气泡射流形成指向刚性柱的脉冲荷载,刚性柱附近脉冲荷载增加。在刚性柱距炸药一个气泡半径距离时,刚性柱附近的气泡脉冲峰值增幅最大,气泡射流荷载最大的测点深度较爆心上移;随着炸药与刚性柱距离的增加,气泡射流现象减弱,刚性柱与炸药距离大于2倍气泡半径时,刚性柱附近不同测深的脉冲荷载相差不大。

参考文献 (19)

目录

    /

    返回文章
    返回