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

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

孟令存 闫明 杜志鹏 张磊

引用本文:
Citation:

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

    作者简介: 孟令存(1994— ),男,硕士研究生,E-mail:1605733924@qq.com;
    通讯作者: 闫明, yanming7802@163.com
  • 中图分类号: O381;U663.85

Research of Underwater Implosion Mechanism of PMT Area Reduction Equivalent Model

    Corresponding author: YAN Ming, yanming7802@163.com ;
  • CLC number: O381;U663.85

  • 摘要: 光电倍增管(photomultiplier tube,PMT)是中微子探测器的核心部件,是由玻璃材料制成的内部真空的薄壳结构,排列在深水中工作,若一个PMT被压溃会产生内爆冲击波,会引起周围PMT发生殉爆。针对PMT内爆,建立了PMT内爆数值计算简化模型,并将计算与试验结果进行对比,验证简化模型的合理性。在此基础上,提出了基于面积折减等效模型的PMT内爆计算方法,通过等效模型分析了防护装置破口面积对PMT内爆的影响,得出随着防护装置破口面积的减小,水流碰撞发生PMT内爆的时刻相应提前,内爆产生的冲击波脉宽基本保持不变,冲击波压力峰值明显减小。该研究有利于找到有效的PMT内爆防护方法。
  • 图 1  光电倍增管实物图

    Figure 1.  The picture of photomultiplier

    图 2  内爆试验装置示意图

    Figure 2.  The picture of tester for implosion

    图 3  PMT内爆过程

    Figure 3.  The process of the PMT implosion

    图 4  内爆数值计算模型

    Figure 4.  Implosion simulation model

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

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

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

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

    图 7  水域流场变化过程

    Figure 7.  Evolution of the water field

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

    Figure 8.  Evolution of the waters’ front velocity

    图 9  PMT防护装置示意图

    Figure 9.  Schematic diagram of PMT protection device

    图 10  PMT防护装置等效模型

    Figure 10.  Equivalent model of the PMT protection device

    图 11  水域测点分布示意图

    Figure 11.  The distribution of measuring points in water field

    图 12  水域流场分布

    Figure 12.  The distribution of water field

    图 13  A1测点冲击波压力峰值

    Figure 13.  The peak pressure of shock wave on measuring point A1

    图 14  冲击波压力峰值分布

    Figure 14.  The distribution of peak pressure

    图 15  水流前锋速度

    Figure 15.  The velocity of the waters front

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

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

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

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

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

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

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

    测点压力峰值/(MPa·ms)误差/%
    试验数值模拟
    11.581.3514.6
    21.210.9124.8
    30.870.83 4.6
    40.320.2812.5
    下载: 导出CSV

    表 3  破碎面积具体值

    Table 3.  The value of break area

    折减系数α10.90.80.7
    破口表面积Sb/m20.7850.7070.6280.55
    下载: 导出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 Zuyu, LONG Yuan, 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 material sunder 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] 吉村太彦. 超级神冈事故原因究明等委员会报告[R]. 宇宙科学研究所, 2001年.
    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.
  • [1] 郝鹏程冯其京胡晓棉 . 内爆加载金属界面不稳定性的数值分析. 爆炸与冲击, doi: 10.11883/1001-1455(2016)06-0739-06
    [2] 庞伟宾李永池何翔 . 化爆冲击波在T型通道内到时规律的实验研究. 爆炸与冲击, doi: 10.11883/1001-1455(2007)01-0063-05
    [3] 余尚江李科杰 . 混凝土结构内冲击波应力传感器设计及其行为. 爆炸与冲击, doi: 10.11883/1001-1455(2005)04-0350-05
    [4] 王启睿张晓忠孔福利张福明 . 多级穿廊结构坑道口部内爆炸冲击波传播规律的实验研究. 爆炸与冲击, doi: 10.11883/1001-1455(2011)05-0449-06
    [5] 梁化彬胡春明徐胜利王昌建费立森 . 爆燃波在含内构件管道中传播现象的实验研究. 爆炸与冲击, doi: 10.11883/1001-1455(2007)06-0487-06
    [6] 高康华李斌刘宇都孙松 . 外爆条件下冲击波与组合壳结构的相互作用. 爆炸与冲击, doi: 10.11883/bzycj-2017-0370
    [7] 赵新颖王伯良李席 . 温压炸药在野外近地空爆中的冲击波规律. 爆炸与冲击, doi: 10.11883/1001-1455(2016)01-0038-05
    [8] 薛冰马宏昊沈兆武余勇 . 爆炸容器内小药量实验动态标定压力传感器. 爆炸与冲击, doi: 10.11883/1001-1455(2015)03-0437-05
    [9] 李雪梅金孝刚李大红 . 钢圆管在内爆加载下的层裂特性研究. 爆炸与冲击, doi: 10.11883/1001-1455(2005)02-0107-05
    [10] 孙可明辛利伟吴迪 . 超临界CO2气爆煤体致裂机理实验研究. 爆炸与冲击, doi: 10.11883/bzycj-2016-0230
    [11] 姚成宝李若田宙郭永辉 . 空气自由场中强爆炸冲击波传播二维数值模拟. 爆炸与冲击, doi: 10.11883/1001-1455(2015)04-0585-06
    [12] 周杰陶钢王健 . 爆炸冲击波对肺损伤的数值模拟. 爆炸与冲击, doi: 10.11883/1001-1455(2012)04-0418-05
    [13] 师华强宗智贾敬蓓 . 水下爆炸冲击波的近场特性. 爆炸与冲击, doi: 10.11883/1001-1455(2009)02-0125-06
    [14] 张柱晋艳娟 . 反向起爆模型下的冲击波加载. 爆炸与冲击, doi: 10.11883/1001-1455(2014)02-0223-06
    [15] 刘贵兵侯海量朱锡张国栋 . 液滴对爆炸冲击波的衰减作用. 爆炸与冲击, doi: 10.11883/1001-1455(2017)05-0844-09
    [16] 李梅蒋建伟王昕 . 复合装药空气中爆炸冲击波传播特性. 爆炸与冲击, doi: 10.11883/bzycj-2016-0209
    [17] 许志宇谭永华李小明 . 冲击波的小波数值计算方法. 爆炸与冲击, doi: 10.11883/bzycj-2018-0467
    [18] 周佩杰王坚陶钢周杰 . 泡沫材料对冲击波的衰减特性. 爆炸与冲击, doi: 10.11883/1001-1455(2015)05-0675-07
    [19] 姚春德续晗姚安仁唐超 . 缸内激波对锥顶型燃烧室的活塞破坏机理. 爆炸与冲击, doi: 10.11883/1001-1455(2015)01-0057-08
    [20] 顾文彬郑向平刘建青李丹俊陆鸣 . 浅层水中爆炸冲击波对混凝土墩斜碰撞作用试验研究. 爆炸与冲击, doi: 10.11883/1001-1455(2006)04-0361-06
  • 加载中
图(16)表(3)
计量
  • 文章访问数:  53
  • HTML全文浏览量:  43
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-11-18
  • 录用日期:  2020-02-11
  • 网络出版日期:  2020-07-25

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

    作者简介:孟令存(1994— ),男,硕士研究生,E-mail:1605733924@qq.com
    通讯作者: 闫明, yanming7802@163.com
  • 1. 沈阳工业大学机械工程学院,辽宁 沈阳 110870
  • 2. 海军研究院,北京 100161

摘要: 光电倍增管(photomultiplier tube,PMT)是中微子探测器的核心部件,是由玻璃材料制成的内部真空的薄壳结构,排列在深水中工作,若一个PMT被压溃会产生内爆冲击波,会引起周围PMT发生殉爆。针对PMT内爆,建立了PMT内爆数值计算简化模型,并将计算与试验结果进行对比,验证简化模型的合理性。在此基础上,提出了基于面积折减等效模型的PMT内爆计算方法,通过等效模型分析了防护装置破口面积对PMT内爆的影响,得出随着防护装置破口面积的减小,水流碰撞发生PMT内爆的时刻相应提前,内爆产生的冲击波脉宽基本保持不变,冲击波压力峰值明显减小。该研究有利于找到有效的PMT内爆防护方法。

English Abstract

  • 在深水环境下工作的光电倍增管(photomultiplier tube,PMT),其玻璃外表面承受高静水压力,当压力载荷超过玻璃壳体的结构强度时,PMT被压溃使水流向中心流动,流体的水下压力势能转化为动能,当水流压缩至最小限度时,高速流体碰撞产生水锤型冲击,发生内爆[1-4],产生内爆冲击波。水下内爆冲击波会对周围的PMT造成破坏,引发相邻的PMT发生殉爆。2014年,美国深海科研潜艇“Nereus”号在执行深海作业时,潜艇的中空陶瓷浮力装置在巨大静水压力下被压溃,发生内爆,使其遇难。在高能物理研究领域,脆性材料[5]制成的内部真空的PMT被用来进行粒子探测,工作在深水环境中,容易发生内爆;2001年,日本超级神冈中微子实验站发生PMT殉爆事故,此次事故炸毁将近8000个PMT,直接经济损失高达3000万美元,中微子研究遭受重大挫折[6]。我国将建设世界领先地位的江门中微子实验站,该实验站拟采用16000只直径约为500 mm的PMT作为探测器,其最大工作水深为40 m。根据国外相关事故的经验教训和江门中微子探测器的实际工作要求,为防止PMT殉爆的发生,探索PMT内爆机理具有重要的应用价值和实际意义。

    国内外对水下内爆进行了诸多研究。Diwan等[7]在大型可加压容器罐内进行了2次0.69 MPa静水压下PMT内爆试验,得到了PMT内爆过程高速摄像及内爆冲击波压力时域曲线,并用LS-DYNA对其进行了数值模拟,但计算结果与试验结果差别较大;Gish等[8]在可加压容器罐中进行了金属圆柱壳内爆试验,并结合数值模拟研究了圆柱壳长度、直径和厚度对内爆冲击波的影响;杜志鹏等[9]将水下爆炸气泡动力学与水下内爆相结合,基于能量守恒关系,推导出不可压缩流体中球形容器内爆理论模型,并利用该理论模型分析了容器尺寸、静水压力对内爆冲击波压力峰值、气泡溃灭时间周期的影响规律;黄治新等[10]基于应力波原理,提出一种在大型钢制压力罐内进行PMT内爆试验的试验方法,并通过此方法成功测得PMT内爆发生过程高速摄像及内爆冲击波。

    上述研究主要针对水下内爆特性进行了试验与数值模拟,但对于PMT内爆并没有有效的防护。本文开展了PMT水下内爆试验,用于验证PMT内爆简化模型模拟水下内爆的准确性;并通过有限元计算的方式探究水流涌入面积对内爆冲击波的影响,提出了当PMT发生破碎时,减小水流涌入面积能有效减小内爆冲击波强度的内爆防护原理。

    • PMT实物如图1所示,其形状如灯泡,PMT玻璃外壳厚度为5 mm,头部球体部分外径为0.508 m,尾部直径为0.1 m,总高度约为0.7 m,内部真空体积为6.992×10−2 m3;PMT内爆试验在如图2所示直径为3米的可加压密封钢制容器罐内进行,压力罐侧面透明视窗用于放置照明和高速摄影设备,罐内设有工作平台、液压式挤压装置。试验中共布置4个PCB动压传感器,F1位于PMT正上方,距PMT中心0.41 m处;F2位于PMT赤道平面,靠近压力罐入口处,距PMT中心0.55 m处;F3、F4分别位于PMT赤道平面,远离压力罐入口处,距PMT中心0.55、1.35 m处。试验前,罐中装入水,罐上方预留出适量空气,考虑PMT殉爆防护安全余量,利用空压机对压力罐中施加0.5 MPa压力模拟PMT深水工作环境;试验时,挤压装置模拟外界恶劣环境挤压PMT发生内爆,内爆发生过程通过高速摄像(3 000 Hz)进行记录,压力数据利用杭州亿恒数据采集系统(采样率为1 MHz)进行采集,试验成功测得PMT内爆过程及4组内爆冲击波压力数据。

      图  1  光电倍增管实物图

      Figure 1.  The picture of photomultiplier

      图  2  内爆试验装置示意图

      Figure 2.  The picture of tester for implosion

      完成PMT内爆试验后,对PMT内爆过程进行分析:由内爆试验高速摄影可知整个内爆过程如图3所示,可将其分为4个典型阶段:(1)挤压过程阶段;(2)裂纹传播阶段;(3)整体压溃形成冲击波阶段;(4)冲击波传播与碎片飞散阶段。挤压装置挤压PMT产生局部初始裂纹,初始裂纹瞬间传遍整个玻壳而发生坍塌,外部高压水流迅速内涌发生碰撞产生内爆冲击波。通过高速摄像充分展现了PMT内爆发生过程和冲击波产生过程,对水下内爆数值模拟具有一定的指导意义。

      图  3  PMT内爆过程

      Figure 3.  The process of the PMT implosion

      下面用有限元方法模拟PMT内爆过程,由高速摄影看到PMT产生局部初始裂纹后会瞬间传遍整个玻壳而发生坍塌,故进行简化数值模拟时,忽略PMT玻壳的影响,只单独计算外部水域流场,按体积等效的方式将PMT等效为半径为0.25 m的真空球体。图4为PMT内爆有限元模型中部截面图,整个水域半径为1.7 m,采用六面体欧拉网格,单元类型为EC3D8R,水域节点共计619 393个,单元共计613 872个;蓝色部分欧拉网格赋予水材料属性,水采用Gruneisen状态方程描述,密度为1 000 kg/m3us-up直线截距(声速)为1 484 m/s,;红色部分欧拉网格设置无材料填充,默认为真空状态,真空半径为0.25 m,用于模拟PMT。利用初始应力的方法设置蓝色部分水域初始静水压力为0.5 MPa,水域外边界设置冲击波无反射、流体自由流入、流出边界条件,模型中压力测点S1~S4与内爆试验测点F1~F4位置相对应,最后通过显示动态分析步(explicit)进行求解。将简化数值模型的计算结果与试验结果相对比,验证简化模型的合理性。

      图  4  内爆数值计算模型

      Figure 4.  Implosion simulation model

    • 以冲击波最大峰值时刻为基准,将数值计算所得各测点的冲击波压力数据与对应试验测点压力数据绘制曲线如图5所示。试验中各测点压力数据首先产生内爆冲击波峰值,最大峰值过后经过2 ms,出现第2个压力峰值,峰值约为2 MPa,这是由于压力罐壁反射冲击波导致的;由于压力罐直径为3 m,按冲击波传播速度为1 500 m/s计算,罐壁反射冲击波应经过2 ms返回至压力测点,这与试验所测得的数据相吻合,证明了冲击波测量的准确性。数值模型中冲击波压力数据与试验数据相比,波形较为平滑,振荡较少,脉宽偏小;最大压力峰值过后无反射冲击波峰值产生,这是由于数值模型设置了无反射边界条件,故所得压力数据无反射波形。将数值模型中各测点冲击波压力峰值与相对应试验压力峰值列于表1,计算值与试验值最大相差13.4%。

      测点压力峰值/MPa误差/%
      试验数值模拟
      114.1312.87 8.9
      2 7.68 6.6513.4
      3 7.08 6.65 6.1
      4 3.19 2.8510.7

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

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

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

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

    • 通过对比数值模拟与试验压力冲击波数据可知,与计算相比试验压力数据曲线存在反射冲击波峰值。为消除试验压力数据反射冲击波对比冲量的影响,提高可比性,对计算与试验数据,均取内爆产生的最大波峰及其之前的压力数据进行积分,得到如图6所示的试验与计算内爆冲击波比冲量对比曲线。试验数据比冲量在积分区域内均呈逐渐增大至平稳的趋势,计算数据比冲量出现负值是由于压力数据存在负压导致的;取计算与试验内爆冲击波比冲量数据最大值与最小值的差值作为内爆冲击波比冲量峰值,将试验与计算数据比冲量峰值列于表2,可看出测点1、2通过计算得到的冲击波比冲量峰值与试验偏差较大,这是由于测点1、2试验压力数据最大峰值后出现较大振荡干扰导致的。通过以上分析可知,与内爆试验相比,内爆数值计算所得冲击波压力峰值最大相差13.4%,冲击波比冲量峰值最大相差24.8%,这对内爆数值计算来说,是可以接受的,满足实际工程需求。

      测点压力峰值/(MPa·ms)误差/%
      试验数值模拟
      11.581.3514.6
      21.210.9124.8
      30.870.83 4.6
      40.320.2812.5

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

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

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

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

    • 下面对内爆简化模型的水域流场与流速进行分析。欧拉网格具有与材料完全隔离的特性,欧拉网格中材料的体积分数以不同的颜色区分,红色代表材料体积分数为1,即水填满整个网格;蓝色代表材料体积分数为0,即没有水填充。内爆简化模型计算中,水域流场变化过程如图7所示。图7(a)为内爆初始状态,蓝色代表PMT真空区域,红色代表外部水域,其初始静水压为0.5 MPa;图7(b)(c)为PMT破碎后,外部水流迅速内涌发生内爆的过程,可看到水流从四周均匀涌向球心;图7(d)为内涌水流碰撞产生冲击波的瞬间。PMT内爆简化模型模拟内爆发生过程共经过10.25 ms,与试验高速摄像拍摄的PMT内爆过程及发生内爆的时间相近。数值计算中水流前锋速度变化过程如图8所示,在整个过程中水流前锋速度逐渐增加,内涌高速水流碰撞前速度达324.9 m/s,说明PMT内爆发生前碰撞水流具有很大的动能。

      图  7  水域流场变化过程

      Figure 7.  Evolution of the water field

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

      Figure 8.  Evolution of the waters’ front velocity

      通过以上PMT内爆数值模拟与试验对比可知,PMT内爆数值简化模型在冲击波压力、比冲量、水域流场方面,能够较为全面、准确地模拟PMT内爆真实物理情形。

    • 在PMT水下内爆数值模拟方法正确建立的基础上,为避免PMT殉爆的发生,下面利用数值模拟对内爆冲击波的影响因素进行研究。PMT发生内爆时,玻壳整体破碎向内坍塌,外部水流从四周涌入,发生碰撞产生冲击波。为探究水流涌入面积对内爆的影响,假设PMT瞬间破碎后,外部存在相同直径的有机玻璃防护装置,如图9所示,其中R为等效半径,H为球冠高度,D为防护装置破口直径。

      图  9  PMT防护装置示意图

      Figure 9.  Schematic diagram of PMT protection device

      定义无量纲面积折减系数α来表示防护装置破口程度,如公式(1)所示,破口面积Sb具体数值列于表3

      折减系数α10.90.80.7
      破口表面积Sb/m20.7850.7070.6280.55

      表 3  破碎面积具体值

      Table 3.  The value of break area

      $\alpha {\rm{ = }}\frac{{{S_{\rm{b}} }}}{{{S_{\rm{a}} }}}$

      式中:Sa=4πR2,为PMT等效表面积;Sb=4πRH,为防护装置破口表面积;

    • 当面积折减系数α为0.9、0.8、0.7时,需在内爆简化模型(α=1)中加入防护装置模型,如图10所示,此时应采用欧拉-拉格朗日耦合分析方法,有机玻璃防护装置利用拉格朗日网格描述,采用厚度10 mm的壳单元。其密度为1 180 kg/m3,弹性模量为2.77 GPa,泊松比为0.376,并约束其所有节点自由度(防护装置保持固定),其余初始、边界条件与内爆简化模型相同。

      图  10  PMT防护装置等效模型

      Figure 10.  Equivalent model of the PMT protection device

      为探究不同面积折减系数下冲击波压力在水中传播情况,在PMT防护装置一侧缺口对应的水域轴线上布置一系列压力测点,如图11所示,在测点A1至A29连线上等间距布置29个压力测点(A1~A29),其中A1为PMT防护装置圆周处测点。

      图  11  水域测点分布示意图

      Figure 11.  The distribution of measuring points in water field

    • (1)水域流场分布

      不同面积折减系数下高速水流碰撞瞬间流场分布如图12所示,同理,以颜色区分欧拉网格中水的体积分数。由水域流场分布可直观地看到,与PMT整体破碎内涌水流四周碰撞相比,随着防护装置破口面积的减小,真空区域内涌水流碰撞量随之大幅度减少,只是水流前锋局部碰撞产生冲击波,这可能导致冲击波强度发生变化。

      图  12  水域流场分布

      Figure 12.  The distribution of water field

      (2)内爆冲击波压力分布及强度分析

      下面对不同面积折减系数下冲击波强度进行分析:图13为各面积折减系数下A1测点处的冲击波压力时域曲线。可看出随着防护装置破口面积减小,PMT圆周处测点的冲击波压力峰值迅速减弱。但冲击波脉宽基本不变;同时,随着防护装置破口面积的减小,水流碰撞发生内爆的时刻也相应提前。进一步对冲击波压力峰值进行分析,取各面积折减下测点A1~A29冲击波压力峰值,以冲击波压力峰值为纵坐标,测点距球心的距离为横坐标,绘制曲线如图14所示,为了曲线表达清晰,每间隔一个测点进行标记显示。以PMT球心为内爆中心,可明显观察到内爆压力峰值随与测点距离的增大而逐渐衰减,在0.25~0.4 m范围内,内爆压力峰值迅速衰减,0.4 m以外其内爆压力峰值衰减相对较慢,在1.6 m处内爆压力平稳回归到约2 MPa;最重要的是随着面积折减系数的减小,防护装置破口面积减小,同一测点水流内涌碰撞产生的冲击波压力峰值均大幅度地减弱,这可能与内涌水流碰撞量的减少有关,但也可能与水流前锋碰撞速度的变化有关。

      图  13  A1测点冲击波压力峰值

      Figure 13.  The peak pressure of shock wave on measuring point A1

      图  14  冲击波压力峰值分布

      Figure 14.  The distribution of peak pressure

      下面对内爆冲击波的影响因素进行进一步分析:图15为各面积折减系数下内爆发生前的水流前锋速度随时间变化曲线,可看出水流前锋速度在前期增长缓慢,后期迅速增加;此外,防护装置破口面积越小,水流前锋速度增加越快,但其碰撞前速度相差不明显,说明防护装置破口面积对内涌水流前锋最终碰撞速度值影响较小。各面积折减系数下,A1测点压力峰值与内涌水流前锋速度变化如图16所示。在防护装置破口面积变化30%情况下,A1测点内爆压力峰值变化73.6%,而对应的内涌水流前锋速度变化仅为13.9%,说明随着防护装置破口面积的减小,内爆冲击波压力峰值大幅度减弱主要是内涌水流碰撞量减少导致的,而内涌水流前锋最终碰撞速度对其影响较小。故以此为减小PMT内爆冲击波压力峰值的内爆机理,为PMT内爆防护提供重要思路。

      图  15  水流前锋速度

      Figure 15.  The velocity of the waters front

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

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

    • (1)以PMT内爆试验为基础,通过数值模拟再现了PMT内爆过程,对比数值模拟与试验结果可知:数值模拟得出的PMT内爆发生过程、冲击波压力峰值和冲击波压力比冲量与试验吻合较好,验证了PMT内爆简化模型的合理性和适用性,提出了一套PMT内爆数值模拟方法。

      (2)在此基础之上,提出了基于面积折减等效模型的PMT内爆计算方法,通过等效模型分析了不同面积折减系数对PMT内爆冲击波的影响,得出基于面积折减等效模型的PMT水下内爆机理:随着防护装置破口面积的减小,水流碰撞引起的PMT内爆发生时刻相应提前,内爆产生的冲击波脉宽基本保持不变、冲击波峰值明显减小,这主要是随着防护装置破口面积的减小使内涌水流碰撞量减少导致的,而内涌水流前锋最终碰撞速度对其影响较小。

参考文献 (10)

目录

    /

    返回文章
    返回