第3章：刚体动力学基础——牛顿-欧拉方程、动量与动量矩、惯性张量与质量属性
各位同学，欢迎来到刚体动力学的核心地带。
说实话，很多搞仿真的朋友，一听到「牛顿-欧拉」这几个字就头大。我当年刚入行时也一样，觉得这玩意儿太理论了，离实际工程很远。直到有一次，我负责一个机械臂的动力学建模，算出来的关节力矩总是跟实测对不上，折腾了整整两周。最后发现，问题出在惯性张量上——我把它当成标量处理了。嗯，从那以后，我再也不敢小看这些基础概念了。
3.1 牛顿-欧拉方程：刚体运动的两把钥匙
刚体动力学，说白了就是研究「力怎么让物体动起来」。这里有两套逻辑：
- 牛顿方程：管平移，描述质心的运动。
- 欧拉方程：管转动，描述绕质心的旋转。
我个人习惯把这两者分开想。先看平移：
F = m * a_c其中 F 是合外力，m 是质量，a_c 是质心加速度。这个公式大家太熟了，但要注意——这里的加速度必须是质心的加速度。我在项目中见过有人直接用末端点的加速度去算，结果全错了。
再看转动部分，欧拉方程：
M = I * α + ω × (I * ω)这里 M 是合外力矩，I 是惯性张量，α 是角加速度，ω 是角速度。你想想看，如果物体转得很快，那个 ω × (I * ω) 项就不能忽略。我做过一个高速旋转的飞轮仿真，一开始没加这一项，结果能量不守恒，飞轮越转越快——这显然不对。
核心要点：牛顿方程管平动，欧拉方程管转动。两者通过质心位置耦合在一起。写代码时，我建议把平动和转动分开求解，最后再合成。
3.2 动量与动量矩：守恒才是硬道理
动量，就是质量乘以速度：p = m * v。动量矩（也叫角动量），是惯性张量乘以角速度：L = I * ω。
为什么要提这两个量？因为它们是守恒量。在没有外力（矩）的情况下，动量和动量矩保持不变。这个性质在仿真中特别有用。
我记得有一次做卫星姿态仿真，卫星在太空中不受外力矩，按理说动量矩应该守恒。但我的仿真结果却显示动量矩在缓慢漂移。查了半天，发现是数值积分方法的问题——欧拉法不保能量。后来换了四阶龙格-库塔法，问题就解决了。
避坑指南：我曾经在仿真中直接用动量矩的变化率去算力矩，结果发现数值噪声被放大了。后来我改用「先算力矩，再积分得到动量矩」的方式，稳定性好很多。记住：微分放大噪声，积分平滑噪声。
3.3 惯性张量与质量属性：别把刚体当质点
很多初学者容易犯一个错误：把刚体当成质点来处理。质点只有质量，没有形状。但刚体有形状，有质量分布，这就引出了惯性张量。
惯性张量是一个 3×3 的矩阵：
I = | Ixx Ixy Ixz |
| Iyx Iyy Iyz |
| Izx Izy Izz |对角线元素 Ixx、Iyy、Izz 是绕各轴的转动惯量。非对角线元素 Ixy、Ixz 等是惯性积，反映了质量分布的不对称性。
我建议你记住一个原则：坐标系选得好，计算量少一半。如果坐标系选在质心，且坐标轴与刚体的惯性主轴重合，那么惯性张量就变成对角矩阵：
I = | Ixx 0 0 |
| 0 Iyy 0 |
| 0 0 Izz |这样一来，欧拉方程就简化成了三个独立的标量方程。我在做机器人动力学时，每次都先把连杆的惯性张量转换到质心坐标系下，再对角化处理。省事不少。
注意：惯性张量是相对于某个坐标系定义的。换坐标系时，惯性张量必须通过平行轴定理和旋转矩阵进行变换。我见过有人直接把 CAD 软件里导出的惯性张量拿来用，结果坐标系没对齐，仿真结果一塌糊涂。
3.4 知识体系框架
下面这张图，是我自己总结的刚体动力学知识结构。你可以把它当作一张地图，随时回来看看自己走到哪了。
3.5 实战中的几个关键点
讲到这里，我想分享几个实战中容易踩的坑：
- 坐标系一致性：所有物理量必须在同一个坐标系下表达。我习惯用全局惯性系表达平动量，用体坐标系表达转动量。
- 惯性张量的更新：刚体旋转后，惯性张量在全局坐标系下会变化。需要用旋转矩阵做变换：I_global = R * I_body * R^T。
- 数值稳定性：欧拉方程中的 ω × (I·ω) 项，在高速旋转时可能引起数值发散。我建议用四元数代替欧拉角，避免万向锁问题。
一个小技巧：如果你用 MATLAB 或 Python 做仿真，可以先用符号推导验证公式，再转成数值代码。我每次写新的动力学模型，都会先用符号算一遍，确保没有代数错误。
好了，这一章的内容就到这里。刚体动力学的基础打牢了，后面讲多体系统建模、约束处理、数值积分时，你才能游刃有余。记住：基础不牢，地动山摇。我当年吃过的亏，希望你们别再吃了。