35 kJ YBCO单螺管型超导储能磁体的多目标优化

35 kJ YBCO单螺管型超导储能磁体的多目标优化
李泰来;杨徉;杨正
【摘要】优化设计超导磁体,不仅可以从技术上保证超导磁储能系统运行的安全和可靠性,而且能够最大限度地降低制造成本.利用模拟退火算法优化在一定条件下的单螺管储能磁体的基本参数(长、外径、内径等),得出最优设计方案.在ANSYS中对超导储能磁体进行模拟仿真实验,结果验证了优化设计的优越性.
【期刊名称】《吉首大学学报（自然科学版）》
【年(卷),期】2019(040)004
【总页数】4页(P48-51)
【关键词】高温超导磁体;多目标优化;ANSYS
【作者】李泰来;杨徉;杨正
【作者单位】吉首大学信息科学与工程学院,湖南吉首416000;吉首大学信息科学与工程学院,湖南吉首416000;吉首大学信息科学与工程学院,湖南吉首416000【正文语种】中文
【中图分类】TM153+.1
超导磁体所储存的电磁能量等效于磁体在自由空间中产生的磁场能量.它由2个部分组成：一部分位于磁体绕组所包围的有限空间中，磁场强度强，磁能密度较高；另一部分位于磁体绕组的外部，磁通密度较弱但分布最广，形成超导储能磁体的漏磁场区.超导储能磁体，特别是微型储能磁体，在系统中常运行在脉动状态下.超导
材料在交变磁场中承受交流损耗，从而增大制冷费用.相同储能量级下，超导材料
用量越少，交流损耗也越少，超导磁体的运行越经济.因此，超导储能磁体的材料
利用率是磁体设计中应考虑的一个重要因素[1].近年来，学者们对磁体优化的研究[1-7]表明，磁体优化设计的主要思路是针对受磁场影响而导致临界电流降低的部分，采取一些措施优化磁体的性能，从而达到优化的目的.笔者拟通过改变磁体形
状因数来优化磁体结构，并通过有限元软件对优化后的单螺管磁体进行建模仿真，验证其优越性.
1 高温超导磁体的优化设计
图1 超导磁体线圈截面Fig. 1 Cross Section of Superconducting Magnet Coil 对于未设置铁磁屏蔽的超导磁体，可以构造如图1所示的模型进行优化设计.该超
导磁体的储能表达式为
(1)
其中：l(p，q) 为点p与点q之间的距离；J(p)和J(q)均为电流密度；S为磁体线
圈横截面的面积.由(1)式可知，磁体的体积和形状参数直接影响磁体储能量.因此，将通常使用的如图2所示的矩形截面磁体线圈模型设为参照方案，优化设计阶梯
形截面磁体线圈模型如图3所示.
图2 矩形截面磁体线圈模型Fig. 2 Model of Magnet Coil with Rectangular Cross Section
图3 阶梯形截面磁体线圈模型Fig. 3 Model of Magnet Coil with Ladder Cross Section
设计材料为YBCO超导磁体线圈，储能量35 kJ.储能磁体为薄壁单螺管双饼结构，单饼厚度4.5 mm，单饼间用聚酰亚胺绝缘，双饼间用铝铜线隔开[5].用线圈体积
作为目标函数.矩形截面磁体线圈(方案1)的体积目标函数为
f=2π((0.065+x1)2-0.0652)h1；
阶梯形截面磁体线圈(方案2)的体积目标函数为
f=2π(((0.065+x2)2-0.0652)h2+((0.065+x2+x3)2-(0.065+x2)2)h3+
((0.065+x2+x3+x4)2-(0.065+x2+x3)2)h4)，
其中各量在图2和图3中标出.以此目标进行优化，在保证磁体储能量的同时，使
体积(即超导带材用量)最小.
对磁体的优化还应同时满足如下约束条件：(1)线圈储能误差不超过±0.2；(2)线圈的内径为65 mm；(3)符合YBCO导体的B-I特性曲线.约束条件的数学表达式为：
其中：E为磁体的储能量；R为固定的线圈内径；I0p为磁体的饱和电流.
按照2种方案的相关约束条件设置参数，使用模拟退火算法对磁体目标函数中的
几何参数寻求最优解.所得优化结果见表1和表2.
表1 方案1形状因数的优化结果Table 1 Optimized Results of Form Factor in Scheme 1 mmx1x2x3x4h1h2h3h4313———1 058———
表2 方案2形状因数的优化结果Table 2 Optimized Results of Form Factor in Scheme 2 mmx1x2x3x4h1h2h3h4—11876266—982654428
根据磁体模型长、宽、高等参数计算出2个方案中磁体的体积和自感系数，并求
得磁体电流和对应的储能量，如表3所示.
表3 磁体其他参数的计算结果Table 3 Calculated Results of Other Parameters 方案体积/m3电流/A 储能量/J11.678 0×10-3210.03 499.8021.466 5×10-3213.83 501.25
从表1～3中可以看出，在选取YBCO超导材料、磁体储能量约为35 kJ的情况下，方案1的单螺管型需要使用带材1.678×10-3 m3；方案2需要使用带材1.466
5×10-3 m3，相比前者可节省13%的超导带材，使得超导磁体设计制作的成本大
大减少，从而增加超导磁储能系统的经济效益.
2 ANSYS仿真分析
ANSYS是目前应用最为广泛的有限元分析软件之一.利用ANSYS对储能磁体进行仿真分析，能够清晰、直观地得出2种方案磁体周围的磁感应强度分布情况.根据表2和表3中的各项磁体参数，可进一步计算得出2个方案在ANSYS建模中的参数设置，列于表4和表5.
表4 方案1的参数设置Table 4 Parameter-Setting of Scheme 1绕组内径/mm 绕组外径/mm绕组高度/mm双饼数电流/A6596.3105.812210.0
表5 方案2的参数设置Table 5 Parameter-Setting of Scheme 2绕组外径
r1/mm绕组外径r2/mm绕组外径r3/mm绕组高度r4/mm绕组高度r5/mm绕组高度r6/mm双饼数电流/A76.884.4111.0105.865.442.812213.8
建立有限元模型并进行网格划分，再利用SOLVE命令求解.在求解过程中，按照规定设置磁体外2 m处的空气单元符合第二类狄利克雷边界条件(磁力线平行)，并对磁体作电流密度加载，然后利用MAGSOLV命令对所建立的磁体模型进行求解，得到2个方案的磁感应强度分布，结果如图4和图5所示.
图4 矩形截面磁体磁感应强度分布云图Fig. 4 Cloud Chart of Magnetic Induction Intensity of Rectangular Cross Section
图5 阶梯形截面磁体磁感应强度分布云图Fig. 5 Cloud Chart of Magnetic Induction Intensity of Ladder Cross Section
整理2个方案的仿真结果，列于表6.
表6 2种优化方案结果比较Table 6 Simulation Comparison of Two Optimized Schemes方案电流/A储能量/kJ中心磁感应强度/T2 m外杂散磁场/T1210.035.043.114 50.527×10-32213.835.923.196 60.126×10-3
比较仿真结果可知，优化方案的储能能力与期望值相比大小相差在2%以内，且2
种方案的中心磁感应强度相差不大，符合设计要求.方案2的2 m外最小杂散磁场磁感应强度符合国家标准(小于5×10-4 T)，并小于方案1的杂散磁场，能够达到优化目标.
3 结语
以降低杂散磁场、节约超导带材为目的，对单螺管型超导磁体进行优化设计.以截
面为矩形的螺线管磁体为参照对象，在此基础上设计出截面为阶梯形的螺线管磁体.仿真实验结果表明，阶梯形截面的螺线管符合储能要求，且磁体外2 m处的杂散
磁场更小.同时，根据对体积的计算，阶梯型截面的磁体绕制过程中所需的超导带
材更少.随着超导磁体储能量的不断增大，接下来将进一步验证本研究的优化思路
在储能量更高的磁体上能否取得更好的优化效果，并积极探索计算量更小、带材使用量更少和磁场分布更优越的优化方法.
参考文献：
【相关文献】
[1] 何超峰，张俊峰，丁先庚，等.高温超导磁体试验装置设计[J].低温与超导，2012，40(9):55-57；
75.
[2] 朱光，刘建华，程军胜，等.25 T超导磁体优化中线圈数量影响分析[J].物理学报，2016，
65(5)[058401]:1-7.
[3] 朱晓辉，王海云，刘胜利.单螺线管型高温超导磁体的设计与优化[J].低温与超导，2016，
44(3):1-5；37.
[4] 朱颢东，钟勇.一种改进的模拟退火算法[J].计算机技术与发展，2009，19(6):32-35.
[5] 杨徉，罗映红，蒋宝灿，等.2.5 MJ单螺管型超导储能磁体的主动屏蔽研究[J].郑州大学学报(理学版)，2015，47(1):112-114.
[6] 张昊蒙.浅谈电力储能技术及其在现代电力系统中的应用[J].电子世界，2014(16):276.
[7] 李开彦.单螺管型高温超导储能磁体的静磁场分析[J].中国科技纵横，2015(21):149-152.。