4. 姿态动力学基础
姿态动力学回答的是“力矩如何改变航天器的姿态运动”。在姿态控制中，我们最终要控制的量是姿态角或姿态角速度，而真正被力矩直接改变的是角动量。因此，理解刚体角动量、惯量矩阵与欧拉动力学方程，是设计控制律、分析耦合和判断机动可行性的基础。
一、刚体角动量：从质点到连续体
对绕固定点转动的刚体，角动量定义为各质量微元位置与动量叉积的体积分：
H = ∫ r × (ω × r) dm其中 r 是质量微元相对参考点的位置矢量，ω 是刚体角速度。展开三重叉积后可得：
H = I · ω这里 I 是 3×3 惯量矩阵，ω 是角速度矢量。这个式子说明：角动量方向一般不与角速度方向重合，除非 ω 恰好沿惯量主轴。这是刚体动力学与质点动力学最本质的区别之一。
二、惯量矩阵与主惯量
惯量矩阵在体坐标系下写为：
I = [ Ixx -Ixy -Ixz
-Ixy Iyy -Iyz
-Ixz -Iyz Izz ]对角元是绕各轴的转动惯量，非对角元是惯性积，反映质量分布的不对称性。对任意刚体，总可以找到一组正交坐标轴，使惯性积为零，这组轴称为主轴，对应的对角元称为主惯量 Ix、Iy、Iz。工程上常把体坐标系取在主轴附近，使惯量矩阵近似对角，从而简化动力学方程。
以某中型遥感卫星为例（示例参数，非实测）：Ix=1200 kg·m²，Iy=1500 kg·m²，Iz=900 kg·m²。由于三个主惯量互不相等，该星属于三轴惯量不等体，机动时会出现明显的陀螺耦合。
三、欧拉动力学方程
在体坐标系下，对角化惯量矩阵后，欧拉动力学方程为：
Ix·ω̇x - (Iy - Iz)·ωy·ωz = Mx
Iy·ω̇y - (Iz - Ix)·ωz·ωx = My
Iz·ω̇z - (Ix - Iy)·ωx·ωy = Mz式中 Mx、My、Mz 是外力矩在体轴上的分量。方程左边第一项是角加速度产生的惯性力矩，第二项是陀螺力矩。当三个主惯量相等时，陀螺项为零，三轴运动解耦；当主惯量不等时，绕某一轴的转动会通过陀螺项耦合到另外两轴。
关键判断：若 Iy≈Iz，则绕 x 轴的机动几乎不产生 y、z 轴耦合；若三者差异大，则单轴机动会激发明显的横向角速度，需要三轴协同控制。
四、陀螺力矩与耦合特性
陀螺力矩的本质是角动量方向在惯性空间中变化时产生的效应。当航天器同时存在两个方向的角速度时，角动量矢量会“扫过”一个锥面，从而产生与角速度乘积成正比的力矩。以示例卫星为例，若 ωx=0.5 °/s、ωy=0.3 °/s，则绕 z 轴的陀螺力矩约为：
Mz_gyro = (Ix - Iy)·ωx·ωy
= (1200 - 1500)×(0.5×π/180)×(0.3×π/180)
≈ -0.045 N·m这个量级虽小，但在高精度指向任务中足以引起姿态漂移。若同时绕三轴机动，耦合项会相互叠加，使角速度轨迹偏离单轴设计值。因此，在机动规划时通常采用“先绕惯量最大轴、再绕最小轴”的顺序，或直接使用三轴耦合的轨迹优化。
五、工程含义与机动策略
从上述关系可以提炼出三条工程结论。第一，惯量矩阵必须通过质量特性测量或CAD模型精确获得，误差会直接进入陀螺力矩计算，影响前馈补偿精度。第二，大角度快速机动时，陀螺耦合不可忽略，控制律中应包含耦合补偿项，或采用基于欧拉方程的轨迹规划。第三，当某一轴惯量远小于其他两轴时，该轴对耦合最敏感，应避免在该轴上单独施加长时间力矩。
实用建议：在姿态机动仿真中，先以单轴模型估算机动时间和力矩需求，再用三轴欧拉方程校核耦合引起的横向角速度峰值。若横向角速度超过指向精度允许的阈值，则需要调整机动顺序或增加耦合补偿。
理解欧拉动力学方程的关键，不是记住三个公式，而是建立“角动量变化—外力矩—角速度耦合”之间的因果直觉。后续的姿态控制律设计、动量管理、甚至故障下的姿态安全，都建立在这一组方程之上。