膨胀柱壳恒定应变率的本构关系*

郭昭亮任国武张世文汤铁钢刘仓理

(1.中国工程物理研究院流体物理研究所，四川 绵阳 621999;2.中国工程物理研究院，四川 绵阳 621999)

膨胀柱壳恒定应变率的本构关系*

郭昭亮1，任国武1，张世文1，汤铁钢，刘仓理2

(1.中国工程物理研究院流体物理研究所，四川 绵阳 621999;2.中国工程物理研究院，四川 绵阳 621999)

考虑到金属柱壳膨胀过程中随机裂纹萌生对测试结果的可能干扰，设计了预置有中心穿透圆孔的柱壳样品，采用多普勒光纤探针测量系统获得了柱壳外壁更优的径向速度历史。基于膨胀柱壳实验中固有的非恒定应变率现象，研究了获得恒定应变率下本构方程的方法，并采用改进后的本构方程确定方法，获得了20钢恒定应变率下的应力应变关系。

固体力学；本构方程；多普勒光纤探针；膨胀柱壳；恒定应变率

材料与结构在高应变率下的变形与破坏规律研究，对深入认识材料的动态响应机理以及工程防护有着重要的作用。爆炸膨胀环作为材料在高应变下动态力学性能的实验手段，最早由P.C.Johnson 等[1]建立，后逐步发展为获得高应变率拉伸加载下材料性能的重要工具。近年来，T. Hiroe等[2-3]采用线起爆技术研究了内部爆炸载荷作用下的材料动力学行为。在测试材料力学性能的实验中，对试件的几何设计有着严格的要求，以避免不合适的试件设计带来的非材料因素的影响。汤铁钢等[4-7]讨论了爆炸膨胀环一维应力假定，爆炸膨胀环实验数据处理方法以及爆炸膨胀环界面尺寸效应等问题，并指出在爆炸膨胀环实验研究中，适当增加膨胀环的宽度既可以提高加载应变率，又可以增加膨胀环运动的稳定性。在此理念上，本文中发展了用于研究材料高应变率拉伸本构的膨胀柱壳样品设计。R.H.Warnes等[8]曾指出，膨胀环实验中应变率在实验过程中减小，这并不是一个缺点，可以通过多次实验画出给定材料的应力-应变-应变率曲线。遗憾的是在其分析过程中依然采用了数据窗口中心的数据来计算应变率，这实际上依旧是采用平均应变率的概念，并未在应力应变关系中剥离应变率变化带来的影响。尤其是对于应变率敏感材料，其动态力学性能将会随应变率较小的变化而发生较大的改变，在材料的应力应变关系测试中需要将应变率的影响分离，进而获得更清晰的认识。霍普金森压杆(SHPB)实验中，可以通过控制加载实现恒应变率条件下的测试[9]，但膨胀环实验中的非恒定应变率是由实验的基本原理确定的，变化范围太大，可能已经不能继续选择平均应变率这一概念。

为了获得膨胀环、膨胀柱壳的恒定应变率本构关系，本文中首先分析爆炸柱壳实验中非恒定应变率的由来；随后针对金属柱壳膨胀断裂过程中随机裂纹的萌生，设计预置中心圆孔的样品，成功避免随机裂纹萌生对测速的干扰；最后给出获得金属柱壳恒定应变率下本构关系的方法。

1 膨胀柱壳非恒定应变率起源

中心线起爆加载下的膨胀环、膨胀柱壳实验中(实验装置示意图见文献[4])，一般采用试样外壁的径向速度历史反演应力应变关系，膨胀柱壳的运动分析如图1所示。

图1 膨胀柱壳运动分析示意图Fig.1 Schematic of dynamic analysis of expanding cylinder

在柱壳中心附近沿着z轴选取高度为dz的范围，在自由飞行阶段，金属柱壳仅在环向应力σθ作用下做减速运动，假定材料不可压，在靠近柱壳中心附近局域沿轴向的变形可以忽略不计。环向流动应力、应变、应变率分别为(理论分析类似于膨胀环，见文献[5])：

(1)

(2)

(3)

式中：ρ0为材料密度，r0为初始外半径，r为随时间变化的膨胀半径，vr为瞬时速度。柱壳自由飞行阶段，近似为平面应变状态[10]，σz=(σr+σθ)/2，σr=0，等效应力可写为：

(4)

假定体积不可压，应变之间的关系满足：

εz=0

(5)

εr=-εθ

(6)

因此，等效塑性应变可写为：

(7)

典型的膨胀柱壳外壁径向速度曲线如图2(a)所示，由公式(3)得到对应的应变率随时间的变化，如图2(b)所示。应变率在自由膨胀阶段，从1.0×104s-1递减至5.0s-1左右，在建模过程中，不能简单采用应变率平均值描述材料的动态响应。

图2 膨胀柱壳典型的速度、应变率随时间的变化曲线Fig.2 Curves for typical velocity and strain rate of expanding cylinder

应变率随着时间减小，是由膨胀环、柱壳实验基本原理决定的，在实验中无法避免。直接采用速度历史获得的应力应变关系本质上不同于应变率下的应力应变关系。一般地，我们希望获得某一恒定应变率下的应力应变关系，以便建模；然而，由于爆炸膨胀环、柱壳特殊的几何结构，造就了在整个膨胀变形过程中，应变率的非恒定性。为此，我们将通过一组20钢爆炸膨胀柱壳实验，建立获得恒定应变下应力应变关系的方法。

2 实验设计

膨胀柱壳实验装置示意图如图3(a)所示。诸实验中驱动器为20钢，几何尺寸不变，炸药采用泰安粉末均匀填装于装药套筒，装药高度为56 mm，实验参数如表1所示，表中r、R、h、D、m、ρ分别表示试样内径、试样外径、试样高度、装药直径、装药质量、装药密度。

实验中除了装药量不同外，其他状态完全一致，其中实验2、3为重复性实验。

实验测试布局示意图如图3(b)所示，样品为含有中心穿透圆孔(直径为1 mm)的金属柱壳，测速点位于圆孔正对的位置，以避免在测速点邻近区域萌生裂纹，从而影响速度测量(不含预置缺陷的金属柱壳在内部载荷作用下，裂纹将随机萌生，测速点附近可能会受到局域卸载波的干扰，具体见实验结果分析)。另一方面，我们采用高速摄影技术观测预置缺陷局部的变形、裂纹萌生及裂纹扩展等现象，以此估计预置缺陷对测速点的干扰。

3 实验结果与分析

采用激光干涉测速技术(DPS)，获得了20钢柱壳实验1～4的外壁径向膨胀速度历史，如图4所示。20钢柱壳自由膨胀阶段持续约40～50 μs，由于装药量不同，20钢柱壳的径向速度响应不同，其中实验2与实验3为重复性实验，在结构完全一致、装药量完全相同的情况下，由于装药密度微小的差异(0.6%)，速度差异约为3.1%。从速度曲线可以得到，自由膨胀阶段，20钢实验1～4的应变率变化范围分别为6.14×103～0 s-1、8.62×103～8.1×102s-1、8.76×103～8.1×102s-1、1.01×104～2.6×103s-1。在如此大的应变率变化情况下，如果采用自由膨胀数据窗口中心处的应变率作为平均值，是不合适的。

通过回收样品可以看出，如图5所示，实验1中，圆孔经历大变形之后，局部萌生裂纹，裂纹微小的扩展(约0.5 mm)后，裂纹滞止；实验2与实验3变形较大，但均未完全破坏；实验4中样品完全破碎。

图4 20钢柱壳径向速度曲线Fig.4 Velocity along radial direction of 20 steel expanding cylinder

图5 20钢实验1～4回收样品Fig.5 Recovered samples of 20 steel experiments

高速摄影采用前照明，每幅之间间隔2 μs。图6给出了实验3高速摄影图像，可以清晰地看到预置圆孔局域的膨胀变形，裂纹萌生、扩展及止裂等现象。

图6 实验3膨胀断裂高速分幅图像Fig.6 High-speed image frames of expanding and fracture

在装置的动作过程中，炸药在中心线起爆之后约4 μs冲击波传递至试样，由实验3高速摄影图像可以看出，在20 μs时，裂纹萌生，随后裂纹扩展，当扩展至一定长度后裂纹止裂。随着裂纹的扩展，所发射的Mott卸载波将会绕射至速度的测试点，下面我们将估计Mott卸载波对速度测试的影响。

依据公式(1)～(2)可以得到实验1～4的环向应力与环向应变随时间的变化关系，再由公式(4)、(7)可以获得等效塑性应力应变关系，如图7(a)所示。需要注意的是，此处的应力应变关系曲线上不同点的应变率不同，这种数据对于建模而言，使用比较困难。

图7 应力应变曲线及恒定应变率数据点分布Fig.7 Curve of stress vs. strain and distribution of constant strain rate dates

本轮20钢实验1～4的应变率变化范围分别为6.14×103～0 s-1、8.62×103～8.1×102s-1、8.76×103～8.1×102s-1、1.01×104～2.6×103s-1，其公共应变率范围为6.14×103～2.6×103s-1。为了获得恒定应变率下的流动应力应变关系，原则上可以选择4发实验公共应变率范围内的任意数据点进行处理。不失一般性，选择5.5×103、5.0×103、4.5×103、4.0×103、3.5×103、3.0×103s-1等6组应变率点，尽可能均匀覆盖公共应变率范围，如图7(b)所示。可以看出，随着应变率的线性增加，流动应力准线性增加，实验数据展现出明显的应变率效应。

标准的JC本构方程为

(8)

(9)

并采用修正后的JC本构方程拟合6组应变率点，如图8所示。

图8 恒定应变率下的应力应变曲线Fig.8 Stress-strain curves at constant strain rates

最终确定20钢在拉伸加载下应变率为3.0×103～5.5×103s-1时的JC本构参数为:A=0.25 GPa、B=0.54 GPa、n=0.31、c=0.002、d=5×10-5。从图8可以看出，采用应变率修正后的JC本构方程，可以更好的拟合20钢柱壳在拉伸加载下的宽应变率范围下的实验结果。需要说明的是，公式(9)是一个单调递增函数，在应变较大时无法描述应变软化特征。

上述数据处理方法，通过一组不同载荷强度下膨胀柱壳的实验数据，在其公共的应变率范围内选择了6组应变率下的应力应变数据，通过拟合，获得了恒定应变率下的20钢柱壳拉伸应力应变关系，这种处理手段避免了膨胀环、膨胀柱壳实验中应变率变化范围过大的缺陷。另一方面也发现，利用膨胀环、柱壳实验装置研究材料的拉伸应力应变关系，仅通过单次实验是无法准确获得某一恒定应变率下的本构关系。这就要求在相同的装置结构上，获得载荷强度存在一定差异的一组数据，再通过后期的数据处理得到恒定应变率下的本构关系。虽然为了获得恒定应变率下的本构关系需要的实验数量较多，但回报却是非常丰厚的，理论上可以获得对应应变率变化范围内任意恒定应变率下的应力应变关系。

4 结 论

针对膨胀环、柱壳实验中非恒定应变率现象，从理论和实验两方面进行研究，初步得到如下结论：

(1)分析了膨胀柱壳非恒定应变率的来源，指出这一现象是实验基本原理导致，在实验中无法避免；

(2)设计了预置中心圆孔的金属柱壳实验，可以有效避免随机裂纹对速度测试的影响；

(3)采用不同载荷下相同结构的膨胀柱壳实验，建立了获得恒定应变率下本构方程的方法，获得了20钢柱壳在恒定应变率(3.0×103～5.5×103s-1)下的本构方程。

感谢张振涛、金山、陈浩玉在实验开展中的帮助，感谢刘明涛、范诚在本文撰写过程中的建议。

[1] Johnson P C, Stein B A, Davis R S. Measurement of dynamic plastic flow properties under uniform stress[C]∥Symposium on the Dynamic Behavior of Materials. ASTM Special Publication, 1963:195-198.

[2] Hiroe T, Matsuo H, Fujiwara K, et al. Dynamic behavior of materials induced by explosive loadings initiated using wire explosion techniques[J]. Journal of Materials Processing Technology, 1999,85(1/3):56-59.

[3] Hiroe T, Fujiwara K, Hata H, et al. Deformation and fragmentation behavior of exploded metal cylinders and the effects of wall materials, configuration, explosive energy and initiated locations[J]. International Journal of Impact Engineering, 2008,35(12):1578-1586.

[4] 汤铁钢,李庆忠,陈永涛,等.实现材料高应变率拉伸加载的爆炸膨胀环技术[J].爆炸与冲击,2009,29(5):546-549. Tang Tiegang, Li Qingzhong, Chen Yongtao, et al. An improved technique for dynamic tension of metal ring by explosive loading[J]. Explosion and Shock Waves, 2009,29(5):546-549.

[5] 汤铁钢,桂毓林,李庆忠,等.爆炸膨胀环实验数据处理方法讨论[J].爆炸与冲击,2010,30(5):505-510. Tang Tiegang, Gui Yulin, Li Qingzhong, et al. A discussion of data processing techniques for expanding ring tests[J]. Explosion and Shock Waves, 2010,30(5):505-510.

[6] 汤铁钢,李庆忠,刘仓理,等.爆炸膨胀环的截面尺寸效应[J].爆炸与冲击,2010,30(1):39-44. Tang Tiegang, Li Qingzhong, Liu Cangli, et al. Size effects of expanding ring by explosive loading[J]. Explosion and Shock Waves, 2010,30(1):39-44.

[7] 汤铁钢,李庆忠,陈永涛,等.爆炸膨胀环一维应力假定的分析与讨论[J].爆炸与冲击,2010,30(6):577-582. Tang Tiegang, Li Qingzhong, Chen Yongtao, et al. Discussion about one-dimensional stress presume for explosion expanding ring test[J]. Explosion and Shock Waves, 2010,30(6):577-582.

[8] Warnes R H, Karpp R R, Follansbee P S. The freely expanding ring test-A test to determine material strength at high strain rates[J]. Journal De Physiques IV, 1985,46(C5):583-590.

[9] 宋力,胡时胜.SHPB测试中的均匀性问题及恒应变率[J].爆炸与冲击,2005,25(3):207-216. Song Li, Hu Shisheng. Stress uniformity and constant strain rate in SHPB test[J]. Explosion and Shock Waves, 2005,25(3):207-216.

[10] 奥尔连科.爆炸物理学[M].孙承纬,译.北京:科学出版社,2011:822-825.

[11] Grady D E, Kipp M E. The growth of unstable thermoplastic shear with application to steady-wave shock compression in solids[J]. Journal of the Mechanics and Physics of Solids, 1987,35(1):95-119.

[12] Grady D E. Fragmentation of Rings and Shells: The Legacy of N. F. Mott[M]. New York: Springer, 2006.

[13] 刘旭红,黄西成,陈裕泽,等.强动载荷下金属材料塑性变形本构模型评述[J].力学进展,2007,37(3):361-374. Liu Xuhong, Huang Xicheng, Cheng Yuze, et al. A review on constitutive models for plastic deformation of metal materials under dynamic loadings[J]. Advances in Mechanics, 2007,37(3):361-374.

(责任编辑 曾月蓉)

Constitutive equation of expanding cylindrical shell at constant strain rate

Guo Zhaoliang1, Ren Guowu1, Zhang Shiwen1,Tang Tiegang1, Liu Cangli2

(1.InstituteofFluidPhysics,ChinaAcademyofEngineeringPhysics,Mianyang621999,Sichuan,China;2.ChinaAcademyofEngineeringPhysics,Mianyang621999,Sichuan,China)

In this work, taking into account of the interference of random cracks in the expansion process, we designed a cylindrical sample made with a preset hole. Then, using the Doppler pins system (DPS), we obtained some better curves for the cylinder’s radial speed histories. Based on the non-constant strain rate phenomena in the expanding ring and expanding cylinder experiments, we conducted experiments and studied a method for obtaining constitutive equation at constant strain rate. The improved constitutive equation method was used to obtain the modified constitutive equation at constant strain rate for 20 steel.

solid mechanics; constitutive equation; Doppler pins system(DPS); expanding cylindrical shell; constant strain rate

10.11883/1001-1455(2016)05-0583-07

2015-01-12； < class="emphasis_bold">修回日期：2015-05-06

2015-05-06

国家自然科学基金项目(11172279,11102191)

郭昭亮(1984— )，男，博士研究生,助理研究员,glogos@caep.cn。

O347.3 <国标学科代码：13015 class="emphasis_bold"> 国标学科代码：13015 文献标志码：A国标学科代码：13015

A