SHPB试验应变率设计.docVIP

SHPB试验的应变率设计 朱耀，庞宝君，盖秉政，王立闻 （哈尔滨工业大学航天学院，哈尔滨，150001） 摘 要：根据一维应力波理论和二波法计算公式推导出SHPB试验的应变率预算公式，根据该公式可由试验条件和试件尺寸预先计算出试验的应变率，并分析了SHPB试验时应变率限制范围。利用LS-DYNA进行了二维SHPB数值仿真试验，试验结果证明该应变率计算公式准确可靠，对SHPB试验的方案设计具有指导意义。 关键词： SHPB；应变率；数值仿真 0 引言 分离式霍普金森压杆（Split Hopkinson Pressure Bar，简称SHPB）是用于研究材料在102 - 104s-1应变率下力学行为的最主要的试验手段。在Davies[1]完成了该项技术的关键分析后，Kolsky[2]对该技术进行了革命性改进，之后SHPB相关的理论分析与试验手段基本成熟。但是，关于影响试验的精度的因素还是得到了许多研究者的关注，而这些因素中能够调整与改进的主要都与与试样尺寸有关，包括以下几个方面：试件与压杆横截面不匹配引起的二维效应[3,4]；压杆与试件之间的端面摩擦效应[5,6]；试件的动态惯性效应[7-10]等。在进行动态试验时，除试件尺寸外，应变率是一个关键的试验参数。本文根据一维应力波理论和二波法计算公式推导出SHPB试验的应变率预算公式，对SHPB试验方案设计具有的指导意义。 1 SHPB试验的应变率 图1 分离式霍普金森压杆SHPB系统示意图 Fig.1 Schematic illustration of split Hopkinson pressure bar system 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 材料参数 Table 1 Material parameters 材料 E / GPa μ ρ / kg/m3 σ0 / MPa Ep / MPa 试件 71 0.3 2780 370 1000 材料 E / GPa μ ρ / kg/m3 压杆 210 0.3 7850 压杆直径为12mm，试件尺寸见表2，因为在此数值模拟仅为了检验应变率预测公式式（9），所以对数值模拟的细节方面不做介绍，仅列出仿真试验结果与式（9）计算结果的对比，见表2，计算结果与仿真结果都是平均应变率,其中计算仿真结果平均应变率时去掉了开始阶段的上升沿以及结