杨育梅 李志鹏

(兰州理工大学理学院,兰州 730050)

引言

YBCO 高温超导带材具有较高的临界转变温度、较高的载流能力及较低的交流损耗,在高磁场下拥有良好的载流能力,因而备受科学界及超导工程技术领域的关注[1],使其成为在液氮温区强电应用最具希望的超导材料,可望在未来超导电缆、超导发电机、超导变压器、超导限流器、STATCOM 等电力设备中得到广泛应用[2-5].然而,超导带材为多层结构,由金属基带、缓冲层、超导层和保护层构成,各材料层性质的差异、复杂的外部场环境及自身的多层结构特性等会造成局部脱黏而极易引起其输电性能及力学性能的退化[6-10].因此,为了推进超导带材的应用,针对其局部脱黏引起的电学及力学性能的研究就显得至关重要.Chen 等[11]提出了基于感应原理对超导带材临界电流进行脉冲测量的方法,并通过实验验证了该方法的可行性.Amemiya 等[12]对RE-123 涂层超导体进行了重复磁场激励实验,研究了其对磁场的稳定性和再现性的影响.而考虑到超导带材所处的复杂外部场环境对其力学性能的影响,Yang 和Wang[13-14]分析了矩形超导薄板仅受载流或仅受外部磁场激励及受承载电流和外部磁场同时作用时,超导薄板内部的电流和磁场的分布情况,给出了薄板内的电流和磁场分布的表述,获得了薄板中的应力、应变及磁致伸缩的解析表述.黑颖顿等[15]对国内外几种典型高温超导带材在77 K 温度、0～1 T/360°环境下的临界电流衰减情况进行了测量,并总结了其对磁场角度依赖性变化的规律.Wang 等[16]对自场下REBCO 超导带材中临界电流密度的变化进行了研究,并分析了在外加磁场、仅载流和外加磁场与传输电流同步施加3 种情况下的电流密度和磁场分布.李腾腾等[17]研究了侧向压力下BSCCO 多芯带材和REBCO 涂层导体载流性能的变化情况,指出,相较于77 K 温度下,BSCCO 在室温下力学性能更强,而REBCO 带材的情形则相反.而随着研究的不断深入,层间剥离及裂纹扩展等引起了超导带材临界电流退化及力学性能丧失,这成为一个不可忽略的问题.鉴于此,Osamura等[18]首次对YBCO 涂层超导带材在室温下超导晶格的变化进行了准确测量,并得到了超导层的内部应变与外加应变之间的变化关系.指出,超导层在热处理后产生-0.2%左右的残余压应变,且该压应变与外加应变呈线性变化直至超导带材损坏.随后,Takematsu 等[19]通过扫描电子显微镜和能谱分析对树脂浸渍过的双饼超导线圈在冷却过程中因残余热应力而引起的剥离分层现象进行了观察,确定了层间剥离发生在超导层与缓冲层之间,同时也发生在超导层内部.Shin 和Gorospe[20]则对YBCO 高温超导带材进行了横向拉伸剥离实验,指出加载面的大小对剥离强度产生直接影响,加载面越大,剥离强度越小,并且对剥离面位置与临界电流退化的关系进行了总结.Shi 等[21]测量了含缺陷的REBCO 超导带在失超过程中的电压和温度曲线,分析对比了两条无缺陷带、一条无缺陷带和一条有缺陷带以及两条带缺陷的叠层REBCO 带的最小失超能(MQE)和正常区沿带材的传播速度(NZPV)的失超特性.指出有缺陷的叠层带的最小失超能小于无缺陷带的最小失超能,有缺陷的叠层带的法向传播速度大于无缺陷带的法向传播速度.为了进一步揭示超导带材的分层机理,刘君等[22]研究了具有不同缺陷高温超导带材中的电磁特性,研究了高温超导带材在磁化过程中的应力分布和磁致伸缩行为,指明,缺陷尺寸和位置会对其磁化及力学性能产生较大影响.Jing 等[23]从理论上分析了存在局部脱黏的超导带材的应力和应变状态,讨论了温度以及承载电流对超导薄膜中的正应力和界面处的剪切应力的影响.而针对局部脱黏后的YBCO 超导带材在外部磁场激励下的力学响应尚未研究.

因此,本文基于超导临界态Bean 模型和弹性力学平面应变方法,考虑局部脱黏后的YBCO 超导带材在外部垂直磁场下的情形,研究了外部磁场上升和下降过程中超导薄膜内电流和磁通密度的分布,并得到了超导体内各应力的具体数值,详细讨论了局部脱黏的存在、磁场大小的变化及基体材质的变化对超导薄膜中正应力及界面处切应力的大小及分布产生的影响.

1 模型建立及理论推导

为简化起见,本研究中考虑的超导带材结构仅包含两层:超导薄膜层和基底层.基于此,考虑一个无限长YBCO 高温超导带材涂层结构受到外部垂直磁场的作用,如图1 所示,薄膜和基底在y方向的总宽度均为2l,厚度分别为tf和ts,杨氏模量分别为Ef和Es(其中f 和s 分别代表薄膜层和基底层).假设薄膜与基底局部发生脱黏,其分离区域的宽度为2ld,(ld=ll,lr).其余部分完美贴合.假设超导薄膜呈各向同性,忽略弯曲刚度,基底材料为非磁性弹性材料.当结构处于外加大小Ba的垂直磁场中时,由于高温超导体俘获磁通线的特性[24],超导薄膜层内将会产生感应电流J和感应磁场B,接下来首先讨论结构内的电流及磁场的分布特征.

图1 超导带材涂层结构横截面示意图Fig.1 Schematic diagram of cross section of superconducting coated conductor tape with local detachment

1.1 超导薄膜中的电流及磁场分布

当结构处于如图1 所示的外部垂直磁场中时,超导薄膜中会产生相应的感应电流以及磁场.随着外加磁场Ba从零开始逐渐增大,磁通线开始从结构的两侧逐渐渗入,当外加磁场未达到完全渗透场时,超导薄膜内的电流密度及磁感应强度的分布如下[25]

1.2 超导体内的磁通钉扎力及应力分布

随着磁通线逐渐渗透到超导薄膜内,在其体内产生相应的磁通钉扎力及基底薄膜界面处的剪切力,如图2(a)所示.基于临界态Bean 模型假定,该钉扎力大小等于洛伦兹力,即:f=J×B.由薄膜中体力f(ζ)与基底薄膜界面处剪切应力g(ζ)在y方向的平衡关系,如图2(b)所示,可以得到超导体内的平衡方程为

图2 超导涂层结构力学等效模型Fig.2 Mechanical equivalent model of superconducting coated conductor tape

其中,cln和crn(n=0,1,2,···,N),为未知系数.考虑到第一类Chebyshev 多项式的性质[29-32],通过联立方程(38)～(39) 可以求得超导薄膜中的正应力及基底薄膜界面处的剪切应力.

2 结果与讨论

为了更直观的描述超导薄膜及基底薄膜界面处的应力分布特征,本节结合以上计算结果进行数值模拟,并讨论了外部磁场上升及下降过程中薄膜内正应力及基底薄膜界面处切应力的分布及随外部磁场的变化规律.薄膜和基体所选取的参数均为典型参数,其非常接近YBCO 超导带材的实际值[23].分别如下:Ef=157 GPa,νf=0.3,Es=200 GPa,νs=0.3,tf=10 µm,l=2 cm,ll=0.4l,ld=0.01l和Jc=3.3 MA/cm2.

2.1 超导薄膜中的正应力分布

图3 给出了当超导薄膜结构在分别受到承载电流作用、外部磁场作用以及承载电流和外部磁场共同作用时,薄膜内的正应力分布.从图3 可以看出,在上述3 种情形下,超导薄膜内的正应力的变化趋势较为相似.由于受到基体材料的约束,当承载电流和外加磁场逐渐增大时,如图3(a)所示,超导薄膜在两端受到不同程度的拉应力作用,而在其他部位受到压应力作用.不同之处在于,当超导薄膜结构在受到承载电流和外加磁场共同作用时,薄膜内的正应力不再是对称分布的,这是由于超导薄膜内的磁通线发生了改变而造成的.需要注意的是,应力奇异性发生在超导薄膜中分离区域的两端,其压应力远大于未分离区域,并且呈线性变化,和未分离区域相比,存在较大差异.而当承载电流和外加磁场从最大值开始下降时,如图3(b)所示,超导薄膜两侧的压应力先增大,而后逐渐减小,直至变为拉应力,在其他部位,又逐渐变为压应力.应力奇异仍然发生在分离部分的左右端.

图3 不同物理场下超导薄膜内的正应力分布Fig.3 Normal stress distribution in superconducting film under different physical fields

图4 描述了超导薄膜中正应力随外部磁场变化的分布规律.当外加磁场从零逐渐上升时,如图4(a)所示,薄膜两侧边缘部分出现拉应力,该拉应力随着外部磁场的增大而显著增大.而在薄膜内部,应力奇异发生在分离区域的末端,该区域压应力值显著增大,意味着该区域薄膜层极易从基底层上脱黏.当外部磁场开始下降时,如图4(b)所示,薄膜两侧边缘出现压应力作用,该压应力随着外部磁场的减小而逐渐减小.随着外部磁场从最大值B0逐渐减小至零,薄膜内的正应力整体呈减小趋势,而当外加磁场进一步从零减小至-B0,薄膜内的正应力整体呈增大趋势.这是因为当外部磁场从最大值下降时,薄膜内俘获的磁通线开始从薄膜边缘排出,导致薄膜内的磁通钉扎力减小,进而引起结构内正应力的逐渐减小.然而,当外部磁场从零持续下降时,磁通线开始反向进入超导薄膜内,此时,俘获在超导薄膜内的磁通线数量开始增多,结构内的磁通钉扎力增大,正应力随之增大.应力奇异现象仍然出现在分离区域的两端,由于超导薄膜中局部磁通密度的分布不同而导致呈现出不同的受力状态.

图4 超导薄膜中的正应力分布Fig.4 Normal stress distribution in superconducting film

图4 超导薄膜中的正应力分布(续)Fig.4 Normal stress distribution in superconducting film(continued)

图5 给出了基底弹性模量对薄膜内正应力的影响规律.无论上升场还是下降场情形,在具有较硬基底的超导薄膜两侧都出现了微小的拉应力作用,而在其他部位,均呈压应力状态.当基底杨氏模量较小甚至为零时,超导薄膜在整体上受压,且随着杨氏模量的减小,压应力逐渐增大.此外,在超导薄膜的分离区域,由于不受基底约束而出现应力奇点.

图5 基底杨氏模量对薄膜内正应力的影响Fig.5 Normal stress distribution in superconducting thin films with different substrates

2.2 基底薄膜界面处的剪切应力分布

图6 描述了基底薄膜界面处切应力随外部磁场变化的分布规律.当外加磁场从零逐渐上升时,如图6(a)所示,在基底薄膜界面两侧及分离区域的末端分别呈现出不同的受力状态,且随着外部磁场的增大而显著增大,这与磁通线从两侧向中心部位渗透的特点有关.当外部磁场开始下降时,如图6(b)所示,基底薄膜界面左侧出现拉应力,并与界面内侧原有的压应力相互作用而形成一个应力尖角,其大小随着外部磁场的减小而逐渐减小.随着外部磁场从最大值B0逐渐减小至零,基底薄膜界面处的切应力整体呈减小趋势,而当外加磁场进一步从零减小至-B0,基底-薄膜界面处的切应力整体呈增大趋势.

图6 基底-薄膜界面处的切应力分布Fig.6 Shear stress distribution on the interface of the substrate

图6 基底-薄膜界面处的切应力分布(续)Fig.6 Shear stress distribution on the interface of the substrate(continued)

图7 给出了基底弹性模量对基底-薄膜界面处的剪切力的影响规律.可以看出,无论是上升场还是下降场情形,基底-薄膜界面处切应力整体呈现出左侧受拉,右侧受压的状态,且对于杨氏模量较大的基底,由于对薄膜具有更强的约束作用,其界面处的切应力呈现出更加剧烈的变化.而当杨氏模量为零时,表示无基底情形,切应力为零.

图7 基底杨氏模量对的超导体界面处的剪切应力分布的影响Fig.7 Shear stress distribution at the interface of superconductors with different substrates

3 结论

本文基于超导临界态Bean 模型,结合弹性力学解决平面应变问题的方法,研究了处于垂直外加磁场下,存在局部脱黏的YBCO 高温超导薄膜中的正应力及基底-薄膜界面处切应力的大小及分布特征.结果显示:

(1)磁场强度的变化对薄膜中的正应力及基底薄膜界面处的切应力大小及分布产生直接影响.对于较高的磁场环境,薄膜承受较大的正应力及界面切应力,且应力极值始逐现在薄膜及基底薄膜界面两侧.因此,薄膜边缘是极易发生脱黏的区域.

(2)在脱黏区域末端会出现应力集中现象,这将加剧带材间裂纹的扩展及层间剥离的发生.

(3)基底材质的变化对薄膜中正应力及基底-薄膜界面处切应力产生显著影响.当基底的杨氏模量较小时,薄膜经受较大的正应力;而对于基底较硬的情形,基底-薄膜界面处的剪切应力在分离区域附近及两侧的变化将更加剧烈.因此,通过综合考虑超导带材所处的外部工作环境、生产工艺及基底材料的性质,可以有效避免超导带材层间剥离的发生,保障其在实际中的安全应用.

附录