Processing math: 100%

Sc、Ti、V、Zr掺杂Cr2B3高压力学性质的第一性原理计算

张畅 孙小伟 宋婷 田俊红 刘子江

刘雨佳, 侯海量, 李茂, 金键, 戴文喜. 前舱物对低速大质量平头弹侵彻金属薄板的影响[J]. 高压物理学报, 2020, 34(1): 015104. doi: 10.11858/gywlxb.20190830
引用本文: 张畅, 孙小伟, 宋婷, 田俊红, 刘子江. Sc、Ti、V、Zr掺杂Cr2B3高压力学性质的第一性原理计算[J]. 高压物理学报, 2022, 36(4): 042201. doi: 10.11858/gywlxb.20210916
LIU Yujia, HOU Hailiang, LI Mao, JIN Jian, DAI Wenxi. Influence of Nose Cabin on Low Speed Blunt Projectile during Penetration of Metal Plate[J]. Chinese Journal of High Pressure Physics, 2020, 34(1): 015104. doi: 10.11858/gywlxb.20190830
Citation: ZHANG Chang, SUN Xiaowei, SONG Ting, TIAN Junhong, LIU Zijiang. First-Principles Study on Mechanical Properties of Sc, Ti, V, Zr-Doped Cr2B3 at High Pressure[J]. Chinese Journal of High Pressure Physics, 2022, 36(4): 042201. doi: 10.11858/gywlxb.20210916

Sc、Ti、V、Zr掺杂Cr2B3高压力学性质的第一性原理计算

doi: 10.11858/gywlxb.20210916
基金项目: 甘肃省高等学校创新基金(2020A-039);甘肃省自然科学基金重点项目(20JR5RA427);甘肃省杰出青年科学基金(145RJDA323)
详细信息
    作者简介:

    张 畅(1996-),女,硕士研究生,主要从事功能新材料高压物性研究.E-mail:zhangch_lzjtu@126.com

    通讯作者:

    孙小伟(1979-),男,博士,教授,主要从事材料高压物性、声子晶体和声学超材料研究.E-mail:sunxw_lzjtu@yeah.net

  • 中图分类号: O521.2; O482.1

First-Principles Study on Mechanical Properties of Sc, Ti, V, Zr-Doped Cr2B3 at High Pressure

  • 摘要: 采用基于密度泛函理论的第一性原理方法,计算了掺杂Sc、Ti、V和Zr的Cr2B3在零压下的晶体结构和电子结构及其在0~150 GPa压力范围内的弹性常数和维氏硬度。结果表明:Cr2B3及其掺杂化合物均具备力学稳定性;在零压下,添加Sc、Ti、V和Zr元素均可提高Cr2B3的维氏硬度,其中Ti掺杂Cr2B3的硬度由26.3 GPa提高至40.2 GPa,提高52.9%,达到超硬材料标准,且 Ti和V掺杂Cr2B3的剪切模量分别提高14.3%和16.2%,杨氏模量分别提高8.2%和12.0%;由电子结构分析可知,Sc、Ti、V和 Zr元素可以加强B与B之间的电子局域化程度,从而增强共价键结合强度,使Cr2B3的硬度升高;Cr2B3的弹性常数、体积弹性模量、剪切模量、杨氏模量以及硬度随着压力的增加而增加,但其硬度仍较低,150 GPa下仅为28.3 GPa,而掺杂V的Cr2B3的硬度在整个压力范围内约为37 GPa。研究结果可为Cr2B3在高压等特殊条件下的应用提供理论参考。

     

  • 半穿甲导弹是在穿甲弹基础上发展起来的用于攻击舰船目标的武器,兼具穿甲和爆破毁伤能力,是舰船防护领域的重点设防对象[1]。前舱物作为半穿甲导弹的重要组成部分,在导弹的制导和控制方面起着至关重要的作用。半穿甲导弹侵彻靶板的过程中,前舱物先于战斗部接触靶板,对靶板产生影响。因此,研究前舱物对战斗部侵彻靶板的影响具有实际意义。迄今为止,对于半穿甲导弹侵彻舰船结构等问题,研究人员通常忽略前舱物的影响,即只考虑战斗部的动能穿甲作用,将其抽象为低速大质量弹体的穿甲问题。陈斌等[2]利用忽略前舱物的半穿甲模拟弹,开展了侵彻陶瓷/钢复合装甲试验,明确了弹着角对半穿甲弹侵彻效能有明显影响;朱锡等[3]采用与实际半穿甲导弹战斗部几何形状相似的弹丸,开展了侵彻多种装甲结构试验,得到了不同装甲结构形式抗半穿甲导弹动能穿甲的特性。目前,针对前舱物对半穿甲导弹侵彻靶板影响的研究相对较少,这是由于前舱物的组成部件复杂,无法完全模拟所有器件,因此,通常对前舱物进行力学等效。

    对前舱物的等效主要有两种方法。第一种方法为采用质量块对前舱物进行等效。楼建锋等[4]通过有限元方法,分别使用金属材料和脆性材料质量块等效模拟前舱物,得到了高速(675 m/s)和低速(250 m/s)两种速度条件下前舱物对战斗部侵彻过程的影响。陈刚[5]利用质量块模拟前舱物,对比分析了前舱物对战斗部高速撞击靶板的影响,得到了前舱物改进战斗部弹体在斜穿甲过程中的荷载环境条件。Chen等[6]考虑前舱物作用导致的靶板结构预响应,将前舱物剩余部分等效为附着在战斗部前部的质量块,并对剪切冲塞模型进行修正,得到了与实验结果相吻合的终点弹道公式。该等效方法的关键是质量块材料的选取,只有选择与前舱物名义力学特性相近的材料,才能正确表征前舱物对导弹弹体动能穿甲的影响。第二种方法为采用薄壁圆筒结构对前舱物进行等效。徐钰巍等[7]采用硬铝圆筒结构对前舱物进行等效,建立了变截面刚度/含缺陷的自由梁模型,并结合模型试验得到了前舱物在弹体侵彻靶板初期对弹体姿态变化的影响规律。该等效方法忽略了前舱物内部仪器设备强度的影响,仅考虑前舱段外壳在侵彻过程中所起的作用,在研究半穿甲导弹正侵彻时,与实际情况的差异较大。对于正侵彻而言,将前舱物等效为质量块的方法是可行的,但现有的等效方法为了体现前舱物稀松多孔、中空易碎的结构特点,所选取的等效材料模型和前舱物在撞击过程中所表现出的力学特性存在不相符的情况。

    本研究采用轻质泡沫铝材料模拟半穿甲导弹前舱物,建立变截面刚度弹体侵彻金属薄板的有限元模型,分析带前舱物平头弹侵彻金属薄板的动态力学过程,研究不同靶板厚度、弹速条件下前舱物对平头战斗部侵彻金属薄板的影响。

    前舱物主要由雷达导引头、导航系统、电源系统、高度表等仪器设备组成,其特点是质量较大(约为战斗部总质量的1/4~1/3)且分布均匀,整体结构中空、易碎、延展性差,结构强度远低于战斗部。前舱物由高分子聚合物、铝合金、单晶硅等多种材料构成,平均密度低(小于1 g/cm3)。前舱物在撞击靶板时,会经历“接触—屈服—压实—破碎飞散”的破坏过程,其中屈服阶段持续时间长,名义应力基本保持不变,即存在较长的平台应力,随后前舱物被压实,名义应力大幅增加,前舱段整体长度大幅减小,这与泡沫铝材料在冲击荷载作用下的位移变化曲线类似[8-10],如图1所示。因此,本研究中前舱物材料模型采用低强度、低密度、有长应力平台的泡沫铝材料进行等效。

    图  1  泡沫铝材料载荷-位移曲线
    Figure  1.  Loading-displacement curve for aluminum foam material

    半穿甲反舰导弹通常由前舱物、战斗部和推进系统3部分组成。战斗部是弹体中强度最高、刚度最大的部分,有较大的药室,装填炸药量较多,大多为钝头[1]。本研究建立的前舱物战斗部等效计算模型如图2所示。战斗部壳体模型外径为428.40 mm,内径为384.79 mm,内置装药模型直径为384.79 mm,前舱物等效为与战斗部截面尺寸相等的圆柱体。

    图  2  带前舱物战斗部等效计算模型
    Figure  2.  Equivalent calculation model of projectile with nose cabin

    采用有限元软件LS-DYNA模拟不同工况下带前舱物战斗部侵彻靶板的过程。计算中,结构采用拉格朗日实体单元进行网格划分,时间积分采用显式中心差分法。靶板平面尺寸为3 000 mm × 3 000 mm,边界条件设置为固支约束。计算模型设置为沿中轴面对称,建立1/2模型,战斗部与前舱物采用六面体网格划分。靶板自中心采用过渡网格,中心区域网格尺寸为2.5 mm × 2.5 mm,距离撞击区域越远,网格尺寸越大。

    数值仿真中主要涉及两种材料模型。战斗部壳体和金属靶板采用金属材料常用的Johnson-Cook材料模型,战斗部内置装药和前舱物等效材料采用弹塑性材料常用的弹塑性硬化模型。

    Johnson-Cook材料模型考虑了塑性强化、应变率强化、温升软化效应,能够反映材料在高温、高应变率下力学特性的变化,其具体形式为

    s=(A+Bˉepn)(1+Cln˙ε˙ε0)[1(TT0TmT0)m] (1)

    式中:s为应力;A为静态屈服强度;Bn通过拟合等效应力-应变曲线数据获得;ˉep为有效塑性应变;C为应变率系数;˙ε为有效塑性应变率;˙ε0为参考塑性应变率,一般取˙ε0 = 1 s–1T0为参考温度(取室温);Tm为材料熔点;m为热软化指数。

    战斗部壳体和金属靶板的材料参数根据文献[11-12]并结合工程实践调整得到,如表1所示,其中cp为定压比热容,ρ为密度,μ为泊松比。

    表  1  弹体和靶板的J-C模型参表
    Table  1.  J-C model parameters of projectile body and target
    Materialρ/(g·cm–3E/GPaμcp/(J·kg–1·K–1T0/KTm/K
    Warhead shell7.802050.28400.903001 765
    Target plate7.802000.30452.002981 881
    Material˙ε/s–1A/MPaB/MPanCm
    Warhead shell17605000.530.0141.13
    Target plate13554500.360.0221.00
    下载: 导出CSV 
    | 显示表格

    弹塑性硬化模型考虑了材料应变的等向或随动塑性强化。其应变率效应由Cowper-Symonds模型描述

    σy=[1+(˙εc)1/p](σ0+βEEtEEtεeff)

    式中:σy为屈服强度,˙ε为等效塑性应变率;cpβ为材料常数;σ0为静态屈服强度,E为杨氏模量,Et为硬化模量,εeff为有效塑性应变。内置装药和前舱物等效泡沫材料的模型参数如表2所示,其中fs为失效应变。考虑到不同导弹的前舱物内包含仪器设备的数量和种类不同,前舱物等效材料的屈服强度无法准确给出,本研究只给出取值范围。

    表  2  弹塑性硬化模型材料参数
    Table  2.  Material parameters of the Elastoplastic hardening model
    Materialρ/(kg·m–3E/GPaμσy/MPaEt/MPaβc/s–1pfs
    Charge1 7505 0.3203000005
    Nose cabin 7000.50.30.5–10.0 00000.4
    下载: 导出CSV 
    | 显示表格

    图3为带前舱物战斗部侵彻金属薄板的过程。图4为侵彻过程中典型时刻的靶板内应力分布。典型工况下,带前舱物战斗部与不计前舱物战斗部侵彻金属薄板的速度时程曲线如图5所示。根据金属薄板受力和变形的特点,可以将带前舱物战斗部侵彻金属薄板的过程分为3个阶段。

    图  3  侵彻过程
    Figure  3.  Progress of perforation
    图  4  靶板应力分布
    Figure  4.  Stress distribution of target
    图  5  速度时程曲线
    Figure  5.  Velocity-time curve

    (1)金属薄板结构预变形阶段。从前舱物接触靶板开始,接触区域内最大压强远大于前舱物等效材料的极限强度和金属薄板的屈服强度。因此,前舱物发生变形破坏,金属薄板产生大范围碟形塑性变形,并获得一定的横向速度,战斗部与靶板之间的速度差减小,弹靶接触区边缘的应力集中程度减轻。由图4可知,在金属薄板结构预变形阶段,前舱物与靶板接触区域内正应力较大,且随着与中心距离的增大而降低;剪应力很小,且只存在于接触区域及其附近区域。在此阶段,战斗部的速度变化可以分为两个时期,如图5所示。AB时期,由于前舱物与金属薄板撞击产生的压缩应力波尚未传递到战斗部,战斗部以初速度向前运动;BC时期,随着前舱物与金属薄板撞击产生的压缩应力作用于战斗部,战斗部的速度开始下降,但由于作用于战斗部的压缩应力强度较小,战斗部速度下降较慢。

    (2)冲塞破坏阶段。从战斗部接触金属薄板开始,接触区域由于速度差产生剪切变形,在拉应力与剪应力的联合作用下,接触区域周围产生裂纹,裂纹直径略大于战斗部直径。如图4所示,此阶段战斗部与靶板接触区域内正应力最大,且由于裂纹拓展释放压缩波,接触区域内出现拉伸应力与压缩应力混合的情况;接触区域内剪应力强度远大于结构预变形阶段,剪应力在接触边缘附近达到峰值,并随着与中心距离的增大迅速衰减。如图5所示,在此阶段,战斗部直接侵彻金属薄板,弹体受到的应力强度大,速度下降明显,可近似视为匀减速运动。

    (3)弹体贯穿阶段。随着弹体的进一步运动,裂纹贯穿整个金属薄板厚度方向,形成与战斗部直径相似的帽形失效块,同时,金属薄板由于惯性作用继续发生碟形变形。如图4所示,该阶段由于战斗部与靶板接触的区域形成失效块与靶板脱离,接触区域内正应力与剪应力均为零。靶板动能转化为变形能,挠度增加,靶板仍有较大正应力,由于接触区域附近靶体材料存在速度差,因此靶板接触区域附近存在一定的剪应力。如图5所示,在该阶段,战斗部的速度变化呈现两个时期:DE阶段,由于冲塞块附着在战斗部头部,战斗部推动冲塞块前进,战斗部的一部分动能转化为冲塞块的动能,导致战斗部速度略有下降;EF阶段,冲塞块与战斗部速度达到一致,战斗部向前匀速运动。

    图5中带前舱物战斗部与不计前舱物战斗部的速度时程曲线对比可知,不计前舱物战斗部侵彻靶板的贯穿时间为90 μs,带前舱物战斗部的侵彻贯穿时间为2 180 μs,侵彻作用时间大幅增加,但增加的侵彻时间主要是前舱物被压溃破坏的时间。在金属薄板预变形阶段,金属薄板与前舱物接触区域获得的速度较小,与战斗部之间的速度差较大,因此,战斗部侵彻贯穿金属薄板所需要的时间并未明显增加。此外,相比于不计前舱物战斗部速度变化情况,前舱物的存在使得战斗部在侵彻过程中的受力更加复杂,战斗部速度变化的阶段性和多样性更加明显。

    从塑性变形范围来看,通过有限元仿真软件的计算结果可知,不计前舱物战斗部侵彻靶板过程[13] 中金属薄板塑性变形的直径为883.8 mm,带前舱物战斗部侵彻靶板过程中,由于存在金属薄板预变形阶段,金属薄板迎弹面塑性变形范围大幅增加,直径达到2 629.7 mm。其原因在于:带前舱物战斗部侵彻金属薄板的过靶时间长,弹靶接触区域产生的塑性波传播范围远大于不计前舱物战斗部。因此,前舱物的“缓冲垫”作用明显,在大幅延长战斗部过靶时间的同时大幅扩大了金属薄板塑性变形范围。

    由于不同导弹的前舱物内包含仪器设备的种类和数量不同,因此,无法定量给出前舱物等效材料的屈服强度。本研究给定前舱物等效材料的屈服强度范围为0.5~10.0 MPa,在此屈服强度范围内研究前舱物对平头战斗部侵彻金属薄板的影响。

    图6给出了初速度为340 m/s的不同前舱物等效屈服强度战斗部侵彻20 mm厚金属薄板的应变云图。由图6可知,前舱物等效屈服强度在0.5~10.0 MPa范围内时,金属薄板的塑性应变分布相似,战斗部贯穿金属薄板形成的帽形失效块直径相等。表3列出了不同前舱物等效屈服强度战斗部侵彻靶板时金属薄板的塑性变形范围和战斗部剩余速度。由表3可知,前舱段等效屈服强度增加时,金属薄板塑性变形范围并未出现明显变化,而是在一定范围内波动,最大波动幅度为0.994%。其原因在于,模型计算结果的输出时间间隔为20 μs,侵彻作用时间选取最接近战斗部贯穿靶板的时间,时间选取存在不准确的情况。随着前舱物等效屈服强度的增加,战斗部剩余速度增大,但增大幅度很小,仅为0.147%。因此,可以认为前舱物等效屈服强度对带前舱物平头弹侵彻金属薄板的影响很小,可以忽略。这是由于前舱物存在破坏强度低、延展性差的特点,变形破坏所耗散的变形能远小于带前舱物战斗部所具有的动能,因此当前舱物等效模型的屈服强度远小于靶板时,前舱物的惯性力起主要作用,结构强度的影响可以忽略不计。

    表  3  不同等效屈服强度下侵彻结果的比较
    Table  3.  Comparison of penetration results with different yield stress
    Equivalent yield strength/MPaDistortion range/mmResidual velocity/(m·s–1
    0.52 629.70325.88
    1.02 624.41326.30
    2.02 629.31326.37
    5.02 625.23326.75
    10.02 659.20326.86
    下载: 导出CSV 
    | 显示表格
    图  6  靶板应变云图
    Figure  6.  Strain distribution of the target

    为了分析不同靶板厚度时前舱物对战斗部侵彻的影响,分别模拟了靶板厚度在10~40 mm范围内的7种工况。不同靶板厚度条件下,战斗部剩余速度如图7所示。由图7可知:带前舱物战斗部剩余速度与不计前舱物战斗部的剩余速度均随靶板厚度的增大而减小,且变化规律一致,均近似呈线性变化;靶板厚度相同时,两者的剩余速度相差很小,最大偏差为0.67%。因此,靶板厚度在10~40 mm范围内时,前舱物对战斗部侵彻的影响很小。

    图  7  不同靶板厚度时的剩余速度对比
    Figure  7.  Comparison of residual velocity with different target thickness

    常见反舰半穿甲导弹的飞行马赫数Ma 在 0.8~1.0之间[14],为了分析不同侵彻速度下前舱物对平头战斗部侵彻金属薄板的影响,分别模拟了侵彻速度在238~340 m/s范围内的7种情况。不同侵彻速度条件下,战斗部的剩余速度如图8所示。由图8可知,随着侵彻速度的增大,不计前舱物平头战斗部和带前舱物平头战斗部的剩余速度均近似呈线性增长,且带前舱物平头战斗部的剩余速度略大于不计前舱物平头战斗部(最大偏差为1.15%)。Chen等[6]提出的带前舱物钝头弹侵彻金属靶的计算公式可得到与本模型相似的结论:侵彻速度相同时,带前舱物平头战斗部的剩余速度略大于不计前舱物平头战斗部。

    图  8  不同侵彻速度下的剩余速度对比
    Figure  8.  Comparison of residual velocity with different initial velocity

    通过对不同靶板厚度和侵彻速度条件下战斗部的剩余速度进行分析可知,带前舱物战斗部的剩余速度均大于不计前舱物战斗部,原因在于:与不计前舱物战斗部相比,带前舱物战斗部的质量更大,侵彻速度相同时动能更大。前舱物变形破坏和金属靶板结构预变形消耗部分前舱物具有的动能,使带前舱物战斗部以更大动能侵彻靶板。带前舱物战斗部对靶板的侵彻破坏模式为冲塞破坏,与不带前舱物战斗部侵彻靶板的破坏模式相同,在侵彻过程中消耗的能量相差不大。因此,带前舱物战斗部的剩余速度略大于不计前舱物战斗部。

    本研究建立了变截面刚度弹体模型,对带前舱物战斗部侵彻金属薄板进行仿真分析。通过对比不同前舱物等效屈服强度战斗部的侵彻性能,以及不同靶板厚度与初速度条件下带前舱物战斗部与不计前舱物战斗部的剩余速度,得出如下主要结论。

    (1)根据靶板的受力变形特征,带前舱物平头战斗部侵彻金属薄板可以分为3个阶段,即靶板预变形阶段、冲塞破坏阶段和弹体贯穿阶段。虽然带前舱物平头战斗部与不计前舱物平头战斗部侵彻靶板的过程存在明显差异,但靶板破坏模式均为剪切冲塞破坏。

    (2)当前舱物等效材料的强度远小于靶板的屈服强度时,前舱物等效泡沫铝材料的屈服强度对战斗部侵彻性能的影响可以忽略不计。

    (3)在导弹侵彻目标的过程中,一方面,前舱物的存在增大了弹体的总动能,且先于战斗部撞击目标,造成目标结构变形甚至破坏,对战斗部侵彻靶板起到了积极作用;另一方面,由于前舱物的强度远小于战斗部,前舱物的存在延长了弹体的过靶时间,降低了靶板过载,起到了“缓冲垫”的作用,不利于平头战斗部对目标的侵彻。在这两方面综合影响下,当飞行马赫数Ma 在 0.8~1.0之间,侵彻厚度为10~40 mm的靶板时,带前舱物平头战斗部的剩余速度比不计前舱物平头战斗部均有所增加,但增加幅度小于2%。

  • 图  (a) Cr2B3 和 (b) CrMB3(M = Sc, Ti, V, Zr)的晶体结构(蓝色、绿色和粉色小球分别代表 Cr、B 和过渡金属原子)

    Figure  1.  Crystal structures of (a) Cr2B3 and (b) CrMB3 (M=Sc, Ti, V, Zr), where the blue, green and pink spheres represent Cr, B and transition metal atoms, respectively

    图  Cr2B3的弹性常数随压力的变化

    Figure  2.  Pressure dependence of the elastic constants for Cr2B3

    图  CrMB3(M=Sc, Ti, V, Zr)的弹性常数随压力的变化

    Figure  3.  Pressure dependence of the elastic constants for CrMB3 (M=Sc, Ti, V, Zr)

    图  Cr2B3在零压(a)和 150 GPa (b)下的声子色散曲线

    Figure  4.  Phonon-dispersion curves of Cr2B3 at (a) 0 GPa and (b) 150 GPa

    图  (a) CrScB3、(b) CrTiB3、(c) CrVB3、(d) CrZrB3在零压下的声子色散曲线

    Figure  5.  Phonon-dispersion curves of (a) CrScB3, (b) CrTiB3, (c) CrVB3 and (d) CrZrB3 at 0 GPa

    图  Cr2B3和 CrMB3(M=Sc, Ti, V, Zr)的体积弹性模量(a)、剪切模量(b)和杨氏模量(c)随压力的变化

    Figure  6.  Pressure dependence of (a) the bulk moduli, (b) the shear moduli, (c) the Young’s moduli for Cr2B3 and CrMB3 (M=Sc, Ti, V, Zr)

    图  Cr2B3和 CrMB3(M=Sc, Ti, V, Zr)的维氏硬度随压力的变化

    Figure  7.  Pressure dependence of the Vicker’s hardness for Cr2B3 and CrMB3 (M=Sc, Ti, V, Zr)

    图  Cr2B3及 CrMB3(M=Sc, Ti, V, Zr)的B/G 随压力的变化

    Figure  8.  Pressure dependence of the B/G for Cr2B3 and CrMB3 (M=Sc, Ti, V, Zr)

    图  零压下 Cr2B3的体积弹性模量的三维表示(a)及其在 xyxzyz 面上的投影(b)

    Figure  9.  (a) 3D representation and (b) 2D projections on xy, xz and yz planes of the bulk modulus for Cr2B3 at 0 GPa

    图  10  零压下 Cr2B3的杨氏模量的三维表示(a)及其在 xyxzyz 面上的投影(b)

    Figure  10.  (a) 3D representation and (b) 2D projections on xy, xz and yz planes of the Young’s modulus for Cr2B3 at 0 GPa

    图  11  零压下 CrMB3(M=Sc, Ti, V, Zr)的体积弹性模量在 xyxzyz 面的投影

    Figure  11.  2D projections of the bulk modulus for CrMB3 (M=Sc, Ti, V, Zr) on xy, xz and yz planes at 0 GPa

    图  12  零压下 CrMB3(M=Sc, Ti, V, Zr)的杨氏模量在 xyxzyz 面的投影

    Figure  12.  2D projections of the Young’s modulus for CrMB3 (M=Sc, Ti, V, Zr) on xy, xz and yz planes at 0 GPa

    图  13  零压下 Cr2B3的总态密度以及分波态密度

    Figure  13.  Total density of states and partial density of states for Cr2B3 at 0 GPa

    图  14  零压下 CrMB3(M=Sc, Ti, V, Zr)的总态密度以及分波态密度

    Figure  14.  Total density of states and partial density of states for CrMB3 (M=Sc, Ti, V, Zr) at 0 GPa

    图  15  零压下(a) Cr2B3、(b) CrScB3、(c) CrTiB3、(d) CrVB3、(e) CrZrB3在(100)平面的 ELF 以及零压下(f) Cr2B3、(g) CrScB3、(h) CrTiB3、(i) CrVB3和(j) CrZrB3在(00¯1)平面的 ELF

    Figure  15.  Electronic local functions contours for (a) Cr2B3, (b) CrScB3, (c) CrTiB3, (d) CrVB3 and (e) CrZrB3 in plane (100) at 0 GPa, and electronic local functions contours for (f) Cr2B3, (g) CrScB3, (h) CrTiB3, (i) CrVB3 and (j) CrZrB3 in plane (00¯1) at 0 GPa

    表  1  零压下Cr2B3及掺杂结构CrMB3(M=Sc, Ti, V, Zr)的晶格常数、形成焓及掺杂结构的形成能

    Table  1.   Lattice constants, formation enthalpy of Cr2B3 and CrMB3 (M = Sc, Ti, V, Zr), and impurity formation energy ofCrMB3 (M=Sc, Ti, V, Zr) at zero pressure

    Compounds Doping-site positionSpace groupLattice constantsEf/eVΔH/(eV∙atom−1)
    abc
    Cr2B3 Cmcm2.898318.04642.9286−0.4731
    CrScB3 Cr1Cmcm3.220718.52223.0293−0.9713−0.6653
    Cr2
    CrTiB3 Cr1Cmcm3.064818.16692.9746−1.9192−0.8561
    Cr2Cmcm3.054118.59002.9883−0.6876−0.6118
    CrVB3 Cr1Cmcm2.954618.08742.9436−1.1856−0.7092
    Cr2Cmcm2.954718.22542.9489−0.6797−0.6101
    CrZrB3 Cr1Cmcm3.297218.58723.0702−1.2236−0.7181
    Cr2Cmcm3.176320.01363.0524 0.1705−0.4424
    下载: 导出CSV

    表  2  Cr2B3中B―B键的键长以及布居数随压力的变化

    Table  2.   Pressure dependence of B―B bond length and population for Cr2B3

    Pressure/GPaBond length/Å Population
    B1―B1B1―B2B2―B3B1―B1B1―B2B2―B3
    01.695861.719711.75978 1.590.631.36
    251.661811.683551.723691.620.641.40
    501.634741.655591.694621.650.651.44
    751.612371.632141.670121.680.661.47
    1001.592981.612251.649311.700.671.50
    1251.575891.595061.630981.730.681.53
    1501.560581.579501.614421.750.691.56
    下载: 导出CSV
  • [1] CAO A H, ZHAO W J, ZHOU Q Y, et al. A superhard allotrope of carbon: ibam-C and its BN phase [J]. Chemical Physics Letters, 2019, 714: 119–124. doi: 10.1016/j.cplett.2018.10.079
    [2] FENG S Q, YANG Y, GUO F, et al. Structural, elastic, electronic and hardness properties of osmium diboride predicted from first principles calculations [J]. Journal of Alloys and Compounds, 2020, 844: 156098. doi: 10.1016/j.jallcom.2020.156098
    [3] WANG C C, TAO Q, DONG S S, et al. Synthesis and mechanical character of hexagonal phase δ-WN [J]. Inorganic Chemistry, 2017, 56(7): 3970–3975. doi: 10.1021/acs.inorgchem.6b03041
    [4] CAI Y X, XIONG J M, LIU Y B, et al. Electronic structure and chemical hydrogen storage of a porous sp3 tetragonal BC2N compound [J]. Journal of Alloys and Compounds, 2017, 724: 229–233. doi: 10.1016/j.jallcom.2017.06.343
    [5] MOHAMMADI R, LECH A T, XIE M, et al. Tungsten tetraboride, an inexpensive superhard material [J]. Proceedings of the National Academy of Sciences of the United States of America, 2011, 108(27): 10958–10962. doi: 10.1073/pnas.1102636108
    [6] KANER R B, GILMAN J J, TOLBERT S H. Designing superhard materials [J]. Science, 2005, 308(5726): 1268–1269. doi: 10.1126/science.1109830
    [7] GOU H Y, LI Z P, NIU H, et al. Unusual rigidity and ideal strength of CrB4 and MnB4 [J]. Applied Physics Letters, 2012, 100(11): 111907. doi: 10.1063/1.3692777
    [8] CHONG X Y, JIANG Y H, ZHOU R, et al. Elastic properties and electronic structures of Cr xB y as superhard compounds [J]. Journal of Alloys and Compounds, 2014, 610: 684–694. doi: 10.1016/j.jallcom.2014.05.010
    [9] ANDERSSON S, LUNDSTRÖM T. The crystal structure of CrB4 [J]. Acta Chemica Scandinavica, 1968, 22(10): 3103–3110.
    [10] KOTZOTT D, ADE M, HILLEBRECHT H. Synthesis and crystal structures of α- and β- modifications of Cr2IrB2 containing 4-membered B4 chain fragments, the τ-boride Cr7.9Ir14.1B6 and orthorhombic Cr2B [J]. Solid State Sciences, 2008, 10(3): 291–302. doi: 10.1016/j.solidstatesciences.2007.09.014
    [11] OKADA S, ATODA T, HIGASHI I. Structural investigation of Cr2B3, Cr3B4, and CrB by single-crystal diffractometry [J]. Journal of Solid State Chemistry, 1987, 68(1): 61–67. doi: 10.1016/0022-4596(87)90285-4
    [12] GIANOGLIO C, PRADELLI G, VALLINO M. Solid state equilibria in the Cr-Fe-B system at the temperature of 1 373 K [J]. Metallurgical Science and Tecnology, 1983, 1(2): 51–57.
    [13] WONG-NG W, MCMURDIE H F, PARETZKIN B, et al. Reference X-ray diffraction powder patterns of fifteen ceramic phases [J]. Powder Diffraction, 1987, 2(4): 257–265. doi: 10.1017/S0885715600012926
    [14] NIU H Y, WANG J Q, CHEN X Q, et al. Structure, bonding, and possible superhardness of CrB4 [J]. Physical Review B, 2012, 85(14): 144116. doi: 10.1103/PhysRevB.85.144116
    [15] ZHANG Y K, WU L L, WAN B, et al. Structural variety beyond appearance: high-pressure phases of CrB4 in comparison with FeB4 [J]. Physical Chemistry Chemical Physics, 2016, 18(4): 2361–2368. doi: 10.1039/C5CP06745F
    [16] WANG S, YU X, ZHANG J, et al. Crystal structures, elastic properties, and hardness of high-pressure synthesized CrB2 and CrB4 [J]. Journal of Superhard Materials, 2014, 36(4): 279–287. doi: 10.3103/S1063457614040066
    [17] OKADA S, KUDOU K, IIZUMI K, et al. Single-crystal growth and properties of CrB, Cr3B4, Cr2B3 and CrB2 from high-temperature aluminum solutions [J]. Journal of Crystal Growth, 1996, 166(1/2/3/4): 429–435.
    [18] MIAO N H, SA B S, ZHOU J, et al. Theoretical investigation on the transition-metal borides with Ta3B4-type structure: a class of hard and refractory materials [J]. Computational Materials Science, 2011, 50(4): 1559–1566. doi: 10.1016/j.commatsci.2010.12.015
    [19] XING W D, MENG F Y, YU R. Strengthening materials by changing the number of valence electrons [J]. Computational Materials Science, 2017, 129: 252–258. doi: 10.1016/j.commatsci.2016.12.037
    [20] ZHANG Y M, LIU D, ZHAO Y H, et al. Physical properties and electronic structure of Cr2B under pressure [J]. Physica Status Solidi (B), 2021, 258(2): 2000212. doi: 10.1002/pssb.202000212
    [21] DOVALE-FARELO V, TAVADZE P, VERSTRAETE M J, et al. Exploring the elastic and electronic properties of chromium molybdenum diboride alloys [J]. Journal of Alloys and Compounds, 2021, 866: 158885. doi: 10.1016/j.jallcom.2021.158885
    [22] OKADA S, ATODA T, HIGASHI I, et al. Preparation of single crystals of a new boride Cr2B3 by the aluminium-flux technique and some of its properties [J]. Journal of the Less Common Metals, 1985, 113(2): 331–339. doi: 10.1016/0022-5088(85)90289-9
    [23] WATANABE K, SAKAIRI M, TAKAHASHI H, et al. Formation of Al-Zr composite oxide films on aluminum by sol-gel coating and anodizing [J]. Journal of Electroanalytical Chemistry, 1999, 473(1/2): 250–255.
    [24] PERDEW J P, RUZSINSZKY A, CSONKA G I, et al. Restoring the density-gradient expansion for exchange in solids and surfaces [J]. Physical Review Letters, 2008, 100(13): 136406. doi: 10.1103/PhysRevLett.100.136406
    [25] VANDERBILT D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism [J]. Physical Review B, 1990, 41(11): 7892–7895. doi: 10.1103/PhysRevB.41.7892
    [26] MONKHORST H J, PACK J D. Special points for Brillouin-zone integrations [J]. Physical Review B, 1976, 13(12): 5188–5192. doi: 10.1103/PhysRevB.13.5188
    [27] VAN DE WALLE C G, NEUGEBAUER J. First-principles calculations for defects and impurities: applications to Ⅲ-nitrides [J]. Journal of Applied Physics, 2004, 95(8): 3851–3879. doi: 10.1063/1.1682673
    [28] WU Z J, ZHAO E J, XIANG H P, et al. Crystal structures and elastic properties of superhard IrN2 and IrN3 from first principles [J]. Physical Review B, 2007, 76(5): 054115. doi: 10.1103/PhysRevB.76.054115
    [29] HILL R. The elastic behaviour of a crystalline aggregate [J]. Proceedings of the Physical Society: Section A, 1952, 65(5): 349–354. doi: 10.1088/0370-1298/65/5/307
    [30] CHEN X Q, NIU H Y, LI D Z, et al. Modeling hardness of polycrystalline materials and bulk metallic glasses [J]. Intermetallics, 2011, 19(9): 1275–1281. doi: 10.1016/j.intermet.2011.03.026
    [31] PUGH S F. XCII. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals [J]. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1954, 45(367): 823–843. doi: 10.1080/14786440808520496
    [32] SILVI B, SAVIN A. Classification of chemical bonds based on topological analysis of electron localization functions [J]. Nature, 1994, 371(6499): 683–686. doi: 10.1038/371683a0
  • 加载中
图(17) / 表(2)
计量
  • 文章访问数:  504
  • HTML全文浏览量:  662
  • PDF下载量:  89
出版历程
  • 收稿日期:  2021-12-10
  • 修回日期:  2022-02-06
  • 录用日期:  2022-02-06
  • 网络出版日期:  2022-06-30
  • 刊出日期:  2022-07-28

目录

/

返回文章
返回