上一篇我们把静态场下的麦克斯韦方程组拆成了两组静电场方程组和静磁场方程组。很多朋友看完后问我为什么明明是四个方程非要分开讲其实道理很简单静态场本质上就是时间不变化电场是电场、磁场是磁场两者勉强靠变化率耦合着一旦加上“静态”两个字耦合的桥梁就直接断了。电场守着电荷、磁场守着电流各干各的互不打扰。这也是为什么工程上能把电磁场问题拆成独立的静电场题目和静磁场题目来处理。这篇续文我打算把重点放在磁场这一半上磁通连续性方程、安培环路方程、矢量磁位 ( \mathbf{A} ) 的引入以及那些直接影响工程仿真和电感计算的边界条件与求解套路。适合正在学电磁场理论的电类专业学生、做电机变压器或磁元件设计的工程师以及用有限元软件算磁场却又总觉得“差口气”的仿真党。1. 静态场方程组到底删掉了哪些东西1.1 完整麦克斯韦方程组里的“动态项”先把完整方程组摆出来看一眼[ abla \cdot \mathbf{D} \rho, \quad abla \cdot \mathbf{B} 0, \quad abla \times \mathbf{E} -\frac{\partial \mathbf{B}}{\partial t}, \quad abla \times \mathbf{H} \mathbf{J} \frac{\partial \mathbf{D}}{\partial t} ]我当年学的时候最要命的就是分不清到底哪个方程是“独立的”哪个是“从属的”。这里有一个非常实用的理解方式散度方程是给场的“源”定规矩的旋度方程是给场的“旋转”定规矩的。电磁波能跑起来靠的是两个旋度方程里的时间导数项互相喂能量一个是 ( \partial \mathbf{B}/\partial t ) 生电场一个是 ( \partial \mathbf{D}/\partial t ) 生磁场。但在静态场下所有物理量对时间的偏导数全部等于零。于是方程组直接断成两截静电场[ abla \cdot \mathbf{D} \rho,\quad abla \times \mathbf{E} 0 ]静磁场[ abla \cdot \mathbf{B} 0,\quad abla \times \mathbf{H} \mathbf{J} ]注意看两个原本“连接”电场和磁场的旋度方程里时间项全部消失了。第一篇文章我们处理的静电场只剩两个方程高斯定律和电场无旋。这一篇要对付的静磁场同样只剩两个磁通连续性定律和安培环路定律。1.2 一个常见的理解误区静态场不等于没有电流很多人一看到“静态”两个字第一反应是都没变化了哪来的电流这是个非常普遍的误解。静态场说的是场量不随时间变化而不是说电荷不能运动。举个最直观的例子一根通着恒定直流电的导线周围产生的磁场是稳定不变的这当然属于静态场的范畴。电流虽然是电荷在流动但只要电流的大小和方向不随时间变它产生的磁场就是静态磁场。所以静磁场方程右端那个 ( \mathbf{J} ) 必须保留它是磁场的源也是静磁场存在的基础。意识到这一点特别重要因为我在给读者答疑时至少有五六个人问过“恒定电流算不算静态场”。答案是肯定的。只要你把电流源一接通导线周围的磁场会从零跳变到稳定值然后一直保持下去。这个稳定后的状态就是静态场问题。电场和磁场没有依赖关系了——恒定电流产生的磁场不会反过来再生成一个电场因为 ( \partial \mathbf{B}/\partial t0 )。1.3 本文要解决的核心问题静磁场的两个方程看起来很简洁但实际使用中坑特别多。比如为什么 ( abla \cdot \mathbf{B}0 ) 就一定意味着磁力线是闭合的它跟“磁单极子不存在”是什么关系( abla \times \mathbf{H}\mathbf{J} ) 为什么能直接推出安培环路定律的积分形式求解静电场时我们有电势 ( \varphi )求解静磁场时对应的那个“势”到底是什么为什么叫矢量磁位铁芯和气隙交界面上磁场到底是怎么“分家”的工程上算磁路时用的那几个边界公式是怎么来的这些问题的答案一半在数学推导里另一半在物理直觉里。我尽量把两边都讲透。2. 磁通连续性方程磁力线为何永远首尾相连2.1 从 (abla \cdot \mathbf{B}0 ) 到“无磁荷”世界观静磁场第一个方程是 ( abla \cdot \mathbf{B}0 )。高斯定理告诉我们一个矢量场的散度代表的是单位体积内“源”的净流出量。散度为零意味着这个场在任何微小体积内没有净流出也就是说磁感应强度线既没有起点、也没有终点它自己必须闭环。对比一下静电场( abla \cdot \mathbf{D}\rho )电荷就是电场的“源”和“汇”。正电荷发出电场线负电荷接收电场线。而磁场没有对应的“磁荷”所以磁场线只能一圈一圈地绕。这个结论放在任何尺度下都成立小到原子核周围的磁场大到地球磁层再到银河系的磁场结构磁力线全都闭合。这个方程是“磁单极子不存在”这一物理事实的数学表达。虽然某些大统一理论预言了磁单极子的存在但至今没人真正观测到。所以工程计算中默认 ( abla \cdot \mathbf{B}0 ) 是铁律。2.2 积分形式与磁通守恒把 ( abla \cdot \mathbf{B}0 ) 在任意体积 ( V ) 上积分再用高斯散度定理转化为面积分就得到[ \oint_S \mathbf{B} \cdot d\mathbf{S} 0 ]它的含义是穿出任意闭合曲面的磁通量恒等于零。进一步可以推出同一个磁通管想象一根由磁力线围成的管道任意截面上的磁通量都相等。这就是变压器、电机里“主磁通守恒”的底层逻辑铁芯里的磁通不会凭空消失它要么在铁芯内部闭合要么通过气隙穿出去形成漏磁。让我给你一个特别实际的工程经验算变压器时如果你发现铁芯两个截面上算出来的磁通量不相等那不用怀疑一定是网格剖分出了问题或者材料属性设置有误。这不是物理规律的问题是模型的问题。因为磁通连续性方程决定了一个铁芯柱上各截面磁通必须相等。2.3 为啥还要单列一个方程有人可能会问既然磁场是电流产生的那我把电流产生的磁场算出来不就行了为什么还要专门强调 ( abla \cdot \mathbf{B}0 )因为求解磁场时这个散度约束是保证解“合法”的硬性条件。后面要讲的矢量磁位就是从它出发构造出来的。换句话说它看起来只是个小小的散度方程实际上直接决定了静磁场求解的整个数学框架。3. 安培环路方程理解电流如何制造旋涡3.1 旋度的直觉让场“转”起来的源静磁场第二个方程是 ( abla \times \mathbf{H}\mathbf{J} )。旋度的物理直觉是场在一个点附近的旋转趋势。电流密度 ( \mathbf{J} ) 就是产生这种旋转趋势的“原动力”。怎么理解呢想象一个平静的湖面你把手伸进水里搅动湖面上就会出现漩涡。电流就好比这只手磁场就是水面形成的旋涡流场。电流越大磁场越强“旋转”越剧烈。注意这里磁场不是真的在“旋转流动”而是一种数学上的旋转趋势——你用右手握住导线四指的方向就是磁场方向。这个方程是“电生磁”的微分表述积分形式就是初中就接触过的安培环路定理[ \oint_L \mathbf{H} \cdot d\mathbf{l} \iint_S \mathbf{J} \cdot d\mathbf{S} I_{enc} ]它说的是磁场强度沿任意闭合回路的线积分等于穿过该回路所围曲面的总电流。这个形式在工程上极好用尤其是碰到对称结构。3.2 为什么静态场下没有位移电流完整的安培环路方程其实是[ abla \times \mathbf{H} \mathbf{J} \frac{\partial \mathbf{D}}{\partial t} ]多出来的 ( \partial \mathbf{D}/\partial t ) 被称为位移电流密度。它在电磁波传播、电容器充放电暂态分析中扮演核心角色。但在静态场里电场不变化位移电流直接归零。这是静态场和时变场最本质的区别之一。我实测过一个对比案例用有限元仿真一个平行板电容器当给它加恒定电压静态场磁场完全由外部导线中的电流决定极板之间的位移电流项为0磁场强度极小。而如果换成高频交流电压极板间的交变电场会产生明显的位移电流磁场分布跟静态时完全不同。这就是为什么很多高频器件仿真不能偷懒用静态求解器。3.3 对称性红利用积分形式手算磁场工程手算时积分形式的安培环路定律往往比微分形式好用得多。遇到无限长直导线、无限长螺线管、环形铁芯这类对称结构时可以直接选合适的安培环路把 ( \mathbf{H} ) 从积分号里提出来。举一个经典例子环形铁芯线圈。设铁芯平均磁路长度为 ( l )线圈匝数为 ( N )电流为 ( I )则取铁芯中心线为积分回路对称性保证 ( H ) 沿回路大小处处相等于是[ H \cdot l NI ]这个公式就是磁路欧姆定律里磁动势 ( NI ) 的直接来源。电机学里算励磁磁动势、变压器里算空载电流全都要靠它。别小看这个公式它是整个磁路计算的基石。4. 矢量磁位A静磁场求解的第一把钥匙4.1 引入矢量磁位的数学动机前面说过静磁场必须满足 ( abla \cdot \mathbf{B}0 )。而在矢量微积分里有一个非常重要的恒等式任意矢量场的旋度的散度必为零。反过来就是说如果某个矢量场的散度为零那它一定能写成另一个矢量场的旋度。[ abla \cdot (abla \times \mathbf{A}) 0 ]于是我们可以定义[ \mathbf{B} abla \times \mathbf{A} ]这里的 ( \mathbf{A} ) 就叫矢量磁位。为什么非要引入它因为直接解 ( \mathbf{B} ) 的约束太严格了既要满足旋度方程又要保证散度为零而且 ( \mathbf{B} ) 本身有三个分量直接解非常头疼。引入 ( \mathbf{A} ) 之后自动满足磁通连续性方程我们就只需要处理安培环路方程了。做一个类比静电场里我们引入电势 ( \varphi )让 ( \mathbf{E} -abla\varphi )自动满足无旋条件静磁场里矢量磁位干的是一模一样的活——通过数学构造把一个耦合方程组变成一个更简单的独立方程。不过这里有个重要的区别电势是标量好处理矢量磁位是矢量有三个分量仍然不简单。但好在很多实际场景比如二维平面场可以大大简化。4.2 库仑规范与拉普拉斯/泊松方程现在把 ( \mathbf{B} abla \times \mathbf{A} ) 代入 ( abla \times \mathbf{H} \mathbf{J} )假设介质线性且各向同性( \mathbf{B} \mu \mathbf{H} )于是[ abla \times ( abla \times \mathbf{A}) \mu \mathbf{J} ]利用矢量恒等式 ( abla \times ( abla \times \mathbf{A}) abla(abla \cdot \mathbf{A}) - abla^2 \mathbf{A} )得到[ abla(abla \cdot \mathbf{A}) - abla^2 \mathbf{A} \mu \mathbf{J} ]这里 ( \mathbf{A} ) 并不唯一——给 ( \mathbf{A} ) 加上任意一个标量场的梯度 ( abla \chi )它的旋度不变( \mathbf{B} ) 也不变。这被称为规范自由度。为了把方程化简我们可以自由选择一个规范条件。工程计算里最常用的就是库仑规范[ abla \cdot \mathbf{A} 0 ]于是方程变成[ abla^2 \mathbf{A} -\mu \mathbf{J} ]这就是矢量形式的泊松方程。在没有电流的区域它退化成矢量拉普拉斯方程[ abla^2 \mathbf{A} 0 ]到这里静磁问题就转化成了数学上相对成熟的偏微分方程求解问题。我读研那会儿第一次看到这个推导时恍然大悟原来所谓“求解磁场”最核心的步骤是求解一个矢量泊松方程。4.3 从矢量磁位到电感一个算例矢量磁位不只是中间变量它跟磁通有直接关系。根据斯托克斯定理[ \Phi \iint_S \mathbf{B} \cdot d\mathbf{S} \iint_S ( abla \times \mathbf{A}) \cdot d\mathbf{S} \oint_L \mathbf{A} \cdot d\mathbf{l} ]也就是说穿过某个曲面的磁通等于矢量磁位沿该曲面边界的环量。这个关系在电感计算里非常好用因为它让磁通的计算从面积分变成了线积分有时候反而更简单。举个例子——同轴电缆的电感计算。设内导体半径为 ( a )外导体内半径为 ( b )电流 ( I ) 均匀分布。根据对称性在内导体内部和内外导体之间的区域矢量磁位只与半径 ( r ) 有关。通过解圆柱坐标系下的泊松方程和拉普拉斯方程配合边界条件可以求出各区域的 ( \mathbf{A} )再沿回路积分得到磁通最后用 ( L \Phi/I ) 得到单位长度电感[ L \frac{\mu}{2\pi} \ln\frac{b}{a} \frac{\mu_0}{8\pi} ]前一项来自内外导体之间的磁场后一项来自内导体内部的磁场。如果你自己去翻教材发现公式长得不太一样大概率是差在是否含内导体的内自感上。我做电力电子高频变压器损耗分析时就踩过这个坑拿公式算电感发现算出来的和高频实测差很多最后才发现是没算内导体的内自感贡献而且高频下还要考虑趋肤效应导致的内自感变化。很多所谓电磁场公式“没用”其实是没搞清楚适用条件。4.4 数值仿真中库仑规范怎么处理如果你用有限元软件算静磁场会发现软件后处理里可以输出矢量磁位 ( \mathbf{A} ) 的分布。但大多数软件并不会直接求解原始方程而是会强制加上库仑规范来保证解的唯一性。这里有个实操要点:在二维静磁分析中比如电机横截面磁场分析,矢量磁位只保留轴向分量 ( A_z )方程直接退化成标量泊松方程。这正是很多商用电磁仿真软件二维静磁求解器的理论基础它名义上叫“求解磁场”实际上内部解的就是 ( A_z )。建模边界条件时通常还会设置某个边界上的 ( A_z ) 为常数或零。比如电机仿真里常见的“平行边界条件”本质就是让磁力线沿边界走。理解这一点你就不至于在设置边界条件时两眼一抹黑——你其实是在给磁位定参考值。5. 分界面的规矩边界条件的工程面孔5.1 法向分量连续为什么磁力线穿过铁芯分界面不断裂把 ( abla \cdot \mathbf{B}0 ) 应用到一个薄的跨界面扁盒子上让盒子厚度趋于零可以得到[ B_{1n} B_{2n} ]含义是磁感应强度在分界面两侧的法向分量相等。换句话说磁力线斜着穿到另一种材料里时它在法向上的“密度”不会突变。磁力线不会在分界面上凭空多出来或消失。这个边界条件的工程意义非常大。在变压器铁芯和气隙交界面上磁通从铁芯进入气隙时法向分量连续意味着“气隙磁通等于铁芯磁通”。只要铁芯截面积和气隙截面积相同气隙磁密就等于铁芯磁密。很多刚做电机设计的人会纠结气隙磁密和铁芯磁密是不是要乘个系数之类的其实在截面积相同时它们就是完全相等的。真正要注意的是切向分量气隙里的磁场会比铁芯里“歪”很多导致磁力线弯曲产生边缘效应。5.2 切向分量跳变面电流密度说了算对 ( abla \times \mathbf{H} \mathbf{J} ) 应用面积分并取极限可以得到切向边界条件[ \mathbf{n} \times (\mathbf{H}_1 - \mathbf{H}_2) \mathbf{J}_s ]( \mathbf{J}_s ) 是分界面上的自由面电流密度。如果界面上没有面电流那么磁场强度的切向分量连续[ H_{1t} H_{2t} ]实际工程里哪里会碰到面电流电机定子铁芯和绕组所在的槽区域如果电流被建模为铁芯表面的电流片这时候面电流密度就出现在边界条件里。还有电磁屏蔽分析中理想导体表面会感应出电流导致磁场被“挤”出去。这些场景下如果忽略面电流项算出来的磁场分布会明显不对。5.3 气隙里的“大头账”磁路设计必看边界条件直接决定了磁路里“磁动势降在哪”。考虑一个铁芯-气隙-铁芯结构铁芯磁导率 ( \mu_{Fe} )、气隙磁导率 ( \mu_0 )。根据切向 ( H ) 连续条件气隙中的 ( H ) 等于铁芯表面附近的 ( H )忽略边缘效应。但由于 ( B \mu H )气隙里磁感应强度 ( B_0 \mu_0 H )铁芯里 ( B_{Fe} \mu_{Fe} H )。因为 ( \mu_{Fe} \gg \mu_0 )同样的 ( H ) 下铁芯里磁通密度远大于气隙。但在开气隙的磁路中真正决定总磁动势大头的是气隙。因为磁路欧姆定律[ NI H_{Fe}l_{Fe} H_gl_g ]气隙虽短但 ( \mu_0 ) 极小所以 ( H_gl_g ) 往往占了大头。这就是为什么带气隙的电感器磁动势大部分消耗在气隙上铁芯的磁导率再高也没用气隙一开整个磁路的主要瓶颈就在那儿。我调过好几个反激变压器的气隙对这一点体会特别深:气隙开大一点电感量掉得很明显就是因为气隙这个“磁阻大头”在起作用。5.4 铁磁材料的非线性μ 不是常数前面所有推导都假设材料线性、各向同性( \mu ) 是常数。但真实铁磁材料有磁饱和和磁滞。硅钢片的 ( B-H ) 曲线在低场区斜率大、高场区平坦一进饱和区( \mu ) 急剧下降。边界条件公式本身是普适的但用它推导电感、磁阻时要小心。如果铁芯已经饱和等效磁导率 ( \mu ) 不是常数磁阻计算就变成非线性问题了。商用有限元软件求解时用的都是迭代方法牛顿-拉夫逊迭代、欠松弛迭代等每步迭代更新 ( \mu )。如果你只是手算近似最好先查一下材料的 ( B-H ) 曲线估出工作点的 ( \mu )再代入公式。这也是为什么很多经验不足的人拿线性公式去算饱和电感结果误差大到离谱。6. 静电静磁大类比一套思路解两类场6.1 标量势与矢量势的统一视角静电场和静磁场求解路径高度相似项目静电场静磁场源电荷密度 ( \rho )电流密度 ( \mathbf{J} )场量( \mathbf{E}, \mathbf{D} )( \mathbf{B}, \mathbf{H} )无旋/无散约束(abla \times \mathbf{E}0 )(abla \cdot \mathbf{B}0 )引入的势标量电势 ( \varphi )矢量磁位 ( \mathbf{A} )本构关系( \mathbf{D}\varepsilon \mathbf{E} )( \mathbf{B}\mu \mathbf{H} )泊松方程(abla^2\varphi-\rho/\varepsilon )(abla^2\mathbf{A}-\mu\mathbf{J} )边界量电容 ( C )电感 ( L )这张表是电磁场入门的“武功秘籍”。你可以把静磁问题完全翻译成静电场问题来记反之亦然。当年备考时我就是靠这张表把两大类题目打通了。碰到一个陌生结构先想它对应的是电容还是电感然后用类比直接套公式。6.2 电路磁路对应从电阻到磁阻电路和磁路的类比也非常经典电路磁路电动势 ( \mathcal{E} )磁动势 ( FNI )电流 ( I )磁通 ( \Phi )电阻 ( R\rho l/S )磁阻 ( \mathcal{R}l/(\mu S) )电导 ( G )磁导 ( \Lambda )这个类比在工程上极其好用变压器、电感器的磁路估算基本都是这个思路。但千万注意一个本质差异电路里导线约束电流沿导线流动断路就没有电流磁路里没有“磁绝缘体”磁通不只在铁芯里走还会漏到空气中。气隙越大漏磁越严重。所以磁路计算永远是工程近似精确结果还得靠数值求解。6.3 类比教你少背公式真正掌握这个类比之后你会发现很多公式根本不用背。比如磁场储能公式 ( W_m \frac{1}{2}L I^2 )你可以从电场储能 ( W_e \frac{1}{2}C V^2 ) 类比过来磁通与磁动势的关系 ( \PhiF/\mathcal{R} )和电流与电动势的关系 ( I\mathcal{E}/R ) 完全同构。但是有一个地方类比会失效静电场有正负电荷电场线有起有落磁场没有磁荷磁力线永远闭合。所以碰到“磁荷”模型某些永磁体简化模型会用等效磁荷代替一定要谨慎它只是数学等效不代表真实存在磁单极子。永磁体建模时等效磁荷模型算外部磁场很快但一旦涉及铁磁材料耦合误差就上来了。7. 实操中的坑从公式到仿真的常见问题7.1 常见错误速查表错误表现根本原因解决办法算电感老是对不上实测值忽略导体内部自感或漏磁路径算全磁通而不是只算主磁通铁芯磁通在仿真里“不守恒”网格剖分太粗或边界条件不正确加密气隙附近网格检查矢量磁位边界气隙磁密算出来比理论大忽略边缘效应或用了线性 ( \mu )设置饱和B-H曲线建三维模型看边缘误以为“静态没有电流”概念混淆记住恒定电流是静态场的源无限大模型仿真结果不对求解域太小截断边界太近远场加渐变吸收边界或扩大求解域面电流边界条件设置后结果跳变方向或大小正负取错用右手定则重新核对电流方向7.2 仿真软件设置里的几点体会我最早用有限元软件做静磁分析时总是急着导入模型、加电流、跑求解结果出来的磁力线图好看到不行一对比磁通数值发现差得离谱。后来才明白静磁仿真最关键的从来不是求解器而是边界条件和材料属性。第一个要点求解域截断边界。静磁问题的磁力线是闭合的它们会延伸到无穷远。有限元法必须把域截断边界远一些更准但代价是网格量大增。工程折中是对开放电磁结构比如螺线管外部磁场用“气球边界条件”或者“渐开远场单元”对封闭结构变压器铁芯则直接以铁芯外表面附近的空气层作为边界适当加厚即可。第二个要点电流源必须是闭合回路。静磁场求解器里如果一个电流源是“断”的只有一段导体没有回流路径解出来磁场会非常奇怪。因为 ( abla \times \mathbf{H} \mathbf{J} ) 要求电流密度是无散的断开的电流意味着电荷在堆积那就不是“静态”问题了。很多新手画完一根导线“从外面接进模型再从模型接出去”容易忽略模型内部的回流路径导致求解器警告电流不闭合。遇到这个问题别慌把端部用理想导体连起来形成闭合回路即可。第三个要点二维近似一定确认对称性。电机和电感横截面通常可以用二维仿真但前提是结构沿轴向无限长或者忽略端部效应。对短绕组、小型铁氧体磁芯端部漏感占比大二维结果明显偏高。我做小型共模电感时对比过二维结果往往比实测电感高10%到30%就是因为端部漏磁不计。如果没有条件建三维模型至少用二维结果乘以一个端部修正系数经验值在0.7到0.9之间具体看结构长宽比。7.3 手算与仿真如何互相验证现在工具发达了很多人直接拉个软件模型就跑边界条件设完也不细想。但我还是建议拿到一个静磁问题先用手算把关键量估一遍再开软件。手算的速度虽然慢但能让你建立“量级感”。比如一个环形铁芯电感先按 ( L N^2 \mu A / l ) 粗算得到大概1mH还是100mH心里有个数仿真出来后看是不是同一个量级。如果差了一个数量级说明哪里肯定错了——要么是有效截面积给错要么是材料 ( \mu ) 设错要么是磁路长度取错。手算的价值在于校验方向而不在于替代精确解。我自己带过几个做电源的工程师做电感仿真发现他们最喜欢直接抄别人模型的 B-H 曲线参数抄完换材料牌号也不改数据仿真结果自然离谱。B-H 曲线设置错了别的都白搭。材料属性是最该较真的地方。8. 最后分享一点实用小经验做这一行久了我越来越觉得学麦克斯韦方程组不在于会背那四个方程而在于你能不能快速判断一个实际问题该用哪一半的方程组。拿到一个电磁场问题第一件事永远是问场是时变的吗如果是老老实实上全套方程如果不是再问是电场主导还是磁场主导。静态安培环路方程最容易被误用的场景就是把它用到高频电流上——高频变压器磁芯里的磁场虽然也服从安培环路定律但频率稍高集肤效应、涡流损耗全跑出来了模型完全不同。我手边常备一张“静态场问题快速判断卡”恒定电流回路→静磁场固定电荷分布→静电场低频交流且电流均匀→暂态场或似稳场准静态场高频电磁波→全波场。每次开仿真前过一遍这张卡能省下大量试错时间。下一次有空我打算把准静态场和涡流问题的求解心得整理一下那是在静态场和全波场之间的“灰色地带”也是工程上绕不开的难题。有需要的朋友可以先在留言区说说自己遇到的具体案例我挑有代表性的补充进去。