惯性矩权威指南：惯性、矩与计算全解析

理解物体如何转动，是机器人动力学、机械设计乃至物理学的核心课题。惯性矩正是解开这一谜题的关键参数。下面我们从“惯性”和“矩”这两个最基础的概念入手，逐步推导到惯性矩在机器人实际工程中的计算与应用。

惯性：物体抗拒状态改变的“刚性”

惯性是物体与生俱来的一种“惰性”——它倾向于维持当前的运动状态，无论静止还是匀速直线运动。这种抵抗状态变化的属性是物质固有的，与是否受力或运动与否无关。

牛顿第一定律精准描述了这种特性：任何物体都保持匀速直线运动或静止状态，除非外力迫使它改变这一状态。

需要明确几个关键点：

• 惯性大小唯一取决于质量：质量越大，惯性越大，启动或停止所需的力就越大。
• 惯性不是力：我们常说“物体因惯性继续运动”，但不能说“物体受到惯性力的作用”。惯性是属性，而非作用力。
• 惯性与速度无关：高速飞驰的汽车，其惯性依然仅由质量决定。高速带来的是更大的动量和动能，但并未改变它抵抗状态改变的固有属性。

矩：描述分布的“加权”工具

“矩”源自拉丁语，最初用于力学中描述使物体旋转的效应。随后演化为一个通用的数学概念，核心在于量化物理量的分布特征。

其逻辑很简单：以某一参考点或轴为基准，用“距离”作为权重来衡量物理量（如质量）的分布。更本质地说，矩是「距离的某次幂 × 物理量」的积分或求和。

根据距离的幂次不同，矩分为不同阶次，各自描述不同的分布特性：

• 零阶矩：距离的0次幂（即1）乘以物理量。就是物理量的总和，比如物体的总质量。
• 一阶矩：距离的1次幂乘以物理量。描述物理量分布的“位置中心”，如质心就是质量的一阶矩除以总质量。
• 二阶矩：距离的2次幂乘以物理量。描述物理量相对于参考点或轴的分散程度。而惯性矩正是“质量的二阶矩”。

听起来抽象？我们通过具体例子来感受。

零阶矩：就是总和

假设两个小球A和B粘在一起构成物体。小球A位于x=2，质量4kg；小球B位于x=3，质量8kg。

该物体的零阶矩（总质量）为：

直接相加，印证了零阶矩就是物理量的总和。

一阶矩：寻找“重心”

仍以上述两球为例。它们的一阶矩计算如下：

一阶矩 = (2 × 4) + (3 × 8) = 32 kg·m

距离在这里充当了“加权”角色：越远的点，其质量对一阶矩的贡献越大。这就像跷跷板，重物离支点越远，转动效果越明显。

如何用一阶矩找出质心（重心）？

想象将两球放在无质量的跷跷板两端，支点放哪里才能平衡？

• 两球质量相等时，支点自然在中间。
• 质量不等时，支点需向质量大的一侧移动。

这个支点位置就是系统质心的坐标。计算公式：质心位置 = 一阶矩 / 总质量。

代入例子：质心位置 = 32 / 12 = 2.67 m。即支点应放在距A球0.67m、距B球0.33m处才能平衡。

二阶矩：衡量“离散”程度

现在看二阶矩。同样两球，二阶矩计算如下：

二阶矩 = (2² × 4) + (3² × 8) = (4×4) + (9×8) = 16 + 72 = 88 kg·m²

注意距离变成了平方，权重被急剧放大。越远的点对二阶矩的贡献呈平方级增长。

因此，二阶矩描述物理量相对于参考点的离散或分散程度。想象一列等距的点绕圆心旋转：

离圆心越远的点，分布越“散”，二阶矩值越大。

惯性矩：转动惯性的度量

铺垫至此，核心概念呼之欲出：惯性矩就是质量的二阶矩（距离² × 质量）。

它的核心逻辑是用“距离平方”作为权重，刻画质量相对于参考轴的分散程度。该数值越大，物体绕该轴旋转时，其惯性（抵抗转动状态改变的能力）就越大。

一个直观例子：两个质量相同的长方体，绕不同轴旋转时惯性矩不同。

图中右侧惯性矩大于左侧，因为质量离旋转轴的平均距离b大于a。距离平方的加权效应使右侧“转动惯性”更大。

惯性矩是刚体转动分析（如机器人、机械结构设计）的核心参数。根据对象不同，分为两类：

• 多零件组合的惯性矩：离散质点系求和。
• 单一零件的惯性矩：连续刚体积分计算。

多零件组合的惯性矩计算

以一个简单机器人关节为例：

关节D末端连接三个零件A、B、C。已知：

• 零件A：质量0.2kg，质心到旋转轴距离0.5m
• 零件B：质量0.1kg，质心到旋转轴距离0.4m
• 零件C：质量0.3kg，质心到旋转轴距离0.2m

求组合体绕关节D旋转轴的惯性矩。

公式即各零件质量乘其质心到转轴距离的平方，再求和：

其中I是惯性矩（kg·m²），m_i是第i个零件质量，r_i是其质心到参考轴的垂直距离（注意是质心距离，不是零件边缘）。

计算如下：

I = (0.2 × 0.5²) + (0.1 × 0.4²) + (0.3 × 0.2²) = 0.05 + 0.016 + 0.012 = 0.078 kg·m²

单一零件（连续刚体）的惯性矩计算

对于单个零件，通常视为密度均匀的刚体。计算惯性矩需用微积分：将刚体看作无数个微小质量元dm的集合。

刚体绕某轴的惯性矩，等于每个质量元的（距离平方×质量）在整个体积上的积分：

其中r是微元dm到参考轴的垂直距离。对于均质刚体，dm = ρ·dV，ρ为密度，dV为微元体积。

幸运的是，常见几何形状的公式前人已推导好，直接套用即可。

常见几何体的惯性矩公式

1. 均质细杆

质量M，长度L（半径忽略）。

• 转轴过质心且垂直于杆：

I = (1/12) M L²

• 转轴过杆一端且垂直于杆：

I = (1/3) M L²

可见，转轴移至端点时，质量分布离轴更“散”，惯性矩从1/12增大到1/3。

2. 均质圆柱

质量M，半径R，高度L。

• 转轴沿圆柱中心轴（长度方向）：

I = (1/2) M R²

质量主要集中在中心轴附近，分散程度小，惯性矩也小。

• 转轴过质心且垂直于圆柱中心轴：

I = (1/12) M (3R² + L²)

质量沿长度方向分布更广，分散程度大，惯性矩更大。

3. 均质球体

质量M，半径R。

• 转轴过球心：

I = (2/5) M R²

4. 均质正方体

质量M，边长a。

• 转轴过质心且平行于棱：

I = (1/6) M a²

• 转轴过棱边且平行于棱：

I = (2/3) M a²

同样，转轴外移导致惯性矩增大。

5. 均质长方体

质量M，长、宽、高分别为a, b, c。

转轴过质心且平行于长度(a)方向：

I = (1/12) M (b² + c²)

这个公式在机器人领域应用极广。我们常将一个不规则零件绕x、y、z轴的惯性矩近似等效为均质长方体的惯性矩，然后根据此数据反推等效长方体的尺寸并可视化。

这样做的好处是：通过观察等效长方体的形状（长宽高比例），可直观判断原始零件的惯性矩数据是否合理。例如，下图展示了宇树H1机器人部分零件的惯性矩可视化结果（蓝色长方体）：

蓝色长方体的宽高比基本与其代表的零件外形比例一致，说明该零件的惯性矩参数设置较为合理。

鉴于均质长方体惯性矩公式如此重要，同时也为了加深对惯性矩计算本质的理解，我们接下来以均质正方体为例，详细走一遍其惯性矩公式的积分推导过程。

均质正方体惯性矩公式推导

推导最常见情况：均质正方体，绕通过质心且平行于棱的轴旋转时的惯性矩公式。

根据定义，需计算积分：I = ∫ r² dm。其中r是质量元dm到转轴的垂直距离。

原理与步骤分析

将正方体视为无数个微小质量元dm的集合。建立坐标系：原点在正方体质心，x、y、z轴分别平行于三条棱。假设转轴为x轴。

空间中任一点(x, y, z)处的质量元dm到x轴的垂直距离r，即该点在y-z平面上的投影到原点的距离：r² = y² + z²。

整个推导分两步：

第一步，计算一个垂直于x轴的薄片（截面）的惯性矩。 在该薄片内对y和z进行二重积分。

第二步，沿x轴方向对所有薄片积分。 完成在整个正方体体积上的三重积分。

公式推导过程

已知正方体边长a，总质量M，密度ρ常数。积分区域：x, y, z 均从 -a/2 到 a/2。

1. 建立积分式：
   惯性矩 I = ∫ (y² + z²) dm。
   对于均质体，dm = ρ dV = ρ dx dy dz。
   所以 I = ∭ (y² + z²) ρ dx dy dz。

2. 代入密度关系：
   总质量 M = ρ * V = ρ * a³，故 ρ = M / a³。
   代入：I = (M / a³) ∭ (y² + z²) dx dy dz。

3. 分离变量积分：
   被积函数(y²+z²)与x无关，三重积分拆分为：
   I = (M / a³) [ ∫_-a/2^a/2 dx ] · [ ∬ (y² + z²) dy dz ]
   第一个积分即x方向长度a。
   所以 I = (M / a²) ∬ (y² + z²) dy dz，y和z积分区间均为[-a/2, a/2]。

4. 计算二重积分：
   ∬ (y² + z²) dy dz = ∬ y² dy dz + ∬ z² dy dz。
   由于积分区域对称，y和z独立，拆分为：
   = [ ∫_-a/2^a/2 y² dy ] · [ ∫_-a/2^a/2 dz ] + [ ∫_-a/2^a/2 z² dz ] · [ ∫_-a/2^a/2 dy ]。
   计算各一维积分：
   ∫_-a/2^a/2 y² dy = [ y³/3 ]_-a/2^a/2 = (a³/24) - (-a³/24) = a³/12。
   ∫_-a/2^a/2 dz = a。
   同理 ∫_-a/2^a/2 z² dz = a³/12，∫_-a/2^a/2 dy = a。
   故二重积分结果 = (a³/12 * a) + (a³/12 * a) = a⁴/12 + a⁴/12 = a⁴/6。

5. 得到最终结果：
   代入：I = (M / a²) * (a⁴ / 6) = (M a²) / 6。
   即 I = (1/6) M a²。

推导完毕。此过程清晰展示了如何从微积分定义得到常用公式。

作为练习，读者可用类似思路推导均质长方体（质量M，边长a, b, c）绕其过质心且平行于边长a的轴的惯性矩公式，结果为：

I = (1/12) M (b² + c²)

均质长方体惯性矩在机器人中的应用

在机器人建模（如URDF文件）中，常需处理零件的惯性矩参数，并将其可视化以校验合理性。

半边长表示法

为方便，机器人领域常用“半边长”表示法：将长方体尺寸定义为2a, 2b, 2c。这里的a, b, c即从质心到各面的距离。

其优势在于质心到面的距离一目了然（a, b, c），且可简化从惯性矩反推尺寸的公式。此时，绕x轴（平行于2a边）的惯性矩公式变为：

I_xx = (1/3) M (b² + c²)

绕y轴和z轴的公式依此类推。

根据惯性矩逆推尺寸与可视化校验

在机器人URDF文件中，惯性矩参数通常如下表示：

name="base">
  
    rpy="0.0 0.0 0.0" xyz="0.0 0.0 0.0"/>
    value="0.01"/>
    ixx="0.0001" ixy="0.0" ixz="0.0" iyy="0.0002" iyz="0.0" izz="0.0003"/>
  

其中ixx, iyy, izz即绕x, y, z轴的惯性矩。可将其等效为一个均质长方体，反推出半边长a, b, c，从而实现可视化：

通过对比可视化长方体的比例与实际零件的大致外形，可快速判断惯性矩数据是否合理。具体步骤如下：

1. 物理合理性校验：
   - 基础校验：质量和所有惯性矩必须为非负数。
   - 三角不等式校验：刚体的惯性矩必须满足“任意两轴惯性矩之和大于等于第三轴惯性矩”，即 I_xx + I_yy ≥ I_zz，以及另外两个类似关系。这是刚体的物理约束。

2. 构建并求解方程组：
   根据半边长方体的惯性矩公式，得方程组：
   I_xx = (1/3) M (b² + c²)
   I_yy = (1/3) M (a² + c²)
   I_zz = (1/3) M (a² + b²)
   解方程组，得半边长a, b, c：
   a = √[ (3/(2M)) (I_yy + I_zz - I_xx) ]
   b = √[ (3/(2M)) (I_xx + I_zz - I_yy) ]
   c = √[ (3/(2M)) (I_xx + I_yy - I_zz) ]

总结

我们从“惯性”和“矩”的底层概念出发，逐步深入到“惯性矩”的本质——它是质量的二阶矩，用距离平方加权衡量质量绕轴的分散程度，从而决定物体的转动惯性。

惯性只由质量决定，而惯性矩同时受质量和质量分布（距离）的双重影响。我们探讨了离散系统和连续刚体的计算方法，并推导了关键的正方体惯性矩公式。

最后，我们看到了这一理论在机器人领域的落地：通过将复杂零件的惯性矩等效为长方体并可视化，来校验物理参数的合理性。理解这个过程，不仅有助于读懂模型文件，更能加深对物体转动特性的本质认识。

来源：互联网