An Instantaneous Optical Pyrometer with Six Channels for the Shock Temperature Measurement in Materials

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

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

1.   实　验

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

图  1  反向碰撞实验示意图

Figure  1.  Schematic of backward-impact experimental configuration

下载: 全尺寸图片 幻灯片

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

图  2  反向碰撞实验p-u图

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

下载: 全尺寸图片 幻灯片

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

式中：C_0w和${\lambda _{\rm{w}}}$为窗口材料的Hugoniot参数，则（1）式可表示为

因此，只需要测得飞片（样品）击靶速度W和窗口的波后粒子速度u_w，就可以获得样品的粒子速度u_s（u_s=W–u_w）和对应的冲击波速度D_s。随后，由波系作用（见图1（b））的几何关系可知，样品材料Hugoniot状态的拉格朗日纵波声速C_L为

相应的欧拉纵波声速C_l为

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

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

2.   实验结果与分析

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

式中：u_w代表修正后的窗口界面粒子速度，u代表实测界面粒子速度，二者单位均为km/s。其中LiF密度为2.638 g/cm³，冲击波速度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/cm³。冲击波速度-粒子速度数据如图4所示，图中还显示了美国洛斯阿拉莫斯国家实验室（LASL）发表的实验数据^[5]。从图4中可以看出，本研究获得的实验数据与已有实验数据具有较好的一致性。

表  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 

| 显示表格

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

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

式中：y_i为输入量Y_i的直接测量值，z为被测量的测量值，f为z和y_i的函数关系，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

下载: 全尺寸图片 幻灯片

3.   结　论

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

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