永久磁铁的磁场分布
永久磁铁是一种能够在无外加磁场的情况下保持磁化状态的材料。其内部的磁畴在外部磁化后趋于一致排列，形成宏观的磁偶极子。磁场 B（磁感应强度）描述了磁铁对空间中其他磁性物体的作用力，其方向由磁感线指向，大小用特斯拉（T）或高斯（G）表示。
- 磁感应强度 B：矢量量，单位 T（1 T = 10⁴ G）。
- 磁场强度 H：与磁介质的磁化强度 M 关系为 B = μ₀(H + M)，其中 μ₀ 为真空磁导率。
- 磁偶极子：永久磁铁在宏观上等效为一个磁偶极子，产生的磁场可用偶极子模型近似描述。
永久磁铁的磁化方式
- 热磁化：在高温下将材料加热至居里温度以上，随后在外加磁场冷却，使磁畴重新排列。
- 机械磁化：通过冲击、压制等机械手段改变晶格缺陷，诱导磁畴取向。
- 电磁磁化：在强磁场中对材料进行磁化，是最常用的工业方法。
不同的磁化方式会影响磁畴的稳定性，从而影响磁场分布的均匀性和最大磁感应强度。
磁场分布的理论模型
1. 磁偶极子模型
对远离磁铁的空间，磁场近似为磁偶极子产生的场：
[
\mathbf{B}(\mathbf{r}) = \frac{\mu_0}{4\pi}\frac{3(\mathbf{m}\cdot\hat{r})\hat{r}-\mathbf{m}}{r^{3}}
]
其中 m 为磁偶极矩，(\hat{r}) 为指向观察点的单位向量，(r) 为距离。该模型适用于 r ≫ 磁体尺寸 的情况。
2. 磁单极子（磁荷）模型
在磁体内部或靠近表面时，可将磁体视为等效的正负磁荷分布。对矩形块磁铁，正负磁荷分别集中在两端面：
[
\sigma_m = \mathbf{M}\cdot\mathbf{n}
]
其中 (\sigma_m) 为面磁荷密度，(\mathbf{M}) 为磁化强度，(\mathbf{n}) 为外法向量。利用库仑型积分可求得任意点的磁感应强度。
3. 有限元数值模拟
当磁体形状复杂或材料非均匀时，解析解难以获得。此时采用 有限元法（FEM），在软件（如 ANSYS Maxwell、COMSOL）中建立几何模型、设定材料磁化曲线，求解 ∇·B = 0 与 ∇×H = 0 的耦合方程，得到高精度的磁场分布图。
常用测量方法
| 方法 | 原理 | 适用范围 | 优缺点 |
|---|---|---|---|
| 霍尔探头 | 利用霍尔效应产生与磁感应强度成比例的电压 | 直流/低频磁场，局部测量 | 体积小，分辨率受温度影响 |
| 磁通门（Fluxgate） | 通过磁芯磁化/去磁产生交变感应电压 | 中低磁场（µT~mT） | 高灵敏度，需校准 |
| 磁力计（Gaussmeter） | 直接读取磁感应强度 | 现场快速检测 | 精度有限，需靠近磁体 |
| 核磁共振（NMR） | 核自旋共振频率与磁场强度线性关系 | 高精度实验室测量 | 设备昂贵，受环境限制 |
测量时应注意探头的空间分辨率、温漂以及磁体本身的磁化退化。
实际应用中的磁场分布分析
1. 永久磁铁电机
在永磁同步电机（PMSM）中，转子磁铁的磁场分布决定了电机的电磁转矩波形。通过 FEM 计算转子表面的 B_r（径向磁感应）与 B_t（切向磁感应），可优化磁铁排布、槽形设计，降低转矩波动。
2. 磁悬浮列车
磁悬浮系统利用永磁体产生的强磁梯度实现无接触支撑。关键在于 磁场梯度 ∂B/∂z 的均匀性和峰值，过大的梯度会导致悬浮不稳，过小则支撑力不足。通过磁场分布仿真，可确定磁体的尺寸、间距以及补偿线圈的电流。
3. 磁共振成像（MRI）
MRI 的主磁场通常由超导磁体提供，但局部的 梯度磁体 多采用永磁体。梯度线圈的磁场均匀性直接影响成像分辨率，需要在设计阶段对永磁体的磁场分布进行精细调校。
示例计算：矩形块磁铁的中心磁感应强度
已知磁铁尺寸：长 a = 30 mm，宽 b = 20 mm，厚 c = 10 mm；磁化方向为 z 轴，磁化强度 M = 1.2 × 10⁶ A/m。求磁铁中心点（即 z = 0）的 B_z。
- 计算等效面磁荷密度
[
\sigma_m = M \quad (\text{正面});,; -M \quad (\text{负面})
]
- 对称性表明中心点只受 z 方向分量贡献，使用积分公式：
[
B_z = \frac{\mu_0}{4\pi},2M\iint_{-a/2}^{a/2}\iint_{-b/2}^{b/2}
\frac{c}{\bigl(x^{2}+y^{2}+ (c/2)^{2}\bigr)^{3/2}} ,dx,dy
]
- 采用数值积分（MATLAB / Python）可得：
import numpy as np
mu0 = 4*np.pi*1e-7
M = 1.2e6
a, b, c = 0.03, 0.02, 0.01
def integrand(x, y):
return c / ( (x**2 + y**2 + (c/2)**2)**1.5 )
xs = np.linspace(-a/2, a/2, 400)
ys = np.linspace(-b/2, b/2, 400)
dx = xs[1]-xs[0]; dy = ys[1]-ys[0]
X, Y = np.meshgrid(xs, ys)
Bz = mu0/(4*np.pi) * 2*M * np.sum(integrand(X, Y))*dx*dy
print(f"Bz = {Bz:.3e} T")
运行结果约为 B_z ≈ 0.48 T，与实际测量值相符，说明简化的面磁荷模型在中心区域具有良好精度。
小结
永久磁铁的磁场分布是磁学设计的核心问题。通过磁偶极子模型、面磁荷模型以及有限元数值仿真，可以在不同尺度上准确描述磁场特性。结合合适的测量手段，能够对磁体的实际磁场进行验证与校正，从而在电机、磁悬浮、医学成像等领域实现高效、可靠的应用。掌握上述理论与方法，是进行磁悬浮技术及其他磁场相关工程设计的必备基础。