第2章：电磁场理论基础
各位同学，大家好。我是你们的老朋友，一个在电机设计圈摸爬滚打十几年的工程师。今天咱们聊聊电磁场理论，这玩意儿是有限元仿真的“内功心法”。你软件操作得再溜，不懂这个，仿真结果就是空中楼阁。
说实话，我刚入行那会儿，也觉得理论枯燥。但后来吃过亏，才明白它的重要性。咱们不搞学术研究，就讲怎么用这些理论指导仿真。
2.1 麦克斯韦方程组：仿真的“宪法”
麦克斯韦方程组，说白了就是电磁场的基本法则。所有电磁现象，都逃不过这四条定律。在电机仿真里，我们主要关心它的积分形式和微分形式。
我个人习惯，在脑海里把电机简化成一个“场”的问题。电流产生磁场，磁场变化又感应出电场。这背后的逻辑，全在方程组里。
核心方程（微分形式）：
- 安培环路定律： ∇ × H = J + ∂D/∂t （电流和变化的电场都能产生磁场）
- 法拉第电磁感应定律： ∇ × E = -∂B/∂t （变化的磁场产生电场，这是电机工作的根本）
- 高斯磁定律： ∇ · B = 0 （磁感线是闭合的，没有起点和终点）
- 高斯电定律： ∇ · D = ρ （电荷是电场的源）
你想想看，我们仿真永磁同步电机，算反电动势、算转矩，哪一步离得开这些？我在做一台高速电机项目时，发现仿真结果和实测差了10%。后来排查，就是忽略了位移电流∂D/∂t的影响。嗯，高频下这玩意儿不能省。
2.2 磁矢位与标量位：求解的“捷径”
直接求解麦克斯韦方程组，计算量太大。所以我们需要引入“位函数”来简化问题。这就好比解方程组，先设个未知数，再代入求解。
- 磁矢位 A： 满足 B = ∇ × A。引入它，可以自动满足高斯磁定律。在二维仿真中，A只有z方向分量，问题大大简化。我常用的二维静磁场仿真，核心就是求解A的泊松方程。
- 标量位 φ： 在无电流区域，可以引入标量磁位 φ_m，满足 H = -∇φ_m。这能进一步减少未知量，加快计算速度。
我的经验： 在Ansys Maxwell里做二维仿真，默认就是求解磁矢位A。你不需要手动设置，但理解它有助于你判断边界条件。比如，你给模型加“矢量势边界条件”，其实就是指定了A的值。
2.3 边界条件：给“场”画个圈
仿真计算必须在有限的区域内进行。边界条件，就是告诉求解器，在这个区域的边界上场量是什么状态。边界设错了，结果肯定不对。
我曾经犯过一个低级错误：做电机外转子仿真时，把最外层的边界设成了“自然边界”（Neumann边界），结果磁力线全跑出去了，漏磁大得离谱。后来改成“平行边界”（Dirichlet边界），结果才正常。
| 边界类型 | 物理含义 | 仿真应用场景 |
|---|---|---|
| 狄利克雷边界 (Dirichlet) | 指定边界上的磁矢位A值 | 模拟磁屏蔽、对称面（A=0） |
| 诺伊曼边界 (Neumann) | 指定边界上A的法向导数 | 模拟开边界、磁通垂直边界 |
| 周期边界 (Periodic) | 边界上场量呈周期变化 | 电机极槽配合、简化模型 |
| 主从边界 (Master/Slave) | 强制两个边界上的场量相等或相反 | 模拟电机的一个极或一个齿槽 |
避坑指南： 我曾经在设置周期边界时，把“主边界”和“从边界”的方向搞反了，导致磁力线不连续，转矩脉动算出来是错的。记住，主从边界的网格节点必须一一对应，否则会报错或结果异常。
2.4 材料本构关系：给“场”赋予灵魂
麦克斯韦方程组是通用的，但材料不同，电磁特性就不同。本构关系，就是描述材料在电磁场中如何响应的方程。
- 线性材料： B = μH，J = σE。μ是磁导率，σ是电导率。空气、铜、铝都属于这类。
- 非线性材料： B和H不是线性关系，比如硅钢片。它的B-H曲线是弯曲的，有饱和特性。仿真时必须输入真实的B-H曲线，否则铁耗、电感都算不准。
- 永磁材料： 有剩磁Br和内禀矫顽力Hcj。在仿真中，通常用“等效电流面”或“剩余磁通密度”来建模。
我记得有一次，客户要求优化一台电机的效率。我用了默认的B-H曲线，结果算出来的铁耗比实测低30%。后来换了供应商提供的实测曲线，结果就对上了。所以，材料数据一定要用真实的，别偷懒。
2.5 能量与力：仿真的“最终目标”
我们做仿真，最终是为了算出力、转矩、损耗这些宏观量。这些量都可以从电磁场能量中推导出来。
- 磁场能量： W = 1/2 ∫ (B·H) dV。存储在磁场中的能量。
- 虚功原理： 力 F = -∂W/∂x。通过计算能量对位移的偏导，得到电磁力。这是有限元软件计算转矩的常用方法。
- 麦克斯韦应力张量： 通过计算包围物体的闭合曲面上的应力积分，得到总力和转矩。这种方法不需要网格变形，计算速度快。
实战要点： 在Maxwell中计算转矩，有两种方法：
- 虚功法： 精度高，但需要设置“参数化”或“移动网格”。适合计算静态转矩。
- 麦克斯韦应力张量法： 速度快，适合瞬态仿真。但要注意，积分路径必须选在气隙中间，且网格要足够密，否则误差大。
我个人习惯，在瞬态场中算转矩，用麦克斯韦应力张量法。但在做齿槽转矩优化时，我会用虚功法，因为它对网格的微小变化更敏感，能捕捉到更细微的波动。
知识体系总览
下面这张图，是我自己总结的本章知识脉络。你可以把它当作一个“地图”，随时回来对照。
好了，这一章的内容就到这里。理论是枯燥的，但它是你成为高手的基础。下一章，我们会把这些理论用到实际的电机模型中去。到时候，你会发现，原来这些公式和概念，都是活的。