卫红B用无别直流电机齿槽转矩仰制方法研究

刘金刚 郑剑云 张聪悦 傅兵 王高升 肖绍熙

摘要：由于齿槽结构的影响，无刷直流电机（brushless direct current motor，BLDCM）中存在齿槽转矩，会导致转矩脉动增加，从而影响无刷直流电机的性能.为此，以电子液压制动系统（Electro-hydraulic brake system，EHB）用BLDCM为研究对象，对结构参数进行分析及改进，来抑制BLDCM的齿槽转矩.首先，通过对齿槽转矩表达式进行推导和分析，发现其数值由基本结构参数以及傅里叶分解系数B（rnz/2p）和Gn所决定;其次，分别推导了偏心距、极弧系数和辅助槽与傅里叶分解系数的关系式，进而得到三者与齿槽转矩的变化关系;最后，通过有限元方法研究了偏心距改变、极弧系数变更和辅助槽添加三种情况下齿槽转矩的变化情况.结果表明：当偏心距、极弧系数和辅助槽形状分别为7.90mm、0.74、椭圆形时，对应的齿槽转矩峰值分别降低96.65%、90.48%、86.74%，齿槽转矩都得到明显抑制.

关键词：无刷直流电机;结构参数;转矩抑制;有限元分析

中图分类号：TM351

文献标志码：A

收稿日期：2021-07-06

基金项目：湖南省杰出青年基金资助项目（2019JJ20015）;国家自然科学基金资助项目（52075465）;湖南省战略性新兴产业科技攻关与重大科技成果转化项目（2019GK4025）;湖南省科技创新计划资助项目（2020RC4038）;湖南省自然科学基金青年项目（2020JJ5540）

作者简介：刘金刚（1979—），男，山东聊城人，湘潭大学教授，博士

†通信联系人，E-mail：2682833218@qq.com

无刷直流电机（BLDCM）具有体积小、质量轻、功率密度高、可靠性高等优点，非常适合用作电子液压制动系统（EHB）的助力电机[1-3].然而，由于转子永磁体与定子铁心的齿槽相互作用，包括BLDCM在内的齿槽结构永磁电机都会产生齿槽转矩[4].齿槽转矩的存在将导致转矩脉动增加，带来振动与噪声，影响无刷直流电机的性能.

如何降低永磁电机齿槽转矩，至今都是国内外学者研究热点之一，众多学者在电机本体结构设计改進方面做了大量工作[5]，在定子侧的结构参数改进上，WangD等人[6]在保证槽口宽度相同前提下，采用非均匀齿冠宽度的方法，并通过解析分析确定了齿冠宽度比，大大降低齿槽转矩;黄守道等人[7]根据齿槽转矩解析表达式，分析了单个槽对齿槽转矩增加规律的影响，提出通过槽口偏移来减弱齿槽转矩的方法;WangK等人[8]利用定子齿的偏心结构来减小齿槽转矩，并在一台外转子PMSM上验证了所提方法的有效性.在转子侧的结构参数改进上，安治国等人[9]研究了转子辅助槽添加位置、开槽深度和宽度等参数对齿槽转矩和漏磁的作用情况，获得最佳开槽结构参数;ZhaoJ等人[10]对永磁体边缘的结构参数进行了研究，通过改变参数减少了转子端部效应的影响，抑制了齿槽转矩;李松生等人[11]进行了转子永磁体对齐安装和错位安装的对比，发现分段斜极安装可以大幅降低齿槽转矩.上述提出的结构改进方法能降低电机齿槽转矩，但鲜有同时兼顾永磁体形状参数和电枢参数改进的方法，也很少综合考虑齿槽转矩和电机其他基本性能.

本文通过对齿槽转矩的表达式进行推导和分析，发现其数值由基本结构参数以及傅里叶分解系数B（rnz/2p）和Gn所决定，以此为依据，推导了傅里叶分解系数和结构参数的关系式.再利用有限元方法分别研究了BLDCM偏心距改变、极弧系数变更和添加辅助槽三种情况下齿槽转矩、气隙磁密和反电动势的变化情况.结果表明，对结构参数进行合理地改进，可以在保持电机性能基本不变的情况下抑制齿槽转矩.

1结构参数与齿槽转矩的关系分析

1.1齿槽转矩分析

由于齿槽结构的影响，当转子永磁体处于不同位置时，磁路中的磁阻会发生变化，从而引起磁场能量变化，导致齿槽转矩的产生.因此，将齿槽转矩定义为磁场能量W相对于位置角a的函数[12]，即：

为便于分析，做以下假设：定子铁心磁导率为无穷大;不同永磁体之间结构及磁性能处于理想状态，无加工制造误差，磁导率与空气相同;硅钢片的叠装系数为1.因此，永磁电机磁场中存储的磁能可近似[13]为气隙与永磁体磁能之和，可表示为：

式中，B为气隙磁通密度;μ0为磁导率;V为气隙体积.其中气隙磁密沿电枢表面分布可表示为：

式中，B（rθ）为永磁体剩磁沿圆周方向的分布;hm（θ）、δ（θ，α）分别为沿圆周、磁化方向的气隙长度.

将式（2）和式（3）代入式（1），进行傅里叶分解展开并求导，可得齿槽转矩表达式如下：

式中，R1、R2分别为转子外径和定子内径;Z为槽数;L为铁心长度;Gn和B（rnz/2p）为傅里叶分解所对应n次和nz/2p次谐波系数;p0为极对数.

由式（4）可知，齿槽转矩由基本结构参数、分解系数B（rnz/2p）和Gn决定.电机设计完成后，R1、R2、Z、L等基本结构参数已经确定，因此抑制齿槽转矩可从分解系数B（rnz/2p）和Gn上着手.B（rnz/2p）的数值主要由永磁体偏心距、极弧系数等决定，而Gn的幅值主要与槽口宽度和辅助槽添加等有关.

1.2偏心距与齿槽转矩的关系分析

BLDCM转子采用面包形表贴式结构，永磁体形状如图1所示，为方便分析将面包形永磁体分为a、b两部分，b部分的圆心与转子圆心重合.

永磁磁体a、b两部分充磁厚度之和为h′以表示为：

式中，Ra、Rb分别表示a、b部分的半径;Δh表示偏心距，是永磁体外部弧面圆心Oa与转子中心Ob的距离;lδ表示永磁体底部平面与转子中心Ob的距离.

面包形永磁体偏心时，气隙长度随之产生变化，但实际气隙长度与永磁体的充磁方向厚度之和保持不变，因此径向气隙磁密B′（θ，α）可以表示为：

进一步可以得到B’（0），其表达式为：

式中，hmax表示永磁体最大充磁厚度.将hm′（θ）的表达式（5）代入式（7），得到B2（0）的表达式：

永磁体偏心设计的目的是使得气隙不均匀，从而获取预想的气隙磁通密度分布，其程度可通过偏心距Δh的数值大小来体现.以Δh为变量时，不考虑极弧系数km的影响，将km设置为初始值0.86，则Br′2（θ）可转化为关于偏心距Δh和位置θ的函数，即：

经过推导，得到B2（0）在分解区间[-πk2p，πk/2p]上的傅里叶各次分解系数Bm：

由式（10）和式（4）可知，随着Δh的增大，Br′2（θ）的傅里叶分解系数Bm降低，进而减小齿槽转矩.

1.3极弧系数与齿槽转矩的关系分析

永磁体极弧的宽度bm和极距τ之比被称之为极弧系数km，其表达式为：

以极弧系数km为改进变量，研究其对齿槽转矩Tc的影响时，不考虑偏心情况，设置Δh取值为零，则Br′2（θ）可以表示为关于位置角θ和极弧系数km的函数，其表达式为：

B，2（0）在区间[-Tkn/2p，πkn/2p]上的傅里叶各次分解系数B的函数，可以转换为与极弧系数飞相关的函数，其表达式如下：

由式（13）和式（4）可知，调整极弧系数km可以使得分解系数Bm取较小值，从而可以对齿槽转矩进行抑制.

1.4辅助槽与齿槽转矩的关系分析

电枢是电机实现能量转化中起到关键枢纽作用的部分，在BLDCM中电枢指的是定子部分，其设计的关键在于槽形参数选取.槽形参数的设计会影响

在定子内表面的分布，改变其傅里叶展开系数Gn的取值，进而改变齿槽转矩Tc的大小.

在分析电枢参数对Tc影响情况时，可忽略转子

其中，槽口宽度β（以弧度制表示）是重要的槽形参数之一，其数值会直接改变

的分解区间，槽口宽度与积分区间的关系如图2所示.

考虑槽口宽度影响时，

式中：β为用弧度表示的槽口宽度，

b0表示定子槽口宽度.

由式（14）可知，适当地选取槽口宽度，可以改变

Gn，从而抑制齿槽转矩.然而，槽口过小时可能会引起漏磁，过大时容易造成定子齿部饱和，因此槽口取值范围要合理選取.

槽口宽度β的改变常与辅助槽添加配合使用，添加辅助槽等于增加BLDCM开槽数目，更改极槽配比，使得新添加的槽与原先的齿槽转矩进行抵消，进而对其进行削弱.矩形的定子辅助槽均匀添加在齿顶上，宽度为β，与原槽口的宽度一致，具有一定开槽深度，辅助槽添加示意图如图3所示.

傅里叶展开区间会发生改变，当每齿添加1个辅助槽的情况下，∫得∫到的Gn如下所示：

由式（15）和式（4）可知，适当地选取辅助槽，可以改变傅里叶分解系数Gn，从而抑制齿槽转矩.

2无刷直流电机有限元建模

理论分析了偏心距、极弧系数和辅助槽与齿槽转矩傅里叶分解系数之间的关系，发现改变三者的数值大小可以减小齿槽转矩.但是，这是基于假设推导得到的，而实际情况十分复杂，需要仿真来进一步分析.有限元方法可以全面考虑电机漏磁、磁饱和等复杂情况，相较于理论分析更为精确.因此，为准确分析齿槽转矩敏感结构参数对齿槽转矩和电机其他性能的影响，本文基于电机结构参数，如表1所示，搭建了无刷直流电机三维模型，如图4所示.并使用Ansoft Maxwell软件建立了8极12槽BLDCM的有限元分析模型，对其进行材料定义、边界设置、激励添加和网格剖分，建立的有限元模型如图5所示.

在有限元求解设置中，将电机设置为恒转速工况，绕组电流设置为零，以永磁体磁场为激励，对该电机进行有限元仿真计算.

3仿真结果分析

3.1偏心距改变对齿槽转矩的影响为获取能有效抑制齿槽转矩的偏心距数值，在其取值范围0.00～9.00mm内，以1.00mm为步长进行仿真分析.得到齿槽转矩峰值Tcmax随偏心距Δh的变化情况，如图6所示.

由图6可知，随偏心距Δh的增加齿槽转矩峰值Tcmax呈现先减小后增大的趋势，这与理论分析中的Bm变化趋势相一致.同时，发现当偏心距Δh取值为8.00mm时，齿槽转矩峰值取到最小值3.51mN·m;当Δh取值为9.00mm时次之，峰值为5.74mN·m.

为更准确地分析偏心距改变对齿槽转矩的影响，对上述取值进行再次细分后仿真.在有限元分析中，偏心距取值以8.00mm为中心，在7.00～9.00mm范围内，以0.10mm为步长进行仿真，Tcmax随Δh变化情况，如图7所示.由图7可知，Δh=7.9mm时齿槽转矩出现最小值，峰值为3.48mN·m.

偏心距Δh改变前后齿槽转矩Tc的曲线波形对比、气隙磁密对比和反电动势曲线的波形对比分别如图8～10所示.由图8可知，改变前永磁体不偏心，即Δh=0mm，齿槽转矩峰值Tcmax为103.91mN·m;改变后，Δh为7.90mm，Tcmax减少为3.48mN·m，相比改变前降低了96.65%.由图9可知，改变前，气隙磁密Bδ的峰值为0.88T;改变后，由于永磁体最大厚度没有改变，所以气隙磁密Bδ的峰值基本不变，但是由于Δh的增加，导致曲线波形变尖.由图10可知，改变前，反电动势E的峰值为5.27V;改变后，波形平顶宽度变窄，反电动势E的峰值减小为4.71V，相较于改变前降低了10.66%.

3.2极弧系数变更对齿槽转矩的影响在极弧系数取值范围0～1内，以0.1为步长进行仿真，得到齿槽转矩峰值Tcmax随极弧系数km变化情况，如图11所示.

由图11可知，Tcmax随km变化而波动，km取值为0.40时，Tcmax出现最小值为10.84mN·m;km取值为0.70时次之，峰值为46.40mN·m.为更准确地分析极弧系数变更对齿槽转矩的影响，对上述取值再次细分后仿真.在有限元仿真分析中，极弧系数取值以0.70为中心，在0.60～0.80范围内，以0.01为步长进行仿真，得到齿槽转矩峰值Tcmax随极弧系数km的变化情况，如图12所示.从图12可以看出，当时齿槽转矩出现最小值，峰值为9.89mN·m.李辉等人[14]观察到齿槽转矩峰值与极弧系数之间有相似的趋势.

极弧系数km变更前后齿槽转矩Tc的曲线波形对比、气隙磁密波形对比和反电动势曲线分别如图13～15所示.由图13可知，极弧系数变更前，齿槽转矩峰值T为103.91mN·m;变更后，Tcmax减小到9.89mN·m，相比于变更前降低了90.47%，齿槽转矩峰值得到了明显抑制.由图14可知，变更前，气隙磁密Bδ的峰值为0.88T;变更后，由于永磁体最大厚度不变，所以气隙磁密的峰值基本不变，但极弧系数的减小会导致曲线顶部的宽度略微变窄.由图15可知，变更前，反电动势E的峰值为5.27V;变更后，反电动势E的峰值降低了0.76%，变为5.23V.

3.3添加辅助槽对齿槽转矩的影响

在齿槽槽口取值范围0～5mm内进行仿真分析，得到了齿槽转矩峰值Tcmax随齿槽槽口宽度b0的变化情况，如图16所示.

由于b0的取值受到限制，在现有范围内并没有出现最佳的取值，可以通过添加辅助槽等手段进一步对Gn进行抑制，进而减少齿槽转矩.但是，添加辅助槽后定子铁心的结构发生变化，磁导率随之改变，从而改变齿槽转矩，而不同形状辅助槽形产生的影响也会不同，因此需要对辅助槽的形状进行分析.本文针对矩形、三角形和椭圆形3种不同形状的定子辅助槽进行分析，不同槽形比较时，将槽宽和槽深设置为不变量，以便比较形状对齿槽转矩的影响，添加不同形状辅助槽时齿槽转矩Tc的曲线波形对比、气隙磁密对比和反电动势曲线分别如图17～19所示.

由图17可知，添加的辅助槽形状为矩形时，齿槽转矩峰值Tcmax为23.10mN·m;为三角形辅助槽时，Tcmax增大为29.59mN·m，相较于矩形辅助槽增加了28.10%;为椭圆形辅助槽时，Tcmax降低到13.78mN·m，相较于矩形辅助槽降低了40.35%，因此椭圆形辅助槽对齿槽转矩的削弱效果最佳.由图18可知，添加辅助槽前，气隙磁密Bδ的曲线波形由于受到原槽口作用，顶部存在1处凹陷;而添加辅助槽后，凹陷数量增加至3个，相比于椭圆形辅助槽，添加矩形辅助槽和三角形辅助槽时气隙磁密波形会变得更尖.由图19可知，添加辅助槽前，反电动势E的峰值为5.27V;而添加不同形状辅助槽后，对应的反电动势E的峰值也有着不同程度地降低，其中添加椭圆形辅助槽时变化最小，其反电动势E的峰值为4.95V，相较于开槽之前仅降低了6.07%.

4结论

本文以EHB用无刷直流电机为研究对象，通过对齿槽转矩的表达式进行推导和分析，发现其数值由基本结构参数以及傅里叶分解系数B（rnz/2p）和Gn所决定，以此为依据，推导了结构参数与傅里叶分解系数的关系式.并分别研究了偏心距改变、极弧系数变更和添加辅助槽对齿槽转矩和电机其他性能的影响，得到了以下结论.

1）当偏心距、极弧系数和辅助槽形状分别为7.90mm、0.74、椭圆形时，对应的齿槽转矩峰值分别为3.48mN·m、9.89mN·m、13.78mN·m，相比于改变前分别降低96.65%、90.48%、86.74%，齿槽转矩都得到了明显抑制.

2）偏心距改变、极弧系数变更和添加辅助槽会增加气隙磁密波形的凹陷数量，也会使反电动势曲线变窄，说明其对电机性能产生了一定的影响，相较而言，添加椭圆形辅助槽的影响最小.

参考文献

[1]时培成，王锁，张荣芸，等.无刷直流电机两管导通零矢量直接转矩控制[J].中国机械工程，2020，31（6）：670-678.

[2]王伟达，武燕杰，史家磊，等.基于驾驶员意图识别的电子液压制动助力系统控制策略[J].吉林大学学报（工学版），2021，51（2）：406-413.

[3]谷贺冲.乘用车电机直驱线控制动系统设计与控制研究[D].长春：吉林大学，2019.

[4]张昌錦，陈志辉，梅庆枭，等.基于转子齿偏移的四相横向磁通永磁电机齿槽转矩削弱方法[J].中国电机工程学报，2018，38（12）：3680-3687.

[5]石松宁，王大志.消弱永磁驱动器齿槽转矩的斜极优化方法[J].电工技术学报，2015，30（22）：27-32.

[6]WANG D H，WANG X H，QIAO D W，et al.Reducing cogging torque in surface-mounted permanent-magnet motors by nonuni⁃formly distributed teeth method[J].IEEE Transactions on Mag⁃netics，2011 ，4 7（ 9 ）：2 2 3 1 - 2 2 3 9 .

[7]黄守道，刘婷，欧阳红林，等.基于槽口偏移的永磁电机齿槽转矩削弱方法[J].电工技术学报，2013，28（3）：99-106.

[8] WANG K，LIANG Y P，WANG D M，et al.Cogging torque reduc⁃tion by eccentric structure of teeth in external rotor permanent magnet synchronous motors[J].IET Electric Power Applications， 2 0 1 9 ，1 3（ 1 ）：5 7 - 6 3 .

[9]安治国，刘高朋，高尉.转子辅助槽对IPM性能的影响研究[J].重庆交通大学学报（自然科学版），2019，38（4）：123-127.

[10] ZHAO J，HU H Z，LIU X D，et al.Influence of edge permanent- magnet shape on the performance of an arc-linear permanent- magnet synchronous machine[J].IEEE Transactions on Magnet⁃i c s ，2 0 1 5 ，5 1（ 1 1 ）：1 - 4 .

[11]李松生，杨焕钊，何国庆，等.基于齿槽转矩的永磁同步电主轴动态性能试验研究[J].振动与冲击，2020，39（20）：216-222.

[12] WU J，ZHANG H R，HE R，et al. A mechatronic brake booster for electric vehicles：design，control，and experiment[J].IEEE Transactions on Vehicular Technology，2020，69（7）：7040-7053.

[13]李彬，马同凯，贺媛，等.弧形永磁同步电机定位转矩补偿控制方法研究[J].中国电机工程学报，2020，40（24）：7909-7917.

[14]李辉，许艮，杨超，等.汽车天窗电机设计参数对齿槽转矩影响及优化[J].电机与控制学报，2012，16（1）：18-22.