我们考虑一个半径为 $a$、总长度为 $L$ 的密绕螺线管。导线沿圆柱面以等距螺旋线方式缠绕，螺距为 $h$（相邻两匝之间的轴向距离），总匝数 $N = L/h$。设导线中通有恒定电流 $I$，方向沿螺旋线参数增加的方向。我们采用等距螺旋线模型，通过毕奥-萨伐尔定律直接计算轴线上任意点的磁场轴向分量，并推广到无限长密绕螺线管内部均匀磁场的情况。
需要注意的是：一根有限长的螺旋导线本身并不具有严格的连续旋转对称性，因此不能简单地断言轴线上磁场只有 $z$ 分量。本文严格推导的是轴向分量 $B_z$。在密绕近似、连续面电流模型或周期平均意义下，横向分量可以认为相互抵消，此时轴上磁场可近似看作只沿轴向。
1. 螺旋线参数化
取圆柱坐标系，使螺线管轴线沿 $z$ 轴。螺旋线的参数方程可写为（$\varphi$ 为绕轴转角）：
$$
\mathbf{r}(\varphi) = (a\cos\varphi, a\sin\varphi, p\varphi), \quad 0 \le \varphi \le 2\pi N.
$$
其中
$$
p = \frac{h}{2\pi}
$$
是螺旋线沿 $z$ 方向的前进系数。也就是说，当 $\varphi$ 增加 $2\pi$ 时，导线绕轴一周，同时 $z$ 坐标增加
$$
2\pi p = h.
$$
因此 $h$ 是螺距，而 $p$ 是单位弧度对应的轴向前进距离。电流方向对应 $\varphi$ 增加的方向。
2. 毕奥-萨伐尔定律
轴线上一点 $P(0,0,z)$ 处的磁感应强度为
$$
\mathbf{B}(z) = \frac{\mu_0 I}{4\pi}\int \frac{d\mathbf{l}\times\mathbf{R}}{R^3}.
$$
式中 $d\mathbf{l}$ 是电流元矢量，$\mathbf{R}$ 是从源点到场点的矢量，$R = |\mathbf{R}|$。若场点位置矢量为
$$
\mathbf{r}_0 = (0,0,z),
$$
则
$$
\mathbf{R} = \mathbf{r}_0 - \mathbf{r}(\varphi).
$$
电流元
由螺旋线参数方程可得
$$
\frac{d\mathbf{r}}{d\varphi} = (-a\sin\varphi, a\cos\varphi, p).
$$
因此
$$
d\mathbf{l} = \frac{d\mathbf{r}}{d\varphi}d\varphi = (-a\sin\varphi, a\cos\varphi, p)d\varphi.
$$
相对位置矢量
从源点指向场点的矢量为
$$
\mathbf{R} = (0,0,z) - (a\cos\varphi, a\sin\varphi, p\varphi).
$$
即
$$
\mathbf{R} = (-a\cos\varphi, -a\sin\varphi, z-p\varphi).
$$
其模长为
$$
R = \sqrt{a^2 + (z-p\varphi)^2}.
$$
叉乘 $d\mathbf{l} \times \mathbf{R}$
下面直接写出叉乘各分量。由
$$
d\mathbf{l} = (-a\sin\varphi, a\cos\varphi, p)d\varphi
$$
和
$$
\mathbf{R} = (-a\cos\varphi, -a\sin\varphi, z-p\varphi)
$$
可得 $x$ 分量：
$$
(d\mathbf{l}\times\mathbf{R})_x = a\cos\varphi(z-p\varphi)d\varphi + ap\sin\varphi d\varphi.
$$
$y$ 分量：
$$
(d\mathbf{l}\times\mathbf{R})_y = a\sin\varphi(z-p\varphi)d\varphi - ap\cos\varphi d\varphi.
$$
$z$ 分量：
$$
(d\mathbf{l}\times\mathbf{R})_z = a^2\sin^2\varphi d\varphi + a^2\cos^2\varphi d\varphi.
$$
由于
$$
\sin^2\varphi + \cos^2\varphi = 1,
$$
所以
$$
(d\mathbf{l}\times\mathbf{R})_z = a^2 d\varphi.
$$
因此
$$
d\mathbf{l}\times\mathbf{R} = (a\cos\varphi(z-p\varphi)+ap\sin\varphi, a\sin\varphi(z-p\varphi)-ap\cos\varphi, a^2)d\varphi.
$$
3. 轴上磁场积分
本文重点计算轴向磁场分量 $B_z$。由毕奥-萨伐尔定律可得
$$
B_z(z) = \frac{\mu_0 I}{4\pi}\int_0^{2\pi N}\frac{a^2}{[a^2+(z-p\varphi)^2]^{3/2}}d\varphi.
$$
作变量代换
$$
u = z - p\varphi.
$$
则
$$
du = -p d\varphi.
$$
即
$$
d\varphi = -\frac{du}{p}.
$$
积分限随之变化：
- 当 $\varphi = 0$ 时，$u=z$；
- 当 $\varphi = 2\pi N$ 时，
$$
u = z - 2\pi Np.
$$
又因为
$$
p = \frac{h}{2\pi},
$$
所以
$$
2\pi Np = Nh = L.
$$
因此上限对应
$$
u = z - L.
$$
于是
$$
B_z(z) = \frac{\mu_0 I a^2}{4\pi}\int_z^{z-L}\frac{-du/p}{(a^2+u^2)^{3/2}}.
$$
调整积分上下限，得到
$$
B_z(z) = \frac{\mu_0 I a^2}{4\pi p}\int_{z-L}^{z}\frac{du}{(a^2+u^2)^{3/2}}.
$$
利用积分公式
$$
\int \frac{du}{(a^2+u^2)^{3/2}} = \frac{u}{a^2\sqrt{a^2+u^2}} + C.
$$
得到
$$
B_z(z) = \frac{\mu_0 I a^2}{4\pi p}\left[\frac{u}{a^2\sqrt{a^2+u^2}}\right]_{z-L}^{z}.
$$
化简得
$$
B_z(z) = \frac{\mu_0 I}{4\pi p}\left(\frac{z}{\sqrt{a^2+z^2}} - \frac{z-L}{\sqrt{a^2+(z-L)^2}}\right).
$$
代入
$$
p = \frac{h}{2\pi},
$$
有
$$
\frac{1}{p} = \frac{2\pi}{h}.
$$
故
$$
\boxed{B_z(z) = \frac{\mu_0 I}{2h}\left(\frac{z}{\sqrt{a^2+z^2}} - \frac{z-L}{\sqrt{a^2+(z-L)^2}}\right)}.
$$
这就是有限长等距螺旋线在轴线上产生的磁场轴向分量。
若引入单位长度匝数
$$
n = \frac{1}{h},
$$
则上式也可写为
$$
\boxed{B_z(z) = \frac{\mu_0 n I}{2}\left(\frac{z}{\sqrt{a^2+z^2}} - \frac{z-L}{\sqrt{a^2+(z-L)^2}}\right)}.
$$
4. 无限长螺线管的极限
当螺线管无限长时，$L \to \infty$，且考察内部远离端点的区域。按照本文的坐标约定，这相当于
$$
z \to +\infty,
$$
以及
$$
z-L \to -\infty.
$$
此时
$$
\frac{z}{\sqrt{a^2+z^2}} \to 1,
$$
并且
$$
\frac{z-L}{\sqrt{a^2+(z-L)^2}} \to -1.
$$
于是
$$
B_z = \frac{\mu_0 I}{2h}[1 - (-1)].
$$
因此
$$
B_z = \frac{\mu_0 I}{h}.
$$
由于单位长度匝数
$$
n = \frac{1}{h},
$$
最终得到
$$
\boxed{B_z = \mu_0 n I}.
$$
在无限长密绕螺线管的理想模型中，磁场方向沿轴线，因此可以写成矢量形式
$$
\boxed{\mathbf{B} = \mu_0 n I\hat{\mathbf{z}}}.
$$
这一结果与教科书中的无限长密绕螺线管内部均匀磁场公式一致。
5. 讨论
上述推导基于螺线管导线的实际螺旋形状，而不是一开始就把螺线管看作许多独立圆环的叠加。计算表明，对于轴向分量 $B_z$，螺旋线模型与圆环叠加模型给出的结果一致。
单根有限长螺旋导线并不严格具有轴对称性，因此其横向磁场分量一般不严格为零。只有在密绕近似、连续面电流模型或周期平均意义下，横向分量才可以认为相互抵消。
对于无限长密绕螺线管，端部效应可以忽略，内部磁场近似均匀，并且大小为 $\mu_0 n I$。对于有限长螺线管，内部磁场只是近似均匀，两端附近会出现明显的端部效应。
当螺距 $h$ 足够小时，螺旋线模型趋近于理想连续螺线管模型。若进一步考虑导线粗细、有限螺距和真实绕制方式，磁场会出现微小修正，但在通常的电磁学教学模型中可以忽略。
结论
通过将通电螺线管视为等距螺旋线，并应用毕奥-萨伐尔定律进行积分，我们推导出了轴线上磁场轴向分量的解析表达式：
$$
\boxed{B_z(z) = \frac{\mu_0 I}{2h}\left(\frac{z}{\sqrt{a^2+z^2}} - \frac{z-L}{\sqrt{a^2+(z-L)^2}}\right)}.
$$
其中 $h$ 为螺距，$n=1/h$ 为单位长度匝数。
在无限长密绕螺线管的理想极限下，该结果化为
$$
\boxed{\mathbf{B} = \mu_0 n I\hat{\mathbf{z}}}.
$$
这一推导过程完整展现了从螺旋电流元到宏观螺线管磁场公式的数学联系，同时也说明了理想螺线管公式成立所依赖的物理近似。
