惯性张量

贡献者： addis

本文存在未完成的内容。

预备知识 1　常见几何体的转动惯量，张量，瞬时转轴

　　¹我们已经初步学习了刚体绕固定轴的转动，现在我们做一个更一般的讨论，也就是惯性系中刚体绕某个固定点的转动，注意这个固定点未必是质心。刚体绕固定轴转动只有一个自由度 $θ$ ，而绕固定点转动却可以有 3 个自由度。为了讨论方便，本文以旋转的固定点作为坐标原点。

1. 惯性张量

　　刚体绕固定点旋转的角动量等于惯性张量 $I$ 乘以瞬时角速度矢量 $ω$

\begin{matrix} (1) & L = I ω . \end{matrix}

该式推导见下文。其中惯性张量

I

是一个 3 维方矩阵，把三维角速度矢量线性映射到三维动量矢量。惯性张量与物体的质量分布、朝向、以及固定点的位置有关。也就是说，在刚体旋转时，

I

一般也会随时间变化。

　　惯性张量的矩阵元一般记为

\begin{matrix} (2) & I = (\begin{matrix} I_{x x} & I_{x y} & I_{x z} \\ I_{y x} & I_{y y} & I_{y z} \\ I_{z x} & I_{z y} & I_{z z} \end{matrix}) . \end{matrix}

如果将

x, y, z

分别记为

x_{1}, x_{2}, x_{3}

，则

I

的矩阵元可以通过体积分计算

\begin{matrix} (3) & I_{i j} = δ_{i j} \int r^{2} ρ (r) d V - \int x_{i} x_{j} ρ (r) d V (i, j = 1, 2, 3) . \end{matrix}

其中

r

是刚体上某点的位置矢量，

r = | r |

，

ρ (r)

是位置

r

处的密度，与刚体的位置和朝向有关。也可以表示成矩阵积分的形式，即对每个矩阵元分别积分

\begin{matrix} (4) & I = \int ρ (r) (\begin{matrix} y^{2} + z^{2} & - x y & - x z \\ - y x & x^{2} + z^{2} & - y z \\ - z x & - z y & x^{2} + y^{2} \end{matrix}) d V . \end{matrix}

若把刚体看作质点系，令第

i

个质点的质量为

m_{k}

，位置为

r_{k} = (x_{1, k}, x_{2, k}, x_{3, k})

，式 3 也可记为

\begin{matrix} (5) & I_{i j} = δ_{i j} \sum_{k} m_{k} r_{k}^{2} - \sum_{k} m_{k} x_{i k} x_{j k} (i, j = 1, 2, 3) . \end{matrix}

　　从定义容易看出， $I$ 是一个对称矩阵，即 $I_{i j} = I_{j i}$ 。而对称矩阵必定有三个正交归一的本征矢，详见 “刚体的主轴”。

例 1　长方体的惯性张量

　　令质心在原点，密度为常数，质量为 $M$ 的长方体延 $x, y, z$ 方向的长度分别为 $a, b, c$ ，那么有

\begin{matrix} (6) & \begin{aligned} ∭ x y d x d y d z = ∭ y z d x d y d z = ∭ z x d x d y d z = 0, \\ ∭ x^{2} d x d y d z = \frac{1}{12} a^{3} b c, \\ ∭ y^{2} d x d y d z = \frac{1}{12} a b^{3} c, \\ ∭ z^{2} d x d y d z = \frac{1}{12} a b c^{3} . \end{aligned} \end{matrix}

使用

M = ρ a b c

可得

\begin{matrix} (7) & I = \frac{1}{12} M (\begin{matrix} b^{2} + c^{2} & 0 & 0 \\ 0 & a^{2} + c^{2} & 0 \\ 0 & 0 & a^{2} + b^{2} \end{matrix}), \end{matrix}

可以发现这个是一个对角矩阵。这就意味着，如果长方体分别绕

x, y, z

轴转动，它的角动量也分别延

x, y, z

轴。由式 12 可以验证三个方向的转动惯量是对的。当转轴与坐标轴不重合时，可见转动惯量的方向与转动方向一般不同。

　　特殊地，当 $a = b = c$ （正方体）时，惯性张量变为单位矩阵乘以常数

\begin{matrix} (8) & I = \frac{1}{6} M a^{2} (\begin{matrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{matrix}) . \end{matrix}

由于单位矩阵乘以任何矢量都得到相同的矢量，正方体延任意转轴转动时，角动量都和角速度矢量成正比

\begin{matrix} (9) & L = I ω = \frac{1}{6} M a^{2} ω, \end{matrix}

这可以验证式 13 的结果。

例 2　细杆的角动量

　　假设细杆沿 $x$ 轴放置，长度为 $L$ 。要计算它的惯性张量，我们可以从头积分（习题 1 ），也可以假设细杆是一个细长的长方体，直接用例 1 的结论（令 $a = L, b \to 0, c \to 0$ ）得

\begin{matrix} (10) & I = \frac{1}{12} M (\begin{matrix} 0 & 0 & 0 \\ 0 & L^{2} & 0 \\ 0 & 0 & L^{2} \end{matrix}) . \end{matrix}

可以看到，如果转轴延

x

方向，即

ω = (ω, 0, 0)^{T}

，角动量为零，如果延

y, z

方向，则角动量大小都等于

\begin{matrix} (11) & L = \frac{1}{12} M L^{2} ω . \end{matrix}

对应的转动惯量为

I = M L^{2} / 12

，这与式 3 一致。

习题 1　

　　直接使用积分计算细杆的角动量（式 10 ）。

推导

　　我们把刚体看作质点系刚体，刚体绕固定点转动时，刚体上的每个点都在做瞬时圆周运动。根据角动量的定义式 1 以及圆周运动的速度式 5

\begin{matrix} (12) & L = \sum_{i} m_{i} r_{i} v_{i} = \sum_{i} m_{i} r_{i} \times (ω \times r_{i}) = \sum_{i} m_{i} r_{i}^{2} ω - \sum_{i} m_{i} (ω \cdot r_{i}) r_{i}, \end{matrix}

其中最后一步使用了 “连续叉乘的化简” 中的式 1 。将上式右边写成分量的形式，并将求和表示为刚体密度

ρ (r)

的积分得

\begin{matrix} (13) & \begin{aligned} (\begin{array}{c} L_{x} \\ L_{y} \\ L_{z} \end{array}) & = \int ρ r^{2} (\begin{array}{c} ω_{x} \\ ω_{y} \\ ω_{z} \end{array}) d V - \int ρ (\begin{array}{c} x x & x y & x z \\ y x & y y & y z \\ z x & z y & z z \end{array}) (\begin{array}{c} ω_{x} \\ ω_{y} \\ ω_{z} \end{array}) d V \\ = [\int ρ (\begin{array}{c} r^{2} - x x & - x y & - x z \\ - y x & r^{2} - y y & - y z \\ - z x & - z y & r^{2} - z z \end{array}) d V] (\begin{array}{c} ω_{x} \\ ω_{y} \\ ω_{z} \end{array}) \end{aligned} . \end{matrix}

对比式 1 可知，方括号中的积分就是惯性张量

I

（式 4 ）。

2. 坐标旋转变换

预备知识 2　空间旋转矩阵，相似变换

　　一般情况下当刚体相对于选取的参考系发生平移或任意方向的转动，（相对于原点的）惯性张量都会发生改变。我们先来讨论惯性张量在两个原点相同朝向不同的坐标系之间如何变换，这叫做惯性张量的旋转变换。在计算刚体运动时，通常先建立一个与刚体始终相对静止的参考系叫做体坐标系（body frame），另外选一个原点相同的惯性系，称为实验室坐标系（lab frame）。实际问题中，刚体以及体坐标系绕坐标原点随时间相对实验室系转动。令同一个矢量在实验室系的坐标为 $r_{L}$ ，在体坐标系的坐标为 $r_{B}$ ，令旋转变换矩阵为 $R$ ，满足

\begin{matrix} (14) & r_{L} = R r_{B}, \end{matrix}

令体坐标系和实验室系中的惯性张量分别为

I_{0}

和

I

，则二者的关系为

\begin{matrix} (15) & I = R I_{0} R^{- 1} = R I_{0} R^{T}, \end{matrix}

其中使用了酉矩阵的性质

R^{- 1} = R^{T}

（式 4 ）。注意该式是一个相似变换，所以

I

和

I_{0}

互为相似矩阵。

　　证明：要计算 $L = I ω$ ，把式 15 右边三个矩阵从右到左依次作用在 $ω$ 上。先计算 $R^{T} ω$ 把实验室系中的角速度矢量 $ω$ 变换到体坐标系中的 $ω_{0}$ 。再计算 $I_{0} ω_{0}$ 得到体坐标系中的角动量 $L_{0}$ ，最后计算 $R L_{0}$ 得到实验室系中的角动量 $L$ 。证毕。

例 3　

　　如果转轴和轻杆法线的夹角为 $θ$ ，（未完成）

3. 坐标平移变换

　　若刚体原来的惯性张量（关于原点）为 $I$ ，总质量为 $M$ ，质心为 $r_{c} = (x_{1, c}, x_{2, c}, x_{3, c})$ ，把它平移 $r_{0} = (x_{1, 0}, x_{2, 0}, x_{3, 0})$ 后令新的惯性张量（关于原点）为 $I^{'}$ ，则

\begin{matrix} (16) & I_{i j}^{'} = I_{i j} + 2 δ_{i j} M r_{c} \cdot r_{0} + M (δ_{i j} r_{0}^{2} - x_{i c} x_{j 0} - x_{j c} x_{i 0} - x_{i 0} x_{j 0}) (i, j = 1, 2, 3) . \end{matrix}

特殊地，若

r_{c} = 0

，则

\begin{matrix} (17) & I_{i j}^{'} = I_{i j} + M (δ_{i j} r_{0}^{2} - x_{i 0} x_{j 0}) (i, j = 1, 2, 3) . \end{matrix}

　　推导：把式 5 中 $r_{k}^{2}$ 替换为 $(r_{k} + r_{0})^{2}$ ， $x_{i k}$ 替换为 $x_{i k} + x_{i 0}$ ， $x_{j k}$ 替换为 $x_{j k} + x_{j 0}$ ，再利用质心的定义（式 6 ） $r_{c} = \sum_{k} m_{k} r_{k} / M$ ，就可以得到式 16 。

习题 2　长方体的惯性张量 2

　　若把例 1 中的长方体平移 $r_{0} = (a / 2, b / 2, c / 2)$ ，使其一个角处于原点，求新的惯性张量。

　　答案：

\begin{matrix} (18) & \begin{aligned} I^{'} & = \frac{1}{12} M (\begin{array}{c} 4 b^{2} + 4 c^{2} & - 3 a b & - 3 a c \\ - 3 a b & 4 a^{2} + 4 c^{2} & - 3 b c \\ - 3 a c & - 3 b c & 4 a^{2} + 4 b^{2} \end{array}) . \end{aligned} \end{matrix}

习题 3　陀螺模型

　　显然，球体的惯性张量就是其转动惯量（式 10 ） $2 M R^{2} / 5$ 乘以单位矩阵。若把它沿 $z$ 轴平移 $a$ ，就得到一个陀螺模型，求陀螺的惯性张量。

　　答案：

\begin{matrix} (19) & I = \frac{2}{5} M R^{2} (\begin{matrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{matrix}) + M a^{2} (\begin{matrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 0 \end{matrix}) . \end{matrix}

未完成：引用陀螺的数值模拟

1. ^ 参考 Wikipedia 相关页面。

惯性张量

1. 惯性张量

例 1 长方体的惯性张量

例 2 细杆的角动量

习题 1

推导

2. 坐标旋转变换

例 3

3. 坐标平移变换

习题 2 长方体的惯性张量 2

习题 3 陀螺模型

例 1　长方体的惯性张量

例 2　细杆的角动量

习题 1　

例 3　

习题 2　长方体的惯性张量 2

习题 3　陀螺模型