搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

体心立方Fe中${ \langle 100 \rangle}$位错环对微裂纹扩展影响的分子动力学研究

梁晋洁 高宁 李玉红

体心立方Fe中${ \langle 100 \rangle}$位错环对微裂纹扩展影响的分子动力学研究

梁晋洁, 高宁, 李玉红
PDF
HTML
导出引用
导出核心图
  • 在辐照环境下, 载能粒子与材料相互作用导致材料中原子移位, 造成辐照损伤. 其中, 由辐照形成的过饱和自间隙原子团簇形成的间隙型位错环, 是体心立方Fe为基的材料中常见的辐照缺陷之一, 其与材料中其他缺陷之间相互作用, 是导致辐照硬化、脆化、肿胀及蠕变等辐照损伤的原因之一. 除此相互作用外, 在材料表面或内部沿晶界、沉积相、惰性气体形成的气泡所导致的微裂纹, 是诱发辐照促进应力腐蚀开裂的重要原因. 因此, 理解辐照条件下间隙型位错环与微裂纹之间的相互作用, 是理解辐照促进应力腐蚀开裂微观机制的重要一步. 在本研究中, 利用分子动力学方法, 模拟了原子尺度微裂纹与间隙型位错环之间的相互作用, 研究了位错环与微裂纹之间的距离、相对位置及位错环尺寸对二者相互作用的影响, 揭示了位错环对微裂纹是否沿滑移面扩展的影响, 发现当二者的相互作用起主导作用时(如在临界水平或垂直距离之内), 形成的以$ \langle 100 \rangle $为主或高密度的1/2$ { \left\langle {111} \right\rangle} $位错网络可以抑制微裂纹沿滑移面的扩展. 当位错环尺寸发生变化时, 只有当位错环位错核与微裂纹尖端相互作用时, 才能抑制微裂纹沿滑移面的扩展. 这些结果为进一步理解辐照应力开裂提供了新的参考.
      通信作者: 李玉红, liyuhong@lzu.edu.cn
    • 基金项目: 国家级-科技部国家磁约束核聚变能发展研究专项项目(2018YFE0308101)
    [1]

    Rana M A 2012 Radit. Meas. 47 50

    [2]

    Ascheron C, Schindler A, Flagmeyer R, Otto G 1989 Nucl. Instrum. Methods B 36 163

    [3]

    Ascheron C, Schindler A, Flagmeyer R, Otto G 1986 Phys. Status Solidi A 96 555

    [4]

    Mikkelson R C, Miller J W, Holland R E, Gemmell D S 1973 J. Appl. Phys. 44 935

    [5]

    Kohyama A, Hishinuma A, Gelles D S, Klueh R L, Dietz W, Ehrlich K 1996 J. Nucl. Mater. 233–237 138

    [6]

    Jung P 2002 J. Nucl. Mater. 301 15

    [7]

    Beyerlein I J, Caro A, Demkowicz M J, Mara N A, Misra A, Uberuaga B P 2013 Mater. Today 16 443

    [8]

    Litwa P, Kurpaska L, Varin R A, Perkowski K, Jagielski J, Jozwiak S, Czujko T 2017 J. Alloys Compd. 711 111

    [9]

    Wei Y P, Liu P P, Zhu Y M, Wang Z Q, Wan F R, Zhan Q 2016 J. Alloys Compd. 676 481

    [10]

    Masters B C 1965 Philos. Mag. 11 881

    [11]

    Scott P 1994 J. Nucl. Mater. 211 101

    [12]

    Was G S, Busby J T 2005 Philos. Mag. 85 443

    [13]

    Fukuya K, Nakano M, Fujii K, Torimaru T, Kitsunai Y 2004 J. Nucl. Sci. Technol. 41 1218

    [14]

    梁晋洁, 高宁, 李玉红 2020 物理学报 69 036101

    Liang J J, Gao N, Li Y H 2020 Acta Phys. Sin. 69 036101

    [15]

    黄依娜, 万发荣, 焦治杰 2011 物理学报 60 036802

    Huang Y N, Wan F R, Jiao Z J 2011 Acta Phys. Sin. 60 036802

    [16]

    杜玉峰, 崔丽娟, 万发荣 2019 工程科学学报 41 1016

    Du Y F, Cui L J, Wan F R 2019 Chin. J. Eng. 41 1016

    [17]

    崔丽娟, 高进, 杜玉峰, 张高伟, 张磊, 龙毅, 杨善武, 詹倩, 万发荣 2016 物理学报 65 066102

    Cui L J, Gao J, Du Y F, Zhang G W, Zhang L, Long Y, Yang S W, Zhan Q, Wan F R 2016 Acta Phys. Sin. 65 066102

    [18]

    邓平, 孙晨, 彭群家, 韩恩厚, 柯伟 2015 中国腐蚀与防护学报 35 479

    Deng P, Sun C, Peng Q J, Han E H, Ke W 2015 J. Chin. Soc. Corros. Prot. 35 479

    [19]

    Gurrutxaga-Lerma B 2018 Int. J. Eng. Sci. 128 165

    [20]

    Gupta J, Hure J, Tanguy B, Laffont L, Lafont M C, Andrieu E 2018 J. Nucl. Mater. 501 45

    [21]

    Decelis B, Argon A S, Yip S 1983 J. Appl. Phys. 54 4864

    [22]

    董杰, 陈学东, 范志超, 江慧峰, 陆守香 2008 金属学报 44 1167

    Dong J, Chen X D, Fan Z C, Jiang H F, Lu S X 2008 Acta Metall. Sin. 44 1167

    [23]

    赵芳欣, 张瑛洁, 张松, 赵祖欣, 尹绍奎, 李德成 1996 金属学报 32 1056

    Zhao F X, Zhang Y J, Zhang S, Zhao Z X, Yin S K, Li D C 1996 Acta Metall. Sin. 32 1056

    [24]

    徐永波, 白以龙, 沈乐天, 薛青, 李环, 孔丹 1995 金属学报 31 485

    Xu Y B, Bai Y L, Shen T L, Xue Q, Li H, Kong D 1995 Acta Metall. Sin. 31 485

    [25]

    陈奇志, 褚武扬, 乔利杰, 王燕斌, 肖纪美 1994 金属学报 30 248

    Chen Q Z, Chu W Y, Qiao L J, Wang Y B, Xiao J M 1994 Acta Metall. Sin. 30 248

    [26]

    Ackland G J, Mendelev M I, Srolovitz D J, Han S, Barashev A V 2004 J. Phys. Condens. Matter 16 S2629

    [27]

    Mendelev M I, Han S, Srolovitz D J, Ackland G J, Sun D Y, Asta M 2003 Philos. Mag. 83 3977

    [28]

    Plimpton S 1995 J. Comput. Phys. 117 1

    [29]

    Stukowski A, Albe K 2010 Modell. Simul. Mater. Sci. Eng. 18 085001

  • 图 1  位错环与微裂纹扩展相互作用模型示意图 (a)无位错环时的微裂纹自由扩展; (b)存在不同位置的位错环时与裂纹扩展相互作用的模型

    Fig. 1.  Schematic of interaction between the interstitial dislocation loop and expansion of micro-crack: (a) The model of micro-crack expansion without the effect of dislocation loop; (b) the model after including the effect of interstitial dislocation loop located at different positions.

    图 2  外应力作用下自由微裂纹扩展的过程 (a)−(c) 1.0, 40.0 和122.4 ps时的形貌

    Fig. 2.  Expansion of free micro-crack under the effect of external stress: (a)−(c) Results at time of 1.0, 40.0, and 122.4 ps, respectively.

    图 3  d = 1.5 nm, R = 1.5 nm, 位错环与微裂纹扩展的相互作用过程, (a)−(c)分别为15, 20和120 ps时的形貌 (a) 反应初始阶段在微裂纹尖端形成凹凸不平; (b)二者反应形成位错网络阶段; (c) 裂纹尖端被位错网络钉扎并保持在XOZ平面

    Fig. 3.  Interaction between interstitial dislocation loop with radius R = 1.5 nm and micro-crack with distance to loop d = 1.5 nm. (a) to (c) are results at time of 15, 20, and 120 ps: (a) The initial stage with the formation of rugged crack tip; (b) the state with the formation of dislocation network after the interaction between loop and crack; (c) the final state with crack tip pinned by dislocation network located on XOZ plane.

    图 4  d = 3.0 nm, R = 1.5 nm, 位错环与微裂纹扩展的相互作用过程, (a)−(c)分别为46, 75和122.6 ps时的形貌 (a)微裂纹在滑移面上开始扩展的初始阶段; (b)形成以1/2$\left\langle {111} \right\rangle $为主的位错网络阶段; (c)微裂纹沿滑移面扩展导致材料断裂的阶段

    Fig. 4.  Interaction between interstitial dislocation loop with radius R = 1.5 nm and micro-crack with distance d = 3.0 nm to loop. (a) to (c) are results at time of 46, 75, and 122.6 ps: (a) The stage at the beginning of crack on slip plane by changing the expansion direction; (b) the state with the formation of dislocation network and 1/2$\left\langle {111} \right\rangle $ segments dominate the network; (c) the final state with crack expansion on slip plane, resulting in the fracture of material.

    图 5  d = 1.5 nm, R = 1.5 nm, l = 1.5 nm, 位错环与微裂纹扩展的相互作用过程 (a)−(e)分别为36.8, 48.6, 50.5, 95.0和122.6 ps时的形貌

    Fig. 5.  Interaction between dislocation loop and micro-crack with d = 1.5 nm, R = 1.5 nm, and l = 1.5 nm: (a) to (e) are results at time of 36.8, 48.6, 50.5, 95.0, and 122.6 ps, respectively.

    图 6  d = 1.5 nm, R = 3.0 nm, 位错环与微裂纹扩展的相互作用过程 (a)−(d)分别为29, 32, 80 和123 ps时的形貌

    Fig. 6.  Interaction between dislocation loop and micro-crack with d = 1.5 nm and R = 3.0 nm: (a) to (d) are results at time of 29, 32, 80, and 123 ps.

    图 7  d =5.5 nm, R = 3.0 nm, 计算122.5 ps后的相互作用结果

    Fig. 7.  Results of interaction between dislocation loop and micro-crack with d = 5.5 nm, R = 3.0 nm after 122.5 ps.

  • [1]

    Rana M A 2012 Radit. Meas. 47 50

    [2]

    Ascheron C, Schindler A, Flagmeyer R, Otto G 1989 Nucl. Instrum. Methods B 36 163

    [3]

    Ascheron C, Schindler A, Flagmeyer R, Otto G 1986 Phys. Status Solidi A 96 555

    [4]

    Mikkelson R C, Miller J W, Holland R E, Gemmell D S 1973 J. Appl. Phys. 44 935

    [5]

    Kohyama A, Hishinuma A, Gelles D S, Klueh R L, Dietz W, Ehrlich K 1996 J. Nucl. Mater. 233–237 138

    [6]

    Jung P 2002 J. Nucl. Mater. 301 15

    [7]

    Beyerlein I J, Caro A, Demkowicz M J, Mara N A, Misra A, Uberuaga B P 2013 Mater. Today 16 443

    [8]

    Litwa P, Kurpaska L, Varin R A, Perkowski K, Jagielski J, Jozwiak S, Czujko T 2017 J. Alloys Compd. 711 111

    [9]

    Wei Y P, Liu P P, Zhu Y M, Wang Z Q, Wan F R, Zhan Q 2016 J. Alloys Compd. 676 481

    [10]

    Masters B C 1965 Philos. Mag. 11 881

    [11]

    Scott P 1994 J. Nucl. Mater. 211 101

    [12]

    Was G S, Busby J T 2005 Philos. Mag. 85 443

    [13]

    Fukuya K, Nakano M, Fujii K, Torimaru T, Kitsunai Y 2004 J. Nucl. Sci. Technol. 41 1218

    [14]

    梁晋洁, 高宁, 李玉红 2020 物理学报 69 036101

    Liang J J, Gao N, Li Y H 2020 Acta Phys. Sin. 69 036101

    [15]

    黄依娜, 万发荣, 焦治杰 2011 物理学报 60 036802

    Huang Y N, Wan F R, Jiao Z J 2011 Acta Phys. Sin. 60 036802

    [16]

    杜玉峰, 崔丽娟, 万发荣 2019 工程科学学报 41 1016

    Du Y F, Cui L J, Wan F R 2019 Chin. J. Eng. 41 1016

    [17]

    崔丽娟, 高进, 杜玉峰, 张高伟, 张磊, 龙毅, 杨善武, 詹倩, 万发荣 2016 物理学报 65 066102

    Cui L J, Gao J, Du Y F, Zhang G W, Zhang L, Long Y, Yang S W, Zhan Q, Wan F R 2016 Acta Phys. Sin. 65 066102

    [18]

    邓平, 孙晨, 彭群家, 韩恩厚, 柯伟 2015 中国腐蚀与防护学报 35 479

    Deng P, Sun C, Peng Q J, Han E H, Ke W 2015 J. Chin. Soc. Corros. Prot. 35 479

    [19]

    Gurrutxaga-Lerma B 2018 Int. J. Eng. Sci. 128 165

    [20]

    Gupta J, Hure J, Tanguy B, Laffont L, Lafont M C, Andrieu E 2018 J. Nucl. Mater. 501 45

    [21]

    Decelis B, Argon A S, Yip S 1983 J. Appl. Phys. 54 4864

    [22]

    董杰, 陈学东, 范志超, 江慧峰, 陆守香 2008 金属学报 44 1167

    Dong J, Chen X D, Fan Z C, Jiang H F, Lu S X 2008 Acta Metall. Sin. 44 1167

    [23]

    赵芳欣, 张瑛洁, 张松, 赵祖欣, 尹绍奎, 李德成 1996 金属学报 32 1056

    Zhao F X, Zhang Y J, Zhang S, Zhao Z X, Yin S K, Li D C 1996 Acta Metall. Sin. 32 1056

    [24]

    徐永波, 白以龙, 沈乐天, 薛青, 李环, 孔丹 1995 金属学报 31 485

    Xu Y B, Bai Y L, Shen T L, Xue Q, Li H, Kong D 1995 Acta Metall. Sin. 31 485

    [25]

    陈奇志, 褚武扬, 乔利杰, 王燕斌, 肖纪美 1994 金属学报 30 248

    Chen Q Z, Chu W Y, Qiao L J, Wang Y B, Xiao J M 1994 Acta Metall. Sin. 30 248

    [26]

    Ackland G J, Mendelev M I, Srolovitz D J, Han S, Barashev A V 2004 J. Phys. Condens. Matter 16 S2629

    [27]

    Mendelev M I, Han S, Srolovitz D J, Ackland G J, Sun D Y, Asta M 2003 Philos. Mag. 83 3977

    [28]

    Plimpton S 1995 J. Comput. Phys. 117 1

    [29]

    Stukowski A, Albe K 2010 Modell. Simul. Mater. Sci. Eng. 18 085001

  • [1] 梁晋洁, 高宁, 李玉红. 表面效应对铁\begin{document}${\left\langle 100 \right\rangle} $\end{document}间隙型位错环的影响. 物理学报, 2020, 69(3): 036101. doi: 10.7498/aps.69.20191379
    [2] 姜少宁, 万发荣, 龙毅, 刘传歆, 詹倩, 大貫惣明. 氦、氘对纯铁辐照缺陷的影响. 物理学报, 2013, 62(16): 166801. doi: 10.7498/aps.62.166801
    [3] 崔丽娟, 高进, 杜玉峰, 张高伟, 张磊, 龙毅, 杨善武, 詹倩, 万发荣. 氢离子辐照纯钒中形成的位错环. 物理学报, 2016, 65(6): 066102. doi: 10.7498/aps.65.066102
    [4] 汪俊, 张宝玲, 周宇璐, 侯氢. 金属钨中氦行为的分子动力学模拟. 物理学报, 2011, 60(10): 106601. doi: 10.7498/aps.60.106601
    [5] 崔振国, 勾成俊, 侯氢, 毛莉, 周晓松. 低能中子在锆中产生的辐照损伤的计算机模拟研究. 物理学报, 2013, 62(15): 156105. doi: 10.7498/aps.62.156105
    [6] 贺新福, 杨文, 樊胜. 论FeCr合金辐照损伤的多尺度模拟. 物理学报, 2009, 58(12): 8657-8669. doi: 10.7498/aps.58.8657
    [7] 朱勇, 李宝华, 谢国锋. 质子对BaTiO3薄膜辐照损伤的计算机模拟. 物理学报, 2012, 61(4): 046103. doi: 10.7498/aps.61.046103
    [8] 冉琴, 王欢, 钟睿, 伍建春, 邹宇, 汪俊. 钨中不同构型的双自间隙原子扩散行为研究. 物理学报, 2019, 68(12): 126701. doi: 10.7498/aps.68.20190310
    [9] 杨斌, 魏烁, 史开元. 基于等效弹性模量的微裂纹-超声波非线性作用多阶段模型. 物理学报, 2017, 66(13): 134301. doi: 10.7498/aps.66.134301
    [10] 曹莉霞, 王崇愚. α-Fe裂纹的分子动力学研究. 物理学报, 2007, 56(1): 413-422. doi: 10.7498/aps.56.413
    [11] 徐威, 兰忠, 彭本利, 温荣福, 马学虎. 微液滴在不同能量表面上润湿状态的分子动力学模拟. 物理学报, 2015, 64(21): 216801. doi: 10.7498/aps.64.216801
    [12] 于涛, 谢红献, 刘波. 温度对镍基单晶高温合金γ/γ'相界面上错配位错运动影响的分子动力学研究. 物理学报, 2011, 60(4): 046104. doi: 10.7498/aps.60.046104
    [13] 许爱国, 卢果, 张广财, 方步青. 有限温度下位错环的脱体现象. 物理学报, 2009, 58(11): 7934-7946. doi: 10.7498/aps.58.7934
    [14] 高云亮, 朱芫江, 李进平. Al辐照损伤初期的第一性原理研究. 物理学报, 2017, 66(5): 057104. doi: 10.7498/aps.66.057104
    [15] 姚宝殿, 胡桂青, 于治水, 张慧芬, 施立群, 沈皓, 王月霞. H,He对Ti3SiC2材料力学性能影响的第一性原理研究. 物理学报, 2016, 65(2): 026202. doi: 10.7498/aps.65.026202
    [16] 吴宜勇, 岳龙, 肖景东, 杨德庄, 何世禹, 蓝慕杰, 胡建民, 张忠卫, 王训春, 钱勇, 陈鸣波. 位移损伤剂量法评估空间GaAs/Ge太阳电池辐照损伤过程. 物理学报, 2011, 60(9): 098110. doi: 10.7498/aps.60.098110
    [17] 齐佳红, 胡建民, 盛延辉, 吴宜勇, 徐建文, 王月媛, 杨晓明, 张子锐, 周扬. 电子辐照下GaAs/Ge太阳电池载流子输运机理研究. 物理学报, 2015, 64(10): 108802. doi: 10.7498/aps.64.108802
    [18] 常晓阳, 尧舜, 张奇灵, 张杨, 吴波, 占荣, 杨翠柏, 王智勇. 基于分布式布拉格反射器结构的空间三结砷化镓太阳能电池抗辐照研究. 物理学报, 2016, 65(10): 108801. doi: 10.7498/aps.65.108801
    [19] 尼启良, 刘世界, 王晓光, 陈 波, 范鲜红, 李 敏. Mo/Si多层膜在质子辐照下反射率的变化. 物理学报, 2008, 57(10): 6494-6499. doi: 10.7498/aps.57.6494
    [20] 王铁山, 张多飞, 陈亮, 律鹏, 杜鑫, 袁伟, 杨迪. 辐照导致硼硅酸盐玻璃机械性能变化. 物理学报, 2017, 66(2): 026101. doi: 10.7498/aps.66.026101
  • 引用本文:
    Citation:
计量
  • 文章访问数:  227
  • PDF下载量:  10
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-29
  • 修回日期:  2020-03-31
  • 刊出日期:  2020-06-05

体心立方Fe中${ \langle 100 \rangle}$位错环对微裂纹扩展影响的分子动力学研究

  • 1. 兰州大学核科学与技术学院, 兰州 730000
  • 2. 中国科学院近代物理研究所, 兰州 730000
  • 3. 山东大学前沿交叉科学青岛研究院, 青岛 266237
  • 通信作者: 李玉红, liyuhong@lzu.edu.cn
    基金项目: 国家级-科技部国家磁约束核聚变能发展研究专项项目(2018YFE0308101)

摘要: 在辐照环境下, 载能粒子与材料相互作用导致材料中原子移位, 造成辐照损伤. 其中, 由辐照形成的过饱和自间隙原子团簇形成的间隙型位错环, 是体心立方Fe为基的材料中常见的辐照缺陷之一, 其与材料中其他缺陷之间相互作用, 是导致辐照硬化、脆化、肿胀及蠕变等辐照损伤的原因之一. 除此相互作用外, 在材料表面或内部沿晶界、沉积相、惰性气体形成的气泡所导致的微裂纹, 是诱发辐照促进应力腐蚀开裂的重要原因. 因此, 理解辐照条件下间隙型位错环与微裂纹之间的相互作用, 是理解辐照促进应力腐蚀开裂微观机制的重要一步. 在本研究中, 利用分子动力学方法, 模拟了原子尺度微裂纹与间隙型位错环之间的相互作用, 研究了位错环与微裂纹之间的距离、相对位置及位错环尺寸对二者相互作用的影响, 揭示了位错环对微裂纹是否沿滑移面扩展的影响, 发现当二者的相互作用起主导作用时(如在临界水平或垂直距离之内), 形成的以$ \langle 100 \rangle $为主或高密度的1/2$ { \left\langle {111} \right\rangle} $位错网络可以抑制微裂纹沿滑移面的扩展. 当位错环尺寸发生变化时, 只有当位错环位错核与微裂纹尖端相互作用时, 才能抑制微裂纹沿滑移面的扩展. 这些结果为进一步理解辐照应力开裂提供了新的参考.

English Abstract

    • 在聚变或裂变堆辐照环境下, 载能粒子与服役材料中的晶格原子相互作用过程中存在能量传递过程. 当二者之间传递的能量超过材料原子的离位阈能时, 将导致材料中原子移位一定的距离形成自间隙缺陷, 并在原来晶格点阵位置形成空位缺陷. 在实际辐照过程中(如裂变和聚变反应堆中), 大量的载能粒子与材料相互作用, 产生足够多的空位和间隙原子等点缺陷, 空位和间隙原子之间发生湮灭反应后, 剩余的空位和间隙原子的浓度也远超过给定温度下的平衡态缺陷浓度, 因此会形成过饱和的间隙原子和空位缺陷. 除了湮灭反应, 这些间隙原子与间隙原子之间及空位与空位之间也会同时发生团簇反应, 形成间隙原子团簇和空位团簇, 随着团簇在辐照过程中尺寸的增加, 逐渐形成间隙型位错环和空洞诱发材料辐照损伤, 如辐照硬化、脆化、肿胀、蠕变及疲劳等[1-9], 影响材料的服役性能. 在目前反应堆结构材料中, 以体心立方(bcc结构) Fe为基的合金已经被广泛应用或者作为重要的候选材料. 这些合金经过长时间服役后形成的间隙型位错环主要有1/2$\left\langle {111} \right\rangle $$\left\langle {100} \right\rangle $两种[10], 其中1/2$\left\langle {111} \right\rangle $位错环具有较低的迁移势垒, 当其尺寸在纳米量级时具有高扩散性; 而$\left\langle {100} \right\rangle $位错环具有较高的迁移势垒, 其迁移过程需要较长的时间, 因此, 在一定时间内可以把$\left\langle {100} \right\rangle $位错环看作固定的二维缺陷, 这种类似沉积相的性质使得$\left\langle {100} \right\rangle $位错环能够对其他缺陷的运动起到一定的钉扎阻止作用, 例如对位错线的钉扎作用, 也可能对辐照诱发的微裂纹产生作用, 如果微裂纹开裂过程中与$\left\langle {100} \right\rangle $位错环相互作用, 其扩展过程可能会受到影响, 因此需要做相应的研究.

      除上述情形外, 位错环与微裂纹相互作用也与辐照促进应力腐蚀开裂相关[11-13], 在裂纹开裂过程中, 在外部应力作用下裂纹及腐蚀介质会与辐照缺陷(如间隙型位错环)相互作用, 因此, 对这些过程的理解同样需要首先理解微裂纹与位错环或其他缺陷的相互作用过程.

      在材料内部, 由氦泡诱发的微裂纹的形成及扩展是造成材料氦脆的主要机制之一[6-9], 在外应力作用下微裂纹在扩展过程中也可能与位错环相互作用, 因此, 需要研究应力作用下材料内部微裂纹扩展与位错环相互作用的相关机制. 人们对材料中微裂纹或位错环单独存在时的演化、微裂纹与简单缺陷的相互作用及对材料的影响做了大量的研究[14-25], 但是, 对应力作用下的微裂纹扩展与间隙型位错环相互作用的机制, 到目前还未见详细的报道, 基于上述二者相互作用的重要性, 本研究工作通过构建外应力作用下微裂纹扩展及间隙型位错环的原子尺度模型, 采用分子动力学方法, 研究Fe中$\left\langle {100} \right\rangle $位错环与微裂纹扩展的相互作用过程, 对其物理机理进行详细研究, 并对辐照损伤的影响做相应的讨论.

    • 本项工作主要研究在bcc结构Fe中低温条件下常见的$\left\langle {100} \right\rangle $间隙型位错环与微裂纹扩展的相互作用过程. 首先对$\left\langle {100} \right\rangle $间隙型位错环与微裂纹的相互作用过程开展研究. 为了计算分析方便, 计算胞的三个方向分别设定为X [100], Y [010]及Z [001], 长度分别约为17.1, 11.4和11.4 nm, 每个计算基胞中含有192100个原子, 根据位错环大小不同总原子数会有相应的变化. 通过构建位错环与外应力作用下的裂纹, 形成本项工作的计算模型, 其示意图如图1所示. 在此模型中, $\left\langle {100} \right\rangle $间隙型位错环伯格斯矢量方向沿X轴, 由于bcc结构中沿$\left\langle {100} \right\rangle $方向的原子层堆垛顺序为AB, 根据伯格斯矢量的大小和位错环半径, 可以沿X方向添加一圆形的AB原子面, 并通过原子面平均位移的方式对构建的位错环结构进行松弛, 使其能够满足之后的分子静力学和分子动力学弛豫. 通过构建此原子模型, 可以使$\left\langle {100} \right\rangle $位错环平面(平行于YOZ面)与微裂纹开裂方向(沿着X方向)相互垂直, 在二者相互作用过程中, 既可以保证二者之间充分作用也能摒弃其他因素的影响; 同时考虑到微裂纹形成时基本都在纳米尺度, 因此, 选择X方向长度约为20 nm, 而YZ方向长度的选择可以保证位错环及微裂纹与其镜像之间没有相互作用. 其中描述微裂纹的参数主要包括微裂纹尖端的位置和微裂纹开裂大小, 在本工作中, 裂纹尖端位置(p)位于XOZ平面上, 上下对称, 裂纹开裂大小由裂纹对应的斜率(k)描述; $\left\langle {100} \right\rangle $位错环伯格斯矢量沿X方向, 其中心位置或者沿X平移, 或者沿Y平移, 因此, 通过改变位错环与裂纹尖端的位置, 可以研究位错环与裂纹尖端之间X方向水平距离(d )、相对位置即Y方向垂直距离(l ) 对位错环与裂纹尖端相互作用的影响. 在模拟中, 位错环尺寸(R)也是影响因素之一. 在本研究中, 这些参数范围如下: R 取值为1.5, 2.0和3.0 nm; 位错环与裂纹尖端的水平距离d取值为1.5—7.5 nm(每隔0.5 nm取一个值, 共13个值); 相对位置l的取值为1.5 nm; 初始斜率 k设定为0.15. 另外, 为了确保裂纹能够在给定的外力下开裂, 我们事先进行了相应的试算, 得到适合本项工作在X面上下两个面每个原子上施加的力. 需要注意的是, 本工作不侧重外力大小对相互作用的影响, 所有的计算都采用了相同的外力大小, 以避免外力大小对相互作用的影响. 这些参数的选取主要基于确保二者之间在原子尺度及分子动力学时间尺度内能够发生反应, 并且能够研究不同参数范围对二者相互作用的影响, 得到的结果可以为实际情况下二者相互作用的理解提供参考. 在模拟过程中, Fe-Fe之间的相互作用由Ackland-2004[26]势函数描述, 此势函数能够准确地反映不同点缺陷之间、不同晶体结构之间的能量差、系统的弹性常数以及位错性质, 并且此类型的势函数也已经被应用到微裂纹和间隙型位错环的研究中. 同时通过更换随机数种子对模拟结果进行了重复计算, 得到了相同的结果, 为了检验势函数的影响, 我们也测试了Mendelev势函数[27], 得到类似的结果. 计算首先采用共轭梯度弛豫, 之后利用分子动力学(MD)弛豫的方法进行, 其中MD过程中温度设定为300 K, 时间步长为1 fs, 总弛豫时间为124—136 ps. 针对材料内部形成的微裂纹, 模拟中考虑了周期性边界条件, 整个计算利用LAMMPS[28]完成, 其中温度控制采用速度调整方法, 采用Verlet方法对运动方程积分, 得到的结果利用Ovito[29]完成原子尺度结果分析.

      图  1  位错环与微裂纹扩展相互作用模型示意图 (a)无位错环时的微裂纹自由扩展; (b)存在不同位置的位错环时与裂纹扩展相互作用的模型

      Figure 1.  Schematic of interaction between the interstitial dislocation loop and expansion of micro-crack: (a) The model of micro-crack expansion without the effect of dislocation loop; (b) the model after including the effect of interstitial dislocation loop located at different positions.

    • 通过分析位错环与裂纹尖端的距离、二者之间的相对位置以及位错环半径的大小对微裂纹扩展的影响, 揭示辐照条件下, 材料在位错环形成后裂纹扩展的机制变化及对力学性能的影响. 需要注意的是, 本工作固定了施加到微裂纹上的剪切应力的大小和方向, 以减少变量影响而集中研究位错环性质对微裂纹扩展的影响.

    • 将材料在非辐照条件下服役时微裂纹的扩展定义为自由微裂纹的扩展, 通过对比此种微裂纹的扩展, 研究辐照形成的间隙型位错环对相同微裂纹扩展的影响.

      计算模型如图1(a)所示, 根据第2节的方法, 得到了自由微裂纹的扩展过程及结果. 研究发现, 当自由微裂纹在外应力下开始扩展, 初始状态如图2(a)所示, 裂纹的尖端在外力作用下逐渐变得凹凸不平, 随着裂纹进一步扩展, 裂纹尖端开始沿着系统的滑移面开始向上或向下扩展, 从而使得微裂纹的对称面从XOZ面变为滑移面(如图2(b)); 随着扩展在滑移面上进行, 裂纹尖端沿其伯格斯矢量方向在裂纹扩展的滑移面上运动, 使得裂纹很容易扩展到材料表面, 造成材料的断裂, 最终结果如图2(c)所示. 对上述过程采用变换随机数种子的多次模拟统计得到相同的结果, 证明自由微裂纹的扩展与相应的滑移面紧密相关, 从而造成了裂纹从原来的对称面改变方向变为沿滑移面扩展. 当裂纹尖端有缺陷时, 如本项工作中的间隙型位错环与微裂纹的相互作用, 可能会影响微裂纹的扩展是否会沿着滑移面, 这些结果在3.23.4节中阐述.

      图  2  外应力作用下自由微裂纹扩展的过程 (a)−(c) 1.0, 40.0 和122.4 ps时的形貌

      Figure 2.  Expansion of free micro-crack under the effect of external stress: (a)−(c) Results at time of 1.0, 40.0, and 122.4 ps, respectively.

    • 当材料中存在辐照形成的间隙型位错环后, 模型如图1(b)所示. 首先研究间隙型位错环位于计算胞的中心时微裂纹的扩展情况. 由于微裂纹的对称平面均是位于XOZ平面上, 因此, 微裂纹与位错环均是以XOZ为对称面. 当R = 1.5 nm, d = 1.5 nm时, 随着模拟的进行, 裂纹在应力作用下裂口逐渐变大, 导致裂纹尖端处变得凹凸不平(如图3(a)所示). 对此过程研究表明, 在微裂纹的初始扩展过程中, 尖端没有形成新的位错线片段, 主要原因可能在于本项工作中, 主要研究微裂纹形核结束后的扩展阶段, 而没有包含形核阶段, 因此, 在裂纹扩展的初始阶段, 没有新的位错线片段的形成. 由于位错环为$\left\langle {100} \right\rangle $间隙型位错环, 在给定的温度和模拟时间内, 其状态基本保持不变没有发生位错环的位移. 在应力作用下, 随着时间的延长, 微裂纹尖端逐渐沿着–X方向扩展, 当时间达到17.7 ps时, 位错环与裂纹尖端发生接触反应, 生成$\left\langle {100} \right\rangle $和1/2$\left\langle {111} \right\rangle $共存的位错网络. 如果没有间隙型位错环, 如3.1节的自由微裂纹扩展, 不会形成$\left\langle {100} \right\rangle $和1/2$\left\langle {111} \right\rangle $共存的位错网络. 这些位错网络的形成, 位错密度的升高, 使得$\left\langle {100} \right\rangle $位错环与形成的位错网络反应加速, 最后导致$\left\langle {100} \right\rangle $位错环完全被吸收, 与裂纹尖端形成包含$\left\langle {100} \right\rangle $和1/2$\left\langle {111} \right\rangle $片段的位错网络, 如图3(b)所示, 图中红色线段为$\left\langle {100} \right\rangle $片段, 长度为5.79 nm, 绿色线段为1/2$\left\langle {111} \right\rangle $片段, 长度为8.34 nm. 同时, 正是位错环与裂纹尖端形成位错网络导致裂纹尖端的扩展速度变慢, 主要原因是形成的位错网络对裂纹尖端扩展的钉扎作用, 抑制了微裂纹的快速扩展, 使得材料表现出更好的抗裂纹扩展的性能, 意味着需要更大的外应力才能使裂纹进一步扩展. 随着反应在外应力下的进行, 裂纹尖端在越过位错网络施加的势垒后, 再次不断向内扩展, 同时, 尖端形貌发生严重形变, 但是$\left\langle {100} \right\rangle $位错片段和1/2$\left\langle {111} \right\rangle $共存的位错片段仍然存在, 如图3(c)所示, $\left\langle {100} \right\rangle $片段为5.99 nm, 1/2$\left\langle {111} \right\rangle $片段为3.18 nm, 裂纹尖端依然处于较中心位置, 而没有像自由微裂纹扩展一样, 发生沿着滑移面的裂纹开裂.

      图  3  d = 1.5 nm, R = 1.5 nm, 位错环与微裂纹扩展的相互作用过程, (a)−(c)分别为15, 20和120 ps时的形貌 (a) 反应初始阶段在微裂纹尖端形成凹凸不平; (b)二者反应形成位错网络阶段; (c) 裂纹尖端被位错网络钉扎并保持在XOZ平面

      Figure 3.  Interaction between interstitial dislocation loop with radius R = 1.5 nm and micro-crack with distance to loop d = 1.5 nm. (a) to (c) are results at time of 15, 20, and 120 ps: (a) The initial stage with the formation of rugged crack tip; (b) the state with the formation of dislocation network after the interaction between loop and crack; (c) the final state with crack tip pinned by dislocation network located on XOZ plane.

      以上分析发现, 位错环与微裂纹尖端发生反应形成位错网络, 会在一定程度上抑制微裂纹的扩展. 因此, 可以通过增加微裂纹和位错环之间的距离, 来研究二者的反应及对微裂纹扩展的影响. 当d由1.5增加到2.0及2.5 nm时, 模拟同样揭示出类似的反应的过程, 当d = 3.0 nm时, 模拟发现随着裂纹的开裂, 在位错环与裂纹尖端发生接触反应之前, 裂纹尖端已经开始沿滑移面开裂 (如图4(a)); 发生接触之后尖端形成了$\left\langle {100} \right\rangle $和1/2$\left\langle {111} \right\rangle $片段, 随着反应的进行, $\left\langle {100} \right\rangle $片段逐渐减少, 整个位错网络转化为以1/2$\left\langle {111} \right\rangle $位错线段为主, 如图4(b), 其中$\left\langle {100} \right\rangle $线段长度为0.59 nm, 1/2$\left\langle {111} \right\rangle $线段长度为11.08 nm. 以1/2$\left\langle {111} \right\rangle $片段为主的位错网络会导致在滑移面上容易发生滑移, 当密度比较低时难以形成钉扎结构, 裂纹尖端在应力场下不断向侧上方沿滑移面方向开裂, 随着位错网络不断反应, 形成一段两端终结于裂纹尖端侧表面、长度为14.42 nm的1/2$\left\langle {111} \right\rangle $位错线段(如图4(c)), 未能抑制裂纹尖端进一步沿滑移面开裂, 造成材料断裂. 可以得出当位错环与裂纹尖端反应形成的位错网络中$\left\langle {100} \right\rangle $片段减少时, 位错网络对裂纹尖端开裂的抑制作用减弱, 因此, 可以推断出$\left\langle {100} \right\rangle $位错片段及密度高的1/2$\left\langle {111} \right\rangle $位错片段在抑制裂纹尖端扩展中可能起到主要作用.

      图  4  d = 3.0 nm, R = 1.5 nm, 位错环与微裂纹扩展的相互作用过程, (a)−(c)分别为46, 75和122.6 ps时的形貌 (a)微裂纹在滑移面上开始扩展的初始阶段; (b)形成以1/2$\left\langle {111} \right\rangle $为主的位错网络阶段; (c)微裂纹沿滑移面扩展导致材料断裂的阶段

      Figure 4.  Interaction between interstitial dislocation loop with radius R = 1.5 nm and micro-crack with distance d = 3.0 nm to loop. (a) to (c) are results at time of 46, 75, and 122.6 ps: (a) The stage at the beginning of crack on slip plane by changing the expansion direction; (b) the state with the formation of dislocation network and 1/2$\left\langle {111} \right\rangle $ segments dominate the network; (c) the final state with crack expansion on slip plane, resulting in the fracture of material.

      d从3.5增加到5.0 nm时, 位错环与裂纹尖端或者形成$\left\langle {100} \right\rangle $和1/2$\left\langle {111} \right\rangle $共存的位错片段抑制微裂纹的扩展, 或形成以低密度的1/2$\left\langle {111} \right\rangle $为主的位错网络, 对微裂纹开裂的抑制作用减弱; 而当d ≥ 5.5 nm后, 二者形成的均是以低密度的1/2$\left\langle {111} \right\rangle $为主的位错网络, 导致微裂纹沿着滑移面扩展而不是保持在XOZ面上. 多次的计算统计得到了类似的结果, 因此得出位错环与裂纹尖端之间的距离影响了二者之间的反应, 当d < 3.0 nm或d > 5.0 nm时, 位错环对微裂纹的抑制或增强或减弱, 但是这种影响在一定范围内(3.0—5.0 nm)并不是完全呈线性关系的, 还需要取决于二者在外应力下反应形成的位错网络的性质及位错网络的密度值是否能够钉扎裂纹的开裂.

    • 除了上述的理想情形外, 位错环与微裂纹的相对位置更多可能存在的情形是二者对称面并不重合, 由于这种不重合, 位错环的位错核与微裂纹的反应也会处在非对称条件下, 二者的相互作用可能会与对称面重合的情形不同. 根据3.1节和3.2节的计算结果, 得出位错环与裂纹尖端的反应不同会导致不同的位错网络形成, 最终影响裂纹尖端的扩展性质. 因此, 本节研究位错环与微裂纹尖端相对位置对微裂纹扩展的影响.

      R = 1.5 nm, d = 1.5 nm, l = 1.5 nm时, 模型如图1(b)所示, 根据第2节中的模拟方法, 得到如下结果. 裂纹在外应力作用下开口逐渐变大, 裂纹尖端周围变得凹凸不平, 分析结果同样表明没有新的位错线片段形成, 原因如3.2节所述. 随着模拟时间的延长, 裂纹尖端沿–X方向扩展, 当模拟时间达到36.8 ps时, 位错环与裂纹尖端发生接触反应, 逐渐由两端终结于裂纹尖端的$\left\langle {100} \right\rangle $位错线段演变为1/2$\left\langle {111} \right\rangle $$\left\langle {100} \right\rangle $共存的位错网络, 随着裂纹尖端的扩展位错网络又演化为只包含$\left\langle {100} \right\rangle $的位错片段, 在95 ps时被尖端完全吸收, 同时, 裂纹尖端继续向侧上方更深处扩展, 最终与自由微裂纹扩展结果类似, 反应过程如图5所示. 当d = 2.0, 2.5 nm直到d = 4.0 nm时, 总能得到类似的过程. 当d增大到4.5 nm之后, 由于微裂纹尖端与位错环距离较远, 微裂纹在外应力作用下的扩展首先是沿着滑移面进行, 而没有与位错环反应, 因此, 整个演化过程与自由微裂纹的扩展类似.

      图  5  d = 1.5 nm, R = 1.5 nm, l = 1.5 nm, 位错环与微裂纹扩展的相互作用过程 (a)−(e)分别为36.8, 48.6, 50.5, 95.0和122.6 ps时的形貌

      Figure 5.  Interaction between dislocation loop and micro-crack with d = 1.5 nm, R = 1.5 nm, and l = 1.5 nm: (a) to (e) are results at time of 36.8, 48.6, 50.5, 95.0, and 122.6 ps, respectively.

      对以上过程和结果分析表明, 当位错环与微裂纹相对位置发生改变后, 二者的相互反应或导致反应形成的位错网络被裂纹扩展所吸收, 没有起到相应的钉扎作用, 并且由于反应的非对称性, 导致微裂纹的扩展最终沿着某一滑移面进行, 造成材料中微裂纹的扩展, 使材料服役性能降低; 或二者距离较远导致反应减弱, 最终使得微裂纹的扩展类似于自由微裂纹扩展的情形. 因此, 位错环是否对微裂纹扩展产生影响, 二者之间的相对位置也是一个影响参量.

    • 除了位错环与微裂纹之间的距离及相对位置等影响因素外, 本工作还模拟研究了位错环尺寸对微裂纹开裂的影响. 根据上面的研究, 如果二者之间的距离比较远, 位错环对微裂纹的扩展影响很小, 微裂纹的扩展同自由微裂纹扩展情形类似, 因此, 本节主要研究二者具有相同的XOZ对称面, 也就是位错环位于计算胞的中心, 并且d在3.0 nm内的情形, 以确保能够通过二者之间的反应揭示位错环尺寸对微裂纹扩展的影响.

      d = 1.5 nm, R = 2.0 nm时, 与3.2节所述类似, 裂纹在应力作用下开口逐渐增大, 裂纹尖端周围变得凹凸不平, 与R = 1.5 nm情形相比, 位错环尺寸的增加使得位错环与裂纹尖端发生接触反应的开始时间更早(约10 ps), 主要原因是位错环尺寸的增加使得其位错核周围应力场增强, 导致与微裂纹反应更早发生. 二者反应后同样形成1/2$\left\langle {111} \right\rangle $$\left\langle {100} \right\rangle $共存的位错网络, 随着反应的进行, 位错网络钉扎在裂纹尖端, 阻隔了尖端向更深处扩展. 当R = 3.0 nm时, 大尺寸的位错环在微裂纹应力场下, 首先发生演化形成1/2$\left\langle {111} \right\rangle $$\left\langle {100} \right\rangle $组成的闭合环, 而之后微裂纹在应力场下扩展然后与形成的1/2$\left\langle {111} \right\rangle $$\left\langle {100} \right\rangle $组成的闭合环相互发生反应, 形成1/2$\left\langle {111} \right\rangle $$\left\langle {100} \right\rangle $共存的位错网络, 随后的演化过程与之前类似, 均形成1/2$\left\langle {111} \right\rangle $$\left\langle {100} \right\rangle $位错网络钉扎裂纹尖端的扩展, 同时保持裂纹尖端在XOZ平面上扩展而不是沿其他滑移面扩展, 当扩展进一步增加时, 裂纹会逐渐沿滑移面开裂, 其相互作用过程如图6所示.

      图  6  d = 1.5 nm, R = 3.0 nm, 位错环与微裂纹扩展的相互作用过程 (a)−(d)分别为29, 32, 80 和123 ps时的形貌

      Figure 6.  Interaction between dislocation loop and micro-crack with d = 1.5 nm and R = 3.0 nm: (a) to (d) are results at time of 29, 32, 80, and 123 ps.

      d = 2.0 nm时, R 由1.5逐渐增加到3.0 nm时, 位错环在与微裂纹相互作用之前, 同样先演化为由两小段1/2$\left\langle {111} \right\rangle $$\left\langle {100} \right\rangle $组成的闭合环, 然后再与微裂纹尖端相互作用, 微裂纹尖端被位错网络钉扎且处于XOZ平面上, 同样, 随着扩展的进行裂纹还会沿滑移面扩展. 需要注意的是, 随着位错环尺寸的增加, 最终形成的钉扎位错网络也越来越复杂, 形成多个1/2$\left\langle {111} \right\rangle $$\left\langle {100} \right\rangle $共存的位错网络, 抑制微裂纹的扩展. 当d ≥ 5.5 nm后, 即使R增加到3.0 nm, 此时微裂纹未与位错环发生接触反应而是沿着滑移面扩展, 如图7所示, 主要原因是当二者之间的距离增加到一定值后, 二者之间的反应减弱, 导致位错环对微裂纹的影响可以忽略, 因此, 微裂纹的扩展与自由微裂纹类似.

      图  7  d =5.5 nm, R = 3.0 nm, 计算122.5 ps后的相互作用结果

      Figure 7.  Results of interaction between dislocation loop and micro-crack with d = 5.5 nm, R = 3.0 nm after 122.5 ps.

    • 利用分子动力学方法, 本项工作在原子尺度模拟了微裂纹与辐照形成的间隙型$\left\langle {100} \right\rangle $位错环之间的相互作用, 研究了不同影响参数, 如位错环与微裂纹尖端的水平距离, 二者之间的相对位置以及位错环尺寸等, 对二者之间的相互作用的影响, 并与没有位错环影响下的自由微裂纹的扩展相比较, 揭示出只有当位错环与微裂纹尖端的相互作用起主导作用时, 通过形成的$\left\langle {100} \right\rangle $和1/2$\left\langle {111} \right\rangle $位错网络抑制微裂纹沿滑移面的扩展, 从而使得微裂纹扩展减弱; 当位错环与微裂纹尖端距离较远时, 二者之间的相互作用较弱, 微裂纹在外部应力场作用下主要以沿着滑移面扩展为主; 当位错环与微裂纹不具有相同的对称面时, 在一定的距离内, 位错环的存在使得微裂纹沿靠近位错环的滑移面扩展增强. 同样, 当位错环尺寸发生变化时, 只有当位错环位错核与微裂纹尖端相互作用时, 才能抑制微裂纹沿滑移面的扩展. 这些结果表明, 间隙型位错环对微裂纹扩展的影响取决于不同的条件, 因此, 对于位错环对辐照促进应力腐蚀开裂的影响, 同样需要考虑在不同局部微观条件下的影响, 才能更好理解辐照促进应力腐蚀开裂.

      感谢中国科学院近代物理研究所王志光研究员的讨论.

参考文献 (29)

目录

    /

    返回文章
    返回