Influence of Inertia Effect on Adiabatic Shear Instability
-
摘要:
基于绝热剪切失稳和流体边界层的相似性,运用广义边界层理论方法,建立了绝热剪切失稳的控制方程,并结合量纲分析的方法对控制方程进行了无量纲化,得到了可以表征惯性效应的广义雷诺数,同时运用线性摄动的方法,对无量纲化后的控制方程组进行了求解,分析了惯性效应对绝热剪切失稳的影响.分析果构表明:当惯性效应达到一定临界值以后,绝热剪切失稳会被抑制,且这种抑制作用会随着惯性效应的增强而增强.
Abstract:Based on the similarity between adiabatic shear instability and fluid boundary layer, the governing equations of adiabatic shear instability are established by the generalized boundary layer theory. The governing equations are dimensionless by the method of dimensional analysis, and the generalized Reynolds number which can represent the inertial effect is obtained. At the same time, the dimensionless governing equations are solved by the method of linear perturbation. The influence of inertia effect on adiabatic shear instability is analyzed. The results showed that the adiabatic shear instability would be restrained when the inertia effect reaches a certain critical value, and the restraining effect would increase with the increase of inertia effect.
-
开放科学(资源服务)标志码(OSID):
0. 引言
ZENER和HOLLOMON[1]最早研究了金属材料绝热剪切失稳问题,首次通过试验分析得到, 当应力应变曲线的斜率为负值时,塑性变形已经从均匀变形转变成局部集中的塑性变形.且这种负的斜率是由两种机制相互竞争决定的, 即加工硬化效应和热软化效应.当加工硬化的作用大于热软化的作用时,此时塑性变形是稳定的,即变形是均匀的.反之,则塑性变形会失去稳定性,就会形成绝热剪切失稳.随后,BAI [2]运用线性摄动近似求解的方法,建立简单剪切加载条件下(一维剪切变形)的绝热剪切失稳的判据. BAI的判据从理论上量化了加工硬化效应、应变率硬化效应和热软化效应在绝热剪切失稳中的作用.当硬化的作用小于软化的作用的时候,系统就会发生塑性失稳.MA和LIU [3-4]研究了在复杂应力加载条件下发生塑性失稳的条件,这些塑性失稳包括颈缩塑性失稳、绝热剪切塑性失稳,同时结合塑性失稳相图来判断加工硬化效应、应变率效应和热软化效应在塑性失稳中的作用,由此建立了二维复杂变形条件下的统一塑性失稳判据.MOLINARI[5]通过数值模拟的方法,研究了加工硬化效应、应变率硬化效应、热软化效应、惯性效应对金属材料绝热剪切失稳的影响.但是,到目前为止,还没有文献报道惯性效应在塑性失稳中作用的理论研究.
本文基于广义边界层[6]理论,运用Navier-Stokes方程、连续性方程、应变演化方程、变形协调方程和能量方程建立可以描述加工硬化效应、应变率硬化效应、热软化效应、惯性效应对金属材料绝热剪切失稳影响的理论模型.并运用线性摄动理论对此理论模型进行了近似求解,结合塑性失稳相图研究了惯性效应对绝热剪切失稳的影响.
1. 理论模型的建立
金属材料在高应变率的加载条件下,通常会形成一条窄的剧烈塑性变化区域(见图 1).这种剪切带与流体力学中的边界层相似,是一个很薄的、具有大的变形梯度的区域,其厚度的量级通常在10~100 μm[7],同时伴有局部温度的剧烈变化 [8-9],并且带的传播速度很快,有时能超过1000 ms-1 [10-11].这种问题的建模通常需要考虑塑形大变形、率敏感性、加工硬化、热传导、热软化和惯性效应等因素.纯粹通过固体力学中定理方程已经很难将这些影响因素考虑进去,必须适当借助流体力学边界层的相关理论进行研究.
N-S方程具有如式(1)的数学表达形式:
{ρ(∂v1∂t+v1∂v1∂x1+v2∂v1∂x2)=∂τ11∂x1+∂τ12∂x2+∂p∂x1,ρ(∂v2∂t+v1∂v2∂x1+v2∂v2∂x2)=∂τ21∂x1+∂τ22∂x2+∂p∂x2, (1) 式中: ρ表示材料的密度,v1、v2表示欧拉速度场,p为静水压力(正压力),τij(i, j=1, 2) 表示应力偏量;t为时间变量; x1代表x方向;x2代表y方向.在高速动态加载条件下,假定材料的弹性变形相对塑形变形可以忽略不计,且认为材料密度不变(不可压缩性),则二维条件下的连续性方程为:
∂v1∂x1+∂v2∂x2=0. (2) 同时,材料的变形必须满足变形协调条件:
∂2ε11∂x22+∂2ε22∂x21=2∂2ε12∂x1∂x2, (3) 其中: εij(i, j=1, 2) 为应变场应变值的大小,由于是在欧拉场研究材料的稳定性,则应变随时间的演化方程为:
˙εij=∂εij∂t+v1∂εij∂x1+v2∂εij∂x2,(i,j=1,2), (4) 温度的分布则由能量方程控制:
ρCP(∂T∂t+v1∂T∂x1+v2∂T∂x2)= K(∂2T∂x21+∂2T∂x22)+ β(σ11˙ε11+σ22˙ε22+2σ12˙ε12), (5) 其中:CP为比定压热容;K为热传导系数;β为泰勒-奎尼系数,用来表征功热转化的百分率,其数值大小一般取0.90~0.95.
2. 控制方程的无量纲化
采用广义边界层理论研究材料的稳定性,运用流体力学中的无量纲化方法:
ˆx1=x1L0,ˆx2=x2εL0,^vi=viU0,ˆt=U0tL0,ˆσij=L0˙ε0σijσ0U0,ˆT=TTm, (6) 式中:L0为材料在流动方向的特征长度, ε为小量,通过εL0表征垂直于流动方向的特征长度;U0为外力的加载速度;σ0为材料的初始屈服应力;Tm为材料的熔点.完成无量纲化后,可以得到广义雷诺数Re,表征惯性效应在绝热剪切失稳中的作用.
3. 控制方程的摄动分析
在进行线性摄动分析的时候,采用式(7)中所示的通用描述形式:
{ˆφij}={ˆφij0}+{ˆδˆφij}eˆωˆt+i(ˆk1ˆx1+ˆk2ˆx2),(i,j=1,2), (7) 其中:{ˆφij0}={ˆvi0,ˆσij0,ˆεij0,ˆ˙εij0,ˆT0},且满足{ˆφij0}≫{δˆφij};ˆvi0、ˆσij0、ˆεij0、ˆεij0、ˆT0分别为控制方程组的稳定解;而δˆvi0、δˆσij0、δˆεij0、δˆ˙εij0、δˆT0为控制方程稳定解的微量变化.如果扰动增长率ˆω具有正值,由式(7)可知δˆvi0、δˆσij0、δˆεij0、δˆ˙εij0、δˆT0等微量的变化不是衰减的,而是发散的,此时系统就会因为这种发散而失稳;如果扰动增长率ˆω具有负值,此时各微量是衰减的,系统则处于稳定的状态.基于此原理,通过线性摄动分析来求解扰动增长率ˆω的近似值,从而得到整个系统的稳定性判据.将式(7)代入无量纲化后的控制方程可得到:
{−C10ˆP0ˆk20δˆT+(C10C11ˆQ0ˆk20+ C10C11ˆR0ˆk20ˆω)δˆεN+(2C10ˆQ0ˆk20+ 2C10ˆR0ˆk20ˆω+2Re2ˆω2)δˆε12+ C12δˆvs=0−(C13Λˆk20+C14ΦˆP0+ˆω)δˆT+ (Φ(η11ˆσ11(0)+η22ˆσ22(0))ˆω+ C11C14ΦˆR0ˆω+C11C14ΦˆQ0)δˆεN+ (2ΦReˆσ12(0)ˆω+2C14ΦˆR0ˆω+ 2C14ΦˆQ0)δˆε12+C15δˆvs=0,C16δˆεN−C17δˆε12=0,C18δˆvs=0, (8) 式中:Q为应变硬化系数; η为各应变分量在总应变中所占的比重,由数值模拟结果给出;Φ为无量纲的热传导系数;Λ为无量纲的功热转换系数.要使齐次方程组(8)有解,则必须保证方程组(8)的四阶系数行列式的值为零,由此可得:
2C23Re2ˆω3+(2C16C23Re2ΦˆP0+2C17C23Re2Λˆk20+(C14C23+C12C24)ˆR0ˆk20)ˆω2+((C14C17C23+C12C17C24)ΛˆR0ˆk40+(C14C16C23+C12C16C24−C10C18C24−C10C20C23)ΦˆP0ˆR0ˆk20−C10C24(η11ˆσ11+η22ˆσ22)ΦˆP0ˆk20−2C10C23Φˆσ12ˆP0ˆk20+(C13C23+C11C24)ˆQ0ˆk20)ˆω+(C13C17C23+C11C17C24)ΛˆQ0ˆk40+(C13C16C23+C11C16C24−C10C19C24−C10C21C23)ΦˆP0ˆQ0ˆk20=0. (9) 显然式(9)是一个扰动增长率的三次方程,想要直接求解是比较困难的,采用以波数ˆk20作为未知量来重组方程,最后达到求解系统稳定性的判据.求解过程如下:
aˆk40+bˆk20+c=0, (10) 由式(10)可知,重组后的方程组为波数ˆk20作为未知量的一元二次方程组,利用一元二次方程组的求根公式可得:
ˆk20=−b±√b2−4ac2a, (11) 对方程(10)两边同时对ˆk20求导,即∂ˆω∂ˆk20=0可得:
ˆk20=−b2a. (12) 比较式(11)与式(12)可得:b2-4ac=0.令G1(ˆω)=b2和G2(ˆω)=4ac可得(对于绝热剪切,热传导效应可以忽略不计),
G1(ˆω)=((2C17C23Re2+(C14C23+C12C24)ˆR0/Λ)ˆω+((C14C16C23+C12C16C24−C10C18C24−C10C20C23)ΦˆP0ˆR0(ρU2/Q0)−C10C24(η11ˆσ11+η22ˆσ22)ΦˆP0(ρU2/Q0)−2C10C23Φˆσ12ˆP0(ρU2/Q0)+(C13C23+C11C24)ˆQ0(ρU2/Q0)))2, (13) G2(ˆω)=8C23Re2ˆω((C14C17C23+C12C17C24)(ˆR0/Λ)ˆω+(C13C17C23+C11C17C24)ˆQ0(ρU2/Q0)). (14) 要得到扰动增长率ˆω的正根,必须满足G1(0)>G2(0),可得:
ΦˆP0(C10C24(η11ˆσ11+η22ˆσ22)+2C10C23ˆσ12−(C14C16C23+C12C16C24−C10C18C24−C10C20C23)ˆR0)>(C13C23+C11C24)ˆQ0. (15) 4. 绝热剪切失稳相图分析
式(15)已经得到了绝热剪切失稳的理论判据,本节通过数值模拟和绝热剪切失稳相图来说明惯性效应对绝热剪切失稳的影响.金属切削过程中发生的塑性失稳作为一种典型的绝热剪切失稳模式,本节将通过金属的切削失稳,结合第3节得到失稳判据来说明惯性效应对绝热剪切失稳的影响.数值模拟中采用Johnson-Cook本构关系模型, 数值模拟中采用的材料参数可以参见文献[5].在金属的切削速度分别为50、70、100、150、170、190、210、250、290 m/s时,临界控制函数(13)和(14)在不同切削速度下,扰动增长率与切削速度的变化关系曲线如图 2所示.
由图 2可知,在低速切削阶段,扰动增长率随着切削速度的增大而增大;当切削速度增大到一定程度时,扰动增长率随着切削速度的增大而减小.这和由数值模拟所得的切削力的波动与切削速度的变化趋势一致.
5. 结论
在低应变率的加载条件下,惯性效应在绝热剪切失稳中作用很小,扰动增长率的数值大小为零值,材料不会发生绝热剪切失稳;随着应变率的增大,惯性效应的作用不断增强,但是此时热软化还是主控因素,此时材料会发生绝热剪切失稳;当应变率超过一定的数值,惯性效应成为主控因素,扰动增长率的数值开始不断减小,随着惯性效应的继续增强,扰动增长率将趋向于零,此时将不易发生绝热剪切失稳.本文通过理论研究,揭示了金属材料中惯性效应对绝热剪切失稳的影响,但是不同材料的物理性质和本构关系不相同,因此对于具体的金属材料,需要结合实验和本构关系进一步研究.
-
[1] ZENER C, HOLLOMON J H. Effect of strain rate on plastic flow of steel[J]. Journal of Applied Physics, 1944, 15(1): 22-32. doi: 10.1063/1.1707363
[2] BAI Y L. Thermal-plastic instability in simple shear[J]. Journal of the Mechanics and Physics of Solids, 1982, 30(4): 195-207. doi: 10.1016/0022-5096(82)90029-1
[3] MA W, LI X W, DAI L H, et al. Instability criterion of materials in combined stress states and its application to orthogonal cutting process[J]. International Journal of Plasticity, 2012, 30/31(10): 18-40.
[4] LIU J, WANG G. Effect of inertial on the adiabatic shear instability: numerical analysis[J]. IOP Conference Series: Materials Science and Engineering, 2020, 740: 012145. doi: 10.1088/1757-899X/740/1/012145
[5] MOLINARI A, MUSQUAR C, SUTTER G. Adiabatic shear banding in high speed machining of Ti-6Al-4V: Experiments and modeling[J]. International Journal of Plasticity, 2002, 18(4): 443-459. doi: 10.1016/S0749-6419(01)00003-1
[6] 谷骁勇. 金属高速及超高速切削的物质点法数值模拟技术研究[D]. 北京: 北京理工大学, 2016. GU Xiaoyong. Numerical simulation of high speed and ultra-high speed machining of metal using the material point method[D]. Beijing: Beijing Institute of Technology, 2016.
[7] TIMOTHY S P. The structure of adiabatic shear bands in metals: A critical review[J]. Acta Metallurgica, 1987, 35(2): 301-306. doi: 10.1016/0001-6160(87)90238-0
[8] LEECH P W. Observations of adiabatic shear band formation in 7039 aluminium alloy[J]. Metallurgical Transactions A, 1985, 16: 1900-1903. doi: 10.1007/BF02670378
[9] MRCHAND A, DUFFY J. An experimental study of the formation process of adiabatic shear bands in structural steel[J]. Journal of the Mechanics and Physics of Solids, 1988, 36(3): 261-283.
[10] ZHOU M, ROSAKIS A J, RAVICHANDRAN G. Dynamically propagating shear bands in impact-loaded prenotched plates: Ⅰ. Experimental investigations of temperature signatures and propagation speed[J]. Journal of the Mechanics and Physics of Solids, 1996, 44(6): 981-1006. doi: 10.1016/0022-5096(96)00003-8
[11] MASON J J, ROSAKIS A J, RAVICHANDRAN G. Full field measurements of the dynamic deformation field around a growing adiabatic shear band at the tip of a dynamically loaded crack or notch[J]. J Mech Phys Solids, 1994, 42(11): 1679-1697. doi: 10.1016/0022-5096(94)90067-1