镁铝合金的冲击熔化行为实验研究

谭叶 肖元陆 薛桃 李俊 金柯

引用本文:
Citation:

镁铝合金的冲击熔化行为实验研究

    作者简介: 谭 叶(1986-),男,硕士,副研究员,主要从事材料高压物性研究. E-mail: yetan@caep.cn;
    通讯作者: 金柯, jinke102@caep.cn
  • 中图分类号: O521.2

Melting of MB2 Alloy under Shock Compression

    Corresponding author: JIN Ke, jinke102@caep.cn
  • CLC number: O521.2

  • 摘要: 采用反向碰撞实验技术,结合具有高时空分辨率的全光纤激光干涉测速技术,对镁铝合金开展了极端动态压缩条件下的动力学行为实验研究,获得了镁铝合金在30~73 GPa压力范围内的Hugoniot和声速实验数据。深入的数据分析表明,所获得的Hugoniot数据与早期的Hugoniot数据一致,但是纵波声速却呈现出明显的向体波声速转变的趋势,对应冲击加载下镁铝合金的固-液熔化相变,相变压力区间为40~57 GPa。
  • 图 1  反向碰撞实验示意图

    Figure 1.  Schematic of backward-impact experimental configuration

    图 2  反向碰撞实验p-u

    Figure 2.  p-u relation for backward-impact experiment

    图 3  MB2/LiF界面粒子速度剖面

    Figure 3.  Particle velocity profile of MB2/LiF interface

    图 4  冲击波速度与粒子速度的关系

    Figure 4.  Shock velocity vs. particle velocity

    图 5  声速与冲击压力的关系

    Figure 5.  Sound velocity vs. shock pressure

    图 6  反碰撞法测量声速实验中不确定度

    Figure 6.  Uncertainties for backward-impact experiment based on law of propagation

    表 1  MB2样品冲击实验参数及结果

    Table 1.  Shock experiment parameters and results of MB2

    Exp.No. hs/mm W/(km·s–1) uw/(km·s–1) us/(km·s–1) Ds/(km·s–1) p/GPa Cl/(km·s–1)
    1 1.980±0.004 3.949±0.020 1.589±0.016 2.360±0.026 7.303±0.166 30.6±0.4 7.983±0.297
    2 1.997±0.004 4.358±0.020 1.763±0.018 2.595±0.027 7.607±0.175 35.0±0.5 9.101±0.394
    3 1.978±0.004 5.379±0.027 2.195±0.022 3.184±0.030 8.317±0.191 47.0±0.7 9.167±0.402
    4 1.981±0.004 5.928±0.030 2.435±0.024 3.493±0.039 8.746±0.214 54.2±0.9 8.912±0.380
    5 1.981±0.004 6.100±0.030 2.514±0.025 3.586±0.039 8.907±0.213 56.7±0.9 8.601±0.348
    6 1.987±0.004 7.220±0.036 3.003±0.030 4.217±0.050 9.747±0.245 73.0±1.2 9.723±0.463
    下载: 导出CSV
  • [1] GAO C Y, ZHANG L C, GUO W G, et al. Dynamic plasticity of AZ31 magnesium alloy: experimental investigation and constitutive modeling [J]. Materials Science and Engineering A, 2014, 613: 379–389. doi: 10.1016/j.msea.2014.06.112
    [2] 胡昌明, 李英雷, 胡时胜, 等. 高温-高应变率下MB2合金的动态力学性能及变形机理 [J]. 兵器材料科学与工程, 2009, 32(5): 8–11 doi: 10.3969/j.issn.1004-244X.2009.05.003
    HU C M, LI Y L, HU S S, et al. Dynamic mechanical properties and deformation mechanical of MB2 alloy under high temperature and high strain rates [J]. Ordnance Material Science and Engineering, 2009, 32(5): 8–11 doi: 10.3969/j.issn.1004-244X.2009.05.003
    [3] SCHMIDT R M, DAVIES F W, LEMPRIERE B M, et al. Temperature dependent spall threshold of four metal alloys [J]. Journal of Physics and Chemistry of Solids, 1978, 39(4): 375–385. doi: 10.1016/0022-3697(78)90079-3
    [4] MILLETT J C F, STIRK S M, BOURNE N K, et al. On the behaviour of the magnesium alloy, AZ61 to one-dimensional shock loading [J]. Acta Materialia, 2010, 58(17): 5675–5682. doi: 10.1016/j.actamat.2010.06.042
    [5] MARSH S P. LASL shock Hugoniot data [M]. Berkeley: University of California Press, 1980: 208.
    [6] YU Y Y, TAN H, HU J B, et al. Shear modulus of shock-compressed LY12 aluminium up to melting point [J]. Chinese Physics B, 2008, 17(1): 264–269. doi: 10.1088/1674-1056/17/1/046
    [7] 胡建波, 谭华, 俞宇颖, 等. 铝的动态屈服强度测量及再加载弹性前驱波的形成机理分析 [J]. 物理学报, 2008, 57(1): 405–410 doi: 10.3321/j.issn:1000-3290.2008.01.063
    HU J B, TAN H, YU Y Y, et al. Measurements of dynamic yield strength of aluminum alloy and mechanism analysis of elastic precursor during reloading [J]. Acta Physica Sinica, 2008, 57(1): 405–410 doi: 10.3321/j.issn:1000-3290.2008.01.063
    [8] 谭叶, 俞宇颖, 戴诚达, 等. 反向碰撞法测量Bi的低压Hugoniot数据 [J]. 物理学报, 2011, 60(10): 106401 doi: 10.7498/aps.60.106401
    TAN Y, YU Y Y, DAI C D, et al. Measurement of low-pressure Hugoniot data for bismuth with reverse-impact geometry [J]. Acta Physica Sinica, 2011, 60(10): 106401 doi: 10.7498/aps.60.106401
    [9] 俞宇颖, 谭叶, 戴诚达, 等. 钒的高压声速测量 [J]. 物理学报, 2014, 63(2): 026202
    YU Y Y, TAN Y, DAI C D, et al. Sound velocities of vanadium under shock compression [J]. Acta Physica Sinica, 2014, 63(2): 026202
    [10] WENG J D, TAN H, WANG X, et al. Optical-fiber interferometer for velocity measurements with picosecond resolution [J]. Applied Physics Letters, 2006, 89(11): 111101. doi: 10.1063/1.2335948
    [11] 谭华. 实验冲击波物理导 [M]. 北京: 国防工业出版社, 2006: 5.
    TAN H. Introduction to experimental shock-wave physics [M]. Beijing: National Defense Industry Press, 2006: 5.
    [12] RIGG P A, KNUDSON M D, SCHARFF R J, et al. Determining the refractive index of shocked [J]. Journal of Applied Physics, 2014, 116(3): 033515. doi: 10.1063/1.4890714
    [13] URTIEW P A, GROVER R. The melting temperature of magnesium under shock loading [J]. Journal of Applied Physics, 1977, 48(3): 1122–1126. doi: 10.1063/1.323789
    [14] 徐锡申, 张万箱. 实用物态方程理论导引 [M]. 北京: 科学出版社, 1986: 527.
    XU X S, ZHANG W X. Introduction to equation of state theory [M]. Beijing: National Defense Industry Press, 1986: 527.
    [15] WU Q, JING F Q. Thermodynamic equation of state and application to Hugoniot predictions for porous materials [J]. Journal of Applied Physics, 1996, 80(8): 4343–4349. doi: 10.1063/1.363391
    [16] 国家质量技术监督局计量司. 测量不确定度评定与表示指南 [M]. 北京: 中国计量出版社, 2000.
  • [1] 俞宇颖谭华胡建波戴诚达马云陈大年 . 利用VISAR测量LY12铝在冲击压缩下的声速. 高压物理学报, 2006, 20(2): 145-152 . doi: 10.11858/gywlxb.2006.02.006
    [2] 张修路袁长迎孙悦韦建军董石吴国栋 . 壬基酚聚氧乙烯醚的冲击压缩实验研究. 高压物理学报, 2005, 19(2): 164-168 . doi: 10.11858/gywlxb.2005.02.011
    [3] 金柯吴强李加波周显明叶素华李俊 . 冲击加载下NaCl单晶高压声速与温度的同步测量. 高压物理学报, 2017, 31(6): 707-717. doi: 10.11858/gywlxb.2017.06.005
    [4] 陈永涛李庆忠唐小军刘文韬 . FeMnNi合金冲击相变与卸载逆相变历程研究. 高压物理学报, 2011, 25(2): 153-158 . doi: 10.11858/gywlxb.2011.02.011
    [5] 李庆忠陈永涛胡海波徐永波 . 冲击加载下FeMnNi合金相变和层裂特性实验研究. 高压物理学报, 2010, 24(2): 107-112 . doi: 10.11858/gywlxb.2010.02.005
    [6] 王金贵李欣竹王为傅秋卫王翔 . 三种钨合金的冲击压缩性. 高压物理学报, 1998, 12(4): 258-263 . doi: 10.11858/gywlxb.1998.04.004
    [7] 杨金科龚自正邓力维张莉费英伟 . 顽火辉石(Mg0.92,Fe0.08)SiO3的冲击相变和高压状态方程及其地球物理意义. 高压物理学报, 2007, 21(1): 45-54 . doi: 10.11858/gywlxb.2007.01.008
    [8] 陈永涛洪仁楷汤铁钢陈浩玉 . 熔化状态下锡样品微喷射现象的诊断. 高压物理学报, 2016, 30(4): 323-327. doi: 10.11858/gywlxb.2016.04.009
    [9] 陈栋泉李茂生蔡灵仓 . 高压熔化的自由体积理论. 高压物理学报, 2001, 15(2): 81-86 . doi: 10.11858/gywlxb.2001.02.001
    [10] 冯磊黄海军冷春蔚杨刚 . 高压下Fe92.5O2.2S5.3的熔化温度. 高压物理学报, 2017, 31(6): 698-706. doi: 10.11858/gywlxb.2017.06.004
    [11] 郑伟涛张瑞林 . 铝、铜、铅压力熔化曲线的研究. 高压物理学报, 1993, 7(1): 42-46 . doi: 10.11858/gywlxb.1993.01.006
    [12] 熊大和B. C. MingM. H. Manghnani . 钙钛矿(GaTiO3)高压相变及等温压缩. 高压物理学报, 1988, 2(1): 1-9 . doi: 10.11858/gywlxb.1988.01.001
    [13] 李庆忠陈永涛汤铁钢 . FeMnNi合金高压加卸载历程测量和相变层裂分析. 高压物理学报, 2009, 23(2): 117-122 . doi: 10.11858/gywlxb.2009.02.007
    [14] 李厚深杨世源向芸蔡灵仓王军霞 . 平面冲击波作用下TiO2晶粒的相变纳米化. 高压物理学报, 2010, 24(5): 395-400 . doi: 10.11858/gywlxb.2010.05.013
    [15] 徐康贺红亮谭华刘建军 . 冲击波温度和压力对二氧化钛相变的影响. 高压物理学报, 1997, 11(1): 27-31 . doi: 10.11858/gywlxb.1997.01.005
    [16] 马云张毅李加波蔡灵仓李剑峰柳雷敬秋民王志刚王翔翁继东 . DAC加载下材料高压声速的激光超声测量技术. 高压物理学报, 2011, 25(5): 416-420 . doi: 10.11858/gywlxb.2011.05.006
    [17] 李敏崔启良李芳菲周强武晓鑫邹广田 . 高温高压条件下液态氧的声速及弹性系数研究. 高压物理学报, 2008, 22(3): 286-290 . doi: 10.11858/gywlxb.2008.03.011
    [18] 马云张毅李加波蔡灵仓李剑峰王翔翁继东李军敬秋民李俊 . 利用激光超声技术测量介观尺度材料的声速. 高压物理学报, 2011, 25(2): 128-132 . doi: 10.11858/gywlxb.2011.02.006
    [19] 陈军曾代朋谭多望王荣波 . 超压爆轰产物声速以及热力学CJ点的实验研究. 高压物理学报, 2010, 24(3): 201-205 . doi: 10.11858/gywlxb.2010.03.007
    [20] 龚自正华劲松经福谦谢鸿森郭捷 . 高压下声速温度系数的一种新算法. 高压物理学报, 2000, 14(4): 241-246 . doi: 10.11858/gywlxb.2000.04.001
  • 加载中
图(7)表(1)
计量
  • 文章访问数:  1597
  • 阅读全文浏览量:  1450
  • PDF下载量:  14
出版历程
  • 收稿日期:  2019-02-26
  • 录用日期:  2019-03-14
  • 网络出版日期:  2019-06-25
  • 刊出日期:  2019-04-01

镁铝合金的冲击熔化行为实验研究

    作者简介:谭 叶(1986-),男,硕士,副研究员,主要从事材料高压物性研究. E-mail: yetan@caep.cn
    通讯作者: 金柯, jinke102@caep.cn
  • 中国工程物理研究院流体物理研究所冲击波物理与爆轰物理重点实验室,四川 绵阳 621999

摘要: 采用反向碰撞实验技术,结合具有高时空分辨率的全光纤激光干涉测速技术,对镁铝合金开展了极端动态压缩条件下的动力学行为实验研究,获得了镁铝合金在30~73 GPa压力范围内的Hugoniot和声速实验数据。深入的数据分析表明,所获得的Hugoniot数据与早期的Hugoniot数据一致,但是纵波声速却呈现出明显的向体波声速转变的趋势,对应冲击加载下镁铝合金的固-液熔化相变,相变压力区间为40~57 GPa。

English Abstract

  • 镁铝合金(MB2)作为一类特殊的合金材料,具有低密度、高强度、易机械加工、耐腐蚀等特点,广泛应用于车辆工程、航空航天等领域,其动态加载下的力学和物理特性对相关结构设计等具有重要意义。国内外对MB2合金的早期研究主要集中于力学特性方面,通过开展低压动态响应特性实验研究,获得了材料的弹塑性响应[1]、动态损伤[2]及层裂行为[3]的初步认识,近年来Millett等[4]通过一维冲击加载研究了MB2合金在早期变形和位错条件下随加载应力、加载脉宽而变化的弹塑性和剪切强度行为。另外,对动态加载下MB2合金的物态方程及相变研究已开展了一系列工作,主要集中于冲击Hugoniot数据测量[5],然而与动态加载下物理、力学特性紧密相关的相变研究尚处于起步阶段。声速作为应力扰动在材料中传播的定量表征参量,是获知材料动态响应特性(如相变、屈服强度、剪切模量)的主要途径之一[67],然而相关的研究工作却鲜见报道。

    本研究拟采用反向碰撞实验技术[89],结合具有高时空分辨率的全光纤激光干涉测速技术DPS(Doppler Pin System)[10],对MB2合金开展30~73 GPa压力范围内的冲击Hugoniot及声速测量实验,并与早期实验数据进行对比验证,分析MB2合金的冲击熔化行为。

    • 本研究涉及低压力区的声速测量,为此采用反向碰撞实验设计,即将样品材料制作成飞片,撞击LiF透明窗口,其原理如图1所示。在拉氏坐标下,飞片以速度W直接撞击窗口(t=t1),在飞片和窗口中分别产生左行冲击波和右行冲击波,引起飞片/窗口界面粒子速度的突跃。当飞片中的左行冲击波到达后界面时,将反射中心稀疏波,该稀疏波的波速就是材料在冲击压缩下的声速。如果样品材料发生冲击熔化,则中心稀疏波将以单一塑性波的形式在样品内传播,并在飞片/窗口界面处(t=t2)发生卸载,引起粒子速度的下降。如果样品材料没有发生冲击熔化,则中心稀疏波包括传播速度较快的弹性波和传播速度相对较慢的塑性波。当传播速度较快的弹性卸载波到达飞片/窗口界面时(t=t2),界面粒子速度下降,在速度剖面上形成第1个拐点;当塑性卸载波到达飞片/窗口界面时(t=t3),界面粒子速度再次突变,在速度剖面上形成第2个拐点;当后续稀疏波陆续到达飞片/窗口界面时,会导致界面粒子速度的连续下降;如果在卸载过程中样品材料发生相变,则也会在速度剖面上形成拐点。

      图  1  反向碰撞实验示意图

      Figure 1.  Schematic of backward-impact experimental configuration

      根据图2所示的界面连续性条件(其中p为冲击压力,u为粒子速度),结合Rankine-Hugoniot关系[11],可以得到样品内的冲击波速度Ds

      图  2  反向碰撞实验p-u

      Figure 2.  p-u relation for backward-impact experiment

      ${D_{\rm{s}}} = \frac{{{\rho _{0{\rm{w}}}}{D_{\rm{w}}}{u_{\rm{w}}}}}{{{\rho _{0{\rm{s}}}}(W - {u_{\rm{w}}})}}\quad\quad\quad\quad\quad\quad\quad(1)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad$

      式中:D$\rho $W分别为冲击波速度、密度和飞片速度,下标s和w分别对应样品和窗口。如果窗口材料的D-u曲线满足线性关系

      ${D_{\rm{w}}} = {C_{0{\rm{w}}}} + {\lambda _{\rm{w}}}{u_{\rm{w}}}\quad\quad\quad\quad\quad\quad\quad(2)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad$

      式中:C0w${\lambda _{\rm{w}}}$为窗口材料的Hugoniot参数,则(1)式可表示为

      $ {D_{\rm{s}}} = \frac{{{\rho _{0{\rm{w}}}}({C_{0{\rm{w}}}} + {\lambda _{\rm{w}}}{u_{\rm{w}}}){u_{\rm{w}}}}}{{{\rho _{0{\rm{s}}}}(W - {u_{\rm{w}}})}} $

      因此,只需要测得飞片(样品)击靶速度W和窗口的波后粒子速度uw,就可以获得样品的粒子速度usus=Wuw)和对应的冲击波速度Ds。随后,由波系作用(见图1(b))的几何关系可知,样品材料Hugoniot状态的拉格朗日纵波声速CL

      ${C_{\rm{L}}} = \frac{{{D_{\rm{s}}}{h_{\rm{s}}}}}{{{D_{\rm{s}}}{t_{12}} - {h_{\rm{s}}}}} $

      相应的欧拉纵波声速Cl

      $ {C_{\rm{l}}} = \frac{{{D_{\rm{s}}}{h_{\rm{s}}}}}{{{D_{\rm{s}}}{t_{12}} - {h_{\rm{s}}}}}\frac{{{D_{\rm{s}}} - {u_{\rm{s}}}}}{{{D_{\rm{s}}}}} $

      式中:hs为样品厚度,下标1和2对应波剖面上不同的时间点。

      可以看到,反向碰撞法以波剖面测量为基础,波系作用简单,通过波剖面的粒子速度及时间信息得到高压声速,实验数据具有较高的精度,但是由于可供选择的窗口材料种类较少,目前使用的LiF窗口的阻抗较低,导致实验压力范围有限。

    • 实验在中国工程物理研究院流体物理研究所$\varnothing$30 mm二级轻气炮上进行。将MB2合金飞片安装在弹丸上,将弹丸发射至稳定的弹道速度W,并撞击LiF单晶窗口,通过DPS测量飞片击靶速度以及飞片/窗口界面粒子速度,获得待测材料的冲击波速度和声速。为了提高测试界面对DPS入射光的反射效率,避免靶室残留气体对测试的干扰,窗口击靶面上镀1 $ {\text{μ}}{\rm{m}}$厚的铝膜并贴8 $ {\text{μ}}{\rm{m}}$厚铝箔。LiF窗口折射率修正采用Rigg等[12]的公式

      $ {u_{\rm{w}}} = 0.7895 \times {u^{0.9918}} $

      式中:uw代表修正后的窗口界面粒子速度,u代表实测界面粒子速度,二者单位均为km/s。其中LiF密度为2.638 g/cm3,冲击波速度D和粒子速度u的关系为D=5.150+1.352u(单位km/s)[5]

      图3给出了6发实验测得的MB2/LiF界面粒子速度剖面(平台高度随加载压力的增大而升高)。通过界面粒子速度剖面,由(1)式~(6)式可得MB2样品在30~73 GPa冲击压力范围内的冲击波速度-粒子速度和声速-压力数据,实验结果列于表1,其中密度采用排水法测量,实测值为(1.775±0.004) g/cm3。冲击波速度-粒子速度数据如图4所示,图中还显示了美国洛斯阿拉莫斯国家实验室(LASL)发表的实验数据[5]。从图4中可以看出,本研究获得的实验数据与已有实验数据具有较好的一致性。

      Exp.No. hs/mm W/(km·s–1) uw/(km·s–1) us/(km·s–1) Ds/(km·s–1) p/GPa Cl/(km·s–1)
      1 1.980±0.004 3.949±0.020 1.589±0.016 2.360±0.026 7.303±0.166 30.6±0.4 7.983±0.297
      2 1.997±0.004 4.358±0.020 1.763±0.018 2.595±0.027 7.607±0.175 35.0±0.5 9.101±0.394
      3 1.978±0.004 5.379±0.027 2.195±0.022 3.184±0.030 8.317±0.191 47.0±0.7 9.167±0.402
      4 1.981±0.004 5.928±0.030 2.435±0.024 3.493±0.039 8.746±0.214 54.2±0.9 8.912±0.380
      5 1.981±0.004 6.100±0.030 2.514±0.025 3.586±0.039 8.907±0.213 56.7±0.9 8.601±0.348
      6 1.987±0.004 7.220±0.036 3.003±0.030 4.217±0.050 9.747±0.245 73.0±1.2 9.723±0.463

      表 1  MB2样品冲击实验参数及结果

      Table 1.  Shock experiment parameters and results of MB2

      图  3  MB2/LiF界面粒子速度剖面

      Figure 3.  Particle velocity profile of MB2/LiF interface

      图  4  冲击波速度与粒子速度的关系

      Figure 4.  Shock velocity vs. particle velocity

      图3可以看到:6发实验测得的界面粒子速度剖面质量良好,粒子速度剖面对应的冲击波、卸载波到达样品/窗口界面的特征信号清晰;加载压力为30.6和35.0 GPa时卸载剖面的弹-塑性特征明显,表明MB2合金在该冲击压力下尚未完全熔化;随着加载压力的升高,卸载剖面的弹-塑性特征逐渐消失,当加载压力达到73.0 GPa时,弹-塑性卸载特征完全消失,表明MB2已完全进入熔化相区,与图5所示的不同加载压力下声速转变特征一致。

      图  5  声速与冲击压力的关系

      Figure 5.  Sound velocity vs. shock pressure

      图5所示,当加载压力由30.6 GPa增加到35.0 GPa时,纵波声速逐渐增大,由7.983 km/s增大至9.101 km/s;但是,当加载压力增大至47.0 GPa时,纵波声速逐渐向体波声速偏转,跃变为9.167 km/s,预示着随着加载压力的升高,材料内部的剪切效应减小,冲击熔化发生;直至压力达到56.7 GPa时,纵波声速转变为体波声速,由此进一步确认了冲击加载下MB2合金发生熔化。Urtiew等[13]通过理论预测MB2合金在57 GPa附近开始熔化,该结果与本研究根据声速判定的冲击熔化区域基本一致,从而进一步证实了MB2合金在该压力范围内发生冲击熔化,只是Urtiew等预估的理论压力略高。图5中的实线是根据以下公式[11]计算得到的体波声速曲线

      $ C_{\rm{b}}^2 = \frac{1}{2}{V^2}{\left(\frac{\gamma }{V}\right)_{\rm{H}}}{p_{\rm{H}}} - {V^2}\frac{{{\rm{d}}{p_{\rm{H}}}}}{{{\rm{d}}V}}\left[1 - \frac{1}{2}{\left(\frac{\gamma }{V}\right)_{\rm{H}}}({V_0} - {V_{\rm{H}}})\right] $

      式中:Cb为体波声速;$\rho $$\gamma $分别为密度和Grüneisen系数,$\rho \gamma = {\rho _0}{\gamma _0}$$\rho _0$为材料初始密度,$\gamma _0$=1.43[14]图5中的虚线是基于吴-经方程[15]计算得到的体波声速曲线。可以看出,理论计算的体波声速曲线较冲击加载实验值偏低约10%。这主要是由于合金材料自身成键及结合能的影响,难以采用传统的混合法则或物理模型准确计算物态方程的基本参数,如Grüneisen系数等,由此导致理论预测结果与实验结果出现差异(通常理论值偏低)。

      根据不确定度传递关系[16],当全部直接测量量(输入量)Yi彼此独立不相关时,由其确定的间接测量量z的合成不确定度${\text{Δ}} {\textit{z}}$由下式确定

      ${\text{Δ}} {{\textit{z}}^2} = \sum\limits_{i = 1}^N {{{\left( {\frac{{\partial f}}{{\partial {y_{{i}}}}}} \right)}^2}{\text{Δ}} y_{{i}}^2} $

      式中:yi为输入量Yi的直接测量值,z为被测量的测量值,fzyi的函数关系,N为输入量的总个数,${\text{Δ}} {y_{{i}}}$为直接测量值的不确定度。因而,当密度为${\rho _{0{\rm{s}}}}$、厚度为${h_{\rm{s}}}$的飞片撞击窗口时,样品/窗口界面处粒子速度跳跃,实验测得飞片速度$W$、界面粒子速度${u_{\rm{w}}}$、时间间隔(剖面平台)${t_{12}}$,相应的测量不确定度为${\text{Δ}} {\rho _{0{\rm{s}}}}$${\text{Δ}} {h_{\rm{s}}}$${\text{Δ}} W$${\text{Δ}} {u_{\rm{w}}}$${\text{Δ}} {t_{12}}$,而窗口Hugoniot参数(${\rho _{0{\rm{w}}}}$${C_{0{\rm{w}}}}$${\lambda _{\rm{w}}}$)的不确定度(${\text{Δ}} {\rho _{0{\rm{w}}}}$${\text{Δ}} {C_{0{\rm{w}}}}$${\text{Δ}} {\lambda _{\rm{w}}}$)已知。由于各测量量是独立测量的,(${C_{0{\rm{w}}}}$${\lambda _{\rm{w}}}$)相互作用项较小,可忽略不计,据此可根据不确定度传递律确定声速测量的不确定度。

      图6所示,就该例反碰撞实验(No.2)而言,影响声速测量不确定度的因素很多,包括样品初始密度、厚度、飞片速度、界面粒子速度、追赶时间、窗口材料冲击Hugoniot参数等。样品内部冲击压缩状态(如粒子速度、冲击波速度、冲击压力等)均通过飞片速度、界面粒子速度及窗口Hugoniot参数计算获得,影响声速测量不确定度的主要因素在于飞片速度、界面粒子速度及稀疏波追赶时间(平台时间),所占比例(即对(8)式中各平方项求和,下同)约为声速测量不确定度的99%,其余如初始密度、厚度等参量在当前诊断水平下对声速测量不确定度的贡献较小,约占总体的1%。在现有诊断条件下,飞片速度采用DPS直接测量,测量扩展不确定度不大于0.5%;界面粒子速度剖面的测量不确定度主要受平台区速度及稀疏波追赶时间的影响,影响因素包含干涉信号数据转换精度、窗口折射率修正、起跳及卸载时刻的判断等,综合而言,平台区速度测量的扩展不确定度不大于1%,追赶时间测量的扩展不确定度约6 ns。总体而言,传递至声速的测量扩展不确定度不超过5%。

      图  6  反碰撞法测量声速实验中不确定度

      Figure 6.  Uncertainties for backward-impact experiment based on law of propagation

    • 采用反向碰撞实验技术,结合具有高时空分辨率的DPS,获得了MB2合金在30~73 GPa压力范围内的冲击Hugoniot及声速数据。随着加载压力的升高,MB2合金纵波声速呈现出明显的向体波声速转变的趋势,预示着材料内部的剪切效应逐渐减小,冲击熔化发生,其相变压力区间为40~57 GPa。该实验结果与不同加载压力下卸载波剖面对应的弹-塑性转变特征完全一致,由此进一步确认了冲击加载下MB2合金熔化相变的发生。

      感谢中国工程物理研究院流体物理研究所黄金、康强、叶素华、方茂林、向曜民、陈志云等在实验过程中给予帮助。

参考文献 (16)

目录

    /

    返回文章
    返回