搜索

文章查询

x

留言板

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

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

H对Mg2Si力学性能影响的第一性原理研究

付正鸿 李婷 单美乐 郭糠 苟国庆

H对Mg2Si力学性能影响的第一性原理研究

付正鸿, 李婷, 单美乐, 郭糠, 苟国庆
PDF
HTML
导出引用
导出核心图
  • Al-Mg-Si系铝合金在氢环境服役时, 因遭受氢损伤而导致力学性能退化. Mg2Si是Al-Mg-Si系铝合金主要的热处理强化相, 其力学性能直接决定着Al-Mg-Si系铝合金的强度. 本文采用基于密度泛函数理论的第一性原理计算方法, 研究了间隙H原子对Mg2Si力学性能的影响. 首先计算了Mg2Si的单晶体弹性常数C11, C12C44, 在通过Hill模型计算了多晶体的弹性模量、泊松比及硬度, 并进一步计算了H对Mg2Si晶粒拉伸性能的影响, 最后对H原子掺杂引起的Mg2Si晶体态密度变化进行了分析. 结果表明: H原子的引入显著降低Mg2Si晶体的剪切模量和弹性模量, 从而使得体系强度、硬度降低, 韧性提高. 拉伸性能计算也表明H原子引起Mg2Si晶体断裂强度降低、断裂延伸率提高. 态密度分析表明, Mg2Si晶体中引入氢原子, 将引起Mg2Si晶体由半导体性质向金属性质的转变, 从而造成体系强度、硬度的降低, 韧性增加. 本文计算结果可以为揭示以Mg2Si为增强相的材料在氢环境中强度降低机理提供参考依据.
      通信作者: 苟国庆, gouguoqing1001@163.com
    • 基金项目: 西南交通大学博士创新基金(批准号: D-CX201831)和国家重点研发计划(批准号: 2016YFB1200602-16)资助的课题.
    [1]

    李建国, 谭红艳, 史子木, 何迁 2008 中国有色金属学报 18 1819

    Li J G, Tan H Y, Shi Z M, He Q 2008 Chin. J. Nonferrous. Met. 18 1819

    [2]

    Zeng F L, Wei Z L, Li J F, Li C X, Tan X, Zhang Z, Zheng Z Q 2011 T. Nonferr. Metal. Soc. 21 2559

    [3]

    Qin Q D, Li W X, Zhao K W, Qiu S L, Zhao Y G 2010 Mat. Sci. Eng. A 527 2253

    [4]

    Tong X, Zhang D, Wang K, Lin J, Liu Y, Shi Z, Li Y, Lin J, Wen C 2018 Mat. Sci. Eng. A 733 9

    [5]

    任玉艳, 刘桐宇, 李英民 2016 中国科学: 物理学 力学 天文学 46 084611

    Ren Y Y, Liu T Y, Li Y M 2016 Sci. Sin.: Phys. Mech. Astron. 46 084611

    [6]

    余本海, 刘墨林, 陈东 2011 物理学报 60 087105

    Yu B H, Liu M L, Chen D 2011 Acta Phys. Sin. 60 087105

    [7]

    韩秀丽, 王清, 孙东立, 张红星 2008 中国有色金属学报 18 523

    Han S L, Wang Q, Sun L D, Zhang H X 2008 Chin. J. Nonferrous. Met. 18 523

    [8]

    张凤春, 李春福, 文平, 罗强, 冉曾令 2014 物理学报 66 227101

    Zhang F C, Li C F, Wen P, Rang Z L 2014 Acta Phys. Sin. 66 227101

    [9]

    马明光, 亢世江, 张红玲, 徐红彬 2015 热加工工艺 44 96

    Ma G M, Kang S J, Zhang H L, Xu H B 2015 Hot Working Technology 44 96

    [10]

    姚宝殿, 胡桂青, 于治水 2016 物理学报 65 026202

    Yao B D, Hu G Q, Yu Z S 2016 Acta Phys. Sin. 65 026202

    [11]

    韩秀丽 2010 博士学位论文 (哈尔滨: 哈尔滨工业大学)

    Han X L 2010 Ph. D. Dissertation (Harbin: Harbin Institute of Technology) (in Chinese)

    [12]

    饶建平, 欧阳楚英, 雷敏生 2012 物理学报 61 047105

    Rao J P, Ouyang C Y, Lei M S 2012 Acta Phys. Sin. 61 047105

    [13]

    石瑜, 白洋, 莫丽玢, 向青云, 黄亚丽, 曹江利 2015 物理学报 64 116301

    Shi Y, Bai Y, Mo L B, Xiang Q Y, Huang Y L, Cao J L 2015 Acta Phys. Sin. 64 116301

    [14]

    刘娜娜, 孙翰英, 刘洪生 2009 材料导报: 纳米与新材料专辑 23 278

    Liu N N, Sun H Y, Liu H S 2009 Mat. Rev. 23 278

    [15]

    Li L Y, Yu W, Jin C Q 2005 J. Phys.: Condens. Mat. 17 5965

    [16]

    Hill R 2002 Proc. Phys. Soc. 65 349

    [17]

    Jang S H, Chichibu S F 2012 J. Appl. Phys. 112 073503

    [18]

    Yu W Y, Wang N, Xiao X B 2009 Solid State Sci. 11 1400

    [19]

    Li Y, Gao Y, Xiao B 2010 J. Alloy Compd. 502 28

    [20]

    Bystricky J, LaFrance P, Lehar F, Perrot F, Winternitz P 1985 Phys. Rev. D 32 575

    [21]

    Parks G S 1973 J. Chem. Educ. 30 82A

    [22]

    Sung C M, Sung M 1996 Mater. Chem. Phys. 43 1

    [23]

    Leger J M, Djemia P, Ganot F 2001 Appl. Phys. Lett. 79 2169

    [24]

    Miao N, Sa B, Zhou J 2011 Comp. Mater. Sci. 50 1559

    [25]

    Senkov O N, Dubois M, Jonas J J 1996 Metall. Mater. Trans. A 27 3963

    [26]

    Maziyar A, Ahad M, Herbert D, Christian G M 2015 Metall. Mater. Trans. B 46 1471

  • 图 1  H原子位于八面体间隙的Mg2Si-H体系超晶胞结构模型

    Fig. 1.  Super-cell structural model of Mg2Si-H system.

    图 2  弹性常数C11, C12, C44C'随氢原子浓度增加的变化趋势

    Fig. 2.  Relationships between C11, C12, C44 and C' as a function of H atoms content.

    图 3  剪切模量(GH)、体模量(BH)、弹性模量(E)、泊松比(ν)、脆性(BH/GH)及硬度(H)与体系H原子浓度的变化关系

    Fig. 3.  Shear moduli G, bulk moduli B, Young's moduli E, Poisson ratios ν, brittleness BH/GH and Hardness H as a function of H atoms content.

    图 4  Mg2Si 1 × 1 × 2超晶胞和4 at.% Mg2Si-H超晶胞沿[001]方向进行拉伸的应力-应变曲线

    Fig. 4.  Stress vs. strain curves of 1 × 1 × 2 super cell with 4 at.% and without H during stretching along [001] direction

    图 5  (a) Mg2Si 和 (b) Mg2Si-7.69 at.%H 总态密度和分波态密度

    Fig. 5.  DOS and PDOS of (a) Mg2Si and (b) Mg2Si-7.69 at.% H.

    表 1  Mg2Si-H体系晶体弹常数计算结果

    Table 1.  Elastic constant of Mg2Si-H system.

    PhaseSourceElastic constants/GPa
    C11C12C44C'C12C44
    Mg2SiThis work113.222.245.245.5–23.0
    Calculated [20]113.722.843.5
    Calculated [6]114.521.545.6
    Experimental [21]126.026.048.5
    Mg2Si-2.05 at.%HThis work110.325.138.842.6–13.7
    Mg2Si-4 at.%HThis work109.026.137.841.4–11.7
    Mg2Si-7.69 at.%HThis work90.635.534.927.50.6
    下载: 导出CSV

    表 2  Mg2Si-H体系模量

    Table 2.  Modulus of Mg2Si-H system.

    PhaseSourceModulus
    GH/GPaBH/GPaE/GPaνBH/GHH/GPa
    Mg2SiThis work45.352.5105.60.161.1610.31
    Calculated [20]46.252.5107.10.16
    Calculated [6]44.353.4104.10.17
    Experimental [21]59.0120.0
    Mg2Si-2.05 at.%HThis work40.253.496.50.171.329.07
    Mg2Si-4 at.%HThis work39.253.794.70.181.378.56
    Mg2Si-7.69 at.%HThis work31.753.879.50.281.694.55
    下载: 导出CSV
  • [1]

    李建国, 谭红艳, 史子木, 何迁 2008 中国有色金属学报 18 1819

    Li J G, Tan H Y, Shi Z M, He Q 2008 Chin. J. Nonferrous. Met. 18 1819

    [2]

    Zeng F L, Wei Z L, Li J F, Li C X, Tan X, Zhang Z, Zheng Z Q 2011 T. Nonferr. Metal. Soc. 21 2559

    [3]

    Qin Q D, Li W X, Zhao K W, Qiu S L, Zhao Y G 2010 Mat. Sci. Eng. A 527 2253

    [4]

    Tong X, Zhang D, Wang K, Lin J, Liu Y, Shi Z, Li Y, Lin J, Wen C 2018 Mat. Sci. Eng. A 733 9

    [5]

    任玉艳, 刘桐宇, 李英民 2016 中国科学: 物理学 力学 天文学 46 084611

    Ren Y Y, Liu T Y, Li Y M 2016 Sci. Sin.: Phys. Mech. Astron. 46 084611

    [6]

    余本海, 刘墨林, 陈东 2011 物理学报 60 087105

    Yu B H, Liu M L, Chen D 2011 Acta Phys. Sin. 60 087105

    [7]

    韩秀丽, 王清, 孙东立, 张红星 2008 中国有色金属学报 18 523

    Han S L, Wang Q, Sun L D, Zhang H X 2008 Chin. J. Nonferrous. Met. 18 523

    [8]

    张凤春, 李春福, 文平, 罗强, 冉曾令 2014 物理学报 66 227101

    Zhang F C, Li C F, Wen P, Rang Z L 2014 Acta Phys. Sin. 66 227101

    [9]

    马明光, 亢世江, 张红玲, 徐红彬 2015 热加工工艺 44 96

    Ma G M, Kang S J, Zhang H L, Xu H B 2015 Hot Working Technology 44 96

    [10]

    姚宝殿, 胡桂青, 于治水 2016 物理学报 65 026202

    Yao B D, Hu G Q, Yu Z S 2016 Acta Phys. Sin. 65 026202

    [11]

    韩秀丽 2010 博士学位论文 (哈尔滨: 哈尔滨工业大学)

    Han X L 2010 Ph. D. Dissertation (Harbin: Harbin Institute of Technology) (in Chinese)

    [12]

    饶建平, 欧阳楚英, 雷敏生 2012 物理学报 61 047105

    Rao J P, Ouyang C Y, Lei M S 2012 Acta Phys. Sin. 61 047105

    [13]

    石瑜, 白洋, 莫丽玢, 向青云, 黄亚丽, 曹江利 2015 物理学报 64 116301

    Shi Y, Bai Y, Mo L B, Xiang Q Y, Huang Y L, Cao J L 2015 Acta Phys. Sin. 64 116301

    [14]

    刘娜娜, 孙翰英, 刘洪生 2009 材料导报: 纳米与新材料专辑 23 278

    Liu N N, Sun H Y, Liu H S 2009 Mat. Rev. 23 278

    [15]

    Li L Y, Yu W, Jin C Q 2005 J. Phys.: Condens. Mat. 17 5965

    [16]

    Hill R 2002 Proc. Phys. Soc. 65 349

    [17]

    Jang S H, Chichibu S F 2012 J. Appl. Phys. 112 073503

    [18]

    Yu W Y, Wang N, Xiao X B 2009 Solid State Sci. 11 1400

    [19]

    Li Y, Gao Y, Xiao B 2010 J. Alloy Compd. 502 28

    [20]

    Bystricky J, LaFrance P, Lehar F, Perrot F, Winternitz P 1985 Phys. Rev. D 32 575

    [21]

    Parks G S 1973 J. Chem. Educ. 30 82A

    [22]

    Sung C M, Sung M 1996 Mater. Chem. Phys. 43 1

    [23]

    Leger J M, Djemia P, Ganot F 2001 Appl. Phys. Lett. 79 2169

    [24]

    Miao N, Sa B, Zhou J 2011 Comp. Mater. Sci. 50 1559

    [25]

    Senkov O N, Dubois M, Jonas J J 1996 Metall. Mater. Trans. A 27 3963

    [26]

    Maziyar A, Ahad M, Herbert D, Christian G M 2015 Metall. Mater. Trans. B 46 1471

  • [1] 胡洁琼, 谢明, 陈家林, 刘满门, 陈永泰, 王松, 王塞北, 李爱坤. Ti3AC2相(A = Si,Sn,Al,Ge)电子结构、弹性性质的第一性原理研究. 物理学报, 2017, 66(5): 057102. doi: 10.7498/aps.66.057102
    [2] 余本海, 刘墨林, 陈东. 第一性原理研究Mg2 Si同质异相体的结构、电子结构和弹性性质. 物理学报, 2011, 60(8): 087105. doi: 10.7498/aps.60.087105
    [3] 王斌, 刘颖, 叶金文. 高压下TiC的弹性、电子结构及热力学性质的第一性原理计算. 物理学报, 2012, 61(18): 186501. doi: 10.7498/aps.61.186501
    [4] 胡玉平, 平凯斌, 闫志杰, 杨雯, 宫长伟. Finemet合金析出相-Fe(Si)结构与磁性的第一性原理计算. 物理学报, 2011, 60(10): 107504. doi: 10.7498/aps.60.107504
    [5] 石瑜, 白洋, 莫丽玢, 向青云, 黄亚丽, 曹江利. H掺杂α-Fe2O3的第一性原理研究. 物理学报, 2015, 64(11): 116301. doi: 10.7498/aps.64.116301
    [6] 林俏露, 李公平, 许楠楠, 刘欢, 王苍龙. 金红石TiO2本征缺陷磁性的第一性原理计算. 物理学报, 2017, 66(3): 037101. doi: 10.7498/aps.66.037101
    [7] 范开敏, 杨莉, 孙庆强, 代云雅, 彭述明, 龙兴贵, 周晓松, 祖小涛. 六角相ErAx (A=H, He)体系弹性性质的第一性原理研究. 物理学报, 2013, 62(11): 116201. doi: 10.7498/aps.62.116201
    [8] 周平, 王新强, 周木, 夏川茴, 史玲娜, 胡成华. 第一性原理研究硫化镉高压相变及其电子结构与弹性性质. 物理学报, 2013, 62(8): 087104. doi: 10.7498/aps.62.087104
    [9] 付现凯, 陈万骐, 姜钟生, 杨波, 赵骧, 左良. Ti3O5弹性、电子和光学性质的第一性原理研究. 物理学报, 2019, 68(20): 207301. doi: 10.7498/aps.68.20190664
    [10] 赵荣达, 朱景川, 刘勇, 来忠红. FeAl(B2) 合金La, Ac, Sc 和 Y 元素微合金化的第一性原理研究. 物理学报, 2012, 61(13): 137102. doi: 10.7498/aps.61.137102
  • 引用本文:
    Citation:
计量
  • 文章访问数:  111
  • PDF下载量:  4
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-04
  • 修回日期:  2019-06-15
  • 上网日期:  2019-11-26
  • 刊出日期:  2019-09-01

H对Mg2Si力学性能影响的第一性原理研究

  • 西南交通大学, 先进材料教育部重点实验室, 成都 610031
  • 通信作者: 苟国庆, gouguoqing1001@163.com
    基金项目: 西南交通大学博士创新基金(批准号: D-CX201831)和国家重点研发计划(批准号: 2016YFB1200602-16)资助的课题.

摘要: Al-Mg-Si系铝合金在氢环境服役时, 因遭受氢损伤而导致力学性能退化. Mg2Si是Al-Mg-Si系铝合金主要的热处理强化相, 其力学性能直接决定着Al-Mg-Si系铝合金的强度. 本文采用基于密度泛函数理论的第一性原理计算方法, 研究了间隙H原子对Mg2Si力学性能的影响. 首先计算了Mg2Si的单晶体弹性常数C11, C12C44, 在通过Hill模型计算了多晶体的弹性模量、泊松比及硬度, 并进一步计算了H对Mg2Si晶粒拉伸性能的影响, 最后对H原子掺杂引起的Mg2Si晶体态密度变化进行了分析. 结果表明: H原子的引入显著降低Mg2Si晶体的剪切模量和弹性模量, 从而使得体系强度、硬度降低, 韧性提高. 拉伸性能计算也表明H原子引起Mg2Si晶体断裂强度降低、断裂延伸率提高. 态密度分析表明, Mg2Si晶体中引入氢原子, 将引起Mg2Si晶体由半导体性质向金属性质的转变, 从而造成体系强度、硬度的降低, 韧性增加. 本文计算结果可以为揭示以Mg2Si为增强相的材料在氢环境中强度降低机理提供参考依据.

English Abstract

    • Mg2Si金属间化合物是Al-Mg-Si系铝合金主要的热处理强化相, 高弹性模量、高硬度的Mg2Si颗粒弥散分布于铝合金基体中, 通过阻碍位错运动在Al-Mg-Si系铝合金形成弥散强化作用, 使铝合金强度硬度大幅度提高[1,2]. 由于Mg2Si的高熔点(1085 ℃)[3], 被用于替代SiC陶瓷制备铝基复合材料, 且具有良好的切削加工性和成形性. 因此, Mg2Si金属间化合物的研究备受关注[4].

      对于Mg2Si金属间化合物的理论研究, 主要采用基于密度泛函数理论的第一性原理计算方法. 任玉艳等[5]对Mg2Si金属间化合物的结构稳定性、热力学和力学性能进行了系统的第一性原理计算研究, 理论计算值与实验值十分接近. 余本海等[6]采用第一性原理的方法对Mg2Si的同质异相体的晶体结构、电子结构及弹性性质进行了全面研究.

      氢由于其原子半径极小, 容易在材料中扩散, 并引起材料性能的变化. 氢原子对晶体结构和力学性质影响的第一性原理计算已经成为目前研究的热点. 韩秀丽等[7]采用第一性原理的方法对Ti-H体系的晶体结构和能量进行了研究, 氢原子的加入引起Ti晶体的晶格畸变, 降低了晶体的结合能. 氢在α-Fe和γ-Fe中掺杂均使晶体结构产生不同程度的畸变, 并降低了材料的塑性[8,9]. 姚宝殿等[10]的研究显示氢对Ti3SiC2的力学性质影响不大. 对α-Ti—H键重叠集居数的研究发现[11], 氢减弱了近邻原子间的成键强度, 存在弱键效应; 而在β-Ti—H中, 氢增加了近邻钛原子间的键重叠集居数, 存在强键效应. 氢的加入也会使Nb晶体发生类似的电子结构变化[12]. α-Fe2O3—H中, H和O之间成键, 费米面附近出现了新的杂化能级, 价带顶变宽, 同时带隙由于导带底下移而减小, 表明氢可以用来调节能带结构[13].

      本文采用基于密度泛函数理论的第一性原理方法, 构建Mg2Si-H体系, 研究了氢原子浓度变化对Mg2Si金属间化合物弹性性质的影响, 并对其力学性能的变化进行了深入分析.

    • Mg2Si金属间化合物具有反萤石结构, 空间群号为Fm3m (No. 255), 为面心立方 (FCC)结构, 晶格常数a = b = c = 0.676 nm, α = β = γ = 90° [14]. Si原子占据顶角和面心位置, Mg原子位于Si原子组成的四面体间隙中形成边长为0.5a的简单立方结构, 余留体心位置的八面体间隙, 成为H等小半径原子优先固溶位置.

      本文采用基于密度泛函数理论(DFT)和平面波赝势的CASTEP(Cambridge sequential total energy package)模块, 对不同H原子浓度的Mg2Si-H体系进行模拟计算. 首先构建Mg2Si晶胞, 几何优化后构建1 × 1 × 1, 1 × 1 × 2, 1 × 2 × 2的Mg2Si超晶胞, 在体心八面体间隙中掺杂H原子, 形成H原子掺杂浓度分别为7.69 at.%, 4 at.%和2.05 at.%的Mg2Si-H体系, 如图1所示. 对掺杂H原子后的超晶胞再进行几何优化, 得到稳定的Mg2Si-H体系, 然后对Mg2Si-H体系进行弹性性质计算. 为更直观地反映H原子掺杂对Mg2Si力学性的影响, 分别计算了1 × 1 × 2的Mg2Si超晶胞和1 × 1 × 2的Mg2Si-H (4 at.%)超晶胞沿[001]方向进行拉伸的应力-应变曲线. 拉伸应变通过改变晶胞[001]方向的尺寸进行施加, 几何优化后获得对应的应力值.

      图  1  H原子位于八面体间隙的Mg2Si-H体系超晶胞结构模型

      Figure 1.  Super-cell structural model of Mg2Si-H system.

      计算过程采用周期性边界条件, 结构优化和弹性性质计算过程中, 采用广义梯度近似(genaral gradient approximation, GGA)中的质子平衡方程(proton balance equation, PEB)来处理电子间的交互关联能, 势函数选用倒异空间中表述的超软赝势(ultrasoft pseudopotential). 计算中平面波截断能量(cut off energy)取400 eV, 最大内应力为1 × 10–3 GPa, 原子最大位移为2 × 10–4 nm, 能量收敛标准为1 × 10–6 eV/atom. 第一布里渊区k-point矢量分别取9 × 9 × 9, 9 × 9 × 4和9 × 4 × 4.

      在三向应力下, 材料共有36个弹性常数Cij (i = 1—6, j = 1—6). Mg2Si晶体为立方结构, 因此存在三个独立的弹性常数, 即C11, C12C44, 分别体现体系在[100], [111]和[011]三个方向的抗变形能力. 当C11, C12C44同时满足C11C12 > 0, C11 > 0, C44 > 0和C11+2C12 > 0四个条件时, 合金体系才可以稳定存在, 满足力学性能稳定性条件[15]. 立方晶系剪切模量C'可由C11C12计算:

      $ C' = \frac{{{C_{11}} - {C_{12}}}}{2}. $

      材料的力学性质在微观上可以通过弹性常数、体模量、杨氏模量、剪切模量、泊松比和体模量/剪切模量之比等体现. 多晶弹性模量通常可由单晶弹性常数计算获得, 计算方法可用Voigt模型和Reuss模型描述. 经验证, Voigt模型和Reuss模型所得模量通常为多晶模量的上限和下限, 因此Hill[16,17]对两者进行数学平均, 发现所得值与实测值一致性更好. 因此, 实际应用中Hill模型更为常用. 因此本文中Mg2Si体模量BH, 剪切模量GH和杨氏模量E的计算如下:

      $ {B_{\rm{V}}} = {B_{\rm{R}}} = \frac{1}{3}\left( {{C_{11}} + 2{C_{12}}} \right), $

      $ {G_{\rm{V}}} = \frac{1}{5}\left( {{C_{11}} - {C_{12}} + {C_{44}}} \right), $

      $ {G_{\rm{R}}} = \frac{{5{C_{44}}\left( {{C_{11}} - {C_{12}}} \right)}}{{4{C_{44}} + 3\left( {{C_{11}} - {C_{12}}} \right)}}, $

      其中BVGV分别为Voigt模型的体积模量和剪切模量, BRGR分别为 Reuss模型的体积模量和剪切模量. Hill模型中体积模量B和剪切模量G分别为:

      $ {B_{\rm{H}}} = \frac{{{B_{\rm{V}}} + {B_{\rm{R}}}}}{2}, $

      $ {G_{\rm{H}}} = \frac{{{G_{\rm{V}}} + {G_{\rm{R}}}}}{2}. $

      杨氏模量E和泊松比ν可由BG计算得到[18,19]:

      $ E = \frac{{9{B_{\rm{H}}}{G_{\rm{H}}}}}{{3{B_{\rm{H}}} + {G_{\rm{H}}}}}, $

      $ \nu = \frac{{3{B_{\rm{H}}} - 2{G_{\rm{H}}}}}{{2\left( {3{B_{\rm{H}}} + {G_{\rm{H}}}} \right)}}. $

      通常情况下, B反映材料抗体积压缩能力; G反映材料抗剪切应变的能力, 且与硬度相关[19]; E反映材料抗正应变的能力, 也反映材料的硬度, E越大, 材料的硬度越高[20]. 剪切应力下, ν越小则稳定性越好. BG的比值是材料工程应用中必须考虑的重要参数, 能反映材料的脆韧性, 以1.75作为临界值. 当 B/G > 1.75时材料呈韧性, 否则呈脆性[20,21].

      材料的硬度与剪切模量、体模量之间存在直接的对应关系[22,23]. Miao等 [24]在通过第一性原理研究材料的力学性质中, 提出了硬度H与杨氏模量E和泊松比ν之间计算的经验公式:

      $ H = \frac{{\left( {1 - 2\nu } \right)E}}{{6\left( {1 + \nu } \right)}}. $

    • Mg2Si-H体系晶体弹性常数C11, C12, C44C'计算结果如表1所列. 从表1中可以看出, 所有体系中C11, C12C44均同时满足C11C12 > 0, C11 > 0, C44 > 0和C11+2C12 > 0四个条件, 说明本文所建立的Mg2Si-H体系可以稳定存在, 满足力学性能稳定性条件[15]. 本工作中的Mg2Si的弹性常数C11, C12C44计算值分别为113.2, 22.2和45.2 GPa, 与其他学者的理论计算结果[6,20]相比, 最大计算误差分别为1.15%, 2.70%和0.88%, 计算精度较高; 与其他学者的试验结果[21]相比, 本工作中的, C12C44的最大计算误差分别为10.16%, 14.61和6.80%, 说明本文计算方法具有较高的可靠性.

      PhaseSourceElastic constants/GPa
      C11C12C44C'C12C44
      Mg2SiThis work113.222.245.245.5–23.0
      Calculated [20]113.722.843.5
      Calculated [6]114.521.545.6
      Experimental [21]126.026.048.5
      Mg2Si-2.05 at.%HThis work110.325.138.842.6–13.7
      Mg2Si-4 at.%HThis work109.026.137.841.4–11.7
      Mg2Si-7.69 at.%HThis work90.635.534.927.50.6

      表 1  Mg2Si-H体系晶体弹常数计算结果

      Table 1.  Elastic constant of Mg2Si-H system.

      图2给出了弹性常数C11, C12, C44C'随氢原子浓度增加的变化趋势. C11值反映了体系在[100]方向的抗变形能力. 从图2中可以看出, 随氢原子浓度从0增加到7.69 at.%过程中, C11从113.2 GPa逐步减小到90.6 GPa, 降低了20%, 说明H原子的引入显著降低了Mg2Si在[100]方向抵抗变形的能力. C12反映了体系[111]方向的抗变形能力. 从图2中可看出, H原子浓度从0增加到7.69 at.%过程中, Mg2Si的C12值从22.2 GPa增加到35.5 GPa, 增加了59.9%, 说明氢原子的引入使得Mg2Si抵抗[111]方向的抗变形能力明显提高. C44反映了体系[011]方向的抗变形能力, 随H原子浓度从0增加到7.69 at.%, Mg2Si的C44从45.2 MPa降低到34.9 GPa, 降低了22.8%, 说明Mg2Si在[011]方向抗变形能力随氢原子浓度的增加而有所降低. 经(1)式计算得到的单晶体剪切模量C'则随着氢浓度从0增加到7.69 at.%, 从45.5 GPa降低到27.5 GPa, 降低了39.5%, 说明H原子具有显著降低Mg2Si单晶体的抗剪切变形能力. 此外, C12C44值还可表征晶体的脆性, 当C12C44 > 0时晶体为塑性, 当C12C44 < 0晶体为脆性[17]. 从表1可以看出, 当H原子浓度达到7.69 at.%时, C12C44值从由负值增加到正值, Mg2Si晶体发生了由脆性向塑性的转变.

      图  2  弹性常数C11, C12, C44C'随氢原子浓度增加的变化趋势

      Figure 2.  Relationships between C11, C12, C44 and C' as a function of H atoms content.

    • 经(5)—(9)式计算的剪切模量(GH)、体模量(BH)、弹性模量(E)、泊松比(ν)、脆性(BH/GH)及硬度(H)如表2所列. 本工作中的GH, BH, Eν值分别为45.3, 52.5, 105.6 GPa和0.16, 与其他学者的理论计算结果[6,20]相比, 本工作的最大计算误差分别为2.21%, 1.73%, 1.42%和4.29%, 计算精度较高; 与其他学者试验结果[21]相比, 本工作中的BHE的最大计算误差分别为12.38%和13.63%, 与试验结果同样十分接近, 说明本文计算方法具有较高的可靠性.

      PhaseSourceModulus
      GH/GPaBH/GPaE/GPaνBH/GHH/GPa
      Mg2SiThis work45.352.5105.60.161.1610.31
      Calculated [20]46.252.5107.10.16
      Calculated [6]44.353.4104.10.17
      Experimental [21]59.0120.0
      Mg2Si-2.05 at.%HThis work40.253.496.50.171.329.07
      Mg2Si-4 at.%HThis work39.253.794.70.181.378.56
      Mg2Si-7.69 at.%HThis work31.753.879.50.281.694.55

      表 2  Mg2Si-H体系模量

      Table 2.  Modulus of Mg2Si-H system.

      图3给出了剪切模量(GH)、体模量(BH)、弹性模量(E)、泊松比(ν)、脆性(BH/GH)及硬度(H)随氢原子含量增加的变化趋势. 弹性模量E是表征材料刚度的重要参数, E值越高, 材料刚度越大. 本工作计算结果表明, 未引入氢原子时, Mg2Si的弹性模量E 为105.6 GPa, 随着体系内H原子浓度从0增加到7.69 at.%, 弹性模量E逐步降低到79.5 GPa, 降低了24.7%. 切变模量G与弹性模量E类似, 随着体系氢原子浓度从0增加到7.69 at.%, 切变模量G从45.3 GPa降低到31.7 GPa, 降低了30%. 体模量B则随体系氢浓度的改变未表现出明显的变化. 弹性模量E和切变模量G的大幅度降低, 表明Mg2Si-H体系内氢原子浓度的增加, 显著降低了材料的强度、硬度. 与α-Ti相似[25], Mg2Si也表现出H致软化现象. 硬度计算结果显示, 体系的硬度从不含H原子时的10.31 GPa逐步降低到7.69 at.% H浓度时的4.55 GPa, 降低了55.8%, 印证了弹性模量E和切变模量G的计算结果. 随着H原子浓度增加, 体系硬度降低, B/G值逐步增加, 泊松比ν也逐步增大. 当H原子浓度达到7.69 at.%时, 体系的泊松ν为0.28, 与普通低合金钢基本相当[26]. 说明当大量H原子溶于Mg2Si时, Mg2Si将逐步失去其本身具有的高强度、高硬度和高脆性, 可能严重降低以Mg2Si为增强相材料的强度.

      图  3  剪切模量(GH)、体模量(BH)、弹性模量(E)、泊松比(ν)、脆性(BH/GH)及硬度(H)与体系H原子浓度的变化关系

      Figure 3.  Shear moduli G, bulk moduli B, Young's moduli E, Poisson ratios ν, brittleness BH/GH and Hardness H as a function of H atoms content.

    • 图4给出了1 × 1 × 2的Mg2Si超晶胞和1 × 1 × 2的Mg2Si-H(4 at.%)超晶胞沿[001]方向进行拉伸的应力-应变曲线. 从图中可以看出, 当Mg2Si中掺杂4 at.%的H原子时, 其抗拉强度从15.21 GPa降低到14.12 GPa, 延伸率则从0.3增加到0.36, 说明H原子的掺杂使Mg2Si的强度降低, 塑性增大. 此外, 掺杂4 at.% H原子的Mg2Si在弹性阶段的曲线斜率明显比不掺杂时要小, 说明H原子的掺杂使Mg2Si的弹性模量降低. 拉伸过程的应力-应变曲线结果与弹性模量计算结果得到很好地相互应证.

      图  4  Mg2Si 1 × 1 × 2超晶胞和4 at.% Mg2Si-H超晶胞沿[001]方向进行拉伸的应力-应变曲线

      Figure 4.  Stress vs. strain curves of 1 × 1 × 2 super cell with 4 at.% and without H during stretching along [001] direction

    • 为探索H原子对Mg2Si力学性能影响微观机制, 对体系的态密度和分波态密度进行了分析, 如图5所示. 图5(a)为Mg2Si晶体的电子态密度和分波态密度. 由图可见, Mg2Si晶体的能带主要密集分布于费米能级两侧, 具有典型的半导体特性. 分波态密度显示, 在–9 — –6 eV能量范围内, Mg2Si的态密度主要由Mg的3s和Si的3s态电子贡献, Si的3p态电子有少量贡献. 在–5—0 eV能量范围内, Mg2Si的态密度主要由Mg的3p和Si的3s态电子贡献, Si的3p态电子有少量贡献. 在0—20 eV能量范围内, Mg2Si的态密度主要由Si的3p和3s态电子贡献, Mg的3p态电子有少量贡献. 当体系中引入7.79 at.%的H原子时, H原子显著增加了Mg和Si原子S和P态电子的贡献率, 总态密度明显提高; 同时, 体系在0—20 eV能量范围内的态密度峰缩减到0—6 eV范围内, 变现为体系的禁带变窄, 如图5(b)所示. 于此同时, H原子的引入, 还引起了费米能级处态密度的增加, 态密度跨过费米能级区域. 因此, H原子的引入使体系具有了一定的金属性质, 从而引起材料强度、硬度的降低, 韧性增加.

      图  5  (a) Mg2Si 和 (b) Mg2Si-7.69 at.%H 总态密度和分波态密度

      Figure 5.  DOS and PDOS of (a) Mg2Si and (b) Mg2Si-7.69 at.% H.

    • 1) Mg2Si-H体系内H原子浓度从0增加到7.69 at.%, 剪切模量和弹性模量分别下降24.7%和31.7%, 从而造成体系强度硬度降低, 韧性提高.

      2)态密度分析表明, Mg2Si晶体中引入氢原子, 将引起Mg2Si晶体由半导体性质向金属性质的转变, 从而造成体系强度、硬度的降低, 韧性增加.

参考文献 (26)

目录

    /

    返回文章
    返回