SHPB试验的应变率设计.doc

宇航学报第卷
宇航学报
哈尔滨工业大学航天学院博士生学术论坛
--
SHPB试验的应变率设计
朱耀，庞宝君，盖秉政，王立闻
（哈尔滨工业大学航天学院，哈尔滨，150001）
摘要：根据一维应力波理论和二波法计算公式推导出SHPB试验的应变率预算公式，根据该公式可由试验条件和试件尺寸预先计算出试验的应变率，并分析了SHPB试验时应变率限制范围。利用LS-DYNA进行了二维SHPB数值仿真试验，试验结果证明该应变率计算公式准确可靠，对SHPB试验的方案设计具有指导意义。
关键词：SHPB；应变率；数值仿真

第卷
SHPB试验的应变率设计
--
0引言
分离式霍普金森压杆（SplitHopkinsonPressureBar，简称SHPB）是用于研究材料在102-104s-1应变率下力学行为的最主要的试验手段。在Davies[1]完成了该项技术的关键分析后，Kolsky[2]对该技术进行了革命性改进，之后SHPB相关的理论分析与试验手段基本成熟。但是，关于影响试验的精度的因素还是得到了许多研究者的关注，而这些因素中能够调整与改进的主要都与与试样尺寸有关，包括以下几个方面：试件与压杆横截面不匹配引起的二维效应[3,4]；压杆与试件之间的端面摩擦效应[5,6]；试件的动态惯性效应[7-10]等。在进行动态试验时，除试件尺寸外，应变率是一个关键的试验参数。本文根据一维应力波理论和二波法计算公式推导出SHPB试验的应变率预算公式，对SHPB试验方案设计具有的指导意义。
1SHPB试验的应变率
子弹
入射杆
试件
透射杆
应变片
超动态应变仪
波形存储器
计算机
图1分离式霍普金森压杆SHPB系统示意图
Fig.1SchematicillustrationofsplitHopkinsonpressurebarsystem
SHPB装置如图1所示。首先，子弹以一定的速度撞击入射杆，在入射杆中产生一个压缩脉冲，即入射波。该压缩脉冲通过入射杆传播到试件处，试件在该压缩载荷作用下高速变形,同时压缩脉冲一部分反射回入射杆中形成反射波，另一部分通过试件传入透射杆中形成透射波。三个波可以通过贴在杆上的应变片记录下来。所有的SHPB试验都是通过入射波、反射波和透射波这三个波的数据来进行计算处理，包括三波法、二波法以及它们的改进算法。相对于三波法，二波法更为常用，且形式较三波法简单易用，二波法具体形式如下
（4）
（5）
（6）
式中，是试件中的平均应变率，是试件中的平均应变，为试件中的平均应力；和分别为反射波和透射波对应的应变，由贴在入射杆和透射杆上的应变片可以直接读出；C0和A0分别是压杆的声速与横截面积；Ls和As分别是试件的初始长度与横截面积。
由于子弹与入射杆材料、横截面积都相同，子弹撞击入射杆产生的入射脉冲应力幅值为
（7）
式中，v0为子弹撞击速度，可以通过试验测得。
根据SHPB试验的试件应力/应变均匀化假定再根据一维应力波理论有
即（8）
根据式（4）—式（8）有

（9）

式中，分别为压杆的密度、波速、弹性模量及压杆直径，分别为试件的长度、直径与工程应力，为子弹撞击速度。当压杆试验装置以及测试材料确定之后，这几个参数就已经确定了，只能通过调节子弹撞击速度和试件的尺寸来实现不同的应变率。
1.1应变率限制
在实际的SHPB试验操作中，子弹的撞击速度v0存在以下几个限制条件：
1、压杆的弹性限制。为确保试验过程中，压杆始终处于弹性状态，子弹撞击速度必须低于屈服速度vY，即
（10）
2、反射波强度不能过小。原因有两个：一是确保入射波强度能够使试件进入塑性屈服阶段；另外就是在利用反射波计算应变的时候确保误差不至于过大。在此假设反射波强度必须达到入射波强度的1/10：
即（11）
将式（8）代入式（11）得
（12）
将式（6）和式（7）带入式（12）得
（13）
（14）
联立式（10）和式（14）可得子弹的撞击速度调节范围为
（15）
将上式带入（9）式可得
（16）
2数值模拟仿真实验
数值模拟采用Ansys/LS-DYNA有限元程序建立二维轴对称模型，压杆材料采用线弹性材料，而试件材料采用双线性各向同性硬化模型，表1给出了具体的材料参数。
表1材料参数
Table1Materialparameters材料E/GPaμρ/kg/m3σ0/MPaEp/MPa试件710.327803701000材料E/GPaμρ/kg/m3压杆2100.37850压杆直径为12mm，试件尺寸见表2，因为在此数值模拟仅为了检验应变率预测公式式（9），所以对数值模拟的细节方面不做介绍，仅列出仿真试验结果与式（9）计算结果的对比，见表2，计算结果与仿真结果都是平均应变率,其中计算仿真结果平均应变率时去掉了开始阶段的上升沿以及结尾的下降沿。
表2仿真结果与计算结果对