Maxwell 2D FEA 仿真助手#

专注电机工程领域的 Maxwell 2D 仿真，覆盖静磁场、瞬态场、参数化扫描等典型仿真场景。

仿真类型与适用场景#

仿真类型求解器适用场景静磁场（Magnetostatic）Maxwell 2D永磁体工作点、气隙磁密、空载反电动势瞬态场（Transient）Maxwell 2D力矩特性、负载响应、齿槽转矩、纹波涡流场（Eddy Current）Maxwell 2D高速电机套筒损耗、制动盘静电力场（Electrostatic）Maxwell 2D电容传感器、静电电机

仿真建模标准流程#

1. 几何建模#

常用参数（单位：mm）

参数说明示例定子外径Dso90定子内径Dsi54转子外径Dri53转子内径Dri_shaft16气隙长度g0.5铁心长度L50极数2p8槽数Q36

建模原则

• 优先使用 1/4 或 1/2 对称模型减少计算量
• 定子外圆和转子内圆建议加空气包域（边界层外扩≥20%模型尺寸）
• 斜槽：用多截面或斜极方式建模

2. 材料定义#

定子/转子铁心

text

材料类型：Silicon Steel（硅钢片）
常用牌号：
  - 35YT310：B_H曲线 + 损耗曲线
  - M19_29Gauss：美标对应
  - DW310-35：国产牌号
设置项：BH曲线、层压方向（0°/90°叠压）、电导率

永磁体

text

材料类型：Nonlinear Magnetic（非线性磁性）
牌号：N42SH、N35UH、35EH 等
关键参数：
  - 剩磁 Br（T）：如 N42SH → Br≈1.25T
  - 矫顽力 Hc（kA/m）：如 N42SH → Hc≈955kA/m
  - 温度系数：Br温度系数约 -0.1%/℃
设置：各向异性充磁方向（Radial / Circumferential）

铜线/绕组

text

材料：Copper（退火铜）
电导率：5.8e7 S/m（20℃）
温度修正：每℃约 +0.4%

3. 边界条件#

边界类型设置方法适用场景Balloon边界Set balloon开放边界（默认）Vector Potential = 0主从边界副边界周期对称模型狄利克雷边界指定 A=0建模对称轴

典型设置（1/4对称模型）

text

上边界（α=0°）：Vector Potential = 0
右边界（r=R_max）：Balloon（或 A=0）
模型外空气包外圆：Balloon边界

4. 激励设置#

永磁体激励

text

方法1：等效面电流
  面电流密度 Jc = Hc / δ（δ为等效厚度）

方法2：直接赋材料
  材料类型选永磁体牌号
  充磁方向：Radial（径向充磁）
    8极：每极充磁方向交替（NSNS...）

绕组激励

text

绕组类型：Stranded（绞线，忽略涡流）
匝数设置：每个线圈单元的串联匝数
电流赋值：瞬态场用随时间变化的电流表达式
  示例：I_phaseA = Imax * sin(2*pi*freq*time)

绕组建模方式

text

方案A：绕组 band 方式
  定子槽内画 coil 区域 → 赋电流密度 J = N*I/SlotArea
  槽内填充率（槽满率）通常 0.45~0.55

方案B：外电路耦合（External Circuit）
  定义绕组端子 → 连接 voltage source
  适合反电动势测试

5. 求解设置#

静磁场求解器

text

求解类型：Magnetostatic
收敛标准：Energy Error < 0.001（或 0.0001 高精度）
自适应求解：开启， refinement per pass = 5
网格加密：永磁体/气隙区域加密
  气隙内侧：Maximum Length < 0.1mm
  磁钢区域：Maximum Length < 0.2mm

瞬态场求解器

text

时间步长：Tstep = 1/（freq × 360×P）/ 20
  8极电机 @ 3000rpm @ 50Hz → Tstep = 1/(50×360×4)/20 ≈ 0.00007s
  最小步长：建议 < 周期/180
求解时间：至少 1 个完整电周期
后处理：提取 T_N（力矩）和 EMF（反电动势）

参数化扫描

text

扫描变量：电流 I、角度 theta、永磁体厚度 hm、槽深 h
扫描方式：Optimetrics → 添加 Parametric Sweep
  变化范围：如 I = 0~20A，步长 2A
结果提取：Map Fields 或 Table

典型仿真案例#

案例1：空载反电动势（静磁场）#

text

目标：计算 8极36槽电机 3000rpm 时的空载反电动势
方法：静磁场 + 角度参数化

设置步骤：
1. 定义转子角度参数 theta（0°~360°，步长10°）
2. 设置静磁场求解
3. 提取气隙磁通密度 B(theta)
4. 后处理计算：E = N * dΦ/dt = N * ω * dΦ/dθ
5. 或用瞬态场直接测量线间电压波形

案例2：齿槽转矩（Cogging Torque）#

text

目标：计算无电流激励下的齿槽转矩
方法：瞬态场，绕组开路

设置步骤：
1. 绕组电流设为 0（开路）
2. 转速设为额定转速（如 3000rpm）
3. 求解 1 个完整机械周期
4. 提取 Torque，转子角度为横坐标
结果：齿槽转矩峰值通常 1%~3% 额定转矩

案例3：额定力矩（瞬态场）#

text

目标：计算 FOC 控制下的额定力矩
方法：瞬态场 + 外电路

设置步骤：
1. 定义 d轴电流 Id = 0 控制
2. q轴电流 Iq = 额定相电流峰值
3. 设置初始角度为 d轴对齐位置
4. 求解 1 个电周期后取平均值

网格质量评估#

区域建议网格尺寸优先级气隙≤0.08mm最高永磁体≤0.15mm高槽口/槽楔≤0.1mm高定子/转子轭部≤0.5mm中空气包外层≤2mm低

网格质量检查

• 三角形单元（Default），三角形比例 Aspect Ratio < 3:1
• 禁用梯形单元（inaccurate）
• 气隙至少有 5 层单元

结果提取与后处理#

常用后处理量

物理量提取方法单位气隙磁密 Bg路径积分（B_vector）T反电动势 EVoltage_transient 探针V电磁转矩 TTorque 探针N·m铁损 P_feLosses → Core LossW铜损 P_cuI²R 计算W效率 ηP_out / (P_out + P_fe + P_cu)%

力矩提取关键点

• 瞬态场 Torque 探针放在转子区域
• 转矩脉动 = (T_max - T_min) / T_avg × 100%
• 关注 6k 次谐波（k=1,2,3...）

常见问题排查#

问题原因解决方案收敛困难气隙网格太粗、边界条件错误加密气隙网格至 <0.1mm力矩为负电流方向错误检查绕组相序磁密过高磁钢过厚/规格选错核实 Br 和 Hc反电动势偏低匝数不够/气隙过大检查绕组连接和 N网格报错模型有间隙/重叠用 Audit 检查几何

参考脚本与资源#

• scripts/maxwell_simulation_checklist.py - 仿真检查清单
• scripts/maxwell_post_processor.py - 仿真结果后处理脚本（齿槽转矩、反电动势、气隙磁密、力矩转速分析）
• references/material_library.md - 常用材料参数库
• references/mesh_guidelines.md - 网格设置规范
• references/simulation_templates.md - 仿真参数配置模板库（静磁场/瞬态场配置、求解器参数、Optimetrics 扫描）

快速使用 maxwell_post_processor.py#

bash

# 齿槽转矩分析（演示数据）
python scripts/maxwell_post_processor.py --mode cogging_torque --pole_pairs 4

# 反电动势谐波分析
python scripts/maxwell_post_processor.py --mode back_emf --fundamental_freq 200 --pole_pairs 4

# 气隙磁密分析
python scripts/maxwell_post_processor.py --mode air_gap_flux

# 力矩转速特性（给定电流矢量）
python scripts/maxwell_post_processor.py --mode torque_speed --pole_pairs 4

# 实际数据导入（CSV格式）
python scripts/maxwell_post_processor.py --mode cogging_torque \
    --file "torque_data.csv" --angle "angle_deg" --torque "torque_nm"

注意： maxwell_post_processor.py 无需 Maxwell 软件，仅用于后处理 CSV/数组数据。 如需与 Maxwell 直接交互（建模、求解控制），需安装 pyedb（ANSYS Python API）。