基于非均匀气隙的永磁电机气隙磁密优化

《微电机》投稿收稿日期：××××-××-××基金项目：北京市优秀人才基金项目作者简介：刘明基(1969)，博士/副教授，研究方向：新型电机理论及运行控制；韩志(1988),硕士研究生，研究方向：新型电机及现代节能控制技术；裴亚伟(1988)，硕士研究生，研究方向：新型电机及现代节能控制技术。基于非均匀气隙的永磁电机气隙磁密优化刘明基，韩志，裴亚伟（华北电力大学，新能源电力系统国家重点实验室，北京，102206）摘要：传统方法设计的永磁电机同步电机气隙磁密谐波含量较大，对电机性能造成不利影响。本文采用不均匀气隙方法对内置“V”型转子磁路结构的永磁电机进行优化。首先，利用Taguchi方法及有限元法，以永磁电机的4个结构参数为优化变量，2种电机特性为约束条件，进行多目标多变量优化，得到了电机优化的空载气隙磁密波形；然后，将该方法与遗传算法进行对比分析，验证Taguchi方法优化结果的优越性；最后，运用时步有限元方法计算得到了考虑定子斜槽的感应电动势波形。关键词：永磁电机；磁极优化；Taguchi方法；有限元法；不均匀气隙中图分类号：×××××；×××××文献标志码:AMagneticFieldWaveformOptimizationofPermanentMagnetSynchronousMotorwithNo-uniformAirGapLIUMingji,HANZhi,PEIYawei（StateKeyLaboratoryofAlternateElectricalPowerSystemwithRenewableEnergySources,NorthChinaElectricPowerUniversity,Beijing,102206,China）Abstract:Harmonicofairgapfluxdensitywaveformofthepermanentmagnetsynchronousmotor(PMSM)designedbytraditionalmethodisalittlelarge,itisharmfultothemotorperformance.AkindofV-typerotorstructurepermanentmagnetmotorwasoptimizedbyNo-Uniformairgapmethodinthispaper.First,4parametersofthepermanentmagnetmotorwereoptimizedthroughTaguchimethodandfiniteelementmethod(FEM).Twomotorcharacteristicsweretakenasconstraintconditions.Bymultipleobjectiveandmultiplevariablesoptimization,aoptimizedairgapfluxdensitywaveformwasget;Then,Taguchimethodwerecomparedwithgeneticalgorithm,andTaguchimethodhasadvantages;Intheend,inductionelectromotiveforceconsideringskewedstatorslotwascalculatedbytime-steppingFEM.Keywords:PMSM；OptimizationofMagneticPole；TaguchiMethod；FEM；No-UniformAirGap中图分类号：×××××；×××××文献标志码：A文章编号：1001-6848（××××）××-××××-×××0引言在用传统方法设计的永磁电机中，空载气隙磁密的谐波含量较大，导致定子绕组感应电动势波形发生畸变，谐波电流含量较大，增加了谐波损耗，并导致转矩波动，引起机械振动及噪声，因此，对空载气隙磁密波形进行优化具有重要意义[1]。在永磁电机气隙磁密优化方面，国内外已有大量文献。文献[2-3]对表面式磁钢结构永磁电机包括偏心距在内的几个结构参数进行了优化；西南交通大学徐英磊采用有限元方法分析了偏心设计的V形内永磁电机的磁场，得到了偏心距与磁场的关系，用曲线拟合找到最小的THD点[4]；华北电力大学陈超采用不均匀气隙结构，利用遗传算法以及响应面法对内置V型转子磁路结构自起动永磁电机磁极形状进行优化，得到了转子磁极最优偏心距，从而较好的优化了空载气隙磁密，提高了电机的稳态运行性能。但该方法存在优化周期长，易陷入局部最优解的问题[5]；文献[6-7]利用Taguchi方法，将等效磁路法和有限元法相结合，对永磁电机内置式V型永磁体的磁钢厚度、磁极夹角、磁钢距转子边界及转轴外边界距离进行多目标优化设计，通过改变各目标函数的权重得到不同的设计结果。该方法被认为在初步优化磁极阶段较简便，节约实验成本，优化结果具有较高的鲁棒性，从而为产品的大规模生产提供保障。《微电机》投稿本文将Taguchi方法与不均匀气隙法相结合，对内置“V”型转子磁路结构的永磁电机磁密波形进行优化，并通过与遗传算法优化的永磁电机性能进行对比分析验证Taguchi方法优化结果的优越性。1永磁同步发电机的转子结构及电机参数1.1永磁同步发电机转子原始结构本文采用内置“V”型磁路结构的永磁电机作为优化对象，具体转子结构示意图如图1所示，该电机的结构参数如表1所示。图1.永磁同步发电机转子结构表1永磁同步发电机的主要参数参数参数值参数参数值额定功率[kW]11极数4定子外径[mm]260定子内径[mm]170定子槽数36定子联接方式星形永磁矫顽力[kA·1m]772永磁磁化长度[mm]5.12优化方法2.1Taguchi方法Taguchi方法的参数设计是在应用中用得最多、效果比较好的一种设计方法，日本著名的统计学家田口玄一将正交试验选择的因子水平组合成表格，称为正交表。Taguchi方法中的正交试验设计根据正交性从全面试验中挑选出部分有代表性的点进行试验，能够高效、快速、经济地得到满意的实验结果，是研究多因子多水平的一种有效设计方法[8-9]。在正交实验设计方法中，影响品质特性的条件称为因子，因子的取值叫因子的水平。2.1.1Taguchi方法的一般步骤1.选定品质特性2.判定品质特性的理想机能3.列出所有影响品质特性的因子4.确定因子的水平5.根据因子及水平确定正交表6.把变量值映射到正交表中7.进行试验，记录数据8.分析各因子水平对品质特性的影响，找出因子的最优组合9.确认实验。2.2磁场优化的不均匀气隙法传统方法设计的永磁电机多为均匀气隙，转子冲片为标准的圆形结构。不均匀气隙法采用偏心距设计，假设转子半径为R，在一个磁极范围内，圆心向外移动了距离D，以该点为新的圆心，以R-D为半径画弧，所得转子冲片结构将不再是圆形，而是梅花形，此时定转子间为不均匀气隙，如图2所示。3转子结构的优化设计3.1实验设计本文以磁极夹角、磁钢距转轴距离、气隙以及偏心距4个参数为优化变量，如图2所示。表2给出了各优化参数受转子结构限制的变化范围。以空载气隙磁密谐波含量（THD）及空载气隙磁密基波幅值（Bm1）为约束条件（品质特性），本文要求Bm1在0.8T左右。其中THD的计算公式为：22BmnTHD100%Bm1n（1）RABDR-DC图2.V型内永磁转子一个极的结构示意图表2各设计变量的变化范围因素数值最小值最大值A:永磁体夹角[]110140B:永磁体距转子距离[mm]814C:气隙[mm]0.81.1D:偏心距[mm]118针对以上4个待优化的变量，分别用A、B、C、D表示，本文采用16L(44)正交表[8-9]，如表3所示。将A、B、C、D四个变量分别在取值范围内分为四组，表4为依据表3列出的16组实验参数组合。表316L(44)正交表实验变量变量变量变量111432123231314414215211162224723428243393134103241113323123412134122144213154331164444表44变量4水平正交表实际值实验A[][]Bmm[]Cmm[]Dmm111081.112211010163110120.8184110140.91512080.816120100.9187120121.168120141129130811810130101.1111130120.91212130140.861314080.9614140100.81215140121116140141.118第1组实验所对应的空载气隙磁密波形及其谐波含量如图3所示。由图3（a）可知，由于大量谐波的影响空载气隙磁密波形近似为尖顶波；由图3（b）可知，3、5、7、17、19次谐波含量较大，其中17、19次谐波为齿谐波，通过改变转子结构及定子绕组联接方式均难以削弱，将在后文采用定子斜槽法削弱其对电动势的影响。当电机对称星形连接时，线电压中不含3次及其倍数次谐波，也不存在3次及其倍数次谐波电流[10]。而3次及其倍数次谐波还可以起到降低磁密波形幅值、提高基波幅值的作用，通常认为3次谐波含量在15%左右是合理的。因此本文仅对去除3次及其倍数次谐波和齿谐波后的谐波进行约束，此时THD的计算公式改写为：622,4Bm(21)THDTHD_3x_17_19Bm1iii（2）6012018000.050.10.150.20.250.30.350.40.451357911131517190102030405060708090100谐波含量[%]磁密[T]空间电角度[º]谐波次数（a）气隙磁密波形（b）谐波分析图3.第1组实验所得空载气隙磁场对表4中的16组电机结构分别进行有限元计算，得到各参数组合下的THD及Bm1如表5所示。表5各变量组合试验结果实验1234THD[%]13.2319.1611.4825.69Bm1[T]0.58760.61820.55890.6174实验5678THD[%]14.899.1217.2912.72Bm1[T]0.73730.60310.64010.5991实验9101112THD[%]8.3810.656.5912.28Bm1[T]0.64170.72870.67930.7217实验13141516THD[%]8.026.1411.246.99Bm1[T]0.78560.75320.77430.6368由表5可知第14组实验结果较好。为进一步优化，下面对各因子水平对品质特性的影响进行研究。表6列出了各因子水平对Bm1的影响，以因子《微电机》投稿B在水平2下对THD的影响为例，计算公式为[6]：B2261014114(1111)BmBmBmBmBm（3）表6不同因子水平对应的空载气隙基波含量因子水平Bm1ofA[T]Bm1ofB[T]Bm1ofC[T]Bm1ofD[T]10.59550.68810.69280.714420.64490.67580.67140.691430.69260.66320.65830.654840.73750.64380.64830.6101将表6中数据进行图形表示，如图4(a)所示。可以清楚的看出不同因子水平对应得Bm1的大小。由图4(a)可知使Bm1最大的组合是A4B1C1D1。同理可得空载气隙磁密谐波含量随各因子水平的变化情况，如图4(b)所示，使空载气隙磁密基波幅值最小的组合是A4B3C1D4。这两组数据将作为调整结构的重要依据。0.580.60.620.640.660.680.70.720.740.76Bm1[T]A4A3A2A1B4B3B2B1C4C3C2C1D4D3D2D1A,B,C,D参数设置(a)(b)图4.各因子水平对Bm1及THD的影响由表5可以计算得到Bm1和THD的平均值，分别表示为Bm1及THD。为确定决定各特性的关键因子，对各参数对Bm1和THD的影响进行方差分析。以因子B对Bm1的影响为例，其计算公式为[6]：24BBii1S