• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

基于面积折减等效模型的光电倍增管水下内爆机理研究

孟令存 闫明 杜志鹏 张磊

张凯, 奚进一, 高均波. 多层爆炸复合中的成波研究[J]. 爆炸与冲击, 1991, 11(1): 37-45. doi: 10.11883/1001-1455(1991)01-0037-9
引用本文: 孟令存, 闫明, 杜志鹏, 张磊. 基于面积折减等效模型的光电倍增管水下内爆机理研究[J]. 爆炸与冲击, 2020, 40(8): 082102. doi: 10.11883/bzycj-2019-0436
MENG Lingcun, YAN Ming, DU Zhipeng, ZHANG Lei. Underwater implosion mechanism of PMT area reduction equivalent model[J]. Explosion And Shock Waves, 2020, 40(8): 082102. doi: 10.11883/bzycj-2019-0436
Citation: MENG Lingcun, YAN Ming, DU Zhipeng, ZHANG Lei. Underwater implosion mechanism of PMT area reduction equivalent model[J]. Explosion And Shock Waves, 2020, 40(8): 082102. doi: 10.11883/bzycj-2019-0436

基于面积折减等效模型的光电倍增管水下内爆机理研究

doi: 10.11883/bzycj-2019-0436
基金项目: 国家自然科学基金(10672181)
详细信息
    作者简介:

    孟令存(1994— ),男,硕士研究生,1605733924@qq.com

    通讯作者:

    闫 明(1978— ),男,博士,教授,yanming7802@163.com

  • 中图分类号: O381;U663.85

Underwater implosion mechanism of PMT area reduction equivalent model

  • 摘要: 光电倍增管(photomultiplier tube,PMT)是中微子探测器的核心部件,是由玻璃材料制成的内部真空的薄壳结构,排列在深水中工作,若一个PMT被压溃会产生内爆冲击波,会引起周围PMT发生殉爆。针对PMT内爆,建立了PMT内爆数值计算简化模型,并将计算与试验结果进行对比,验证简化模型的合理性。在此基础上,提出了基于面积折减等效模型的PMT内爆计算方法,通过等效模型分析了防护装置破口面积对PMT内爆的影响,得出随着防护装置破口面积的减小,水流碰撞发生PMT内爆的时刻相应提前,内爆产生的冲击波脉宽基本保持不变,冲击波压力峰值明显减小。该研究有利于找到有效的PMT内爆防护方法。
  • 炸药在密闭空间内爆炸,冲击波在空间内往返运动,并在该过程中逐渐减弱并消失,之后空间内的压力分布较均匀,随时间变化缓慢,这种压力状态被称为准静态压力状态。准静态气压状态下空间内的压力因气体与外界热交换而缓慢变化,而一般受关注的是最初的准静态压力值或者准静态压力峰值,下文将准静态压力峰值简称为准静态压力。

    准静态压力受到的关注远不及爆炸冲击波形成的脉冲压力,但其仍是非常重要的力学参数,比如:准静态压力是爆炸容器安全设计的重要参考数据[1-2];另外,对于地下核爆炸,爆室内的准静态压力是核爆炸形成的剧毒物质在岩土孔隙和裂缝中扩散的动力[3]。国外实验研究工作较丰富[4-8],Anderson等[7]基于文献中的实验数据,利用相似理论拟合得到容器内爆炸准静态气压无量纲峰值的经验公式,Anthony等[8]通过实验得到了炸药LX-13在容器内爆炸形成的准静态压力经验公式,Bultman[1]在容器设计时也引用了基于实验数据的经验公式,但这些经验公式差异较大。国内的研究工作不多,特别是实验工作比较缺乏。王等旺等[9]、林俊德[10]开展了柱形容器内准静态压力的实验研究,钟方平[11]通过数值计算研究了柱形容器内的准静态压力,钟巍等[12]研究了考虑化学反应动力学影响的约束爆炸准静态压力的计算方法。

    本文中,在球形爆炸容器内开展了4种不同当量的爆炸加载实验,实验中压力传感器采用了2种安装方式,均获得了理想数据;还理论推导了准静态压力的表达式,并通过拟合实测数据得到具体的经验公式。

    实验在内径523 mm的球形爆炸容器内进行,容器如图 1所示。炸药采用球形装药,装药包括主装药和起爆微型药球两部分,如图 2所示。主装药由2个半球组成,为60%RDX和40%TNT浇铸而成。起爆药球的当量为1 g TNT,直径为10 mm,由微米级的PETN粉压制而成。起爆药球连接一根柔爆索,其中心为柔爆索的一端。柔爆索直径为1 mm,线装药密度为0.5 g RDX/m,为铅皮柔爆索。

    图  1  实验用球形爆炸容器
    Figure  1.  The spherical explosion containment vessel used in experiment
    图  2  实验装药
    Figure  2.  Explosive installation

    实验时采用细绳吊挂装药,装药置于容器中心。容器内除了炸药、导爆索和细绳,再无其他物体,这样的做法可以尽量减少容器内除炸药以外的多余物体。雷管处于雷管罩内,与容器内部体积隔开。由雷管引爆较细的柔爆索,柔爆索在中心起爆微型药球,微型药球再起爆球形主装药,该起爆方式能保证形成较理想的球面冲击波。采用27、60、185、370 g TNT 4种当量,对应的比距离范围为0.36~0.87 m/kg1/3

    容器赤道面上均布4个用于传感器安装的接口。压力传感器安装包括齐平安装和导孔安装2种方式,如图 3所示。齐平安装方式中传感器敏感面与容器内壁齐平;导孔安装方式设计了一个空腔,传感器敏感面处于空腔内,该空腔通过6个直径1 mm的导孔连接容器。导孔安装方式可以避免传感器遭受强冲击波或高速破片的破坏,但也导致整个测试的频响降低[13],因此该方式仅适合准静态压力这种变化较缓慢的参数测试。

    图  3  压力传感器安装结构示意图
    Figure  3.  Diagrams of pressure sensor installations

    实验中采用了某型压电式压力传感器和某型压阻式压力传感器,参数见表 1。由于实验中压阻式传感器在采用齐平安装方式下容易损坏或实测信号质量不高,最终压阻式传感器采用导孔安装方式,压电式传感器采用齐平安装方式。2种压力测量系统通过同步机实现时间同步。

    表  1  压力传感器参数
    Table  1.  Performance parameters of pressure sensors
    类型 量程/MPa 频响/kHz 上升时间/μs 非线性度 精度
    压电式传感器 34.48, 68.95 500 ≤1.0 ≤1.0%Fs 0.14, 0.28 kPa
    压阻式传感器 15, 40 600 ≤0.3 ≤1.0%Fs 0.5%Fs
    下载: 导出CSV 
    | 显示表格

    图 4为齐平安装压力传感器测得的典型波形。可见压力波形的前期(爆炸后1 ms之前)存在3次较明显的脉冲,表明炸药在容器中心爆炸后,形成的冲击波在容器壁和容器中心之间往返运动,在容器壁上反射了3次;在3个脉冲之后,压力波形变化较缓,表明容器内的压力趋于平静,进入准静态压力状态。当然,在比距离比本文的更小的工况下,也可能出现冲击波往返次数多于3次,容器内压力进入准静态状态的现象。准静态压力的取值是对压力波形的准静态压力部分进行平均处理,图 4中冲击波反射压力达到约16 MPa,准静态压力仅为1.5 MPa。

    图  4  齐平安装传感器获得典型压力波形(当量60 gTNT)
    Figure  4.  Pressure-time curve obtained by the pressure sensor of flush installation (when explosion equivalent is 60 gTNT)

    图 5为导孔安装压力传感器获得的典型压力波形。可见波形早期有一个压力最大值(约3.5 MPa),之后压力缓慢变化,且压力在初期时间内(爆炸后0.2 s之前)近似地随时间线性下降。

    图  5  导孔安装传感器获得的典型波形(当量60 gTNT)
    Figure  5.  Pressure-time curve obtained by the pressure sensor of guide-hole installation (when explosion equivalent is 60 gTNT)

    虽然导孔安装传感器未与容器内部空间直接接触,而是通过6个直径1 mm的小孔连接,但爆炸冲击波通过小孔仍会对传感器加载一个较弱的冲击载荷。因此,简单地读取图 5中压力波形的最大值作为准静态压力值是不合适的,该峰值应该是冲击波压力在低频测试系统中的表现。准静态气压是由冲击波压力连续变化而来的,如何确定准静态气压值显得比较困难,这可能也是不同文献中实测准静态气压差别较大的原因之一。本文的取值方法如图 6所示,按照线性下降斜率进行反向回推,与波形上升沿交于一个点,该点的压力值即为准静态气压。这种取值方法也是Anderson等[7]推荐的方法。

    图  6  准静态压力值的选取方法
    Figure  6.  The method for determining quasi-static pressure

    表 2为两种安装方式的传感器测量到的准静态压力值,可见两种测量方式得到的准静态压力值相近。

    表  2  压力传感器参数
    Table  2.  Quasi-static pressures obtained by pressure sensors
    当量/(gTNT) 准静态压力/MPa
    压电式传感器 压阻式传感器
    27 0.7, 0.7, 0.6
    60 1.4, 1.5 1.5, 1.2
    185 2.4 3.7
    370 5.5, 4.9, 5.3 5.4
    下载: 导出CSV 
    | 显示表格

    值得注意的是,由于冲击波的作用,压力波形的早期存在一个远高于准静态压力值的压力峰,如图 4~5所示。传感器的量程必须根据该压力峰值来进行选取,但传感器量程过高必然导致准静态压力值的测量精度下降。

    密闭空间内的准静态气压主要是爆炸气体产物在爆炸释放热量的作用下温度升高而引起的气体热运动压力。影响准静态气压的主要因素包括两部分:一是由爆炸形成的大量气体产物使得密闭空间内气体量增大;二是爆炸释放的热量使气体的温度升高。下面基于理想气体的假设,简单地推导密闭空间内准静态压力的表达式。

    理想气体的状态方程为:

    p=ρRT
    (1)

    式中:ρ为气体密度;R为理想气体参数,为普适气体参数R0除以气体摩尔质量M,即R=R0/MT为绝对温度。

    忽略气体与外界热交换带来的能量损失,且忽略空间内介质气体的质量(因为其质量远小于炸药的质量),同时假设爆炸释放的能量全部用来加热球形容器内的气体物质(因为大部分爆炸能量转换为了气体内能),则气体的温升为:

    ΔT=EQcV
    (2)

    式中:E为爆炸释放的总能量,Q为爆炸当量,cV为容器内气体的平均定容比热。

    实验初始状态为常温、常压,因此爆炸后密闭空间内的准静态压力为:

    pq=ρRT=QVR0M(T0+ΔT)=QVR0M(T0+EQcV)=QV(R0T0M+eR0McV)
    (3)

    式中:T0为初始温度,V为容器体积,e为炸药释放的质量能量,即比内能。对于特定的炸药来说,气体摩尔质量M、定容比热cV、比内能e皆为常数,从式(3)可以看出,空间内准静态压力与当量容积比Q/V呈正比例关系。该关系式与Pastrnak等[5]采用的表达式形式相同。

    根据表 2中数据,绘制准静态压力与当量容积比Q/V的变化关系,如图 7所示。可见,2种安装方式传感器获得的数据变化规律一致。另外,也可以看出,球形爆炸容器内准静态压力与当量容积比Q/V呈近似的线性关系。

    图  7  球形容器内准静态压力与当量容积比的关系
    Figure  7.  Relation between quasi-static pressure and explosion equivalent-to-vessel volume ratio in spherical vessel

    基于式(3)和表 2中的实测数据,可以拟合得出准静态压力经验公式:

    pq=1.10Q/V0.36kgTNT/m3Q/V4.94kgTNT/m3
    (4)

    式中pq单位为MPa,当量单位为kg TNT,容器体积单位为m3

    在理论推导准静态压力表达式—式(3)的过程中,忽略了容器内介质气体的质量,因为大部分情况下介质气体的质量远小于炸药质量,但当炸药质量较小,与容器内介质气体质量相当时,这种处理方式是不合适的,此时经验公式(4)不再适用。另外,文献[10]表明,容器内的氧气环境对准静态压力有影响,在爆炸后的高温环境中,爆轰产物与氧气化合燃烧,增大容器内气压,如果爆炸当量比本文实验的爆炸当量更大,可能导致氧气相对不足,爆轰产物与氧气的化合燃烧不充分,此时经验公式(4)也不再适用。因此,基于现有的实验结果,当量容积比范围0.36~4.94 kg TNT/m3是式(4)成立的重要前提条件,如果当量容积比不在该范围,式(4)可能不再适用。

    再者,式(3)显示,炸药种类不变的情况下准静态压力和当量容积比之间关系的线性系数才是一个常数。因此,式(4)可能不适用于性质与本文实验所用炸药相差较大的炸药。

    在球形容器内开展了4种当量的爆炸加载实验,实验中压力传感器采用了齐平和导孔两种安装方式,均获得了准静态压力数据,且两组数据一致。研究结果表明:

    (1) 在比距离0.36~0.87 m/kg1/3、当量容积比0.36~4.94 kg TNT/m3范围内,爆炸冲击波在球形容器内往返3次后,容器内气体进入准静态压力状态。

    (2) 球形容器内准静态压力与当量容积比呈正比例关系,系数为1.10 MPa·m3/ kg TNT。

  • 图  1  光电倍增管实物图

    Figure  1.  Picture of photomultiplier tube

    图  2  内爆试验装置示意图

    Figure  2.  Schematic of test device for implosion

    图  3  PMT内爆过程

    Figure  3.  Process of the PMT implosion

    图  4  内爆数值计算模型

    Figure  4.  Implosion simulation model

    图  5  数值模拟与试验测点压力对比

    Figure  5.  Comparison of the simulation and test pressures of the measuring points

    图  6  数值模拟与试验测得的比冲量对比

    Figure  6.  Comparison of the simulation and test impulse of the measuring points

    图  7  水域流场变化过程

    Figure  7.  Evolution of the water field

    图  8  水流前锋速度变化过程

    Figure  8.  Evolution of the water front velocity

    图  9  PMT防护装置示意图

    Figure  9.  Schematic diagram of PMT protection device

    图  10  PMT防护装置等效模型

    Figure  10.  Equivalent models of the PMT protection device

    图  11  水域测点分布示意图

    Figure  11.  Distribution of measuring points in water field

    图  12  水域流场分布

    Figure  12.  Distribution of water field

    图  13  A1测点冲击波压力变化

    Figure  13.  Shock wave pressures varied with time at measuring point A1

    图  14  冲击波压力峰值分布

    Figure  14.  Distribution of peak pressure

    图  15  水流前锋速度

    Figure  15.  Velocities of the water fronts

    图  16  压力峰值与水流前锋速度变化

    Figure  16.  Variation of the peak pressure and the velocity of water front

    表  1  数值模拟与试验所得的测点压力峰值对比

    Table  1.   Difference between the simulation and test peak pressures of the measuring points

    测点压力峰值/MPa误差/%测点压力峰值/MPa误差/%
    试验数值模拟试验数值模拟
    114.1312.87 8.93 7.08 6.65 6.1
    2 7.68 6.6513.44 3.19 2.8510.7
    下载: 导出CSV

    表  2  数值模拟与试验所得的测点比冲量峰值对比

    Table  2.   Difference between the simulation and test peak impulse of the measuring points

    测点比冲量峰值/(kPa·s)误差/%测点比冲量峰值/(kPa·s)误差/%
    试验数值模拟试验数值模拟
    11.581.3514.630.870.83 4.6
    21.210.9124.840.320.2812.5
    下载: 导出CSV

    表  3  破碎面积具体值

    Table  3.   The value of break area

    αSb/m2αSb/m2
    1.00.7850.80.628
    0.90.7070.70.550
    下载: 导出CSV
  • [1] LING J J, MARY B, MILIND D, et al. Implosion chain reaction mitigation in underwater assemblies of photomultiplier tubes [J]. Nuclear Instruments and Methods in Physics Research, 2013, 729: 491–499. DOI: 10.1016/j.nima.2013.07.056.
    [2] IMAEDA H, SUN Mingyu. Dynamic characteristics of underwater objects after shock wave loading [C] // AIAA Aerospace Sciences Meeting. Kissimmee Florida, 2018. DOI: 10.2514/6.2018-0579.
    [3] SONG G, CHEN Z Y, LONG Y, et al. Experimental and numerical investigation of the centrifugal model for underwater explosion shock wave and bubble pulsation [J]. Ocean Engineering, 2017, 142: 523–531. DOI: 10.1016/j.oceaneng.2017.04.035.
    [4] NAVAL U. Underwater implosion of cylindrical metal tubes [J]. Journal of Applied Mechanics, 2013, 80: 1–11.
    [5] 包亦望. 脆性材料在双向应力下的断裂实验与理论分析 [J]. 力学学报, 1998, 30(6): 682–690. DOI: 10.3321/j.issn:0459-1879.1998.06.007.

    BAO Y W. Experiments and theoretic analysis for the fracture of brittle materials under biaxial stress [J]. Journal of Theoretical and Applied Mechanics, 1998, 30(6): 682–690. DOI: 10.3321/j.issn:0459-1879.1998.06.007.
    [6] YOSHIMURA M. Report on the Super-Kamiokande accident[R]. Institute of Space and Astronautical Science, 2001.
    [7] MILIND D, JEFFREY D, LING J J, et al. Underwater implosions of large format photo-multiplier tubes [J]. Nuclear Instruments and Methods in Physics Research, 2012, 670: 61–67. DOI: 10.1016/j.nima.2011.12.033.
    [8] GISH L A, WIERZBICKI T. Estimation of the underwater implosion pulse from cylindrical metal shells [J]. International Journal of Impact Engineering, 2015, 77: 166–175. DOI: 10.1016/j.ijimpeng.2014.11.018.
    [9] 杜志鹏, 杜俭业, 李营, 等. 不可压缩流体中球型容器内爆理论模型研究 [J]. 兵工学报, 2015, 36(S1): 92–96.

    DU Z P, DU J Y, LI Y, et al. An implosion theory for the spherical hollow vessel in the incompressible fluid [J]. Acta Armamentarii, 2015, 36(S1): 92–96.
    [10] 黄治新, 喻敏, 杜志鹏, 等. 水下中空结构物内爆试验方法研究 [J]. 振动与冲击, 2017, 36(3): 27–31. DOI: 10.13465/j.cnki.jvs.2017.03.005.

    HUANG Z X, YU M, DU Z P, et al. Implosion test method for underwater hollow structures [J]. Journal of Vibrtion and Shock, 2017, 36(3): 27–31. DOI: 10.13465/j.cnki.jvs.2017.03.005.
  • 期刊类型引用(16)

    1. 陈青华,陶彦光,梁振刚,杨佳岐. 双舱室模型内爆炸临舱压力载荷特性分析. 舰船科学技术. 2024(11): 75-79 . 百度学术
    2. 李营,杜志鹏,陈赶超,王诗平,侯海量,李晓彬,张攀,张伦平,孔祥韶,李海涛,郭君,姚术健,王志凯,殷彩玉. 舰艇爆炸毁伤与防护若干关键问题研究进展. 中国舰船研究. 2024(03): 3-60 . 百度学术
    3. 胡俊华,董奇,胡八一,任逸飞,黄广炎. 抗爆容器的内部爆炸效应和动态力学行为研究进展. 含能材料. 2024(09): 986-1008 . 百度学术
    4. 王昭,吴祖堂,杨军,李焰,刘文祥. 新型薄膜式压力传感器的参数设计. 爆炸与冲击. 2023(07): 174-184 . 本站查看
    5. 冯蕴雯,林心怡,薛小锋,杨祥,刘佳奇. 高可靠单向爆破的民机防爆结构设计. 航空学报. 2023(18): 191-205 . 百度学术
    6. 刘文祥,张德志,程帅,马艳军. 爆炸容器研究进展. 现代应用物理. 2023(03): 61-69 . 百度学术
    7. 刘欣,顾文彬,蔡星会,王涛,刘建青,王振雄,沈慧铭. 圆柱形爆炸容器的内壁爆炸载荷. 爆炸与冲击. 2022(02): 19-30 . 本站查看
    8. 刘正,聂建新,徐星,朱英中,刘攀,郭学永,闫石,张韬. 密闭空间内六硝基六氮杂异伍兹烷基复合炸药能量释放特性. 兵工学报. 2022(03): 503-512 . 百度学术
    9. 张鹏宙,董奇,杨沙. 爆炸载荷特征参数对无限长圆柱壳弹性动态响应的影响. 爆炸与冲击. 2021(06): 48-57 . 本站查看
    10. 孙琦,董奇,杨沙,张刘成. 内爆炸准静态压力对球形容器弹塑性动态响应的影响. 含能材料. 2020(01): 25-31 . 百度学术
    11. 张翔,崔春生,刘双峰,刘嘉颖,崔光强. 子弹在爆炸容器内的爆燃压力测试技术研究. 中国测试. 2020(01): 39-43 . 百度学术
    12. 张龙,邹虹,张宝国,张继军,张东亮,孔德骞. 有限空间爆炸静态压力的温度补偿方法. 爆炸与冲击. 2020(03): 100-109 . 本站查看
    13. 徐景林,顾文彬,刘建青,王振雄,陆鸣,徐博奥. 圆柱形爆炸容器内爆炸载荷的分布规律. 振动与冲击. 2020(18): 276-282 . 百度学术
    14. 夏彬伟,高玉刚,刘承伟,欧昌楠,彭子烨,刘浪. 缝槽水压爆破破岩载荷实验研究. 工程科学学报. 2020(09): 1130-1138 . 百度学术
    15. 张明明,张连生,王鑫. TNT内爆准静态压力实验和数值模拟研究. 兵器装备工程学报. 2019(05): 195-199 . 百度学术
    16. 孙琦,董奇,杨沙,张刘成. 内爆炸准静态压力对球形容器弹性动态响应的影响. 含能材料. 2019(08): 698-707 . 百度学术

    其他类型引用(10)

  • 加载中
图(16) / 表(3)
计量
  • 文章访问数:  4014
  • HTML全文浏览量:  1300
  • PDF下载量:  41
  • 被引次数: 26
出版历程
  • 收稿日期:  2019-11-18
  • 修回日期:  2020-02-11
  • 网络出版日期:  2020-07-25
  • 刊出日期:  2020-08-01

目录

/

返回文章
返回