Nb3Sn高场复合超导体临界性能力学变形效应的多尺度模拟

A15型Nb₃Sn超导体是制造高场( > 10 T)超导磁体线圈的主要材料, 被广泛应用于磁约束可控核聚变、高能物理等强磁场超导磁体装备制造领域. 力学变形诱导的Nb₃Sn超导临界性能退化给高场超导磁体装备的电磁性能指标和安全稳定运行造成了极其不利的影响. 鉴于Nb₃Sn超导体具有复杂的多尺度结构特征, 不同尺度下变形与超导电性能耦合行为是相互关联的, 本文建立了考虑微/细/宏观关联的非线性力电磁耦合本构模型, 提出了从原子尺度A15晶体结构到超导体微结构到宏观非均质Nb₃Sn复合超导体的多尺度模拟模型. 基于多晶体有限元方法, 对静水压加载条件下Nb₃Sn多晶体超导临界温度衰退和单轴拉压加载条件下Nb₃Sn复合多晶体临界性能衰退行为进行了模拟预测, 预测结果与实验观测结果定性吻合. 该模型揭示了Nb₃Sn复合超导体变形-超导电性能多尺度耦合机理, 实现对高场超导体力、电、磁、热耦合行为的预测, 有助于提高对A15型金属间化合物高场超导复合材料力、电、磁、热多尺度耦合行为的认识和描述能力, 为强磁场超导磁体的设计与制造提供有力的理论支撑.      

超导线圈在强磁场作用下的弹塑性变形分析

根据超导线圈受电磁力、热膨胀、弯曲应力作用的受力机制,首先理论研究了超导线圈在发生弹塑性变形时内部的应力、应变、位移,然后使用两段直线之间由过渡曲线相连接的简化模型作为超导线的应力-应变曲线的近似特性。利用有限元方法,数值模拟了电磁力作用下典型的Nb3Sn超导线材在不同边界条件下的应力、应变、位移分布情况。计算结果定量显示了超导线的弹塑性变形特性同电流密度的关系,初步预测了超导线圈发生塑性变形的区域和扩展以及支撑结构对超导磁体力学行为的影响。     

考虑纤维断裂特征的Nb3Sn股线统计拉伸强度和失效概率

Nb₃Sn 超导复合股线在强磁场工程中有着重要的应用,其拉伸强度为保证和评估其长期服役安全的一个关键参数。本文从 “剪滞”理论出发,利用基于Weibull/Possion统计理论的复合材料纤维断裂的Curtin-Zhou模型,较好地描述复合股线中Nb₃Sn超导纤维碎片化过程,建立了分析超导复合股线统计拉伸强度和失效概率的模型。计算结果表明,青铜法Nb₃Sn超导复合股线的拉伸强度随着初始损伤参数的增大而迅速减小;在4.2K服役温度下,当Weibull模量为8时,随着初始损伤参数的增加,股线的拉伸强度约从900MPa衰减至480MPa,与已有实验结果吻合良好;初始损伤参数约为1时,正规化方差最大。初始损伤和Weibull模量对Nb₃Sn复合股线统计拉伸强度和失效概率函数的分布有着显著影响。     

静水压作用下Nb3Sn多晶体超导临界温度退化的耦合模型

Nb3Sn超导材料主要用于强磁场超导磁体的制造。力学变形诱导的其超导临界性能退化给强磁场超导磁体装置的电磁性能指标和安全运行造成了极其不利的影响。针对静水压作用下,Nb3Sn单晶体和多晶体表现出的不同退化行为,本文基于Maki De Gennes（MDG）关系式,建立了描述Nb3Sn单晶体变形-超导临界温度耦合响应的本构关系,并借助于多晶体有限元方法,对静水压作用下Nb3Sn多晶体超导临界温度退化响应进行了预测,预测结果与实验结果定性吻合。模型实现了从Nb3Sn单晶体到Nb3Sn多晶体变形-超导临界温度退化响应曲线的一致性预测。研究结果有助于提高对Nb3Sn高场超导材料变形-超导电性能耦合行为的认识,为发展描述运行工况下Nb3Sn超导材料力-电磁-热多物理场多尺度耦合行为的建模与数值计算方法提供了一定的基础;同时,相关结果对于特殊工况下高场超导磁体性能的评估和高应变耐受性超导材料的制备也具有一定的指导作用。     

静水压作用下Nb3Sn多晶体超导临界温度退化的耦合模型

Nb3Sn超导材料主要用于强磁场超导磁体的制造。力学变形诱导的其超导临界性能退化给强磁场超导磁体装置的电磁性能指标和安全运行造成了极其不利的影响。针对静水压作用下，Nb3Sn单晶体和多晶体表现出的不同退化行为，本文基于Maki De Gennes（MDG）关系式，建立了描述Nb3Sn单晶体变形-超导临界温度耦合响应的本构关系，并借助于多晶体有限元方法，对静水压作用下Nb3Sn多晶体超导临界温度退化响应进行了预测，预测结果与实验结果定性吻合。模型实现了从Nb3Sn单晶体到Nb3Sn多晶体变形-超导临界温度退化响应曲线的一致性预测。研究结果有助于提高对Nb3Sn高场超导材料变形-超导电性能耦合行为的认识，为发展描述运行工况下Nb3Sn超导材料力-电磁-热多物理场多尺度耦合行为的建模与数值计算方法提供了一定的基础；同时，相关结果对于特殊工况下高场超导磁体性能的评估和高应变耐受性超导材料的制备也具有一定的指导作用。     

基于迭代法的非线性弹性均质化研究

微观结构对复合材料的宏观力学性能具有至关重要的影响,通过合理设计复合材料微观结构可以得到期望的宏观性能.均质化方法作为一种有效的设计方法,它从微观结构的角度出发,利用均匀化的概念,实现了对复合材料宏观力学性能的预测和设计.而当考虑非线性因素,均质化的实现就非常困难.本文利用双渐近展开方法,将位移按照宏观位移和微观位移展开,推导了非线性弹性均质化方程.通过直接迭代法,对非线性弹性均质化方程进行了求解,并给出了具体的迭代方法和实现步骤.本文基于迭代步骤和非线性弹性均质化方程编写MATLAB程序,对3种典型本构关系的周期性多孔材料平面问题进行了计算,对比细致模型的应变能、最大位移和等效泊松比,对程序及迭代方法的准确性进行了验证.之后对一种三元橡胶基复合材料进行多尺度均质化,将其分为芯丝尺度和层间尺度.用线弹性的均质化方法得到了芯丝尺度的等效弹性参数,并将其作为层间尺度的材料参数.在层间尺度应用非线性弹性均质化方法对结构进行计算,得到材料的宏观等效性能,并以实验结果为基准进行评价.     

基于迭代法的非线性弹性均质化研究

微观结构对复合材料的宏观力学性能具有至关重要的影响, 通过合理设计复合材料微观结构可以得到期望的宏观性能. 均质化方法作为一种有效的设计方法, 它从微观结构的角度出发, 利用均匀化的概念, 实现了对复合材料宏观力学性能的预测和设计. 而当考虑非线性因素, 均质化的实现就非常困难. 本文利用双渐近展开方法, 将位移按照宏观位移和微观位移展开, 推导了非线性弹性均质化方程. 通过直接迭代法, 对非线性弹性均质化方程进行了求解, 并给出了具体的迭代方法和实现步骤. 本文基于迭代步骤和非线性弹性均质化方程编写MATLAB 程序, 对3种典型本构关系的周期性多孔材料平面问题进行了计算, 对比细致模型的应变能、最大位移和等效泊松比, 对程序及迭代方法的准确性进行了验证. 之后对一种三元橡胶基复合材料进行多尺度均质化, 将其分为芯丝尺度和层间尺度. 用线弹性的均质化方法得到了芯丝尺度的等效弹性参数, 并将其作为层间尺度的材料参数. 在层间尺度应用非线性弹性均质化方法对结构进行计算, 得到材料的宏观等效性能, 并以实验结果为基准进行评价.      

Nb3Sn单晶超导相转变时的弹性性能变化

Nb3Sn超导体失超，是超导磁体装备运行过程中的重要现象。失超，即超导体从超导相转变为正常相的过程；在强磁场超导磁体工程中，由于超高的储能密度，失超伴随着力/热/电等物理参量在瞬时的剧烈变化。失超瞬时，超导相转变的同时伴随着弹性力学性能的突变，研究超导相转变时弹性性能的变化是失超诱发应力跨尺度分析的关键。本文首先采用第一性原理方法计算了弹性常数随温度的变化规律，结果表明由于未考虑A15结构的超导材料在环境温度变化的情况下产生的特殊电子能带结构，基于准静态近似方法的材料弹性常数计算从0K外推至有限温度时，会导致模拟结果与实验观测结果出现定性上的差异；之后，基于晶格自由能函数，给出了描述立方相Nb3Sn单晶弹性性能随温度变化的解析模型，模型预测结果与实验观测结果定性吻合，初步实现了对Nb3Sn单晶超导相转变时弹性性能变化的理论描述和预测。研究结果对于超导体失超应力的跨尺度模拟及超导磁体的安全分析具有一定的理论参考价值。     

考虑面层约束时蜂窝芯弹性常数的确定

本文研究了面板约束效应对蜂窝芯弹性性能的影响并基于六边形蜂窝芯夹层板的代表性体单元建立了考虑面板约束效应的蜂窝芯位移模式.根据真实蜂窝芯代表性单元与其等效材料的代表性单元在宏观均匀应变作用下的总弹性应变能相等,导出了蜂窝芯材料等效弹性常数的公式,进一步探讨了蜂窝几何参数对其弹性常数的影响规律.通过与现有理论结果、有限元结果和实验数据进行对比证实本文导出的弹性常数公式具有较好的预报精度.结果表明面板约束对蜂窝芯等效弹性常数有显著影响.     

低温超导核磁共振陀螺仪模型

建立了单工作物质的三自由度3He低温超导核磁共振陀螺仪结构;利用了量子力学和经典动力学,经过严密的力学分析和数学演算,给出了三轴陀螺仪的工作原理和结构示意图;针对陀螺仪的交叉轴角速率耦合问题,给出了附加磁场线圈解耦法和冗余设计结构,最后建立了基于超导量子干涉仪探测磁矩,并采用最小二乘估计法来推导陀螺进动频率.陀螺仪测速范围可以达到10-9～103 rad/s,漂移为10-4 (°)/h.该陀螺仪结合低温超导技术具有高精度的前景.需要进一步对超导量子干涉磁矩检测仪的精度与陀螺性能进行研究.     

Electric-magnetic-force characteristics of rare earth barium copper oxide superconductor high-field coils based on screening effect and strain sensitivity

Rare earth barium copper oxide (REBCO) is the most researched and commercialized second-generation high-temperature superconducting material. Due to the anisotropic structure, strong deformation sensitivity, and central field errors caused by screening current effects, it is still a challenge for commercialization applications. In this study, the transversely isotropic constitutive relationship is selected as the mechanical model based on the structural characteristics of REBCO tapes, and suitable microelements are selected to equate the elastic constants using their average stress-strain relationships. Then, a two-dimensional axisymmetric model for coils wound by single-layer tapes is constructed to analyze the dependence of the electric-magnetic-force distribution in the tape on the strain. Finally, the anisotropic approximation of the homogenized bulk method is used to equate large-turn high-field coils, and the electric-magnetic-force distribution characteristics of the coils with/without screening effects and mechanical strain conditions are investigated, respectively. The results reveal that the mechanical strain has a weakening effect on the electromagnetic field distribution of superconducting tapes, but causes a significant enhancement in the force field distribution. In the presence of 0.5% mechanical strain, the maximum weakening of the peak value of the current density and the critical current density inside the high-field coil can reach about 8% and 13%, respectively, with a nearly 5 times increase in the peak stress. The screening current makes the current field distribution inside the coil improve by about 10 times. The screening current induced magnetic field can reach up to 0.8 T, making the relative error of the high-field coil center up to 7.8%.     

Analytic calculation of magnetic force between two current-carrying coils

Current-carrying coils are basic elements in electromagnetic equipments, for example, in high field magnets from high temperature superconducting wires or tapes. In the assembly of these systems and their current-carrying operation, unavoidable mis-alignment and shift from the original position can be induced by disturbances such as the imbalance of magnetic force due to safety problems. For two current-carrying coils with non-coplanar axes, the analytic expression of the magnetic force between the two coils is presented according to the rule of Ampere circulation and the Biot-Savart law. Based on the expression, the dependence of the magnetic force on the size and the relative position of each other is further investigated, and the variation of the magnetic force is obtained with the above parameters.     

Quench characteristics and mechanical responses during quench propagation in rare earth barium copper oxide pancake coils

Quench and mechanical behaviors are critical issues in high temperature superconducting(HTS)coils.In this paper,the quench characteristics in the rare earth barium copper oxide(REBCO)pancake coil at 4.2K are analyzed,and a two-dimensional(2D)axisymmetric electro-magneto-thermal model is presented.The effects of the constituent materials,background field,and coil size are analyzed.An elastoplastic mechanical model is used to study the corresponding mechanical responses during the quench propagation.The variations of the temperature and strain in superconducting layers are compared.The results indicate that the radial strain evolutions can reflect the transverse quench propagation and the tensile hoop and radial stresses in superconducting layers increase with the quench propagation.The possible damages are discussed with the consideration of the effects of the background field and coil size.It is concluded that the high background field significantly increases the maximum tensile hoop and radial stresses in quenching coils and local damage may be caused.     

超导磁体等效材料参数的有限元预测

HT—7U超导托卡马克装置中纵场超导磁体是由多种材料组成的且具有周期性分布的大型复杂的重要结构，对其宏观等效材料参数的准确掌握在其结构设计阶段是十分重要的。本文着重介绍了利用小参数展开的渐近均匀化方法对其进行等效处理的过程，并采用有限元方法预测出其等效弹性模量与热膨胀系数等，这些都为纵场磁体结构的力学计算与传热分析提供了必要的参数。     

超导电磁悬浮力的数值模拟

利用超导体宏观电磁场的Kim模型对高温超导体－永磁铁悬浮系统中二才之间的电磁悬浮力进行了定量计算与模拟。通过与已有实验结果比较表明：对超导悬浮力的数值模拟能较好地反映悬浮力的主要特征。在此基础上，给出了超导体内屏蔽电流，铅直悬浮力和面内电磁力的分布情况，数值仿真结果显示了悬浮力的大小与永磁铁的表磁强度和尺寸有关。这些结果为提高超导悬浮力和改进超导宏观电磁场模型提供了理论依据。     

耦合运动多级感应线圈炮三维机电耦合模型

对于多级同步感应线圈炮来说,在以前的文献中比较成熟的模型是基于电流丝方法建立起来的机电模型,它的优点计算速度快,缺点无法考虑线圈炮的外围结构和材料特性对其产生的影响,针对这种情况,推导给出了总体采用矢量磁位 、电路中的驱动线圈中电流 作为求解对象的场路耦合瞬态电磁场模型,利用节点力法计算电枢所受的电磁力,根据电枢运动学模型,进一步推导出耦合运动的多级感应线圈炮的机电模型,并进行有限元离散求解。应用上述模型对5级感应线圈炮进行分析,计算结果与实验结果相吻合,验证了该模型的正确性。     

HT-7U超导托卡马克装置极向场线圈的等效弹性常数及其残余刚度

HT－7装置的超导极向场线圈由导管，超导材料和绝缘体组合而成，具有复杂的正交各向异性的性质。本文从微观力学的观点出发，将极向场线圈整体作为一种复合材料来考虑，采用微观力学的复合律和损伤力学的细观损伤模型，求出了极向场线圈的等效弹性常数和残余刚度；并且，残余刚度的计算结果已得到了实验验证。这些材料特性已被用于有限元分析的输入数据。     

Transient multi-physics behavior of an insert high temperature superconducting no-insulation coil in hybrid superconducting magnets with inductive coupling

A transient multi-physics model incorporated with an electromagneto-thermomechanical coupling is developed to capture the multi-field behavior of a single-pancake (SP) insert no-insulation (NI) coil in a hybrid magnet during the charging and discharging processes. The coupled problem is resolved by means of the finite element method (FEM) for the magneto-thermo-elastic behaviors and the Runge-Kutta method for the transient responses of the electrical circuits of the hybrid superconducting magnet system. The results reveal that the transient multi-physics responses of the insert NI coil primarily depend on the charging/discharging procedure of the hybrid magnet. Moreover, a reverse azimuthal current and a compressive hoop stress are induced in the insert NI coil during the charging process, while a forward azimuthal current and a tensile hoop stress are observed during the discharging process. The induced voltages in the insert NI coil can drive the currents flowing across the radial turns where the contact resistance exists. Therefore, it brings forth significant Joule heat, causing a temperature rise and a uniform distribution of this heat in the coil turns. Accordingly, a thermally/mechanically unstable or quenching event may be encountered when a high operating current is flowing in the insert NI coil. It is numerically predicted that a quick charging will induce a compressive hoop stress which may bring a risk of buckling instability in the coil, while a discharging will not. The simulations provide an insight of hybrid superconducting magnets under transient start-up or shutdown phases which are inevitably encountered in practical applications.     

Buckling behavior of a superconducting magnet coil

The buckling behavior of a superconducting magnet coil has been investigated for out of plane perturbations. It has been shown that such a coil becomes unstable due to its own magnetic field for out of plane bending at the operating values of electric currents. A technique devised by the author that computes the inductance change due to deformation, has been incorporated in this study. The evaluation shows that although the magnetic fields external to the coil are important and that there is additional structural stiffness provided by the surrounding structure, the effects of the self magnetic field should be evaluated and incorporated in the design of such coils. However, based on previous studies, both analytical and experimental, it has been established that the out-of-plane deformation is the most critical mode from a stability standpoint. Therefore, this mode has been exclusively addressed in this study.