Bi在低压段区Hugoniot数据的精确测量.docxVIP

Bi在低压段区Hugoniot数据的精确测量 1. 引 言 　　雨贡纽( Hugoniot) 数据是材料从同一初始状态出发、经过不同的冲击压缩达到的终态的集合，描述了冲击波后任一压缩态下热力学状态量之间的内在联系． 一方面，通过 Hugoniot 数据的测量，能够研究材料在高压下的相变和物态方程; 另一方面，已知 Hugoniot 数据也是进行冲击波实验设计的基本材料参数．  　　  　　Bi 在冲击波加 载 下 会 出 现 较 复 杂 的 相 变 序列，是进行冲击相变研究的典型材料; 而其低压区复杂的多形相变和反常熔化，使得该压力段内Hugoniot 数据的测量对研究材料弹塑性响应和冲击相变具有重要意义． 早期 Duff 等就利用基于电探针技术的自由面速度法测量了 Bi 的低压 Hugoniot数据，并在 2. 5 GPa 附近观测到了 BiⅠ( 斜方六面体结构)rarr;BiⅡ( 底心单斜结构) 相变，相变压力较静压的结果偏高． 之后 Romain利用电探针技术测量了 Bi 在 20 GPa 压力下的 Hugoniot 数据，实验发现在 7 GPa 附近的冲击压缩下 Bi 发生了 BiⅢ( 四方晶系结构)rarr;BiⅤ( 体心立方结构) 相变． 随后，国外相继利用多种加载手段对 Bi 的冲击相变热力学、动力学以及多相物态方程进行了研究． 但作为上述研究基础工作的 Bi 的低压区 Hugoniot 数据，却存在较大的分散性，由此导致获得的相边界的不确定度较大，对 Bi 低压段多相物态方程的建立存在较大的影响． 因而有必要对 Bi 在低压段的 Hugoniot数据进行精确测量．  　　  　　目前人们还不能完全从理论的方法确定材料的 Hugoniot 数据，必须借助于实验测量． 而现有的Hugoniot 数据测量中，冲击波速度 D 主要利用电探针技术测量冲击波在已知厚度样品中的走时( 又称渡越时间) 得到． 由于飞片在加速过程中不可避免地会产生倾斜和弯曲变形，而样品内也存在复杂的波系作用，因此必须要对由此引起的冲击波走时进行修正，该修正量不确定度的引入，导致了冲击波速度合成不确定度的加大; 同时，电探针法作为一种接触式的测量方法，要求了电探针的有效触发和快时间响应，然而电探针的低压触发响应较慢，这使得材料低压区 Hugoniot 数据的精密测量变得更加困难．  　　  　　本研究工作采用反向碰撞的实验技术，探索了通过样品 /窗口界面粒子速度剖面的测量获得材料低压区 Hugoniot 数据的可行性． 利用该方法，获得了 Bi 在 1045 GPa 压 力 范 围 内 的 高 精 度Hugoniot 数据，并与已有的文献报道数据进行了比较．  　　  　　2. 实验方法 　　对于静止靶初始状态为常压的情形，由平面冲击波阵面前后状态满足质量、动量与能量守恒关系得到 Rankine-Hugoniot公式【1-3】    　　  　　式中 rho;，D，u，p 和 E 分别表示冲击波阵面后的密度、冲击波速度、粒子速度、压力和比内能，下标 0 表示初始状态．目前还没有直接测量比内能 E 的方法，而压力p 及密度 rho; 测量的精度较低，因此，Hugoniot 数据测量实验通常是对冲击波速度 D 和粒子速度 u 进行测量． 对于大多数凝聚介质，D-u 曲线在一定范围内满足下列线性关系:【4】    　　  　　式中 C0和 lambda; 为常数，即 Hugoniot 参数． 当前发展的实验方法主要有绝对测量法和对比测量法． 在上述两种方法中，冲击波速度 D 都是通过测量冲击波在一定厚度样品中的走时得到，其精度主要取决于冲击波走时的精度． 为提高冲击波走时的测量精度，需要控制入射冲击波的变形程度( 弯曲和倾斜等) ，并在样品前后界面设计合理的探针布局，通过一定的统计分析模型对冲击波走时进行修正． 在绝对测量法中，由于飞片和样品为同种材料，而且初始状态也完全相同，冲击波后粒子速度 u 通过测量飞片速度 W 得到( u = W /2) ; 在对比测量法中，冲击波后粒子速度 u 无法直接测量得到，必须通过相关量的间接测量，利用阻抗匹配法获得． 由于飞片和样品材料阻抗不匹配，粒子速度的不确定度还需考虑飞片材料 Hugoniot 镜像线近似二次冲击绝热线或卸载线而引入的不确定度．本研究采用反向碰撞法测量 Hugoniot 数据，其原理如图 1． 如图 1( a) 的 X-t 图所示，在拉格朗日坐标系中，当作为飞片的样品以速度 W 撞击静止的透明窗口时，将分别在样品和窗口中产生左、右行冲击波，并相应地使样品和窗口中的压力、粒子速度等状