当年第一次把“磁感应强度B、磁场强度H、矢量磁位A”这三个符号并排写在笔记本上时,我整个人是懵的。明明是同一个磁场,为什么要搞出三个物理量?H的单位是A/m,B的单位是T,A的单位又是Wb/m,三者之间还互相用公式绑着,稍不留神就会把概念记串。后来做了几年电磁场相关的工作,再回头看这段学习经历,才慢慢明白:B、H、A不是三个互相独立的场,而是描述同一个磁场的三种不同视角。这篇笔记,就是想把这三者的关系、来历、用途,以及我踩过的坑一次说清楚。
不管你是电气工程、物理专业的学生,还是做电机、变压器、电磁兼容、感应加热的工程师,只要需要和磁场打交道,这篇文章都值得看完。我会从头梳理B、H、A各自的物理意义,用几个手算案例把它们串起来,再讲一讲有限元仿真里到底该看哪个量。内容偏基础,但里面有不少是我实际干活时才悟出来的细节。
1. 为什么同一个磁场,要同时定义三个量?
1.1 三个量的“分工”完全不同
先说结论:B回答“磁场对电荷和电流产生多大的力”,H回答“电流这个源在空间中激发出多强的磁场”,A回答“这个磁场如何用一个位函数来表达,方便计算”。这三件事,恰好对应电磁场理论里的三个需求层次。
- B是受力意义上的磁场。运动电荷在磁场中受洛伦兹力F=qv×B,载流导线受安培力,磁通量Φ=∫B·dS,这些都是B在起作用。所以B是“结果”,是你最终关心的物理量。
- H是源意义上的磁场。安培环路定理的微分形式是∇×H=J,也就是说H的旋度直接等于传导电流密度。给你一套线圈,通上电流I,绕一圈算一下∮H·dl就等于电流,这事跟介质材料是什么关系不大,H直接由源决定。
- A是数学意义上的辅助量。由于B永远是无散场(∇·B=0),可以写成某个矢量场的旋度,这个矢量场就是A。它不是直接可测量的量,但在计算和仿真里特别好用。
这三个量之间由本构关系B=μH连着,而A和B之间则是B=∇×A。理解了这个“三角关系”,后面所有内容都顺了。
1.2 历史账里藏着为什么容易搞混
很多教材会提到,H这个量其实比B更早被定义。早期研究磁学时,人们把磁极之间的相互作用力类比于电荷之间的库仑力,定义了一个“磁场强度”,单位用奥斯特,也就是今天H的前身。后来安培发现磁场的本源是电流,法拉第又提出磁通、磁力线的概念,麦克斯韦统一电磁理论后,B作为“磁感应强度”承担了受力计算的重任。
但历史留下来的名称却让人很头疼:中文里“磁场强度”这四个字,直觉上应该指磁场本身的强弱,可它偏偏给了H;真正决定磁场强弱的B反而叫“磁感应强度”。我在给新人讲课时经常强调一个心得:不要按名字去猜物理意义,直接按公式去记忆。H负责“算源”,B负责“算力”,这个分工比名字重要得多。
1.3 三种视角的适用场景
工程里最常见的分工是这样的:
- 计算磁路的励磁安匝数、线圈电流、气隙磁动势时,几乎都用H。因为H沿着磁路的线积分就是安匝数N×I,直接对应设计输入。
- 判断铁芯会不会饱和、计算受力、计算磁通密度限制时,用B。因为B直接决定铁芯材料的磁化状态和磁通量。
- 求解电磁场分布、数值仿真、计算感应电场时,A反而是主力。尤其涡流场问题,电流密度和A直接相关,软件内部求解的未知量也经常就是A。
所以你看,这三个量不是“谁替代谁”的关系,而是各自管辖一片区域。下面我把每个量单独拆开细讲。
2. 三个量的定义、单位与物理意义,别再靠背公式理解
2.1 磁感应强度B:真正告诉你磁场“有多强”的量
B的定义来自洛伦兹力:一个电荷q以速度v在磁场中运动,受到的磁力为F=qv×B。这个力垂直于速度和磁场构成的平面。从力的意义上,B才是那个“磁场实力值”。单位特斯拉(T),1T=1N/(A·m),形象点说,1T的磁场对一根1米长、通1A电流的导线产生1N的力。
在SI单位制下,磁通量Φ=∫B·dS,单位是韦伯(Wb),因此B也可写成Wb/m²。高速磁悬浮、MRI、粒子加速器里面算磁场强弱,说的都是B。比如地球表面地磁约为30-60μT,普通永磁体表面约0.2-0.5T,电机气隙磁密通常设计在0.7-1.2T,磁共振成像设备能达到1.5T甚至3T。这些数值级别的概念,用的全是B。
B还有一个极其重要的微分性质:∇·B=0。这意味着磁力线永远闭合,没有“磁单极子”,这也是后面引入A的原因所在。任何关于B的边界条件、磁路计算,源头都能追到这个散度为零。
2.2 磁场强度H:把“媒介影响”剥离出来的助手
H的定义表面上没有B那么直观,但它解决了一个实际问题:同样一个电流线圈,放在真空里和套在铁芯上,B差了上千倍,因为铁芯内部会产生很强的磁化场。这样一来,用B去反推“电流到底贡献了多少”就很麻烦。H的出现把这些介质贡献剥离开来,让电流这个源和磁场响应之间的关系变得干净。
在均匀线性介质中,本构关系是B=μH,其中μ=μ0μr,μ0=4π×10⁻⁷ H/m是真空磁导率,μr是相对磁导率。更一般地,H=B/μ0-M,M是磁化强度。铁磁性材料里μr可以从几百到几万,这也是为什么变压器铁芯能大幅增强磁通。
H的微分性质是∇×H=J(静磁情况下),积分形式也就是安培环路定理∮H·dl=ΣI。用H去计算一个线圈在磁路里产生的“磁势”,实在太顺手了。磁路欧姆定律里的磁动势F=NI,磁阻Rm=l/(μS),最后磁通Φ=NI/Rm,链路里的核心变量就是H。
2.3 矢量磁位A:被低估的“中间人”
矢量磁位A的引入,纯粹是数学上的一步棋。既然B的散度恒等于零,而任意矢量场旋度的散度也恒等于零,那么只要令B=∇×A,自动就能满足∇·B=0。这就把一个矢量场的求解问题,转化成另一个矢量场的求解问题。听起来像原地打转,但在很多场景里是降维打击。
先看A的单位:Wb/m,也就是T·m。它在静磁问题中满足的方程是∇²A=-μJ(在库仑规范和均匀介质下)。这个形式和静电场的泊松方程∇²φ=-ρ/ε完全同构。也就是说,只要你熟悉电荷如何产生电位,就能照着理解电流如何产生矢量磁位,只是这里每个方向的分量都要单独解一遍。
A不是直接可测量,它的物理意义藏在环量里:根据斯托克斯定理,∫∫B·dS=∮A·dl。也就是说,A沿着某个闭合回路的线积分,等于穿过这个回路的总磁通。这个关系在计算电感、磁通耦合时特别关键,后面第4节的实例里会看到。
2.4 一张表说清三个量的关系
我整理了一张自己一直在用的对照表,基本上把新手需要记忆的核心信息都塞进去了:
| 物理量 | 磁感应强度B | 磁场强度H | 矢量磁位A |
|---|---|---|---|
| 中文俗称 | 磁通密度 | 磁场强度 | 磁矢位 |
| 单位 | T(特斯拉),即Wb/m² | A/m | Wb/m |
| 主要作用 | 计算受力、磁通、饱和 | 计算电流源、磁动势、磁路 | 位函数,求场、算磁通、数值仿真 |
| 微分关系 | ∇·B=0 | ∇×H=J | B=∇×A |
| 积分关系 | Φ=∫B·dS | ∮H·dl=I | ∮A·dl=Φ |
| 边界条件 | 法向连续 | 切向连续(无面电流时) | 切向连续(规范选择有讲究) |
| 本文中角色 | “力”视角 | “源”视角 | “计算”视角 |
这张表建议收藏,做题前先看一遍,能把很多公式记混的问题提前消灭掉。
3. 矢量磁位A为什么要引入?这里藏着最大的一处理解门槛
3.1 从“散度为零”凑出来的旋度场
很多初学者不理解,既然B已经能描述磁场,为什么还要自找麻烦定义A?我换个角度解释你就能通:在实际求磁场时,你经常需要面对复杂的电流分布和边界形状,直接解B需要同时满足两个方程——散度方程和旋度方程,还是挺麻烦的。但如果我先定义一个辅助量A,令B=∇×A,那么B的散度方程自动满足,剩下只需要解旋度方程,而旋度方程在库仑规范下又能变成泊松方程。
这个过程,就相当于把一个二阶偏微分方程组的求解,转化为几个独立的泊松方程,把一个矢量场B的三个未知分量,转化为A的三个未知分量。看起来未知数没变,但数学结构简单了,边界条件的处理方式也更统一。更重要的是,在二维平面场问题里,A只有一个方向的分量,形式大大简化。
3.2 规范自由度:A不唯一,但磁场唯一
这里有个容易忽略的坑:给定B,A并不是唯一的。如果你把A换成A'=A+∇ψ,其中ψ是任意标量函数,那么因为旋度的旋度为零(∇×∇ψ=0),B一点都不会变。这种变换叫规范变换,A的取值有一定的“自由发挥”空间。
为了让求解结果唯一,我们需要人为加一个约束,最常用的是库仑规范∇·A=0。加了它之后,A的泊松方程才能从∇×(∇×A)=μJ整理成∇²A=-μJ。实际做解析计算时,A的参考点怎么定是一个很自由的选择,只要保证最终算出的B不变就行。
我在第一次算无限长直导线的A时,就因为参考点选择不当导致结果里出现发散项,折腾了半天。后来经验是:选哪个半径为参考点本身不重要,关键是你的问题边界允许你选哪里。数值计算中,软件会自己处理规范,但在手算时你得心里有数,否则会得出“同一个B对应两个不同的A”这种让自己困惑的结果。
3.3 A的边界条件为什么很少被讲透
B和H的边界条件大家多少都能记住,A的边界条件却常被一笔带过。但实际做解析解或有限元时,A的边界条件直接决定求解成败。
在两种介质分界面上,如果不存在面电流,A的切向分量必须连续。这是为了保证磁通连续。而在库仑规范下,A的法向分量也连续。许多教材直接写“A连续”,背后其实已经默认选好了规范。如果二维问题中A只有轴向分量,比如z轴方向,那么在垂直于z轴的截面内,A就是一个标量函数,边界条件退化为“A本身连续”。这也是为什么很多电磁场数值教材会用“A_z”举例。
3.4 二维场景下A的等值线就是磁力线
A最大的“隐藏福利”出现在二维静磁问题里。当电流都沿z轴方向流动,磁场分布在x-y平面内时,A只有z方向分量A_z(x,y),B_x=∂A_z/∂y,B_y=-∂A_z/∂x。这种形式下,A_z的等值线就是磁力线。
这是我看懂A的转折点。以前总觉得A是个抽象得没法想象的量,后来在有限元软件后处理里看到,那些漂亮的“磁力线图”本质上就是A_z的等值线图。磁力线密集的地方,对应A_z的梯度大,也就是B强;A_z相同的点连成的线,正好是磁场方向。每当你疑惑A到底有什么用,就想想那些磁力线云图——它的底子就是A。
4. 三道手算题,把B、H、A从公式变成手感
4.1 无限长直导线的B/H/A:同一根导线,三个量三种画法
先来最经典的例子。一根半径为a的无限长直导线,通有电流I,沿z轴方向。真空中,求导线外部和内部的B、H、A。
- 外部(r>a):根据安培环路定理,∮H·dl=I,得H=I/(2πr)(方向沿圆周切向),B=μ0I/(2πr)。A只有z分量,由Bφ=-∂A_z/∂r可得A_z=-(μ0I/2π)ln(r/a)+C。取导线表面r=a处A_z=0,则C=0,外部A_z=-(μ0I/2π)ln(r/a)。
- 内部(r<a):电流密度J=I/(πa²),∮H·dl=I(r²/a²),得H=Ir/(2πa²),B=μ0Ir/(2πa²)。再积分一次,并让r=a处内外A_z连续,可得A_z=-(μ0I/4π)(r²/a²-1)。
算完之后你会发现:H和B在导线表面处都是最大的(外部随1/r衰减),而A在导线表面取零、轴线处为正,这说明A的取值高低本身没有绝对意义,有意义的是A的梯度。这比任何语言都直观地解释了“A是位函数”这个说法。
4.2 螺线管内部的B/H:均匀场不等于均匀A
第二个例子是无限长密绕螺线管,单位长度匝数n,电流I。理想情况下,管内磁场均匀,管外磁场为零。用安培环路定理取一个矩形回路,可得管内H=nI,B=μ0nI。
很多人到此就停了,但我建议再算一步A。取柱坐标,因为管内B只有z方向的均匀值B0,由B=∇×A,可知A只有圆周方向分量Aφ,并且满足(1/r)∂(rAφ)/∂r=B0,积分得Aφ=(B0r)/2(加上常数可以忽略)。结果是:哪怕B完全均匀,A仍然随半径线性变化。这再次说明A的绝对大小没有物理意义,它的空间变化率才有。
这个例子在做感应电机、加速器等分析时特别有用,因为很多工程公式都涉及“在均匀磁场中运动的导体感受到的感应电场”,那里面的矢量磁位A就会出现。
4.3 变压器的电感计算链路:从电流到H到B再到电感
最贴近工程的一道题:环形铁芯线圈,铁芯平均磁路长度l=0.2m,截面积S=4×10⁻⁴m²,匝数N=100,电流I=1A,相对磁导率μr=2000,算电感。
计算链路是这样的:
- 磁场强度:H=NI/l=100×1/0.2=500A/m;
- 磁感应强度:B=μ0μrH=4π×10⁻⁷×2000×500≈1.256T;
- 磁通:Φ=BS=1.256×4×10⁻⁴≈5.02×10⁻⁴Wb;
- 磁链:Ψ=NΦ=100×5.02×10⁻⁴≈0.0502Wb;
- 电感:L=Ψ/I=0.0502/1≈50.2mH。
如果铁芯抽掉,μr=1,电感大概只有25μH,相差超过2000倍。这条链路很好地展示了为什么设计磁路时要先用H:H完全由NI和磁路长度决定,先算出H,再由材料磁导率得到B,最后才得到磁通和电感。B和H在这一步步里各司其职,谁也不能替代谁。
4.4 算完这几道题后,我对三者关系的理解
这三道题算下来,我的体会有两点。第一,H和B的关系就像“电压”和“电流”:电压(H)是源驱动,电流(B)是介质响应,中间夹着材料的“电导率”(磁导率)。第二,A不是用来替代B的,它是B的“原函数”,就像电位是电场的“原函数”一样。做解析计算时,有些场分布先求A比直接求B容易;做仿真时,软件内部更是直接把A当未知量。
5. 工程应用里,哪些场合该信B,哪些场合该信H,哪些问题离不开A
5.1 磁路设计:磁动势NI和H,比B更接近“设计输入”
设计变压器、电感、电机磁路时,工程师拿到的最初参数一定是电流和匝数。这时H的价值立刻凸显:沿着磁路把H积分一圈,就等于总的安匝数NI。气隙、铁芯段、永磁体段可以像电阻串联一样处理,每段磁路的磁阻R=l/(μS),总磁动势等于Φ×ΣR。
即便你在做永磁电机,也需要通过H描述永磁体的工作点。永磁体的退磁曲线上横坐标就是H,纵坐标是B。设计时要在B-H曲线上找到合适的负载线,使永磁体工作在最大磁能积附近。如果只盯着B,不考虑H,你根本没法判断永磁体会不会退磁。
5.2 磁滞回线与损耗计算:B和H必须成对看
材料手册里的磁滞回线,横轴是H,纵轴是B。这是因为H代表外加励磁,B代表材料响应。铁磁材料在交变磁化时,B的变化滞后于H,形成回线,回线包围的面积就代表每个周期的磁滞损耗,数学上对应∮H dB。
所以你在做铁损计算时,绝对不能只取一个量。励磁电流确定H的幅值,材料饱和特性决定B的幅值。用正弦波近似时,铁损公式P=k f^α B_m^β中的B_m是磁感应强度峰值,但算励磁电流时又要回头用H。这种来回切换,正是B和H并存的实际意义。
5.3 涡流场和感应问题:A的时间导数直接出场
到了时变场,A的重要性彻底暴露。电磁感应定律的完整形式是E=-∂A/∂t-∇φ。也就是说,感应电场并不仅仅由电位梯度∇φ决定,还有A的时间导数。在涡流场问题中,导体里的电流密度J=σE,其中很大一部分来源就是-σ∂A/∂t。
这也是为什么感应加热、无损检测、电机涡流损耗分析等高阶应用,几乎都建立在求解A的基础上。你答不出A不代表算不了静态场,但一旦涉及时变场,A就绕不过去了。那些声称“B直接就能搞定一切”的直觉,到这里基本失效。
5.4 有限元软件内部:求解变量、后处理量各自是什么
用Ansys Maxwell、COMSOL、JMAG这类低频电磁场软件时,你可能发现求解器设置里有个自由度叫“Vector Potential”或A_z。二维静磁问题里,软件求解的就是A_z这个标量,然后自动后处理成B和H云图。三维问题常用A-V(矢量磁位加标量电位)公式,或者T-Ω公式,这背后的原因就是第3章说的:你不直接解B的方程,而是解一个更容易加边界条件的位函数方程。
所以,你在后处理里看到的磁力线图,本质上是A_z等值线;看到的B云图,本质上是A_z梯度再乘以材料属性;看到的H云图,又是在B基础上除以磁导率得到的。理解这个链条后,你就能明白为什么网格质量差时B云图会出现“锯齿”:因为软件需要对A做数值微分,一阶单元里B是分段常数,云图当然不平滑。
6. 学习电磁场量时最常踩的坑与几点实操建议
6.1 坑1:边界条件记反
最常见的错误是把界面条件记成“B的切向连续、H的法向连续”。真实情况:磁感应强度B的法向分量连续,磁场强度H的切向分量连续(不存在面电流时)。一旦记反,遇到铁芯-空气界面就会得出完全错误的结果。
一个特别实用的验证方法是:铁芯和空气界面上,如果磁场垂直穿入铁芯表面,那么B法向连续意味着两边B的垂直分量一样,但铁芯里H的垂直分量小得可怜;如果磁场平行于表面,H切向连续意味着两边H的平行分量一样,但铁芯里B的平行分量大得惊人。正因为这样,铁芯表面附近的磁力线总是“歪”的,B会在界面上突然转向。画场图时,这个转向方向就能反过来检验你有没有记反。
6.2 坑2:看到旋度就怕,其实旋度只是“绕圈程度”
很多初学者被∇×B、∇×H这些符号吓住。可以把旋度理解成“场线绕某点旋转的倾向”。安培环路定理∮H·dl=I就是旋度定理的积分形式,它说得很直白:H沿着一个圈绕一圈,积分结果等于穿过这个圈的电流。电流越大,H的“绕圈”越明显。
矢量磁位A之所以经常被说“难”,也是因为大家一看到B=∇×A就心里打鼓。但如果你已经算过第4.1节的直导线,会发现这个旋度运算就是一次求导:Bφ=-∂A_z/∂r。理解了旋度是“导数家族”的一员,很多公式就不那么吓人了。
6.3 坑3:仿真云图很漂亮,却不知道软件是拿什么算出来的
我见过不少新人直接把软件生成的B云图当真理,完全不看求解设置。实际上,软件里B和H的精度高度依赖网格密度和单元阶次。一阶三角形/四面体单元中,A_z是线性变化的,对它求导得到的B在单元内是常数,云图会呈现明显的色块分层。这时候加密网格,色块分层会减轻,但计算时间也会暴涨。
正确习惯是:先看A云图是否平滑,再看B云图是否有异常跳变。求解完成后,又在关键路径上画一条线,查看B的法向和切向分量变化。如果法向分量在介质界面处发生跳变,而你明明没有设置面电流,那基本可以断定是网格不够密或边界条件设置错了。
6.4 学习顺序与自检清单
如果你正在学这部分内容,我的建议是不要一上来就死磕麦克斯韦方程组,先按下面这个顺序走:
- 从洛伦兹力公式理解B的物理含义,记住单位T。
- 用安培环路定理理解H,记住H的环量等于电流。
- 用B=μH把B和H联系起来,做几道磁路计算题。
- 再接触A,理解B=∇×A和库仑规范,做无限长直导线和螺线管两道解析题。
- 最后用有限元软件跑两个例子,分别看B、H、A、磁力线四种结果的差异。
每走完一步,可以拿第2.4节那张表做一次自检:能不能不看资料说出每个量的单位、边界条件、核心公式?如果有一项卡壳,就回头重看对应章节。
对我来说,真正打通B、H、A这三个量,不是在某一次课堂上学到的,而是在无数次“算错了又回头查哪里记混”的过程中磨出来的。现在做工程分析,我先用H估算励磁和磁动势,再用B判断饱和和受力,遇到时变场和仿真时想到A。这个思维链条一旦形成,后面再学磁路设计、电磁兼容、无线充电,都会有稳稳的地基。