3. 刚体动力学基础：转动惯量、欧拉方程与姿态动力学建模

各位同学，欢迎来到《卫星姿态控制仿真环境搭建指南》的第三讲。今天咱们聊聊刚体动力学。说实话，这部分内容是整个姿态控制仿真的“地基”。地基没打牢，后面盖的楼再漂亮也得塌。我当年刚入行时，就因为转动惯量算错了一个符号，导致仿真出来的卫星在天上翻跟头，那叫一个尴尬。

好，废话不多说，咱们直接进入正题。

3.1 转动惯量：卫星的“惯性脾气”

先问大家一个问题：为什么推一辆空购物车很轻松，推一辆装满货物的购物车就很费劲？

答案很简单：质量越大，惯性越大。但卫星在太空中，除了平动惯性，还有转动惯性。转动惯量就是描述物体“抗拒”转动的物理量。说白了，就是卫星的“转动脾气”——你给它一个力矩，它转得快还是慢，全看这个量。

对于刚体卫星，转动惯量是一个3×3的矩阵，我们叫它惯量张量：

I = [[Ixx, Ixy, Ixz],
     [Iyx, Iyy, Iyz],
     [Izx, Izy, Izz]]

其中，Ixx、Iyy、Izz是主转动惯量，而Ixy、Ixz这些是惯性积。如果卫星的质量分布对称，惯性积就为零。嗯，这里要注意：实际卫星几乎不可能完全对称，所以惯性积不能随便忽略。

3.2 欧拉动力学方程：卫星的“运动法则”

有了转动惯量，接下来就要描述卫星怎么动了。这就引出了欧拉动力学方程。你想想看，牛顿第二定律告诉我们F = ma，那转动版本就是M = I·α（力矩等于转动惯量乘角加速度）。但事情没那么简单——因为卫星在三维空间里转，坐标轴也在跟着转。

欧拉动力学方程的矢量形式如下：

M = I·ω̇ + ω × (I·ω)

其中：

• M：作用在卫星上的外力矩（比如喷气推力、飞轮反作用力矩）
• I：惯量张量
• ω：卫星角速度
• ω̇：角加速度
• ω × (I·ω)：陀螺力矩项，这是三维转动的“特色产物”

为什么会多出这一项？因为卫星在转动时，角动量方向会变化，产生一个“虚拟”的力矩。我在项目中遇到过，有人忘了加这一项，结果仿真出来的卫星姿态完全不符合物理规律——它自己就莫名其妙地加速了。记住：陀螺力矩不是可选项，是必选项。

3.3 角动量定理：守恒才是王道

欧拉方程是从力矩角度看的。换个角度，从角动量来看，事情会更清晰。

角动量H = I·ω。角动量定理说：外力矩等于角动量的变化率。

M = dH/dt

如果外力矩为零（比如卫星在太空中没有喷气、没有飞轮作用），那么角动量守恒：

H = 常数

这个性质非常有用。举个例子：卫星上装一个反作用飞轮，当飞轮加速时，卫星本体就会向相反方向转动——这就是姿态控制的底层原理。角动量在卫星本体和飞轮之间转移，但总角动量保持不变。

我个人习惯在仿真中，每一步都检查一下总角动量是否守恒。如果发现它变了，那一定是代码有bug。这个方法帮我抓出过不少低级错误。

3.4 姿态动力学建模：把理论变成代码

好了，理论讲完了，咱们看看怎么把这些东西写成仿真代码。下面是一个简化版的姿态动力学仿真框架：

import numpy as np

class RigidBodyDynamics:
    def __init__(self, inertia_matrix):
        # 惯量张量
        self.I = inertia_matrix
        # 惯量张量的逆（用于计算角加速度）
        self.I_inv = np.linalg.inv(self.I)
        # 角速度（初始为零）
        self.omega = np.array([0.0, 0.0, 0.0])
        
    def update(self, torque, dt):
        """
        根据欧拉动力学方程更新角速度
        torque: 外力矩（3维向量）
        dt: 时间步长
        """
        # 计算角加速度：ω̇ = I⁻¹ · (M - ω × (I·ω))
        omega_dot = self.I_inv @ (torque - np.cross(self.omega, self.I @ self.omega))
        # 欧拉积分更新角速度
        self.omega += omega_dot * dt
        return self.omega

这段代码虽然简单，但核心逻辑都在了。实际工程中，我们还会加入姿态表示（四元数或欧拉角）、干扰力矩模型等。但万变不离其宗——欧拉方程是骨架。

下面这张图总结了本章的知识体系，方便大家对照理解：

好了，这一章的内容就到这里。记住：刚体动力学是姿态控制仿真的基石，转动惯量、欧拉方程、角动量定理这三者缺一不可。下次咱们会在此基础上，引入姿态表示方法——四元数和欧拉角，到时候就能真正让卫星“动起来”了。

---