2、卫星姿态控制基础：姿态描述、姿态动力学方程、姿态控制执行机构概述

卫星姿态控制是动量轮卸载策略设计与工程实现的物理基础。在讨论卸载算法、推力器点火逻辑或磁力矩器卸载律之前，必须首先建立三个核心认知：姿态如何数学描述（运动学）、姿态如何受力矩作用而演化（动力学）、以及执行机构如何产生力矩并影响角动量（控制接口）。本章将围绕这三条主线展开，为后续章节中动量轮饱和机理与卸载策略分析奠定严格的理论与工程框架。

图2-1展示了姿态控制系统的完整逻辑链路。从工程视角看，姿态描述解决“卫星现在朝向哪里”的问题，姿态动力学解决“卫星为什么以及如何改变朝向”的问题，而执行机构则是“用什么手段改变朝向”的物理载体。动量轮卸载策略本质上是在执行机构层面，通过外部力矩（磁力矩、推力）来管理动量轮内部存储的角动量，使其不超出轮子转速的物理极限。因此，理解本章内容是掌握卸载策略的前提。

2.1 姿态描述：从欧拉角到四元数

2.1.1 核心原理：姿态的数学表征

卫星姿态描述的核心是建立本体坐标系（Body Frame，记为B）相对于参考坐标系（如地心惯性系ECI、轨道系LVLH）的旋转关系。工程中常用的描述方法有三种：

• 欧拉角（Euler Angles）：通过三次连续旋转定义姿态。常见转序有3-1-3（偏航-滚动-偏航）和3-2-1（偏航-俯仰-滚动，即Yaw-Pitch-Roll）。欧拉角直观但存在万向节锁（Gimbal Lock）问题，且三角函数计算存在奇异性。
• 四元数（Quaternion）：定义为 q = [q0, q1, q2, q3]，其中 q0 为标量部分，满足归一化约束 q0² + q1² + q2² + q3² = 1。四元数无奇异、计算效率高，是星载计算机的标准选择。
• 姿态矩阵（DCM）：3×3正交矩阵，直接用于坐标变换，但参数冗余（9个元素），不适合实时递推。

在动量轮卸载策略中，四元数误差 δq 是控制律计算的关键输入。例如，当星敏测得当前姿态四元数 q_meas 与目标姿态 q_ref 存在偏差时，误差四元数计算为：

δq = q_ref⁻¹ ⊗ q_meas
其中 ⊗ 表示四元数乘法，q_ref⁻¹ 为共轭四元数。

2.1.2 判断方法：如何选择姿态描述方式

工程选型时需根据任务阶段和计算资源权衡：

2.1.3 案例拆解：某遥感卫星姿态确定中的四元数处理

某低轨遥感卫星（轨道高度650km）使用星敏感器+陀螺组合定姿。星敏输出四元数频率5Hz，陀螺输出角速度频率100Hz。在姿态确定滤波中，采用扩展卡尔曼滤波（EKF）融合两者数据。初始对准时，地面注入目标四元数 q_ref = [1, 0, 0, 0]（即本体与轨道系对齐）。星敏实测 q_meas = [0.9998, 0.012, -0.008, 0.005]，则误差四元数计算为：

δq = q_ref⁻¹ ⊗ q_meas = [0.9998, 0.012, -0.008, 0.005]（因q_ref为单位四元数）
对应欧拉角（3-2-1转序）：滚动≈1.37°，俯仰≈-0.92°，偏航≈0.57°

该误差值作为控制律输入，驱动动量轮产生恢复力矩。若误差超过预设阈值（如滚动>0.5°），则触发姿态机动模式，动量轮输出饱和力矩。此案例说明，姿态描述的精度直接影响卸载策略的触发时机判断。

2.1.4 常见误区与总结

总结：姿态描述是控制系统的“眼睛”，四元数因其数学优越性成为星载主流选择。但需注意归一化、转序一致性以及与欧拉角的转换校验。

2.2 姿态动力学方程：动量轮耦合与干扰力矩

2.2.1 核心原理：欧拉方程的扩展形式

带动量轮的卫星姿态动力学方程可写为：

I·ω̇ + ω × (I·ω + h_w) = T_ext + T_act
其中：
- I：卫星主惯量矩阵（3×3）
- ω：本体角速度矢量（rad/s）
- h_w：动量轮总角动量矢量（N·m·s）
- T_ext：外部干扰力矩（重力梯度、气动、太阳光压、剩磁等）
- T_act：执行机构产生的控制力矩（动量轮自身力矩或磁力矩等）

该方程揭示了动量轮卸载的物理本质：当外部干扰力矩长期存在时，动量轮角动量 h_w 会持续累积。例如，重力梯度力矩在低轨卫星上通常为10⁻⁵~10⁻⁴ N·m量级，若持续作用数小时，动量轮转速将逼近其饱和极限（通常为3000-6000 rpm）。卸载策略的核心就是通过外部力矩（磁力矩或推力）来“吸收”这部分累积角动量，使 h_w 回到安全范围。

2.2.2 判断方法：干扰力矩分析与卸载需求评估

工程上，判断是否需要卸载以及何时卸载，需进行以下分析：

1. 干扰力矩预算：列出各干扰源量级，计算总角动量累积速率 ΔH = ∫T_ext dt。
2. 动量轮容量裕度：定义安全阈值（如额定角动量的80%），超过则触发卸载。
3. 卸载能力匹配：评估磁力矩器或推力器能提供的最大卸载力矩，确保在有限时间内完成卸载。

2.2.3 案例拆解：某对地观测卫星动量轮饱和事件复盘

背景：某太阳同步轨道卫星（轨道高度700km，倾角98°）搭载三台反作用轮（偏航轮+俯仰轮+滚动轮），额定角动量4 N·m·s，最大转速6000 rpm。入轨后第3天，俯仰轮转速异常升高至5800 rpm，接近饱和。

初始现象：遥测显示俯仰轮转速持续上升，且姿态控制精度下降（俯仰角误差从0.05°恶化至0.2°）。地面分析发现，卫星处于当地正午附近，太阳光压与重力梯度力矩叠加，且卫星质心偏差导致气动力矩增大。

关键参数：

俯仰轴干扰力矩：T_d = 2.3×10⁻⁵ N·m（重力梯度1.2e-5 + 气动0.8e-5 + 光压0.3e-5）
累积速率：ΔH = T_d × 轨道周期(5800s) ≈ 0.133 N·m·s/圈
预计饱和时间：从3000rpm（2 N·m·s）到5800rpm（3.87 N·m·s）需约14圈

分析判断：地面团队确认动量轮即将饱和，需立即启动卸载。由于卫星配备磁力矩器（最大磁矩100 A·m²），且当地磁场强度约3×10⁻⁵ T，最大卸载力矩约3×10⁻³ N·m，远大于干扰力矩，具备卸载条件。

操作步骤：

1. 切换至安全模式，暂停成像任务。
2. 计算卸载律：根据当前磁场方向与动量轮角动量方向，生成磁力矩指令 M = -k × (B × h_w)。
3. 执行卸载：持续30分钟，俯仰轮转速从5800rpm降至3200rpm。
4. 恢复姿态控制，重新进入对地指向模式。

结果对比：卸载后俯仰角误差恢复至0.05°以内，动量轮转速裕度恢复至45%。

风险与代价：卸载期间卫星姿态扰动增大（约0.3°），且磁力矩器功耗增加（约15W）。复盘发现，若提前2圈启动卸载，可避免进入饱和预警状态。

2.2.4 常见误区与总结

总结：姿态动力学方程是卸载策略的“物理引擎”。理解干扰力矩的累积效应、动量轮的角动量交换机制，是设计有效卸载律的前提。

2.3 姿态控制执行机构概述：动量轮、磁力矩器与推力器

2.3.1 核心原理：执行机构的力矩产生机制

卫星姿态控制执行机构主要分为三类，其力矩产生原理和卸载适用性差异显著：

• 动量轮/反作用轮（RW）：通过电机加速/减速飞轮，产生反作用力矩。特点是控制精度高、响应快，但角动量存储有限（典型2-20 N·m·s）。卸载时，轮子本身不能产生净力矩，必须借助外部力矩。
• 磁力矩器（MTQ）：利用地磁场与线圈磁矩的相互作用产生力矩 T = M × B。特点是质量轻、功耗低、无限寿命，但力矩方向受限于当地磁场方向，且力矩量级小（10⁻⁴~10⁻³ N·m）。
• 推力器（推进系统）：通过喷气产生推力，形成力矩。特点是力矩大、方向可控，但消耗推进剂，寿命受限。通常用于大角度机动或紧急卸载。

2.3.2 判断方法：卸载执行机构选型依据

2.3.3 案例拆解：磁力矩器卸载的工程约束分析

场景：某低轨卫星（轨道高度500km，倾角97.6°）采用三轴动量轮+三轴磁力矩器方案。磁力矩器最大磁矩为80 A·m²，地磁场强度约2.5×10⁻⁵ T，则最大卸载力矩为：

T_max = M_max × B = 80 × 2.5e-5 = 2×10⁻³ N·m

该力矩仅能产生约0.002 N·m的卸载力矩，而动量轮累积速率约0.1 N·m·s/圈。因此，完成一次卸载（降低0.5 N·m·s）需要约250秒的持续卸载时间。工程上，通常采用“分段卸载”策略：在每个轨道周期的特定弧段（如磁场方向有利时）开启磁力矩器，避免长时间连续工作导致功耗过高。

关键约束：磁力矩器卸载存在不可控自由度——力矩方向始终垂直于磁场方向，因此无法独立控制三轴角动量。工程上通过交叉耦合控制律（如B-dot或LQR）实现近似解耦。

2.3.4 常见误区与总结

总结：执行机构是卸载策略的“手脚”。动量轮负责精细控制，磁力矩器负责日常卸载，推力器负责紧急情况。选型时需综合考虑力矩量级、资源消耗、寿命和安全性。