赫姆霍兹线圈与三线圈均匀磁场原理

赫姆霍兹线圈是一种用于产生近似均匀磁场的经典结构。它由两个完全相同的同轴圆形线圈组成，两个线圈通以相同方向、相同大小的电流。当两个线圈中心间距等于线圈半径时，中心区域的磁场会变得非常平坦。

如果进一步加入第三个线圈，合理选择线圈位置和安匝比，可以比普通赫姆霍兹线圈获得更大的均匀磁场区域。本文从单个圆线圈的轴向磁场出发，推导赫姆霍兹双线圈的均匀性条件，并说明三线圈结构的改进原理。

单线圈轴向磁场

设一个圆形线圈半径为 $R$，匝数为 $N$，电流为 $I$。线圈位于 $xy$ 平面，圆心在原点，轴线为 $z$ 轴。

由毕奥-萨伐尔定律：

$d\mathbf{B} = \frac{\mu_0 I}{4\pi} \frac{d\boldsymbol{\ell}\times \mathbf{r}}{r^3} $

对于圆形线圈，在轴线上距离圆心为 $z$ 的点，磁场方向沿 $z$ 轴。积分后得到单个圆线圈的轴向磁场：

$B(z)=\frac{\mu_0 N I R^2}{2(R^2+z^2)^{3/2}}$

其中 $\mu_0$ 为真空磁导率，$N$ 为匝数，$I$ 为电流，$R$ 为半径，$z$ 为观察点到线圈中心的轴向距离。在线圈中心 $z=0$ 处：

$ B(0) = \frac{\mu_0 N I}{2R} $

这是分析赫姆霍兹线圈和多线圈系统的基础。

赫姆霍兹双线圈

设两个相同线圈同轴放置，半径均为 $R$，匝数均为 $N$，电流均为 $I$，中心分别位于 $ z = -\frac{s}{2} $ 与 $ z = +\frac{s}{2} $，$s$ 为两线圈中心间距。轴线上总磁场为：

$ B(z) = \frac{\mu_0 N I R^2}{2} \left[ \frac{1}{\left(R^2+(z-s/2)^2\right)^{3/2}}+ \frac{1}{\left(R^2+(z+s/2)^2\right)^{3/2}}\right] $

系统关于 $z=0$ 对称，中心附近可展开为 $ B(z) = B_0 + A_2z^2 + A_4z^4 + A_6z^6 + \cdots $。一阶项因对称性为零；令 $ A_2 = 0 $ 可得经典条件 $ s = R $（间距等于半径）。

当 $ s=R $ 时，两线圈位于 $ z=-\frac{R}{2} $、$ z=+\frac{R}{2} $。中心磁场经化简为：

$ B_0 = \left(\frac{4}{5}\right)^{3/2} \frac{\mu_0 N I}{R} \approx 0.7155 \frac{\mu_0 N I}{R} $

此时 $ A_2=0 $，展开主导项变为四阶：$ B(z) = B_0 + A_4z^4 + A_6z^6 + \cdots $，中心附近磁场曲线更平坦，即赫姆霍兹线圈近似均匀磁场的根本原因。

引入无量纲坐标 $ u=z/R $，标准双线圈位于 $ u=\pm \frac{1}{2} $，轴向场可写为：

$ B(u) = \frac{\mu_0 N I}{2R} \left[ \frac{1}{\left(1+(u-1/2)^2\right)^{3/2}} + \frac{1}{\left(1+(u+1/2)^2\right)^{3/2}} \right] $

中心场 $ B(0) = \left(\frac{4}{5}\right)^{3/2} \frac{\mu_0 N I}{R} $；归一化 $ B(u)/B(0) $ 在中心附近应尽量接近 1。

三线圈系统

任意多个同轴圆线圈可用叠加原理。第 $i$ 个线圈位于 $z=z_i$、半径 $R_i$、匝数 $N_i$、电流 $I_i$ 时：

$ B_z(z) = \sum_i \frac{\mu_0 N_i I_i R_i^2} {2\left[R_i^2+(z-z_i)^2\right]^{3/2}} $

在赫姆霍兹双线圈中间加入中心线圈，即常见的同半径三线圈：外侧位于 $z=\pm hR$，中心位于 $z=0$，外侧安匝 $NI$，中心安匝 $mNI$（$m$ 为安匝比）。轴向场为：

$ B(u) = \frac{\mu_0 N I}{2R} \left[ \frac{1}{\left(1+(u-h)^2\right)^{3/2}} + \frac{1}{\left(1+(u+h)^2\right)^{3/2}} + m \frac{1}{\left(1+u^2\right)^{3/2}} \right] $

双线圈仅可调间距，通常只能令 $ C_2=0 $；三线圈另有 $ h $、$ m $，可同时令 $ C_2=0 $ 且 $ C_4=0 $，使展开变为 $ B(u) = B_0 + C_6u^6 + \cdots $，均匀区更宽。

一种常见优化参数为 $ h\approx 0.76005 $、$ m\approx 0.5315 $，即 $ z=\pm 0.76005R $ 与中心 $ z=0 $，匝数比约 $59:111$。中心场约为 $ B_0 \approx 0.7704 , \mu_0 N_o I / R $（$N_o$ 为外侧匝数；比较场强时还需考虑总匝数、电阻与发热）。

结构	线圈位置	安匝比	主导误差	特点
单线圈	$z=0$	1	二阶	中心场强，均匀区小
赫姆霍兹双线圈	$z=\pm0.5R$	$1:1$	四阶	结构简单
优化三线圈	$z=0,\pm0.76005R$	$0.5315:1:1$	六阶	均匀区更大

工程与数值计算

磁场与安匝数成正比：$ B\propto NI $。功耗 $ P=I^2 R_c $（$R_c$ 为线圈电阻），长时间运行需考虑导线直径、散热、温升与电源能力。几何对准误差（间距、半径、匝数比、同轴度）会明显影响均匀性；实验中还须注意地磁场与铁磁干扰。

任意多线圈轴向场可用 Python 计算，例如：

import numpy as np

mu0 = 4 * np.pi * 1e-7

def axial_field(z, coils):
    B = 0.0
    for coil in coils:
        R = coil["radius"]
        N = coil["turns"]
        I = coil["current"]
        z0 = coil["z"]
        B += mu0 * N * I * R**2 / (2 * (R**2 + (z - z0)**2)**1.5)
    return B

R, I, N = 1.0, 1.0, 100
helmholtz = [
    {"radius": R, "turns": N, "current": I, "z": -0.5 * R},
    {"radius": R, "turns": N, "current": I, "z": +0.5 * R},
]
triple = [
    {"radius": R, "turns": N, "current": I, "z": -0.76005 * R},
    {"radius": R, "turns": 0.5315 * N, "current": I, "z": 0.0},
    {"radius": R, "turns": N, "current": I, "z": +0.76005 * R},
]

归一化 $ B(z)/B(0) $ 越接近 1，该处磁场越接近中心值。交互可视化中轴线场可用 JavaScript 的 axialB(z, coils) 累加各线圈贡献；空间任意点则用毕奥-萨伐尔数值积分：

$ \mathbf{B}(\mathbf{r}) \approx \frac{\mu_0 I}{4\pi} \sum_j \frac{ \Delta\boldsymbol{\ell}_j \times (\mathbf{r}-\mathbf{r}_j) } { |\mathbf{r}-\mathbf{r}_j|^3 } $

小结

赫姆霍兹线圈通过对称叠加使 $ A_2=0 $，在 $ s=R $ 时 $ B_0 = \left(\frac{4}{5}\right)^{3/2} \mu_0 N I / R $。三线圈在 $ z=0,\pm0.76005R $、安匝比 $ m \approx 0.5315 $ 下可进一步压低低阶误差，获得更宽的均匀磁场平台。

参考资料