电磁学基础 — 麦克斯韦方程组到电机仿真完整解析
老师，我做电机设计，听说要用电磁仿真软件，但我本科学的机械，电磁场理论基本忘了。能从头讲讲吗？
没问题！电磁学的核心就是麦克斯韦方程组——4个方程统治了所有经典电磁现象，从手机天线到电动汽车电机。不同的是，电磁分析的场变量是向量（$\mathbf{E}$、$\mathbf{B}$），而不是标量（温度）或位移向量，数值处理上有独特的挑战。
1. 基本电磁量
| 符号 | 名称 | 单位 | 物理含义 |
|---|---|---|---|
| $\mathbf{E}$ | 电场强度 | V/m | 单位正电荷所受的电场力 |
| $\mathbf{D}$ | 电位移矢量 | C/m² | $\mathbf{D} = \varepsilon\mathbf{E}$，含介质极化效应 |
| $\mathbf{B}$ | 磁感应强度（磁通量密度） | T（特斯拉） | 磁场对运动电荷的力 $\mathbf{F} = q\mathbf{v}\times\mathbf{B}$ |
| $\mathbf{H}$ | 磁场强度 | A/m | $\mathbf{H} = \mathbf{B}/\mu$，不含磁化效应 |
| $\mathbf{J}$ | 电流密度 | A/m² | 单位面积的电流 |
| $\rho_e$ | 自由电荷体密度 | C/m³ | 单位体积的自由电荷 |
| $\varepsilon$ | 介电常数（电容率） | F/m | $\varepsilon = \varepsilon_r \varepsilon_0$，$\varepsilon_0 = 8.854\times10^{-12}$ F/m |
| $\mu$ | 磁导率 | H/m | $\mu = \mu_r \mu_0$，$\mu_0 = 4\pi\times10^{-7}$ H/m |
| $\sigma$ | 电导率 | S/m | $\mathbf{J} = \sigma\mathbf{E}$（欧姆定律微分形式） |
2. 麦克斯韦方程组
麦克斯韦方程组到底有几个？它们分别说的是什么？
4个方程，每个都对应一个基本物理定律。简单说：高斯电场定律（电荷产生电场）、高斯磁场定律（没有磁单极子）、法拉第电磁感应定律（变化磁场产生电场）、安培-麦克斯韦定律（电流和变化电场产生磁场）。
2.1 微分形式
麦克斯韦的天才贡献在于第四个方程里加入了位移电流 $\frac{\partial \mathbf{D}}{\partial t}$——正是这一项使方程组预言了电磁波的存在，赫兹在1888年实验验证了这一预言。
2.2 本构关系
在线性各向同性介质中：
2.3 麦克斯韦方程的"统一"意义
这4个方程各自对应什么"已知定律"？麦克斯韦的贡献是什么？
前三个方程的核心已经在麦克斯韦之前建立（高斯、法拉第等人）。麦克斯韦的关键贡献有两个：① 在安培定律里加入了"位移电流"项（$\partial\mathbf{D}/\partial t$），使方程组自洽；② 从这4个方程推导出波动方程，预言光是电磁波，并计算出波速 $c = 1/\sqrt{\varepsilon_0\mu_0}$——和当时测量的光速完全吻合！这是19世纪物理学最伟大的成就之一。
3. 静电场分析
静电场（$\partial/\partial t = 0$，无磁场）中，法拉第定律退化为 $\nabla\times\mathbf{E}=0$，意味着电场是保守场，可以引入标量电位 $\varphi$：
代入高斯电场定律，得到泊松方程：
无自由电荷时（大多数介质内部），退化为拉普拉斯方程：
3.1 静电场的应用
- 电容计算：给定几何形状的电极，计算电容 $C = Q/V$（电荷/电压）
- 绝缘体耐压分析：高压电力设备（变压器、电缆接头）中找电场最强点，判断击穿风险
- 静电放电（ESD）分析：电子产品中评估芯片静电损伤风险
- MEMS传感器：静电驱动微执行器的力-位移关系
高压电力设备的绝缘分析，具体怎么做？
典型流程：建立几何（电极、绝缘体、空气间隙）→ 设定材料介电常数 → 在电极上施加电压边界条件 → 求解泊松方程得到电位分布 → 计算电场强度 $\mathbf{E} = -\nabla\varphi$ → 找出最大电场点，与材料的击穿场强（kV/mm级别）比较。关键注意：尖端、凹陷、气泡处电场会急剧集中——这些地方最容易放电。
4. 静磁场分析
静磁场（$\partial/\partial t = 0$，无时变电场）中，$\nabla\cdot\mathbf{B}=0$ 意味着可以引入磁矢位 $\mathbf{A}$：
代入安培定律，得到控制方程：
其中 $\mathbf{J}_s$ 是源电流密度（线圈电流）。
4.1 静磁场分析的应用
- 永磁体设计：永磁电机、扬声器、MRI磁体的磁场分布
- 电感计算：$L = \Phi / I$（磁通量/电流），用于变压器、电感器设计
- 力和转矩计算：对电流或磁化材料施加的磁力，是电机转矩的基础
- 磁屏蔽设计：高导磁材料（坡莫合金）用于屏蔽外部磁场干扰
4.2 磁力计算——麦克斯韦应力张量
作用在磁性材料上的力密度可通过麦克斯韦应力张量计算：
将应力张量在闭合曲面上积分，可得作用在该曲面内所有物体上的总力。这是电机转矩计算的理论基础。
| 材料 | 相对磁导率 μᵣ | 特点 | 典型应用 |
|---|---|---|---|
| 空气/真空 | 1.0 | 基准 | — |
| 铜（导体） | ~1.0 | 顺磁性，几乎不影响磁场 | 线圈绕组 |
| 硅钢（电工钢） | 2000～10000 | 高导磁，低铁损 | 电机铁芯、变压器 |
| 坡莫合金（80Ni-20Fe） | ~80000 | 极高导磁 | 磁屏蔽、传感器 |
| 钕铁硼永磁体（NdFeB） | ~1.05（磁导率低） | 高矫顽力，强剩磁（B_r~1.2T） | 永磁电机、硬盘 |
5. 涡流分析与集肤效应
做变压器铁芯设计时，提到"铁损"，里面有"涡流损耗"，这是什么？怎么仿真？
涡流是变化磁场在导体中感应出的循环电流（法拉第定律）。这些电流会在导体中产生焦耳热，就是涡流损耗。变压器铁芯叠片（0.3～0.5 mm薄片）就是为了切断涡流路径，减少损耗。频率越高、导体越厚，涡流越严重。
5.1 涡流方程
准静态（忽略位移电流，$\omega \ll \sigma/\varepsilon$）的时谐电磁场（角频率 $\omega$）满足：
第二项 $j\omega\sigma\mathbf{A}$ 就是涡流项——包含了导体中感应电流对磁场的反作用。
5.2 集肤效应与集肤深度
交变电流在导体中不均匀分布，倾向于集中在表面附近，这就是集肤效应。电流密度从表面向内按指数衰减：
其中集肤深度 $\delta$（电流密度降至表面值的 $1/e \approx 37\%$ 的深度）：
$\rho = 1/\sigma$ 是电阻率，$f$ 是频率。频率越高，集肤深度越浅，有效导电截面积越小，电阻越大。
| 材料 | 电阻率 ρ (Ω·m) | 集肤深度 δ（50Hz） | 集肤深度 δ（1MHz） |
|---|---|---|---|
| 铜 | 1.68×10⁻⁸ | 9.3 mm | 66 μm |
| 铝 | 2.82×10⁻⁸ | 12 mm | 85 μm |
| 钢（低碳钢） | ~1×10⁻⁷，μᵣ~100 | ~0.7 mm | ~50 μm |
哦，所以高频电源的线圈用细导线编织的利兹线，就是为了解决集肤效应问题？
完全正确！利兹线（Litz wire）由许多根绝缘的细导线绞合而成，每根细线的直径都远小于集肤深度，使所有细线都能有效导电。比如100kHz的感应加热线圈，单根铜线的集肤深度约0.2 mm，你得用直径 <0.1 mm 的细线才有意义。
5.3 近邻效应
集肤效应的"升级版"——当多导体靠近时，相邻导体的磁场会影响电流分布，使有效截面积进一步减小。多层线圈、高频变压器绕组设计时必须考虑。
6. 电磁波与高频分析
天线和射频（RF）电路的仿真，用的还是麦克斯韦方程吗？
对，高频电磁分析（RF/微波、天线）也是麦克斯韦方程，但不能忽略位移电流项了——正是位移电流使电磁波能在真空中传播。频率高到一定程度（比如GHz级别），导线的分布电感/电容不能再忽略，整个电路都得用场的概念来分析，而不是集总电路模型。
6.1 波动方程
在无源、无损耗介质中，从麦克斯韦方程组导出电场的波动方程：
电磁波传播速度：
真空中 $c_0 = 3\times10^8$ m/s。介质中的波速更慢，折射率 $n = \sqrt{\mu_r\varepsilon_r}$。
6.2 高频分析方法
| 频率范围 | 典型应用 | 仿真方法 | 软件 |
|---|---|---|---|
| DC～10 kHz | 电机、变压器、感应加热 | 准静态FEM（涡流） | ANSYS Maxwell, Flux |
| 10 kHz ～ 100 MHz | 开关电源、EMC分析 | 时域FEM/FDTD | ANSYS HFSS, CST |
| 100 MHz ～ 30 GHz | 天线、微波电路、5G | 频域FEM/矩量法（MoM） | ANSYS HFSS, FEKO |
| 光频（THz～PHz） | 光波导、光子器件 | FDTD, RCWA | Lumerical FDTD, MEEP |
6.3 S参数
高频电路用散射参数（S参数）描述端口的反射和传输特性：
$a_i$、$b_i$ 是端口入射/反射波幅。$|S_{11}|^2$ 是反射功率比，$|S_{21}|^2$ 是传输功率比。天线的阻抗匹配要求 $|S_{11}| < -10$ dB（反射功率 <10%）。
7. FEM电磁分析
电磁场的有限元和结构分析的有限元，在数值处理上有什么不同？
最大的区别是场变量类型。结构FEM的未知量是标量（温度）或位移向量（在节点上定义）；电磁FEM的核心未知量是矢量场——用棱边元（Edge Element，Nédélec元）在单元边上定义切向分量，这样才能正确处理界面上的切向连续条件，避免"伪解"（虚假的非物理模式）。
7.1 棱边有限元（Nédélec元）
电磁矢量场的正确离散方式：
- 节点元（H1 元）：在节点上定义，适合标量场（电位 $\varphi$），但用于矢量场会产生伪解
- 棱边元（H(curl) 元）：在单元边上定义切向分量，自动满足 $\nabla\times\mathbf{A}$ 的正则性，是电磁场FEM的标准选择
- 面元（H(div) 元）：在面上定义法向分量，适合 $\mathbf{B}$、$\mathbf{D}$ 等场量
7.2 非线性磁性材料——B-H曲线
铁磁材料（硅钢、铸铁）的磁导率不是常数，而是随磁场强度变化：$\mu = B/H$，通过B-H曲线（磁化曲线）描述。仿真时需要输入B-H曲线数据，求解器在每个迭代步用Newton-Raphson迭代处理非线性：
磁饱和时 $\mu$ 急剧下降，电机铁芯的工作点通常设在B-H曲线的"膝部"（饱和开始前），以获得最大磁通量密度。
7.3 时谐分析与瞬态分析
- 时谐（AC稳态）分析：假设所有场量以 $e^{j\omega t}$ 变化，将时间导数换成 $j\omega$，得到复数代数方程。适合稳态电机、变压器、感应加热分析
- 瞬态分析：需要直接时间积分，处理开关瞬态、脉冲电流、启动过程。计算量比时谐分析大得多
8. 电机与执行器CAE
电动汽车的永磁同步电机，用CAE怎么分析？需要考虑哪些物理场？
永磁同步电机（PMSM）是典型的多物理场耦合问题。核心是电磁分析（转矩、反电动势、铁损计算），但还要考虑热分析（绕组铜损和铁损引起的温升，影响永磁体退磁）、结构分析（转子高速旋转的离心应力、振动噪声），以及转子旋转引起的运动边界处理。
8.1 永磁电机关键仿真内容
- 磁通量分布：气隙磁通量密度波形（决定转矩纹波大小）
- 铁损计算：Steinmetz方程 $P_{core} = k_h f B^n + k_e f^2 B^2$（磁滞损耗+涡流损耗）
- 转矩计算：通过虚功原理或麦克斯韦应力张量积分
- 反电动势（Back-EMF）：评估波形谐波含量，影响控制性能
- 齿槽转矩：永磁体与定子齿相互作用产生的脉动转矩，影响低速平顺性
8.2 旋转机械的滑动网格
电机转子旋转时，气隙两侧（定子侧和转子侧）相对运动。处理方法：
- 滑动网格（Sliding Mesh）：气隙处设置滑动界面，转子网格随时间旋转，每步重新计算界面插值。最准确但最耗时。
- 运动带（Motion Band）：气隙区域用特殊的运动区域处理，减少网格重划分。
- 冻结渗透法：转子在若干固定位置求解，然后叠加，适合线性分析快速评估。
8.3 电磁-热-结构耦合流程
- 电磁分析 → 计算铜损（$P_{Cu} = \int J^2/\sigma\, dV$）和铁损（Steinmetz方程）
- 热分析 → 以损耗为热源，计算温度场（绕组、铁芯、永磁体温度）
- 更新材料参数 → 铜的电阻率随温度变化（$\sigma(T) = \sigma_0/(1+\alpha\Delta T)$）；永磁体剩磁随温度降低
- 返回步骤1重新迭代，直至温度和损耗收敛
- 结构分析 → 以温度场为载荷，评估热应力和振动噪声（NVH）
| 分析类型 | 主要求解量 | 软件选择 | 网格要求 |
|---|---|---|---|
| 静磁场（永磁体） | 磁通量密度B、转矩 | ANSYS Maxwell, Flux, JMAG | 气隙需足够细（至少3层网格） |
| 涡流（铁损） | 涡流损耗、集肤效应 | 同上（时谐模式） | 集肤深度内至少2～3层网格 |
| 瞬态（启动/换相） | 时变电流、转矩波形 | Maxwell瞬态，Flux Transient | 滑动网格，运算时间长 |
| 振动噪声（NVH） | 谐波力、声辐射 | Maxwell + ANSYS Mechanical + Fluent | 多物理场耦合 |