王浩玉 农智升† 王继杰 朱景川
1) (沈阳航空航天大学材料科学与工程学院, 沈阳 110136)
2) (哈尔滨工业大学材料科学与工程学院, 哈尔滨 150001)
(2018 年10 月23日收到; 2018 年12 月4日收到修改稿)
为了探索AlxCrFeNiTi系高熵合金组成成分和弹性性质的关系, 结合固溶体特征参数和第一性原理计算, 研究Al元素含量对AlxCrFeNiTi (x = 0, 0.5, 1, 2, 3, 4)合金结构和弹性性质的影响, 并分析合金固溶体特征参数与弹性性质之间的关系. 结果表明: AlxCrFeNiTi系合金的价电子浓度随着Al含量的增加逐渐减小,合金在体心立方结构下的形成焓均低于面心立方结构, 说明研究的AlxCrFeNiTi系合金会形成单一的体心立方结构固溶体; 合金的晶格常数和形成能力强弱随着Al含量的增加而增大, 但合金的结构稳定性略有下降;当合金元素按照等原子比进行成分配比时, 合金的原子尺寸差异最大; AlxCrFeNiTi系合金中不同原子之间除了金属键结合外, 还表现出一定的共价和离子键结合特征; 对于AlxCrFeNiTi系合金而言, 随着热力学熵焓比的增大, 合金体弹模量和韧性随之增大; 随着合金混合焓的增加, 合金在压缩方向的各向异性程度明显降低. 热力学熵焓比和混合焓可作为AlxCrFeNiTi系高熵合金成分设计的重要参数.
一般地, 将五种及五种以上金属元素以等摩尔比或近摩尔比混合得到的合金称为高熵合金[1−3].高熵合金通常为体心立方(body centered cubic,BCC)或面心立方(face centered cubic, FCC)结构的固溶体, 一般认为高熵效应是形成这种独特固溶体结构的主要原因[2]. 基于这样的固溶体结构,高熵合金具有传统合金不能比拟的高硬度、高强度、耐腐蚀性、耐磨性和电磁性等综合性质[4,5]. 目前对于高熵合金的研究多数集中在含AlCrFeNi元素的体系, 并通过添加不同的合金元素对合金的性能进行调控, 以获得具有不同特性的高熵合金体系, 如高耐磨的AlCo1.5CrFeNi1.5Ti合金[6]、高强的 AlCoCrFeNiV合金[7]以及抗高温氧化的AlCrFeMn0.5NiTi合金[8]体系等. 研究表明, 高熵合金的结构和性能与合金成分有着密切的联系. 刘瑞文[9]研究了Al元素对AlxCoCrFeNiTi0.5高熵合金结构和性能的影响, 发现随着Al元素的增加,合金的相结构由FCC固溶体逐渐向BCC固溶体转变, 并且合金的硬度逐渐增加.
对高熵合金成分的表征一般通过合金的固溶体特征参数进行描述, 并结合统计方法进行结构预测. 原子尺寸差δ和热力学熵焓比Ω通常被用于预测多组元合金分别形成无序固溶体、有序固溶体、金属间化合物和块状非晶时所处的参数区域[10]. 研究发现当合金的固溶体特征参数在δ≤ 6.6%和Ω≥ 1.1时最易形成稳定的固溶体. 同时, 多主元合金的价电子浓度(valance electron concentration,VEC)在预测高熵合金形成单一BCC、单一FCC或BCC+FCC的混合结构上起到很好的效果. 高熵合金的结构稳定性和弹性性质通常采用基于密度泛函理论(density functional theory, DFT)的第一性原理计算进行研究. 对单相TiZrNbMoVx(x= 0—1.5)高熵合金的弹性性质的第一性原理计算表明[11], 当合金的VEC在4.72左右时, 该合金具有最优弹性各向异性性质. 尽管在理论上已经有这些对高熵合金结构和弹性性质的研究, 但成分和性能之间的直接关系对于高熵合金体系设计而言仍是需关注的重点.
为了探索高熵合金的组成成分和弹性性质之间的关系, 本文选择AlxCrFeNiTi (x= 0, 0.5, 1,2, 3, 4)系列单相的合金体系进行研究. 该新型无Co的合金体系既包含了高熵合金最常采用的合金元素, 又满足了降低高熵合金成本的设计原则. 合金的组成成分采用固溶体特征参数进行表征, 理论的弹性性质通过第一性原理计算获得. 本文首先研究Al含量变化对合金结构和弹性性质的影响, 最后深入探讨固溶体特征参数与弹性性质之间的关系, 为更好地进行高熵合金体系设计提供理论参考.
为了表征不同Al元素含量下AlxCrFeNiTi (x=0, 0.5, 1, 2, 3, 4)合金的组成成分, 引入原子尺寸差异δ、混合焓 ∆Hmix、VEC 和热力学熵焓比Ω这四个固溶体特征参数. 其中, 原子尺寸差异δ与置换固溶体的形成有关, 计算式为[12]
式中,ci是第i个元素的原子百分比,ri是该原子的半径. 通常来说, 稳定的置换固溶体的原子尺寸差异δ在1—6之间. 混合焓 ∆Hmix是另一个影响固溶体形成的重要参数, 当混合焓的值在-0.17 ≤∆Hmix≤ 0.1 eV区间时[10], 高熵合金将形成稳定的固溶体结构. 混合焓的计算式为[13]
其中(VEC)i表示第i个元素的VEC. 热力学熵焓比Ω是一个结合热力学熵和焓的复合固溶体特征参数, 其表达式为[10]
式中, ∆Sconf为固溶体形成时的混合熵,Tm为合金的熔点, ∆Hmix是由(2)式所计算的合金的混合焓.本文对AlxCrFeNiTi系合金固溶体参数的计算所用到的相关元素的物理参数如表1所列.
表1 AlxCrFeNiTi系合金中元素的物理参数[15]Table 1. Physical parameters of alloying elements in AlxCrFe NiTi alloys[15].
为了研究AlxCrFeNiTi (x= 0, 0.5, 1, 2, 3, 4)高熵合金的结构稳定性和弹性性质, 本文采用基于DFT的第一性原理对AlxCrFeNiTi系高熵合金进行研究. AlxCrFeNiTi系合金的晶体结构在基于BCC和FCC结构的超晶胞上建立, 其中所有原子进行随机占位, 并且满足最邻近原子不能为同类型原子的原则[8]. 合金元素的电子和离子之间的相互作用采用超软贋势进行描述, 电子交换关联函数采用广义梯度近似中的Perdew−Burke−Ernzerhof形式[16]. 对建立的BCC结构的AlxCrFeNiTi (x=0, 0.5, 1, 2, 3, 4)合金, DFT计算中所使用的k点分别设置为 9 × 9 × 4, 8 × 8 × 1, 9 × 9 × 2, 9 ×9 × 3, 9 × 9 × 1 和 8 × 8 × 2; 对建立的 FCC 结构的AlxCrFeNiTi合金,k点分别设置为7 × 6 × 7,6 × 8 × 1, 7 × 7 × 1, 6 × 6 × 3, 7 × 6 × 1 和6 × 8 × 3. 为了兼顾精度和计算速度, 将平面波截断能设置为400 eV. 第一性原理自洽计算的收敛条件为: 每个原子上的力低于0.01 eV/Å, 公差偏移小于 5.0 × 10−4Å, 最大压力小于 0.02 GPa 和总能量小于5.0 × 10−6eV/atom.
材料的多晶体弹性性质, 包括体弹模量B、剪切模量G、杨氏模量E和泊松比ν, 可以根据Voigt−Reuss−Hill (V−R−H)近似结合第一性原理计算获得的弹性常数进行预测, 计算公式分别为[17]
其中多晶体体弹模量B和剪切模量G分别是根据Voigt和Reuss模型获得的模量的算数平均值.体弹模量B和剪切模量G的计算公式分别为[17]
其中Cij代表材料的单晶弹性常数;Sij为弹性柔度, 是弹性矩阵Cij的逆矩阵.
材料的弹性各向异性是一个重要的弹性性质.基于Voigt和Reuss模型获得的弹性模量, 在压缩方向上的各向异性AB和剪切方向上各向异性的AG通常被应用于研究材料的弹性各向异性, 计算式分别为[18]
AB和AG的值能够描述合金的弹性各向异性性质,值越大表明合金在该方向上的各向异性程度越大.当AB和AG的值为0%时, 结构表现为各向同性特征; 而当AB和AG的值达到100%时就意味着合金的各向异性特征最大.
根据经典的多组元合金的VEC理论[15], 不同的晶体结构有不同的VEC特征, 当VEC< 6.87时形成单一的BCC结构固溶体; 当VEC≥ 8.00时形成单一的FCC结构固溶体; 在6.87 ≤VEC< 8.00时形成BCC和FCC结构的混合固溶体.不同Al元素含量下AlxCrFeNiTi (x= 0, 0.5, 1,2, 3, 4)高熵合金的VEC如图1所示. 除了x= 0时的CrFeNiTi四元合金的VEC略高于6.78外,其余AlxCrFeNiTi (x= 0.5, 1, 2, 3, 4)五种合金的VEC均小于6.87, 说明当Al含量的原子百分比为11.11 at.%, 20 at.%, 33.33 at.%, 42.86 at.%和50 at.%时, AlxCrFeNiTi系合金有可能形成单一BCC结构的固溶体. 此外, 从能量角度描述合金形成不同结构的能力强弱可以通过形成焓Hf进行表征,AnBm合金形成焓的计算式为[19]
式中Etot是单相AlxCrFeNiTi系合金经晶格优化后得到的平衡结构的总能量,n和m分别是不同合金元素在结构中对应的原子数目,是X元素在其初始平衡结构下的能量.
图1 Al含量对AlxCrFeNiTi系合金的VEC、形成焓和结构的影响, 其中当x = 0, 0.5, 1, 2, 3, 4时, 分别对应Al含量的原子百分比为0, 11.11 at.%, 20 at.%, 33.33 at.%,42.86 at.%, 50 at.%Fig.1. Effect of Al contents on valence electron concentr−ation, formation enthalpy and structure of AlxCrFeNiTi alloys. When x = 0, 0.5, 1, 2, 3, 4, the atomic percentages corresponding to Al content are 0, 11.11 at.%, 20 at.%,33.33 at.%, 42.86 at.%, 50 at.% respectively.
形成焓的数值代表化合物形成前后的能量差异, 拥有更低负值的形成焓的合金晶体结构更容易形成. 采用 Broyden−Fletcher−Goldfarb−Shanno 方法[20]对建立的BCC和FCC结构的AlxCrFeNiTi(x= 0, 0.5, 1, 2, 3, 4)合金进行晶格结构优化, 获得BCC和FCC结构下AlxCrFeNiTi系合金的平衡结构, 并根据(8)式计算得到不同Al含量下AlxCrFeNiTi系合金形成不同结构的形成焓, 结果如图1所示. 可以看出, 在相同的Al含量下, BCC结构AlxCrFeNiTi合金的形成焓均低于FCC结构下的合金, 这说明从能量角度上分析, 相对于FCC结构, 合金形成BCC结构固溶体所需能量更少, 由此可以推断, BCC结构为AlxCrFeNiTi系合金优先形成的结构. 另外, 随着Al含量的增加,AlxCrFeNiTi系合金的形成焓明显下降, 说明合金的形成能力随着Al含量的增加而增强. 与此同时,由图1可以发现, AlxCrFeNiTi系合金的VEC随着Al含量的增大也逐渐下降. 固溶体的VEC值越低, 说明合金越易于形成低原子堆垛特征的结构, 即BCC结构, 这是因为这种原子堆垛方式的能量更稳定[14]. 这与形成焓的计算结果相符合.
VEC理论和形成焓计算结果对AlxCrFeNiTi(x= 0.5, 1, 2, 3, 4)这五种合金形成结构倾向性的结论完全符合, 但对于四元CrFeNiTi合金, 由于其区别于传统意义上的五元及以上的高熵合金, 虽然其VEC略高于判据中的6.87, 但根据有关四元等摩尔比合金的文献报道[21]以及形成焓计算可知,四元CrFeNiTi合金同样形成单一的BCC结构固溶体. 因此, 结合VEC理论和形成焓计算结果, 可以预测AlxCrFeNiTi (x= 0, 0.5, 1, 2, 3, 4)合金的晶体结构均为单相的BCC结构. 本文其余有关第一性原理的计算, 包括晶格常数、电子结构和弹性性质等, 对AlxCrFeNiTi系合金的研究均基于BCC结构.
图2为Al含量对AlxCrFeNiTi系合金平衡结构下晶格常数和原子尺寸差异的影响. 由于所研究的AlxCrFeNiTi系合金的晶体结构是基于BCC结构上建立的超晶胞, 在计算充分弛豫后, 组成超晶胞的各个独立的单胞具有相似的晶格常数a,b和c. 为了便于与实验数据进行对比, 采用晶格常数的算术平均数aavg对AlxCrFeNiTi系合金的晶格参数进行描述, 计算公式为[8]
图2 Al含量对AlxCrFeNiTi系合金的晶格常数和原子尺寸差异的影响Fig.2. Effect of Al contents on lattice constant and atomic size difference of AlxCrFeNiTi alloys.
式中ai,bi,ci, 分别表示结构优化后超晶胞中第i个单胞的晶格常数;n是超晶胞中单胞的数目. 根据含AlCrFeNi元素高熵合金的相关文献[1−5], 形成稳定BCC结构时合金的晶格常数约为2.92 Å.由图2的晶格常数曲线可知, 当Al元素的原子百分比为0, 11.11%和20%时, AlxCrFeNiTi (x= 0,0.5, 1)合金的晶格常数在2.92 Å附近, 形成的BCC结构固溶体更为稳定; 当Al元素的原子百分比为33.33%, 42.86%和50%时, AlxCrFeNiTi (x=2, 3, 4)合金的晶格常数在2.98—3.00 Å之间. 这是由于Al元素的原子半径较大, Al元素在合金体系中的比重增大使得AlxCrFeNiTi系合金的晶格参数明显增大.
从图2中原子尺寸差异δ的曲线可以看出, 虽然合金的晶格常数随着Al含量增大逐渐增加, 但表征晶格畸变程度的原子尺寸差异并非持续增大.当Al元素原子百分比含量低于20 at.%时, 随着大原子半径的Al元素的增加, 合金的原子尺寸差异逐渐增大; 当Al元素的含量超过20 at.%时, 此时随着Al含量的增加, 合金的原子尺寸差异却逐渐减小; 当AlxCrFeNiTi系合金中各元素的配比为等摩尔比时, 即合金为AlCrFeNiTi时, 合金的原子尺寸差异到达最大. 这是由于初始较大原子半径的Al原子的加入提高了CrFeNiTi四元合金的晶格畸变程度, 但随着Al元素含量在合金中的比重增加, 使得合金趋向于形成以Al为基体的固溶体,因此转而使合金的原子尺寸差异即晶格畸变降低.
图3 CrFeNiTi, Al2CrFeNiTi和Al4CrFeNiTi合金的总态密度Fig.3. Density of states (DOS) of CrFeNiTi, Al2CrFeNiTi and Al4CrFeNiTi alloys.
为了研究AlxCrFeNiTi系合金的结构稳定性来源, 计算了CrFeNiTi合金、Al2CrFeNiTi合金和Al4CrFeNiTi合金在平衡结构下的总态密度,结果如图3所示. 费米能级处(0 eV)较高的态密度值说明了这三种合金均表现出明显的金属特征.随着合金体系中Al元素含量的不断增加, 合金的总态密度在低能级处的峰位(−72, −57, −44和−24 eV附近)向高能级方向移动, 同时, 总态密度在高能级处(费米能级附近)呈现同样的规律, 说明随着Al含量增加, AlxCrFeNiTi系合金晶体结构的稳定性有所下降.
除了分析总态密度的变化趋势之外, 合金的结构稳定性还可以通过计算合金在平衡结构下的结合能Ec进行定性描述. 结合能Ec反映了晶体被分解为单个原子所需要的能量,AnBm合金结合能的计算公式为[19]
其中Etot是合金在平衡结构下的总能量,n和m分别是不同合金元素在结构中对应的原子数目,是孤立X原子的能量. 更低的结合能Ec说明分解该晶体结构需要更大的能量, 此时该晶体结构的稳定性更高. 计算获得的AlxCrFeNiTi (x= 0,0.5, 1, 2, 3, 4)合金在平衡结构下Ec分别为−5.64,说明随着Al含量增加, AlxCrFeNiTi合金的结合能不断升高, 分解合金结构所需能量减小, 合金的结构稳定性逐渐下降. 这与态密度的分析结果完全一致.
图4 AlxCrFeNiTi系合金的电荷分布 (a) Al2CrFeNiTi在(110)面的差分电荷密度; (b) AlxCrFeNiTi系合金的原子布居数Fig.4. Charge distribution of AlxCrFeNiTi alloys: (a) Elec−tron density difference map on (110) plane for Al2CrFe NiTi; (b) the atomic Mulliken population of AlxCrFeNiTi alloys.
差分电荷密度和原子布居数能够描述合金键合之后原子电荷分布的变化情况, 可以用来研究不同元素对合金中成键的贡献. 以Al2CrFeNiTi合金作为例子, 图4(a)为选取的包含所有组成原子的部分(110)截面上的差分电荷密度分布图. 原子之间电荷密度的明显增加说明成键后合金内形成离域的价电子云, 这是原子间形成金属键结合的显著特点. 此外, Fe和Ni原子在得电子区域以及Al和Ti原子在失电子区域均表现出一定的方向性, 说明这些原子的键合可能具有离子结合的特征;而Cr−Cr原子之间存在一定共价键特征的差分电荷密度分布. 图4(b)为BCC结构的AlxCrFeNiTi系合金在平衡结构下不同元素的原子布居数分析结果. 由于在计算过程中同种元素的每一个原子都具有极其相似的得失电子的能力, 因此在图4(b)中不同元素的原子布居数是该种元素所有原子的原子布居数的平均值. 当原子布居数的平均值为正时, 意味着该元素在键合之后失去电子, 而布居数为负值则意味着得到电子; 更小的负值(或更大的正值)说明该元素具有更强的得(失)电子的能力[22].由图4(b)可知, Al原子和Ti原子在成键后均是失去电子, 而Ni原子和Fe原子均得到电子. Cr原子除了在四元CrFeNiTi合金时失去电子, 在其余AlxCrFeNiTi (x= 0.5, 1, 2, 3, 4)合金中均得到电子. 价电子主要从Al原子和Ti原子的s轨道上转移到了Cr原子、Ni原子和Fe原子上, 并且Ti原子展现出比Al原子更强的失电子的能力. 随Al元素含量的增加, Cr和Ni原子得电子能力呈现出增强的趋势, 而Fe原子基本保持不变. 另外, 除四元CrFeNiTi合金以外, 其余合金中Fe元素与Ti元素的布居数随Al元素含量变化的曲线上呈一定的对称分布特征, 可以推测Ti原子失去的电子有可能主要迁移到Fe原子上. 在合金中存在元素电荷迁移现象意味着有不同元素原子之间存在一定的离子键作用力, 这有利于合金形成稳定的结构. 结合差分电荷密度和原子布居数分析, 在AlxCrFeNiTi系合金中原子之间的键和作用既有金属键结合作用又存在一定的共价和离子键结合作用, 这样的多种键合相互作用, 使得合金具有较好的结构稳定性.
材料的弹性常数是直观表述材料弹性性质的重要参数, 反映了晶体结构的力学稳定性, 对研究合金的弹性性质如弹性模量和弹性各向异性等具有重要意义. 表2为通过第一性原理计算获得的不同Al含量下AlxCrFeNiTi系合金的弹性常数Cij.
合金的机械稳定性可通过弹性常数进行判断.对于本文研究的基于BCC结构下的AlxCrFeNiTi系合金, 基态下(0 K, 0 GPa)机械稳定性判据为[23]
结合表2中获得的AlxCrFeNiTi系合金的弹性常数, 可知所有数据均满足弹性稳定性判据, 说明AlxCrFeNiTi系合金在基态下均具有良好的机械稳定性, 同样说明了AlxCrFeNiTi系合金具有稳定的晶体结构.
为了预测AlxCrFeNiTi系合金的弹性性质, 采用(5)和(6)式中的 V−R−H近似方法, 计算得到了多晶体AlxCrFeNiTi系合金的弹性常数, 结果如图5所示. 体弹模量B反映晶体在外部压力下抵抗体积变形的能力, 剪切模量G是晶体在剪切应力下抵抗可逆变形的度量, 杨氏模量E是表征晶体在弹性极限内抗压或抗拉的物理量[23]. 从图5可以看出, 合金的体弹模量B分布在130—200 GPa之间, 剪切模量G大小从50到80 GPa, 杨氏模量E在130—240 GPa之间. 随着AlxCrFeNiTi系合金中Al含量的增多, 合金的体弹模量B整体呈下降趋势, CrFeNiTi的体弹模量B最大; 剪切模量G整体变化不大, Al4CrFeNiTi的剪切模量G最大; 杨氏模量E呈先升高后下降最后上升的趋势, Al4CrFeNiTi的杨氏模量E最大. 因此,Al含量的增加使得AlxCrFeNiTi合金抵抗体积变形的能力逐渐减弱, 而抵抗可逆变形以及在弹性范围内抗压的能力有所提高.
表2 不同Al含量下AlxCrFeNiTi系合金的弹性常数CijTable 2. Elastic constant Cij (in GPa) of AlxCrFeNiTi alloys with different Al contents.
图5 Al含量对AlxCrFeNiTi系合金弹性模量的影响Fig.5. The effect of Al contents on elastic constants AlxCrFeNiTi alloys.
脆韧性是材料一个非常重要的性质, 用来描述材料在断裂前吸收能量和进行塑性变形的能力, 通常综合使用体弹模量和剪切模量的比值B/G、泊松比ν和(C12-C44)这三个参数对合金的脆韧性进行预测. 当合金同时满足B/G> 1.75,ν> 0.26和(C12-C44) > 0时, 合金表现出韧性特征, 否则呈脆性[24]. 图6为Al含量对AlxCrFeNiTi系合金的B/G、泊松比ν和(C12-C44)的影响. 除了Al4CrFeNiTi呈现脆性(其中,B/G大约为1.5, 泊松比ν约为2.3, (C12-C44)为−18 GPa)外, 其余的CrFeNiTi, Al0.5CrFeNiTi, Al1CrFeNiTi, Al2Cr FeNiTi和Al3CrFeNiTi合金均表现明显的韧性特征. 并且随着Al含量的增加, 这三个参数均呈缓慢下降趋势, 其中B/G曲线和泊松比ν曲线变化一致. 在AlxCrFeNiTi系合金中CrFeNiTi合金的B/G、泊松比ν和 (C12-C44)均为最大, 说明CrFeNiTi合金的韧性最好, AlxCrFeNiTi系合金的韧性随着Al含量的增加而降低, 而当Al的原子百分比含量增大至50%时, 合金由韧性转变为脆性.
图6 Al含量对AlxCrFeNiTi系合金的B/G、泊松比ν和 (C12− C44)的影响Fig.6. Effect of Al contents on B/G, Poisson’s ratioν and (C12−C44 ) of AlxCrFeNiTi alloys.
高熵合金组成成分和弹性性质之间的关系是高熵合金体系成分设计所关注的重点, 结合研究获得的Al含量对AlxCrFeNiTi系合金弹性性质的影响, 通过固溶体特征参数对合金的组成成分进行表征, 并作为横坐标, 以计算获得的合金弹性性质作为纵坐标, 探索AlxCrFeNiTi系合金成分和弹性性质的关系. 图7为AlxCrFeNiTi系合金固溶体特征参数和弹性性质之间的关系. 从图7(a)可以看出,剪切模量G、杨氏模量E与热力学熵焓比Ω虽然没有明显的关系, 但体弹模量B和B/G都随着热力学熵焓比Ω的升高呈上升趋势. 这说明通过调整元素成分改变AlxCrFeNiTi系合金的热力学熵焓比Ω, 可以控制合金抵抗压缩时体积变形和韧脆性的能力. 当提高合金的热力学熵焓比Ω时, 合金将获得更好的抵抗压缩时体积变形和韧性的能力. 对于合金的弹性各向异性而言, 仅存在与混合焓∆Hmix有明显的关系, 如图7(b)所示. 在压缩方向上的各向异性程度AB随着混合焓 ∆Hmix的增加而下降, 这说明若想在压缩方向上获得各向同性特性较好的AlxCrFeNiTi系合金, 可以通过改变成分提高合金的混合焓 ∆Hmix来实现.
为了进一步验证AlxCrFeNiTi系合金弹性各向异性的性质, 根据图7(b)选取了弹性各向异性较小的CrFeNiTi和各向异性较大的Al2CrFeNiTi的两个合金, 对其杨氏模量E在三维空间里各向异性的曲面表示进行研究, 结果如图8所示. 由于建立的AlxCrFeNiTi系合金的晶体结构属于四方晶系晶体结构, 相应的杨氏模量E在三维空间里各向异性的曲面可以根据(12)式进行计算[25]:
其中l1,l2和l3分别是X,Y和Z轴方向上的余弦;Sij是弹性柔度. 在三维空间上的图形越接近球形说明所研究的结构越趋近于各向同性特征, 反之表现为各向异性特征. 从图8可以看出, CrFeNiTi,Al2CrFeNiTi合金的弹性模量三维曲面图均不是球形, 说明这两个合金都表现出一定的各向异性特征, 但CrFeNiTi合金的三维曲面图相对于Al2CrFeNiTi合金的三维曲面图更接近球形特征,表明Al2CrFeNiTi合金显示了更为明显的弹性各向异性特征. 该结论与图7(b)中合金的弹性各向异性大小结论相符合. 尽管可以通过热力学熵焓比Ω和混合焓 ∆Hmix对AlxCrFeNiTi系合金的弹性性质进行设计, 但仍需要在理论和实验方面进行进一步的探索, 为高熵合金的成分设计提供参考.
图7 AlxCrFeNiTi系合金固溶体特征参数和弹性性质之间的关系 (a) 热力学熵焓比 Ω 与B, G, E, B/G的关系; (b) 混合焓ΔHmix与AG, AB的关系Fig.7. Correlativity between solid solution characteristics and elastic properties of AlxCrFeNiTi alloys: (a) Ratio of thermodynamic entropy to enthalpy and B, G, E, B/G; (b) formation entahlpy Δ Hmix and AG, AB.
图8 合金杨氏模量E的各向异性 (a) CrFeNiTi; (b) Al2CrFeNiTiFig.8. Anisotropy of Young modulus E of CrFeNiTi, Al2CrFeNiTi alloys: (a) CrFeNiTi; (b) Al2CrFeNiTi.
1) 研究的AlxCrFeNiTi (x= 0.5, 1, 2, 3, 4)系合金的VEC均低于6.87, 且BCC结构下的形成焓均要低于FCC结构下的形成焓, 合金将形成单一的BCC结构固溶体. 随着Al元素含量的增加, AlxCrFeNiTi系合金的VEC逐渐降低, 但合金的形成能力不断增强; 晶格参数随着Al含量增加不断上升, 合金的原子尺寸差异(反映合金晶格畸变的程度)先上升后下降.
2) AlxCrFeNiTi系合金的总态密度随着Al含量的增加往高能级方向移动, Al和Ti在成键后均是失去电子, 而Ni和Fe均得到电子. AlxCrFeNiTi系合金的原子之间形成复合的键和作用.
3) 随着Al含量的增加, AlxCrFeNiTi系合金的体弹模量B呈下降趋势而剪切模量G整体变化不大, 韧性则逐渐降低. 其中Al4CrFeNiTi合金呈明显的脆性特征.
4) 热 力 学 熵 焓 比Ω和 混 合 焓 ∆Hmix可 作 为AlxCrFeNiTi系合金弹性性质预测的参数. 随着热力学熵焓比Ω的升高, 合金的体弹模量B和B/G明显升高; 随着混合焓 ∆Hmix的增加, 合金在压缩方向的各向异性程度AB逐渐下降.