2. 航天器姿态动力学基础回顾

安全模式的所有判据、阈值和恢复动作，最终都要落到姿态动力学方程上。若对运动学与动力学的基本关系没有清晰认识，安全模式设计就容易变成“凭经验设阈值”：既说不清为什么某个角速度门限是 0.5 °/s，也说不清为什么在重力梯度力矩占主导的轨道上，某些姿态反而更稳定。本节把与安全模式最相关的动力学基础重新梳理一遍，重点不是推导，而是建立“方程—参数—环境力矩—安全判据”之间的对应关系。

一、姿态描述：三种方法各管一段

描述航天器姿态有三套常用工具，它们在安全模式中承担不同角色，不能互相替代。

方向余弦矩阵（DCM）是最基础的表示。设本体系 b 相对参考系 n 的姿态矩阵为 C_bn，则任意矢量满足 v^b = C_bn vⁿ。它的优点是物理意义直接、可无奇异地进行矢量变换，缺点是 9 个元素存在 6 个约束，直接用于滤波或积分会累积正交性误差。安全模式中，DCM 主要用于把太阳方向、地磁矢量、气动来流方向从轨道系或惯性系变换到本体系，供敏感器和执行机构使用。

欧拉角最直观，工程上常用 3-1-2 或 3-2-1 转序。以 3-2-1 为例，偏航 ψ、俯仰 θ、滚动 φ 对应本体系相对参考系的三次转动。它的致命问题是万向节锁定：当俯仰角 θ 接近 ±90° 时，ψ 与 φ 耦合，无法唯一确定。安全模式若采用欧拉角判据，必须显式处理这一奇异区，否则在接近奇异姿态时会出现“角度正常但控制发散”的假安全。

四元数用 4 个参数表示姿态，无奇异，适合姿态误差计算和闭环控制。设误差四元数 q_e = q_c^-1 ⊗ q_b，则标量部分 q_e0 = cos(Φ/2)，矢量部分给出旋转轴。安全模式中，四元数最适合做“姿态偏差是否超过门限”的判据，例如要求 q_e0 ≥ cos(Φ_max/2)。需要注意四元数的双覆盖：q 与 -q 表示同一姿态，误差计算时必须先取符号一致，否则会得到 360° 的错误偏差。

二、运动学与动力学：两套方程的分工

姿态运动学只描述姿态参数如何随角速度变化，不涉及力和力矩。以四元数为例：

dq/dt = 0.5 * q ⊗ [0, ω_b]^T

其中 ω_b 是本体系相对参考系的角速度。运动学是“几何关系”，它决定了姿态偏差如何随时间累积。安全模式中，若角速度持续存在而姿态控制失效，姿态偏差会按上式积分增长，这就是“角速度门限”和“姿态偏差门限”必须联合使用的原因。

姿态动力学描述角速度如何被力矩改变。刚体转动方程在本体系下为：

J * dω_b/dt + ω_b × (J * ω_b) = T_ext + T_ctrl

其中 J 为惯量矩阵，T_ext 为环境力矩，T_ctrl 为控制力矩。展开后可以看到，即使没有外力矩，ω_b × (J ω_b) 这一陀螺项也会让角速度在三个轴之间交换能量。对于惯量不等的主轴，自由刚体的运动是绕最大惯量轴的稳定锥形运动与绕最小惯量轴的不稳定运动。安全模式设计必须知道本星惯量分布：若最小惯量轴与太阳帆板法线接近，一旦该轴出现扰动，角速度可能快速发散，此时“等待自然稳定”是不可行的。

把两套方程放在一起看：动力学决定角速度的演化，运动学决定姿态的演化。安全模式的控制逻辑通常是“先阻尼角速度，再恢复姿态”，因为角速度是姿态偏差的源头，也是陀螺、动量轮等执行机构最直接的负载。

三、环境力矩：安全模式中最容易被低估的扰动源

安全模式往往发生在控制能力下降之后，此时环境力矩的相对重要性显著上升。四类主要环境力矩对安全模式的影响各不相同。

重力梯度力矩在轨道系下可写为 T_gg = 3μ/r³ · (n̂ × (J n̂))，其中 n̂ 为地心指向航天器的单位矢量。它与轨道半径三次方成反比，在 500 km 轨道上量级约为 10^-5～10^-4 N·m。重力梯度有一个重要特性：它会让最小惯量轴趋向于指向地心，最大惯量轴趋向于轨道法线。这意味着如果安全模式把航天器置于“自由漂移”状态，姿态不会完全随机，而是会向重力梯度稳定姿态缓慢收敛。设计时应判断：这个自然收敛姿态是否满足太阳帆板对日、天线对地等安全约束。若不满足，就不能依赖自然稳定。

气动力矩与大气密度、来流速度和迎风面积有关，在 300～600 km 低轨是主要扰动之一。其大小可近似为 T_aero ≈ 0.5 ρ V² C_d A · l，其中 l 为压心与质心的距离。安全模式中，若航天器进入翻滚，迎风面积和压心位置都会变化，气动力矩可能从阻尼变为激励。因此低轨安全模式通常要求尽快把角速度降到较低水平，避免长时间大迎风面积翻滚。

磁力矩由星上剩磁矩 m 与地磁场 B 叉乘产生：T_m = m × B。它的特点是随轨道位置周期变化，在磁赤道和极区方向不同。安全模式中，磁力矩既可能是扰动，也可以被磁力矩器主动利用来阻尼角速度。若星上剩磁较大，即使不主动控制，磁力矩也会让角速度出现周期性波动，判据设计时要避免把这种正常波动误判为故障。

太阳光压力矩与太阳光压、受照面积和压心-质心偏移有关，量级通常最小，但对大型太阳帆板或大面质比航天器不可忽略。它的方向随太阳矢量变化，在安全模式中可能造成缓慢的姿态漂移。若安全模式持续时间较长（数小时到数天），太阳光压力矩的累积效应必须纳入姿态偏差预算。

四、把基础关系落到安全模式设计上

综合来看，安全模式设计需要回答三个问题：第一，当前角速度和姿态偏差是否超过门限，这由运动学和四元数误差判据给出；第二，角速度发散的风险有多大，这由惯量分布和动力学陀螺项决定；第三，环境力矩是帮助收敛还是加剧发散，这需要按轨道高度、姿态和航天器构型分场景评估。

一个实用的做法是建立“环境力矩占优表”：按轨道高度段（如 300～500 km、500～800 km、800 km 以上）和典型姿态（对地、对日、翻滚），估算四类力矩的量级，找出主导项。例如在 400 km 低轨对地姿态下，气动力矩和重力梯度力矩通常同量级；在 800 km 以上，气动力矩迅速衰减，重力梯度和磁力矩成为主导。安全模式的阻尼策略和等待时间应据此调整，而不是全轨道段使用同一套阈值。

最后需要强调参照系一致性。运动学方程中的 ω_b 是本体系相对参考系的角速度，动力学方程中的 J 和 T 也都在本体系下定义。若敏感器输出在轨道系，必须先变换到本体系再代入方程。安全模式中常见的“角速度看着不大但姿态持续偏”的问题，往往就是参照系混用或四元数符号未统一造成的。把这些基础关系理清，后续的安全模式判据、阈值和恢复策略才有可靠的依据。