极低温下Nb3Sn超导体单晶裂纹动态扩展模拟

王豪阳 卫颖 乔力

王豪阳, 卫颖, 乔力. 极低温下Nb3Sn超导体单晶裂纹动态扩展模拟[J]. 高压物理学报, 2022, 36(3): 034201. doi: 10.11858/gywlxb.20210884
引用本文: 王豪阳, 卫颖, 乔力. 极低温下Nb3Sn超导体单晶裂纹动态扩展模拟[J]. 高压物理学报, 2022, 36(3): 034201. doi: 10.11858/gywlxb.20210884
WANG Haoyang, WEI Ying, QIAO Li. Simulation of Dynamic Crack Propagation in Superconducting Nb3Sn at Extreme Low Temperature[J]. Chinese Journal of High Pressure Physics, 2022, 36(3): 034201. doi: 10.11858/gywlxb.20210884
Citation: WANG Haoyang, WEI Ying, QIAO Li. Simulation of Dynamic Crack Propagation in Superconducting Nb3Sn at Extreme Low Temperature[J]. Chinese Journal of High Pressure Physics, 2022, 36(3): 034201. doi: 10.11858/gywlxb.20210884

极低温下Nb3Sn超导体单晶裂纹动态扩展模拟

doi: 10.11858/gywlxb.20210884
基金项目: 国家自然科学基金(11772212,11402159)
详细信息
    作者简介:

    王豪阳(1998-),男,硕士研究生,主要从事电磁固体力学研究. E-mail:why18734584608@163.com

    通讯作者:

    乔 力(1984-),男,教授,主要从事微纳米尺度材料力学和电磁固体力学研究.E-mail:qiaoli@tyut.edu.cn

  • 中图分类号: O511.3; O521.2

Simulation of Dynamic Crack Propagation in Superconducting Nb3Sn at Extreme Low Temperature

  • 摘要: 研究Nb3Sn超导体的损伤断裂行为对于揭示超导临界性能弱化背后的力学机制具有重要的意义。采用分子动力学模拟方法,研究了极低温下不含裂纹和含中心裂纹的Nb3Sn单晶在力学拉伸变形作用下的断裂机制和裂纹扩展行为,同时分析了应变率效应对Nb3Sn单晶断裂机制与裂纹扩展行为的影响。结果表明:不含裂纹的Nb3Sn单晶在结构受力后出现滑移,滑移带上位错塞积导致应力集中,应力集中使原子键断裂从而萌生裂纹致使Nb3Sn单晶断裂;而含中心裂纹的Nb3Sn单晶则由于裂纹尖端应力集中使得原子键断裂形成微裂纹,裂纹扩展致使Nb3Sn单晶断裂。Nb3Sn单晶在不同的应变率下表现出不同的断裂机制,在低应变率下表现为脆性断裂,而在高应变率下表现为韧性断裂。

     

  • Nb3Sn超导体在磁约束核聚变以及高能物理超导磁体领域有着广泛的应用需求[1-4]。Nb3Sn超导磁体在极低温(4.2 K)、高电流(68 kA)、强磁场(12 T)的极端工况下运行,由于Nb3Sn超导体的制备温度约为943 K,而使用环境温度约为4.2 K,巨大的温度差异导致残余应力的产生;另外,强磁场对载流股线产生洛伦兹力,在温度载荷与电磁载荷的共同作用下,Nb3Sn超导体会发生损伤断裂。Nb3Sn超导体的损伤断裂是其临界性能不可逆退化的主要原因,研究Nb3Sn超导体的损伤断裂行为对于揭示超导临界性能弱化背后的力学机制具有重要的意义,同时也为Nb3Sn强磁场超导磁体装置的设计和安全运行提供一定的理论指导。

    Sanabria等[5]利用金相图像分析方法,对国际热核聚变实验堆(international thermonuclear experimental reactor, ITER)管内电缆导体(cable in conduit conductor, CICC)进行分析,探究了洛伦兹力引起的股线运动对Nb3Sn细丝的影响,结果表明局部Nb3Sn细丝断裂是由高洛伦兹力引起的。Sheth等[6]对ITER环向磁场股线在单轴循环加载下的纤维开裂进行了研究,发现青铜工艺股线经过104次加载循环后,在1%的应变下开始出现显著的长丝开裂,而内锡工艺股线经过103次加载循环后,在0.7%的应变下开始出现显著的长丝开裂,说明制造工艺对股线性能产生明显影响。Lu等[7]在对Nb3Sn股线的不可逆应变研究中,利用直拉装置测量了由Luvata Waterbury公司(Luvata)和牛津仪器超导技术公司(Oxford Instrument Superconducting Technology,OST)制造的ITER导线的不可逆应变,结果表明:与OST细丝相比,Luvata细丝的不可逆应变较大。Nijhuis等[8]在4.2 K下测量了2种股线由完好直至发生破坏的轴向应力-应变关系,发现内锡股线在应变为0.90%(300 MPa)时发生破坏,青铜股线在应变为1.35%(380 MPa)时发生破坏。这种股线与股线之间的性能差异,与其几何构型以及微观结构缺陷有关[7, 9]。目前大部分研究集中在Nb3Sn超导复合体的断裂上,并将Nb3Sn复合体(基体、二次相和孔洞等缺陷)作为一个完整的结构单元,然而这种处理难以描述Nb3Sn的固有断裂性能,也很难将这些结果应用于Nb3Sn超导丝在使用条件下的临界电流随应变变化的预测模型中。为了研究Nb3Sn的固有断裂性能,Dylla等[10]从熔合型Nb3Sn复合线材中提取出单独的Nb3Sn细丝,利用一种测量小直径长丝拉伸性能的新方法对Nb3Sn细丝进行拉伸测试,结果表明:在室温下,当应变小于0.4%时,Nb3Sn细丝发生断裂。

    在数值模拟研究方面,Zhai等[11]发展了新的超导电缆数学模型,研究了热载荷和电磁力所致的弯曲变形以及超导丝断裂对超导电缆超导性能的影响,结果表明超导丝断裂是超导电缆超导性能退化的主要原因。Wang等[12]通过对超导线缆的多层级建模及力-电-磁行为研究,探讨了Nb3Sn超导股线在轴向拉伸和纯弯曲作用下超导丝断裂对股线临界性能的影响,计算结果表明,股线中超导丝的断裂对股线的临界电流具有非常重要的影响。人们还从微观角度对Nb3Sn的力-电和力-热耦合行为进行了模拟研究[13-16],得到了极低温环境及宽温度域环境下超导单晶体的基本弹性力学性能参数,然而受第一性原理计算方法在模拟体系上的限制,这些研究并未涉及超导体单晶的损伤演化行为。

    实验结果表明,Nb3Sn超导体的超导电性对力学响应(如应力、应变)非常敏感,外载作用下超导临界性能退化存在可逆变化向不可逆变化的转变。随着拉伸应变的增加,存在某应变极限值,当应变超过该极限值时,轴向变形导致的Nb3Sn超导体临界电流退化具有不可逆性,进而引起超导体性能的永久退化,并伴随着超导体中基体材料的塑性变形和脆性Nb3Sn超导丝的断裂。当Nb3Sn股线处于服役状态时,复杂应力状态下的超导丝断裂会导致超导临界性能的不可逆退化,进而引发失超。鉴于Nb3Sn超导体复杂的多尺度结构特征,很难观察其变形和失效过程全貌,加之测量环境苛刻,极低温下裂纹萌生和演化的精确观测难度很大,因此需要从微观尺度对Nb3Sn超导体的断裂行为进行模拟,以补充宏观实验的观测结果。此外,在Nb3Sn超导体制备过程中会不可避免地引入各种缺陷,如位错、空位、孔洞、晶界、裂纹等,这些缺陷的存在会改变原子间结合能,直接影响材料的力学性能和超导临界性能,另外实际工程中Nb3Sn超导丝断裂引起的超导临界性能退化也需要以极低温下含缺陷Nb3Sn超导体的变形和断裂行为为基础进行研究。Ren等[17]研究了Nb3Sn单晶和多晶样品在不同压力下的电阻率和超导性能变化,发现随着压力的增大,Nb3Sn单晶和多晶样品的超导临界温度均降低,但是压力响应不同,单晶样品比多晶样品的抗压性更强。Nb3Sn单晶和多晶的损伤断裂对其超导临界性能具有一定的影响,但是对这种影响尚没有定量阐述,为此需要通过微观尺度模拟细致地刻画裂纹萌生和扩展过程。

    分子动力学方法可以从微观层次通过原子间的相互作用力研究材料的变形和断裂过程。在给定的初始条件和边界条件下,对利用经典力学方法建立的体系运动方程进行数值求解,得到体系原子的位置、速度等信息,进而得到粒子系统在空间中随时间演化的轨迹。与量子力学计算方法相比,分子动力学方法具有节省计算资源、模拟体系大等特点[18],在超导体研究方面应用广泛。例如:Chudinov 等[19]采用分子动力学方法模拟了Nb3Sn化合物中的级联碰撞过程;陈伟[20]借助分子动力学方法研究了第Ⅰ类超导磁体系统的静态磁畴结构及其输运性质,在一定程度上解释了实验现象。

    本研究通过建立Nb3Sn单晶模型,采用分子动力学方法研究不含裂纹和含中心裂纹的Nb3Sn单晶在力学拉伸作用下的断裂机制和裂纹扩展行为,同时考虑应变率效应对Nb3Sn单晶断裂机制和裂纹扩展行为的影响,研究结果可为Nb3Sn强磁场超导磁体装置的设计与安全运行提供一定的理论指导。

    Nb3Sn具有A15相晶体结构(见图1(a)),通过扫描透射电镜实验,得到Nb3Sn晶格结构的高角度环形暗场(high-angle annular dark field,HAADF)图像[21],如图1(b)所示,从中可以清晰地观察到原子排列特征。基于实验图像,利用LAMMPS软件建立不含裂纹和含中心裂纹的Nb3Sn单晶模型,本研究采用60a×60a×5a计算模型(a为Nb3Sn的晶格常数)。图2(a)为Nb3Sn单晶立体模型;图2(b)为不含裂纹的Nb3Sn单晶模型,体系中有144000个原子;图2(c)为含中心裂纹的Nb3Sn单晶模型,模型大小也为60a×60a×5a,通过删除原子法预制长度为10a的裂纹,体系中有143015个原子。需要说明的是,受裂纹密度的影响,理论上应该选取边界更大的计算模型,以避免边界条件对裂纹扩展模拟结果的影响,为了清楚地显示裂纹的成核和扩展过程,并综合考虑计算效率,本研究选择上述模型。xyz轴分别对应晶体的[100][010][001]晶向,为了减小边界效应的影响,在xyz方向上采用周期性边界条件。

    图  1  Nb3Sn的晶格结构
    Figure  1.  Lattice structure of Nb3Sn
    图  2  Nb3Sn单晶的分子动力学计算模型
    Figure  2.  Molecular dynamics simulation model of Nb3Sn single crystal

    采用LAMMPS程序对不含裂纹和含中心裂纹的Nb3Sn单晶进行拉伸模拟,模拟过程主要包括弛豫阶段和加载阶段。在模拟过程中,采用共轭梯度法对不含裂纹和含中心裂纹的Nb3Sn单晶模型进行能量最小化,以获得初始平衡结构,之后在等温等压的NPT系综下进行弛豫,弛豫结束后,在NVT系综下沿x轴拉伸。选取应变率时,借鉴了众多学者利用分子动力学方法模拟金属材料及金属间化合物时的应变率(加载应变率基本在108~1011 s−1[22-24],结合数值试验结果,应变率在109 s−1以下时Nb3Sn的计算结果对应变率不敏感,为此设定5×108 s−1作为常温下Nb3Sn单晶拉伸模拟和4.2 K下Nb3Sn超导体单晶裂纹扩展模拟时的应变率;选取5×108、1×109、5×109、1×1010、5×1010 s−1 5组应变率研究Nb3Sn超导体单晶动态断裂的应变率效应,模拟时间步长设为0.001 ps,模拟过程中每隔一段时间输出体系的应变、应力等力学参数,模拟结果通过后处理软件OVITO进行可视化分析。

    原子之间的相互作用势是所有原子水平计算模拟的基础,它直接影响分子动力学模拟结果的准确性。本研究采用组合势函数的方法描述Nb3Sn单晶体系中原子之间的相互作用力,Nb与Nb之间的势函数采用Zhang等[25]开发的用于Ni62Nb38的原子嵌入势,Sn与Sn之间的势函数为MEAM势函数[26],Nb与Sn之间的势函数为L-J势,其中的参数来自文献[19]。为了验证势函数的可靠性,计算了在4.2 K、0 GPa下Nb3Sn单晶的弹性常数(C11C12C44)和晶格常数,并与第一性原理模拟结果[13-15, 27]进行比较,如表1所示,误差在可接受范围内,表明该势函数在表征力学性能参数上准确可靠。

    表  1  Nb3Sn单晶的弹性常数和晶格常数
    Table  1.  Elastic constants and lattice constant of Nb3Sn single crystal
    MethodC11/GPaC12/GPaC44/GPaa
    This work284.10 95.8453.765.21
    First principle284.32107.7067.075.32
    下载: 导出CSV 
    | 显示表格

    在分子动力学模拟中,原子尺度下的应力计算与宏观尺度下的应力计算不同。以下从能量原理出发,给出原子应力的Virial表达式[28]。对于一个由N个原子组成的体积为Ω的体系,其总应力等于所有原子维里之和除以总体积

    Παβ = 1ΩNi=1wαβi
    (1)

    i个原子维里为

    wαβi=mivαivβi12ji(1rijϕij(rij)r)rαijrβij+...
    (2)

    式中:αβ为笛卡尔坐标,ϕij为原子间相互作用势,rij为第 i 个原子与第j个原子之间的距离,mivi 为第 i 个原子的质量和速度,省略号代表分子模拟中可能出现的高阶原子维里贡献。选取合适的原子体积 ωi,如Wigner-Seitz/Voronoi体积,则第i个原子的维里应力可定义为

    παβi=wαβi/ωi
    (3)

    本研究对不含裂纹模型的144000个原子和含中心裂纹模型的143015个原子的单原子应力进行了计算,原子体积采用Voronoi体积,其中Nb原子体积为17.7 Å3,Sn原子体积为17.2 Å3

    常温300 K下,对不含裂纹和含中心裂纹的Nb3Sn单晶进行拉伸模拟,得到的应力-应变曲线如图3所示。 目前还没有Nb3Sn单晶在极低温下的应力-应变曲线报道,为了从侧面间接验证势函数和模拟结果的可靠性,在图3中给出了常温下Nb3Sn细丝的拉伸实验曲线[10]

    图  3  300 K下不含裂纹和含中心裂纹的Nb3Sn单晶拉伸应力-应变曲线
    Figure  3.  Tensile stress-strain curves of Nb3Sn single crystal without crack and with central crack at 300 K

    通过比较可以发现,模拟得到的应力-应变曲线与实验测量结果定性一致,表明本研究所用的势函数具有一定的可行性。常温下模拟得到的断裂应变和拉伸强度均高于实验值,其原因是模拟采用了理想的Nb3Sn单晶,而实验采用的样品为含有裂纹、孔洞、杂质等缺陷的Nb3Sn多晶体。对比不含裂纹与含中心裂纹的Nb3Sn单晶的拉伸结果,可以发现,裂纹的存在会降低断裂应变和拉伸强度。

    为了研究极低温下Nb3Sn单晶的断裂,在4.2 K下对不含裂纹和含中心裂纹的Nb3Sn单晶进行拉伸模拟,得到的应力-应变曲线如图4所示。对于不含裂纹的Nb3Sn单晶,在初始弹性阶段,应力随应变呈线性增加,表明在微观尺度下Nb3Sn单晶的弹性变形仍符合胡克定律,在此过程中原子间的距离发生微小的变化,晶格排列完整。弹性阶段后出现“拐点”(图4中红色箭头所指处),“拐点”的产生是由于在极低温下原子动能较小,且在较低的应变率下晶格只能发生微量塑性变形。当应变ε为0.110时,应力达到最大值,随着应变继续增加,应力急剧下降,原子间作用力减弱直至消失,Nb3Sn发生断裂。对于含中心裂纹的Nb3Sn单晶,在初始弹性阶段,应力随应变线性增加,当应变为0.045时,应力达到最大值,此后应力急剧下降,Nb3Sn发生断裂。在极低温下,通过对比不含裂纹和含中心裂纹的Nb3Sn单晶的拉伸结果,可以发现,裂纹的存在会降低断裂应变和拉伸强度。

    图  4  4.2 K下Nb3Sn单晶的拉伸应力-应变曲线
    Figure  4.  Tensile stress-strain curves of Nb3Sn single crystal at 4.2 K

    为了观察极低温下Nb3Sn单晶的裂纹萌生和扩展过程,借助OVITO软件进行可视化处理,得到了温度为4.2 K、应变率为5×108 s−1时不含裂纹和含中心裂纹的原子结构演化图像,如图5图6所示。从图5可以看出:不含裂纹的Nb3Sn单晶在初始拉伸时发生弹性变形,在此过程中,原子排列整齐,结构内部有微量塑性变形;当应变ε>0.110时,位错开始出现(由脆性断裂理论[29-30]可知,脆性裂纹的成核与切应力导致的滑移是分不开的,对于单晶体,在拉伸作用力下,沿着最大切应力的滑移面滑移,单轴拉伸时最大切应力方向一般与拉伸轴成45°角),与拉伸轴成45°方向上出现滑移带;当应变ε=0.122时,滑移带上出现孔洞,随着应变的增加,孔洞增多并逐渐变大,随后孔洞融合形成裂纹,裂纹扩展最终使Nb3Sn单晶完全断裂。从图5还可以看到,裂纹两边的晶格结构仍然保持有序,没有明显的塑性变形,Nb3Sn单晶发生脆性断裂,模拟结果与实验观测到的Nb3Sn的脆性断裂特征[31-32]吻合。从图6可以看出:含中心裂纹的Nb3Sn单晶在初始拉伸时(0<ε<0.045)同样只发生弹性变形,晶格结构未发生变化;随着应变的增加,裂纹逐渐变长变宽,裂纹尖端发生塑性变形,原子乱序排列,同时裂纹尖端发射位错,结构内部出现少量滑移,之后裂纹扩展形成孔洞,孔洞与之前的中心裂纹融合形成长裂纹,最终导致Nb3Sn单晶发生断裂。

    图  5  4.2 K下不含裂纹的Nb3Sn单晶的原子结构演化
    Figure  5.  Atomic structure evolution of Nb3Sn single crystal without crack at 4.2 K
    图  6  4.2 K下含中心裂纹的Nb3Sn单晶的原子结构演化
    Figure  6.  Atomic structure evolution of Nb3Sn single crystal with central crack at 4.2 K

    为了进一步解释不含裂纹和含中心裂纹的Nb3Sn单晶断裂过程的不同,导出拉伸方向上的原子应力云图,如图7图8所示。观察图7发现,对于不含裂纹的Nb3Sn单晶,当应变ε>0.110时,位错在滑移带上前进,当位错运动到滑移带的交接处和边界时,由于边界阻塞和位错交互作用,出现应力集中现象。从图7中还可以发现,滑移带上的原子应力高于其他原子,应力集中使原子键发生断裂,从而萌生缺陷,出现裂纹。而对于含中心裂纹的Nb3Sn单晶,如图8所示,裂纹尖端的原子应力较高,这是由于随着应变的增加,裂纹开始扩展,裂纹尖端发生塑性变形,塑性区出现位错塞积,当塞积数目较大时,塞积群的应力集中导致原子键发生断裂,因此在裂纹扩展方向上出现孔洞,孔洞与之前的中心裂纹融合形成长裂纹,长裂纹继续扩展直至单晶体完全断裂。对比不含裂纹和含中心裂纹的Nb3Sn单晶断裂过程发现:不含裂纹的Nb3Sn单晶是结构受力后出现滑移,滑移带上的位错塞积导致应力集中,进而使原子键断裂萌生裂纹,Nb3Sn单晶断裂;而含中心裂纹的Nb3Sn单晶则是裂纹尖端应力集中,使得原子键断裂形成微裂纹,微裂纹扩展致使材料断裂。对于含中心裂纹的Nb3Sn单晶,当应变ε>0.045时,裂纹尖端出现滑移带,在滑移带上可以发现应力集中现象,但是在材料断裂前滑移带上没有萌生缺陷,裂纹自身扩展是Nb3Sn单晶断裂的主要原因,说明中心裂纹的存在加速了Nb3Sn单晶的断裂,同时缺陷的存在使材料更容易失效。

    图  7  4.2 K下不含裂纹的Nb3Sn单晶在拉伸方向的原子应力云图
    Figure  7.  Atomic stress distribution of Nb3Sn single crystal without crack at 4.2 K
    图  8  4.2 K下含中心裂纹的Nb3Sn单晶在拉伸方向的原子应力云图
    Figure  8.  Atomic stress distribution of Nb3Sn single crystal with central crack at 4.2 K

    为了分析极低温下不同应变率对不含裂纹和含中心裂纹的Nb3Sn单晶断裂机制和裂纹扩展行为的影响,设定5组应变率,分别为5×108、1×109、5×109、1×1010、5×1010 s−1,控制体系温度在4.2 K,对不含裂纹和含中心裂纹的Nb3Sn单晶进行不同应变率下的拉伸模拟,得到的应力-应变曲线如图9所示。可以看出:在初始弹性阶段,斜率基本相同,表明应变率变化对Nb3Sn单晶弹性变形行为几乎没有影响;随着应变率的增加,不含裂纹和含中心裂纹的Nb3Sn单晶的拉伸曲线均逐渐表现出塑性特性;在不同的应变率下,不含裂纹的Nb3Sn单晶的拉伸强度没有太大差异,而含中心裂纹的Nb3Sn单晶的拉伸强度则随着应变率的增加而增大。

    图  9  4.2 K不同应变率下Nb3Sn单晶的拉伸应力-应变曲线
    Figure  9.  Tensile stress-strain curves of Nb3Sn single crystal under different strain rates at 4.2 K

    为了观察不含裂纹和含中心裂纹的Nb3Sn单晶在不同应变率下的裂纹萌生和扩展过程,选取应变率为5×109、1×1010 s−1时的原子结构演化图像进行分析,如图10图11图12图13所示。

    图  10  应变率为5×109 s−1时不含裂纹的Nb3Sn单晶的原子结构演化
    Figure  10.  Atomic structure evolution in Nb3Sn single crystal without crack at strain rate of 5×109 s−1
    图  11  应变率为1×1010 s−1时不含裂纹的Nb3Sn单晶的原子结构演化
    Figure  11.  Atomic structure evolution in Nb3Sn single crystal without crack at strain rate of 1×1010 s−1
    图  12  应变率为5×109 s−1时含中心裂纹的Nb3Sn单晶的原子结构演化
    Figure  12.  Atomic structure evolution in Nb3Sn single crystal with central crack at strain rate of 5×109 s−1
    图  13  应变率为1×1010 s−1时含中心裂纹的Nb3Sn单晶的原子结构演化
    Figure  13.  Atomic structure evolution in Nb3Sn single crystal with central crack at strain rate of 1×1010 s−1

    可以看出,在不同的应变率下,不含裂纹的Nb3Sn单晶在初始阶段均发生弹性变形,随着应变的增加,结构内部出现大量滑移带,位错堆积,之后滑移带上出现孔洞,孔洞交汇形成裂纹,最后Nb3Sn单晶断裂,材料断裂前结构内部出现大量塑性变形,表现出韧性断裂行为。

    对于含中心裂纹的Nb3Sn单晶,在不同的应变率下,其初始阶段均发生弹性变形,随着应变的增加,裂纹变长变宽,裂纹尖端出现位错形核,在裂纹扩展的同时,裂尖发射大量位错,位错堆积,随着塞积数目的变大,塞积群的应力集中使得原子键发生断裂,之后该区域出现孔洞,孔洞变大后与中心裂纹融合,最终使Nb3Sn单晶断裂。对比2.2节的模拟结果可以发现,低应变率(5×108 s−1)下裂纹自身扩展致使材料断裂,而高应变率下裂纹尖端发射的位错是材料断裂的主要原因。实验研究[29-30]指出,若位错在应力场中的运动速度低于裂纹本身的传播速度,就会产生脆性断裂,反之产生韧性断裂。从图12图13可以看出,裂尖发射位错的速度比自身扩展速度快,可见在高应变率下材料发生韧性断裂。

    为了解释不同应变率下含中心裂纹的Nb3Sn单晶的拉伸强度随应变率增加而增大的现象,模拟得到不同应变率下应力峰值处的原子应力云图,如图14所示。随着应变率的增加,裂纹尖端的原子应力集中现象越来越明显,裂纹尖端原子应力越大,裂纹扩展需要克服的阻力也越大,拉伸强度就越大。另外图14表明,在高应变率下,由于裂纹扩展的阻力变大,裂纹扩展速度变慢,因此裂纹尖端发射的位错是材料发生断裂的主要原因。

    图  14  不同应变率下应力峰值处的原子应力云图
    Figure  14.  Atomic stress distribution at stress peak under different strain rates

    工程中应用的Nb3Sn复合超导体具有复杂的多尺度、多级结构特征,大量宏观实验揭示出在极低温区超导体的力学及力-电耦合行为具有多尺度耦合特征[33]。为了更加明确各个尺度上独立的力学及力-电耦合行为,近年来的实验[17]将研究对象调整为单晶体和多晶体。单晶体的结构简单,影响因素明确,有利于揭示背后的力学及力-电耦合机制,缺点是样品制备和加载困难(目前仅有静水压加载实验);而多晶体在制备过程中容易形成Sn原子积聚,产生新的相结构。微观模拟获得的裂纹扩展结果为揭示具有复杂结构特征的Nb3Sn超导体的力学及力-电耦合背后的物理机理提供了一定的参考,也有助于高应变耐受型超导材料的制备和开发。

    采用分子动力学模拟方法,研究了不含裂纹和含中心裂纹的Nb3Sn单晶在拉伸作用下的断裂机制与裂纹扩展行为,分析了应变率效应对Nb3Sn超导体单晶断裂机制与裂纹扩展行为的影响,得到以下结论。

    (1) 4.2 K下,不含裂纹和含中心裂纹的Nb3Sn单晶表现出不同的裂纹扩展行为。不含裂纹的Nb3Sn单晶是结构受力后出现滑移,滑移带上位错塞积导致应力集中,应力集中使原子键断裂,从而萌生裂纹,致使Nb3Sn单晶断裂;而含中心裂纹的Nb3Sn单晶则是裂纹尖端应力集中,使得原子键断裂形成微裂纹,微裂纹扩展致使材料断裂。

    (2)应变率变化对Nb3Sn单晶的弹性变形行为几乎没有影响,随着应变率的增加,Nb3Sn单晶表现出不同的断裂机制,低应变率下表现为脆性断裂,而高应变率下表现为韧性断裂。

    本研究结果为发展实用超导体中裂纹的多尺度扩展模型奠定了一定的基础,同时为工程超导体断裂特征的解释提供参考。

  • 图  Nb3Sn的晶格结构

    Figure  1.  Lattice structure of Nb3Sn

    图  Nb3Sn单晶的分子动力学计算模型

    Figure  2.  Molecular dynamics simulation model of Nb3Sn single crystal

    图  300 K下不含裂纹和含中心裂纹的Nb3Sn单晶拉伸应力-应变曲线

    Figure  3.  Tensile stress-strain curves of Nb3Sn single crystal without crack and with central crack at 300 K

    图  4.2 K下Nb3Sn单晶的拉伸应力-应变曲线

    Figure  4.  Tensile stress-strain curves of Nb3Sn single crystal at 4.2 K

    图  4.2 K下不含裂纹的Nb3Sn单晶的原子结构演化

    Figure  5.  Atomic structure evolution of Nb3Sn single crystal without crack at 4.2 K

    图  4.2 K下含中心裂纹的Nb3Sn单晶的原子结构演化

    Figure  6.  Atomic structure evolution of Nb3Sn single crystal with central crack at 4.2 K

    图  4.2 K下不含裂纹的Nb3Sn单晶在拉伸方向的原子应力云图

    Figure  7.  Atomic stress distribution of Nb3Sn single crystal without crack at 4.2 K

    图  4.2 K下含中心裂纹的Nb3Sn单晶在拉伸方向的原子应力云图

    Figure  8.  Atomic stress distribution of Nb3Sn single crystal with central crack at 4.2 K

    图  4.2 K不同应变率下Nb3Sn单晶的拉伸应力-应变曲线

    Figure  9.  Tensile stress-strain curves of Nb3Sn single crystal under different strain rates at 4.2 K

    图  10  应变率为5×109 s−1时不含裂纹的Nb3Sn单晶的原子结构演化

    Figure  10.  Atomic structure evolution in Nb3Sn single crystal without crack at strain rate of 5×109 s−1

    图  11  应变率为1×1010 s−1时不含裂纹的Nb3Sn单晶的原子结构演化

    Figure  11.  Atomic structure evolution in Nb3Sn single crystal without crack at strain rate of 1×1010 s−1

    图  12  应变率为5×109 s−1时含中心裂纹的Nb3Sn单晶的原子结构演化

    Figure  12.  Atomic structure evolution in Nb3Sn single crystal with central crack at strain rate of 5×109 s−1

    图  13  应变率为1×1010 s−1时含中心裂纹的Nb3Sn单晶的原子结构演化

    Figure  13.  Atomic structure evolution in Nb3Sn single crystal with central crack at strain rate of 1×1010 s−1

    图  14  不同应变率下应力峰值处的原子应力云图

    Figure  14.  Atomic stress distribution at stress peak under different strain rates

    表  1  Nb3Sn单晶的弹性常数和晶格常数

    Table  1.   Elastic constants and lattice constant of Nb3Sn single crystal

    MethodC11/GPaC12/GPaC44/GPaa
    This work284.10 95.8453.765.21
    First principle284.32107.7067.075.32
    下载: 导出CSV
  • [1] 王秋良. 高磁场超导磁体科学 [M]. 北京: 科学出版社, 2008: 24–37.

    WANG Q L. High-field superconducting magnets science [M]. Beijing: Science Press, 2008: 24–37.
    [2] 周又和, 王省哲. ITER超导磁体设计与制备中的若干关键力学问题 [J]. 中国科学(物理学 力学 天文学), 2013, 43(12): 1558–1569. doi: 10.1360/132013-166

    ZHOU Y H, WANG X Z. Review on some key issues related to design and fabrication of superconducting magnets in ITER [J]. Scientia Sinica (Physica, Mechanica & Astronomica), 2013, 43(12): 1558–1569. doi: 10.1360/132013-166
    [3] 梁明, 张平祥, 卢亚锋, 等. 磁体用Nb3Sn超导体研究进展 [J]. 材料导报, 2006, 20(12): 1–4. doi: 10.3321/j.issn:1005-023X.2006.12.001

    LIANG M, ZHANG P X, LU Y F, et al. Advances in Nb3Sn superconductor for magnet application [J]. Materials Review, 2006, 20(12): 1–4. doi: 10.3321/j.issn:1005-023X.2006.12.001
    [4] DEVRED A, BACKBIER I, BESSETTE D, et al. Challenges and status of ITER conductor production [J]. Superconductor Science and Technology, 2014, 27(4): 044001. doi: 10.1088/0953-2048/27/4/044001
    [5] SANABRIA C, LEE P J, STARCH W, et al. Evidence that filament fracture occurs in an ITER toroidal field conductor after cyclic Lorentz force loading in SULTAN [J]. Superconductor Science and Technology, 2012, 25(7): 075007. doi: 10.1088/0953-2048/25/7/075007
    [6] SHETH M K, LEE P J, MCRAE D M, et al. Study of filament cracking under uniaxial repeated loading for ITER TF strands [J]. IEEE Transactions on Applied Superconductivity, 2012, 22(3): 4802504. doi: 10.1109/TASC.2011.2174554
    [7] LU J, HAN K, GODDARD R E, et al. The IC irreversible strain of some ITER high JC Nb3Sn wires [J]. IEEE Transactions on Applied Superconductivity, 2009, 19(3): 2637–2640. doi: 10.1109/TASC.2009.2019597
    [8] NIJHUIS A, MIYOSHI Y, JEWELL M C, et al. Systematic study on filament fracture distribution in ITER Nb3Sn strands [J]. IEEE Transactions on Applied Superconductivity, 2009, 19(3): 2628–2632. doi: 10.1109/TASC.2009.2018082
    [9] CHEGGOUR N, LEE P J, GOODRICH L F, et al. Influence of the heat-treatment conditions, microchemistry, and microstructure on the irreversible strain limit of a selection of Ti-doped internal-tin Nb3Sn ITER wires [J]. Superconductor Science and Technology, 2014, 27(10): 105004. doi: 10.1088/0953-2048/27/10/105004
    [10] DYLLA M T, SCHULTZ S E, JEWELL M C, et al. Fracture strength distribution of individual Nb3Sn filaments [J]. IEEE Transactions on Applied Superconductivity, 2016, 26(8): 6001907. doi: 10.1109/TASC.2016.2602819
    [11] ZHAI Y H, BIRD M D. Florida electro-mechanical cable model of Nb3Sn CICCs for high-field magnet design [J]. Superconductor Science and Technology, 2008, 21(11): 115010. doi: 10.1088/0953-2048/21/11/115010
    [12] WANG X, GAO Y W. Tensile behavior analysis of the Nb3Sn superconducting strand with damage of the filaments [J]. IEEE Transactions on Applied Superconductivity, 2016, 26(4): 6000304. doi: 10.1109/TASC.2015.2509601
    [13] PAPADIMITRIOU I, UTTON C, TSAKIROPOULOS P. Ab initio investigation of the intermetallics in the Nb-Sn binary system [J]. Acta Materialia, 2015, 86: 23–33. doi: 10.1016/j.actamat.2014.12.017
    [14] SUNDARESWARI M, RAMASUBRAMANIAN S, RAJAGOPALAN M. Elastic and thermodynamical properties of A15 Nb3X (X = Al, Ga, In, Sn and Sb) compounds-first principles DFT study [J]. Solid State Communications, 2010, 150(41/42): 2057–2060. doi: 10.1016/j.ssc.2010.08.004
    [15] ZHANG R, GAO P F, WANG X Z, et al. First-principles study on elastic and superconducting properties of Nb3Sn and Nb3Al under hydrostatic pressure [J]. AIP Advances, 2015, 5(10): 1–9. doi: 10.1063/1.4935099
    [16] DE MARZI G, MORICI L, MUZZI L, et al. Strain sensitivity and superconducting properties of Nb3Sn from first principles calculations [J]. Journal of Physics: Condensed Matter, 2013, 25(13): 135702. doi: 10.1088/0953-8984/25/13/135702
    [17] REN Z, GAMPERLE L, FETE A, et al. Evolution of T2 resistivity and superconductivity in Nb3Sn under pressure [J]. Physical Review B, 2017, 95(18): 184503. doi: 10.1103/PhysRevB.95.184503
    [18] 严六明, 朱素华. 分子动力学模拟的理论与实践 [M]. 北京: 科学出版社, 2013: 1−6.

    YAN L M, ZHU S H. Theory and practice of molecular dynamics simulation [M]. Beijing: Science Press, 2013: 1−6.
    [19] CHUDINOV V G, GOGOLIN V P, GOSHCHITSKII B N, et al. Simulation of collision cascades in intermetallic Nb3Sn compounds [J]. Physica Status Solidi (A), 1981, 67(1): 61–67. doi: 10.1002/pssa.2210670103
    [20] 陈伟. 第Ⅰ类超导体中的平衡磁结构及其动力学研究 [D]. 南京: 南京大学, 2016: 1−20.

    CHEN W. Equilibrium intermediate-state patterns and dynamics in a type-Ⅰ superconducting slab [D]. Nanjing: Nanjing University, 2016: 1−20.
    [21] HALL D L. New insights into the limitations on the efficiency and achievable gradients in Nb3Sn SRF cavities [D]. New York: Cornell University, 2017: 1−60.
    [22] ZHANG Y Q, JIANG S Y. Molecular dynamics simulation of crack propagation in nanoscale polycrystal nickel based on different strain rates [J]. Metals, 2017, 7(10): 432. doi: 10.3390/met7100432
    [23] LI L, CHEN H T, FANG Q H, et al. Effects of temperature and strain rate on plastic deformation mechanisms of nanocrystalline high-entropy alloys [J]. Intermetallics, 2020, 120: 106741. doi: 10.1016/j.intermet.2020.106741
    [24] WU W P, YAO Z Z. Influence of a strain rate and temperature on the crack tip stress and microstructure evolution of monocrystalline nickel: a molecular dynamics simulation [J]. Strength of Materials, 2014, 46(2): 164–171. doi: 10.1007/s11223-014-9531-0
    [25] ZHANG Y, ASHCRAFT R, MENDELEV M I, et al. Experimental and molecular dynamics simulation study of structure of liquid and amorphous Ni62Nb38 alloy [J]. The Journal of Chemical Physics, 2016, 145(20): 204505. doi: 10.1063/1.4968212
    [26] KO W S, KIM D H, KWON Y J, et al. Atomistic simulations of pure tin based on a new modified embedded-atom method interatomic potential [J]. Metals, 2018, 8(11): 900. doi: 10.3390/met8110900
    [27] 石震天, 杨绪佳, 王豪阳, 等. 高压下Nb3Sn单晶的超导相转变 [J]. 高压物理学报, 2021, 35(2): 021102. doi: 10.11858/gywlxb.20200615

    SHI Z T, YANG X J, WANG H Y, et al. Superconducting transition of Nb3Sn single crystal under high-pressure [J]. Chinese Journal of High Pressure Physics, 2021, 35(2): 021102. doi: 10.11858/gywlxb.20200615
    [28] DARKINS R, SUSHKO M L, LIU J, et al. Stress in titania nanoparticles: an atomistic study [J]. Physical Chemistry Chemical Physics, 2014, 16(20): 9441–9447. doi: 10.1039/C3CP54357A
    [29] 钟群鹏, 赵子华. 断口学 [M]. 北京: 高等教育出版社, 2006: 54–77, 176−241.

    ZHONG Q P, ZHAO Z H. Fractography [M]. Beijing: Higher Education Press, 2006: 54−77, 176−241.
    [30] 上海交通大学《金属断口分析》编写组. 金属断口分析 [M]. 北京: 国防工业出版社, 1979: 4−26.
    [31] WEST A W, RAWLINGS R D. The microstructure and mechanical properties of Nb3Sn filamentary superconducting composites [J]. Journal of Materials Science, 1979, 14(5): 1179–1186. doi: 10.1007/BF00561303
    [32] TAKEUCHI T, TSUCHIYA K, SAEDA M, et al. Electron backscatter diffraction analysis on Nb3Sn and Nb3Al multifilaments [J]. IEEE Transactions on Applied Superconductivity, 2011, 21(3): 2541–2545. doi: 10.1109/TASC.2010.2083625
    [33] GODEKE A. Performance boundaries in Nb3Sn superconductors [D]. Enschede: University of Twente, 2005: 1−30.
  • 加载中
图(14) / 表(1)
计量
  • 文章访问数:  1334
  • HTML全文浏览量:  441
  • PDF下载量:  30
出版历程
  • 收稿日期:  2021-09-28
  • 修回日期:  2021-12-02
  • 刊出日期:  2022-05-30

目录

/

返回文章
返回