Multi-Objective Optimization for Flywheel Motors Based on Parameter Priority Division
-
摘要:
为改善电机转矩性能,同时降低永磁材料成本,首先,根据飞轮电池的运行模式,提出飞轮电机的设计要求,并选取输出转矩、转矩脉动和永磁体用量为设计目标;在此基础上,采用等效磁路法和有限元法计算电机的主要尺寸和性能,同时确定合适的槽极配合以及初始结构参数;其次,使用相关系数法分析该电机各结构参数与优化目标之间的相关性,并合理划分转子尺寸和定子尺寸的优先级;最后,采用有限元代理模型和多目标优化算法(非支配排序遗传算法Ⅱ,NSGA-Ⅱ)相结合的方法对设计参数进行逐级寻优,并通过场路联合仿真和样机实验进行验证. 研究结果表明:采用参数优先级划分的方式可以减少优化过程中数据点的采样,缩短了整个过程中有限元分析的时间;优化后,电机反电势梯形化程度更高,转矩脉动减小了40%,永磁体用量降低了8%.
Abstract:This paper aims to improve the torque performance of a flywheel motor and reduce its permanent magnet cost. Firstly, according to the operating mode of the flywheel battery, the design requirements of the flywheel motor were proposed, and the output torque, torque ripple, and consumption of permanent magnets were selected as the design objectives. On this basis, the equivalent magnetic circuit method and finite element method were used to calculate the main size and performance of the motor, and the appropriate number of slots and poles and initial structural parameters were determined. Secondly, the correlation coefficient method was utilized to analyze the correlation between the structural parameters of the motor and the optimization objectives, and the priority of rotor sizes and stator sizes was divided reasonably. Finally, the finite element surrogate model and the multi-objective optimization algorithm, namely the non-dominated sorting genetic algorithm Ⅱ (NSGA-Ⅱ) were used to optimize the design parameters step by step, and the correctness of optimization results was verified by the field-circuit co-simulation method and prototype experiment. The results show that the proposed method of parameter priority division can reduce the sampling of data points in the optimization process and shorten the time of finite element analysis in the whole process. After optimization, the back electromotive force of the motor is more trapezoidal. The torque ripple is reduced by 40%, and the consumption of permanent magnets is reduced by 8%.
-
地铁具有交通量大、速度快、污染程度小、占地面积小等优点,能在很大程度上缓解城市地面交通压力,刺激城市经济,带动城市经济发展. 但由于施工成本高、工期长、资金回收率高,地铁大都会处于亏损状态,不得不寻求发展新模式. “地铁 + 上盖物业”可以有效利用和节约土地资源,带动相关产业的联合发展,从而提高地铁建设和运营的科学性和经济性.
地铁振动通过道床、隧道结构和基础传递到上层建筑,进而影响上层建筑和人类活动[1-2]. 国内外科学家对此做了大量研究. 邹超[3]以地铁车辆段及其上部建筑为研究对象,分析地铁车辆段及其上部建筑的振动传播规律,建立建筑振动预测模型,提出了有效的振动控制措施. 袁葵[4]分析试车线上盖建筑的车致振动机理、振动传播规律,在此基础上,提出减振降噪措施. 吕文婷[5]利用有限元软件建立车辆段上盖物业模型,并对其振动情况进行预测. 谢伟平等[6]结合2个实际工程,研究地铁列车进出车辆段对上部建筑特性振动的影响,并分析车辆段的振动荷载特性. Qu等[7]研究了地铁列车和公路车辆同时通过时的地面各测点振动响应. Real等[8]建立三维有限元模型,并用现场数据进行验证,预测铁路列车的振动. 赵彦辉[9]基于现场振动测试,探讨空沟和填充沟隔振措施的隔振效果. Yang等[10]通过相似模型试验,研究不同隧道断面对隧道及周围土体动力响应的影响. Zhang等[11]在模型试验中施加了扫频载荷,研究积分横向通道(ITP)的影响机理,提出一种新型的组装横向通道(ATP),并研究其减振效果. 郭治岳等[12]建立一个精细的数值分析模型,研究地铁附近建筑物的振动发展特征. 兰凯等[13]采用两步分析方法,预测列车以不同速度通过时引起的上盖建筑振动,并判断其是否满足标准要求. 孟坤等[14]认为在列车通行过程中,公路桥梁结构构件均存在明显的振动,但未发生扩散性共振效应,振动幅值很小.
国内外科学家也对地铁荷载进行了大量研究:陈行[15]使用国际公认的时间拟合公式确定高速列车的振动载荷;马龙祥等[16]基于无限-周期结构理论,建立车轨动力耦合模型,分析振动源车轨系统的竖向振动,计算列车在频率范围内通过钢轨传递到隧道基底的竖向振动激振力. 基于轮轨关系理论,刘凯[17]建立车辆-轨道系统的动力学分析模型,模拟拉伸载荷;刘维宁等[18]构建高精度的车轨耦合解析模型,并将计算结果与实测数据进行了对比.
综上所述,以往的研究主要集中在地铁振动对车辆段[19]或相邻建筑物的影响[20],而对不同结合类型地铁车站及其上方的动力反应特性的研究很少. 本文基于轨道车辆垂向耦合动力学理论,确定地铁B列列车的振动荷载,建立三维精细数值模拟模型,从时域、频域以及动应力角度,研究不同结合型式下地铁上盖物业的动力响应规律. 研究结果旨在为地铁车站上盖物业的结构设计和动力分析提供参考.
1. 地铁车站上盖物业结合型式
地铁车站与上盖物业的结合型式主要有层间夹土和板、柱等结构直接相连2种型式,在此将其分别归为“软结合”“硬结合”. 其中,“软结合”是上盖物业与地铁车站结构不直接相连,二者之间存在回填或原状夹层,两结构间无刚性连接. “硬结合”是上盖物业与地铁车站间采用刚性连接,具体又分为:1) 硬结合Ⅰ,物业首层作为地铁车站顶板(板共用);2) 硬结合Ⅱ,柱网结合型式(板不共用). 2种结合型式如图1所示.
2. 动力响应分析方法
2.1 列车荷载的确定
列车-轨道振动系统动力学模型分析法是利用解析法或者有限元数值法求解列车荷载,以轮轨不平顺视为车辆振动的激励源,并通过轮轨间赫兹接触关系将其传递到轨道结构中,从而形成整个车轨系统的动力学耦合过程. 车辆-轨道垂向耦合振动的动力模型[21]包括车辆模型、轨道模型和轮轨模型[22]. 车辆模型包括车体、2个转向架以及2组轮对,共计10个自由度.
2.2 数值仿真模型及参数
依托工程为佛山金融高新地铁站,该站为地下2层双柱双跨岛式站台,上盖物业上方有商场及裙楼,车站于下方斜穿上盖建筑,轨道为无砟形式. 采用FLAC3D建立数值仿真模型,鉴于动力计算的难度和耗时,仅考虑受列车振动影响大的上部物业结构. 为减小模型边界效应,确定模型长121.4 m、宽64 m、高72.5 m,地铁车站位于整体模型纵向中心位置. 为精确模拟波的传播,模型网格尺寸应小于输入波形最高频率对应波长的1/8~1/10[23]. 在轨行区、站厅层(A0)及上盖物业(A1~A4)楼板设置25个测点,模型、测点位置及编号见图2.
周边岩土体采用理想弹塑性模型并服从Mohr-Coulomb屈服准则[24],车站结构与上盖物业[25]均为线弹性模型[26]. 动力分析中选用静态边界条件,可有效减小模型边界上的入射波影响[27]. 具体做法是在模型的法向与切向分别设置阻尼器,从而实现吸收入射波的目的.
阻尼器所提供的法向黏性力tn和切向黏性力ts分别为
{tn=−ρCpvn,ts=−ρCsvs, (1) 式中:vn、vs分别为模型边界上法向和切向速度分量,m/s;ρ为介质密度,kg/m3;Cp、Cs分别为P波和S波的波速,m/s.
根据地铁站周边地层勘察资料,车站采用C30混凝土,上盖物业采用C25混凝土,地层、轨道和相关结构的物理力学参数取值如表1和表2所示.
表 1 地层及结构物理力学参数Table 1. Physical and mechanical parameters of strata and structures序号 名称 重度/(KN·m−3) 弹性模量/MPa 泊松比 厚度/m 内摩擦角/(°) 黏聚力/kPa 1 素填土 16.5 18.00 0.33 3.5 5.7 9.4 2 淤泥质粉细砂 17.1 30.00 0.25 6.0 23.0 3 粉质黏土 19.8 105.00 0.32 2.5 22.6 22.2 4 全风化粉砂岩 19.3 225.00 0.29 38 49.3 10.3 5 钢轨 78.5 2.01×105 0.30 6 地铁车站 24.0 3.00×104 0.20 7 上盖物业 24.0 2.80×104 0.20 表 2 轨道部件物理力学参数Table 2. Physical and mechanical parameters of track components钢轨 扣件 道床板 质量/
(kg·m−1)密度/
(kg·m−3)弹性
模量/GPa泊松比 垂向刚度/
(MN·m−1)扣件间距/m 弹性模量/
GPa泊松比 密度/
(kg·m−3)60 7850 205.9 0.30 59.2 0.6 32.5 0.24 2400 2.3 阻尼设置
瑞利阻尼C能够较好减弱系统自振模式的振幅,其矩阵计算如式(2)所示.
C=αM+βK, (2) 式中:α、β分别为与质量、刚度成比例的阻尼常数,M、K分别为质量、刚度矩阵.
阻尼常数可通过最小临界阻尼比和最小中心频率确定. 岩土材料临界阻尼比一般是2%~5%[24],而结构系统的临界阻尼比一般是2%~10%[28]. 模型基本构架大多为矩形块体,相对简单,可采用自振频率作为中心频率. 最终,通过试算得到中心频率为2.1 Hz,土体的质量和刚度阻尼常数分别为0.105、0.024,结构系统的质量和刚度阻尼常数分别为0.042、0.010.
2.4 列车荷载移动化
移动荷载与静载对结构的动力响应存在很大区别[29],在列车运行时,经过轨道上不同位置扣件的反力时程曲线在波形上除相位存在差异外,荷载时程曲线及频谱特性基本一致,该相位差的存在是由于列车通过相邻两观测点的时刻不同所导致的. 因此,将前述计算得出的列车荷载作用在钢轨上,施加移动的轮轨作用力. 首先,根据车速和沿列车行进方向单元的长度,计算出列车通过一个单元所需要的时间Δt,再将Δt除以动力时间步长,得到相邻单元所间隔的时间步数,进而实现对不同位置单元施加移动列车荷载. 荷载施加如图3所示,图中,ft为列车振动荷载.
3. 计算方法验证
3.1 现场测试
为验证列车荷载确定方法和构建数值模型的准确性,对依托工程的地铁列车振动进行了现场测试,将数值仿真结果与该地铁站测试数据进行时域和频域的对比分析. 振动信号采集仪器主要有加速度传感器及TST5912动态信号测试分析系统. 测试站点位于该站负1层的站厅层(A0)及上盖物业1~4层(A1~A4),安排专人负责采集列车经过时的振动信号.
3.2 计算模型与验证
该地铁站与上方裙楼及物业的连接方式是以物业首层楼板作为地铁站顶板,故采用2.2节“网格划分”给出的数值模型,仅去除掉软结合夹层,边界条件与材料阻尼参数均与2.2节和2.3节相同,数值仿真模型如图4所示.
结合列车速度信息,选取现场实测上盖物业1层(A1)加速度的时程数据,与列车速度40 km/h时数值仿真结果进行对比,并将两者加速度数据进行1/3倍频程转换,对现场实测与数值仿真的加速度级进行对比,结果如图5所示.
通过两者的对比分析可以看出:现场实测A1最大振动加速度为3.38 mm/s2,数值仿真计算站厅层最大振动加速度达到3.09 mm/s2,二者的误差为0.29 mm/s2;现场实测和数值仿真的1/3倍频程谱曲线贴合度较高,仅部分频率存在一定差异. 鉴于现场实测不可避免地会受到环境、设备及行人等影响,少量偏差在所难免,但二者数据分布基本一致,可以认为数值仿真方法与计算参数具有良好的可靠性.
4. 仿真结果分析
在数值仿真方法可行、计算参数合理的基础上,建立软结合、硬结合Ⅰ型、硬结合Ⅱ型3种型式下地铁车站上盖物业数值仿真模型,从时域、频域以及动应力出发,研究地铁列车以最大设计速度80 km/h运行时,上盖物业车致振动影响及差异.
4.1 振动加速度时程分析
限于篇幅,以断面3的测点振动数据为例进行分析,其振动加速度时程曲线如图6~8所示.
由图6~8可看出:1) 随着列车接近、穿越并驶离建筑物,振动加速度呈现明显的先增大、后减小的规律. 2) 软结合型式下,站厅层A0振动加速度峰值最大,为41.20 mm/s2,上盖物业A1~A4振动加速度峰值分别为12.73、15.61、16.78、13.78 mm/s2,A0到A1振动加速度峰值衰减了69.10%;硬结合Ⅰ型下,A0振动加速度峰值为43.84 mm/s2,A1~A4分别为42.93、50.36、33.56、44.26 mm/s2,A0到A1振动加速度峰值减小了2.08%;硬结合Ⅱ型下,A0振动加速度峰值为32.68 mm/s2,A1~A4分别为33.64、43.47、27.13、30.92 mm/s2,A0到A1振动加速度峰值增大了2.94%. 3) 软结合型式下,上盖物业各楼层振动加速度整体小于硬结合型式,且软结合型式A0到A1加速度衰减较大,而硬结合型式下变化很小,说明软结合型式的间隔土对列车振动有削减作用. 4) 硬结合Ⅰ型各楼层振动加速度略大于硬结合Ⅱ型,原因是硬结合Ⅱ型的车站-上盖物业刚度更大.
4.2 频域下动力响应分析
依据断面3各测点振动加速度的计算结果,采用快速傅里叶变换将时域结果转换到频域,通过MATLAB得到1/3倍频程各中心频率的振动加速度级. 3种结合型式下各测点加速度级峰值曲线如图9所示.
3种结合型式下各测点加速度级峰值曲线如图9所示. 以断面3为例,软结合型式下地铁车站上盖物业加速度级最大值位于A0,为79.5 dB,A1的值为68.2 dB,较A0减小了11.3 dB,说明地铁车站与上盖物业结构间的间隔土对列车振动有较大的衰减作用;除A1外,各楼层加速度级最大值均位于断面3,加速度级从断面3往两侧逐渐减小,说明对于同一楼层,随着距振源距离的增大振动减小. 硬结合Ⅰ型与硬结合Ⅱ型地铁车站上盖物业加速度级水平与竖向分布规律基本相同;对于上盖物业同一楼层,均是断面3加速度级最大,往两侧逐渐减小,说明随着距振源距离的增大振动减小;地铁车站上盖物业加速度级最大值均位于A4,硬结合Ⅰ型最大值为83.4 dB,硬结合Ⅱ型最大值为79.4 dB,均较A0有所增大. 硬结合型式下,上盖物业加速度级整体上较软结合大,并且硬结合Ⅰ型大于硬结合Ⅱ型. 因此,在条件允许的情况下,地铁车站上盖物业优先选择软结合型式,其次选择硬结合Ⅱ型.
图10为3种结合型A1各断面振动加速度频谱,图11为3种结合型断面3各楼层的振动加速度频谱. 由图10、11可以看出:3种结合型,上盖物业车致振动响应频率的分布规律基本一致,上盖物业车致振动频率主要集中在40~90 Hz;对于同一断面,各楼层振动响应频率范围及振幅变化不大.
4.3 动应力分析
3种结合型式下,断面3各楼层测点第一主应力时程曲线如图12所示,各主应力峰值如表3所示.
从图12、表3可以看出:软结合、硬结合Ⅰ型、硬结合Ⅱ型的第一主应力最大值均位于A0,分别为899.33、849.54、952.66 Pa,A1~A4第一主应力最大值远小于A0最大值;软结合型式第一主应力从A0到A1衰减了85.81%,衰减较大,主要是车站与上盖物业间隔土削弱了列车振动,进而减小了列车荷载作用下上盖物业的附加应力;硬结合Ⅰ、Ⅱ型第一主应力从A0到A1分别衰减了63.46%、72.27%,比软结合型衰减慢,主要是因为车站与上盖物业刚性连接,列车振动可以通过连接柱向上传播.
表 3 不同结合型断面3各楼层的第一主应力峰值Table 3. First principal stress peak value of the third section of each floor under different combination types楼层 软结合 硬结合Ⅰ 硬结合Ⅱ A0 899.33 849.54 952.66 A1 127.60 310.39 264.16 A2 113.31 264.24 180.05 A3 109.44 200.76 172.24 A4 120.30 233.98 192.38 衰减率/% 85.81 63.46 72.27 注:衰减率指A0到A1的第一主应力衰减率. 5. 结 论
采用现场测试和数值仿真手段对提出的地铁车站-上盖物业间的软结合、硬结合Ⅰ型和硬结合Ⅱ型3种型式的车致振动响应差异进行了研究,得到如下结论:
1) 依托现场实测数据验证了数值仿真方法与计算参数的可靠性,两者的A1最大振动加速度差值仅为0.29 mm/s2,且1/3倍频程谱曲线分布也基本一致.
2) 从A0到A1,振动加速度峰值软结合型减小了69.10%,硬结合Ⅰ型减小了2.08%,硬结合Ⅱ型增大了2.94%;硬结合型式下上盖物业振动加速度较软结合型式大.
3) 3种结合型式下,上盖物业振动的频率主要在40~90 Hz,且对于上盖物业同一楼层,车致振动随着距振源距离的增大而逐渐减小;软结合型A1加速度级最大值为68.2 dB,较A0减小11.3 dB;硬结合型上盖物业加速度级整体较软结合大,硬结合Ⅰ型上盖物业加速度级最大值为83.4 dB,硬结合Ⅱ型为79.4 dB.
4) 地铁列车振动造成上盖物业附加第一主应力很小,且在向上传播过程中衰减很快;从A0到A1,软结合型、硬结合Ⅰ型、硬结合Ⅱ型分别衰减85.81%、63.46%、72.27%,间隔土对附加应力有明显衰减作用.
5) 在相同列车荷载作用下,列车振动对软结合型地铁车站上盖物业影响最小,其次为硬结合Ⅱ型,最后为硬结合Ⅰ型;实际工程中,如果现场条件允许,建议优选软结合型式.
-
表 1 SPMSM主要尺寸参数
Table 1. Main size parameters of SPMSM
参数 电机 12 极 14 槽 24 极 28 槽 转速/(r·min−1) 4000 驱动方式 无刷直流 有效轴长/mm 8 转子外径/mm 100 定子外径/mm 73 定子内径/mm 40 气隙长度/mm 0.5 永磁体用量/kg 0.02 导体线径/mm 0.6 (8 股并绕) 电流密度/(A·mm−2) 8 槽满率/% 75 槽总面积/mm2 912 1248 每相串联匝数/匝 36 48 反电势平均值/V 6.52 8.70 输出转矩/(N·m) 0.78 1.03 表 2 电机设计变量初始值和优化范围
Table 2. Initial value and optimal ranges of motor design variables
mm 设计变量 初值 范围 Hs0 1.0 [0.8,1.2] Hs1 0.5 [0.4,0.6] Hs2 10 [8,12] Bs0 2.0 [1.6,2.4] Bs1 6.0 [4.8,7.2] aPM 0.90 [0.64,0.96] hPM 1.6 [1.2,1.8] hg 0.5 [0.2,0.8] 表 3 参数优先级
Table 3. Priority of parameters
优先级 目标 参数变量 优先级 1 Tavg, Tr Hs0,Hs1,Hs2,Bs1 优先级 2 Tavg,Tr,VPM aPM,hPM,hg 表 4 NSGA-Ⅱ算法的参数设定值
Table 4. Parameter setting values of NSGA-II algorithm
参数 设定值 初始种群/个 100 迭代数目/个 80 交叉概率/% 90 变异概率/% 5 最大迭代次数/次 20 最大允许帕累托百分比/% 70 表 5 优先级1优化结果
Table 5. Optimization results in in priority 1
设计变量 初始点 优化结果 Hs0/mm 1.0 0.8 Hs1/mm 0.5 0.4 Hs2/mm 10.0 11.9 Bs1/mm 6.00 5.58 Tavg /(N·m) 1.04 1.08 Tr/% 20.9 16.1 Kf/% 76 75 表 6 优先级2优化结果
Table 6. Optimization results in priority 2
设计变量 优先级 1 结果 优先级 2 优化结果 aPM /mm 0.9 0.8 hPM/mm 1.60 1.25 hg/mm 0.5 0.6 Tavg/(N·m) 1.08 1.02 Tr/% 16.1 14.7 VPM /mm3 2430 2234 -
[1] 张维煜,张林东,于焰均. 磁悬浮支承-飞轮系统稳定运行关键技术综述[J]. 西南交通大学学报,2022,57(3): 627-639.ZHANG Weiyu, ZHANG Lindong, YU Yanjun. Review on key technologies of stable operation for magnetic suspension support-flywheel system[J]. Journal of Southwest Jiaotong University, 2022, 57(3): 627-639. [2] 邓自刚,刘宗鑫,李海涛,等. 磁悬浮列车发展现状与展望[J]. 西南交通大学学报,2022,57(3): 455-474,530.DENG Zigang, LIU Zongxin, LI Haitao, et al. Development status and prospect of maglev train[J]. Journal of Southwest Jiaotong University, 2022, 57(3): 455-474,530. [3] LIU K, YIN M, HUA W, et al. Design and optimization of an external rotor ironless BLDCM used in a flywheel energy storage system[J]. IEEE Transactions on Magnetics, 2018, 54(11): 1-5. [4] LIU Z Q, WANG K, LI F. Design and analysis of permanent magnet homopolar machine for flywheel energy storage system[J]. IEEE Transactions on Magnetics, 2019, 55(7): 8105706.1-8105706.6. [5] 汤平华,漆亚梅,黄国辉,等. 定子无铁心飞轮电机绕组涡流损耗分析[J]. 电工技术学报,2010,25(3): 27-32.TANG Pinghua, QI Yamei, HUANG Guohui, et al. Eddy Current loss analysis of ironless flywheel electric machine’s winding[J]. Transactions of China Electrotechnical Society, 2010, 25(3): 27-32. [6] 付兴贺,江政龙,吕鸿飞,等. 电励磁同步电机无刷励磁与转矩密度提升技术发展综述[J]. 电工技术学报,2022,37(7): 1689-1702.FU Xinghe, JIANG Zhenglong, LÜ Hongfei, et al. Review of the blushless excitation and torque density improvement in wound field synchronous motors[J]. Transactions of China Electrotechnical Society, 2022, 37(7): 1689-1702. [7] LIU G H, ZENG Y, ZHAO W X, et al. Permanent magnet shape using analytical feedback function for torque improvement[J]. IEEE Transactions on Industrial Electronics, 2018, 65(6): 4619-4630. doi: 10.1109/TIE.2017.2772142 [8] 胡鹏飞,王东,靳栓宝,等. 弓形磁极永磁电机谐波削极技术的优化研究[J]. 中国电机工程学报,2020,40(18): 5987-5997.HU Pengfei, WANG Dong, JIN Shuanbao, et al. Research on the optimization of harmonic pole shaping technologies for arcuate pole permanent magnet motors[J]. Proceedings of the CSEE, 2020, 40(18): 5987-5997. [9] 赵品志,杨贵杰,李勇. 三次谐波注入式五相永磁同步电机转矩密度优化[J]. 中国电机工程学报,2010,30(33): 71-77.ZHAO Pinzhi, YANG Guijie, LI Yong. Torque density optimization for five-phase PMSM with third harmonic injection[J]. Proceedings of the CSEE, 2010, 30(33): 71-77. [10] 殷明,王森,陶侃凯,等. 双圈磁钢飞轮电机的设计及优化[J]. 飞控与探测,2020,3(2): 75-80.YIN Ming, WANG Sen, TAO Kankai, et al. Design and optimization of a double-circle permanent magnet motor used in reaction flywheel system[J]. Flight Control & Detection, 2020, 3(2): 75-80. [11] KIM K H, PARK H I, JANG S M, et al. Comparison of characteristics of double-sided permanent-magnet synchronous motor/generator according to magnetization patterns for flywheel energy storage system using an analytical method[J]. IEEE Transactions on Magnetics, 2015, 51(3): 1-4. [12] 罗成,张昆仑,王滢. 永磁电磁混合Halbach阵列电动悬浮的稳定性控制[J]. 西南交通大学学报,2022,57(3): 574-581.LUO Cheng, ZHANG Kunlun, WANG Ying. Stability control of electrodynamic suspension with permanent magnet and electromagnet hybrid Halbach array[J]. Journal of Southwest Jiaotong University, 2022, 57(3): 574-581. [13] ZHU X Y, WU W Q, QUAN L, et al. Design and multi-objective stratified optimization of a less-rare-earth hybrid permanent magnets motor with high torque density and low cost[J]. IEEE Transactions on Energy Conversion, 2019, 34(3): 1178-1189. doi: 10.1109/TEC.2018.2886316 [14] XU G H, ZHAO W X, LIU G H, et al. Torque performance improvement of consequent-pole PM motors with hybrid rotor configuration[J]. IEEE Transactions on Transportation Electrification, 2021, 7(3): 1561-1572. doi: 10.1109/TTE.2020.3041194 [15] 张邦富,程明,王飒飒,等. 基于改进型代理模型优化算法的磁通切换永磁直线电机优化设计[J]. 电工技术学报,2020,35(5): 1013-1021.ZHANG Bangfu, CHENG Ming, WANG Sasa, et al. Optimal design of flux-switching permanent magnet linear machine based on improved surrogate-based optimization algorithm[J]. Transactions of China Electrotechnical Society, 2020, 35(5): 1013-1021. [16] XU L, ZHAO W X, LIU G H, et al. Design optimization of a spoke-type permanent-magnet vernier machine for torque density and power factor improvement[J]. IEEE Transactions on Vehicular Technology, 2019, 68(4): 3446-3456. doi: 10.1109/TVT.2019.2902729 [17] 刘国海,王艳阳,陈前. 非对称V型内置式永磁同步电机的多目标优化设计[J]. 电工技术学报,2018,33(增2): 385-393.LIU Guohai, WANG Yanyang, CHEN Qian. Multi-objective optimization of an asymmetric V-shaped interior permanent magnet synchronous motor[J]. Transactions of China Electrotechnical Society, 2018, 33(S2): 385-393. [18] 赵玫,于帅,邹海林,等. 聚磁式横向磁通永磁直线电机的多目标优化[J]. 电工技术学报,2021,36(17): 3730-3740.ZHAO Mei, YU Shuai, ZOU Hailin, et al. Multi-objective optimization of transverse flux permanent magnet linear machine with the concentrated flux mover[J]. Transactions of China Electrotechnical Society, 2021, 36(17): 3730-3740. [19] 张河山,邓兆祥,妥吉英,等. 基于改进人工蜂群算法的轮毂电机多目标优化[J]. 西南交通大学学报,2019,54(4): 671-678.ZHANG Heshan, DENG Zhaoxiang, TUO Jiying, et al. Multi-objective optimum design for in-wheel motor based on improved artificial bee colony algorithm[J]. Journal of Southwest Jiaotong University, 2019, 54(4): 671-678. [20] 陈云云,朱孝勇,全力,等. 基于参数敏感度的双凸极永磁型双定子电机的优化设计和性能分析[J]. 电工技术学报,2017,32(8): 160-168. doi: 10.19595/j.cnki.1000-6753.tces.2017.08.019CHEN Yunyun, ZHU Xiaoyong, QUAN Li, et al. Parameter sensitivity optimization design and performance analysis of double-salient permanent-magnet double-stator machine[J]. Transactions of China Electrotechnical Society, 2017, 32(8): 160-168. doi: 10.19595/j.cnki.1000-6753.tces.2017.08.019 [21] LEI G, LIU C C, ZHU J G, et al. Techniques for multilevel design optimization of permanent magnet motors[J]. IEEE Transactions on Energy Conversion, 2015, 30(4): 1574-1584. doi: 10.1109/TEC.2015.2444434 [22] ZHAO W X, YAO T, XU L, et al. Multi-objective optimization design of a modular linear permanent-magnet vernier machine by combined approximation models and differential evolution[J]. IEEE Transactions on Industrial Electronics, 2021, 68(6): 4634-4645. doi: 10.1109/TIE.2020.2988233 [23] 李培强,李慧,李欣然. 基于灵敏度与相关性的综合负荷模型参数优化辨识策略[J]. 电工技术学报,2016,31(16): 181-188.LI Peiqiang, LI Hui, LI Xinran. Optimized identification strategy for composite load model parameters based on sensitivity and correlation analysis[J]. Transactions of China Electrotechnical Society, 2016, 31(16): 181-188. [24] SUN X D, SHI Z, LEI G, et al. Multi-objective design optimization of an IPMSM based on multilevel strategy[J]. IEEE Transactions on Industrial Electronics, 2021, 68(1): 139-148. doi: 10.1109/TIE.2020.2965463 [25] CROUX C, DEHON C. Influence functions of the Spearman and Kendall correlation measures[J]. Statistical Methods & Applications, 2010, 19(4): 497-515. [26] 黄允凯,周涛. 基于等效磁路法的轴向永磁电机效率优化设计[J]. 电工技术学报,2015,30(2): 73-79.HUANG Yunkai, ZHOU Tao. Efficiency optimization design of axial flux permanent magnet machines using magnetic equivalent circuit[J]. Transactions of China Electrotechnical Society, 2015, 30(2): 73-79. [27] 周扬,周瑾,张越,等. 基于RBF近似模型的磁悬浮轴承结构优化设计[J]. 西南交通大学学报,2022,57(3): 682-692.ZHOU Yang, ZHOU Jin, ZHANG Yue, et al. Optimum structural design of active magnetic bearing based on RBF approximation model[J]. Journal of Southwest Jiaotong University, 2022, 57(3): 682-692. [28] LIU G H, WANG Y Y, CHEN Q, et al. Multiobjective deterministic and robust optimization design of a new spoke-type permanent magnet machine for the improvement of torque performance[J]. IEEE Transactions on Industrial Electronics, 2020, 67(12): 10202-10212. doi: 10.1109/TIE.2019.2962472 [29] LI Y W, ZHU C S, WU L J, et al. Multi-objective optimal design of high-speed surface-mounted permanent magnet synchronous motor for magnetically levitated flywheel energy storage system[J]. IEEE Transactions on Magnetics, 2019, 55(7): 8202708.1-8202708.8. -