1. 项目背景与核心价值
盘式电机作为一种特殊结构的永磁同步电机,凭借其轴向磁通、薄饼式结构的特点,在电动汽车、航空航天和精密伺服领域展现出独特优势。这次要解析的"双定单转24槽20极"结构,是盘式电机中典型的分数槽集中绕组设计方案,在转矩密度和齿槽转矩抑制方面具有显著特点。
Maxwell作为电磁场仿真领域的标杆软件,其脚本化建模能力能够大幅提升复杂电机模型的创建效率。这个项目提供的核心代码,实际上封装了从几何建模、材料定义、边界条件设置到求解参数配置的全流程自动化方案。对于电机设计工程师而言,掌握这套代码意味着:
- 建模时间从小时级缩短到分钟级
- 可快速进行参数化扫描研究
- 确保不同设计迭代间的建模一致性
- 为后续优化算法集成提供接口基础
2. 模型架构解析
2.1 拓扑结构特点
24槽20极的配合属于q=0.4的分数槽绕组,这种设计通过选择适当的槽极组合,能够有效抑制齿槽转矩。双定子单转子的对称结构,则实现了磁路的双向利用,其核心优势在于:
- 转矩密度提升约30-50% compared to conventional designs
- 轴向磁通路径缩短,铁损降低
- 自然冷却表面积增加
在Maxwell中建模时需要特别注意:
- 定子铁芯的扇形分段处理
- Halbach阵列永磁体的磁化方向定义
- 气隙区域的网格加密策略
2.2 坐标系与参数化体系
代码中建立了三层参数化系统:
-
几何参数层(G_开头变量):
- G_Stator_OD = 320mm // 定子外径
- G_Airgap = 1.2mm // 单边气隙长度
-
电磁参数层(E_开头变量):
- E_CurrentDensity = 6A/mm²
- E_MagnetGrade = "N38SH"
-
求解控制层(S_开头变量):
- S_MeshSize = 0.5 // 全局网格系数
- S_TimeStep = 0.1ms // 瞬态求解步长
这种命名规范使得代码可维护性大幅提升,也便于进行DOE实验设计。
3. 核心代码模块详解
3.1 几何建模模块
python复制def create_stator_sector():
# 定子扇形创建逻辑
points = [
[0, G_Stator_ID/2],
[G_Tooth_Width/2, G_Stator_ID/2],
[G_Tooth_Width/2 + G_Slot_Open, G_Stator_OD/2],
[0, G_Stator_OD/2]
]
stator_sector = modeler.create_polygon(points)
modeler.rotate(stator_sector, 15) # 360/24=15°
这段代码展示了定子单个扇区的建模方法,关键点在于:
- 采用底部向上(bottom-up)的坐标点构建方式
- 通过参数化控制齿宽(G_Tooth_Width)和槽口宽(G_Slot_Open)
- 旋转复制时使用精确的角度计算(360°/24槽)
3.2 材料属性定义
永磁体定义采用正交各向异性材料模型:
python复制magnet_props = {
"mu_x": 1.05,
"mu_y": 1.05,
"mu_z": 1.0,
"Hc": -890e3, # 矫顽力(A/m)
"Br": 1.25 # 剩磁(T)
}
set_material("N38SH", magnet_props)
特别注意磁化方向的设置技巧:
- 对于20极结构,每个磁极占18°机械角(360°/20)
- 采用Halbach阵列时需要计算每个磁块的方向角
- 通过方向余弦矩阵实现斜极效果
3.3 边界条件设置
轴对称边界和主从边界的正确设置是保证求解精度的关键:
python复制# 周期性边界条件
setup.assign_master_slave(
master="Section_0",
slave="Section_1",
phase_diff=150 # 电角度差(24槽20极的机械角30°对应电角度150°)
)
这里150°的电角度差计算来源于:
- 机械角度:360°/24槽 × 极对数10 = 150°
- 这种设置确保了周期性边界条件的相位连续性
4. 求解器配置技巧
4.1 网格剖分策略
盘式电机的三维场求解需要特殊的网格处理:
python复制mesh_settings = {
"air_region": {"method": "taubin", "size": 0.8},
"stator": {"method": "curvature", "min_size": 0.3},
"coil": {"method": "surface", "layers": 3}
}
经验建议:
- 气隙区域至少划分5层网格
- 永磁体边缘使用曲率自适应加密
- 绕组区域采用边界层网格
4.2 瞬态求解设置
python复制transient_setup = {
"stop_time": "20ms",
"time_step": "0.1ms",
"nonlinear_residual": 1e-4,
"skew_angle": 7.5 # 机械角度斜极
}
关键参数选择依据:
- 时间步长应小于1/20电气周期(本例中电气频率≈333Hz)
- 非线性残差设为1e-4平衡精度与速度
- 斜极角度按齿距角的一半设置(15°/2)
5. 后处理与结果验证
5.1 转矩波形分析
python复制torque = post.get_torque("Moving1")
fft_result = post.fft_analyze(torque)
典型问题诊断:
- 若6次谐波突出 → 检查永磁体偏心
- 若12次谐波明显 → 绕组不对称
- 若高频振荡 → 网格过粗
5.2 磁密分布检查
python复制b_plot = post.plot_field(
"Mag_B",
plane="xy",
view="section"
)
安全阈值参考:
- 定子齿部 < 1.8T
- 轭部 < 1.6T
- 永磁体工作点 > 0.3T
6. 常见问题排查指南
| 现象 | 可能原因 | 解决方案 |
|---|---|---|
| 转矩波动超15% | 1. 槽口宽度过大 2. 极弧系数不当 |
1. 调整G_Slot_Open 2. 优化磁极包角 |
| 求解不收敛 | 1. 材料非线性设置错误 2. 网格畸变 |
1. 检查BH曲线数据 2. 重构关键区域网格 |
| 效率偏低 | 1. 涡流损耗未考虑 2. 绕组端部效应 |
1. 启用涡流求解 2. 添加端部修正系数 |
7. 进阶优化方向
基于现有模型可扩展:
- 多物理场耦合:添加Thermal模块进行温升分析
- 结构优化:参数化齿形和磁极形状
- 控制系统联合仿真:导出Lookup Table供Simulink使用
实际工程应用中,我们曾通过调整极弧系数使转矩脉动从12.3%降至6.8%,同时保持平均转矩不变。这需要编写循环优化脚本,自动修改参数并提取关键指标。