2. 空间磁场环境与磁控物理基础

磁控执行的前提，是航天器必须“有磁场可用、有力矩可算、有方向可判”。本章围绕四个问题展开：近地轨道磁场用什么模型描述、沿轨道磁场如何变化、磁矩与磁力矩如何产生、以及工程上可用的力矩量级与方向特性。理解这四点，才能判断某次磁控任务是否可行、需要多大磁矩、在什么弧段执行最有效。

一、地磁场模型：从偶极子到高阶球谐

工程上最常用的近似是倾斜偶极子模型：把地磁场等效为一个位于地心、与自转轴约成 11.5° 夹角的磁偶极子。在轨道高度 400～800 km 范围内，偶极子模型给出的磁场强度误差通常在 10%～20%，足以支撑初步力矩估算与磁矩方向规划。若需要更精确的磁场矢量（例如用于磁强计定标或高精度姿态确定），则采用球谐模型，如 IGRF（国际地磁参考场）或 WMM（世界地磁模型），其标量磁位可写为：

V(r,θ,φ) = a · Σ(n=1..N) (a/r)^(n+1) · Σ(m=0..n)
           [ g_n^m·cos(mφ) + h_n^m·sin(mφ) ] · P_n^m(cosθ)

B = -∇V
其中：a 为地球参考半径，g_n^m、h_n^m 为高斯系数，
      P_n^m 为施密特准归一化勒让德函数，
      N 为截断阶数（IGRF-13 取 N=13）。

实际使用中，输入为 UTC 时间、地理经纬度与海拔高度，输出为北东地（NED）或地心惯性系下的磁场三分量。需要特别说明坐标系约定：本章及后续章节中，磁矩 m 与磁感应强度 B 均在同一坐标系下表达，磁力矩按 T = m × B 计算，方向由右手定则确定。若在轨道系与本体系之间切换，必须先做坐标变换再计算叉乘，否则方向会整体反号。

二、轨道磁场分布特征

低轨航天器绕地飞行时，磁场呈现明显的周期性变化。以 500 km 太阳同步轨道为例，轨道周期约 94.6 min，磁场强度在赤道附近约 2.2×10⁻⁵ T（约 22000 nT），在极区附近可升至 5.5×10⁻⁵ T 以上。磁场方向随纬度连续旋转，在轨道系下表现为近似正弦的矢量轨迹。

更关键的是“磁场变化率”。磁力矩控制本质上是利用磁场方向变化产生平均力矩，因此 dB/dt 直接决定控制效率。在 500 km 轨道上，磁场方向在约 1/4 轨道周期内可旋转近 90°，对应角速率约 0.02 rad/s 量级。这意味着磁控的带宽天然受限，通常只能用于慢速姿态调整、动量卸载和阻尼控制，而不适合快速机动。

三、磁矩与磁力矩产生机理

磁控执行机构的核心是磁矩线圈（磁力矩器）。当线圈通电流 I、匝数为 N、包围面积为 A 时，产生的磁矩大小为 |m| = N·I·A，方向由右手定则确定：四指沿电流方向弯曲，拇指指向磁矩方向。三轴正交布置的磁力矩器可合成任意方向的磁矩矢量。

磁矩与地磁场相互作用产生磁力矩：T = m × B。其大小 |T| = |m|·|B|·sinθ，方向同时垂直于 m 和 B。这意味着：当 m 与 B 平行时力矩为零，垂直时力矩最大。工程上正是利用这一特性，通过调节三轴电流使合成磁矩与当地磁场保持接近垂直，从而获得最大控制效率。

以某 500 kg 级遥感卫星为例，若配置三轴磁力矩器，单轴最大磁矩 50 A·m²，在 |B| = 3×10⁻⁵ T 的轨道位置，当 m 与 B 垂直时单轴最大力矩约为 50 × 3×10⁻⁵ = 1.5×10⁻³ N·m。三轴合成后，在最优方向下可接近 2.6×10⁻³ N·m。这个量级足以在数小时内完成飞轮动量卸载，但无法在秒级完成姿态机动。

四、可用力矩量级与方向特性

磁力矩的可用性受三个因素制约：磁矩上限、当地磁场强度、以及 m 与 B 的夹角。由于磁场方向沿轨道持续变化，即使磁矩大小恒定，力矩方向也在不断旋转。因此磁控通常采用“平均力矩”策略：在一个轨道周期内，通过周期性切换电流方向，使力矩在某一方向上的分量累积，而交变分量相互抵消。

方向特性上还需注意：磁力矩只能产生垂直于 B 的力矩分量，沿 B 方向的分量恒为零。这意味着磁控对绕磁场方向轴的姿态控制能力最弱，工程上常与飞轮或推力器配合使用。在动量卸载任务中，通常选择磁场方向变化较大的弧段（如极区附近）执行卸载，以获得更好的方向覆盖。