SHPB试验的应变率设计

哈尔滨工业大学航天学院博士生学术论坛

SHPB试验的应变率设计

朱耀，庞宝君，盖秉政，王立闻

（哈尔滨工业大学航天学院，哈尔滨，150001）

摘 要：根据一维应力波理论和二波法计算公式推导出SHPB试验的应变率预算公式，根据该公式可由试验条件和试件尺寸预先计算出试验的应变率，并分析了SHPB试验时应变率限制范围。利用LS-DYNA进行了二维SHPB数值仿真试验，试验结果证明该应变率计算公式准确可靠，对SHPB试验的方案设计具有指导意义。 关键词： SHPB；应变率；数值仿真

0 引言

分离式霍普金森压杆（Split Hopkinson Pressure Bar，简称SHPB）是用于研究材料在24-1

10 - 10s应变率下力学行为的最主要的试验手

[1]

段。在Davies完成了该项技术的关键分析后，

[2]

Kolsky对该技术进行了革命性改进，之后SHPB相关的理论分析与试验手段基本成熟。但是，关于影响试验的精度的因素还是得到了许多研究者的关注，而这些因素中能够调整与改进的主要都1 SHPB试验的应变率

与与试样尺寸有关，包括以下几个方面：试件与

[3,4]

压杆横截面不匹配引起的二维效应；压杆与试

[5,6]

件之间的端面摩擦效应；试件的动态惯性效应[7-10]

等。在进行动态试验时，除试件尺寸外，应变率是一个关键的试验参数。本文根据一维应力波理论和二波法计算公式推导出SHPB试验的应变率预算公式，对SHPB试验方案设计具有的指导意义。

子弹 入射杆 试件 透射杆 应变片 超动态应变仪 波形存储器 计算机 图1 分离式霍普金森压杆SHPB系统示意图 Fig.1 Schematic illustration of split Hopkinson pressure bar system SHPB装置如图1所示。首先，子弹以一定的速度撞击入射杆，在入射杆中产生一个压缩脉冲，即入射波。该压缩脉冲通过入射杆传播到试件处，试件在该压缩载荷作用下高速变形,同时压缩脉冲一部分反射回入射杆中形成反射波，另一部分通过试件传入透射杆中形成透射波。三个波可以通过贴在杆上的应变片记录下来。所有的SHPB试验都是通过入射波、反射波和透射波这三个波的数据来进行计算处理，包括三波法、二

- 1 -

波法以及它们的改进算法。相对于三波法，二波法更为常用，且形式较三波法简单易用，二波法具体形式如下

?s??.2C0?rLs（4） （5） （6）

?s??2C0Ls??d?0rt?s?EA0?tAsSHPB试验的应变率设计 ?s是试件中的平均应变率，?s是试件中式中，?的平均应变，?s为试件中的平均应力；?r和?t 分别为反射波和透射波对应的应变，由贴在入射杆和透射杆上的应变片可以直接读出；C0和A0分别是压杆的声速与横截面积；Ls和As分别是试件的初始长度与横截面积。

由于子弹与入射杆材料、横截面积都相同，子弹撞击入射杆产生的入射脉冲应力幅值为

1（7） ?i???0C0v02式中，v0为子弹撞击速度，可以通过试验测得。 根据SHPB试验的试件应力/应变均匀化假定再根据一维应力波理论有

9?i?10?t（12）

将式（6）和式（7）带入式（12）得

?Ds?9?Cv?000?10?s??D??2?0?2（13） （14）

20?sDsv0??29?0C0D02联立式（10）和式（14）可得子弹的撞击速度调节范围为

22?Y20?sDs?v??（15） 2?0C009?0C0D0将上式带入（9）式可得

22???2C0?D2C0??sDs?s?????????????s?Ys?2?（16） ????lsElsE?9D0??D0????i??r??t（8） ?i??r?? t即

根据式（4）—式（8）有

2C?2C0(?i??t)2C0(?i??tE)?s??0r???LsLsLsE

2?Ds??2C0?1??????0C0v0??s???LsE?2?D0????2 数值模拟仿真实验

数值模拟采用Ansys/LS-DYNA有限元程序建

立二维轴对称模型，压杆材料采用线弹性材料，而试件材料采用双线性各向同性硬化模型，表1给出了具体的材料参数。 表1 材料参数 Table 1 Material parameters 材料 E / GPa 试件 材料 压杆 71 μ 0.3 E / GPa 210 （9）

式中，?0,C0,E,D0分别为压杆的密度、波速、弹性模量及压杆直径，Ls,Ds,?s分别为试件的长度、直径与工程应力，v0为子弹撞击速度。当压杆试验装置以及测试材料确定之后，?0,C0,E,D0,?s这几个参数就已经确定了，只能通过调节子弹撞击速度v0和试件的尺寸Ls,Ds来实现不同的应变率。 1.1 应变率限制

在实际的SHPB试验操作中，子弹的撞击速度v0存在以下几个限制条件：

1、 压杆的弹性限制。为确保试验过程中，压杆始终处于弹性状态，子弹撞击速度必须低于屈服速度vY，即

2?Yv0?vY?（10）

?0C02、 反射波强度不能过小。原因有两个：一是确保入射波强度能够使试件进入塑性屈服阶段；另外就是在利用反射波计算应变的时候确保误差不至于过大。在此假设反射波强度必须达到入射波强度的1/10：

?r??i/10?r??i/10 即 （11）

ρ / kg/m3 σ0 / MPa Ep / MPa 2780 370 μ 0.3 1000 ρ / kg/m3 7850 压杆直径为12mm，试件尺寸见表2，因为在

此数值模拟仅为了检验应变率预测公式式（9），所以对数值模拟的细节方面不做介绍，仅列出仿真试验结果与式（9）计算结果的对比，见表2，计算结果与仿真结果都是平均应变率,其中计算仿真结果平均应变率时去掉了开始阶段的上升沿以及结尾的下降沿。 表2 仿真结果与计算结果对比 Table 2 Contrast between simulation and calculation 编号 1 2 3 4 5 6 试件直试件长子弹速计算结仿真结误差 径(mm) 度(mm) 度(m/s) 果(s-1) 果(s-1) 5 8 10 3 5 8 5 5 5 5 5 5 30 30 30 20.6 24.5 4804.6 4799.6 0.10% 3380 3349.1 0.93% 2396.1 2334.9 2.62% 3749 3782.6 -0.89% 3833.4 3873.7 -1.04% 将式（8）代入式（11）得

- 2 -

34.15 3997.4 3990.5 0.17% SHPB试验的应变率设计 7 10 5 43 4065.2 3954.3 2.80% 表2中编号1-3是同一撞击速度下不同尺寸试件的仿真试验；编号4-7是预先利用式（9）简单估算得到不同试件尺寸同一应变率下的仿真试验，具体过程为：根据式（1）计算得到直径10mm的试件最大工程应变不能超过0.3056，压缩过程0作用时间t0?2L0/C0，从而计算出最大平均应变率为3951s-1，根据最大工程应变可以计算出塑性段平均工程应力，将平均应变率3951s-1和平均工程应力以及试件尺寸代入应变率预算公式式（9）中，即可计算得到相应的撞击速度，然后建立有限元模型进行仿真试验。

0-1000-2000-3000-4000-5000-60000.00-1000Strain rate (1/s)-2000 Calculation Simulation-3000-4000-50000.0015.00μ30.00μ45.00μ60.00μ75.00μ90.00μStrain rate (1/s) Calculation SimulationTime (s)0-1000-2000 15.00μ30.00μ45.00μ60.00μ75.00μ90.00μTime (s)0-1000 35302520151015.00μ30.00μ45.00μ60.00μ Calculation Simulation-2000 Calculation SimulationStrain rate (1/s)-3000-4000-5000-6000-70000.00Strain rate (1/s)-3000-4000-5000-6000-7000velocity / m/s15.00μ30.00μ45.00μ60.00μ75.00μ90.00μTime (s) -80000.00 300mm strike bar 200mm strike bar75.00μ90.00μ5Time (s) 图2 计算应变率与仿真试验应变率的对比 00.00.2Fig.2 Contrast between calculation and simulation 图2(a)-(d)是4-7号仿真试验的计算应变率与仿真试验应变率的对比图，由图2可以发现计算结果与仿真试验结果符合的比较好。4-7号仿真试验的设计应变率是一样的，仿真结果得到的平均应变率之间的最大误差为4.34%，同时，所有仿真试验结果与计算结果的最大误差为2.80%，说明利用式（9）计算应变率是比较准确的。 3 AA7055铝合金动态压缩试验

AA7055铝合金动态压缩试验要求应变率为1000、2000、3000、4000和5000s-1。根据预试试验得到AA7055屈服强度约为700MPa，且应变硬化比较小，为了计算方便，在根据式（9）进行应变率预算时取?s为700MPa。在SHPB实验中是通过调节充气压力来改变子弹速度的，我们进行试验的Ф12.7SHPB的充气压力与子弹速度关系如图3所示。考虑到充气压力过高对实验装置损害较大，所以限制充气压力为0.8MPa。综上所述，根据第1节的分析，我们确定试验方案如表3所示。

- 3 -

350.40.60.81.01.2Charge Pressure / MPa

velocity / m/s30 20 15 10 300mm strike bar 200mm strike bar5 00.00.20.40.60.81.01.2

Charge Pressure / MPa

图2 计算应变率与仿真试验应变率的对比 Fig.2 Contrast between calculation and simulation 25表2 仿真结果与计算结果对比 Table 2 Contrast between simulation and calculation 编号 试件直试件长子弹速计算结仿真结误差 径(mm) 度(mm) 度(m/s) 果(s-1) 果(s-1)