高级搜索

应变诱发四方相小角度对称倾侧晶界位错反应的晶体相场模拟

,  ,  , 

downloadPDF

上一篇

下一篇

夏文强, 赵彦, 刘振智, 鲁晓刚. 应变诱发四方相小角度对称倾侧晶界位错反应的晶体相场模拟[J]. 物理学报, 2022, 71(9): 096102-1. doi: 10.7498/aps.71.20212278
引用本文: 夏文强, 赵彦, 刘振智, 鲁晓刚. 应变诱发四方相小角度对称倾侧晶界位错反应的晶体相场模拟[J]. 物理学报, 2022, 71(9): 096102-1. doi: 10.7498/aps.71.20212278
shu
Wen-Qiang Xia, Yan Zhao, Zhen-Zhi Liu, Xiao-Gang Lu. Phase field crystal simulation of strain-induced square phase low-angle symmetric tilt grain boundary dislocation reaction[J]. Acta Physica Sinica, 2022, 71(9): 096102-1. doi: 10.7498/aps.71.20212278
Citation: Wen-Qiang Xia, Yan Zhao, Zhen-Zhi Liu, Xiao-Gang Lu. Phase field crystal simulation of strain-induced square phase low-angle symmetric tilt grain boundary dislocation reaction[J]. Acta Physica Sinica, 2022, 71(9): 096102-1. doi: 10.7498/aps.71.20212278
shu

应变诱发四方相小角度对称倾侧晶界位错反应的晶体相场模拟

  • 上海大学材料科学与工程学院, 省部共建高品质特殊钢冶金与制备国家重点实验室, 上海　200444

    • 通讯作者: E-mail: zhaoyan8626@shu.edu.cn; 

  • 中图分类号: 61.72.Bb, 61.72.Hh, 07.05.Tp

Phase field crystal simulation of strain-induced square phase low-angle symmetric tilt grain boundary dislocation reaction

  • State Key Laboratory of Advanced Special Steel, School of Materials Science and Engineering, Shanghai University, Shanghai 200444, China

    • Corresponding author: E-mail: zhaoyan8626@shu.edu.cn; 

  • MSC: 61.72.Bb, 61.72.Hh, 07.05.Tp

  • 摘要: 采用晶体相场法研究了外加应变作用下, 不同取向差的四方相对称倾侧小角度晶界的位错运动与反应及反应过程中的位错组态, 通过采用几何相位法对位错周围应变场进行了表征. 结果表明, 凝固弛豫达到稳态后, 晶界两侧位错平行且方向相反, 随晶界两侧晶粒取向差增大, 位错数目增加, 距离减小, 且体系自由能增加. 在外加应变作用下, 晶界位错经历攀移、发射、反应湮灭等阶段, 体系自由能呈现波动. 当取向差增大时, 位错运动方式由攀移向攀滑移转变, 产生更多类型的位错组构型, 并发生相应的位错与位错组之间的反应. 对于不同构型的位错反应, 正切应变驱动位错靠近, 负切应变驱动位错湮灭.
    Abstract: In this paper, the phase field crystal method is used to study the dislocation motion and reaction of the square phase symmetric tilt low-angle grain boundaries, and the dislocation configurations with different misorientation angles are analyzed under the action of applied strain. The geometric phase approach is used to characterize the strain field around the dislocations. The results show that after the solidification relaxation, the interfacial dislocations on both sides of the grain are distributed in parallel but opposite direction. With the increase of misorientation angle between grains, the number of dislocations increases, the spacing between them decreases, and the free energy of the system increases. Imposed by the applied strain, the grain boundary dislocations undergo climbing, launching, and reactive annihilation, with the free energy fluctuating. When the misorientation increases, the dislocation motion mode changes from climbing to climbing-sliping, resulting in more dislocation group configurations, and more reactions between dislocations and dislocation groups. For the dislocation reactions of different configurations, positive shear strain drives dislocations to approach, and negative shear strain drives dislocations to annihilate.
  • 加载中

  • 图 1  二维四方相相图(黑色部分为两相共存区)

    Figure 1.  Phase diagram of the two-dimensional square lattices (the black part is the two-phase coexistence region).

    下载: 全尺寸图片 幻灯片

    图 2  $ \theta ={4} $°时双晶弛豫过程的模拟　(a) t = 100; (b) t = 1300; (c) t = 1700; (d) t = 30000

    Figure 2.  Simulation of the relaxation process of bicrystals with $ \theta ={4} $°: (a) t = 100; (b) t = 1300; (c) t = 1700; (d) t = 30000.

    下载: 全尺寸图片 幻灯片

    图 3  不同取向差对应的位错构型、位错数目与位错间距

    Figure 3.  Dislocation configuration, number of dislocations and dislocation distance corresponding to different misorientations.

    下载: 全尺寸图片 幻灯片

    图 4  取向差对弛豫过程体系自由能变化的影响

    Figure 4.  Effect of misorientations on the change of free energy during relaxation process.

    下载: 全尺寸图片 幻灯片

    图 5  不同取向差下晶界处$ {\varepsilon }_{xx} $应变云图　(a)$ \theta ={4} $°; (b)$ \theta ={6} $°; (c) $ \theta ={8} $°

    Figure 5.  Contours of $ {\varepsilon }_{xx} $ strain at the grain boundary under different misorientations: (a)$ \theta ={4} $°; (b)$ \theta ={6} $°; (c) $ \theta ={8} $°.

    下载: 全尺寸图片 幻灯片

    图 6  应变作用下$ \theta ={4} $°时晶界位错运动模拟图　(a) t = 50; (b) t = 14150; (c) t = 15250; (d) t = 17550; (e) t = 29100; (f) t = 48500. 图6(a)中插图(g)和(h)分别表示刃位错的类型和位错对应的伯氏矢量类型

    Figure 6.  Simulation of grain boundary dislocation motion under strain with $ \theta ={4} $°: (a) t = 50; (b) t = 14150; (c) t = 15250; (d) t = 17550; (e) t = 29100; (f) t = 48500. Inset (g) and (h) represent the type of edge dislocation and the Burgers vector type respectively corresponding to the dislocation in Fig. 6(a).

    下载: 全尺寸图片 幻灯片

    图 7  应变作用下$ \theta ={6} $°时晶界位错运动模拟图　(a) t = 50; (b) t = 15000; (c) t = 16350; (d) t = 19550; (e) t = 23500; (f) t = 26950

    Figure 7.  Simulation of grain boundary dislocation motion under strain with $ \theta ={6} $°: (a) t = 50; (b) t = 15000; (c) t = 16350; (d) t = 19550; (e) t = 23500; (f) t = 26950.

    下载: 全尺寸图片 幻灯片

    图 8  应变作用下$ \theta ={8} $°时晶界位错运动模拟图　(a) t = 100; (b) t = 16350; (c) t = 17500; (d) t = 20400; (e) t = 25000; (f) t = 34750

    Figure 8.  Simulation of grain boundary dislocation motion under strain with $ \theta ={8} $°: (a) t = 100; (b) t = 16350; (c) t = 17500; (d) t = 20400; (e) t = 25000; (f) t = 34750.

    下载: 全尺寸图片 幻灯片

    图 9  不同取向差下体系自由能变化曲线　(a) $ \theta ={4}$°; (b) $ \theta ={6} $°; (c) $ \theta ={8} $°

    Figure 9.  Variations of free energy under different misorientations: (a) $ \theta ={4} $°; (b) $ \theta ={6} $°; (c) $ \theta ={8} $°.

    下载: 全尺寸图片 幻灯片

    图 10  含位错组的晶体位错构型及对应的不同应变取向应变云图　(a) 晶体构型图; (b) ${\varepsilon }_{xx}$; (c) ${\varepsilon }_{xy}$; (d) ${\varepsilon }_{yy}$

    Figure 10.  Dislocation configuration of crystals containing dislocation groups and corresponding strain contours of different strain orientations: (a) Crystal configuration diagram; (b) ${\varepsilon }_{xx}$; (c) ${\varepsilon }_{xy}$; (d) ${\varepsilon }_{yy}$.

    下载: 全尺寸图片 幻灯片

    图 11  不同位错反应构型的演化图　(a)—(d) R1构型; (e)—(h) R2构型; (i)—(l) R3构型; (m)—(p) R4构型; (q)—(t) R5构型

    Figure 11.  Evolution of different dislocation reaction configurations: (a)–(d) R1 configurations; (e)–(h) R2 configurations; (i)–(l) R3 configurations; (m)–(p) R4 configurations; (q)–(t) R5 configurations.

    下载: 全尺寸图片 幻灯片

    图 12  不同位错反应构型演化时应变(${\varepsilon }_{xy}$)云图　(a)—(d) R1构型; (e)—(h) R2构型; (i)—(l) R3构型; (m)—(p) R4构型; (q)—(t) R5构型

    Figure 12.  Strain (${\varepsilon }_{xy}$) contours during the evolution of different dislocation reaction configurations: (a)–(d) R1 configuration; (e)–(h) R2 configuration; (i)–(l) R3 configuration; (m)–(p) R4 configuration; (q)–(t) R5 configuration.

    下载: 全尺寸图片 幻灯片

    表 1  模拟设定的参数

    Table 1.  Parameters in the simulation.

    方案初始原子密度$ {\overline \psi _0} $温度相关参量ε取向差θ
    A$－0.40$$ 0.45 $$ {4} $°
    B$－0.40$$ 0.45 $$ {6} $°
    C$－0.40$$ 0.45 $$ {8} $°
    下载: 导出CSV

    表 2  不同构型的位错分解发射与位错反应

    Table 2.  Dislocation decomposition and dislocation reaction of different configuration.

    编号位错分解发射编号位错反应
    L1$n1\xrightarrow{\text{分解}}n1+n3+n4$R1$ n3+n4\xrightarrow{\text{反应}} 0 $
    $ n2 \xrightarrow{\text{分解}} n2+n3+n4 $R2$n1+\overline{n2 n4}\xrightarrow{\text{反应} }n4,$
    $n2+\overline{n2 n3}\xrightarrow{\text{反应}}n4$
    L2$n1\xrightarrow{\text{分解}}\overline{n1 n3}+n4$R3$n3+\overline{n2 n4}\xrightarrow{\text{反应} }n2,$
    $n4+\overline{n1 n3}\xrightarrow{\text{反应}}n1$
    $n2\xrightarrow{\text{分解}}\overline{n2 n3}+n4$R4$\overline{n1 n3}+\overline{n2 n4}\xrightarrow{\text{反应} }0,$
    $\overline{n1 n4}+\overline{n2 n3}\xrightarrow{\text{反应}}0$
    L3$n1\xrightarrow{\text{分解}}\overline{n1 n4}+n3$R5$ n1+n2\xrightarrow{\text{反应}}0 $
    $n2\xrightarrow{\text{分解}}\overline{n2 n4}+n3$
    下载: 导出CSV
  • [1] Elder K R, Katakowski M, Haataja M, Grant M 2002 Phys. Rev. Lett. 88 245701 doi: 10.1103/PhysRevLett.88.245701
    [2] Elder K R, Grant M 2004 Phys. Rev. E 70 051605 doi: 10.1103/PhysRevE.70.051605
    [3] Elder K R, Provatas N, Berry J, Stefanovic P, Grant M 2007 Phys. Rev. B 75 064107 doi: 10.1103/PhysRevB.75.064107
    [4] 高英俊, 罗志荣, 黄创高, 卢强华, 林葵 2013 物理学报 62 050507 doi: 10.7498/aps.62.050507 Gao Y J, Luo Z R, Huang C G, Lu Q H, Lin K 2013 Acta Phys. Sin. 62 050507 doi: 10.7498/aps.62.050507
    [5] Yang T, Gao Y, Wang D, Shi R P, Chen Z, Nie J F, Li J, Wang Y 2017 Acta Mater. 127 481 doi: 10.1016/j.actamat.2017.01.036
    [6] Greenwood M, Rottler J, Provatas N 2011 Phys. Rev. E 83 031601 doi: 10.1103/PhysRevE.83.031601
    [7] Gao Y J, Deng Q Q, Liu Z Y, Huang Z J, Li Y X, Luo Z R 2020 J. Mater. Sci. Technol. 49 236 doi: 10.1016/j.jmst.2020.01.030
    [8] Gao Y J, Huang L L, Deng Q Q, Zhou W Q, Luo Z R, Lin K 2016 Acta Mater. 117 238 doi: 10.1016/j.actamat.2016.06.021
    [9] 赵宇龙, 陈铮, 龙建, 杨涛 2013 物理学报 62 118102 doi: 10.7498/aps.62.118102 Zhao Y L, Chen Z, Long J, Yang T 2013 Acta Phys. Sin. 62 118102 doi: 10.7498/aps.62.118102
    [10] 赵宇龙, 陈铮, 龙建, 杨涛 2014 金属学报 27 81 doi: 10.1007/s40195-014-0027-5 Zhao Y L, Chen Z, Long J, Yang T 2014 Acta Metall. Sin. 27 81 doi: 10.1007/s40195-014-0027-5
    [11] Gao Y J, Luo Z R, Huang L L, Mao H, Huang C G, Lin K 2016 Modell. Simul. Mater. Sci. 24 055010 doi: 10.1088/0965-0393/24/5/055010
    [12] Hu S, Chen Z, Xi W, Peng Y Y 2017 J. Mater. Sci. 52 5641 doi: 10.1007/s10853-017-0799-x
    [13] Seymour M, Sanches F, Elder K, Provatas N 2015 Phys. Rev. B 92 184109 doi: 10.1103/PhysRevB.92.184109
    [14] Faghihi N, Provatas N, Elder K R, Grant M, Karttunen M 2013 Phys. Rev. E 88 032407 doi: 10.1103/PhysRevE.88.032407
    [15] Berghoff M, Nestler B 2015 Comput. Condens. Matter 4 46 doi: 10.1016/j.cocom.2015.08.002
    [16] Tang S, Wang Z J, Guo Y L, et al. 2012 Acta Mater. 60 5501 doi: 10.1016/j.actamat.2012.07.012
    [17] 高英俊, 卢成健, 黄礼琳, 罗志荣, 黄创高 2014 金属学报 50 110 doi: 10.3724/SP.J.1037.2013.00308 Gao Y J, Lu C J, Huang L L, Luo Z R, Huang C G 2014 Acta Metall. Sin. 50 110 doi: 10.3724/SP.J.1037.2013.00308
    [18] Hu S, Chen Z, Yu G G, Xi W, Peng Y Y 2016 Comput. Mater. Sci. 124 195 doi: 10.1016/j.commatsci.2016.07.030
    [19] Gutkin M Y, Ovid'ko I A 2001 Phys. Rev. B 63 064515 doi: 10.1103/PhysRevB.63.064515
    [20] Stefanovic P, Haataja M, Provatas N 2009 Phys. Rev. E 80 046107 doi: 10.1103/PhysRevE.80.046107
    [21] Bobylev S V, Gutkin M Y, Ovid'ko I A 2004 Phys. Solid State 46 2053 doi: 10.1134/1.1825548
    [22] Ohno T, Ii S, Shibata N, Matsunaga K, Ikuhara Y, Yamamoto T 2004 Mater. Trans. 45 2117 doi: 10.2320/matertrans.45.2117
    [23] Humphreys F J 2001 J. Mater. Sci. 36 3833 doi: 10.1023/A:1017973432592
    [24] Hirouchi T, Takaki T, Tomita Y 2009 Comput. Mater. Sci. 44 1192 doi: 10.1016/j.commatsci.2008.08.001
    [25] 龙建, 王诏玉, 赵宇龙, 龙清华, 杨涛, 陈铮 2013 物理学报 62 218101 doi: 10.7498/aps.62.218101 Long J, Wang Z Y, Zhao Y L, Long Q H, Yang T, Chen Z 2013 Acta Phys. Sin. 62 218101 doi: 10.7498/aps.62.218101
    [26] 鲁娜, 王永欣, 陈铮 2014 物理学报 63 180508 doi: 10.7498/aps.63.180508 Lu N, Wang Y X, Chen Z 2014 Acta Phys. Sin. 63 180508 doi: 10.7498/aps.63.180508
    [27] 杨涛, 陈铮, 董伟平 2011 金属学报 47 1301 doi: 10.3724/SP.J.1037.2011.00111 Yang T, Chen Z, Dong W P 2011 Acta Metall. Sin. 47 1301 doi: 10.3724/SP.J.1037.2011.00111
    [28] Berkels B, Ratz A, Rumpf M, Voigt A 2008 J. Sci. Comput. 35 1 doi: 10.1007/s10915-007-9157-5
    [29] Elsey M, Wirth B 2014 Multiscale Model. Sim. 12 1 doi: 10.1137/130916515
    [30] Singer H M, Singer I 2006 Phys. Rev. E 74 031103 doi: 10.1103/PhysRevE.74.031103
    [31] Wang Z J, Li J J, Guo Y L, Tang S, Wang J C 2013 Comput. Phys. Commun. 184 2489 doi: 10.1016/j.cpc.2013.06.018
    [32] Hytch M J, Snoeck E, Kilaas R 1998 Ultramicroscopy 74 131 doi: 10.1016/S0304-3991(98)00035-7
    [33] Galindo P L, Kret S, Sanchez A M, et al. 2007 Ultramicroscopy 107 1186 doi: 10.1016/j.ultramic.2007.01.019
    [34] Hytch M J, Putaux J L, Penisson J M 2003 Nature 423 270 doi: 10.1038/nature01638
    [35] Guo Y, Wang J, Wang Z, Li J, Tang S, Liu F, Zhou Y 2015 Philos. Mag. 95 973 doi: 10.1080/14786435.2015.1011250
    [36] Wu K A, Adland A, Karma A 2010 Phys. Rev. E 81 061601 doi: 10.1103/PhysRevE.81.061601
    [37] Chen L Q, Shen J 1998 Comput. Phys. Commun. 108 147 doi: 10.1016/S0010-4655(97)00115-X
  • 加载中
图( 12) 表( 2)
计量
  • 文章访问数:  302
  • HTML全文浏览数:  302
  • PDF下载数:  0
  • 施引文献:  0
出版历程
  • 收稿日期:  2021-12-08
  • 刊出日期:  2022-05-05

应变诱发四方相小角度对称倾侧晶界位错反应的晶体相场模拟

    通讯作者: E-mail: zhaoyan8626@shu.edu.cn; 
  • 上海大学材料科学与工程学院, 省部共建高品质特殊钢冶金与制备国家重点实验室, 上海　200444
  • 收稿日期:  2021-12-08

摘要: 采用晶体相场法研究了外加应变作用下, 不同取向差的四方相对称倾侧小角度晶界的位错运动与反应及反应过程中的位错组态, 通过采用几何相位法对位错周围应变场进行了表征. 结果表明, 凝固弛豫达到稳态后, 晶界两侧位错平行且方向相反, 随晶界两侧晶粒取向差增大, 位错数目增加, 距离减小, 且体系自由能增加. 在外加应变作用下, 晶界位错经历攀移、发射、反应湮灭等阶段, 体系自由能呈现波动. 当取向差增大时, 位错运动方式由攀移向攀滑移转变, 产生更多类型的位错组构型, 并发生相应的位错与位错组之间的反应. 对于不同构型的位错反应, 正切应变驱动位错靠近, 负切应变驱动位错湮灭.

English Abstract

    全文HTML

    1.   引　言

    • 基于密度泛函的晶体相场法[1-3](phase-field-crystal, PFC)能够实现扩散时间尺度与原子空间尺度的模拟, 且自洽耦合了晶体点阵结构, 从而被广泛用于相结构转变[4-6]、晶界变形与位错机制[7,8]、反霍尔-佩奇效应[9,10]、微观裂纹扩展[11,12]、铁电磁复合材料的微结构[13,14]、凝固的枝晶生长[15,16]等的模拟, 是探究原子尺度微观组织演化的重要方法之一.

      材料的加工变形行为与晶界、位错等微观缺陷的构型及运动有关. 在外加应力、应变作用下, 位错将发生运动、分解[17]、反应[18], 导致晶界迁移, 甚至湮灭[19-21], 从而影响材料的力学性能. 近年来有许多关于实验方法研究晶界和位错的报道[22,23], 但难以观察到其运动, 而通过PFC方法能够从理论模拟上研究位错运动. Hirouchi等[24]通过引入等体积法假设, 实现了在外加应变作用下PFC方法模拟晶界与位错的运动, 随后的研究者在此基础上进行了进一步研究[25-27]. 但由于其内在的规律及机理较为复杂, 上述研究主要讨论了外在因素对晶界、位错运动的影响.

      同时, 晶界、取向等特征以及缺陷引起的应变的表征对理解晶界和位错运动有着重要意义. 基于Mumford-Shah型泛函[28]和能量泛函最小化的变分方法[29]能够实现对晶界的表征, 并且上述方法还能表征缺陷周围的应变场. 为了捕获多晶变形中的晶粒取向, 通过引入小波函数, 利用小波变换的方法[30]可视化晶粒取向. Wang等[31]成功地将该方法由二维推向了三维. 除此之外, 几何相位法[32] (geometric phase approach, GPA)和峰对法[33] (peak pairs algorithm, PPA)等也能表征应变. PPA是利用原子近邻峰之间的关系获得局部应变值, 但对于原子位置的精确程度较高, 而GPA通过计算相位场获得应变, 并兼顾两个实虚空间的信息, 不易受原子位置的干扰, 从而常用于原子尺度下应变的表征[34,35].

      本文采用PFC方法模拟不同取向差下对称倾侧小角度晶界的结构及其外加应变下位错运动和相互作用. 通过位错运动与体系自由能之间的联系, 对位错的发射和位错反应湮灭的过程进行研究, 分析不同位错反应的构型特征, 并表征了晶界和位错周围的应变场分布, 揭示位错反应的内在机理.

    2.   计算模拟方法

      2.1.   双模晶体相场模型

    • 采用Wu等[36]提出的双模PFC模型构建四方相点阵结构, 其无量纲的自由能密度函数F可以表示为

      式中, $ \psi \left( {r, t} \right) $作为序参量表示原子密度场, 其原子位置不随时间变化; 而对于液相, 时间发生改变, 原子位置将发生变化, 用常数$ {\overline \psi _0} $表示其平均原子密度; ε为表征温度的参数, 控制固-液两相区的大小, 随着ε值减小, 温度逐渐降低; Q为次邻近波矢幅值与最近邻波矢幅值之比, 该值由晶体结构决定; R则控制双模与单模结构的相对稳定性, 随着R的值增大, 双模会逐渐向单模转变. 对于本文构建的四方相$ Q = \sqrt 2 $, $ R = 1 $.

      对于无量纲原子密度场$ \psi \left( {r, t} \right) $, 在体积不变的情况下是一个保守场变量, 其演化过程可通过Cahn-Hilliard[1]动力学方程来描述:

      将(1)式代入(2)式得

      采用半隐式Fourier谱方法[37]对方程(3)进行数值求解.

      采用双模近似, 四方相的原子密度函数的表达式为[36]

      式中, $ {\overline \psi _0} $为平均原子密度值; $ {A_{{\text{sq}}}} $和$ {B_{{\text{sq}}}} $分别表示密度波的振幅; $ {\boldsymbol{q}} $为波矢, $q = 1.$

      将四方相原子密度函数(4)式代入(1)式, 使分别关于$ {A_{{\text{sq}}}} $和$ {B_{{\text{sq}}}} $自由能最小化得到$ {A_{{\text{sq}}}} $和$ {B_{{\text{sq}}}} $. 经过能量最小化后的四方相自由能密度函数表达式为

      利用公切线法则, 通过化学势((6)式)与巨势((7)式)相等, 求解(6)式与(7)式的方程组, 可得到如图1所示的双模二维PFC相图.

      • 2.2.   模拟参数设置

    • 以四方相为对象, 所设置的模拟参数如表1所列, 平均原子密度为$ {\overline \psi _0} = - 0.4 $, 温度为$ \varepsilon = 0.45 $, 空间步长为$\Delta x = \Delta y = {{\text{π }}}/{4}$, 时间步长为$ \Delta t = 0.5 $, 边界条件采用周期性边界条件. 计算时, 模拟区域的大小为$ {L_x} \times {L_y} = 512\Delta x \times 512\Delta y $, 初始区域分为三部分, $ 0 < y < \dfrac{{{L_y}}}{8} + {d_0} $, $ \dfrac{{3{L_y}}}{8} - {d_0} < y < \dfrac{{5{L_y}}}{8} + {d_0} $, $ \dfrac{{7{L_y}}}{8} - {d_0} < y < {L_y}\left( {0 < x < {L_x}} \right) $. 其中$0 < y < \dfrac{{{L_y}}}{8} $$ + {d_0}$与$\dfrac{{7{L_y}}}{8} - {d_0} < y < {L_y}$属于上下两侧晶粒, 取向角为$ - \dfrac{\theta }{2} $; $ \dfrac{{3{L_y}}}{8} - {d_0} < y < \dfrac{{5{L_y}}}{8} + {d_0} $属于中间晶粒, 取向角为${\theta }/{2}$; $ {d_0} $表示两晶粒之间设置的液相区域宽度为${{{L_y}}}/{8}$. 通过液-固相转变, 经过30000时间步长弛豫过程得到稳定的对称倾侧小角度晶界.

      • 2.3.   单轴拉伸模型

    • 采用等体积法[24]模型, 模拟不同取向四方相的形变过程, 则有

      式中, $ \Delta x $和$ \Delta y $为初始空间步长, $ \Delta x' $和$ \Delta y' $为变形后的空间步长, s为模拟区域的面积. 如对体系在x方向上施加一个拉应变, 经过n个时间步长的变形后, 空间步长为

      式中, $ \dot \varepsilon $为应变速率, 应变量表示为$ \varepsilon = n\dot \varepsilon \Delta t $, 本文中采用的应变速率为$ \dot \varepsilon = 5 \times {10^{ - 6}} $.

      • 2.4.   晶体应变表征方法

    • 为了反映晶界与位错等缺陷的应变, 采用GPA[32]方法对其进行表征. 在GPA方程的构造中, 原子图像的密度场$ I\left({\boldsymbol{r}}\right) $可以表示为

      式中, $ {\boldsymbol{g}} $为无畸变的参考点阵的倒格矢, $ {\boldsymbol{r}} $为实空间的位置矢量. 通过对局部傅里叶分量$ {I_{\boldsymbol{g}}}\left( {\boldsymbol{r}} \right) $的计算, 可获得相对于参考点阵的畸变. 经滤波后, 局部傅里叶分量为

      式中, $ {A_{\boldsymbol{g}}}\left( {\boldsymbol{r}} \right) $为振幅场, $ {P_{\boldsymbol{g}}}\left( {\boldsymbol{r}} \right) $为相位场. 而$ {P_{\boldsymbol{g}}}\left( {\boldsymbol{r}} \right) $与位移场$ {\boldsymbol{\mu}} \left( {\boldsymbol{r}} \right) $关系为

      对于二维图像, 要确定位移场需要两个相位场: $ {P_{{{\boldsymbol{g}}_1}}}\left( {\boldsymbol{r}} \right) $和$ {P_{{{\boldsymbol{g}}_2}}}\left( {\boldsymbol{r}} \right) $, 因此通过(13)式可以得到位移场为[32]

      其中$ {{\boldsymbol{a}}_1} $和$ {{\boldsymbol{a}}_2} $分别为傅里叶空间中$ {{\boldsymbol{g}}_1} $与$ {{\boldsymbol{g}}_2} $对应实空间的非空线的矢量. 获得位移场后, 采用经典弹性力学可求得应变场[32]:

    3.   模拟结果分析

    • 为研究不同取向差的对称倾侧小角度晶界位错构型以及外加应变下位错运动和反应, 本文先在图1相图所示的相点模拟四方相双晶的弛豫过程, 得到稳定的晶界位错构型; 其次, 对弛豫后的晶体构型施加应变, 观察对应的位错响应, 分析位错运动行为、体系自由能、位错周围应变场等; 最后分析不同位错反应内在机理及与宏观性质的对应关系.

      • 3.1.   四方相双晶的弛豫过程模拟

    • 通过表1中A方案的参数来模拟双晶凝固并得到稳定晶界构型的过程, 如图2所示. 设置中间晶粒与两侧晶粒的取向差为$ \theta ={4}{°} $, 其经历了晶界逐步接触-晶界弛豫后形成稳定的对称倾侧小角度晶界. 中间晶粒两侧小角度晶界各由4个伯氏矢量b相同的位错组成, 并且两两间距相同如图2(d), 两行位错伯氏矢量b方向相反、大小相同, 这些位错为刃型位错, 其半原子面与晶界平行.

      对于不同的取向差弛豫后其稳定的晶界位错构型, 如图3所示. 图3中不同颜色的点表示两侧晶界的位错在模拟区域内的位置, 蓝色与青色的折线分别表示位错数量和位错间距的变化. 图中看到位错在晶界处呈直线排布, 当$ \theta ={6} $°与$ \theta ={8} $°时, 两侧位错相对于区域中心对称, 但$ \theta ={4} $°时, 位错沿区域下部偏移. 随着取向差的增大, 位错数量增加, 位错间距减小, 这与小角度对称倾侧晶界的位错模型相符合:

      进一步观察凝固弛豫过程中体系自由能的变化, 如图4. 随着液相的不断消耗, 体系自由能减小, 随后保持不变, 表明体系达到了稳定状态, 不同取向差呈现相同的趋势. 但取向差越大, 稳定时体系自由能越高, 可认为, 随着取向差增大, 晶界能越大, 而导致整个体系自由能增大.

      为了分析不同取向差下晶界位错附近的应变分布, 现获得晶界处$ {\varepsilon }_{xx} $分量云图, 如图5. 位错应变场呈蝴蝶状范围直线排布在晶界两侧, 随着取向差的增大, 单个位错的应变范围逐渐减小, 并且晶界位错之间的无应变区域也减小. 这表明小角度晶界的取向差与晶界附近的本征应变有着很大的联系.

      • 3.2.   取向差对外加应变下位错运动的影响

    • 本文采用2.2节设置的A, B和C三种方案, 对弛豫后的晶体结构沿x轴施加拉应变, 同时在y轴施加压应变, 分析在此条件下不同晶界取向差θ对对称倾侧晶界位错运动的影响.

      如图6为$ \theta ={4} $°时晶界位错演变. 图中黑色箭头表示位错的运动方向, 蓝色虚线框中表示位错演化过程中对应的特征构型. 为便于表示晶界发射位错类型及后续的位错反应, 在图6(a)标定了四方相的(11)面, 并用红色位错符号在插图(g)中标识四种位错类型$ n1 $, $ n2 $, $ n3 $, $ n4 $, 分别对应的伯氏矢量为: ${{\boldsymbol{b}}}_{1}=\dfrac{a}{2}[01], {{\boldsymbol{b}}}_{2}=\dfrac{a}{2}[0\overline{1}]$, ${{\boldsymbol{b}}}_{3}=\dfrac{a}{2}[10], {{\boldsymbol{b}}}_{4}= $$ \dfrac{a}{2}[\overline{1}0]$, 这里, a为晶体的点阵常数.

      图中看到, 由于晶界位错攀移克服晶界弹性阻碍所需的能量小于滑移, 晶界位错先沿水平方向攀移, 且上下两列位错运动方向相反(Step 1). 随后, 位错附近原子发生预熔化, 当预熔化区域达到一定宽度时, 两侧晶界位错将垂直于晶界方向同时发射出两个大小相等方向相反的刃位错(Step 2).

      从图6(b)可以观察到晶界位错分别为$ n1 $和$ n2 $, 发射位错为$ n3 $和$ n4 $. 故上述晶界位错分解发射为L1型(见表2)位错反应. 随后, 晶界和发射位错各自进行攀移, 当$ n3 $与$ n4 $型位错相遇时发生R1(见表2)型反应, 直到$ n3 $与$ n4 $型位错全部湮灭(图6(c)中的Step 3).

      随后重复“Step 1→Step 2→Step 3”过程, 如图6(d)—(f)所示. 但从图6(f)可以看到下侧晶界位错不再呈直线排布, 而是在垂直方向发生了滑移, 这是由于R1型反应未能使$ n3 $与$ n4 $型位错全部湮灭, 剩余的发射位错和晶界位错随机形成位错组, 发生攀滑移, 而导致下侧晶界位错位置改变

      $ \theta ={6} $°(图7)与$ \theta ={4} $°的位错运动演化相比, 仍按Step 1和Step 2过程进行(图7(a)和图7(b)). 但晶界位错不再同时分解发射出$ n3 $与$ n4 $型位错, 而是产生如图7(b)中新的位错组态: 发生L2(红色圈)和L3(黄色圈)型(见表2)的分解发射, 形成单个位错和位错组; 褐色圈中虽然也为L1构型, 但比较于$ \theta ={4}$°, 不再具有规则的对称性; 绿色圈中的位错构型, 未分解发射. 可能原因是由于取向差增大, 晶界位错数量增加, 当施加外应变时, 各位错受到的应变不再平均化, 不一致的应变导致分解发射不同的位错组态. 在Step 3过程, 位错组发生分解, 而对于单个位错不再只发生攀移, 而是攀移、滑移一起进行, 如图7(b)与图7(c)中所标出的A与B. 随后仅部分$ n3 $与$ n4 $型位错相遇发生R1构型位错反应, 如图7(c)中红色矩形框中所示.

      随着演化的进行, 不再重复“Step 1→Step 2→Step 3”, 而是保持原有的运动状态进行Step 4, 攀滑移时晶界位错和发射位错相遇时两两结合形成位错组, 例如图7(d)中C与D位错, 结合后形成图7(e)中的E位错组. 这是由于发射的位错只有部分发生反应而湮灭, 模拟区域内位错数量较多, 并部分以位错组存在, 各位错受到的应变不足以驱使位错分解发射. 当演化进行到23500步时, 不同位错或位错组相遇发生如表2中标号为R2, R3, R4, R5的位错反应, 如图7(e)中红色方框中所示, 为Step 5. 通过上述的位错反应, 位错发生湮灭, 如图7(f)所示, 位错数量减少, 体系畸变减小, 取向差趋于统一, 但存在少量位错, 未能形成完美单晶状.

      如图8所示, 当晶界两侧晶粒取向差增大到$ \theta ={8} $°时, Step 1和Setp 2与$ \theta ={6} $°时相似, 但不存在未分解发射的晶界位错, 并且发射的位错不再优先向垂直方向进行攀移, 而是如蓝色实线矩形中与邻近晶界位错结合形成位错组, 如图8(b)所示. 这是由于取向差的进一步增大, 导致晶界位错数与分解发射的位错数量进一步增加, 晶界处位错密度增大, 晶界位错与发射位错之间间距减少, 之间相互作用力增加, 从而发射位错优先于晶界位错形成位错组. Step 3时位错组运动中不再分解, 这表明随着取向差的增大, 晶界位错对发射位错的吸引力更大, 当演化到t = 17500时, 部分发射位错或者与位错组之间发生R1和R3构型的反应, 如图8(c).

      部分位错湮灭后继续保持原有状态运动, 位错组和少数单个位错发生攀滑移, 如图8(d), 为Step 4. Step 5主要为位错组发生R4构型的反应, 同时过程中也伴随着R1构型的反应, 如图8(e). 随着位错之间的反应湮灭, 畸变减小, 最终未能形成单晶.

      • 3.3.   不同取向差下体系自由能变化

    • 通过不同时刻自由能密度积分获得了体系自由能随时间的演化(如图9). 可以看到, 体系自由能随时间演化大致分为6个阶段(S1—S6). 当$ \theta ={4} $°时(图9(a)), 对应3.2节原子构型模拟结果, S1阶段对应Step 1过程, 攀移导致畸变能增加, 体系自由能逐渐升高, 直到最高点a点, 并且此时原子也预熔化完毕; S2阶段对应Step2—Step3过程, 位错的分解发射、反应导致攀移时所存储的畸变能得到释放, 体系自由能降低, 当发射位错全部湮灭时, 体系自由能降到最低点b点. S3—S6阶段为上述过程的重复, 但c点下降到d点的后半段时存在能量的起伏波动且未达到与b点相同的最低点, 是由于发生R1型位错反应时发射未全部湮灭.

      而进一步比较$ \theta ={6} $°和$ \theta ={8} $°, 如图9(b)和图9(c)所示. S1—S2阶段能量变化依旧为其Step 1—Step 3过程导致. 而S3—S4阶段能量上升与下降的幅度相较于$ \theta ={4} $°时明显降低, 这是由于其未重复“Step 1→Step 2→Step 3”, 而是Step 4和Step 5过程位错与位错组攀滑移以及位错之间反应湮灭引起自由能的增加与减少. S5和S6阶段, 未湮灭完的位错继续攀移, 体系自由能上升, 当少数位错湮灭时, 能量下降到最低点.

      • 3.4.   位错反应演化过程及应变表征

    • 本节将通过GPA方法构建位错反应类型与应变之间的关系, 探究诱发位错反应的内在因素. 从晶体构型与应变关系图10看到, $ {\varepsilon }_{xx} $和$ {\varepsilon }_{yy} $应变图中并未反映出位错组构型, 单个位错的应变范围较小, 这是因为晶界位错主应变主要分布在水平方向, 而发射位错分布在垂直方向. 因此本节选用切应变$ {\varepsilon }_{xy} $作为参考, 图11和图12为R1—R5构型位错反应演化过程及对应的应变表征, 其反应式可见表2.

      R1为$ n3 $与$ n4 $型位错之间的反应. 演化过程中, 位错在y方向攀移, x方向滑移, 位错之间原子层发生错配, 并呈正切应变(图12(b)). 在正切应变的驱动下, 位错相互靠近, 如图11(a)和图11(b). 当正切应变消失后, 位错之间仅剩两个原子层厚的预熔区, 呈负切应变(图12(c)). 在负切应变的驱动下, 位错湮灭, 局部应变消失, 如图11(d)和图12(d)所示.

      而对于R2反应, 为晶界位错与位错组之间的反应, 位错组中位错呈$ {90^\circ } $, 负切应变区重合, 位错与位错组在y方向滑移, x轴方向攀移, 之间错配的原子层呈正切应变, 中心处较两边应变值大, 如图11(e)和图12(e). 在正切应变驱动下, 位错与位错组距离减小3个原子层厚度后正切应变消失, 剩下的负切应变与单个位错的负切应变区相连(图12(f)). 随后在负切应变的驱动下, 原子层的畸变消失, $ n1 $和$ n2 $型位错湮灭, 得到含有标准应变分布的$ n3 $型位错, 如图11(h)和图12(h).

      R3构型也为位错与位错组的反应, 除了攀滑移方向以及应变场分布是与R2相反, 其过程和应变表征与R2相似, 最后$ n3 $与$ n4 $位错湮灭, 得到含有标准应变分布的$ n1 $位错, 如图11(i)—(l)和图12(i)—(l).

      R4构型位错反应为两位错组间的演化过程, 如图11(m)—(p). 各位错组间由呈90o的两位错组成. 随时间演化, 位错组沿y方向靠近, 并沿x方向滑移, 随位错组间靠近至接近2层晶面间距, 反向的$ n3 $与$ n4 $型位错湮灭, 并留下较大负切应变区, 最终驱动$ n1 $与$ n2 $型位错全部湮灭, 如图11(o)与图11(p)和图12(o)与图12(p).

      R5反应为$ n1 $与$ n2 $型位错之间的反应, 除了攀滑移方向及应变场分布与R1相反, 其过程与应变表征是与R1相似, 最后$ n1 $与$ n2 $型位错湮灭见图11(q)—(t)与图12(q)—(t).

    4.   结　论

    • 本文通过晶体相场方法模拟在外应变作用下不同取向差的对称倾侧小角度晶界的位错运动与反应, 研究结果如下:

      1)在外应变的作用下, 不同取向差初始时均经历晶界位错的攀移、分解发射以及位错反应等过程. 随着取向差增大, 位错密度增加, 产生更多位错组构型, 以及相应位错组之间的反应.

      2)不同取向差下体系自由能变化分为6个阶段. S1阶段, 位错攀移导致能量积累, 体系自由能增加; S2阶段, 位错分解发射以及反应湮灭使能量得到释放, 体系自由能降低; S3—S6阶段为上述过程的重复. 随着取向差的增大, S3—S6阶段不再重复, 而是继续运动的位错及位错组的攀滑移及之间反应湮灭导致体系自由能升高与下降.

      3)不同的构型位错反应中, 位错或位错组在正切应变的驱动下, 通过攀滑移靠近. 当位错或位错组接近时, 正切应变消失, 剩下的负应变驱动位错湮灭.

    参考文献 (37)

目录

/

返回文章
返回