专利名称：基于应变的管内电缆导体结构模拟设计方法
技术领域：
本发明涉及管内电缆导体，尤其涉及一种基于应变的管内电缆导体结构模拟设计 方法。
背景技术：
面对当今的能源危机和短缺，聚变装置产生的核能是解决它的理想办法，因为聚 变能既清洁又巨大。但聚变能的实现需要激发产生和控制高温等离子体，目前国内外通用 的是磁约束方法，常规磁体由于焦耳热等问题难以通过快变大电流，现在采用超导磁体技 术。管内电缆导体(Cable-in-ConduitConductor，CICC)因具有良好冷却、高电压绝缘、 大电流、低损耗、多股绞缆等的优点，而成为超导磁体的首选导体，如中国大科学工程项目 EAST、韩国的KSTAR以及国际热核聚变反应堆ITER等都采用CICC导体。CICC导体是由超 导组分、稳定体、冷却介质和Jacket构成的复合体，并由超导股线或铜股线经多级扭绞而 成的，迫流氦在它们的空隙中流动，其间的热交换是一个复杂的过程，当导体受到扰动时， 就会有正常区产生，甚至可能会失超，因此，CICC导体的稳定性受到多因素的制约，合理设 计CICC导体是保证其在快变磁体中稳定安全运行的前提和关键。因为受电磁参数和几何参数的影响，CICC导体的设计理论非常复杂，工程实际中 CICC导体的设计和研制是多次反复尝试的烦琐过程，需要花费大量人力、物力和财力，如中 国的EAST上CICC导体的设计就用了多年时间，花费了大量研究人员的精力，但实际中又需 要大量的CICC导体，为了快速和准确地设计CICC导体，需要开展CICC导体数值模拟设计 研究，开发有效的软件设计支持系统，以帮助加快工程上CICC导体的设计进度。目前针对 CICC导体的模拟研究工具主要有Gandalf、CID以及MFEM1D。Gandalf主要是针对已经存 在的导体结构进行失超模拟分析，缺少导体结构设计理念；CID的仿真设计是通过机械式 的复制超导股线以及不同层的绞缆和压缩过程来堆积形成CICC的几何形状；国内二维失 超分析程序MFEMlD用于NMR磁体的设计分析，它通过建立自适应移动网格有限元模型来分 析CICC导体结构。这几种设计方法对磁场变化、电流密度，特别是应变对稳定性干扰因素 的影响没有加以充分考虑。由于Nb3Sn等A15材料构成的CICC导体，若应变的干扰作用得 不到遏制，加上应变对η值(η值是超导体从超导态向正常态转变时退钉扎的同步率)产生 的作用效果，会有正常区的扩散和传播，进而就会出现CICC导体的失超现象，从而影响到 CICC导体的稳定性，显然应变因素制约和影响着CICC导体的稳定运行。目前在导体的数值 模拟设计上，也多是基于稳定性、铜超比优化的空间电流密度理论等，加之事后的交流(AC) 损耗来验证模拟获得相对合理的导体结构形式，没有对Nb3Sn等的应变效果进行考虑。同 时从现有研究工作和文献中看，缺乏针对CICC导体进行数值模拟设计的研究，需要研制基 于应变的管内电缆导体数值模拟设计方法。
发明内容
本发明的目的是提供一种基于应变的管内电缆导体结构模拟设计方法，模拟设计出导体的结构，保证导体的稳定性。本发明采用下述技术方案一种基于应变的管内电缆导体结构模拟设计方法，包 括步骤(1)、获得导体的应变数据具体是假设作用在导体上的机械应变全部传给导体 的丝线，由实验测试数据可以得到导体的有效应变ε = e。p+eth+eextra，其中，ε。ρ为导体 纵向应变，￡th是零电流时的热应变，是从实验数据中估算出来的Sraitea为其它额外原因导 致的应变；(2)、获得在导体应变影响下的临界电流密度具体是根据实验获得的体积钉扎 力公式f=|jxb|，结合给定条件的磁场强度β和应变势能函数，以及用cicc导体焓值的 估计公式，可获得临界电流密度Jc(B,T, ε) =C(e) (Bc2 (Τ, ε )) (l_t2) 2b_"2 (l_b)2其中Jc (B，Τ, ε )是给定运行条件下的临界电流密度，B是运行时的背景场强，T 是 运行温度，ε是股线中Nb3Sn超导丝在纵向上的总应变；C(0 =C0(l-a0 ε |17)1/2;Bc2 (Τ, ε ) = Bc20 ( ε ) (1-tT) (l_t/3)；
, BTο =- .t =-；
BciiT, ε) Tco{s)T⑶(O =TC0M(l-a0| ε |L7)1/3;BC20(O = BC20M(l-aJ ε Γ)；(3)、假设条件；(4)、获得稳定性裕度、获得铜超比优化的空间电流密度、获得温度裕度及临界电 流，从而得出关于导体铜超比、股线直径、股线根数的方程；(5)、求解并进行导体结构数值模拟设计。所述的步骤⑵中临界电流密度在对于国际热核聚变反应堆ITER上的股线，当磁 场B为零时，临界电流密度为Jcc (B, Τ, ε ) = (1/α(Β，Τ，ε)+Jc0 (B, Τ, O))；Jc0 (B, Τ, ε ) = Jc00 (B, Τ, ε ) (l_(t)2)2 ；其中J。qq(B，Τ, ε )是根据一定条件由实验获得的经验值，Jc0(B, Τ，ε )，在低场中 起矫正作用，在高场中J。o(B，T，ε)可以忽略不计。本发明基于应变的管内电缆导体结构模拟设计方法在实验实测的应变数据基础 上，构建应变对导体临界电流密度的数学关系，根据假设条件以及模拟运行条件，获得CICC 导体模拟设计的模型，形成新的导体模拟设计方法，为合理设计股线直径，选取超导体、铜 组分以及获得Cicc导体结构提供了可靠依据，保证了导体的稳定性；本设计方法能快速和 准确地设计CICC导体，通过数值模拟为ITER等超导托卡马克装置上磁体系统中CICC导体 的稳定性研究、工程设计以及探讨应变作用机理提供方便可靠工具；并且还有助于推动国 内聚变能技术的快速发展、使CCIC导体的模拟设计研究进入国际前沿领域，并处于领先地 位，这具有重要意义。
图1为本发明中的方法流程图2为本发明模拟设计得到的导体结构示意图。
具体实施例方式本发明基于应变的管内电缆导体模拟设计方法主要针对Nb3Sn等A15材料构成的 CICC导体的应变对稳定性的影响进行研究，而对Nb3Sn等A15材料构成的CICC导体，其的 稳定性不仅受到温度裕度、能量裕度、磁场梯度以及交流(AC)损耗等的作用，还受应变对 临界电流密度的影响。本发明研究了应变对导体的稳定性影响，获得股线应变对临界电流 密度等参数影响的量化结果，对Nb3Sn等材料的应变与临界电流密度以及η值的关系进行 研究，根据应变对CICC导体临界电流密度影响的机理和经验数据以及某些假设条件，在建 立关于铜超比、股线直径、股线根数的矩 阵方程基础上，进行数值模拟获得CICC导体的结 构，如图1所示，具体方法步骤如下所述(1)获得导体的应变数据，具体过程为假设作用在CICC导管上的机械应变全部 传给导体的丝线，由实验测试数据可以得到Nb3Sn的有效应变ε ε = ε。ρ+ε th+ε extra (1)其中^㈤为导体纵向应变(也称为运行时应变)，￡th是零电流时的热应变，是从 实验数据中估算出来的，^xtra为其它额外原因导致的应变；严格意义上^㈤和Srartra是导 体真正的纵向应变，实际设计中，ε extra与η值是影响导体结构的重要参数，其变化与导体 以及股线的性能有关；(2)获得导体应变影响下的临界电流密度，具体是根据实验获得的体积钉扎力 公式F= I JXBl，结合给定条件的磁场强度B和应变势能函数(它与Nb3Sn等超导材料的 电子、强子耦合系数相关)，以及用CICC导体焓值的估计公式，可获得要研究的Nb3Sn临界 电流密度Jc (B, Τ, ε ) = C ( ε ) (B C2 (Τ，ε)) _1/2 (l_t2) 2b_"2 (l_b)2 (2)其中J。(B，T，ε )是给定运行条件下的临界电流密度，B是运行时的背景场强，T是 运行温度，ε是股线中Nb3Sn超导丝在纵向上的总应变；其中C(0 =C0(l-aJ ε L7)1/2; (3)Bc2 (Τ, ε ) = Bc20 ( ε ) (1-tT) (l_t/3) ； (4) , B = ⑷T⑶(O =Tcom (Ι—α。I ε |L7)1/3; (7)BC20(e) =BC20M(l-aJ ε |L7) ；(8)其中a Q的取值情况如下当ε < 0时(Nb3Sn受压时)，a Q = 900 ；当ε > 0时 (Nb3Sn受拉时)，a Q = 1250。CQ、Tcom和Bc2qm是根据不同条件由实验获得的经验初始值；对于ITER上的股线，当磁场B为零时，临界电流密度为Jcc (B, Τ, ε ) = (1/(JC(B,T, ε)+Jc0 (B, Τ, O)) ； (9)Jc0 (B, Τ, ε) =Jc00 (B, Τ, e)(l-(t)2)2；(10)其中Jc00 (B, Τ, ε )是由实验获得的经验值;Jc0 (B，Τ, ε )在低场中起矫正作用；在高场中J。o(B，Τ, ε )可以忽略不计；(3)假定条件假定条件为(a)、一级子缆采用3根超导股线；(b)、为获得合理的 电流密度，设定导体运行在过度区间，并取Stekly参数小于1 ； (c)、若有纯铜股线则可以处 理成三种情况在计算起稳定作用的铜面积和湿边周长时都包括纯铜股线的作用；仅在计 算起稳定作用的铜面积时考虑铜股线的作用；仅在计算湿边周长时考虑铜股线的作用；(4)稳定性裕度方面公式根据导体运行在过度区，由稳定性裕度可获得Stekly 参数可表示为
其中fhe为CICC导体的空隙率；Δ Etr为过渡区的稳定性裕度；P和Che分 别为液氦的密度和液氦的热容；Τ。和τ。ρ分别是超导体的临界温度和运行温度；由 于Stekly参数α是表示导体产生的焦耳热与该热量和液氦之间热量传递比率，可得 其中ρ。u为铜的电阻率；PW为导体的湿边周长；A。u为铜的面积；h为氦的热传递系 数，它是一个变化量；结合CICC导体的湿边周长、超导组分和铜组分的面积等，对式(12)进行变换，结 合式(11)并利用铜的面积和湿边周长乘积的表达式可获得第一个关于铜超比(Rcu)、股线 直径(dsc)、股线根数(Nsc)的方程 其中乂=2jr .；； τ、(14) π KPah(Tc-Top)Rcu表示铜超比，ds。、Ns。分别代表CICC导体中超导股线的直径和根数，Kp是计算湿 边周长因子。考虑到CICC导体结构，对三角股的湿边周长一般取为5/6 ；(5)铜超比优化及空间电流密度方面公式下、上空间电流密度公式如下 其中J。为由步骤(2)中获得的临界电流密度，COS θ为扭转角；由此获得第二个关于铜超比(RJ、股线直径(ds。)、股线根数(Ns。)的方程 其中 (6)温度裕度及临界电流方面公式
温度裕度 临界电流 由此获得第三个关于铜超比(Reu)、股线直径(ds。)、股线根数(Nse)的方程 其中I。p为运行中电流；(7)求解并进行导体结构数值模拟设计求解过程为在步骤⑷、(5)、(6)分别 获得三个包含铜超比R。u、股线直径ds。、股线根数Ns。未知量的矩阵方程即公式(13)、(17)、 (21)，但当把公式(21)代入到公式(13)时会得到类似于公式(17)的方程，这说明在空间 电流密度上对CICC导体铜组分的优化不是孤立的问题，是与稳定性裕度、温度裕度及临界 电流有着必然的联系，于是把股线直径ds。表示成股线根数Ns。的函数 在数值模拟设计过程中，首先需要求解股线根数Ns。，它的获得是根据ITER等装 置工程的实际需要初步确定CICC导体的层数，然后采用逐步尝试法，由一、二、三、四等层 的结构来循环求出超导股数Ns。，模拟程序同时把导体的层数以及不同层中股线数目记录 在数据库中，还包括运行电流、温度等数据参数，再由股线直径与股线根数的关系求出符 合实际要求的股线直径ds。，最后求出铜超比R。u，这样由数值模拟设计过程中获得的超导 股数Ns。、股线直径ds。以及铜超比1^数据就可以得到CICC导体结构。如图2所示，其中 (3SC+0Cu)X3X3X3X6表示的是第一级子缆由3根超导股线构成，第二级子缆由3根一 级子缆构成，第三级子缆由3根二级子缆构成，第四级子缆由3根三级子缆构成，CICC导体 (第五级缆)由6根四级子缆构成。本发明所述的方法简化了导体结构的设计过程，具有深 远的意义。
权利要求
一种基于应变的管内电缆导体结构模拟设计方法，其特征在于包括步骤(1)、获得导体的应变数据具体是假设作用在导体上的机械应变全部传给导体的丝线，由实验测试数据可以得到导体的有效应变ε＝εop+εth+εextra，其中，εop为导体纵向应变，εth是零电流时的热应变，是从实验数据中估算出来的εextra为其它额外原因导致的应变；(2)、获得在导体应变影响下的临界电流密度具体是根据实验获得的体积钉扎力公式F＝|J×B|，结合给定条件的磁场强度B和应变势能函数，以及用CICC导体焓值的估计公式，可获得临界电流密度Jc(B，T，ε)＝C(ε)(B C2(T，ε))-1/2(1-t2)2b-1/2(1-b)2其中Jc(B，T，ε)是给定运行条件下的临界电流密度，B是运行时的背景场强，T是运行温度，ε是股线中Nb3Sn超导丝在纵向上的总应变；C(ε)＝C0(1-a0|ε|1.7)1/2；BC2(T，ε)＝BC20(ε)(1-tT)(1-t/3)； <mrow><mi>b</mi><mo>=</mo><mfrac> <mi>B</mi> <mrow><msub> <mi>B</mi> <mrow><mi>C</mi><mn>2</mn> </mrow></msub><mrow> <mo>(</mo> <mi>T</mi> <mo>,</mo> <mi>ϵ</mi> <mo>)</mo></mrow> </mrow></mfrac><mo>;</mo> </mrow> <mrow><mi>t</mi><mo>=</mo><mfrac> <mi>T</mi> <mrow><msub> <mi>T</mi> <mrow><mi>C</mi><mn>0</mn> </mrow></msub><mrow> <mo>(</mo> <mi>ϵ</mi> <mo>)</mo></mrow> </mrow></mfrac><mo>;</mo> </mrow>TC0(ε)＝TCOM(1-α0|ε|1.7)1/3；BC20(ε)＝BC20M(1-α0|ε|1.7)；(3)、假设条件；(4)、获得稳定性裕度、获得铜超比优化的空间电流密度、获得温度裕度及临界电流，从而得出关于导体铜超比、股线直径、股线根数的方程；(5)、求解并进行导体结构数值模拟设计。
2.根据权利要求1所述的基于应变的管内电缆导体结构模拟设计方法，其特征在于 所述的步骤(2)中临界电流密度在对于国际热核聚变反应堆ITER上的股线，当磁场B为零 时，临界电流密度为Jcc(B, Τ, ε ) = (1/(JC(B,T, ε )+Jc0(B, Τ, ε )))；jc0 (β, τ, ε) = jc00 (β, τ, οα-α)2)2;其中Jcfltl(B，Τ, ε )是根据一定条件由实验获得的经验值，Jc0(B, Τ，ε )，在低场中起矫 正作用，在高场中J。。(B，Τ, ε )可以忽略不计。
全文摘要
本发明公开了一种基于应变的管内电缆导体结构设计方法，包括步骤(1)获得导体的应变；(2)获得在导体应变影响下的临界电流密度(3)假设条件；(4)获得稳定性裕度、获得铜超比优化的空间电流密度、获得温度裕度及临界电流，从而得出关于导体铜超比、股线直径、股线根数的方程；(5)求解并进行导体结构数值模拟设计。本方法为合理设计股线直径，选取超导体、铜组分以及获得导体结构提供了可靠依据，保证了导体的稳定性；能快速和准确地设计导体，通过数值模拟为ITER等超导托卡马克装置上磁体系统中导体的稳定性研究、工程设计以及探讨应变作用机理提供方便可靠工具；并且还对推动国内聚变能技术的快速发展有重要意义。
文档编号G06F17/50GK101840448SQ20101013839
公开日2010年9月22日 申请日期2010年4月2日 优先权日2010年4月2日
发明者史卫亚, 李战升, 杨凯, 王雪涛, 王高平, 蒋华伟, 赵玉娟, 马浩歌, 魏欣 申请人:河南工业大学