-
Al-Mg-Si系铝合金在氢环境服役时, 因遭受氢损伤而导致力学性能退化. Mg2Si是Al-Mg-Si系铝合金主要的热处理强化相, 其力学性能直接决定着Al-Mg-Si系铝合金的强度. 本文采用基于密度泛函数理论的第一性原理计算方法, 研究了间隙H原子对Mg2Si力学性能的影响. 首先计算了Mg2Si的单晶体弹性常数C11, C12和C44, 在通过Hill模型计算了多晶体的弹性模量、泊松比及硬度, 并进一步计算了H对Mg2Si晶粒拉伸性能的影响, 最后对H原子掺杂引起的Mg2Si晶体态密度变化进行了分析. 结果表明: H原子的引入显著降低Mg2Si晶体的剪切模量和弹性模量, 从而使得体系强度、硬度降低, 韧性提高. 拉伸性能计算也表明H原子引起Mg2Si晶体断裂强度降低、断裂延伸率提高. 态密度分析表明, Mg2Si晶体中引入氢原子, 将引起Mg2Si晶体由半导体性质向金属性质的转变, 从而造成体系强度、硬度的降低, 韧性增加. 本文计算结果可以为揭示以Mg2Si为增强相的材料在氢环境中强度降低机理提供参考依据.The mechanical properties of Al-Mg-Si-type aluminum alloys may degenerate due to the hydrogen damage during servicing in hydrogen environment. The Mg2Si is the main strengthening phases in Al-Mg-Si-type aluminum alloys. Therefore, the mechanical properties of Mg2Si directly determine the strengths of Al-Mg-Si-type aluminum alloys. In this work, the effects of hydrogen atoms on the mechanical properties of Mg2Si are investigated by first principle calculation, which is based on the density function theory. First of all, we calculate the single crystal elasticity constants of C11, C12 and C44. Then the elasticity modulus, Poisson’s ratio and hardness of polycrystalline are calculated by using the crystal elasticity constants. Furthermore, we also calculate the tensile properties of Mg2Si with and without H atoms. The difference between the densities of states with and without H atoms is used to investigate the change of Mg2Si induced by H atoms. The results show that hydrogen atoms significantly reduce the shear modulus and elastic modulus of Mg2Si, resulting in the strength and hardness decreasing, but the toughness increasing. The calculations of tensile properties indicate that H atoms reduce the fracture strength but enhance the fracture elongation of Mg2Si. The analysis of density of states indicates that hydrogen atoms will induce the properties of Mg2Si to transform from semiconductor to metal properties. The calculated results in this paper can provide a reference basis for revealing the mechanism of strength reduction of Mg2Si materials in a hydrogen environment.
-
Keywords:
- the first principles /
- Mg2Si /
- elastic properties /
- density of states
- PACS:
71.20.-b (Electron density of states and band structure of crystalline solids) 71.22.+i (Electronic structure of liquid metals and semiconductors and their Alloys) 43.35.Cg (Ultrasonic velocity, dispersion, scattering, diffraction, and Attenuation in solids; elastic constants) 31.15.E- (Density-functional theory) 1. 引 言
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金属间化合物弹性性质的影响, 并对其力学性能的变化进行了深入分析.
2. 计算方法
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体系超晶胞结构模型Fig. 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, C12和C44, 分别体现体系在[100], [111]和[011]三个方向的抗变形能力. 当C11, C12和C44同时满足C11–C12 > 0, C11 > 0, C44 > 0和C11+2C12 > 0四个条件时, 合金体系才可以稳定存在, 满足力学性能稳定性条件[15]. 立方晶系剪切模量C'可由C11和C12计算:
C′=C11−C122. (1) 材料的力学性质在微观上可以通过弹性常数、体模量、杨氏模量、剪切模量、泊松比和体模量/剪切模量之比等体现. 多晶弹性模量通常可由单晶弹性常数计算获得, 计算方法可用Voigt模型和Reuss模型描述. 经验证, Voigt模型和Reuss模型所得模量通常为多晶模量的上限和下限, 因此Hill[16,17]对两者进行数学平均, 发现所得值与实测值一致性更好. 因此, 实际应用中Hill模型更为常用. 因此本文中Mg2Si体模量BH, 剪切模量GH和杨氏模量E的计算如下:
BV=BR=13(C11+2C12), (2) GV=15(C11−C12+C44), (3) GR=5C44(C11−C12)4C44+3(C11−C12), (4) 其中BV和GV分别为Voigt模型的体积模量和剪切模量, BR和GR分别为 Reuss模型的体积模量和剪切模量. Hill模型中体积模量B和剪切模量G分别为:
BH=BV+BR2, (5) GH=GV+GR2. (6) E=9BHGH3BH+GH, (7) ν=3BH−2GH2(3BH+GH). (8) 通常情况下, B反映材料抗体积压缩能力; G反映材料抗剪切应变的能力, 且与硬度相关[19]; E反映材料抗正应变的能力, 也反映材料的硬度, E越大, 材料的硬度越高[20]. 剪切应力下, ν越小则稳定性越好. B与G的比值是材料工程应用中必须考虑的重要参数, 能反映材料的脆韧性, 以1.75作为临界值. 当 B/G > 1.75时材料呈韧性, 否则呈脆性[20,21].
材料的硬度与剪切模量、体模量之间存在直接的对应关系[22,23]. Miao等 [24]在通过第一性原理研究材料的力学性质中, 提出了硬度H与杨氏模量E和泊松比ν之间计算的经验公式:
H=(1−2ν)E6(1+ν). (9) 3. 计算结果与分析
3.1 晶体弹性常数
Mg2Si-H体系晶体弹性常数C11, C12, C44和C'计算结果如表1所列. 从表1中可以看出, 所有体系中C11, C12和C44均同时满足C11–C12 > 0, C11 > 0, C44 > 0和C11+2C12 > 0四个条件, 说明本文所建立的Mg2Si-H体系可以稳定存在, 满足力学性能稳定性条件[15]. 本工作中的Mg2Si的弹性常数C11, C12和C44计算值分别为113.2, 22.2和45.2 GPa, 与其他学者的理论计算结果[6,20]相比, 最大计算误差分别为1.15%, 2.70%和0.88%, 计算精度较高; 与其他学者的试验结果[21]相比, 本工作中的, C12和C44的最大计算误差分别为10.16%, 14.61和6.80%, 说明本文计算方法具有较高的可靠性.
表 1 Mg2Si-H体系晶体弹常数计算结果Table 1. Elastic constant of Mg2Si-H system.Phase Source Elastic constants/GPa C11 C12 C44 C' C12–C44 Mg2Si This work 113.2 22.2 45.2 45.5 –23.0 Calculated [20] 113.7 22.8 43.5 Calculated [6] 114.5 21.5 45.6 Experimental [21] 126.0 26.0 48.5 Mg2Si-2.05 at.%H This work 110.3 25.1 38.8 42.6 –13.7 Mg2Si-4 at.%H This work 109.0 26.1 37.8 41.4 –11.7 Mg2Si-7.69 at.%H This work 90.6 35.5 34.9 27.5 0.6 图2给出了弹性常数C11, C12, C44和C'随氢原子浓度增加的变化趋势. 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单晶体的抗剪切变形能力. 此外, C12–C44值还可表征晶体的脆性, 当C12–C44 > 0时晶体为塑性, 当C12–C44 < 0晶体为脆性[17]. 从表1可以看出, 当H原子浓度达到7.69 at.%时, C12–C44值从由负值增加到正值, Mg2Si晶体发生了由脆性向塑性的转变.
图 2 弹性常数C11, C12, C44和C'随氢原子浓度增加的变化趋势Fig. 2. Relationships between C11, C12, C44 and C' as a function of H atoms content.3.2 多晶体弹性模量
经(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]相比, 本工作中的BH和E的最大计算误差分别为12.38%和13.63%, 与试验结果同样十分接近, 说明本文计算方法具有较高的可靠性.
表 2 Mg2Si-H体系模量Table 2. Modulus of Mg2Si-H system.Phase Source Modulus GH/GPa BH/GPa E/GPa ν BH/GH H/GPa Mg2Si This work 45.3 52.5 105.6 0.16 1.16 10.31 Calculated [20] 46.2 52.5 107.1 0.16 — — Calculated [6] 44.3 53.4 104.1 0.17 — — Experimental [21] — 59.0 120.0 — — — Mg2Si-2.05 at.%H This work 40.2 53.4 96.5 0.17 1.32 9.07 Mg2Si-4 at.%H This work 39.2 53.7 94.7 0.18 1.37 8.56 Mg2Si-7.69 at.%H This work 31.7 53.8 79.5 0.28 1.69 4.55 图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原子浓度的变化关系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.3.3 拉伸性能
图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]方向进行拉伸的应力-应变曲线\r\nFig. 4. Stress vs. strain curves of 1 × 1 × 2 super cell with 4 at.% and without H during stretching along [001] direction](/fileWLXB/journal/article/wlxb/2019/17/PIC/17-20190368-4.png)
图 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] direction3.4 态密度分析
为探索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 总态密度和分波态密度Fig. 5. DOS and PDOS of (a) Mg2Si and (b) Mg2Si-7.69 at.% H.4. 结 论
1) Mg2Si-H体系内H原子浓度从0增加到7.69 at.%, 剪切模量和弹性模量分别下降24.7%和31.7%, 从而造成体系强度硬度降低, 韧性提高.
2)态密度分析表明, Mg2Si晶体中引入氢原子, 将引起Mg2Si晶体由半导体性质向金属性质的转变, 从而造成体系强度、硬度的降低, 韧性增加.
[1] 李建国, 谭红艳, 史子木, 何迁 2008 中国有色金属学报 18 1819
Google Scholar
Li J G, Tan H Y, Shi Z M, He Q 2008 Chin. J. Nonferrous. Met. 18 1819
Google Scholar
[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
Google Scholar
[3] Qin Q D, Li W X, Zhao K W, Qiu S L, Zhao Y G 2010 Mat. Sci. Eng. A 527 2253
Google Scholar
[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
Google Scholar
[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
Google Scholar
Yu B H, Liu M L, Chen D 2011 Acta Phys. Sin. 60 087105
Google Scholar
[7] 韩秀丽, 王清, 孙东立, 张红星 2008 中国有色金属学报 18 523
Google Scholar
Han S L, Wang Q, Sun L D, Zhang H X 2008 Chin. J. Nonferrous. Met. 18 523
Google Scholar
[8] 张凤春, 李春福, 文平, 罗强, 冉曾令 2014 物理学报 66 227101
Google Scholar
Zhang F C, Li C F, Wen P, Rang Z L 2014 Acta Phys. Sin. 66 227101
Google Scholar
[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
Google Scholar
Yao B D, Hu G Q, Yu Z S 2016 Acta Phys. Sin. 65 026202
Google Scholar
[11] 韩秀丽 2010 博士学位论文 (哈尔滨: 哈尔滨工业大学)
Han X L 2010 Ph. D. Dissertation (Harbin: Harbin Institute of Technology) (in Chinese)
[12] 饶建平, 欧阳楚英, 雷敏生 2012 物理学报 61 047105
Google Scholar
Rao J P, Ouyang C Y, Lei M S 2012 Acta Phys. Sin. 61 047105
Google Scholar
[13] 石瑜, 白洋, 莫丽玢, 向青云, 黄亚丽, 曹江利 2015 物理学报 64 116301
Google Scholar
Shi Y, Bai Y, Mo L B, Xiang Q Y, Huang Y L, Cao J L 2015 Acta Phys. Sin. 64 116301
Google Scholar
[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
Google Scholar
[16] Hill R 2002 Proc. Phys. Soc. 65 349
[17] Jang S H, Chichibu S F 2012 J. Appl. Phys. 112 073503
Google Scholar
[18] Yu W Y, Wang N, Xiao X B 2009 Solid State Sci. 11 1400
Google Scholar
[19] Li Y, Gao Y, Xiao B 2010 J. Alloy Compd. 502 28
Google Scholar
[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
Google Scholar
[23] Leger J M, Djemia P, Ganot F 2001 Appl. Phys. Lett. 79 2169
Google Scholar
[24] Miao N, Sa B, Zhou J 2011 Comp. Mater. Sci. 50 1559
Google Scholar
[25] Senkov O N, Dubois M, Jonas J J 1996 Metall. Mater. Trans. A 27 3963
Google Scholar
[26] Maziyar A, Ahad M, Herbert D, Christian G M 2015 Metall. Mater. Trans. B 46 1471
Google Scholar
-
图 1 H原子位于八面体间隙的Mg2Si-H体系超晶胞结构模型
Fig. 1. Super-cell structural model of Mg2Si-H system.
图 2 弹性常数C11, C12, C44和C'随氢原子浓度增加的变化趋势
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.
Phase Source Elastic constants/GPa C11 C12 C44 C' C12–C44 Mg2Si This work 113.2 22.2 45.2 45.5 –23.0 Calculated [20] 113.7 22.8 43.5 Calculated [6] 114.5 21.5 45.6 Experimental [21] 126.0 26.0 48.5 Mg2Si-2.05 at.%H This work 110.3 25.1 38.8 42.6 –13.7 Mg2Si-4 at.%H This work 109.0 26.1 37.8 41.4 –11.7 Mg2Si-7.69 at.%H This work 90.6 35.5 34.9 27.5 0.6 
下载: 导出CSV
表 2 Mg2Si-H体系模量
Table 2. Modulus of Mg2Si-H system.
Phase Source Modulus GH/GPa BH/GPa E/GPa ν BH/GH H/GPa Mg2Si This work 45.3 52.5 105.6 0.16 1.16 10.31 Calculated [20] 46.2 52.5 107.1 0.16 — — Calculated [6] 44.3 53.4 104.1 0.17 — — Experimental [21] — 59.0 120.0 — — — Mg2Si-2.05 at.%H This work 40.2 53.4 96.5 0.17 1.32 9.07 Mg2Si-4 at.%H This work 39.2 53.7 94.7 0.18 1.37 8.56 Mg2Si-7.69 at.%H This work 31.7 53.8 79.5 0.28 1.69 4.55
下载: 导出CSV
-
[1] 李建国, 谭红艳, 史子木, 何迁 2008 中国有色金属学报 18 1819
Google Scholar
Li J G, Tan H Y, Shi Z M, He Q 2008 Chin. J. Nonferrous. Met. 18 1819
Google Scholar
[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
Google Scholar
[3] Qin Q D, Li W X, Zhao K W, Qiu S L, Zhao Y G 2010 Mat. Sci. Eng. A 527 2253
Google Scholar
[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
Google Scholar
[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
Google Scholar
Yu B H, Liu M L, Chen D 2011 Acta Phys. Sin. 60 087105
Google Scholar
[7] 韩秀丽, 王清, 孙东立, 张红星 2008 中国有色金属学报 18 523
Google Scholar
Han S L, Wang Q, Sun L D, Zhang H X 2008 Chin. J. Nonferrous. Met. 18 523
Google Scholar
[8] 张凤春, 李春福, 文平, 罗强, 冉曾令 2014 物理学报 66 227101
Google Scholar
Zhang F C, Li C F, Wen P, Rang Z L 2014 Acta Phys. Sin. 66 227101
Google Scholar
[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
Google Scholar
Yao B D, Hu G Q, Yu Z S 2016 Acta Phys. Sin. 65 026202
Google Scholar
[11] 韩秀丽 2010 博士学位论文 (哈尔滨: 哈尔滨工业大学)
Han X L 2010 Ph. D. Dissertation (Harbin: Harbin Institute of Technology) (in Chinese)
[12] 饶建平, 欧阳楚英, 雷敏生 2012 物理学报 61 047105
Google Scholar
Rao J P, Ouyang C Y, Lei M S 2012 Acta Phys. Sin. 61 047105
Google Scholar
[13] 石瑜, 白洋, 莫丽玢, 向青云, 黄亚丽, 曹江利 2015 物理学报 64 116301
Google Scholar
Shi Y, Bai Y, Mo L B, Xiang Q Y, Huang Y L, Cao J L 2015 Acta Phys. Sin. 64 116301
Google Scholar
[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
Google Scholar
[16] Hill R 2002 Proc. Phys. Soc. 65 349
[17] Jang S H, Chichibu S F 2012 J. Appl. Phys. 112 073503
Google Scholar
[18] Yu W Y, Wang N, Xiao X B 2009 Solid State Sci. 11 1400
Google Scholar
[19] Li Y, Gao Y, Xiao B 2010 J. Alloy Compd. 502 28
Google Scholar
[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
Google Scholar
[23] Leger J M, Djemia P, Ganot F 2001 Appl. Phys. Lett. 79 2169
Google Scholar
[24] Miao N, Sa B, Zhou J 2011 Comp. Mater. Sci. 50 1559
Google Scholar
[25] Senkov O N, Dubois M, Jonas J J 1996 Metall. Mater. Trans. A 27 3963
Google Scholar
[26] Maziyar A, Ahad M, Herbert D, Christian G M 2015 Metall. Mater. Trans. B 46 1471
Google Scholar
期刊类型引用(2)
1. 陆贺港，金玉，黄海东. 三赞胶-乙基纤维素基包膜缓释肥的制备及特性研究. 西北农业学报. 2022(09): 1194-1201 . 
百度学术2. 宋莉莉，傅高升，李龙泽，王火生. θ-Al_2Cu金属间化合物中氢原子占位的第一性原理研究. 机械工程材料. 2021(12): 67-71+77 . 百度学术其他类型引用(5)
下载:
下载:
计量
- 文章访问数: 10693
- PDF下载量: 103
- 被引次数: 7
