齿槽转矩解析与低减手法

分类:电磁场解析 > 电动机设计 | 更新 2026-04-11
Cogging torque analysis showing magnetic flux distribution between permanent magnets and stator slots
永磁电动机中齿槽转矩的磁通密度分布与低减手法的概念图

齿槽转矩解析与低减手法的理论基础

什么是齿槽转矩

🧑‍🎓

齿槽转矩是指电动机旋转时的卡滞感吗?

🎓

正是。这是永久磁石与槽口之间的磁气力产生的脉动转矩。直接影响电动汽车的静音性和机器人的定位精度。

🧑‍🎓

不通电的情况下也会产生转矩吗? 这很奇怪呢。

🎓

完全正确。即使线圈无通电,永久磁石的磁力线通过槽口时,磁路阻抗(磁阻)随转角变化。磁力线倾向于选择"阻力小的路径",因此当磁石经过槽口时会出现"卡住"的力。这就是齿槽转矩。

🧑‍🎓

像手摇发电机用手转动时的"卡卡卡"的感觉吗?

🎓

非常好的比喻。在电动汽车的转向电动机中,如果齿槽转矩过大,驾驶员会感到方向盘的振动。在工业机器人中,它会降低定位精度。在直线电动机中,甚至会影响停止精度。因此"如何最小化齿槽转矩"成为电动机设计者的重要课题。

发生机制

🧑‍🎓

能详细解释一下为什么会产生脉动吗?气隙的磁通与这有关吗?

🎓

齿槽转矩的本质是气隙中的磁能 $W_{mag}$ 随转角 $\theta$ 的变化而周期性变动。永久磁石产生的起磁力分布与槽口开口产生的磁阻分布的"乘积"决定了磁通密度。

🎓

简单来说,磁阻在槽口处急剧下降。每当磁石通过这个"孔洞"时,能量都会变化,其对转角的微分就表现为齿槽转矩。

$$ T_{cog}(\theta) = -\frac{\partial W_{mag}(\theta)}{\partial \theta}\bigg|_{i=0} $$

其中 $i=0$ 表示无通电状态。这意味着齿槽转矩完全由磁石与槽口的几何关系决定。

🧑‍🎓

磁石经过每个槽口时都会出现"卡住→释放"的循环,对吧?

🎓

原理上是这样,但更准确地说,单个磁石可能同时跨越多个槽口。因此齿槽转矩的周期不是简单地"槽数分",而是由极数与槽数的最小公倍数(LCM)决定的。

支配方程与傅里叶级数展开

🧑‍🎓

用公式如何表达齿槽转矩?

🎓

齿槽转矩是周期函数,自然地用傅里叶级数展开:

$$ T_{cog}(\theta) = \sum_{n=1}^{\infty} T_n \sin\!\bigl(n \cdot N_{LCM} \cdot \theta + \phi_n\bigr) $$

其中 $N_{LCM} = \text{LCM}(N_s,\, 2p)$(槽数 $N_s$ 与极数 $2p$ 的最小公倍数),$T_n$ 是第 $n$ 次谐波的振幅,$\phi_n$ 是相位。

🧑‍🎓

LCM越大,基频就越高,振幅会变小吗?

🎓

完全正确。通常高次谐波的振幅衰减非常快。因此LCM较大的极槽组合,齿槽转矩会明显更小。

🎓

另一个重要公式是利用Maxwell应力张量计算转矩。在气隙封闭曲面 $S$ 上积分:

$$ T = \frac{r}{\mu_0} \oint_S B_r \cdot B_\theta \, dS $$

其中 $B_r$、$B_\theta$ 分别是径向与切向磁通密度分量,$r$ 是气隙半径。在FEM解析中,通过数值积分来计算齿槽转矩。

极数与槽数的LCM理论

🧑‍🎓

能给出一些具体的LCM例子吗? 实际上差异有多大?

🎓

让我列举一些典型的组合:

极数 2p槽数 NsLCM齿槽转矩特征
4612较大(整数槽距)
81224中等
101260非常小(分数槽距)
141284极小
8972小(分数槽集中绕)
🧑‍🎓

10极12槽在电动汽车驱动电动机中很常见,原来是因为齿槽转矩小啊!

🎓

正是!LCM=60意味着每转360度产生60次细微脉动,分散后基本上感觉不到。但分数槽也有缺点,比如绕线系数会下降。"极数与槽数的选择在设计阶段就决定了七成的齿槽转矩问题",这是行业共识。

低减手法体系

🧑‍🎓

除了改变极槽,还有哪些降低齿槽转矩的技术?

🎓

主要的低减手法体系如下:

手法原理低减效果副作用
斜极(倾斜)转子或定子沿轴向倾斜50~90%平均转矩略降(数%)
分数槽距增大LCM进行高次谐波分散60~95%绕线系数可能下降
磁石形状倒角缓和磁石端部磁阻变化30~60%基本无
磁石偏心配置改变磁石圆弧中心位置20~50%设计复杂化
槽口宽度优化缩小磁阻变动范围20~40%漏磁增加
分段斜极沿轴向分段改变斜极角60~85%层叠工序复杂
不等间距磁石改变相邻磁石夹角30~70%影响转矩脉动
🧑‍🎓

仅用斜极就能减少90%? 这太厉害了!

🎓

理论上,一个槽距(360°/Ns)的斜极能使基波齿槽转矩趋零。但这会让平均转矩按sinc函数衰减,所以实务中通常采用半槽距的斜极。"采用哪种手法、达到什么程度"需要通过FEM解析定量判断。

Coffee Break 闲谈

仅对槽口进行"倒角"就能把齿槽转矩降低一半

在齿槽转矩低减的实用手法中,性价比最高的是"对定子槽口进行倒角(chamfering)"。只需将槽口边缘倾斜打磨,就能缓和磁石经过时的急剧磁阻变化,使齿槽转矩降低30~50%。制造成本几乎不增加,磁通量也基本不变。通过FEM参数化扫描槽口倒角角度和深度来找最优值是业界常见做法。"实施复杂多变量优化前,先试试倒角"是设计现场的不成文规则。

齿槽转矩解析与低减手法的数值计算手法

转矩计算的数值手法

🧑‍🎓

用FEM计算齿槽转矩时,具体使用哪些方法?

🎓

主要有三种方法。它们在精度和计算成本间有不同的权衡:

手法原理精度网格敏感度
虚功法(VW法)从FEM场计算 $T = -\partial W_{mag}/\partial \theta$低(鲁棒)
Maxwell应力张量法(MST)在气隙表面积分 $B_r \cdot B_\theta$中~高高(对积分面敏感)
荷载张量法(Arkkio法)在气隙体积进行体积积分中等
🧑‍🎓

MST法对网格敏感,那实务中用哪种最安全?

🎓

JMAG内部采用改进的法线力法,精度很高。Ansys Maxwell默认使用虚功法。COMSOL则让用户选择。如使用MST法,应将积分面设在气隙中心,并用多个积分面的结果平均,这是最佳实践。

磁矢势法

🧑‍🎓

齿槽转矩解析的基础方程是什么?

🎓

电动机的电磁场解析采用矢量磁位 $\mathbf{A}$。通过定义 $\mathbf{B} = \nabla \times \mathbf{A}$,可自动满足磁通守恒 $\nabla \cdot \mathbf{B} = 0$。二维解析的支配方程为:

$$ \nabla \times \left(\frac{1}{\mu} \nabla \times \mathbf{A}\right) = \mathbf{J}_s + \nabla \times \left(\frac{\mathbf{B}_r}{\mu_0}\right) $$

其中 $\mathbf{J}_s$ 是电流密度(齿槽转矩计算中 $=0$),$\mathbf{B}_r$ 是永久磁石的剩余磁通密度。二维问题可简化为 $A_z$ 标量,计算量较少。

🧑‍🎓

齿槽转矩计算是无电流状态,所以右边的电流项为零,只有磁石项,对吗?

🎓

正是。因此齿槽转矩解析相对简单,是纯静磁场问题。但考虑铁芯的非线性B-H特性后,需要牛顿-拉弗逊迭代。而且转子要转过多个位置角进行数十步计算,所以计算时间还是不少的。

非线性B-H特性的处理

🧑‍🎓

铁芯的磁饱和会影响齿槽转矩吗?

🎓

影响很大。齿尖处的磁通容易集中而饱和,饱和后导磁率急剧下降。饱和区域变得"接近空气",相当于槽口变宽,齿槽转矩反而可能增大。

🎓

所以齿槽转矩解析必须考虑非线性B-H曲线。收敛判定通常为残差相对值 $||R||/||R_0|| < 10^{-4}$,一般10~30次迭代即可收敛。

转子网格处理

🧑‍🎓

转子旋转时,网格怎么处理? 每次旋转都重新生成吗?

🎓

重新生成网格计算成本很高,所以实务中标准做法是在气隙中设置滑动边界(旋转带)

  • 滑动网格法:在气隙中建立圆环形网格带,转子和定子的网格独立保持。每次旋转时更新带上节点的对应关系
  • 网格变形法:根据旋转角度改变气隙单元形状。定期重新网格化
  • 多项式插值法:用多项式在带上插值磁位,连接转子/定子间的节点
🧑‍🎓

不同工具实现滑动边界的方式不同,所以同一模型结果会有差异吗?

🎓

完全同意。特别是齿槽转矩这种微小量,对接合精度很敏感。JMAG采用内部的AZ框架法,精度较高。Maxwell需要充分加密带上的网格才能避免数值噪声混入。

角度刻度与周期性的应用

🧑‍🎓

转子转过多少度采样一次合适? 1度刻度的话就是360步,太多了。

🎓

利用齿槽转矩的周期对称性可大幅减少计算量:

  • 周期对称性:利用 $\text{GCD}(N_s, 2p)$ 的对称性,解析域可缩小到 $1/\text{GCD}$ 倍。例如10极12槽,GCD=2,只需半模型
  • 齿槽周期:一个齿槽周期为 $360°/N_{LCM}$,只需计算这个范围。10极12槽只需6度范围
  • 角度分割:至少20~30等分齿槽周期。10极12槽需0.2~0.3度的刻度
🧑‍🎓

6度范围用0.25度刻度,只需24步。二维解析的话几分钟就完了吧。

🎓

二维的话确实可以。三维考虑斜极时计算量大增,这时常用"多片二维切片法"——沿轴向分成多片,各片做二维解析后求和平均。这是很实用的折中方案。

齿槽转矩解析与低减手法的实务应用

齿槽转矩解析步骤

🧑‍🎓

从零开始进行齿槽转矩解析的流程是什么?

🎓

标准的齿槽转矩解析工作流如下:

  1. 建立模型:建立二维截面模型。利用周期对称性缩小至 $1/\text{GCD}(N_s, 2p)$ 模型
  2. 定义材料:铁芯采用B-H曲线(硅钢板35H300等),磁石设定剩余磁通密度 $B_r$ 和保磁力 $H_c$
  3. 网格生成:气隙至少4~6层网格。特别加密槽口开口部分
  4. 边界条件:设置周期对称边界(奇偶性)和外围气球边界
  5. 参数化扫描:转子旋转一个齿槽周期($360°/N_{LCM}$),各角度进行静磁场解析
  6. 后处理:提取转矩波形并进行FFT分析,识别高次谐波成分
🧑‍🎓

只需计算一个齿槽周期就够了。不必跑全圆周。

🎓

对。但在评估低减手法时,最好计算一个电气周期($360°/p$)范围,同时检查齿槽转矩和转矩脉动。只减少齿槽转矩而增加转矩脉动就得不偿失了。

气隙网格设计

🧑‍🎓

气隙网格这么重要,具体怎样设计?

🎓

齿槽转矩解析中最关键的网格设计要点:

区域推荐单元尺寸层数原因
气隙间隙长的1/4~1/64~6层保证MST法积分精度
槽口开口开口宽度的1/5以下3~4层捕捉磁阻急变
磁石表面磁石厚的1/33层以上捕捉表面磁通变动
齿本体齿宽的1/32~3层评估饱和影响
轭部、后铁可粗化2层提高计算效率
🧑‍🎓

气隙0.5mm的话需要0.1mm以下的单元。这太细了!

🎓

就是这样。网格收敛确认是必须的,至少用3个水平(如3层→5层→7层)验证齿槽转矩峰值在5%内收敛。不足的话数值噪声会混入转矩中,看不清真实值。

斜极解析的实现

🧑‍🎓

用FEM计算斜极的效果需要三维解析吗?

🎓

三维最准确,但实务中常用多片二维法。将轴向分为 $M$ 片,各片以不同的斜极角进行二维解析,结果平均:

$$ T_{cog,skew}(\theta) = \frac{1}{M} \sum_{k=0}^{M-1} T_{cog,2D}\!\left(\theta + k \cdot \frac{\theta_{sk}}{M}\right) $$

$M = 5$~$10$ 片就能获得足够精度。JMAG的"斜极解析"功能内置了多片法,二维模型即可自动计算,很方便。

🧑‍🎓

分段斜极(磁石分段改变角度)的情况也用同样的方法?

🎓

基本相同。三段分段斜极可用 $M=3$ 直接计算。但分段斜极的轴向边界效应(漏磁)较为明显,某些情况需要三维解析才能准确评估。

磁石形状优化

🧑‍🎓

改变磁石形状来降低齿槽转矩怎样实施?

🎓

磁石形状的参数化优化是最实用的齿槽转矩低减手法之一。主要参数:

  • 磁石圆弧角(embrace):磁石宽/极距比。在0.7~0.9范围内扫描找最优值
  • 端部倒角角度和深度:对磁石端部斜切以缓和磁阻变化
  • 磁石表面曲率:从平面型到波浪型。改变气隙有效长度的周向分布
  • 磁石中心偏心量:改变磁石圆弧中心相对转子中心的位置
🧑‍🎓

比如embrace从0.7变到0.9,齿槽转矩如何变化?

🎓

通常呈V字形或波形变化。某个特定值(如0.833 = 5/6这样的分数)齿槽转矩最小。最优值随极槽组合而异,用FEM参数扫描最可靠。实务上先用0.05步长做粗搜索掌握趋势,再用0.01步长细化有前景的范围,这样效率最高。

解析结果验证

🧑‍🎓

计算结果是否正确,怎样确认?

🎓

齿槽转矩解析的合理性验证有以下检查点:

  • 波形对称性:齿槽转矩波形应为原点对称(平均值为零)。若平均值不为零,模型有问题
  • 周期整合性:波形峰值数应与每转 $N_{LCM}$ 次相符
  • 网格收敛:三个水平以上的网格,峰值在5%内收敛
  • 计算方法对比:虚功法和MST法结果偏差超过20%需警惕
  • 实验对比:用分辨率0.001 N·m以下的转矩传感器测量对比
🧑‍🎓

实测与计算比较时精度能达到多少?

🎓

网格充分细化的二维解析,偏差通常在10~30%以内。加上三维模型和端部效应考虑,能达到10%以内。但量产品的齿槽转矩离散度本身就很大——着磁不均、层叠精度等因素会引入20~50%的个体差异。所以FEM是用来指导设计方向的,最终品质保证还是要靠实测。

Coffee Break 闲谈

"空气区域不能缩小"是电磁场解析的铁则

初学者在做电动机电磁场解析时常会问"为什么要网格化空气?"。答案很简单——磁力线会散射到铁芯外面。如果解析区域限制在铁芯周围,无处可去的磁力线会在边界反射,导致实际不存在的磁通集中。最低要确保定子外径的1.5~2倍半径作为空气区域。就像在狭小房间内扔篮球会被墙壁反复弹回一样,磁力线也会出现虚假的反射现象。

齿槽转矩解析与低减手法的软件对比

主要工具特点

🧑‍🎓

有哪些商用软件可以进行齿槽转矩解析?

🎓

列举一下电动机设计中常用的主要电磁场解析工具:

工具名开发商特点齿槽转矩解析
JMAG-DesignerJSOL(日本)电动机设计专用。日本汽车业广泛使用内置斜极功能,法线力法高精度
Ansys MaxwellAnsys Inc.(美国)汎用电磁场解析。Electronics Desktop集成虚功法、自适应网格
COMSOL MultiphysicsCOMSOL AB(瑞典)多物理场集成。PDE自定义灵活MST/VW可选,易于多物理耦合
Motor-CADAnsys Inc.电动机专用设计GUI。初期设计特化简化计算(内置FEM或解析式)
FEMMDavid Meeker(开源)二维静磁/交流磁场。免费VW法。Lua/Python脚本自动化
🧑‍🎓

免费软件FEMM也能计算齿槽转矩?

🎓

可以的。但FEMM只支持二维,没有内置斜极功能,需要自己写Python脚本实现多片法。学生研究或初期评估足够,但产品设计用JMAG或Maxwell更高效。

齿槽转矩功能对比

🧑‍🎓

各工具的齿槽转矩功能具体怎样对比?

🎓
功能JMAGMaxwellCOMSOLFEMM
二维齿槽解析
三维齿槽解析×
斜极解析(内置)△(仅3D)△(仅3D)×
多片二维法○(自动)△(脚本)△(脚本)△(脚本)
参数化优化△(Lua/Python)
FFT转矩分析△(后处理)×
振动噪声耦合○(Workbench)×
拓扑优化×

工具选型指导

🧑‍🎓

最后,应该怎样选工具? 考虑预算的话?

🎓

按应用领域分类选择:

  • 汽车制造商、电机Tier1的电动机设计部:基本是JMAG。日本汽车OEM和Tier1的70%以上采用。斜极解析和损耗计算精度高
  • 海外企业、汎用电磁场解析:Ansys Maxwell。与Workbench集成,结构-电磁耦合方便。社内有其他Ansys产品时协同效应大
  • 研究机构、多物理场耦合:COMSOL。一个平台内完成振动、声学、热等耦合。自定义PDE适合验证新方法
  • 学生、初期评估、预算有限:FEMM入门,之后升迁到JMAG或Maxwell的学生版
Coffee Break 闲谈

为什么同一电动机用不同工具计算齿槽转矩结果会有差异

用多个商用工具对同一IPM电动机计算齿槽转矩,结果往往相差几个到十几个百分点。主要原因是气隙网格细度和滑动网格(旋转边界接合)的实现方式差异。JMAG采用改进的法线力法,空隙磁通密度到齿槽转矩的转换精度很高。Ansys Maxwell基于虚功法,但根据参数设置精度会变化。为避免以工具差异为借口,应建立社内的标准,对标实际测量数据。

齿槽转矩解析与低减手法的故障排除

齿槽转矩与实测不符

🧑‍🎓

FEM计算的齿槽转矩和实测差很大,应该从哪开始检查?

🎓

这是齿槽转矩解析中最常见的咨询。按优先度列出检查项:

检查项常见原因对策
气隙网格层数不足(2层以下)增加到最低4~6层,确认收敛
B-H曲线采用线性材料 / 曲线精度差使用厂商实测B-H数据
气隙长公差±0.05mm影响大做上限、下限、标准三工况
磁石特性Br离散度±5%敏感度分析评估影响度
斜极考虑二维未考虑斜极用多片法反映三维效果
制造离散着磁不均、打冲精度差认识到FEM是理想形状,进行离散评估
🧑‍🎓

气隙仅差0.05mm就会让齿槽转矩变化这么大?

🎓

会的。气隙长与磁阻成反比,0.5mm改为0.45mm,齿槽转矩可能增加10~20%。和实测对标时,必须对公差范围进行敏感度分析。

数值噪声与虚假高次谐波

🧑‍🎓

FFT分析中出现理论上不该有的高次谐波,这是什么?

🎓

数值噪声引起的"虚假高次谐波"。主因有三点:

  • 滑动边界网格不匹配:转子和定子网格节点数不同,接合部出现数值脉动。解决办法是增加旋转带上的节点数
  • 角度刻度不足:根据采样定理,最高次谐波需要2倍以上采样。刻度不足会出现走样,虚假低次成分
  • 非线性收敛不充分:各角度步的收敛残差不均,导致角度间的偏差变成噪声
🧑‍🎓

怎样区分虚假高次与真实高次?

🎓

网格加倍细化后重新计算,如果振幅大幅变化就是数值噪声。物理成分对网格不敏感。另外,齿槽转矩理论上仅含 $n \cdot N_{LCM}$ 次($n=1,2,3,...$)分量,其他次数的成分肯定是虚假的。

斜极效果不显著

🧑‍🎓

加了斜极后齿槽转矩没有降多少,是配置错了吗?

🎓

常见的问题及检查:

  • 斜极角的单位:机械角和电气角混淆。JMAG是机械角输入,Maxwell的设置随工程不同而异
  • 切片数不足:多片法用2~3片不够精确。最少需5片以上
  • 二维未激活斜极功能:二维解析本身不反映斜极,需明确启用"多片斜极分析"或"斜极选项"
  • 斜极方向错误:定子和转子斜极效果相当,但与不等间距磁石联合使用时方向很重要

非线性收敛失败

🧑‍🎓

特定转角处会不收敛,为什么?

🎓

磁石正对槽口时,齿尖处磁通密度急变,在非线性B-H曲线的"饱和过渡点"产生振荡,牛顿-拉弗逊法不稳定。

对策:

  • 前一步的解作为初值(连续扫描模式)
  • 收敛残差设定为 $10^{-3}$~$10^{-4}$(不能过严)
  • B-H曲线高磁通部分平滑插值(避免尖锐折点)
  • 阻尼系数(欠松弛因子)设为0.5~0.8
🧑‍🎓

通过这次讨论,齿槽转矩解析的整体框架清楚多了。理论、网格、验证的三位一体很关键!

🎓

完全同意。齿槽转矩虽然"小",但很"讨厌"。网格详细度和结果验证的马虎,会产生完全无意义的数字。但做对的话,在设计阶段就能把齿槽转矩最小化,从根本上防止NVH(振动噪声)问题。从简单的二维模型开始,亲手摸索一遍是最好的学习。

相关模拟器

通过此领域的交互式模拟器来体验理论

模拟器列表

相关领域

多物理耦合结构分析热分析
本文评分
感谢您的反馈!
对我
有帮助
需要
更详细
内容有
对我有帮助
0
需要更详细
0
内容有误
0
撰写:NovaSolver贡献者
匿名工程师与人工智能 — 网站地图
查看个人资料