刚体

本条目中，向量与标量分别用粗体与斜体显示。例如，位置向量通常用 $\mathbf{r}\,\!$ 表示；而其大小则用 $r\,\!$ 来表示。

<image>

在物理学里，理想刚体（Rigid body或Rigid object）是一种有限尺寸，可以忽略形变的固体。不论是否感受到外力，在刚体内部，质点与质点之间的距离都不会改变。这种理想模型适用条件是，运动过程比固体中的弹性波的传播要缓慢得多。根据相对论，这种物体不可能实际存在，但物体通常可以假定为完美刚体，前提是必须满足运动速度远小于光速的条件。

在经典力学里，刚体通常被视为连续质量分布体；在量子力学里，刚体被视为一群粒子的聚集。例如，分子（由假定为质点的电子与核子组成）时常会被视为刚体。

## 运动学

### 位置与取向

刚体是由一群数量超多的质点组成。实际而言，不可能精确地追踪其中每一个质点的运动。为了简化运算，可以利用刚体的“刚性”，即其内部所有质点彼此之间距离不变的性质。假若物体具有刚性，则倚靠设定三个不共线质点的位置，就足以设定此物体的位置。这意味着，在三维空间里，刚体至多只有九个自由度，但由于假定三个质点之间的距离固定不变，所以，刚体只有六个自由度。假设还有其它约束，例如，刚体的运动必需绕着其内部一点旋转（定点转动），或绕着其内部一直轴旋转（定轴转动），则自由度会小于六。

关于其它任意质点P的位置，只要知道质点P对于上述三个质点之中的任意一个质点的相对位置，就可以重建这质点的位置。通常，整个刚体的空间位形可以简易地以参数设定：

- 刚体的“位置”：挑选刚体内部一点G来代表整个刚体，通常会设定物体的质心或形心为这一点。从空间参考系S观测，点G的位置就是整个刚体在空间的位置。位置可以应用向量的概念来表示：向量的起点为参考系S的原点，终点为点G。设定刚体的位置需要三个坐标，例如，采用直角坐标系，这三个坐标为x-坐标、y-坐标、z-坐标。这用掉了三个自由度。
- 刚体的取向：描述刚体取向的方法有好几种，包括方向余弦、欧拉角、四元数等等。这些方法设定一个附体参考系B的取向（相对于空间参考系S）。附体参考系是固定于刚体的参考系。相对于刚体，附体参考系的取向固定不变。由于刚体可能会呈加速度运动，所以附体参考系可能不是惯性参考系。空间参考系是某设定惯性参考系，例如，在观测飞机的飞行运动时，附着于飞机场控制塔的参考系可以设定为空间参考系，而附着于飞机的参考系则可设定为附体参考系。刚体的取向需要用到另外三个自由度。

#### 方向余弦

<image>

方向余弦方法可以用来设定附体参考系B的取向，即刚体的取向。假设沿着参考系S的坐标轴的三个单位向量分别为$\hat{\mathbf{x}}_1$、$\hat{\mathbf{x}}_2$、$\hat{\mathbf{x}}_3$，沿着参考系B的坐标轴的三个单位向量分别为$\hat{\mathbf{e}}_1$、$\hat{\mathbf{e}}_2$、$\hat{\mathbf{e}}_3$。定义$\hat{\mathbf{e}}_i$与$\hat{\mathbf{x}}_j$之间的方向余弦$a_{ij}$为

$a_{ij}\ \stackrel{def}{=}\ \cos{(\theta_{ij})} =\hat{\mathbf{e}}_i\cdot\hat{\mathbf{x}}_j$。

其中，$\theta_{ij}$是$\hat{\mathbf{e}}_i$与$\hat{\mathbf{x}}_j$之间的夹角。

$\hat{\mathbf{e}}_1$、$\hat{\mathbf{e}}_2$、$\hat{\mathbf{e}}_3$与$\hat{\mathbf{x}}_1$、$\hat{\mathbf{x}}_2$、$\hat{\mathbf{x}}_3$之间的关系分别为

$\hat{\mathbf{e}}_1=\cos{(\theta_{11})}\hat{\mathbf{x}}_1+\cos{(\theta_{12})}\hat{\mathbf{x}}_2+\cos{(\theta_{13})}\hat{\mathbf{x}}_3=a_{11}\hat{\mathbf{x}}_1+a_{12}\hat{\mathbf{x}}_2+a_{13}\hat{\mathbf{x}}_3$。

$\hat{\mathbf{e}}_2=\cos{(\theta_{21})}\hat{\mathbf{x}}_1+\cos{(\theta_{22})}\hat{\mathbf{x}}_2+\cos{(\theta_{23})}\hat{\mathbf{x}}_3=a_{21}\hat{\mathbf{x}}_1+a_{22}\hat{\mathbf{x}}_2+a_{23}\hat{\mathbf{x}}_3$。

$\hat{\mathbf{e}}_3=\cos{(\theta_{31})}\hat{\mathbf{x}}_1+\cos{(\theta_{32})}\hat{\mathbf{x}}_2+\cos{(\theta_{33})}\hat{\mathbf{x}}_3=a_{31}\hat{\mathbf{x}}_1+a_{32}\hat{\mathbf{x}}_2+a_{33}\hat{\mathbf{x}}_3$。

两个参考系的坐标轴所形成的矩阵称为“方向余弦矩阵”$A$：

$A= \begin{bmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{bmatrix}$。

采用爱因斯坦求和约定，由于$\hat{\mathbf{e}}_i=a_{ij}\hat{\mathbf{x}}_j$，给定方向余弦矩阵$A$，则可设定附体参考系B的取向，也就是刚体的取向。

反过来，经过一番运算，可以得到$\hat{\mathbf{x}}_j=a_{ij}\hat{\mathbf{e}}_i$。

给定位置向量$\mathbf{r}$

$\mathbf{r}=x_1\hat{\mathbf{x}}_1+x_2\hat{\mathbf{x}}_2+x_3\hat{\mathbf{x}}_3=e_1\hat{\mathbf{e}}_1+e_2\hat{\mathbf{e}}_2+e_3\hat{\mathbf{e}}_3$。

则$\mathbf{r}$与$\hat{\mathbf{e}}_i$的内积为

$\mathbf{r} \cdot \hat{\mathbf{e}}_i=e_i =a_{i1}x_1+a_{i2}x_2+a_{i3}x_3=a_{ij}x_j$。

方向余弦矩阵$A$可以将从空间参考系S观测的位置坐标$(x_1,x_2,x_3)$，变换为从附体参考系B观测的位置坐标$(e_1,e_2,e_3)$，因此又称为“变换矩阵”。

变换矩阵$A$也可以做反变换如下：

$x_j=a_{ij}e_i$。

变换矩阵$A$是一种正交矩阵，满足“正交条件”

$a_{ij}a_{ik}=\delta_{jk}$。

其中，$\delta_{jk}$是克罗内克函数。

注意到$\theta_{ij}$与$\theta_{ji}$不同，夹角$\theta_{ji}$是$\hat{\mathbf{e}}_j$与空间参考系S的坐标轴单位向量$\hat{\mathbf{x}}_i$之间的夹角。变换矩阵$A$通常不是对称矩阵。

<image>

对于二维旋转，变换矩阵$A$可以视为旋转矩阵。例如，将附体参考系B或刚体旋转，从$\hat{\mathbf{e}}_1$、$\hat{\mathbf{e}}_2$、$\hat{\mathbf{e}}_3$旋转$\theta$角弧成为$\hat{\mathbf{e}}'_1$、$\hat{\mathbf{e}}'_2$、$\hat{\mathbf{e}}'_3$；其中，$\hat{\mathbf{e}}_3=\hat{\mathbf{e}}'_3$。对于这旋转，旋转矩阵$A$为

$A= \begin{bmatrix} \cos{\theta} & \sin{\theta} & 0 \\ -\sin{\theta} & \cos{\theta} & 0 \\ 0 & 0 & 1 \end{bmatrix}$。

参考轴$\hat{\mathbf{e}}'_i$与$\hat{\mathbf{e}}_j$之间的关系为

$\hat{\mathbf{e}}'_i=a_{ij}\hat{\mathbf{e}}_j$。

旋转矩阵$A$也可以视为将点P的位置向量$\mathbf{r}=x_i\hat{\mathbf{x}}_i$旋转$-\theta$角弧成为点P'的位置向量$\mathbf{r}'=x'_i\hat{\mathbf{x}}_i$：

$x'_i=a_{ij}x_j$。

#### 欧拉角

<image>

方向余弦矩阵$A$足以设定附体参考系B的取向。但是，矩阵$A$有九个元素，而刚体只能供给三个自由度来设定取向，因为这九个元素不是自变量。欧拉角的三个自变量可以用来设定刚体的取向。

相对于空间参考系S，附体参考系B的取向，可以用三个欧拉角来设定。参阅右图。设定xyz-轴为空间参考系S的坐标轴，XYZ-轴为附体参考系B的坐标轴。称xy-平面与XY-平面的相交为“交点线”，用英文字母（N）代表。按照“zxz顺规”，欧拉角可以这样定义：

- $\alpha$是x-轴与交点线（N）之间的夹角。
- $\beta$是z-轴与Z-轴之间的夹角。
- $\gamma$是交点线（N）与X-轴之间的夹角。

每一个欧拉角的旋转都对应于一个简单的旋转矩阵：

$A_{\alpha}= \begin{bmatrix} \cos \alpha & \sin \alpha & 0 \\ -\sin \alpha & \cos \alpha & 0 \\ 0 & 0 & 1 \end{bmatrix}$。

$A_{\beta}=\begin{bmatrix} 1 & 0 & 0 \\ 0 & \cos \beta & \sin \beta \\ 0 & -\sin \beta & \cos \beta \end{bmatrix}$。

$A_{\gamma}= \begin{bmatrix} \cos \gamma & \sin \gamma & 0 \\ -\sin \gamma & \cos \gamma & 0 \\ 0 & 0 & 1 \end{bmatrix}$。

设定刚体取向的旋转矩阵$A$是由三个简单旋转矩阵$A_{\alpha}$、$A_{\beta}$、$A_{\gamma}$共同合成：

$A=A_{\gamma}A_{\beta}A_{\alpha}$。

单独分开工作，每个矩阵各自代表一种旋转。按照顺序相乘。

- 最里面的（最右的）矩阵代表绕着z轴的旋转。
- 最外面的（最左的）矩阵代表绕着Z轴的旋转。
- 在中间的矩阵代表绕着交点线的旋转。

经过一番运算，可以得到$A$矩阵：

$A= \begin{bmatrix} \cos\alpha\cos\gamma-\cos\beta\sin\alpha\sin\gamma & \sin\alpha\cos\gamma+\cos\beta\cos\alpha\sin\gamma & \sin\beta\sin\gamma \\-\cos\alpha\sin\gamma-\cos\beta\sin\alpha\cos\gamma & -\sin\alpha\sin\gamma+\cos\beta\cos\alpha\cos\gamma & \sin\beta\cos\gamma \\ \sin\beta\sin\alpha & -\sin\beta\cos\alpha & \cos\beta \end{bmatrix}$。

$A$的逆矩阵是：

$A^{-1}= \begin{bmatrix} \cos\alpha\cos\gamma-\cos\beta\sin\alpha\sin\gamma & -\cos\alpha\sin\gamma-\cos\beta\sin\alpha\cos\gamma & \sin\beta\sin\alpha \\ \sin\alpha\cos\gamma+\cos\beta\cos\alpha\sin\gamma & -\sin\alpha\sin\gamma+\cos\beta\cos\alpha\cos\gamma & -\sin\beta\cos\alpha \\ \sin\beta\sin\gamma & \sin\beta\cos\gamma & \cos\beta \end{bmatrix}$。

#### 欧拉旋转定理

欧拉旋转定理表明，在三维空间里，假设约束刚体内部一点固定不动，则其任意位移等价于绕着某固定轴的一个旋转，而这固定轴必包含这固定点。换句话说，设定附体参考系B的原点为这固定点，则附体参考系B不会因为这位移而涉及任何平移运动，再设定附体参考系B的z-轴与固定轴同轴，则这位移对应于绕着附体参考系B的z-轴旋转$\gamma$角弧，而z-轴的方向是由$\alpha$与$\beta$角弧给出。

对于内部有一点被约束固定不动的刚体（或原点固定不动的参考系），欧拉旋转定理将其任意位移等价为绕着某固定轴的一个旋转。这允许使用旋转来表达取向的改变。因此，变换矩阵$A$可以视为三维旋转的旋转矩阵，将附体参考系B或刚体做任意环绕着固定点的旋转，从$\hat{\mathbf{e}}_1$、$\hat{\mathbf{e}}_2$、$\hat{\mathbf{e}}_3$旋转成为$\hat{\mathbf{e}}'_1$、$\hat{\mathbf{e}}'_2$、$\hat{\mathbf{e}}'_3$。参考轴$\hat{\mathbf{e}}'_i$与$\hat{\mathbf{e}}_j$之间的关系为

$\hat{\mathbf{e}}'_i=a_{ij}\hat{\mathbf{e}}_j$。

#### 沙勒定理

<image>

<image>

当刚体移动时，它的位置与取向都可能会随着时间演进而改变。沙勒定理是欧拉旋转定理的一个推论。根据沙勒定理，刚体的最广义位移等价于一个平移加上一个旋转。因此，刚体运动可分为平移运动与旋转运动。刚体的现在位置与现在取向可以视为是从某个初始位置与初始取向经过平移与旋转而成。

如右图所示，从时间$t_1$到时间$t_2$，当刚体在做平移运动时，任意内部两点，点P与点Q的轨迹（以黑色实线表示）相互平行，线段$\overline{PQ}$（以黑色虚线表示）的方向保持恒定。

挑选刚体内部一点G来代表整个刚体，设定附体参考系B的原点于点G（称为“基点”），则从空间参考系S观测，在刚体内部任意一点P的位置$\mathbf{r}_P$为

$\mathbf{r}_P=\mathbf{r}_{G}+\mathbf{r}_{P/G}$。

其中，$\mathbf{r}_{G}$、$\mathbf{r}_{P/G}$分别是基点G的位置、点P对于基点G的相对位置。

从附体参考系B观测，刚体内部每一点的位置都固定不变，但从空间参考系S观测，刚体从时间$t_1$到时间$t_2$的运动，可以分为基点G从$\mathbf{r}_{G}(t_1)$到$\mathbf{r}_{G}(t_2)$的平移运动，与位移$\mathbf{r}_{P/G}$从时间$t_1$到时间$t_2$的旋转运动。

### 平移速度与角速度

从不同的参考系观测刚体运动，可能会获得不同的平移速度和不同的角速度。为了确保测量结果具有实际物理意义，必需先给定参考系。

刚体的平移速度是向量，是其位置向量的时间变化率，是附着于刚体的基点G的速度。对于纯平移运动（没有任何旋转运动），刚体内部所有点的移动速度相同。假设涉及旋转运动，则通常刚体内部任意两点的瞬时速度不相等；只有当它们恰巧处于同一直轴，而这直轴平行于转动瞬轴，则它们的瞬时速度相等。

角速度也是向量，描述刚体取向改变的角速率，和刚体旋转时的瞬时转轴的方向（欧拉旋转定理保证瞬时转轴的存在）。在任意时间，刚体内部每一个质点的角速度相同。

#### 向量的时间变化率

假设一刚体呈纯旋转运动，其附体参考系B也会跟着旋转，因此，对于任意向量$\mathbf{F}$，从这附体参考系B与从空间参考系S观测，会得到不同的结果。假设附体参考系B $(\hat{\mathbf{e}}_1,\hat{\mathbf{e}}_2,\hat{\mathbf{e}}_3)$与空间参考系S $(\hat{\mathbf{x}}_1,\hat{\mathbf{x}}_2,\hat{\mathbf{x}}_3)$同原点。对于这些参考系，三维含时向量$\mathbf{F}(t)$分解为

$\mathbf{F}=f_i\hat{\mathbf{x}}_i=F_i\hat{\mathbf{e}}_i$。

$\mathbf{F}(t)$对于时间的导数为

$\frac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}=\frac{\mathrm{d}f_i}{\mathrm{d}t}\hat{\mathbf{x}}_i=\frac{\mathrm{d}F_i}{\mathrm{d}t}\hat{\mathbf{e}}_i+F_i\frac{\mathrm{d}\hat{\mathbf{e}}_i}{\mathrm{d}t}$。

单独计算附体参考轴对于时间的导数：

$\frac{\mathrm{d}\hat{\mathbf{e}}_i}{\mathrm{d}t} =\frac{\mathrm{d}}{\mathrm{d}t}(a_{ij}\hat{\mathbf{x}}_j) =\dot{a}_{ij}\hat{\mathbf{x}}_j =\dot{a}_{ij}a_{kj}\hat{\mathbf{e}}_k$。

其中，$\dot{a}_{ij}\ \stackrel{def}{=}\ \frac{\mathrm{d}a_{ij}}{\mathrm{d}t}$是方向余弦对于时间的导数。

由于$\frac{\mathrm{d}\hat{\mathbf{e}}_i}{\mathrm{d}t}$垂直于$\hat{\mathbf{e}}_i$，$\frac{\mathrm{d}\hat{\mathbf{e}}_i}{\mathrm{d}t}$只能是其他两个单位向量的线性组合：

$\frac{\mathrm{d}}{\mathrm{d}t}\hat{\mathbf{e}}_i=-\epsilon_{ij\ell}\omega_{ij}\hat{\mathbf{e}}_{\ell}$。

其中，$\epsilon_{ij\ell}$是列维-奇维塔符号，$\omega_{ij}$是系数。

对于任意$\hat{\mathbf{e}}_{m}\ne \hat{\mathbf{e}}_i$ ,

单位向量$\hat{\mathbf{e}}_i$与$\hat{\mathbf{e}}_{m}$的内积对于时间的导数为

$$\begin{align}\frac{\mathrm{d}}{\mathrm{d}t}(\hat{\mathbf{e}}_{i}\cdot\hat{\mathbf{e}}_{m}) & = \left(\frac{\mathrm{d}}{\mathrm{d}t}\hat{\mathbf{e}}_{i}\right) \cdot\hat{\mathbf{e}}_{m} +\hat{\mathbf{e}}_{i}\cdot\left(\frac{\mathrm{d}}{\mathrm{d}t}\hat{\mathbf{e}}_{m}\right) \\ <pre> & =-\epsilon_{ij\ell}\omega_{ij}\hat{\mathbf{e}}_{\ell} \cdot\hat{\mathbf{e}}_{m} </pre> -\hat{\mathbf{e}}_{i} \cdot\epsilon_{mj\ell}\omega_{mj}\hat{\mathbf{e}}_{\ell} \\ <pre> & =-\epsilon_{ijm}\omega_{ij}-\epsilon_{mji}\omega_{mj} \\ & =-\epsilon_{ijm}(\omega_{ij}-\omega_{mj}) \\ & =0 \\ </pre> \end{align}$$

所以，$\omega_{ij}$的下标$i$多余无用，可以删除，变为$\omega_{j}$：

$\frac{\mathrm{d}\hat{\mathbf{e}}_i}{\mathrm{d}t}=-\epsilon_{ij\ell}\omega_{j}\hat{\mathbf{e}}_{\ell}$。

思考$\frac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}$方程式最右边项目$F_i\frac{\mathrm{d}\hat{\mathbf{e}}_i}{\mathrm{d}t}$，对换傀标$i$ $\ell$，可以得到

$F_i\frac{\mathrm{d}\hat{\mathbf{e}}_i}{\mathrm{d}t} =-\epsilon_{ij\ell}F_i\omega_{j}\hat{\mathbf{e}}_{\ell} =\epsilon_{ij\ell}\hat{\mathbf{e}}_{i}\omega_{j}F_{\ell} =\boldsymbol{\omega}\times\mathbf{F}$。

向量$\boldsymbol{\omega}$是由三个系数$\omega_1$、$\omega_2$、$\omega_3$组成，对应于附体参考系的三个参考轴$\hat{\mathbf{e}}_1$、$\hat{\mathbf{e}}_2$、$\hat{\mathbf{e}}_3$，系数数值可以从欧拉角计算求得：

$\omega_1=\dot{a}_{2j}a_{3j} =\dot{\alpha}\sin{\beta}\sin{\gamma}+\dot{\beta}\cos{\gamma}$。

$\omega_2=-\dot{a}_{1j}a_{3j} =\dot{\alpha}\sin{\beta}\cos{\gamma}+\dot{\beta}\sin{\gamma}$。

$\omega_3=\dot{a}_{1j}a_{2j}=\dot{\alpha}\cos{\beta}+\dot{\gamma}$。

试想对应于欧拉角$\alpha$、$\beta$、$\gamma$的三个旋转轴分别为$\hat{\mathbf{z}}$、$\hat{\mathbf{N}}$、$\hat{\mathbf{Z}}$，三个角速度分别为

$\boldsymbol{\omega}_{\alpha}=\dot{\alpha}\hat{\mathbf{z}} =\dot{\alpha}\sin{\beta}\sin{\gamma}\hat{\mathbf{X}}+\dot{\alpha}\sin{\beta}\cos{\gamma}\hat{\mathbf{Y}}+\dot{\alpha}\cos{\beta}\hat{\mathbf{Z}}$。

$\boldsymbol{\omega}_{\beta}=\dot{\beta}\hat{\mathbf{N}} =\dot{\beta}\cos{\gamma}\hat{\mathbf{X}}-\dot{\beta}\sin{\gamma}\hat{\mathbf{Y}}$。

$\boldsymbol{\omega}_{\gamma}=\dot{\gamma}\hat{\mathbf{Z}}$。

这三个角速度的向量和，对于附体参考系B的分量分别为

$\omega_X=\dot{\alpha}\sin{\beta}\sin{\gamma}+\dot{\beta}\cos{\gamma}=\omega_1$。

$\omega_Y=\dot{\alpha}\sin{\beta}\cos{\gamma}+\dot{\beta}\sin{\gamma}=\omega_2$。

$\omega_Z=\dot{\alpha}\cos{\beta}+\dot{\gamma}=\omega_3$。

注意到附体参考系B的$\hat{\mathbf{e}}_1$、$\hat{\mathbf{e}}_2$、$\hat{\mathbf{e}}_3$就是欧拉角的$\hat{\mathbf{X}}$、$\hat{\mathbf{Y}}$、$\hat{\mathbf{Z}}$，所以，向量$\boldsymbol{\omega}$是附体参考系B旋转的角速度。

总结，向量$\mathbf{F}(t)$对于时间的导数为

$\frac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}=\frac{\mathrm{d}f_i}{\mathrm{d}t}\hat{\mathbf{x}}_i=\frac{\mathrm{d}F_i}{\mathrm{d}t}\hat{\mathbf{e}}_i+\boldsymbol{\omega}\times\mathbf{F}$。

设定$\left(\frac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}\right)_{\mathrm{space}}$、$\left(\frac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}\right)_{\mathrm{body}}$分别为从空间参考系S、附体参考系B观测到的向量$\mathbf{F}(t)$对于时间的导数，上述方程式可以表达为

$\left(\frac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}\right)_{\mathrm{space}}=\left(\frac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}\right)_{\mathrm{body}}+\boldsymbol{\omega}\times\mathbf{F}$。

项目 $\boldsymbol{\omega}\times\mathbf{F}$可以想像为，从空间参考系S观测，刚体内部位置向量为$\mathbf{F}$的质点，由于刚体旋转而产生的速度。

向量$\mathbf{F}(t)$是任意向量，因此可以将$\left(\frac{\mathrm{d}}{\mathrm{d}t}\right)_{\mathrm{space}}$、$\left(\frac{\mathrm{d}}{\mathrm{d}t}\right)_{\mathrm{body}}$当作算符，这样，对应的算符方程式的形式为：

$\left(\frac{\mathrm{d}}{\mathrm{d}t}\right)_{\mathrm{space}}=\left(\frac{\mathrm{d}}{\mathrm{d}t}\right)_{\mathrm{body}}+\boldsymbol{\omega}\times$。

这算符方程式可以作用于任意含时向量。

### 运动学方程式

根据沙勒定理，刚体的最广义位移等价于一个平移加上一个旋转。挑选刚体内部一点G来代表整个刚体，设定附体参考系B的原点于基点G，则从空间参考系S观测，在刚体内部任意一点P的位置$\mathbf{r}_P$为

$\mathbf{r}_P=\mathbf{r}_{G}+\mathbf{r}_{P/G}$。

其中，$\mathbf{r}_{G}$、$\mathbf{r}_{P/G}$分别是基点G的位置、点P对于基点G的相对位置。

点P的速度$\mathbf{v}_P=\left(\frac{\mathrm{d}\mathbf{r}_P}{\mathrm{d}t}\right)_{\mathrm{space}}$为

$\mathbf{v}_P=\mathbf{v}_{G}+\mathbf{v}_{P/G}$。

其中，$\mathbf{v}_{G}=\left(\frac{\mathrm{d}\mathbf{r}_G}{\mathrm{d}t}\right)_{\mathrm{space}}$、$\mathbf{v}_{P/G}=\left(\frac{\mathrm{d}\mathbf{r}_{P/G}}{\mathrm{d}t}\right)_{\mathrm{space}}$分别是基点G的速度、点P对于基点G的相对速度。

应用前段推导出的适用于任意含时向量的算符方程式，可以计算出$\mathbf{v}_{P/G}$。由于从附体参考系B观测，刚体内部每一点的位置都固定不变，项目$\left(\frac{\mathrm{d}\mathbf{r}_{P/G}}{\mathrm{d}t}\right)_{\mathrm{body}}$等于零：

$\mathbf{v}_{P/G}=\left(\frac{\mathrm{d}\mathbf{r}_{P/G}}{\mathrm{d}t}\right)_{\mathrm{body}}+\boldsymbol{\omega}\times\mathbf{r}_{P/G}=\boldsymbol{\omega}\times\mathbf{r}_{P/G}$。

其中，$\boldsymbol{\omega}$是刚体的角速度向量。

所以，点P的速度为

$\mathbf{v}_P=\mathbf{v}_{G}+\boldsymbol{\omega}\times\mathbf{r}_{P/G}$。

点P的加速度$\mathbf{a}_P=\left(\frac{\mathrm{d}\mathbf{v}_P}{\mathrm{d}t}\right)_{\mathrm{space}}$为

$\mathbf{a}_P=\mathbf{a}_{G}+\mathbf{a}_{P/G}$。

其中，$\mathbf{a}_{G}=\left(\frac{\mathrm{d}\mathbf{v}_G}{\mathrm{d}t}\right)_{\mathrm{space}}$、$\mathbf{a}_{P/G}=\left(\frac{\mathrm{d}\mathbf{v}_{P/G}}{\mathrm{d}t}\right)_{\mathrm{space}}$分别是基点G的速度、点P对于基点G的相对速度。

再应用前段推导出的算符方程式，可以计算出

$\mathbf{a}_{P/G}=\left(\frac{\mathrm{d}\boldsymbol{\omega}}{\mathrm{d}t}\right)_{\mathrm{space}}\times\mathbf{r}_{P/G}+\boldsymbol{\omega}\times\mathbf{v}_{P/G} =\boldsymbol{\alpha}\times\mathbf{r}_{P/G}+\boldsymbol{\omega}\times(\boldsymbol{\omega}\times\mathbf{r}_{P/G})$。

其中，$\boldsymbol{\alpha}=\left(\frac{\mathrm{d}\boldsymbol{\omega}}{\mathrm{d}t}\right)_{\mathrm{space}}$是附体参考系B旋转的角加速度向量。

## 动力学

主项目：刚体动力学

当描述刚体的动力运动时，必需先处理妥善刚体的平移运动，即先选择一个参考点来代表刚体与其附体参考系B。刚体内部任意一点都可以被选为参考点（附体参考系B的原点）。但是，根据实际应用需要，比较适当的选择为：

- 刚体的质心：对于自由移动于空间的刚体，以其质心为参考点的方法通常会给出最简单的运动。
- 平移运动为零或可以简易研算的点：例如，在轴或铰链上、在万向接头的中心等等。

当质心被选为参考点时：

- 刚体的动量与其旋转运动无关。在任何时间，动量等于刚体的总质量乘以平移速度。
- 不论刚体的平移运动为何，对于质心的角动量皆等同。所以，在计算角动量时，可以忽略平移运动。在任何时间，角动量等于惯性张量乘以角速度。假若知道刚体绕主轴的角速度，那么，角动量对于每一主轴的分量，是对应的主惯性矩乘以对应的角速度；力矩是转动惯量乘以角加速度。
- 在无外力作用下，可能形成的运动为等速直线运动、稳定绕定轴转动、零力矩进动等等。
- 作用于刚体的净外力，等于总质量乘以刚体平移运动的加速度（也就是说，不论净外力矩是否为零，或这刚体是否在作旋转运动，牛顿第二运动定律可以正确地应用于刚体平移运动）。
- 总动能是平移动能与旋转动能的总和。

## 刚体的转动

- 刚体的转动定律：$M=I\alpha$，其中$M$为刚体所受合外力的力矩，$I$为刚体转动惯量，$\alpha$为刚体角加速度。
- 刚体的转动动能定理：$W=\int_{\theta_1}^{\theta_2} M\, d\theta=\frac{1}{2}I\omega_2^2-\frac{1}{2}I\omega_1^2$，其中$\int_{\theta_1}^{\theta_2} M\, d\theta$表示合外力的力矩$M$在角位移$\theta_2-\theta_1$上所作的功，$I$为刚体的转动惯量，$\omega$为刚体角速度。
- 刚体的转动和平动可以合成为刚体的平面运动，由柯尼希定理，其动能为$E_k=\frac{1}{2}mV_c^2+\frac{1}{2}I\omega^2$，其中$V_c$为刚体质心对参考系的速度。

## 参阅

- 无穷小旋转
- 欧拉运动方程式
- 欧拉运动定律
- 刚体动力学
- 玻恩刚性
- 刚性转子
- 平面刚体运动

## 参考文献

- J.L. Meriam, L.G. Kraige, "Engineering Mechanics: Dynamics,"第三版，ISBN 978-0-471-59273-0。

## 外部链接

- 喷气推进实验室DARTS网页：空间算符代数