齿槽转矩解析与低减手法
齿槽转矩解析与低减手法的理论基础
什么是齿槽转矩
齿槽转矩是指电动机旋转时的卡滞感吗?
正是。这是永久磁石与槽口之间的磁气力产生的脉动转矩。直接影响电动汽车的静音性和机器人的定位精度。
不通电的情况下也会产生转矩吗? 这很奇怪呢。
完全正确。即使线圈无通电,永久磁石的磁力线通过槽口时,磁路阻抗(磁阻)随转角变化。磁力线倾向于选择"阻力小的路径",因此当磁石经过槽口时会出现"卡住"的力。这就是齿槽转矩。
像手摇发电机用手转动时的"卡卡卡"的感觉吗?
非常好的比喻。在电动汽车的转向电动机中,如果齿槽转矩过大,驾驶员会感到方向盘的振动。在工业机器人中,它会降低定位精度。在直线电动机中,甚至会影响停止精度。因此"如何最小化齿槽转矩"成为电动机设计者的重要课题。
发生机制
能详细解释一下为什么会产生脉动吗?气隙的磁通与这有关吗?
齿槽转矩的本质是气隙中的磁能 $W_{mag}$ 随转角 $\theta$ 的变化而周期性变动。永久磁石产生的起磁力分布与槽口开口产生的磁阻分布的"乘积"决定了磁通密度。
简单来说,磁阻在槽口处急剧下降。每当磁石通过这个"孔洞"时,能量都会变化,其对转角的微分就表现为齿槽转矩。
其中 $i=0$ 表示无通电状态。这意味着齿槽转矩完全由磁石与槽口的几何关系决定。
磁石经过每个槽口时都会出现"卡住→释放"的循环,对吧?
原理上是这样,但更准确地说,单个磁石可能同时跨越多个槽口。因此齿槽转矩的周期不是简单地"槽数分",而是由极数与槽数的最小公倍数(LCM)决定的。
支配方程与傅里叶级数展开
用公式如何表达齿槽转矩?
齿槽转矩是周期函数,自然地用傅里叶级数展开:
其中 $N_{LCM} = \text{LCM}(N_s,\, 2p)$(槽数 $N_s$ 与极数 $2p$ 的最小公倍数),$T_n$ 是第 $n$ 次谐波的振幅,$\phi_n$ 是相位。
LCM越大,基频就越高,振幅会变小吗?
完全正确。通常高次谐波的振幅衰减非常快。因此LCM较大的极槽组合,齿槽转矩会明显更小。
另一个重要公式是利用Maxwell应力张量计算转矩。在气隙封闭曲面 $S$ 上积分:
其中 $B_r$、$B_\theta$ 分别是径向与切向磁通密度分量,$r$ 是气隙半径。在FEM解析中,通过数值积分来计算齿槽转矩。
极数与槽数的LCM理论
能给出一些具体的LCM例子吗? 实际上差异有多大?
让我列举一些典型的组合:
| 极数 2p | 槽数 Ns | LCM | 齿槽转矩特征 |
|---|---|---|---|
| 4 | 6 | 12 | 较大(整数槽距) |
| 8 | 12 | 24 | 中等 |
| 10 | 12 | 60 | 非常小(分数槽距) |
| 14 | 12 | 84 | 极小 |
| 8 | 9 | 72 | 小(分数槽集中绕) |
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解析定量判断。
仅对槽口进行"倒角"就能把齿槽转矩降低一半
在齿槽转矩低减的实用手法中,性价比最高的是"对定子槽口进行倒角(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$。二维解析的支配方程为:
其中 $\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/\text{GCD}(N_s, 2p)$ 模型
- 定义材料:铁芯采用B-H曲线(硅钢板35H300等),磁石设定剩余磁通密度 $B_r$ 和保磁力 $H_c$
- 网格生成:气隙至少4~6层网格。特别加密槽口开口部分
- 边界条件:设置周期对称边界(奇偶性)和外围气球边界
- 参数化扫描:转子旋转一个齿槽周期($360°/N_{LCM}$),各角度进行静磁场解析
- 后处理:提取转矩波形并进行FFT分析,识别高次谐波成分
只需计算一个齿槽周期就够了。不必跑全圆周。
对。但在评估低减手法时,最好计算一个电气周期($360°/p$)范围,同时检查齿槽转矩和转矩脉动。只减少齿槽转矩而增加转矩脉动就得不偿失了。
气隙网格设计
气隙网格这么重要,具体怎样设计?
齿槽转矩解析中最关键的网格设计要点:
| 区域 | 推荐单元尺寸 | 层数 | 原因 |
|---|---|---|---|
| 气隙 | 间隙长的1/4~1/6 | 4~6层 | 保证MST法积分精度 |
| 槽口开口 | 开口宽度的1/5以下 | 3~4层 | 捕捉磁阻急变 |
| 磁石表面 | 磁石厚的1/3 | 3层以上 | 捕捉表面磁通变动 |
| 齿本体 | 齿宽的1/3 | 2~3层 | 评估饱和影响 |
| 轭部、后铁 | 可粗化 | 2层 | 提高计算效率 |
气隙0.5mm的话需要0.1mm以下的单元。这太细了!
就是这样。网格收敛确认是必须的,至少用3个水平(如3层→5层→7层)验证齿槽转矩峰值在5%内收敛。不足的话数值噪声会混入转矩中,看不清真实值。
斜极解析的实现
用FEM计算斜极的效果需要三维解析吗?
三维最准确,但实务中常用多片二维法。将轴向分为 $M$ 片,各片以不同的斜极角进行二维解析,结果平均:
$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是用来指导设计方向的,最终品质保证还是要靠实测。
"空气区域不能缩小"是电磁场解析的铁则
初学者在做电动机电磁场解析时常会问"为什么要网格化空气?"。答案很简单——磁力线会散射到铁芯外面。如果解析区域限制在铁芯周围,无处可去的磁力线会在边界反射,导致实际不存在的磁通集中。最低要确保定子外径的1.5~2倍半径作为空气区域。就像在狭小房间内扔篮球会被墙壁反复弹回一样,磁力线也会出现虚假的反射现象。
齿槽转矩解析与低减手法的软件对比
主要工具特点
有哪些商用软件可以进行齿槽转矩解析?
列举一下电动机设计中常用的主要电磁场解析工具:
| 工具名 | 开发商 | 特点 | 齿槽转矩解析 |
|---|---|---|---|
| JMAG-Designer | JSOL(日本) | 电动机设计专用。日本汽车业广泛使用 | 内置斜极功能,法线力法高精度 |
| Ansys Maxwell | Ansys Inc.(美国) | 汎用电磁场解析。Electronics Desktop集成 | 虚功法、自适应网格 |
| COMSOL Multiphysics | COMSOL AB(瑞典) | 多物理场集成。PDE自定义灵活 | MST/VW可选,易于多物理耦合 |
| Motor-CAD | Ansys Inc. | 电动机专用设计GUI。初期设计特化 | 简化计算(内置FEM或解析式) |
| FEMM | David Meeker(开源) | 二维静磁/交流磁场。免费 | VW法。Lua/Python脚本自动化 |
免费软件FEMM也能计算齿槽转矩?
可以的。但FEMM只支持二维,没有内置斜极功能,需要自己写Python脚本实现多片法。学生研究或初期评估足够,但产品设计用JMAG或Maxwell更高效。
齿槽转矩功能对比
各工具的齿槽转矩功能具体怎样对比?
| 功能 | JMAG | Maxwell | COMSOL | FEMM |
|---|---|---|---|---|
| 二维齿槽解析 | ○ | ○ | ○ | ○ |
| 三维齿槽解析 | ○ | ○ | ○ | × |
| 斜极解析(内置) | ○ | △(仅3D) | △(仅3D) | × |
| 多片二维法 | ○(自动) | △(脚本) | △(脚本) | △(脚本) |
| 参数化优化 | ○ | ○ | ○ | △(Lua/Python) |
| FFT转矩分析 | ○ | ○ | △(后处理) | × |
| 振动噪声耦合 | ○ | ○(Workbench) | ○ | × |
| 拓扑优化 | ○ | △ | △ | × |
工具选型指导
最后,应该怎样选工具? 考虑预算的话?
按应用领域分类选择:
- 汽车制造商、电机Tier1的电动机设计部:基本是JMAG。日本汽车OEM和Tier1的70%以上采用。斜极解析和损耗计算精度高
- 海外企业、汎用电磁场解析:Ansys Maxwell。与Workbench集成,结构-电磁耦合方便。社内有其他Ansys产品时协同效应大
- 研究机构、多物理场耦合:COMSOL。一个平台内完成振动、声学、热等耦合。自定义PDE适合验证新方法
- 学生、初期评估、预算有限:FEMM入门,之后升迁到JMAG或Maxwell的学生版
为什么同一电动机用不同工具计算齿槽转矩结果会有差异
用多个商用工具对同一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(振动噪声)问题。从简单的二维模型开始,亲手摸索一遍是最好的学习。
有帮助
更详细
误