Lateral release effect in shock-loaded specimens during soft recovery process
-
摘要: 通过数值模拟, 计算冲击加载下样品经历一维应变加载过程和侧向稀疏过程产生的塑性功, 给出试样内部从冲击加载开始到进入回收桶前全过程的应力随时间变化的历程。结果表明:侧向稀疏过程开始后,样品在径向汇聚波的作用下受循环拉、压载荷作用,拉压循环的振幅在中等冲击压力下达到最大。如果振幅超过了材料的层裂强度,样品中心将发生拉伸破坏不能完整回收。侧向稀疏与一维应变加载产生的塑性功之比随冲击速度的增加而减小。在冲击速度为某临界值时,侧向稀疏产生的塑性功与一维应变加载产生的塑性功相等。在一定的冲击速度下,采用低初始屈服应力的材料可减轻侧向稀疏效应。对理想塑性材料的理论分析表明,侧向稀疏与一维应变加载产生的塑性功之比随冲击速度与屈服强度比值的增大而减小,与数值模拟结果一致。Abstract: Under shock loading a specimen undergoes a uniaxial-strain loading process and a lateral release process, both of which have an influence on the residual structure, while the influence of the latter is often underestimated or even totally neglected. The plastic work generated in these two processes is calculated in this paper, and the stress history from the beginning of the shock loading to the specimen entering the recovery bin is given. It is found that after the lateral release process begins, the specimen experiences cyclic tension and compression load and the amplitude of the cyclic load reaches its maximum under moderate impact pressure. If the amplitude of the cyclic load is larger than the spall strength, the center of the specimen will be destroyed and the specimen cannot be recovered successfully. The ratio of the plastic work produced during the lateral release to that produced during the uniaxial-strain loading decreases as the impact velocity increases. When the impact velocity reaches a certain critical value, the plastic work produced during the lateral release is equal to that produced during the uniaxial-strain loading. At a certain impact velocity, decreasing the initial yield stress of the materials reduces the lateral release effects. Theoretical analysis of the ideally plastic material shows that the ratio of the plastic work produced during the lateral release to that produced during the uniaxial-strain loading decreases as the ratio of the impact velocity to the yield strength increases, which is consistent with the numerical results.
-
Key words:
- solid mechanics /
- lateral release /
- residual strain /
- plastic work /
- shock loading /
- soft recovery
-
冲击加载下的软回收实验是研究材料动态力学性能的重要途径。通过对回收试样的金相分析(如X射线衍射、透射电子显微镜和电子背散射衍射等),可得到材料的微结构信息,如点阵结构、位错密度、晶粒取向等,从而提高对冲击过程的认识。一个有效的软回收装置要求样品主要受到一维应变加载的影响而不是侧向稀疏的影响。为减轻侧向稀疏效应,W.F.Hartman[1]提出在样品周围加上保护环,后来的一些回收实验也采用了这样的结构[2-3],如图 1(a)所示。为降低样品的动能,G.T.Gray Ⅲ等[4]提出在样品右侧增加层裂板,如图 1 (b)所示。
飞片撞击样品产生的压缩脉冲进入层裂板后,遇自由面反射拉伸脉冲发生层裂,释放了冲击产生的能量,保证样品以低速进入回收桶,防止了样品高速进入回收桶可能产生的二次损伤[5-10]。图 1所示的2种结构通过在样品周围增加保护环减轻了侧向稀疏效应,但减轻的效果如何却少见报道。A.L.Stevens等[11]指出,虽然保护环吸收了侧向稀疏波的能量,但样品卸载后侧向应力并不为零,因此当保护环同样品分离后,新的侧向稀疏效应(径向汇聚波)依然存在,从而对样品内的应力状态、残余应变等产生影响,但关于该影响的大小如何却少见报道。A.L.Stevens等[11]同时指出,影响材料残余效应的因素是塑性功而不是冲击压力或者残余应变,因此准确计算样品内部的塑性功有重要意义。
本文中,基于LS-DYNA有限元软件,对图 1(b)装置中样品的冲击响应过程进行数值模拟。给出装置内部应力分布及样品中心处从冲击加载开始到进入回收桶前全过程的应力随时间变化的关系,计算样品内部一维应变加载过程和侧向稀疏过程产生的塑性功。
1. 计算模型
采用图 1(b)结构进行计算,其中试样厚度为3 mm, 半径为11 mm, 保护环内径为11 mm, 外径为20 mm, 层裂板厚度为5 mm。单元类型为轴对称Lagrange单元,接触条件采用面面接触。飞片、试样、保护环和层裂板材料同为6061-T6 Al, 采用Grüneisen状态方程和各向同性硬化本构模型。材料参数为:密度ρ0=2.7 g/cm3, 剪切模量G=27.6 GPa, 屈服强度Y0=295 MPa, 硬化模量Ep=1.58 GPa, Grüneisen系数γ=2.1, c0=5.37 km/s, s=1.34[11-12]。
2. 计算结果
2.1 试样内部应力状态和应力时间关系
图 2给出了冲击压力为2.4 GPa,冲击速度为300 m/s时不同时刻样品内部的径向应力状态。从图 2(a)看出,飞片撞击盖板产生的压缩脉冲在t=0.8 μs时进入样品,此时侧向稀疏波在保护环中传播不会对样品造成影响。经过0.8 μs后,压缩脉冲完全进入层裂板,如图 2(b)所示。随后,该脉冲将遇自由面反射拉伸脉冲,使层裂板发生层裂,如图 3所示。此外,t=0.8 μs压缩脉冲波尾和侧向稀疏波相遇后在交界处形成负压区(或拉伸区) [13],该负压区随着时间的推移而逐渐扩大,在t=1.6 μs进入保护环,为后续保护环和样品的分离提供驱动力。
从图 3中可以看出,发生层裂后层裂板分裂成2个层裂片,冲击产生的能量一部分用于形成层裂面(转化成表面能),另一部分转化成了2个层裂片的动能。计算中采用最大拉伸应力断裂准则,材料层裂强度为1.2 GPa,比热容为875 J/(kg·K)。
图 4给出了样品和2个层裂片的速度时间历史。从图 4中可以看出,样品在一维应变加载脉冲过后速度趋近于零,冲击产生的能量被2个层裂片带走。另外,样品在冲击加载后处于高温状态,因此回收时需要对样品进行快速降温才能将微结构固定下来,常用的降温材料(如油、水、液氮等)波阻抗都较高,如果样品以100 m/s或km/s量级的高速撞击这些材料势必会造成二次损伤[6],导致微结构的进一步改变,这样对回收样品进行金相分析得到的微观结构就不是仅由冲击加载引起,给后续分析造成困难。由此可见,保证样品低速进入回收桶是十分重要的。图 5~6给出了冲击压力为2.4 GPa时6061-T6Al样品在一维应变加载脉冲下的纵向应力应变和纵向应力径向应力曲线。
从图 2、5和6可以看出,当t=0.8 μs压缩脉冲波头进入样品时,样品从状态0沿路径0→1→2加载到状态2;当t=1.6 μs压缩脉冲波尾离开样品时,样品从状态2沿路径2→3→4反向塑性加载到状态4。设飞片速度为v0,由于是对称碰撞,则状态2的纵向速度v(2)z=12v0(z方向是冲击方向,r、θ分别是径向和环向)。对相容关系dσz=±ρ0cpdvz和dvz=±cpdεz从状态1到状态2积分可得[14]:
σ(2)z=−ρ0cp(v(2)z−v(1)z)+σ(1)z (1) ε(2)z=−(v(2)z−v(1)z)/cp+ε(1)z (2) 式中:cp是塑性波速,对理想塑性材料而言,cp=cb=√K/ρ0, cb是体波波速,K是体积模量。对相容关系dvz=±cpdεz从状态3到状态4积分可得状态4的纵向速度大小为:
v(4)z=(Y0ρ0cb)a0(a0−1)2φ0 (3) 式中:a0=√1+43φ0,φ0=GK。由于很多金属的泊松比约为13,则剪切模量和体积模量的比值φ0约为38,因此a0(a0−1)2φ0约为0.367。Y0ρ0cb项对很多金属在0~30之间,从而vz(4)在0~11 m/s内,因此样品的残余速度是很低的。
从图 6可以看出,状态4样品内的纵向应力σz=0但径向应力σr≠0,因此当保护环在拉应力的作用下(图 2(b)中所示拉伸区)同样品分离后,将形成径向汇聚波向样品中心传播,使样品在随后的过程中处在拉压交替的状态。图 7给出了冲击压力为2.4 GPa时6061-T6 Al样品中心处应力随时间变化的关系。
从图 7中可知, 样品从初始状态0在一维应变加载脉冲的作用下到达状态4,该状态下纵向应力σz=0,剪切应力τrz= 0,仅径向应力σr和和环向应力σθ不为零,且|σr|=|σθ|=σ=Y1D,其中σ为有效应力,Y1D为一维应变加载硬化后的屈服强度。在保护环和样品分离前,σr和σθ始终保持为常数,如图 7(b)中所示的状态4处的平台。保护环的内外径之差越大,该平台的持续时间也越长。保护环和样品分离后,样品边界形成自由面,而内部侧向应力σr不为零,从而将形成径向汇聚波向样品中心传播,侧向稀疏过程开始,如图 7(b)中所示。此后样品将处在拉压交替的状态,拉压循环的周期T=2dce,其中d为样品直径,ce=√(E/ρ0)(1−ν2),为平面应力柱面波的弹性波速[15],E和ν分别为弹性模量和泊松比。拉压循环的振幅A与Y1D近似相等。对于理想塑性材料而言,拉压循环的振幅A为常数Y0。
图 8~9给出了冲击压力为2.4、13.0和60.0 GPa时6061-T6 Al样品中心处有效应力σ和等效塑性应变εp随时间的变化关系,其中对于13.0和60.0 GPa的计算采用Steinberg-Guinan模型,模型参数取自文献[16-17]。
从图 8~9可以看出,3种冲击压力下有效应力和等效塑性应变在t=1.6 μs后不再改变,说明此时一维应变加载过程已经结束而侧向稀疏过程还未开始,此时样品处于状态4, 见图 5~7。从图 8可以看出,Y1D在冲击压力为13.0 GPa时达到最大, 为0.36 GPa;在冲击压力为2.4 GPa时为0.32 GPa;而在冲击压力为60.0 GPa时由于温度软化效应使Y1D达到最小,为0.24 GPa。A.Molinari等[18]指出,冲击压力为9.0 GPa时6061-T6 Al的温升仅70 K,因此,在冲击压力较低时(2.0~13.0 GPa)应变硬化占主导,Y1D随冲击压力的升高而升高。从图 9中看出,冲击压力在60.0 GPa时等效塑性应变达到了0.47,远远高于冲击压力为2.4 GPa时的0.022,但冲击压力为60.0 GPa时的Y1D却低于冲击压力为2.4 GPa时的Y1D(0.24 GPa<0.32 GPa),因此在冲击压力较高时(13.0~60.0 GPa),温度软化占主导,Y1D随冲击压力的升高而降低。在后续的侧向稀疏过程中,由于拉压循环的振幅A和Y1D近似相等,因此中等冲击压力将导致大的振幅。值得注意的是,如果振幅超过了材料的层裂强度,在拉压循环过程中将导致样品中心发生拉伸破坏不能完整回收。
2.2 一维应变加载和侧向稀疏产生的塑性功
长久以来,学者们试图将冲击加载下的残余效应归结为一维应变加载过程(包括冲击压力、压缩脉冲持续时间等)的影响[19]。塑性变形的本质是原子的重组(包括位错的运动、孪晶的形成等等),原子脱离平衡位置所需激活能源自塑性功,等效塑性应变的变化可以反映塑性功的变化。图 10~11给出了冲击压力为2.4 GPa时6061-T6 Al样品中心处径向应力、等效塑性应变随时间的变化关系。
从图 11看出,在侧向稀疏开始前,等效塑性应变历经2次跳跃:第1次是从初始状态0经历弹塑性加载到达状态2, 然后保持不变直到状态3(状态2到状态3历经弹性卸载,不会改变等效塑性应变,见图 5~6);第2次跳跃从状态3历经反向塑性加载到达状态4,然后保持不变,直到样品与保护环分离。侧向稀疏过程开始后样品在30 μs内共历经4次拉压循环(Tn表示第n个周期)。等效塑性应变在第1个周期T1内从0.023上升到0.055,在第2个周期T2内从0.055上升到0.06,在第3个周期T3内从0.06上升到0.063,此后便不再增加。等效塑性应变(或塑性功)的变化率d¯εpdt随周期数n的增加而减小。当周期数n>nmax(此处nmax=3)时,等效塑性应变保持为常数,不再增加,直至样品进入回收桶。假设一维应变加载过程(图 10~11中0~3 μs)产生的塑性功为W1Dp,侧向稀疏过程(图 10~11中3~30 μs)产生的塑性功为Wrp[20],则
Wp1D=∫(σz−σr)dεpz (4) Wpr=∫σijdεpij=∫σrdεpr+∫σθdεpθ (5) 由于侧向稀疏过程开始后,样品中心处的纵向应力σz和剪切应力τrz趋近于零,仅σr和σθ不为零,见图 7,因此式(5)中的ij指标求和只有σr和σθ两项。从文献[21]可知,等效塑性应变率可表示为:
d¯εpdt=√23(˙εp2z+˙εp2r+˙εp2θ) 利用εkkp=0(无塑性体应变)和εrp=εθp,则
d¯εp=|dεpz| (6) 采用LS-DYNA结合式(4)~(6),计算不同冲击速度(0.2~1.2 km/s)下,1100-O Al、6061-T6 Al、LY12、7039 Al、无氧高导电性铜(oxygen-free high-conductivity copper, OFHC)和黄铜样品内部一维应变加载过程产生的塑性功和侧向稀疏过程产生的塑性功,相应的冲击压力范围为1~10 GPa,如图 12所示,材料参数取自文献[11, 22]。
从图 12可以看出,6种材料的侧向稀疏与一维应变加载产生的塑性功之比Wrp/W1Dp随冲击速度的增加而减小。对于6061-T6 Al、LY12和7039 Al而言(见图 12(b)~(c)),在冲击速度为某临界值时(对于6061-T6 Al和7039 Al, 该临界冲击速度分别为425和580 m/s),侧向稀疏产生的塑性功与一维应变加载产生的塑性功相等,低于该速度侧向稀疏效应的影响占主导,这种情况下将残余效应归结为一维应变加载的影响是不正确的;在一定的冲击速度下,1100-O Al和OFHC的Wrp/W1Dp比值低于其余4种材料, 如冲击速度为1.2和0.6 km/s时,1100-O Al、OFHC的Wrp/W1Dp比值分别为0.22、0.20,而6061-T6 Al、LY12的Wrp/W1Dp比值分别为0.29和0.48。由于1100-O Al和OFHC的初始屈服应力明显小于其余4种材料,从而采用低初始屈服应力的材料可减轻侧向稀疏效应。
3. 塑性功的理论分析
对于屈服强度为Y0的理想塑性材料,利用式(6)和样品中心处|σr|=|σθ|=σ=Y1D关系,则式(4)~(5)分别化简为:
Wp1D=Y0ˉεp1DWpr=Y0ˉεp1 式中:ε1Dp为一维应变加载过程累积的等效塑性应变, εlp为侧向稀疏过程累积的等效塑性应变。
根据式(6),则ε1Dp=|εzp(2)-εzp(1)|+|εzp(4)-εzp(3)|,其中εzp(i)为状态i(i=1, 2, 3, 4,见图 5~7)的纵向塑性应变。由于状态1的纵向塑性应变为0,且状态2和状态3的纵向塑性应变相等,则ε1Dp=|εzp(2)|+|εzp(4)-εzp(2)|。
由于状态1、2的应力已知,根据Hooke定律可得状态2的弹性应变εze(2),结合式(2)给出的总应变εz(2),利用εzp(2)=εz(2)-εze(2)可得εp(2)z=−13[v0cb−a0(Y0G)],各物理量定义参见前文。状态4的总应变为ε(4)z=−2Y03K [1],由于状态1、4的应力已知,根据Hooke定律可得状态4的弹性应变为εe(4)z=(3K−2G)Y09KG,两者相减则状态4的塑性应变为εp(4)z=−(3K+4G)Y09KG,把εzp(2)、εzp(4)代入ε1Dp=|εzp(2)|+|εzp(4)-εzp(2)|得:
ˉεp1D=2v03cb−a1Y0G (7) 式中:a1=13(2a0+1)+49φ0。
由于侧向稀疏过程中拉压循环的振幅A和Y0相等,从而侧向稀疏过程累积的等效塑性应变εlp正比于Y0ρ0c2e,则
WprWp1D∝a22(Gv0cbY0)−3a1 (8) 式中:a2=3(3+4φ0)4(3+φ0)。
从式(8)看出,WprWp1D随冲击速度v0的增加而减小。在一定的冲击速度v0下,屈服强度Y0越小WprWp1D越小,因此,从定性规律上看,式(8)的理论预测与前面的数值模拟结果一致。
4. 结论
冲击加载下样品经历了一维应变加载过程和侧向稀疏过程,2种过程对回收试样的残余结构都有影响,而侧向稀疏的影响常常被低估或忽略。本文通过数值模拟,计算了这2种过程产生的塑性功,给出了样品内部从冲击加载开始到进入回收桶前全过程的应力随时间变化的历程。对理想塑性材料,还给出了侧向稀疏与一维应变加载产生的塑性功之比的理论解。得到结论如下:
(1) 一维应变加载结束后,样品中心的应力状态为σz=0,τrz=0,σr=σθ=-Y1D,保护环的内外径之差越大该状态持续的时间越长。侧向稀疏过程开始后,样品在径向汇聚波的作用下受循环拉压载荷作用。拉压循环的周期T≈2dce。拉压循环的振幅A与Y1D在数值上近似相等。Y1D在中等冲击压力下(对6061-T6 Al而言,这个压力为13.0 GPa)达到最大值。当冲击压力小于该最大值时应变硬化占主导,Y1D随冲击压力的升高而升高;反之,温度软化占主导,Y1D随冲击压力的升高而降低。
(2) 侧向稀疏过程开始后,等效塑性应变(或塑性功)的变化率随拉压循环周期数n的增加而减小。当周期数n大于某个nmax时,等效塑性应变保持为常数,不再增加。
(3) 侧向稀疏与一维应变加载产生的塑性功之比随冲击速度的增加而减小。在冲击速度为某临界值时,侧向稀疏产生的塑性功与一维应变加载产生的塑性功相等,低于该冲击速度侧向稀疏的影响将占主导。在一定的冲击速度下,采用低初始屈服应力的材料可减轻侧向稀疏效应。
(4) 对理想塑性材料,给出了侧向稀疏与一维应变加载产生的塑性功之比的理论解。结果表明,侧向稀疏与一维应变加载产生的塑性功之比随Gv0cbY0的增大而减小,其趋势与数值模拟结果一致。
-
-
[1] Hartman W F.Determination of unloading behavior of uniaxially strained 6061 T6 Aluminum from residual strain measurements[J].Journal of Applied Physics, 1964, 35(7):2090-2096. doi: 10.1063/1.1702796 [2] Koller D D, Hixson R S, Gray Ⅲ G T, et al.Influence of shock-wave profile shape on dynamically induced damage in high-purity copper[J].Journal of Applied Physics, 2005, 98(10):103518. doi: 10.1063/1.2128493 [3] Escobedo J P, Dennis-Koller D, Cerreta E K, et al.Effects of grain size and boundary structure on the dynamic tensile response of copper[J].Journal of Applied Physics, 2011, 110(3):033513. doi: 10.1063/1.3607294 [4] Gray Ⅲ G T, Follansbee P S, Frantz C E.Effect of residual strain on the substructure development and mechanical response of shock-loaded copper[J].Materials Science and Engineering A, 1989, 111(89):9-16. [5] Follansbee P S, Gray Ⅲ G T.Dynamic deformation of shock prestrained copper[J].Materials Science and Engineering A, 1991, 138(1):23-31. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CC026527003 [6] Blumenthal W R, Gray Ⅲ G T, Claytor T N.Response of aluminium-infiltrated boron carbide cermets to shock wave loading[J].Journal of Materials Science, 1994, 29(17):4567-4576. doi: 10.1007/BF00376280 [7] Bourne N K, Gray Ⅲ G T, Millett J C F.On the shock response of cubic metals[J].Journal of Applied Physics, 2009, 106(9):091301. doi: 10.1063/1.3218758 [8] Bourne N K, Millett J C F, Gray Ⅲ G T.On the shock compression of polycrystalline metals[J].Journal of Materials Science, 2009, 44(13):3319-3343. doi: 10.1007/s10853-009-3394-y [9] Bourne N K, Gray Ⅲ G T.Soft-recovery of shocked polymers and composites[J].Journal of Physics D:Applied Physics, 2005, 38(19):3690-3694. doi: 10.1088/0022-3727/38/19/018 [10] Mogilevsky M A, Newman P E.Mechanisms of deformation under shock loading[J].Physics Reports, 1983, 97(6):357-393. doi: 10.1016/0370-1573(83)90030-3 [11] Stevens A L, Jones O E.Radial stress release phenomena in plate impact experiments:Compression-release[J].Journal of Applied Mechanics, 1972, 39(2):359-366. doi: 10.1115/1.3422683 [12] Bertholf L D, Karnes C H.Axisymmetric elastic-plastic wave propagation in 6061-T6 Aluminum bars of finite length[J].Journal of Applied Mechanics, 1969, 36(3):533-541. doi: 10.1115/1.3564713 [13] 王继海.二维非定常流和激波[M].北京:科学出版社, 1994. [14] 王礼立.应力波基础[M].2版.北京:国防工业出版社, 2005. [15] Srinivasan M G, Ting T C T.Cylindrical elastic-plastic waves due to discontinuous loading at a circular cavity[J].International Journal of Solids and Structures, 1975, 11(9):1057-1077. doi: 10.1016/0020-7683(75)90048-7 [16] Steinberg D J, Cochran S G, Guinan M W.A constitutive model for metals applicable at highstrain rate[J].Journal of Applied Physics, 1980, 51(3):1498-1504. doi: 10.1063/1.327799 [17] 彭建祥.Johnson-Cook本构模型和Steinberg本构模型的比较研究[D].绵阳: 中国工程物理研究院, 2006. [18] Molinari A, Ravichandran G.Fundamental structure of steady plastic shock waves in metals[J].Journal of Applied Physics, 2004, 95(4):1718-1732. doi: 10.1063/1.1640452 [19] Murr L E, Kuhlmann-wilsdorf D.Experimental and theoretical observations on the relationship between dislocation cell size, dislocation density, residual hardness, peak pressure and pulse duration in shock-loaded nickel[J].Acta Metallurgica, 1978, 26(5):847-857. doi: 10.1016/0001-6160(78)90034-2 [20] 王肖钧, 胡秀章, 李永池.硬化材料中弹塑性柱面波的数值方法[J].爆炸与冲击, 1991, 11(2):97-105. http://www.bzycj.cn/article/id/10775Wang Xiaojun, Hu Xiuzhang, Li Yongchi.A computational method of cylindrical elastic-plastic waves in strain hardening materials[J].Explosion and Shock Waves, 1991, 11(2):97-105. http://www.bzycj.cn/article/id/10775 [21] 李永池, 谭福利, 姚磊, 等.含损伤材料的热粘塑性本构关系及其应用[J].爆炸与冲击, 2004, 24(4):289-298. http://www.bzycj.cn/article/id/9958Li Yongchi, Tan Fuli, Yao Lei, et al. Thermo-viscoplastic constitutive relation of damaged materials with application[J].Explosion and Shock Waves, 2004, 24(4):289-298. http://www.bzycj.cn/article/id/9958 [22] Johnson G R, Cook W H.A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures[C]//Proceedings of the 7th International Symposium on Ballistics.Hague, Netherlands, 1983. 期刊类型引用(0)
其他类型引用(1)
-