供风量与气泡帷幕层数协同下水中爆炸冲击波的削波效果

杜明燃 陈智凡 陆少锋 梁进 李基锐 王尹军 王天照 陈宇航

牛林耕, 闫栋, 王根伟, 宋辉, 郭美卿. 类向日葵夹芯圆柱壳的准静态轴向吸能特性与优化[J]. 高压物理学报, 2023, 37(4): 044206. doi: 10.11858/gywlxb.20230637
引用本文: 杜明燃, 陈智凡, 陆少锋, 梁进, 李基锐, 王尹军, 王天照, 陈宇航. 供风量与气泡帷幕层数协同下水中爆炸冲击波的削波效果[J]. 高压物理学报, 2024, 38(1): 015103. doi: 10.11858/gywlxb.20230705
NIU Lingeng, YAN Dong, WANG Genwei, SONG Hui, GUO Meiqing. Quasi-Static Axial Energy Absorption Characteristics and Optimization of Sunflower-Like Sandwich Cylindrical Shells[J]. Chinese Journal of High Pressure Physics, 2023, 37(4): 044206. doi: 10.11858/gywlxb.20230637
Citation: DU Mingran, CHEN Zhifan, LU Shaofeng, LIANG Jin, LI Jirui, WANG Yinjun, WANG Tianzhao, CHEN Yuhang. Synergistic Effect of Air Supply Volume and Bubble Curtain Layer on the Shock Wave Attenuation of Underwater Explosion[J]. Chinese Journal of High Pressure Physics, 2024, 38(1): 015103. doi: 10.11858/gywlxb.20230705

供风量与气泡帷幕层数协同下水中爆炸冲击波的削波效果

doi: 10.11858/gywlxb.20230705
基金项目: 安徽省高校科学研究项目(KJ2021A0431);广西重点研发计划(桂科AB22035001);安徽省自然科学基金(1908085QA33);安徽理工大学引进人才科研启动基金(11881);广西重点研发计划(防科AB21014001)
详细信息
    作者简介:

    杜明燃(1987-),男,博士,副教授,主要从事炸药性能和爆破技术研究. E-mail:dumingranaust@163.com

    通讯作者:

    陈智凡(1999-),男,硕士研究生,主要从事水中爆炸研究. E-mail:786923403@qq.com

  • 中图分类号: O389; TD235

Synergistic Effect of Air Supply Volume and Bubble Curtain Layer on the Shock Wave Attenuation of Underwater Explosion

  • 摘要: 气泡帷幕能有效地削弱水中冲击波对周围环境的影响。为了研究气泡帷幕供风量和层数对水中冲击波的协同作用,在供风量为30、60、90 L/min的条件下分别设计了含1、2、3层气泡帷幕的水下爆炸试验。结果表明,气泡帷幕的衰减率随供风量和层数增加而增大。当供风量较小(如30、60 L/min)时,随着气泡帷幕层数的增加,相邻层数之间峰值压力的衰减效率越来越低;当供风量较大(如90 L/min)时,随着气泡帷幕层数的增加,相邻层数之间峰值压力的衰减效率越来越高。结合实际工程的经济效益和水下复杂环境问题对削波效果进行分析,确定在供风量为30 L/min时开启2层气泡帷幕是最优的削波方案,为相关的实际工程问题提供参考和新思路。

     

  • 超硬材料在基础科学和工业应用中具有重要的意义,其设计与合成受到人们的广泛关注[1-2]。由轻元素(B、C、N和O)形成的强共价键化合物是潜在超硬材料的典型来源,近年来取得了丰硕的研究成果[3],并成功合成出B-C[4-5]、B-O[6-7]、B-C-N[8-10]等超硬材料。超硬材料的理论设计和预测也取得了巨大成就,提出了由一元单质[11]、二元[12]及三元[13]化合物组成的新型超硬材料,标志性工作是能够匹敌金刚石硬度的超硬材料β-C3N4的预测[14]

    自1997年采用高压技术通过沃克型多砧装置合成B-C-O化合物(B6C1.1O0.33和B6C1.28O0.31)以来[15],三元化合物B-C-O进入了人们的研究视野。2001年,Bolotina等[16]在5.5 GPa、1400 K的条件下使1∶1混合的B4C和B2O3发生反应,成功制备出B(C,O)0.1555。实验合成的三元B-C-O化合物证实了其存在的物质基础,自此B-C-O体系化合物作为超硬材料候选物质来源引起了广泛的关注。

    以B2CO(具有与金刚石等电子体的最简配比B-C-O化合物)为例,Li等[17]开展了对B-C-O化合物的开创性理论研究,预测出两种四方晶系超硬结构,分别为tP4-B2CO和tI16-B2CO。实际上,在B-C-O三元体系中存在很多金刚石等电子体化合物,如B2CxO(x = 2, 3, 4, ···),随着C含量的增加,B-C-O三元化合物结构中存在C的堆积,进而形成更短更强的sp3杂化C―C共价键。因此,Li等[17]认为B2CxO(x = 2, 3, 4, ···)化合物中极有可能存在比B2CO更硬的结构。随后,Zhang等[18]针对B2CxO化合物开展了研究,通过系统的粒子群优化算法结构搜索,提出了3种四方晶系的超硬B2CxO(x ≥ 2)相,分别为B2C2O、B2C3O和B2C5O,并发现随着碳含量的增加,sp3杂化C―C共价键有利于弹性模量和硬度的提高[18],与Li等的观点吻合。2017年,Liu等[19]从理论上预测出一种类蓝丝黛尔石lonsdalite的斜方晶系B2CO结构(oP8-B2CO),这是首次提出非四方结构的B2CO化合物。基于键阻模型的微观硬度理论分析,预测oP8-B2CO具有高达47.70 GPa的硬度。随后Qiao等[20]系统研究了高压下oP8-、tI16-和tP4-B2CO化合物的结构和电子性质,基于杂化泛函的研究表明,oP8-、tI16-和tP4-B2CO都是宽的间接带隙半导体材料,其带隙均随压力的增大而增大。

    受到近年来所预测的B-C-O结构特征的启发(如tP4-B2CO和tI16-B2CO均具有类金刚石结构[17]oP8-B2CO具有类蓝丝黛尔石结构[19],B2C2O、B2C3O和B2C5O均具有类金刚石结构),可由不同的tP4-B2CO超胞通过部分B和O被C取代产生[18]。基于结构对性质的影响,Liu等[21]以具有超硬特性的碳同素异形体为模型,通过多步骤手动构建结合第一性原理研究,预测出具有超硬性能的两种正交晶系B2CO相oP16-B2CO和oC16-B2CO,它们分别与Cco-C8和Bct-C4具有相似的结构。

    对于B2CO三元化合物,虽然已经成功预测出一些结构,但是考虑到轻元素的灵活性和报道的宽稳定压力范围,寻找新型未知的热力学稳态相仍然值得关注[22]。为此,Yan等[22]采用粒子群优化晶体结构分析软件CALYPSO,探索了B2CO在环境条件下的热力学稳定结构,发现了一种sp2-sp3杂化共存的新型高硬度正交晶系oI16-B2CO,比所有先前提出的B2CO结构在能量上更稳定。受到Yan等研究成果的鼓舞,Chen等[23]提出了一种四方晶系sp2-sp3 B―C键和B―O键共存的亚稳态B2CO化合物tP16-B2CO。研究表明,随着压力的增加,tP16-B2CO的带隙先增大后减小。通过计算tP16-B2CO的硬度、应力-应变关系等力学性能,发现其为硬质材料,硬度为23.19 GPa,且最大应力沿[001]方向的各向异性远高于沿[100]方向[23]

    此外,通过长期探索,非金刚石等电子体的B-C-O化合物研究取得了较大进展。2016年,Wang等[24]通过从头算变组成演化模拟,探究了0~50 GPa压力范围内的B-C-O体系,发现了一种新型四方相B4CO4。B4CO4在高于23 GPa的压力下是热力学稳定的,但在室温条件下保持亚稳态,计算显示其硬度为38~41 GPa,表明B4CO4具有超硬的潜力。Nuruzzaman等[25]研究发现,B4CO4具有较大的体积模量、剪切模量、杨氏模量、维氏硬度以及较小的泊松比,表明B4CO4可能是超硬材料的候选材料。鉴于直观上B4CO4结构中心存在较大的空心空间,B4CO4是否是一种超硬材料还需要进一步澄清[26]。Zheng等[26]采用第一性原理计算方法,详细推导了B4CO4的杨氏模量和剪切模量的方向依赖关系,结果表明B4CO4具有高度的各向异性。基于应变-应力法,计算了B4CO4沿主晶体方向的理想拉伸强度和剪切强度,研究表明,B4CO4的塑性变形主要是沿(001)[100]滑移体系的剪切模式,这与B―O键的断裂有关。最弱的理想剪切强度为27.5 GPa,证明了B4CO4的高硬度属性。

    Wang等[24]还报道了B-C-O体系中两个低焓亚稳相B6C2O5和B2CO2。基于计算的声子色散和弹性常数,B6C2O5和B2CO2在环境条件下都满足动力学和弹性力学稳定性。通过计算证实了B6C2O5和B2CO2具有高硬属性,硬度值分别为36和40 GPa[24]。Qiao等[27]计算了B2CO2和B6C2O5的能带结构,发现B2CO2和B6C2O5是宽带隙半导体材料,其间接能隙分别为5.66和5.24 eV。弹性各向异性结果表明,B6C2O5的剪切模量具有较大的各向异性,而B2CO2的杨氏模量具有较大的各向异性。此外,Lyakhov-Oganov模型的硬度计算表明B2CO2是一种潜在的超硬材料。随后,Liu等[28-29]用第一性原理计算预测了一种四方晶系B-C-O化合物tI12-B6C4O2tI12-B6C4O2具有较大的力学模量和较高的硬度,是一种硬度为21.9 GPa的高硬材料。各向异性研究表明,沿特定方向的杨氏模量存在E[110] < E[100] < E[001] < E[111] < E[011]的关系。

    纵观三元B-C-O化合物体系研究发现,虽然取得了丰富的成果,但是缺乏组分对性能调控机理的研究,尤其对于非金刚石等电子体系列的B-C-O化合物,其组分多样,相关的结构报道不多,研究不够系统,尚未揭示组分对其性能的作用规律。为此,本研究开展了非金刚石等电子体系列的B-C-O化合物结构预测,提出一种新型B4C6O4化合物,并基于B4C6O4、B2CO2和B4CO4,系统研究BCxO化合物的性质与组分之间的关系。

    通过基于粒子群优化算法的晶体结构预测软件CALYPSO[30-32],在常压下通过变组分方式生成非金刚石等电子体B-C-O化合物的潜在结构,然后在CASTEP模块[33]中采用广义梯度近似的Perdewey-Burke-Ernzerhof(PBE)泛函[34]进行全弛豫结构优化及性质研究。其中超软赝势[35]的截止能量设为380 eV,并采用Koelling-Harmon方式进行相对论处理。通过Monkhorst-Pack网格划分产生K点,其中划分密度为2π × 0.04 Å−1。在几何优化过程中,电子步自洽精度为5 × 10−4 meV/atom,离子步收敛精度为每原子能量变化小于5 × 10−3 meV,指定最大应力的收敛阈值为0.02 GPa,指定最大位移小于5 × 10−4 Å。计算电子能带结构和态密度时涉及到的能带能量容差低至1 × 10−5 eV。声子计算采用有限位移法[36],计算弹性常数时每个应变产生9个模式,最大应变幅值为0.009,收敛精度为每原子能量变化不超过1 × 10−3 meV,最大力为0.002 eV/Å,最大位移不超过1 × 10−4 Å。

    此外,给出了计算过程中涉及的模型类型及其布里渊区高对称点信息。B4C6O4单胞:G(0, 0, 0)→Z(0, 0, 0.5)→T(−0.5, 0, 0.5)→Y(−0.5, 0, 0)→S(−0.5, 0.5, 0)→X(0, 0.5, 0)→U(0, 0.5, 0.5)→R(−0.5, 0.5, 0.5);B4C6O4原胞:G(0, 0, 0)→Z(0, 0, 0.5)→T(0.5, 0.5, 0.5)→Y(0.5, 0.5, 0)→G(0, 0, 0)→S(0, 0.5, 0)→R(0, 0.5, 0.5)→Z(0, 0, 0.5);B4CO4原胞:Z(−0.5, 0.5, 0.5)→G(0, 0, 0)→X(0, 0, 0.5)→P(0.25, 0.25, 0.25)→N(0, 0.5, 0)→G(0, 0, 0);B2CO2原胞:L(−0.5, 0, 0.5)→M(−0.5, −0.5, 0.5)→A(−0.5, 0, 0)→G(0, 0, 0)→Z(0, −0.5, 0.5)→V(0, 0, 0.5)。

    通过CALYPSO软件预测一种新型非金刚石等电子体的B-C-O化合物,常压下优化得到的结构模型如图1所示。该B-C-O化合物具有(0.5, 0.5, 0)中心对称的正交晶系结构,空间群Cmm2,晶胞参数abc分别为4.587、5.441和3.812 Å。晶胞中B、C、O的原子数分别为4、6、4,这里记为B4C6O4。常压下B4C6O4的密度为3.130 g/cm3,其中C原子有C1和C2两类,结构中各类原子占位情况如表1所示。

    图  1  常压下B4C6O4结构(2 × 2 × 2)的超胞模型
    Figure  1.  Structure graphs for B4C6O4 (2 × 2 × 2) supercell at ambient pressure
    表  1  常压下B4C6O4的原子坐标
    Table  1.  Atomic Wyckoff positions of B4C6O4 at ambient pressure
    AtomWyckoff sitexyz
    B4d0.2990.5000.267
    C12b00.5000.991
    C24c0.2500.2500.491
    O4e00.2790.746
    下载: 导出CSV 
    | 显示表格

    观察发现:结构中存在B原子之间的成键,B―B键长为1.842 Å;O原子与C1原子形成的O―C1键的键长为1.522 Å,该类键与x轴呈空间垂直关系;O原子与C2原子形成的O―C2键的键长为1.511 Å,该类键沿x轴相互串联。C2与B原子形成的键具有两类键长:一类在空间上呈现沿y轴无限连接,键长为1.622 Å;另一类键则在空间上垂直于y轴,键长1.728 Å,并与B―B键以等键数相间的方式连接,沿x轴延伸。整个结构中C原子以4配位方式与周边2个O原子和2个B原子成键,B原子与周边2个C2原子、1个C1原子和1个近邻B原子连接成键,而O原子只以3配位形式与周边3个C原子成键。

    B4C6O4的弹性稳定性通过计算其独立的弹性常数Cij,并经波恩弹性稳定性判据验证。B4C6O4的晶体结构属于正交晶系,其独立的弹性常数有9个,波恩弹性稳定性判据[37]

    Cii>0(i=1,4,5,6),C11C22>C212,C11C22C33+2C12C13C23>C11C223+C22C213+C33C212
    (1)

    计算得到常压和100 GPa高压下B4C6O4的弹性力学常数见表2。经分析发现,独立的弹性常数Cij满足稳定性判据,证明B4C6O4在常压和100 GPa高压条件下满足弹性稳定性。

    表  2  常压和100 GPa下B4C6O4的弹性常数
    Table  2.  Elastic parameters of B4C6O4 at ambient pressure and 100 GPa  GPa
    PressureC11C12C13C22C23
    Ambient 693.28 15.78 77.13542.70167.35
    1001394.79146.17335.20974.72497.69
    PressureC33C44C55C66
    Ambient 516.07255.45197.31222.21
    1001000.58517.68363.80292.39
    下载: 导出CSV 
    | 显示表格

    虚频会导致晶格畸变,预示晶体的动力学不稳定。图2为常压和100 GPa高压下B4C6O4单胞整个布里渊区的声子色散谱和相关声子态密度(Density of states,DOS)。在其单胞的整个布里渊区没有负声子模式,表明常压和100 GPa高压下B4C6O4动力学稳定。由此可知,预测的B4C6O4同时满足弹性稳定性和动力学稳定性。

    图  2  常压和100 GPa下B4C6O4的声子散射谱和声子态密度
    Figure  2.  Phonon dispersion spectrum and phonon density of states of B4C6O4 at ambient pressure and 100 GPa

    随后基于实际环境开展对B4C6O4的分子动力学模拟:在恒定273 K和标准大气压的系综下进行总时间为3 ps、步长6 fs的模拟研究。模拟结果如图3所示。整个模拟过程中,结构没有明显变化。此外,为了阐明分子动力学模拟过程中结构的微调,这里一并选取了B―B键的键长以及C2与C2、O与O、B与B的原子间距(如图3中红色虚线所示),将B4C6O4的密度ρ作为观察对象,发现在整个分子动力学模拟过程中B4C6O4结构的局部化学键和原子间距(见表3)存在一定程度的拉伸或压缩,但总体变化微小,模拟初态和末态的体系密度几乎无变化。综上所述,即便在考虑温度和压力的真实环境下,B4C6O4依旧保持其结构完整性和稳定性。

    图  3  B4C6O4在恒温273 K和常压下的分子动力学模拟
    Figure  3.  Molecular dynamics simulation for B4C6O4 at 273 K and ambient pressure
    表  3  恒温常压条件下B4C6O4在分子动力学过程中的动态结构信息
    Table  3.  Structural information of B4C6O4 during molecular dynamics under constant temperature and ambient pressure
    Time/psρ/(g·cm−3)B―B bond length/ÅAtom distance/Å
    C2―C2O―OB―B
    03.1301.8422.7203.0382.757
    0.63.1751.7762.6773.0132.775
    1.23.1371.8692.6973.0262.722
    1.83.0661.8272.7413.0102.810
    2.43.1231.8462.7113.0632.806
    3.03.1281.8842.7383.0332.798
    下载: 导出CSV 
    | 显示表格

    基于密度泛函理论,B4C6O4在常压下的物理性质如电学、力学等得以研究。鉴于传统PBE泛函存在低估半导体带隙的不足[38],为了获取B4C6O4精确的电学性质,基于杂化泛函HSE06[39],分别在B4C6O4单胞和原胞模型上进行电学性质研究,结果如图4所示。虽然B4C6O4的单胞和原胞具有不同的倒易空间和布里渊区及高对称点路径,但是研究发现基于原胞和单胞计算所得的价带最高点和导带最低点均位于同一高对称点K上,最高价带(VBM)和最低导带(CBM)之间存在带隙,分别为2.248和2.247 eV。不同结构模型的研究均表明,B4C6O4具有带隙宽度约2.25 eV的直接带隙半导体属性。

    图  4  基于HSE06计算得到的常压下B4C6O4的电子能带结构和态密度:(a)原胞,(b)单胞(水平红线、暗青色曲线和金色曲线分别代表费米能级、VBM和CBM)
    Figure  4.  Calculated electronic band structures and density of states of B4C6O4 phases via HSE06 with primitive cell (a) and unit cell (b) at ambient pressure (The horizontal red line, dark cyan curve and gold curve represent the Fermi energy level, VBM and CBM, respectively.)

    纵观研究提出的B4C6O4以及Wang等[24]报道的B2CO2和B4CO4,三者均具有B、O含量相等的特点,即化学式可统一表达为BCxO(x = 3/2, 1/2, 1/4)。三者的结构中O呈现3配位成键,B和C均呈现4配位成键的诸多结构共性,为此通过理论研究探究三者的电学性质随组分的变化关系。鉴于杂化泛函HSE06对计算设备的硬件要求严苛、计算时间冗长等不足,而PBE泛函可以快速得到电学性质以及外界条件对电学性质的调控趋势,为此这里采用基于原胞模型和PBE交换关联泛函开展对BCxO(x = 3/2, 1/2, 1/4)的电学性质计算以及压力对电学性质的调控规律研究。

    首先研究常压下BCxO(x = 3/2, 1/2, 1/4)的电学性质,如图5所示。BCxO(x = 3/2, 1/2, 1/4)均属于半导体材料,其中B4C6O4为直接带隙,而B2CO2的价带最高点位于高对称点V,导带最低点位于LM之间的非高对称点,B4CO4的价带最高点位于高对称点G,导带最低点位于X,表明B2CO2和B4CO4均具有间接带隙。对比研究发现,随着C含量的降低,体系的带隙增大。

    图  5  基于PBE计算常压下BCxO原胞结构的电子能带结构和态密度(水平红线、暗青色曲线和金色曲线分别代表费米能级EF、VBM和CBM)
    Figure  5.  Calculated electronic band structures and density of states of three BCxO phases with primitive cell via PBE at ambient pressure (The horizontal red line, dark cyan curve and gold curve represent the Fermi energy level EF, VBM and CBM, respectively.)

    一般而言,压力作用下原子间距减小,体积减小,致密度增大。对于3种BCxO(x = 3/2, 1/2, 1/4)化合物,归一化体积随压力的变化如图6(a)所示。随着x的降低,归一化体积减小;B4C6O4、B2CO2和B4CO4三者的归一化体积均随压力的增大而持续降低。在所研究的压力范围内,B4C6O4、B2CO2、B4CO4的归一化体积分别减小了5.221、3.311、3.486 Å3,体系分别被压缩了21.955%、20.197%和22.764%。

    图  6  基于PBE泛函和原胞模型计算的BCxO(x = 3/2, 1/2, 1/4)的归一化体积(a)和带隙(b)与压力的关系
    Figure  6.  Relationship of the normalized volume (a) and the calculated electronic band gaps (b) of BCxO (x = 3/2, 1/2, 1/4) with pressure based on PBE funtional and primitive cell

    凝聚态材料的电学性质会受压力的影响[40]。通常而言,高压下原子间距变小,电子重叠度增大,电子不局域于单个原子或键,即形成离域电子,导致带隙减小乃至金属化,典型代表就是金属氢[41]。然而高压下,电子出现局域化的同时也可能伴随成键态和反键态的形成,导致价带能量更低而导带能量更高,即出现带隙变大现象,典型代表就是钠绝缘化[42]。很多时候两种作用相互竞争,共同影响材料的电学性质。

    为此,研究了BCxO(x = 3/2, 1/2, 1/4)带隙随压力的变化。基于原胞和PEB算法的研究结果如图6(b)所示。可见,B2CO2的带隙有较明显的受压变化趋势,而B4CO4和B4C6O4的带隙受压力的影响较小。在所研究的压力范围内,B2CO2的带隙随压力的升高呈现持续下降现象,100 GPa内带隙减小425 meV,即在高压的作用下B2CO2的晶格参数变小,布里渊区变大,能带展宽,进而导致带隙减小。B4C6O4的带隙也呈现随压力的升高而降低的变化规律,在低于100 GPa的压力范围内带隙降低了168 meV,降幅明显小于B2CO2。对于B4CO4,带隙在常压到40 GPa区间随压力的升高而持续增大,增幅达62 meV,随后保持带隙最大值5.495 eV直至50 GPa,50 GPa之后带隙随压力的升高而缓慢下降,100 GPa时带隙为5.470 eV。纵观整个升压过程,B4CO4的带隙变化微小,即两种作用相互竞争,近乎相持。

    通过模拟特定方向的结构应变获得凝聚态材料的内在应力是理解结构变形和强度的常用方法。材料的理想强度(一个完美晶体在弹性力学上由稳定变为不稳定时的临界应力)代表了材料强度的上限。通过模拟拉伸过程的应变-应力关系和键断裂过程,可以从微观原子层面深入理解结构变形乃至失效的理论机制。为此,基于密度泛函理论模拟研究了BCxO(x = 3/2, 1/2, 1/4)在常见晶向上的拉伸过程,所得应力-应变关系如图7所示。

    图  7  模拟得到的3种BCxO沿特定方向的拉伸应力-应变关系
    Figure  7.  Simulated tensile stress-strain relationship along the specific direction for three BCxO phases

    B4C6O4为正交晶系结构,3个常见方向的应力-应变关系如图7(a)所示。[100]方向上应力先随应变的增大而增大,在应变为0.1125时达到峰值42.1 GPa,随后应力略微减小,呈现快速失效特征。[010]和[001]方向上,在低于0.05的小应变区域内应力相近,应变继续增大时应力差异显现,其中[001]方向的应力在应变为0.10时达到最大值40.4 GPa,此后继续增大拉伸应力会导致应变超过极限变形量而造成结构失效,而[010]方向上的应力随应变的增大持续增大到49.5 GPa(应变0.1625),随后应力随应变的增大而略微下降,结构快速失效。

    B2CO2为单斜晶系结构,常见方向的拉伸模拟研究结果如图7(b)所示。[100]和[001]方向的应力-应变曲线具有相似性,二者的应力先随着应变的增大而增大,均在0.056处达到最大值,分别为21.0和23.5 GPa,随后应力随应变的增大而减小,最大应变均为0.069,应变继续增大会导致结构失效。在[010]方向,低应变区域内应力随着应变的增大而快速增大,高应变区域内应力随应变的增大而缓慢增大,最大应力为66.7 GPa,此时应变也达到极限值0.225,随后若应力继续增大,则结构失效。

    B4CO4为四方晶系,同样选取了[100]、[110]和[001] 3个方向,其应力-应变曲线如图7(c)所示。在[100]和[001]方向上应力随应变的增大而增大,其中:[100]方向的最大应力为42.3 GPa,此时应变达到极限0.138;[001]方向的最大应力31.0 GPa出现在极限应变0.100处,继续增大应力会导致应变超过极限应变,进而使结构失效塌陷;在[110]方向,应力先随应变的增大而显著增大,当应变为0.100时应力达到最大值55.9 GPa,随后应变继续增大,结构中部分B―O键断裂(图7(c)S1S2),应力突降至11.2 GPa,此时结构未完全失效,当应变继续增大时应力持续增大,并在0.2875应变处应力达到38.0 GPa,当应变继续增大,结构完全失效。二次变形的存在导致结构在[110]方向的应变高于[100]和[001]方向。

    基于模拟的拉伸过程,在结构未失效下开展了3种BCxO(x = 3/2, 1/2, 1/4)的带隙-应变关系研究,结果见图8。如图8(a)所示,在[001]方向的拉伸过程中,B4C6O4的带隙随应变的增大单调增大;在[100]方向的拉伸过程中,带隙随应变的增大持续下降;而[010]方向的拉伸则引起带隙先缓慢上升,在应变0.063处达到最大值1.068 eV,随后带隙先均匀降低,在临近失效阶段快速降低。如图8(b)所示,B2CO2在3个方向的拉伸应力作用下,带隙随结构的改变先快速下降,在达到一定的应变后,[001]方向的拉伸应变引起带隙持续下降,而[100]和[010]方向的应变则均先使体系带隙缓慢上升,达到峰值后,随应变的增大而降低,其中[100]方向的带隙峰值4.155 eV出现在应变0.069处,[010]方向的带隙峰值4.173 eV出现在应变0.050处。B4CO4的带隙随[100]、[110]和[001] 3个方向的拉伸应变的变化如图8(c)所示。[001]方向的应变使带隙先升后降,峰值6.218 eV出现在应变0.075处;在[100]方向上带隙随应变先微弱上升后持续降低;而[110]方向上带隙的变化与应力的变化一样,出现两个阶段,即应变在0.100以内带隙快速下降,应变由0.100增大到0.113的过程中带隙出现突降,随后带隙缓慢下降。

    图  8  模拟得到的3种BCxO在沿特定方向拉伸过程中的电阻-应变关系
    Figure  8.  Simulated gap-strain relationship along the specific direction for three BCxO phases during the tensile process

    基于密度泛函理论,系统研究了3种BCxO(x = 3/2, 1/2, 1/4)的力学性质。首先,基于PBE泛函计算了BCxO(x = 3/2, 1/2, 1/4)的弹性常数Cij,通过Cij计算Voigt和Reuss形式的体积模量B和剪切模量G[43],并经改进型硬度计算公式计算维氏硬度HV[44]。改进型硬度计算公式为

    HV=0.92k1.137G0.708,k=G/B
    (2)

    3种BCxO(x = 3/2, 1/2, 1/4)的体积模量B和剪切模量G与压力的关系如图9(a)图9(b)所示,其中Voigt和Reuss形式数值分别代表模量的上限和下限。对于体积模量而言,B4CO4的体积模量的Voigt和Reuss形式数值非常接近,相差不到0.8 GPa,B2CO2和B4C6O4的体积模量的Voigt和Reuss形式数值差值均小于8 GPa。结合图9(a)可知,3种BCxO(x = 3/2, 1/2, 1/4)均具有较大的体积模量。此外,3种BCxO(x = 3/2, 1/2, 1/4)的体积模量均呈现与压力正相关的线性关系。如图9(b)所示,对于剪切模量而言,3种BCxO(x = 3/2, 1/2, 1/4)的Voigt和Reuss形式数值均存在较大差值,且该差值随压力的升高而增大。根据Voigt和Reuss形式的体积模量和剪切模量,计算得到两种形式的维氏硬度HV及其与压力的关系,如图9(c)所示。常压下3种BCxO(x = 3/2, 1/2, 1/4)均具有较高的硬度,尤其是新预测的B4C6O4具有40.3~37.3 GPa的常压硬度,即B4C6O4是潜在的常压准“超硬”材料;B2CO2和B4CO4也分别具有39.2~36.7 GPa、38.9~33.7 GPa的常压硬度,表明二者具有常压高硬性质。除B4C6O4的硬度在常压到10 GPa区间呈现小幅上升外,3种BCxO(x = 3/2, 1/2, 1/4)的硬度随压力的升高呈现下降趋势,在100 GPa下依旧保持高硬特性。

    图  9  3种BCxO(x = 3/2, 1/2, 1/4)的体积模量(a)、剪切模量(b)、硬度(c)与压力的关系
    Figure  9.  Relationship of bulk modulus (a), shear modulus (b), and hardness (c) with pressure for BCxO (x = 3/2, 1/2, 1/4)

    晶体的声子振动频率受其结构中化学键强度的影响,即最大声子振动频率在一定程度上反映了化学键的强弱。基于不同压力下BCxO(x = 3/2, 1/2, 1/4)的结构,计算得到其声子色散曲线和声子态密度,分析发现在常压到100 GPa的压力范围内BCxO(x = 3/2, 1/2, 1/4)均满足动力学稳定性。基于声子色散曲线,详细探讨了3种BCxO(x = 3/2, 1/2, 1/4)的最大振动频率随压力变化的关系。如图10(a)所示,常压下BCxO的最高声子振动频率均高于30 THz,且呈B4CO4 > B2CO2 > B4C6O4的关系。此外,随着压力的增大,原子间距减小,成键原子的相互作用力增强,即化学键增强,导致体系的最高声子振动频率随压力同步递增。

    图  10  BCxO(x = 3/2, 1/2, 1/4)的最大声子振动频率(a)和零点振动能(b)与压力的关系
    Figure  10.  Relationship of maximum vibrational frenquency (a), and zero-point vibration energy (b) with pressure for BCxO (x = 3/2, 1/2, 1/4)

    零点振动能EZP是量子在绝对温度0 K下的振动能量。零点振动的振幅随温度的升高而增大,原子越轻,零点振动越明显。根据声子散射谱和声子态密度,可以计算出零点振动能EZP

    EZP=0.5F(ω)ωdω
    (3)

    式中:ћ为普朗克常数,F(ω)为声子态密度。基于相关的声子散射谱,进一步研究了归一化零点振动能与压力的关系。

    图10(b)所示,常压下B4C6O4的归一化零点振动能最大,其次是B2CO2,B4CO4最小,与三者的归一化相对质量大小密切相关。此外,随着压力的增大,3种BCxO的零点振动能均增大,且增幅也随x的增大而递减,EZP的增幅分别为134.8、87.8和80.1 meV。以上说明C含量对BCxO的最大振动频率和零点振动能具有显著的影响。

    通过粒子群优化算法生成候选结构,研究了具有正交晶系结构的新型B-C-O化合物B4C6O4。结合第一性原理的声子散射谱和独立弹性常数,证明了B4C6O4结构的常压和高压稳定性。B4C6O4具有带隙宽度约2.25 eV的直接带隙半导体属性。研究同属BCxO系列且结构具有相似性的B4C6O4、B2CO2和B4CO4发现,C含量的降低导致带隙增大,三者的分子式体积随C含量的降低而降低,且100 GPa的高压对三者体积均形成高达20%的压缩。高压作用导致B2CO2和B4C6O4的带隙持续降低,而B4CO4的带隙先升后降,整体变化微弱。应力-应变模拟结果表明,B4C6O4的最大应力49.5 GPa出现在[010]方向的0.1625应变处,[100]和[001]方向能承受的最大应力及应变均低于[010]方向;同样,B2CO2在[010]方向能承受的最大拉伸应力66.7 GPa和应变0.225也远大于[100]和[001]方向的应力和应变。B4CO4在最大应力和最大应变的[110]方向上出现二次拉伸特质,应变为0.100时应力达到最大值55.9 GPa,增大应变至0.113时结构中出现部分断键但结构未塌陷,继续增大应变到超过0.288时结构的完整性被破坏。应力引起的应变会影响3种BCxO(x = 3/2, 1/2, 1/4)的带隙。3种BCxO化合物的力学性能研究表明,它们都具有高弹性模量和高硬度,尤其是新预测的B4C6O4具有40.3~37.3 GPa的常压硬度,是潜在的准“超硬”材料。三者在100 GPa下依旧保持高硬特性。声子散射谱研究表明,常压下BCxO的最高声子振动频率均高于30 THz,且B4CO4的最高声子振动频率最大,B2CO2次之,B4C6O4最小,高压作用会引起体系键能的持续增强。

  • 图  气泡帷幕发生器

    Figure  1.  Bubble curtain generator

    图  现场试验中气泡帷幕发生器的作用效果

    Figure  2.  Effect of bubble curtain generator in field test

    图  现场布置示意图

    Figure  3.  Schematic diagram of test arrangement

    图  空白对照组的冲击波压力时程曲线

    Figure  4.  Time histories of shock wave pressure in blank control group

    图  供风量为30、60、90 L/min时的冲击波压力时程曲线

    Figure  5.  Time histories of shock wave pressure with air supply volume of 30,60 and 90 L/min

    图  供风量为30 L/min、不同气泡帷幕层数时的冲击波压力时程曲线

    Figure  6.  Time histories of shock wave pressure under different layers of bubble curtains with the air supply volume of 30 L/min

    图  相邻层数之间的衰减效率Bi

    Figure  7.  Attenuation efficiency Bi between adjacent layers

    图  相邻梯度供风量之间的衰减效率Cj

    Figure  8.  Attenuation efficiency Cj between adjacent air supply volumes

    表  1  试验设计方案

    Table  1.   Design scheme of field test

    Test No. v/(L·min−1) n Test No. v/(L·min−1) n Test No. v/(L·min−1) n
    1 0 0 8 90 1 15 30 3
    2 0 0 9 30 2 16 30 3
    3 30 1 10 30 2 17 60 3
    4 30 1 11 60 2 18 60 3
    5 60 1 12 60 2 19 90 3
    6 60 1 13 90 2 20 90 3
    7 90 1 14 90 2
    下载: 导出CSV

    表  2  不同层数气泡帷幕条件下水下冲击波峰值压力

    Table  2.   Peak pressure of underwater explosion shock wave under different layers of bubble curtains

    v/(L·min−1) n=1 n=2 n=3
    pm/MPa ¯pm/MPa A/% pm/MPa ¯pm/MPa A/% pm/MPa ¯pm/MPa A/%
    30 2.069
    1.991
    2.030 49.67 0.931
    0.944
    0.938 76.74 0.831
    0.870
    0.851 78.90
    60 1.547
    1.501
    1.524 62.21 0.842
    0.801
    0.822 79.62 0.633
    0.681
    0.657 83.71
    90 0.760
    0.796
    0.778 80.71 0.655
    0.684
    0.670 83.39 0.390
    0.392
    0.391 90.30
    下载: 导出CSV
  • [1] 王兴雁, 詹发民, 周方毅, 等. 气泡帷幕削减水击波压力作用因素分析 [J]. 爆破, 2012, 29(4): 23–27. doi: 10.3963/j.issn.1001-487X.2012.04.006

    WANG X Y, ZHAN F M, ZHOU F Y, et al. Effect of bubble curtains on underwater shockwave reducing [J]. Blasting, 2012, 29(4): 23–27. doi: 10.3963/j.issn.1001-487X.2012.04.006
    [2] 寇晓枫. 空气夹层结构的水下爆炸防护效果研究 [D]. 武汉: 武汉大学, 2017.

    KOU X F. Protection effects of sandwich structures with air interlayer subjected to underwater explosion [D]. Wuhan: Wuhan University, 2017.
    [3] 张轶凡, 刘亮涛, 王金相, 等. 水下爆炸冲击波和气泡载荷对典型圆柱壳结构的毁伤特性 [J]. 兵工学报, 2023, 44(2): 345–359. doi: 10.12382/bgxb.2021.0598

    ZHANG Y F, LIU L T, WANG J X, et al. Damage characteristics of underwater explosion shock wave and bubble load on typical cylindrical shell structure [J]. Acta Armamentarii, 2023, 44(2): 345–359. doi: 10.12382/bgxb.2021.0598
    [4] 王高辉, 高政, 卢文波, 等. 考虑初始应力的混凝土重力坝水下爆炸毁伤特性研究 [J]. 振动与冲击, 2022, 41(11): 133–140. doi: 10.13465/j.cnki.jvs.2022.11.017

    WANG G H, GAO Z, LU W B, et al. Damage characteristics of underwater explosion of concrete gravity dam considering initial stress [J]. Journal of Vibration and Shock, 2022, 41(11): 133–140. doi: 10.13465/j.cnki.jvs.2022.11.017
    [5] 张志波, 李春军, 李红勇, 等. 气泡帷幕在水下爆破减震工程中的应用 [J]. 爆破, 2003, 20(2): 75–76, 89. doi: 10.3963/j.issn.1001-487X.2003.02.028

    ZHANG Z B, LI C J, LI H Y, et al. Application of air bubble purdah in the damping measure in the underwater blasting [J]. Blasting, 2003, 20(2): 75–76, 89. doi: 10.3963/j.issn.1001-487X.2003.02.028
    [6] 贾虎, 郑伟花, 罗强, 等. 爆炸气泡帷幕对水中冲击波能量的衰减特性 [J]. 含能材料, 2015, 23(10): 1015–1019. doi: 10.11943/j.issn.1006-9941.2015.10.018

    JIA H, ZHENG W H, LUO Q, et al. Attenuation characteristics of underwater explosion bubble curtain on the shock [J]. Chinese Journal of Energetic Materials, 2015, 23(10): 1015–1019. doi: 10.11943/j.issn.1006-9941.2015.10.018
    [7] 司剑峰. 深水钻孔爆破的冲击波衰减规律及防护研究 [D]. 武汉: 武汉科技大学, 2021.

    SI J F. Research on the attenuation law of shock wave and protection in deep-water drilling and blasting [D]. Wuhan: Wuhan University of Science and Technology, 2021.
    [8] 陆少锋, 梁进, 覃才勇, 等. 供风量对水下爆炸冲击波气泡帷幕消波效应的影响 [J]. 工程爆破, 2022, 28(5): 143–148. doi: 10.19931/j.EB.20210013

    LU S F, LIANG J, QIN C Y, et al. Influence of air supply rate on wave attenuation effect of bubble curtain for underwater explosion shock wave [J]. Engineering Blasting, 2022, 28(5): 143–148. doi: 10.19931/j.EB.20210013
    [9] 谢金怀, 何树斌, 屈科, 等. 气泵法生成气泡帷幕的特性研究 [J]. 海洋技术学报, 2019, 38(1): 12–17. doi: 10.3969/j.issn.1003-2029.2019.01.003

    XIE J H, HE S B, QU K, et al. Study on the characteristics of bubble curtain generated by the air pump method [J]. Journal of Ocean Technology, 2019, 38(1): 12–17. doi: 10.3969/j.issn.1003-2029.2019.01.003
    [10] 范怀斌, 陆少锋, 莫崇勋, 等. 多层差异性气泡帷幕对水下爆破冲击波的衰减效应的试验研究 [J]. 爆破器材, 2023, 52(2): 48–55. doi: 10.3969/j.issn.1001-8352.2023.02.008

    FAN H B, LU S F, MO C X, et al. Experimental study on attenuation effect of multi-layer differential bubble curtain on underwater blasting shock wave [J]. Explosive Materials, 2023, 52(2): 48–55. doi: 10.3969/j.issn.1001-8352.2023.02.008
    [11] 刘欣, 顾文彬, 陈学平. 气泡帷幕对水中冲击波衰减特性的数值模拟研究 [J]. 爆破, 2015, 32(3): 79–84. doi: 10.3963/j.issn.1001-487X.2015.03.014

    LIU X, GU W B, CHEN X P. Numerical simulation study of attenuation characteristics of water shock wave under bubble curtain [J]. Blasting, 2015, 32(3): 79–84. doi: 10.3963/j.issn.1001-487X.2015.03.014
    [12] 张成兴, 王永学, 王国玉, 等. 静水中气泡帷幕产生水平流的数值模拟研究 [J]. 水动力学研究与进展A辑, 2010, 25(1): 59–66. doi: 10.3969/j.issn.1000-4874.2010-01.009

    ZHANG C X, WANG Y X, WANG G Y, et al. Numerical simulation study on the horizontal current generated by air bubbles curtain in still water [J]. Chinese Journal of Hydrodynamics, 2010, 25(1): 59–66. doi: 10.3969/j.issn.1000-4874.2010-01.009
    [13] 鲁天龙. 气泡帷幕周围流场运动特性数值模拟研究 [D]. 长沙: 长沙理工大学, 2020.

    LU T L. Numerical simulation of flow field motion characteristics around bubble curtain [D]. Changsha: Changsha University of Science & Technology, 2020.
    [14] 刘天云, 龚书堂, 胡伟才, 等. 水下钻孔爆破水击波的传播规律及气泡帷幕对水击波的削减作用 [J]. 爆破器材, 2020, 49(2): 16–22. doi: 10.3969/j.issn.1001-8352.2020.02.003

    LIU T Y, GONG S T, HU W C, et al. Propagation law of water hammer wave in underwater drilling blasting and reduction of bubble curtain on water hammer wave [J]. Explosive Materials, 2020, 49(2): 16–22. doi: 10.3969/j.issn.1001-8352.2020.02.003
    [15] 谢达建, 吴立, 洪江, 等. 气泡帷幕对水下爆破冲击波的削弱作用研究 [J]. 人民长江, 2018, 49(8): 72–77. doi: 10.16232/j.cnki.1001-4179.2018.08.014

    XIE D J, WU L, HONG J, et al. Study on weakening effect of bubble curtain on water shock wave in underwater blasting [J]. Yangtze River, 2018, 49(8): 72–77. doi: 10.16232/j.cnki.1001-4179.2018.08.014
    [16] 胡亚峰, 金建峰, 顾文彬, 等. 爆炸实验水池防护性能及动力学响应分析 [J]. 爆炸与冲击, 2017, 37(6): 1001–1009. doi: 10.11883/1001-1455(2017)06-1001-09

    HU Y F, JIN J F, GU W B, et al. Protective performance and dynamic response analysis of explosion testing pool [J]. Explosion and Shock Waves, 2017, 37(6): 1001–1009. doi: 10.11883/1001-1455(2017)06-1001-09
    [17] 胡伟才, 吴立, 舒利, 等. 不同设置方式下气泡帷幕对水中冲击波衰减特性的影响 [J]. 科学技术与工程, 2018, 18(17): 33–38. doi: 10.3969/j.issn.1671-1815.2018.17.006

    HU W C, WU L, SHU L, et al. Influence of water shock wave on attenuation characteristics under bubble curtain with different settings [J]. Science Technology and Engineering, 2018, 18(17): 33–38. doi: 10.3969/j.issn.1671-1815.2018.17.006
    [18] 彭亚雄, 吴立, 李春军, 等. 水下钻孔爆破水击波特性及气泡帷幕削压效果研究 [J]. 爆破, 2019, 36(1): 38–43. doi: 10.3963/j.issn.1001-487X.2019.01.006

    PENG Y X, WU L, LI C J, et al. Characteristics of water shock wave from underwater hole blasting and weakening pressure effect of bubble curtain in water [J]. Blasting, 2019, 36(1): 38–43. doi: 10.3963/j.issn.1001-487X.2019.01.006
    [19] 俞统昌, 王晓峰, 王建灵. 炸药的水下爆炸冲击波性能 [J]. 含能材料, 2003, 11(4): 182–186. doi: 10.3969/j.issn.1006-9941.2003.04.002

    YU T C, WANG X F, WANG J L. Underwater shockwave performance of explosives [J]. Energetic Materials, 2003, 11(4): 182–186. doi: 10.3969/j.issn.1006-9941.2003.04.002
    [20] 冯沐桦. 水下爆炸与水池抗冲击载荷研究 [D]. 南京: 南京理工大学, 2023.

    FENG M H. Study on underwater explosion and methods of anti-shock load for explosive pool [D]. Nanjing: Nanjing University of Science & Technology, 2023.
  • 加载中
图(8) / 表(2)
计量
  • 文章访问数:  150
  • HTML全文浏览量:  64
  • PDF下载量:  26
出版历程
  • 收稿日期:  2023-08-08
  • 修回日期:  2023-09-14
  • 网络出版日期:  2024-01-29
  • 刊出日期:  2024-02-05

目录

/

返回文章
返回