1. 课程导入与数值模拟概览

刚体转动是经典力学中最容易被低估、却在实际工程中无处不在的一类运动。当我们把研究对象从“质点”换成“有大小、有形状的物体”时，描述方式会发生一次本质跃迁：力不再只决定质心的加速度，力对某一点的力矩才决定物体绕轴的转动状态；而描述转动“惯性大小”的量，也从单一的质量变成了依赖质量分布的转动惯量。正是这一跃迁，让刚体转动成为连接静力学、动力学与工程结构分析的关键枢纽。

与转动相伴的核心守恒量是角动量。在没有外力矩（或外力矩之和为零）的条件下，系统的角动量矢量保持不变。这条规律解释了许多看似反直觉的现象：花样滑冰运动员收臂加速旋转、直升机尾桨为何必须存在、陀螺为何能“站住”而不倒。角动量守恒之所以重要，是因为它把复杂的内部相互作用“屏蔽”掉了——无论物体内部如何变形、如何重新分配质量，只要外力矩为零，总角动量就锁定不变。这种“只看总量、不问细节”的特性，使它成为分析碰撞、爆炸、姿态控制等问题时最有力的工具之一。

然而，真实刚体运动往往无法用解析式干净地求解。当转动惯量随姿态变化、当外力矩依赖角度或速度、当多个刚体通过约束耦合在一起时，运动方程会变成一组强非线性的微分方程。这时，数值模拟就从“可选项”变成“必需品”。数值模拟的角色可以概括为三层：第一层是求解器，把牛顿—欧拉方程或拉格朗日方程离散成时间步进格式，逐步推进状态；第二层是实验台，在计算机里反复改变参数、观察极端工况，成本远低于真实实验；第三层是验证器，用守恒量（能量、角动量）是否漂移来检验模型与算法的正确性。工程中的卫星姿态确定、机器人关节动力学、车辆侧倾稳定性、旋转机械的临界转速分析，背后都是同一套“建模—离散—积分—校验”的数值链条。

上图把本课程的知识结构拆成三层并给出递进路径。需要强调的是，三层之间不是单向流水线，而是存在反馈回路：工程应用暴露出的偏差，会反过来要求我们修正物理模型或数值格式。这正是数值模拟区别于“算一道题”的地方——它是一个可迭代、可检验的闭环。

为了让抽象概念落地，这里给出一个贯穿课程的教学示例（以下参数为教学设定，非实测数据）：考虑一个绕固定水平轴转动的均质细杆，质量 m = 2.0 kg，长度 L = 1.2 m，转轴过其一端且垂直于杆。其转动惯量为 I = mL²/3 = 2.0 × 1.44 / 3 = 0.96 kg·m²。若初始时刻杆静止于水平位置，重力矩 τ = mg(L/2)cosθ 会驱动它下摆。取 g = 9.8 m/s²，初始力矩为 2.0 × 9.8 × 0.6 = 11.76 N·m，初始角加速度 α = τ/I ≈ 12.25 rad/s²。这个简单系统没有解析初等解（涉及椭圆积分），却非常适合作为第一个数值实验对象：用四阶龙格—库塔以步长 h = 0.001 s 积分，记录角速度与机械能，观察能量相对漂移是否随步长减小而收敛。

判断模拟是否可信，不能只看曲线“像不像”，而要看守恒量。对上述无摩擦摆杆，总机械能应恒定；若能量在长时间积分后单调漂移，说明格式耗散或增益过大，应减小步长或换用更高阶/辛格式。对角动量守恒问题，则要检查 L 矢量的模与方向是否稳定。这类“用不变量做标尺”的习惯，是本课程希望建立的第一项工程素养。