剛體

本條目中，向量與標量分別用粗體與斜體顯示。例如，位置向量通常用 $\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網頁：空間算符代數