Adiabatic Shear Failure Behavior of 30CrMnMo Steel under Pulse Stress Impact
-
摘要: 为研究30CrMnMo钢在脉冲应力冲击载荷下的绝热剪切失效及演化特性,利用分离式霍普金森压杆对一种轴对称帽型试件进行冲击剪切实验,并运用LS-DYNA动力学有限元软件对不同入射脉冲应力载荷下的剪切失效演化及剪切区温度分布进行数值模拟。结果表明,帽型试件的绝热剪切失效与脉冲应力比冲量相关,对于30CrMnMo钢帽型试件,其绝热剪切失效对应的脉冲应力比冲量近似为常量。数值模拟中,当网格尺寸小于剪切带宽度时,能够有效模拟剪切带内的局部温升热点特性。绝热剪切演化表现为失稳由帽型试件剪切区拐角处启动并同时向中心传播,剪切带内外材料主要经历均匀剪切变形和失稳快速扩展2个阶段。Abstract: In order to study the adiabatic shear failure behavior and evolution characteristics of 30CrMnMo steel under pulse stress impact, a split Hopkinson pressure bar was used to conduct an axisymmetric cap shaped specimen for impact shear experiments. The shear failure evolution and temperature distribution in the shear zone under different incident pulse stress loads were numerically simulated using LS-DYNA dynamic finite element software. The results indicate that the adiabatic shear failure of the cap shaped specimen is related to the specific impulse of the pulse stress. For cap shaped specimen of 30CrMnMo steel, the specific impulse of pulse stress corresponding to the adiabatic shear failure is approximately constant. In numerical simulation, when the grid size is smaller than the width of the shear band, the local temperature rise of hot spot within the shear band can be effectively simulated. The evolution of adiabatic shear instability is characterized by simultaneous propagation from the corner of the shear zone to the center, and the materials inside and outside the shear zone mainly undergo two stages: uniform shear deformation and rapid expansion of instability.
-
软质高分子聚合物材料因生产成本低、质量轻、具有良好的抗冲击性能等优点在军工、汽车、航空、电子[1-4]等行业得到广泛应用,如用于汽车碰撞实验的假人皮肤材料[5]、航空航天仿生材料[6]等。在上述领域,构件/材料常受动态冲击载荷作用,明确其动态力学性能,构建描述其黏弹性力学特征的本构模型,可为软质高分子聚合物材料的设计开发、性能优化及可靠性分析提供理论模型和方法支撑。
学者们曾利用分离式霍普金森压杆(Split Hopkinson pressure bar, SHPB)装置对PVC弹性体等软质高分子聚合物材料开展了动态压缩力学性能研究[7-8]。与SHPB相比,分离式霍普金森拉杆(Split Hopkinson tensile bar, SHTB)装置采用拉伸加载形式,便于开展材料的动态冲击损伤、变形、断裂等行为研究。近年来,研究人员已对丁氰橡胶[9]、聚脲[10]、半晶聚合物高密度聚乙烯(HDPE)[11]、猪皮[12]、聚氨酯弹性体[13]等低阻抗软材料进行了静、动态拉伸力学性能测试。同时,自动网格法应变测量[14]、脉冲整形[15]等技术以及高/低温环境箱[16]、扫描电镜[17-18]、高速摄像机[19]等装置也被用于软材料的拉伸变形、裂纹萌生及扩展机理分析。通过静、动态拉伸实验发现:软质高分子聚合物材料的拉伸强度会随着应变率对数的增大而线性增加,且不同分子链结构的软质高分子聚合物在同一高拉伸应变率下有不同的拉伸强度和拉伸应变[20]。
但是,软质高分子聚合物材料动态拉伸实验过程中仍存在试样连接方式难以确定、胶黏剂黏接强度低、加载应变率不恒定等问题。为揭示PVC弹性体材料的静、动态拉伸力学性能,利用Instron-5943万能材料试验机和改进型SHTB实验装置对PVC弹性体进行静、动态直接拉伸实验,分析应变率对PVC弹性体材料静、动态拉伸力学性能的影响规律。
1. 实验材料及测试方法
1.1 准静态拉伸实验
实验材料为PVC弹性体。准静态拉伸试样根据GB/T 528–2009[21]中的3型试样标准设计由模具冲压而成,具体尺寸见图1,试样呈哑铃状,厚度为2 mm。由于弹性体材料动态拉伸实验试样的设计尚无统一标准,因此本研究中动态拉伸实验试样将依据动态拉伸实验的具体效果进行设计并制备,详见2.3.1节。
PVC弹性体的准静态拉伸力学性能实验是在Instron-5943万能材料试验机上进行的。实验时,将试验机的加载速率设置为96 mm/min,使其应变率为0.1 s−1。拉伸过程中的载荷(F)和位移(ΔL)数据通过传感器输入至计算机中,利用式(1)、式(2)将其转化为应力、应变数据
σ=FAs (1) ε=ΔLLs (2) 式中:
σ 为应力,As为试样截面积,ε 为应变,Ls为试样长度。1.2 动态拉伸实验
采用合肥姜水动态力学实验技术有限公司第一实验室的分离式霍普金森拉杆(SHTB)实验系统开展聚氯乙烯(PVC)弹性体材料在高应变率下的动态拉伸实验,其装置原理及实物见图2。
调节储气瓶阀向储气室充气,待储气室内气压到达预定值时,关闭储气瓶阀,触发发射按钮,使储气室内高压气体流入发射通道,推动管状套筒子弹在入射杆上加速至一定速度,并以该速度撞击入射杆顶端的法兰盘,在入射杆中产生一个压缩应力波,其持续时间由子弹长度决定。套筒子弹与法兰盘撞击瞬间形成的压缩应力波分为两部分:一部分压缩应力波被吸收杆捕获,最终被能量吸收器(阻尼器)吸收;另一部分压缩应力波在入射杆自由端面处反射形成入射杆中的拉伸应力波。当该波传递至试样端面时,由于杆系与试样间的波阻抗不匹配,一部分拉伸应力波会被反射回入射杆形成反射波,另一部分则会穿过试样并传递至透射杆形成透射波。
根据一维弹性应力波的传播理论,由入射杆和透射杆上的半导体应变片分别获取入射波
εi(t) 、反射波εr(t) 和透射波εt(t) 的应变脉冲信号,可计算出试样两端面的力F1 、F2 及位移u1 、u2 F1=E0A0[εi(t)+εr(t)] (3) F2=E0A0εt(t) (4) u1=C0∫[εi(t)−εr(t)]dt (5) u2=C0∫εt(t)dt (6) 当试样内部满足应力均匀性假设[22]时,材料的应力、应变、应变率具有以下关系
σs=F1+F22As (7) εs=u1−u22Ls (8) ˙εs=dεdt (9) 由于二波法和三波法比较时,三波法的可信度较高,更能反映真实的测试结果,故使用三波法处理实验所得数据。整理式(3)~式(9)可得数据处理公式
˙εs=C0Ls[εi(t)−εr(t)−εt(t)] (10) εs=C0Ls∫t0[εi(t)−εr(t)−εt(t)]dt (11) σs=A0E02As[εi(t)+εr(t)+εt(t)] (12) 式中:
˙εs 、εs 和σs 分别为动态拉伸实验中试样的应变率、应变和应力,εi(t) 、εr(t) 和εt(t) 分别为应变片所获取的电压信号经转换之后的入射应变、反射应变和透射应变,C0 、A0 和E0 分别为杆的弹性波波速、横截面积和弹性模量,t 为应力波在试样内的传播时间。1.3 动态拉伸实验关键问题及解决方法
SHTB实验装置的拉杆均为直径20 mm的实心铝杆,杆系总长度为9500 mm,其中,入射杆长度为5000 mm,透射杆长度为3000 mm,吸收杆长度为1500 mm。法兰盘安装于入射杆上,管状铝质套筒子弹长度为600 mm。软质高分子聚合物材料的SHTB实验与SHPB实验不同,实验过程存在动态拉伸试样连接方式及胶黏剂优选、恒应变率加载与入射波完全卸载等问题,可联合波形分析和高速摄像等方法加以解决。
1.3.1 试样连接方式的确定及胶黏剂优选
目前,SHTB实验中试样的连接方式主要有3种:螺纹连接、挂接以及黏接。由于PVC弹性体材料自身特性不宜采用螺纹连接,因而本研究首先尝试了较常用的挂接方式(使用特殊夹具夹持固定试样)。
挂接夹具和哑铃状试样尺寸如图3所示,由于聚氯乙烯(PVC)弹性体材料的波速较低,难以达到应力均匀,故试样不可过长;同时若试样长度过短,试样两端夹持部分会影响实验结果的准确性。根据橡胶材料测试经验[23],本研究动态拉伸试样有效拉伸长度设置为2 mm。
挂接所得波形见图4。由图4可知:入射脉冲为矩形方波,且在100 μs内便可达到最大应变;透射脉冲较宽且透射信号在0~0.2 μs之间有明显下降点,但始终未下降至水平基线位置。
通过分析发现,出现这一现象主要有两方面原因:一是挂接中使用的夹具对波形传播造成了干扰;二是透射杆长度过短,使得通过应变片的透射波沿透射杆返回,与尚未完全通过应变片的透射波在应变片处发生了叠加。
为解决上述问题,对原有SHTB实验装置进行了改进。改进后的SHTB杆系总长度为11500 mm,其中,入射杆长度为5000 mm,透射杆长度为5000 mm,吸收杆长度为1500 mm。试样的连接方式改为黏接(将试样用合适强度的胶黏剂黏接于入射杆与透射杆之间),试样尺寸为
∅ 20 mm× 2 mm。由于胶黏剂完全固化后才可达到最大黏接强度,为加快实验速度,首先选用了强度高、完全固化时间短(约10 min)的502瞬间强力胶进行了实验,试样黏接示意图及加载波形见图5。
由图5可知,入射波仍为矩形方波,且可在100 μs内达到最大应变;透射脉冲在极短时间内上升至最大应变后迅速下降且回到基线位置。为明确透射波上升至最大应变后立即卸载的原因,采用高速摄像机(V1212,深圳约克科技有限公司)记录了动态拉伸实验过程,发现波传递至试样与入射杆黏接面的瞬间,502瞬间强力胶发生强烈抖动,导致试样与入射杆黏接面脱胶,致使透射波上升至最大应变后立即卸载,高速摄像照片见图6。
为解决脱胶问题,对比了多种类型的胶黏剂(如YH-818专用橡胶胶水、固特灵401胶水、环氧树脂AB胶、JL-330橡胶专用胶等)并进行了胶黏剂优选。研究发现,环氧树脂AB胶和JL-330橡胶专用胶可以将试样牢固地黏接于入射杆与透射杆之间。利用SHTB装置对这两种胶黏剂黏接的试样进行了动态拉伸实验,以确定这两种胶黏剂的适用性,实验所得波形见图7。
由图7可知,采用两种胶黏剂所得入射脉冲均为矩形方波,且可在100 μs内达到最大应变,其中,环氧树脂AB胶所得透射脉冲上升至最大应变后迅速下降并可回到基线,JL-330橡胶专用胶所得透射波加载时间变长(约在200 μs左右),且达到最大应变后缓慢下降,说明JL-330橡胶专用胶可在一定程度上缓解透射波上升至最大应变后立即卸载的问题,故后续实验将尝试用JL-330橡胶专用胶进行聚氯乙烯弹性体动态拉伸实验。
1.3.2 恒应变率加载
当SHTB实验采用直接加载时,入射波为方波且在极短时间内即可达到最大应变,难以保证整个加载过程为恒应变率加载。为确保SHTB实验结果的准确性,在法兰盘端面(套筒子弹撞击端面)上放置脉冲整形器,以获得所需的入射波特征(通常保证入射波上升沿所对应的加载时间大于加载脉冲在试样中3个往复所需时间)并实现恒应变率加载。在本SHTB实验中,脉冲整形器材料采用铜版纸,并加工成外径为28.4 mm、内径为20 mm的环形片,涂抹少量凡士林将其黏附在法兰盘的端面上,脉冲整形器的安装见图8。
试样中应力波波速约为70 m/s,达到基本均匀需要约240 μs才能使波在试样中传播3个来回。经过多次实验发现,实验中加入整形器后所得入射波的上升沿为270~300 μs,试样达到应力平衡状态的时间充足,其波形如图9所示。
1.3.3 入射波完全卸载
SHTB实验过程中,常将吸收杆与入射杆间留出一定空隙(2~5 mm),以避免入射杆与吸收杆之间发生二次撞击,但在PVC弹性体的SHTB实验过程中发现,入射杆与吸收杆间的空隙会导致入射波偏离基线,即未完全卸载,见图10(a)。对实验装置进行分析发现,法兰盘厚度远小于子弹长度,使得子弹与法兰盘一起运动,这是引发入射波偏离基线的主要原因。为了解决这一问题,在实验时应消除入射杆与吸收杆之间的空隙,以卸载入射波,进而获得理想实验波形,见图10(b)。
2. 实验结果与分析
2.1 准静态拉伸应力-应变曲线分析
PVC弹性体在准静态加载下的拉伸载荷-位移曲线如图11所示。由图11可以看出,PVC弹性体在准静态下的拉伸力学性能基本呈线性增加趋势。从载荷-位移曲线中可以计算得出,PVC弹性体在0.1 s−1应变率下的拉伸模量为27 MPa。
2.2 动态拉伸应力-应变曲线分析
对PVC弹性体动态拉伸实验数据进行处理,得到应变率为400~1850 s−1时的拉伸应力-应变曲线,如图12所示。由图12可以看出,PVC弹性体在高应变率拉伸加载条件下的应力-应变曲线可分为弹性、塑性以及卸载3个阶段,呈现出明显的非线性特征,且在不同的高应变率载荷作用下具有不同的拉伸弹性模量、峰值应力、峰值应变,表现出一定的应变率敏感性。图12中,应变为0.1时的对应点与原点连线,用其斜率表征杨氏模量Eg
,称之为割线模量。为进一步说明PVC弹性体的应变率效应,取图12中应变为0.1处对应的割线模量、峰值应力、峰值应变进行比较,其具体变化情况见表1。 表 1 动态参数变化情况Table 1. Variations of dynamic propertiesStrain rate/s−1 Secant modulus/MPa Peak stress/MPa Peak strain/% 400 89.78 9.19 14.12 700 115.62 13.25 16.07 950 134.82 17.99 26.85 1 850 166.81 27.63 53.19 由表1可知,当应变率范围为400~1850 s−1时,割线模量、峰值应力及峰值应变均随应变率的增加而增大。
3. ZWT非线性黏弹性模型
高分子聚合物的动态力学性能可由一个非线性弹簧和两个Maxwell单元组成的朱-王-唐(ZWT)本构模型(见图13)描述,其表达式为
σ=E0ε+αε2+βε3+E1∫t0˙εexp(−t−τθ1)dτ+E2∫t0˙εexp(−t−τθ2)dτ (13) 式中:前3项用于表征材料非线性弹性响应且与应变率无关,后两项分别代表低、高应变率下的黏弹性响应,E0、
α 和β 为试验确定的弹性常数,E1、θ 1和E2、θ 2分别为试验确定的低、高应变率下的弹性常数和松弛时间。当材料受到低应变率加载时,高应变率所对应的Maxwell单元始终处于松弛状态;而当材料受到高应变率加载时,低应变率所对应的Maxwell单元则来不及松弛。相应的模型表达式可简化为σ=E0ε+αε2+βε3+E1∫t0˙εexp(−t−τθ1)dτ (14) σ=(E0+E1)ε+αε2+βε3+E2∫t0˙εexp(−t−τθ2)dτ (15) 本实验的加载率可近似看作恒应变率加载,故式(14)和式(15)可写为
σ=E0ε+αε2+βε3+E1θ1˙ε[1−exp(−εθ1˙ε)] (16) σ=(E0+E1)ε+αε2+βε3+E2θ2˙ε[1−exp(−εθ2˙ε)] (17) 采用式(16)和式(17)拟合PVC弹性体在高应变率下的实验数据,相应的拟合参数见表2。
表 2 ZWT模型拟合参数值Table 2. ZWT model fitting parameter valuesStrain rates/s−1 E0 E1 α β E2 θ2 400 −7.364 6.930 0.166 0.009 8.414 0.002 0 700 −0.350 0.630 0.076 0.004 5.074 0.002 0 950 51.249 −51.151 0.051 0.002 5.482 0.001 7 1 850 2.392 −2.550 0.027 0 8.647 0.001 0 图14为实验结果与本构模型拟合结果的对比。可以看出,两者吻合较好,说明该本构能够较好地描述PVC弹性体在高应变率范围内的动态力学性能。
4. 结 论
开展了PVC弹性体材料的动态拉伸实验,联合波形分析和高速摄像等方法对动态拉伸试样连接方式及胶黏剂进行了优选,实现了恒应变率加载与入射波完全卸载,研究了PVC弹性体在高应变率(400~1850 s−1)拉伸载荷作用下的力学性能,并构建了黏弹性本构模型。通过实验结果分析发现:
(1)提出的基于SHTB装置的动态拉伸技术可用于软质高分子聚合物材料,能满足霍普金森实验的两个基本假设,测试所得结果可反映材料的动态力学性能,且装置结构简单,便于操作和推广;
(2)聚氯乙烯(PVC)弹性体在低应变率(0.1 s−1)拉伸载荷作用下具有明显的非线性弹性特征,PVC弹性体在高应变率(400~1850 s−1)载荷作用下具有不同的拉伸弹性模量、峰值应力、峰值应变,具有明显的应变率敏感性;
(3)PVC弹性体的动态拉伸力学性能的黏弹性特征明显,用ZWT本构模型描述特征时误差较小。
研究结果可为软质高分子聚合物材料的静、动态力学性能研究提供参考。为研究温度对PVC弹性体静、动态拉伸力学性能的影响,后期将用高/低温环境箱开展不同环境温度下的PVC弹性体材料的静、动态拉伸力学性能测试实验。
-
表 1 帽型试件的冲击剪切实验
Table 1. Impact shear experiment on cap shaped specimens
Case L/mm v0/(m·s−1) σim/MPa Δt/μs ts/μs As/mm2 I/(MPa·μs·mm−2) 1 100 43.00 900 39 34 54.6 560.4 2 200 19.00 370 77 54.6 528.6 3 200 24.33 506 77 61 54.6 565.3 4 300 19.12 384 115 81 54.6 569.7 5 300 20.12 412 115 75 54.6 565.9 表 2 采用3种网格尺寸模拟得到的剪切区
Table 2. Results of shear zone simulated with three grid sizes
Mesh size/(μm×μm) Calculation time/h Data file/GB CPU number ASB formation 8×8 90.00 180 8 Yes 15×15 1.36 54 8 Yes 30×30 0.24 30 8 Yes -
[1] KURIYAMA S, MEYERS M A. Numerical modeling of the propagation of an adiabatic shear band [J]. Metallurgical Transactions A, 1986, 17(3): 443–450. doi: 10.1007/BF02643951 [2] BONNET-LEBOUVIER A S, MOLINARI A, LIPINSKI P, et al. Analysis of the dynamic propagation of adiabatic shear bands [J]. International Journal of Solids and Structures, 2002, 39(16): 4249–4269. doi: 10.1016/S0020-7683(02)00244-5 [3] MEDYANIK S N, LIU W K, LI S F. On criteria for dynamic adiabatic shear band propagation [J]. Journal of the Mechanics and Physics of Solids, 2007, 55(7): 1439–1461. doi: 10.1016/j.jmps.2006.12.006 [4] TENG X, WIERZBICKI T, COUQUE H. On the transition from adiabatic shear banding to fracture [J]. Mechanics of Materials, 2007, 39(2): 107–125. doi: 10.1016/j.mechmat.2006.03.001 [5] 付应乾, 董新龙. 帽型试样动态绝热剪切破坏演化分析 [J]. 固体力学学报, 2015, 36(5): 392–400. doi: 10.19636/j.cnki.cjsm42-1250/o3.2015.05.004FU Y Q, DONG X L. Study of evolution of adiabatic shear failure in hat-shaped specimen under dynamic loading [J]. Chinese Journal of Solid Mechanics, 2015, 36(5): 392–400. doi: 10.19636/j.cnki.cjsm42-1250/o3.2015.05.004 [6] GUO Y Z, RUAN Q C, ZHU S X, et al. Dynamic failure of titanium: temperature rise and adiabatic shear band formation [J]. Journal of the Mechanics and Physics of Solids, 2020, 135: 103811. doi: 10.1016/j.jmps.2019.103811 [7] ZHU S X, GUO Y Z, CHEN H S, et al. Formation of adiabatic shear band within Ti-6Al-4V: effects of stress state [J]. Mechanics of Materials, 2019, 137: 103102. doi: 10.1016/j.mechmat.2019.103102 [8] 王财源, 余丰. 复杂应力状态下GCr15轴承钢的动态剪切及熔融特性研究 [J]. 力学季刊, 2022, 43(1): 84–92. doi: 10.15959/j.cnki.0254-0053.2022.01.009WANG C Y, YU F. Research on dynamic shearing and melting characteristics of GCr15 bearing steel under complex stress state [J]. Chinese Quarterly of Mechanics, 2022, 43(1): 84–92. doi: 10.15959/j.cnki.0254-0053.2022.01.009 [9] MEYER L W, MANWARING S. Critical adiabatic shear strength of low alloyed steel under compressive loading [M]//MURR L E, STAUDHAMMER K P, MEYERS M A. Metallurgical Applications of Shock Wave and High-Strain-Rate Phenomena. New York: Marcel Dekker Inc., 1985: 657–674. [10] 王礼立. 应力波基础 [M]. 2版. 北京: 国防工业出版社, 2005.WANG L L. Foundation of stress waves [M]. 2nd ed. Beijing: National Defense Industry Press, 2005. [11] 杨扬, 程信林. 绝热剪切的研究现状及发展趋势 [J]. 中国有色金属学报, 2002, 12(3): 401–408. doi: 10.3321/j.issn:1004-0609.2002.03.001YANG Y, CHENG X L. Current status and trends in researches on adiabatic shearing [J]. The Chinese Journal of Nonferrous Metals, 2002, 12(3): 401–408. doi: 10.3321/j.issn:1004-0609.2002.03.001 [12] TRESCA M H. On further applications of the flow of solids [J]. Proceedings of the Institution of Mechanical Engineers, 1878, 29(1): 301–345. doi: 10.1243/PIME_PROC_1878_029_017_02 [13] RULE W K, JONES S E. A revised form for the Johnson-Cook strength model [J]. International Journal of Impact Engineering, 1998, 21(8): 609–624. doi: 10.1016/S0734-743X(97)00081-X [14] BØRVIK T, DEY S, CLAUSEN A H. Perforation resistance of five different high-strength steel plates subjected to small-arms projectiles [J]. International Journal of Impact Engineering, 2009, 36(7): 948–964. doi: 10.1016/j.ijimpeng.2008.12.003 [15] 周刚毅. TA2钛合金绝热剪切破坏特性及应力状态、晶粒度影响 [D]. 宁波: 宁波大学, 2018.ZHOU G Y. Adiabatic shearing behavior of TA2 titanium alloy and its influence of stress state, grain size [D]. Ningbo: Ningbo University, 2018. [16] LI J R, YU J L, WEI Z G. Influence of specimen geometry on adiabatic shear instability of tungsten heavy alloys [J]. International Journal of Impact Engineering, 2003, 28(3): 303–314. doi: 10.1016/S0734-743X(02)00022-2 [17] WEI Z G, YU J L, LI J, et al. Influence of stress condition on adiabatic shear localization of tungsten heavy alloys [J]. International Journal of Impact Engineering, 2001, 26(1): 843–852. [18] RANC N, TARAVELLA L, PINA V, et al. Temperature field measurement in titanium alloy during high strain rate loading-adiabatic shear bands phenomenon [J]. Mechanics of Materials, 2008, 40(4/5): 255–270. doi: 10.1016/j.mechmat.2007.08.002 [19] WRIGHT T W, OCKENDON H. A model for fully formed shear bands [J]. Journal of the Mechanics and Physics of Solids, 1992, 40(6): 1217–1226. doi: 10.1016/0022-5096(92)90013-R [20] BRONKHORST C A, CERRETA E K, XUE Q, et al. An experimental and numerical study of the localization behavior of tantalum and stainless steel [J]. International Journal of Plasticity, 2006, 22(7): 1304–1335. doi: 10.1016/j.ijplas.2005.10.002 -