电磁分析 — 有限元仿真指南
教授，电磁分析和结构分析有什么本质不同？我听说电机设计必须用电磁仿真，但我不太清楚它到底在算什么。
简单说，结构分析算的是力和变形，而电磁分析算的是空间里的电场和磁场分布。比如你设计一个电机，你需要知道转子磁铁产生的磁通量怎么分布在铁芯里、哪里会磁饱和、扭矩是多少——这些都要靠电磁仿真来预测，而不是等到做出实物再测试。
那控制这些现象的基础方程是什么？我记得大学里学过麦克斯韦方程组，但那时候感觉很抽象……
麦克斯韦方程组就是电磁学的"根基"，有限元电磁求解器本质上就是在数值求解这四个方程。你不需要背下来，但要理解每个方程的物理含义：安培定律说电流和变化电场会产生磁场；法拉第定律说变化磁场会产生电场；第三个说没有磁单极子——磁力线总是闭合的；第四个说自由电荷是电场的来源。
那实际工程里是哪类问题最常见？每次都要把四个方程全部算吗？
不，大多数工程问题都是简化版的。低频的电机、变压器用"静磁分析"或"涡流分析"就够了；只有在微波天线、雷达这类高频应用里才需要全波仿真。选对分析类型，计算量可以差几个数量级——这是电磁仿真工程师最重要的判断力之一。
1. 麦克斯韦方程组 — 电磁分析的物理基础
所有电磁现象均由麦克斯韦方程组（Maxwell's Equations）描述。微分形式如下：
安培定律（含位移电流）
$$\nabla \times \mathbf{H} = \mathbf{J} + \frac{\partial \mathbf{D}}{\partial t}$$其中 $\mathbf{H}$ 为磁场强度，$\mathbf{J}$ 为自由电流密度，$\mathbf{D}$ 为电位移矢量。该方程描述电流（含位移电流）产生磁场的机制，是电感、电机和变压器分析的出发点。
法拉第电磁感应定律
$$\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$$变化的磁通量 $\mathbf{B}$ 会感应出电场 $\mathbf{E}$。这是感应加热（涡流加热）和无线充电背后的基本原理。
磁场散度定律（无磁单极子）
$$\nabla \cdot \mathbf{B} = 0$$磁力线永远是闭合曲线，不存在磁单极子。该方程在有限元求解中对规范化（Gauging）处理有重要意义。
高斯定律
$$\nabla \cdot \mathbf{D} = \rho_f$$其中 $\rho_f$ 为自由电荷密度。自由电荷是电位移场的散度源，是静电分析的核心方程。
本构关系
在线性各向同性介质中：
$$\mathbf{B} = \mu \mathbf{H}, \quad \mathbf{D} = \varepsilon \mathbf{E}, \quad \mathbf{J} = \sigma \mathbf{E}$$其中 $\mu$ 为磁导率，$\varepsilon$ 为介电常数，$\sigma$ 为电导率。在铁磁材料中，$\mu$ 是磁场强度的非线性函数（B-H曲线），这是电机仿真必须处理非线性的根本原因。
2. 主要电磁分析类型
分析类型这么多，我怎么知道我的项目该用哪一种？
问自己两个问题：一，有没有时变量？如果是静止的电场或磁场，就用"静态分析"。二，频率高不高？工频（50/60 Hz）电机用涡流分析；GHz级的天线就要全波仿真了。这两个判断基本能把你导向正确的分析类型。
2.1 静磁场分析（Magnetostatics）
适用于：直流电流激励的电磁铁、永磁体磁路、直流电机空载磁场。忽略时间导数，方程简化为：
$$\nabla \times \mathbf{H} = \mathbf{J}_s, \quad \nabla \cdot \mathbf{B} = 0$$引入磁矢量位（Magnetic Vector Potential）$\mathbf{A}$，令 $\mathbf{B} = \nabla \times \mathbf{A}$，转化为关于 $\mathbf{A}$ 的偏微分方程。典型输出：磁通量密度分布、电感、磁饱和区域、扭矩（麦克斯韦张量积分）。
2.2 静电场分析（Electrostatics）
适用于：电容器、高压绝缘件（GIS开关、变压器绝缘设计）、MEMS传感器。忽略磁效应，用电标量位 $\varphi$ 求解：
$$\nabla \cdot (\varepsilon \nabla \varphi) = -\rho_f$$典型输出：电场强度分布、电容值、击穿风险区域（$E$ 场超过材料绝缘强度的位置）。
2.3 涡流分析（Eddy Current / AC Magnetics）
适用于：交流电机（50/60 Hz）、变压器铁芯损耗、感应加热系统、涡流检测（NDT）。时变磁场在导体中感应出涡流 $\mathbf{J}_{eddy}$，控制方程为：
$$\nabla \times \left(\frac{1}{\mu}\nabla \times \mathbf{A}\right) + \sigma \frac{\partial \mathbf{A}}{\partial t} = \mathbf{J}_s$$典型输出：铁损（磁滞损耗 + 涡流损耗）、集肤效应深度 $\delta = \sqrt{2/(\omega \mu \sigma)}$、温升分布（配合热分析）。
2.4 全波仿真（Full-Wave / High Frequency）
适用于：天线辐射、微波滤波器、雷达截面（RCS）、PCB上的信号完整性（SI）。频率足够高时，位移电流项 $\partial \mathbf{D}/\partial t$ 不可忽略，必须求解完整麦克斯韦方程组。常用方法：
- 时域有限差分（FDTD）：对宽频带问题高效，是CST Microwave Studio的核心算法
- 有限元频域（FEM-FD）：对复杂几何精度高，是ANSYS HFSS的核心算法
- 矩量法（MoM）：适合天线辐射等开域问题，如FEKO
3. 主要工程应用
电磁仿真在工业上的应用有哪些？感觉比结构分析更小众……
其实随着电动汽车和新能源的爆发，电磁仿真的需求已经大幅增长了。EV的驱动电机效率图、逆变器的EMC合规测试、车载充电器的感应线圈设计——这些现在都是汽车厂商的刚需。更别说通信、雷达、医疗MRI这些传统大头了。
| 应用领域 | 典型问题 | 分析类型 |
|---|---|---|
| 电机设计 | PMSM/SRM扭矩-转速特性、效率图（Loss Map）、去磁风险 | 静磁 + 涡流瞬态 |
| 变压器 | 铁芯磁通密度分布、绕组漏感、铁损/铜损计算、温升预测 | 涡流（交流） |
| EMC / EMI | 设备辐射发射量测量、印刷电路板（PCB）传导干扰抑制 | 全波（FDTD / FEM-FD） |
| 天线设计 | 5G毫米波阵列增益/方向图、VSWR匹配优化 | 全波（FEM-FD / MoM） |
| 感应加热 | 工件内涡流分布优化、加热均匀性、线圈几何设计 | 涡流 + 热耦合 |
| 高压绝缘 | GIS开关局部放电风险评估、绝缘子电场集中区识别 | 静电场 |
| 无线充电（WPT） | 线圈互感、传输效率、漏磁对周边组件的影响 | 涡流（准静态） |
4. 有限元求解：节点元 vs 棱边元
我看到有些电磁软件里有"节点元"和"棱边元"两种网格选项，有什么区别？用错了会怎样？
这是电磁FEM里非常关键的一点。节点元把未知量（磁矢量位 $\mathbf{A}$）的每个分量定义在节点上，简单但有一个严重问题：它在材料界面上会产生"虚假解"（Spurious modes），即数学上成立但物理上无意义的伪模式。棱边元把未知量定义在单元的边（Edge）上，天然地满足界面切向连续条件，完全消除了伪模式问题。所以现代电磁仿真软件，尤其是做全波仿真的，都默认用棱边元。
4.1 节点元（Nodal Elements）
将未知量（如标量位 $\varphi$ 或矢量位分量）定义在网格节点上。优点是实现简单，与结构FEM体系一致。缺点：对矢量电磁场问题，在材料界面处强制施加了切向连续条件，会引入虚假解（Spurious Modes）。主要用于静电场（标量位）和低频磁标量位问题。
4.2 棱边元（Edge Elements / Nédélec Elements）
将矢量场的切向分量定义在单元的"棱"（Edge）上。自动满足：
- 切向分量在材料界面连续 $\mathbf{n} \times (\mathbf{H}_1 - \mathbf{H}_2) = \mathbf{K}$
- 法向分量允许在界面不连续（对应 $\mathbf{B}$ 的法向连续但 $\mathbf{H}$ 不连续）
完全消除虚假模式，是现代全波仿真（ANSYS HFSS、COMSOL、Elmer FEM）和涡流分析的标准选择。
4.3 规范化（Gauging）
磁矢量位 $\mathbf{A}$ 有自由度（加任意梯度场 $\nabla \psi$ 不改变 $\mathbf{B}$），需施加规范条件来确保方程组有唯一解：
$$\text{库仑规范：} \quad \nabla \cdot \mathbf{A} = 0$$ $$\text{洛伦兹规范：} \quad \nabla \cdot \mathbf{A} + \mu\varepsilon \frac{\partial \varphi}{\partial t} = 0$$多数低频求解器采用树-余树（Tree-Cotree）技术自动处理规范化，用户通常无需手动设置。
5. 多物理场耦合
电磁分析和热分析可以耦合吗？比如电机发热这种问题？
当然可以，而且实际工程里这种耦合非常常见。以感应加热为例：涡流在工件中产生焦耳热 $Q = \mathbf{J} \cdot \mathbf{E}$，热量导致温度升高；而温度升高会改变材料的电导率 $\sigma(T)$，进而影响涡流分布——这就构成了双向耦合。另一个典型是电机里的电磁-结构耦合：定子绕组受到洛伦兹力 $\mathbf{F} = \mathbf{J} \times \mathbf{B}$，这个力会引起定子铁芯的振动，产生电磁噪声（NVH问题）。
5.1 电磁-热耦合（焦耳加热 / 感应加热）
导体中的涡流产生焦耳热，作为热源项输入热分析：
$$Q_{Joule} = \frac{|\mathbf{J}|^2}{\sigma} = \sigma |\mathbf{E}|^2 \quad \text{[W/m}^3\text{]}$$温度反过来影响材料导电率 $\sigma(T)$ 和磁导率 $\mu(T)$，形成双向耦合。应用：感应淬火、感应烹饪、无线充电效率分析、电机绕组寿命预测。
5.2 电磁-结构耦合（洛伦兹力）
载流导体在磁场中受到体积力（洛伦兹力），作为载荷传入结构分析：
$$\mathbf{f}_{vol} = \mathbf{J} \times \mathbf{B} \quad \text{[N/m}^3\text{]}$$对于铁磁体，还需考虑麦克斯韦应力张量（Maxwell Stress Tensor）。应用：大型变压器绕组的机械强度校核、电磁成形（电磁力驱动金属变形）、电机定子铁芯噪声（NVH）分析。
5.3 压电耦合（Piezoelectric）
应变与电场之间的直接/逆压电效应，用于传感器（超声探头）和执行器（压电扬声器）设计。COMSOL和ANSYS均支持结构-电场联立求解。
6. 主要软件对比
电磁仿真软件好像比结构分析的种类更多，各家有什么侧重？我要做电机设计，该用哪款？
做电机的话，日系厂商首选JMAG——它是专为旋转电机设计的，电机效率图（Loss Map）和多工况扫描功能非常完善，在日本汽车行业用得极多。如果你已经在用ANSYS生态，Maxwell和Mechanical可以无缝耦合，做电磁-热-结构三场耦合很方便。COMSOL的优势是多物理场灵活性，但电机专用功能不如JMAG深。免费开源的话，Elmer FEM也能做到不错的涡流和静磁分析，适合学习用。
| 软件 | 开发商 | 核心优势 | 典型用途 |
|---|---|---|---|
| ANSYS Maxwell | ANSYS Inc. | 低频电磁（电机/变压器）专业求解器；与ANSYS Mechanical/Fluent多物理场耦合无缝；自动自适应网格（AML） | 电机、变压器、EMC、感应加热 |
| ANSYS HFSS | ANSYS Inc. | 业界标杆高频全波求解器（FEM-FD）；天线、微波器件精度极高；SBR+射线追踪法处理大型场景 | 天线、雷达、PCB信号完整性 |
| COMSOL Multiphysics | COMSOL AB | AC/DC + RF + 结构 + 热等多场同平台耦合；用户界面友好；支持自定义PDE方程；适合科研 | 多物理场研究、压电、生物电磁 |
| JMAG | JSOL Corporation（日本） | 旋转/线性电机专用；效率图（Loss Map）自动生成；铁损模型精度高；NVH振动噪声分析集成 | EV/HEV驱动电机、发电机 |
| CST Studio Suite | Dassault Systèmes | 时域FDTD（T求解器）宽频带高效；低频FEM+高频FDTD统一平台；EMC/EMI合规测试首选 | EMC、天线系统、SI/PI |
| Elmer FEM | CSC Finland（开源） | 完全开源；支持静磁、涡流、全波；可与OpenFOAM耦合做电磁-流体问题；适合学术研究 | 学术研究、MHD（磁流体） |