卫星磁力矩器姿态控制模拟器

参数设置

卫星质量

惯性矩 I_sat

kg·m²

轨道高度

轨道倾斜角

51.6°=ISS、98°=太阳同步

线圈圈数 N

线圈面积 A

cm²

线圈电流 I

目标角速度 ω_t

°/s

姿态变更要达到的角速度

计算结果

—

轨道磁场 B (μT)

—

磁矩 m (A·m²)

—

最大力矩 τ (μN·m)

—

角加速度 α (°/s²)

—

目标到达时间 (s)

—

消耗功率 (W)

—

卫星·线圈·地磁场向量

卫星本体（中央）的磁线圈产生磁矩 m，与地球磁场向量 B（蓝箭头）相互作用产生力矩 τ=m×B（红箭头）。轨道用绿椭圆显示。

最大力矩 τ 对轨道高度变化

CubeSat 级别角加速度对比

理论·主要公式

$$\vec\tau = \vec m \times \vec B,\quad m = N\,I\,A,\quad \alpha = \tau/I_{sat}$$

磁力矩器基本公式。m=磁矩 [A·m²]，N=线圈圈数，I=电流 [A]，A=线圈面积 [m²]，B=地磁场向量 [T]，τ=力矩 [N·m]，I_sat=卫星惯性矩 [kg·m²]，α=角加速度 [rad/s²]。

$$B(h) \approx B_0 \left(\frac{R_E}{R_E+h}\right)^{3}\!\!\bigl[1+0.3\sin i\bigr]$$

轨道高度 h、倾斜 i 下的偶极子近似地磁场强度。R_E=6371 km，B_0≈25 μT。详细分析用 IGRF 模型。

卫星磁力矩器姿态控制 — B-dot、IGRF、CubeSat

🙋

我第一次听说"磁力矩器"，卫星能用磁石改变姿态吗？太空没有地面，怎么会产生力矩呢？

🎓

你问得好。磁力矩器是"在卫星内部制造电磁铁，与地球磁场相互作用"的装置。在 LEO（高度几百 km）仍有 25～50 μT 的地磁气存在，卫星内线圈通电后就能从 τ=m×B 产生真实力矩。不需要地面。原理和指南针的针指向北一样，只是把它用来"转动整个卫星"。因为没有活动部件，寿命长，体积小、功耗低，特别适合 CubeSat。

🙋

我用默认值算过，最大力矩是 3.94 μN·m... 这超级小啊。有意义吗？

🎓

确实"微牛顿·米"级别感觉像零。但太空没有摩擦，小力积分起来也有效。现在的默认设定（1U CubeSat、惯性 0.06 kg·m²）到达 0.5°/s 目标需要 133 秒。地面上难以相信，但轨道周期 94 分钟的卫星这速度够快了。真正的问题是"初始 detumble"，火箭分离后卫星可能每秒转 5～10°，磁力用几十分钟到几小时才能停止。但没有其他办法，几乎所有 CubeSat 都配备这个。

🙋

刚才提到"B-dot 控制"，这是什么？姿态估计听起来很复杂...

🎓

B-dot 是 detumble 用最简单的控制律，代码就 5 行。就是磁场传感器测地磁 B(t)，求时间导数 dB/dt，反号后通电。公式是 m = −k·dB/dt。卫星转得越快，B(t) 变化越大，就反向产生磁矩抵消。角动量转移到地球磁场，卫星缓慢衰减。完全不需要姿态估计器（卡尔曼滤波之类），只要磁传感器 1 个和磁力矩器，发射直后黑暗的"暂时模式"最强大。

🙋

那磁力矩器"单独"能完美控制姿态吗？我听过反应轮，为什么需要两个？

🎓

好问题！磁力矩器单独其实"欠驱动（underactuated）"。从 τ=m×B 外积，平行于 B 的轴无法产生力矩。所以瞬间无法全方向控制。轨道周期内 B 向改变，可以平均控制 3 轴，但很慢，精密指向困难。实际做法是"精密指向=反应轮，反应轮饱和时=磁力矩器卸载动量"。Bradford Space MTQ-200、CubeSpace CubeTorquer、ISIS iMTQ 是商用品，Planet Labs Doves 卫星也用这个组合。

🙋

最近 CubeSat 业界有什么新动向吗？

🎓

有意思的是"发射器内 detumble"。传统做法是卫星从火箭分离后才开始 B-dot，现在新的发射器（Yaesu 系或小卫星舱）提前在分离前就磁力整定卫星。还有 IGRF（International Geomagnetic Reference Field）更新让轨道磁场预测精度提高，磁力控制律也用"B 模型+陀螺+磁传感器"混合估计改进精度。基础理论 60 年前就有，但这领域今天仍在活跃发展。

常见问题

线圈型磁力矩器用 m = N·I·A 求得。N 是圈数，I 是电流（A），A 是线圈 1 圈的有效面积（m²）。1U CubeSat 用的典型值为 N=400 圈，A=80 cm²=8e-3 m²，I=50 mA 时 m≈0.16 A·m²。本工具用 cm² 和 mA 输入，m 换算为 A·m²，然后从 τ=m×B 算出最大力矩。

地球磁场由内部偶极子磁石近似，B∝1/r³，轨道高度越高磁场越弱。LEO（高度 200～1000 km）约 20～50 μT，地表约 50 μT，500 km 处约 25 μT。高纬度极地通过时磁场向量较强，低纬度（赤道）较弱，可用轨道倾斜角的 sin 修正。详细用 IGRF（International Geomagnetic Reference Field）模型。

不能。从 τ=m×B 外积式，平行于磁场向量 B 的轴无法产生力矩，瞬间处于欠驱动（underactuated）状态。纯磁力控制在轨道周期内平均可控制 3 轴，但反应时间为分～小时量级。精密指向需与反应轮组合：反应轮进行精密控制，磁力矩器进行动量卸载。

B-dot 控制是 detumble（发射后阻止不规则旋转）用的最简单控制律。磁矩设为 m = −k·dB/dt，与地磁场时间导数成正比，将角动量转移到地球磁场，卫星缓慢衰减。不需要姿态估计器，只用磁场传感器 1 个和磁力矩器，火箭分离直后也能可靠减速，几乎所有 CubeSat 都装有此"安全模式"。

实际应用

CubeSat 商用卫星群（Planet Labs Doves 等）：Planet Labs 运营数百颗 3U Dove 卫星，每颗都用磁力矩器+反应轮组合维持地球指向。发射后 detumble 用 B-dot，成像时精密指向用反应轮，长期运营中动量卸载再用磁力矩器。线圈规格数百圈×数十 cm²×50～100 mA 是典型，与本工具默认值一致。

科学观测任务（ICESat-2、CYGNSS）：NASA 冰盖高度测量卫星 ICESat-2 和飓风观测卫星群 CYGNSS（8 颗）都用磁力矩器进行动量卸载。小型化必需且要求精密指向的科学卫星标配反应轮+磁力矩器。极轨道或 50° 左右倾斜轨道运行，磁场向量定期大幅摆动，磁力矩器效率高。

商用 ADCS 单元（Bradford Space、ISIS、CubeSpace）：1U～12U CubeSat 市场有 Bradford Space MTQ-200、CubeSpace CubeTorquer、ISIS iMTQ、Honeywell ACS-Series 等商品。磁矩 0.1～1.0 A·m²，质量 30～200 g，功耗 0.1～2 W 是典型。作为 COTS（商用现成品）购买组装，新创卫星公司很依赖。

CAE 中磁环境仿真：磁力矩器设计通常用 IGRF 模型生成轨道上 B(t)，与 6 自由度刚体方程相结合，在 MATLAB/Simulink 或 Basilisk 等宇航仿真器验证。本工具的闭式近似"线圈规格是否现实"初步评估，详细轨道平均响应和非线性耦合用数值仿真。

常见误解和注意事项

常见误解首先是"磁力矩器能完全精密 3 轴指向"。从 τ=m×B 外积，B 平行轴无力矩，瞬间欠驱动。轨道周期内 B 转向可平均 3 轴动，反应时间分～小时级。秒级精密指向（地观摄像、激光通信、太空望远镜）必须配反应轮/CMG。单磁力宣称精密指向的设计需警惕。

其次是"地磁场 B 处处一定"错觉。本工具用 B(h)=B₀·(R_E/(R_E+h))³·[1+0.3·sin i] 近似，实际经度、纬度、时间（太阳活动、磁暴）大幅变动。南大西洋异常（SAA）磁场局部削弱，磁力矩器效能减半。磁暴时 B 短时强化产生意外力矩。准确设计需 IGRF-13/14 并轨道计算生成 B(t) 时间序列。

最后是"只看线圈电力 P=R·I²"误区。本工具仅计电阻损耗，实际磁力驱动电路有 PWM 驱动、滤波、过流保护，0.1～0.3 W 额外损耗。磁力矩器磁场还与卫星其他器件（磁传感器、霍尔传感器、反应轮 HDD）干扰，需 PWM 管理（典型 50% on/off）。设计要分"平均功率"和"峰值磁矩"评估。

卫星磁力矩器姿态控制 — B-dot、IGRF、CubeSat

常见问题

实际应用

常见误解和注意事项

使用指南

具体计算例

实务注意事项

卫星 磁力矩器姿态控制模拟器

卫星 磁力矩器姿态控制 — B-dot、IGRF、CubeSat

常见问题

实际应用

常见误解和注意事项

使用指南

具体计算例

实务注意事项

卫星磁力矩器姿态控制模拟器

卫星磁力矩器姿态控制 — B-dot、IGRF、CubeSat