基于晶体塑性理论研究铝材料高压高应变率下的强度特性

潘 昊,王升涛,吴子辉,胡晓棉

(北京应用物理与计算数学研究所，北京 100094)

金属材料在军事、航空航天、汽车等领域有广泛的应用，它在冲击加载条件下的弹塑性变形行为一直被关注。冲击加载条件下，材料往往受到强烈的载荷作用，处于高压、高应变率状态，动态力学特性与静态存在显著差别[1-2]。随着检测技术的发展，发现金属材料在高速加载下具有复杂的动力学响应特征，如：单晶Al材料在极端加载条件下可能会发生孪晶[3]，多晶Al材料的强度在高应变率104s-1以上的加载下会大幅增加[4]。另外，也发现在动态加载下金属材料的织构也会随着加载方向及载荷水平发生变化[5]。

为了从理论上描述金属材料在动态加载条件下的力学行为，早期研究基于对实验数据的数学归纳建立了多种宏观本构模型，如MTS模型、JC模型、SG模型等[6-8]，但这些模型中基本没有考虑材料的微结构及其演化，没有从材料塑性变形的物理本质上描述材料的动力学行为，很难普适地描述材料在不同加载条件下的动力学行为。为此，越来越多的本构模型尝试从材料塑性变形的微观机理出发，对位错运动、孪晶演化和硬化机制进行描述，并通过体积平均、自洽场平均等平均化方法得到宏观应力和应变、应变率、温度等参数的关系[9-11]。晶体塑性理论就是其中的代表，且已成为研究材料准静态力学行为的重要方法[12]。该方法从描述晶体的滑移系运动过程出发，将材料的宏观塑性变形和微观的位错运动相联系。但该方法对材料在高速变形下弹性变形过程及位错运动描述不够准确，主要原因包括：(1)大多只考虑应力张量与应变张量的线性关系，这种处理方法适用于体积应变比较小、压力小于1 GPa、温度为常温的情况。在冲击加载下，晶体绝热压缩且有熵增，温升较高，并且体应变相对较大，此时线弹性假设就不再适用，需要考虑温度对弹性变形的影响。(2)位错动力学模型较为简单，指数型经验形式的位错动力学方程只适用于准静态或低应变率加载的情况，对于高压、高应变率加载并不适用。(3)硬化形式考虑不充分，以前的模型大多只考虑了应变硬化，未考虑压力的硬化和温度的软化过程。

针对这些不足，已经开展了一定研究，包括：Clayton[13]基于欧氏应变，通过对三阶弹性常数进行变换，实现了更高压力下的晶体弹性变形过程的描述；Austin等[14]和Mayer等[15]分别提出了适用于冲击加载下位错运动和增值过程的演化方程；Lloyd等[16]采用一维差分方法，对单晶和多晶Al材料的冲击波传播过程进行了模拟。这些研究为基于晶体塑性理论描述材料冲击加载下的力学性能变化提供了有益的帮助，但也存在如位错动力学部分描述过于复杂、只有基于一维加载过程的建模、缺乏对更高压力加载下材料强度及织构演化过程的描述等缺点。

本文中，基于经典晶体塑性理论，针对金属材料在高压、高应变率下的塑性变形特点，采用高阶的欧拉弹性方程和更合理的位错动力学模型和硬化模型，建立充分考虑材料微结构演化的热弹-黏塑性理论模型。和传统的晶体塑性模型相比，本文中的模型不仅能更准确地描述材料在高压下的动态力学行为，也更适用于更大的应变率范围。采用该模型也研究了单晶、多晶Al材料在高压、高应变率加载下强度变化特征和织构演化规律。

1 热弹-黏塑性晶体塑性模型

晶体塑性理论已经有了多种研究框架[17-19]，本文中在Kalidindi等[20]的晶体塑性模型基础上，针对高压、高应变率的加载情况，采用高阶欧拉弹性方程，同时考虑热激活和声子拖曳的位错运动模型，并考虑加工硬化，压力硬化和温度软化，建立了单晶体的热弹-黏塑性晶体塑性模型。

晶体的变形梯度张量可以分解为弹性变形和塑性变形两部分[20]：

F=FeFp

(1)

式中：Fe和Fp分别表示弹性和塑性部分的变形梯度张量。塑性变形梯度的演化可以表示为：

(2)

式中：Lp为塑性速度梯度张量。 只考虑位错滑移对塑性变形的贡献，则：

(3)

1.1 弹性变形方程

欧氏应变的定义为：

(4)

基于欧氏应变，同时考虑熵增Δη对应力的影响，第二Piola-Kirchhoff应力张量为[13]：

(5)

(6)

1.2 位错动力学

1.2.1位错运动

位错运动及演化过程是晶体滑移变形的核心过程，位错滑移运动的控制机制包括短程作用和长程作用。所谓短程作用是指克服由位错交互作用所形成的障碍，晶格的热振动过程有助于克服该障碍，此时的位错运动受热激活过程控制。长程作用是指位错、析出相质点和晶界等引起的长程内应力，此时的位错运动主要受声子的拖曳作用[22]。同时考虑热激活机制和声子拖曳机制时，位错滑移的平均速率可以表示为：

(7)

式中：λ是材料内部长程作用的长度，为0.1～10 μm；tTA是热激活过程特征时间，tD是声子拖曳的特征时间。根据Arrhenius方程：

(8)

式中：ve是有效的越过障碍的频率，k是玻尔兹曼常数，G是活化自由能。G可以表示为：

(9)

(10)

(11)

(12)

(13)

式(7)可以改写为：

(14)

1.2.2位错演化

冲击加载下，材料会产生更多的位错，根据位错理论，并非所有的位错都对滑移系的剪切应变有影响，位错可分为可动和非可动位错两部分。各部分的增值及湮灭过程如下[15,24]：

(15)

(16)

(17)

(18)

(19)

(20)

式中：ρI是非可动位错密度，ρD是可动位错密度，η为经验系数，一般取η=0.1，εD=8 eV/b，vI是位错合并时的特征速度，d表示材料的平均晶粒尺寸，ka是位错湮灭因子。

1.2.3硬化

根据Taylor关系[25]，材料的剪切力与位错满足：

(21)

式中：τβ是初始状态下应力值，与晶格的阻碍和材料的夹杂相关，β表征位错纠缠的长程相互作用系数，G为剪切模量。根据SG本构模型，有：

G=G0[1+dpV1/3-h(T-300)]

(22)

以上模型所采用的参数分别为：λ=0.1～10 μm[22],k=1.38×10-23J/K,ve=1010～1011s-1[22],b=2.86×10-7mm[14],g0=0.65[14],ρ=2.7 g/cm3,θ0=427 K,h=6.626×10-34J·s,cb=5 150 m/s,η=0.1,εD=8 eV/b[16,24],vI=5 m/s[24],ρ0=10-7cm-1[24],ka=10[24],d=10 μm[24],β=1.638[14],G0=25.2 GPa,τβ=27.627 GPa[14]。

1.3 多晶本构模型

由晶体塑性模型可以给出各个晶粒的物理性质(应力应变结果)，将各晶粒的应力应变关系扩展至多晶，则需要满足多晶材料取向分布特征的同时，还要满足晶粒与晶粒之间变形协调条件和晶界处的应力平衡条件。本文中所用求解方法包括Taylor模型和晶体塑性有限元模型，两种模型的具体实现过程如下。

1.3.1Taylor模型[26]

令所有的晶粒在加载过程中拥有同样的变形梯度，根据样品初始的织构分布函数，给出每个晶粒的初始取向(φ0,θ0,ω0)及旋转矩阵Q0，通过旋转变化计算各个晶体的初始弹性常数和Grüneisen系数等，再由单晶体的晶体塑性方程，求解每个晶粒在当前变形梯度下的应力及新的取向(φ,θ,ω)，令总的应力为多个晶粒的平均，就可以得到多晶材料的应力应变关系：

(23)

式中：N为需要统计的晶粒数，V为各晶粒的体积。该模型比较适合于简单加载条件，满足变形协调条件，但并不满足应力平衡条件。相关研究[27-28]表明，对于FCC、HCP多晶材料的简单拉伸及压缩变形问题，在不考虑晶粒间的滑移作用下，Taylor模型中所作的假设仍近似成立，获得的应力应变关系和织构演化与其他模型的结果基本一致。

1.3.2晶体塑性有限元模型(crystalplasticityfiniteelementmethod,CPFEM)[20]

将晶体塑性模型作为有限元程序的一个材料模块， 1个晶粒可以由1个或多个单元构成，通过有限元计算，获得多晶体材料在边界载荷的作用下，整体和局部的变形以及与其他介质的相互作用。在计算过程中，可通过改变界面之间的滑移线(面)类型，模拟晶粒之间的相互作用如滑移、摩擦等。为简单起见，本文中假设晶粒与晶粒是束缚滑移，暂不考虑晶粒之间的滑移作用。有限元方法保证了计算单元之间是满足变形协调和应力平衡的，因此该模型中各晶粒之间也满足变形协调和应力平衡条件。CPFEM能够灵活地施加边界条件，较为适用于多晶材料的复杂冲击-卸载问题。

2 结果与讨论

2.1 位错速度与剪应力的关系

采用上述同时考虑热激活效应和声子拖曳效应的模型计算了Al材料的位错滑移速度和剪应力之间的关系，并将结果和实验结果以及只考虑热激活效应、只考虑声子拖曳效应的模型的计算结果进行了对比，如图1所示。

本文模型的计算得到的位错滑移速度在100～400 MPa的压力下都能和实验结果较好吻合，而只考虑热激活效应或声子拖曳效应，所得位错滑移速度将不能很好吻合实验结果。若只考虑热激活效应，在剪应力较低时，计算得到的位错速度与实验结果较为接近；但当剪应力较大时，模型计算得到的位错速度偏高，和实验结果之间的差距随着剪应力的增加而逐渐变大，甚至可以超过剪切波速，与物理规律不符。若只考虑拖曳效应，在剪应力较低时，模型计算得到的位错速度较实验结果明显偏高。

图1 不同模型的Al材料位错滑移速度和剪应力及与实验结果比较Fig.1 Dislocation glide speed-shear stress curves of Al materials by different theoretical models compared with experimental results

2.2 Al材料剪切强度

采用上述模型，模拟了单晶铝材料在<100>、<110>和<111>加载方向上的剪切强度随压力变化的情况，加载应变率为105s-1。计算得到的剪切强度随加载压力的关系如图2所示，图中散点为文献[29]的实验数据，可见本文中计算的不同压力下单晶Al 材料的剪应力与实验结果符合较好。另外，当Al材料发生塑性变形后，不同加载方向的剪切强度都随压力的增加而单调增加。但在相同加载压力下，不同加载方向的剪切强度大小不同，其中<100>加载方向上剪切强度最低，<110>方向上剪切强度最高，<111>方向上的剪切强度介于两者之间。可见，即使对于对称性较高的FCC 晶体，其单晶材料也具有明显的各向异性。

进一步，在获得单晶材料力学特性的基础上，采用Taylor假设，模拟了500个晶粒(随机取向)的多晶Al材料剪切强度随加载压力的变化，如图3所示。加载条件为一维应变加载，D(2,2)=105s-1，即只在y方向受压，而其余方向为固支约束。由图3可见，本文中模拟的剪切强度-压力的关系与实验结果[30-31]符合较好。另外，图3中还给出了Austin等[32]的计算结果，图中点划线和虚线分别表示其模型预测的Al材料剪切强度的上下界。和该模型相比，本文模型采用了高阶欧拉应变描述材料的弹性变形，能更准确地描述高压下的剪切强度。综合来看，本文发展的热弹-黏塑性晶体塑性模型能比较准确地预测Al材料在高压、高应变率加载条件下的力学行为。

图2 单晶Al材料剪切强度随加载压力的变化Fig.2 Effect of loading pressure on shear strength of single-crystal Al materials

图3 多晶Al材料剪切强度随加载压力的变化Fig.3 Effect of loading pressure on shear strength of polycrystalline Al materials

2.3 织构演化规律

采用上述晶体塑性模型，对一维应变冲击加载下Al 材料织构的变化情况进行了研究，加载条件仍为D(2,2)=105s-1。图4 是500个晶粒{101}方向上织构(初始为随机织构)随加载压力的变化情况，可见随着加载压力的增加，织构发生了一定程度的变化，特别是当加载压力大于20 GPa 后，织构的变化更为明显。这表明，即使只有一个方向能够变形，在冲击加载下，Al 材料的晶体同样会发生旋转，出现晶体择优取向的现象。更深层次上，导致晶体旋转的主要原因是剪切变形，在一维应变冲击加载下，有效塑性变形应约占总变形的2/3，当压力较高、相应的剪切变形较大时，晶体的择优取向效应更明显。

图4 多晶Al材料不同加载压力下的织构演化结果Fig.4 Texture evolution of polycrystalline Al materials under different loading pressures

3 结 论

针对典型金属材料在高温、高压加载条件下的力学响应，基于晶体塑性理论建立了热弹-黏塑性晶体塑性模型。该模型同时考虑了热激活效应和声子拖曳效应，并采用了适用于中高应变率加载的位错运动及硬化方程，能更加准确地描述金属材料在不同加载压力下的弹塑性变形过程。采用该模型研究了加载方向对单晶铝材料剪切强度的影响，结果表明单晶铝材料具有明显的各向异性。最后，研究了多晶Al材料在冲击加载条件下织构的演化过程，发现即使在较低的压力下，织构也会发生择优取向，在更高的压力下，织构的择优取向更加明显。本文中所做的研究主要针对以位错滑移为主要变形机制的Al材料，尚未考虑孪晶等其他塑性变形机制，将在后续研究中予以考虑。