跳转到内容

刚体:修订间差异

维基百科,自由的百科全书
删除的内容 添加的内容
老陳留言 | 贡献
內容擴充
老陳留言 | 贡献
內容擴充
第1行: 第1行:
{{noteTA|G1=物理學}}
{{noteTA|G1=物理學}}
[[File:Flight_dynamics_with_text.png|thumb|200px|right|刚体的空间组态决定于其[[质心]]位置与其[[定向(幾何)|定向]](最少有六个参数)。<ref name=Sciavicco>
[[File:Flight_dynamics_with_text.png|thumb|200px|right|刚体的空间组态决定于其[[质心]]位置与其[[定向 (幾何)|定向]](最少有六个参数)。<ref name=Sciavicco>
{{cite book |title=Modelling and control of robot manipulators |author=Lorenzo Sciavicco, Bruno Siciliano |chapter=§2.4.2 Roll-pitch-yaw angles |page=32 |isbn=1852332212 |year=2000 |edition=2nd |publisher=Springer}}
{{cite book |title=Modelling and control of robot manipulators |author=Lorenzo Sciavicco, Bruno Siciliano |chapter=§2.4.2 Roll-pitch-yaw angles |page=32 |isbn=1852332212 |year=2000 |edition=2nd |publisher=Springer}}
</ref>]]
</ref>]]
第67行: 第67行:
方向餘弦矩陣 <math>A</math> 足以設定附體參考系B的定向。但是,矩陣 <math>A</math> 有九個元素,而剛體只能供給三個自由度來做定向,因為這九個元素不是[[自變量]]({{lang|en|independent variable}}))。歐拉角的三個自變量可以用來設定剛體的定向。
方向餘弦矩陣 <math>A</math> 足以設定附體參考系B的定向。但是,矩陣 <math>A</math> 有九個元素,而剛體只能供給三個自由度來做定向,因為這九個元素不是[[自變量]]({{lang|en|independent variable}}))。歐拉角的三個自變量可以用來設定剛體的定向。


相對於空間參考系S,附體參考系B的定向,可以用三個歐拉角來設定。參右圖。設定 xyz-軸為空間參考系S的坐標軸,XYZ-軸為附體參考系B的坐標軸。稱 xy-平面與 XY-平面的相交為'''交點線''',用英文字母(N)代表。按照「zxz 順規」,歐拉角可以這樣定義:
相對於空間參考系S,附體參考系B的定向,可以用三個歐拉角來設定。參右圖。設定 xyz-軸為空間參考系S的坐標軸,XYZ-軸為附體參考系B的坐標軸。稱 xy-平面與 XY-平面的相交為'''交點線''',用英文字母(N)代表。按照「zxz 順規」,歐拉角可以這樣定義:
*<math>\alpha</math> 是 x-軸與交點線之間的夾角,
*<math>\alpha</math> 是 x-軸與交點線之間的夾角,
*<math>\beta</math> 是 z-軸與Z-軸之間的夾角,
*<math>\beta</math> 是 z-軸與Z-軸之間的夾角,
第90行: 第90行:


單獨分開工作,每個矩陣各自代表一種旋轉。按照順序相乘,
單獨分開工作,每個矩陣各自代表一種旋轉。按照順序相乘,
* 最裏面的(最右的) 矩陣代表繞著 z 軸的旋轉。
* 最裏面的最右的矩陣代表繞著 z 軸的旋轉。
* 最外面的(最左的) 矩陣代表繞著 Z 軸的旋轉。
* 最外面的最左的矩陣代表繞著 Z 軸的旋轉。
* 在中間的矩陣代表繞著交點線的旋轉。
* 在中間的矩陣代表繞著交點線的旋轉。


第111行: 第111行:
[[歐拉旋轉定理]]表明,在三維空間裏,假設約束剛體內部一點固定不動,則其任意位移等價於繞著某固定軸的一個旋轉,而這固定軸必含有這固定點。換句話說,設定附體參考系B的原點為這固定點,則附體參考系B不會因為這位移而涉及任何平移運動,再設定附體參考系B的z-軸與固定軸同軸,則這位移對應於繞著附體參考系B的z-軸旋轉 <math>\gamma</math> 角弧,而z-軸的方向是由 <math>\alpha</math> 與 <math>\beta</math> 角弧給出。<ref name=Whittaker_1917>{{cite book|last=Whittaker|first=Edmund |title=A treatise on the analytical dynamics of particles and rigid bodies; with an introduction to the problem of three bodies |publisher=Cambridge University Press|pages=2-5|date=1917|url= http://www.archive.org/details/cu31924001080294}}</ref>
[[歐拉旋轉定理]]表明,在三維空間裏,假設約束剛體內部一點固定不動,則其任意位移等價於繞著某固定軸的一個旋轉,而這固定軸必含有這固定點。換句話說,設定附體參考系B的原點為這固定點,則附體參考系B不會因為這位移而涉及任何平移運動,再設定附體參考系B的z-軸與固定軸同軸,則這位移對應於繞著附體參考系B的z-軸旋轉 <math>\gamma</math> 角弧,而z-軸的方向是由 <math>\alpha</math> 與 <math>\beta</math> 角弧給出。<ref name=Whittaker_1917>{{cite book|last=Whittaker|first=Edmund |title=A treatise on the analytical dynamics of particles and rigid bodies; with an introduction to the problem of three bodies |publisher=Cambridge University Press|pages=2-5|date=1917|url= http://www.archive.org/details/cu31924001080294}}</ref>


[[Image:Change of axes.svg|thumb|right|200px|改變剛體的定向與旋轉其附體參考系相同。]]
對於內部有一點被約束固定不動的剛體(或原點固定不動的參考系),歐拉旋轉定理將其任意位移等價為繞著某固定軸的一個旋轉。這允許使用旋轉來表達定向的改變。因此,變換矩陣 <math>A</math> 可以視為[[旋转矩阵#三维空间|三維旋轉]]的[[旋轉矩陣]],將附體參考系B或剛體做任意環繞著固定點的旋轉,從 <math>\hat{\mathbf{e}}_1</math> 、<math>\hat{\mathbf{e}}_2</math> 、<math>\hat{\mathbf{e}}_3</math> 旋轉成為 <math>\hat{\mathbf{e}}'_1</math> 、<math>\hat{\mathbf{e}}'_2</math> 、<math>\hat{\mathbf{e}}'_3</math> 。參考軸 <math>\hat{\mathbf{e}}'_i</math> 與 <math>\hat{\mathbf{e}}_j</math> 之間的關係為
對於內部有一點被約束固定不動的剛體(或原點固定不動的參考系),歐拉旋轉定理將其任意位移等價為繞著某固定軸的一個旋轉。這允許使用旋轉來表達定向的改變。因此,變換矩陣 <math>A</math> 可以視為[[旋转矩阵#三维空间|三維旋轉]]的[[旋轉矩陣]],將附體參考系B或剛體做任意環繞著固定點的旋轉,從 <math>\hat{\mathbf{e}}_1</math> 、<math>\hat{\mathbf{e}}_2</math> 、<math>\hat{\mathbf{e}}_3</math> 旋轉成為 <math>\hat{\mathbf{e}}'_1</math> 、<math>\hat{\mathbf{e}}'_2</math> 、<math>\hat{\mathbf{e}}'_3</math> 。參考軸 <math>\hat{\mathbf{e}}'_i</math> 與 <math>\hat{\mathbf{e}}_j</math> 之間的關係為
:<math>\hat{\mathbf{e}}'_i=a_{ij}\hat{\mathbf{e}}_j</math> 。
:<math>\hat{\mathbf{e}}'_i=a_{ij}\hat{\mathbf{e}}_j</math> 。
第120行: 第121行:
其中,<math>\mathbf{r}_{G}</math> 、<math>\mathbf{r}_{P/G}</math> 分別是點G的位置、從點G到點P的位移。
其中,<math>\mathbf{r}_{G}</math> 、<math>\mathbf{r}_{P/G}</math> 分別是點G的位置、從點G到點P的位移。


從附體參考系B觀測,剛體內部每一點的位置都固定不變。剛體從時間 <math>t_0</math> 到時間 <math>t</math> 的運動,可以分為點G從 <math>\mathbf{r}_{G}(t_0)</math> 到 <math>\mathbf{r}_{G}(t)</math> 的平移運動,與位移 <math>\mathbf{r}_{P/G}</math> 從時間 <math>t_0</math> 到時間 <math>t</math> 的旋轉運動。
<math>\mathbf{r}_{P/G}</math> 可以分解
:<math>\mathbf{r}_{P/G}=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</math> ,

從附體參考系B觀測,剛體內部每一點的位置都固定不變,所以,<math>(e_1,e_2,e_3)</math> 是常向量。剛體從時間 <math>t_0</math> 到時間 <math>t</math> 的運動,可以分為點G從 <math>\mathbf{r}_{G}(t_0)</math> 到 <math>\mathbf{r}_{G}(t)</math> 的平移運動,與位移 <math>\mathbf{r}_{P/G}</math> 從時間 <math>t_0</math> 到時間 <math>t</math> 的旋轉運動。


===平移速度與角速度===
===平移速度與角速度===
從不同的參考系觀測剛體運動的平移速度角速度是從某參考系可能不同,為了確保測量結果具有實際物理意義,必需先定參考系。
從不同的參考系觀測剛體運動,可能會獲得不同的平移速度和不同的角速度為了確保測量結果具有實際物理意義,必需先定參考系。


剛體的平移速度是向量,是其位置向量的時間變化率,是附著於剛體的附體參考系的原點的速度。對於純平移運動(沒有任何旋轉運動),剛體內部所有點的移動速度相同。假設涉及旋轉運動,則通常剛體內部任意兩點的瞬時速度不相等;只有當它們恰巧處於同一直軸,而這直軸平行於[[轉動瞬軸]]({{lang|en|instantaneous axis of rotation}}),則它們的瞬時速度相等。
剛體的平移速度是向量,是其位置向量的時間變化率,是附著於剛體的附體參考系的原點的速度。對於純平移運動(沒有任何旋轉運動),剛體內部所有點的移動速度相同。假設涉及旋轉運動,則通常剛體內部任意兩點的瞬時速度不相等;只有當它們恰巧處於同一直軸,而這直軸平行於[[轉動瞬軸]]({{lang|en|instantaneous axis of rotation}}),則它們的瞬時速度相等。


角速度也是向量,描述剛體定向改變的角速率,以及剛體旋轉時的瞬時轉軸([[歐拉旋轉定理]]保證瞬時轉軸的存在)。在任意時間,剛體內部每一個粒子的角速度相同。
角速度也是向量,描述剛體定向改變的角速率,以及剛體旋轉時的瞬時轉軸([[歐拉旋轉定理]]保證瞬時轉軸的存在)。在任意時間,剛體內部每一個粒子的角速度相同。

====向量的時間變化率====
假設一剛體呈純旋轉運動,其附體參考系B也會跟著旋轉,因此,對於任意向量 <math>\mathbf{F}</math> ,從這附體參考系B與從空間參考系S觀測,會得到不同的結果。假設附體參考系B <math>(\hat{\mathbf{e}}_1,\hat{\mathbf{e}}_2,\hat{\mathbf{e}}_3)</math> 與空間參考系S <math>(\hat{\mathbf{x}}_1,\hat{\mathbf{x}}_2,\hat{\mathbf{x}}_3)</math> 同原點。對於這些參考系,三維含時向量 <math>\mathbf{F}(t)</math> 分解為
:<math>\mathbf{F}=f_i\hat{\mathbf{x}}_i=F_i\hat{\mathbf{e}}_i</math> 。

<math>\mathbf{F}(t)</math> 對於時間的導數
:<math>\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}</math> 。

單獨計算附體參考軸對於時間的導數:
:<math>\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</math> ;

其中,<math>\dot{a}_{ij}\ \stackrel{def}{=}\ \frac{\mathrm{d}}{\mathrm{d}t}a_{ij}</math> 是方向餘弦對於時間的導數。

由於 <math>\frac{\mathrm{d}}{\mathrm{d}t}\hat{\mathbf{e}}_i</math> 垂直於 <math>\hat{\mathbf{e}}_i</math> ,<math>\frac{\mathrm{d}}{\mathrm{d}t}\hat{\mathbf{e}}_i</math> 只能是其他兩個單位向量的[[線性組合]]:
:<math>\frac{\mathrm{d}}{\mathrm{d}t}\hat{\mathbf{e}}_i=-\epsilon_{ij\ell}\omega_{ij}\hat{\mathbf{e}}_{\ell}</math> ;

其中,<math>\epsilon_{ij\ell}</math> 是[[列維-奇維塔符號]],<math>\omega_{ij}</math> 是係數。

對於任意 <math>\hat{\mathbf{e}}_{m}\ne \hat{\mathbf{e}}_i</math> ,
單位向量 <math>\hat{\mathbf{e}}_i</math> 與 <math>\hat{\mathbf{e}}_{m}</math> 的內積對於時間的導數為
:<math>\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) \\
& =-\epsilon_{ij\ell}\omega_{ij}\hat{\mathbf{e}}_{\ell} \cdot\hat{\mathbf{e}}_{m}
-\hat{\mathbf{e}}_{i} \cdot\epsilon_{mj\ell}\omega_{mj}\hat{\mathbf{e}}_{\ell} \\
& =-\epsilon_{ijm}\omega_{ij}-\epsilon_{mji}\omega_{mj} \\
& =-\epsilon_{ijm}(\omega_{ij}-\omega_{mj}) \\
& =0 \\
\end{align}</math>

所以, <math>\omega_{ij}</math> 的下標 <math>i</math> 多餘無用,可以刪除,變為 <math>\omega_{j}</math> :
:<math>\frac{\mathrm{d}}{\mathrm{d}t}\hat{\mathbf{e}}_i=-\epsilon_{ij\ell}\omega_{j}\hat{\mathbf{e}}_{\ell} </math> 。

思考 <math>\frac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}</math> 方程式最右邊項目 <math>F_i\frac{\mathrm{d}\hat{\mathbf{e}}_i}{\mathrm{d}t}</math> ,對換傀標 <math>i</math> <math>\ell</math> ,可以得到
:<math>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}</math> 。

向量 <math>\boldsymbol{\omega}</math> 是由三個係數 <math>\omega_1</math> 、<math>\omega_2</math> 、<math>\omega_3</math> 組成,對應於附體參考系的三個參考軸 <math>\hat{\mathbf{e}}_1</math> 、<math>\hat{\mathbf{e}}_2</math> 、<math>\hat{\mathbf{e}}_3</math> ,係數數值可以從歐拉角計算求得:
:<math>\omega_1=\dot{a}_{2j}a_{3j}
=\dot{\alpha}\sin{\beta}\sin{\gamma}+\dot{\beta}\cos{\gamma}</math> 、
:<math>\omega_2=-\dot{a}_{1j}a_{3j}
=\dot{\alpha}\sin{\beta}\cos{\gamma}+\dot{\beta}\sin{\gamma}</math> 、
:<math>\omega_3=\dot{a}_{1j}a_{2j}=\dot{\alpha}\cos{\beta}+\dot{\gamma}</math> 。

試想對應於歐拉角 <math>\alpha</math> 、<math>\beta</math> 、<math>\gamma</math> 的三個旋轉軸分別為 <math>\hat{\mathbf{z}}</math> 、<math>\hat{\mathbf{N}}</math> 、<math>\hat{\mathbf{Z}}</math> ,三個角速度分別為
:<math>\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}}
</math> 、
:<math>\boldsymbol{\omega}_{\beta}=\dot{\beta}\hat{\mathbf{N}}
=\dot{\beta}\cos{\gamma}\hat{\mathbf{X}}-\dot{\beta}\sin{\gamma}\hat{\mathbf{Y}}
</math> 、
:<math>\boldsymbol{\omega}_{\gamma}=\dot{\gamma}\hat{\mathbf{Z}}</math> 。

這三個角速度的向量和,對於附體參考系B的分量分別為
:<math>\omega_X=\dot{\alpha}\sin{\beta}\sin{\gamma}+\dot{\beta}\cos{\gamma}=\omega_1</math> 、
:<math>\omega_Y=\dot{\alpha}\sin{\beta}\cos{\gamma}+\dot{\beta}\sin{\gamma}=\omega_2</math> 、
:<math>\omega_Z=\dot{\alpha}\cos{\beta}+\dot{\gamma}=\omega_3</math> 。

注意到附體參考系B的 <math>\hat{\mathbf{e}}_1</math> 、<math>\hat{\mathbf{e}}_2</math> 、<math>\hat{\mathbf{e}}_3</math> 就是歐拉角的 <math>\hat{\mathbf{X}}</math> 、<math>\hat{\mathbf{Y}}</math> 、<math>\hat{\mathbf{Z}}</math> ,所以,向量 <math>\boldsymbol{\omega}</math> 是附體參考系B旋轉的角速度。

總結,向量 <math>\mathbf{F}(t)</math> 對於時間的導數為
:<math>\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}</math> 。

設定 <math>\left(\frac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}\right)_{\mathrm{space}}</math> 、<math>\left(\frac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}\right)_{\mathrm{body}}</math> 分別為從空間參考系S、附體參考系B觀測到的向量 <math>\mathbf{F}(t)</math> 對於時間的導數,上述方程式可以表達為
:<math>\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}</math> 。

由於向量 <math>\mathbf{F}(t)</math> 是任意向量,可以將 <math>\left(\frac{\mathrm{d}}{\mathrm{d}t}\right)_{\mathrm{space}}</math> 、<math>\left(\frac{\mathrm{d}}{\mathrm{d}t}\right)_{\mathrm{body}}</math> 當作算符,這樣,對應的算符方程式的形式為:
:<math>\left(\frac{\mathrm{d}}{\mathrm{d}t}\right)_{\mathrm{space}}=\left(\frac{\mathrm{d}}{\mathrm{d}t}\right)_{\mathrm{body}}+\boldsymbol{\omega}\times</math> 。

這算符方程式可以作用於任意含時向量。


==动力学==
==动力学==

2011年11月10日 (四) 06:38的版本

刚体的空间组态决定于其质心位置与其定向(最少有六个参数)。[1]

物理学裏,理想刚体rigid body)是一種有限尺寸,可以忽略形变固体。不论是否感受到外力,在刚体內部,与点之间的距离都不会改变。根據相對論,這種物體不可能實際存在,但物體通常可以假定為完美剛體,前提是必須滿足運動速度超小於光速的條件。

经典力学裡,刚体通常被认为是连续质量分佈体;在量子力学裏,刚体被認為是一群粒子的聚集。例如,分子(由假定為質點电子核子组成)时常會被视为刚体(請参阅條目分子的分类为刚性转子classification of molecules as rigid rotors))。

运动学

位置與定向

剛體是由一群數量超多的粒子組成。實際而言,不可能精確地追蹤其中每一個粒子的運動。為了簡化運算,可以利用剛體的剛性,即其內部所有粒子彼此之間距離不變的性質。假若物體具有剛性,則倚靠設定三個不同線粒子的位置,就足以設定此物體的位置。這意味著,在三維空間裏,剛體至多只有九個自由度,但由於假定三個粒子之間的距離固定不變,所以,剛體只有六個自由度。假設還有其它約束,例如,剛體的運動必需繞著其內部一點旋轉,或繞著其內部一直軸旋轉,則自由度會小於六。

關於其它任意粒子P的位置,只要知道粒子P與上述三個粒子之中的任意一個粒子的相對位置,就可以重建這粒子的位置。通常,整個剛體的空間組態可以簡易地以以下參數設定:

  1. 剛體的「位置」:挑選剛體內部一點G來代表整個剛體,通常會設定物體的質心形心為這一點。從空間參考系S觀測,點G的位置就是整個剛體在空間的位置。表示位置可以應用向量的概念。向量的起點為參考系S的原點,終點為點G。設定剛體的位置需要三個坐標,例如,採用直角坐標系,這三個坐標為x-坐標、y-坐標、z-坐標。這用掉了三個自由度。
  2. 剛體的定向:描述剛體定向的方法有好幾種,包括方向餘弦歐拉角四元數等等。這些方法設定一個附體參考系B的定向(相對於空間參考系S)。附體參考系是固定於剛體的參考系。相對於剛體,附體參考系的定向固定不變。由於剛體可能會呈加速度運動,所以附體參考系可能不是慣性參考系。空間參考系是某設定慣性參考系,例如,在觀測飛機的飛行運動時,附著於飛機場控制塔的參考系可以設定為空間參考系,而附著於飛機的參考系則可設定為附體參考系。剛體的定向需要用到另外三個自由度。

方向餘弦

方向餘弦方法可以用來設定,附體參考系B的定向,從而設定剛體的定向。假設沿著參考系S的坐標軸的三個單位向量分別為 ,沿著參考系B的坐標軸的三個單位向量分別為 。定義 之間的方向餘弦

其中, 之間的夾角。

之間的關係分別為

兩個參考系的坐標軸所形成的矩陣稱為「方向餘弦矩陣」

採用愛因斯坦求和約定,由於 ,給定方向餘弦矩陣 ,則可設定附體參考系B的定向,也就是剛體的定向。

反過來,經過一番運算,可以得到

給定位置向量

內積

方向餘弦矩陣 可以將從空間參考系S觀測的位置坐標 ,變換為從附體參考系B觀測的位置坐標 ,因此又稱為「變換矩陣」。

變換矩陣 也可以做反變換如下:

變換矩陣 是一種正交矩陣,滿足「正交條件」

其中,克羅內克函數

注意到 不同,夾角 與空間參考系S的坐標軸單位向量 之間的夾角。變換矩陣 通常不是對稱矩陣

左圖顯示「主動變換」:參考軸不變,點P被旋轉角弧 成為點P'。右圖顯示「被動變換」:參考軸被旋轉角弧 ,而點P保持不變。

對於二維旋轉,變換矩陣 可以視為旋轉矩陣。例如,將附體參考系B或剛體旋轉,從 旋轉角弧 成為 ;其中, 。對於這旋轉,旋轉矩陣

參考軸 之間的關係為

換一方面,旋轉矩陣 也可以視為將點P的位置向量 旋轉角弧 成為點P'的位置向量

歐拉角

按照順序使用最多三個歐拉角 () ,可以從xyz-軸變換到XYZ-軸。交點線 (N)是以綠色表示 。

方向餘弦矩陣 足以設定附體參考系B的定向。但是,矩陣 有九個元素,而剛體只能供給三個自由度來做定向,因為這九個元素不是自變量independent variable))。歐拉角的三個自變量可以用來設定剛體的定向。

相對於空間參考系S,附體參考系B的定向,可以用三個歐拉角來設定。參閱右圖。設定 xyz-軸為空間參考系S的坐標軸,XYZ-軸為附體參考系B的坐標軸。稱 xy-平面與 XY-平面的相交為交點線,用英文字母(N)代表。按照「zxz 順規」,歐拉角可以這樣定義:

  • 是 x-軸與交點線之間的夾角,
  • 是 z-軸與Z-軸之間的夾角,
  • 是交點線與X-軸之間的夾角。

每一個歐拉角的旋轉都對應於一個簡單的旋轉矩陣:

設定剛體定向的旋轉矩陣 是由三個簡單旋轉矩陣 共同合成:

單獨分開工作,每個矩陣各自代表一種旋轉。按照順序相乘,

  • 最裏面的(最右的)矩陣代表繞著 z 軸的旋轉。
  • 最外面的(最左的)矩陣代表繞著 Z 軸的旋轉。
  • 在中間的矩陣代表繞著交點線的旋轉。

經過一番運算,可以得到 矩陣:

逆矩陣是:

歐拉旋轉定律

歐拉旋轉定理表明,在三維空間裏,假設約束剛體內部一點固定不動,則其任意位移等價於繞著某固定軸的一個旋轉,而這固定軸必含有這固定點。換句話說,設定附體參考系B的原點為這固定點,則附體參考系B不會因為這位移而涉及任何平移運動,再設定附體參考系B的z-軸與固定軸同軸,則這位移對應於繞著附體參考系B的z-軸旋轉 角弧,而z-軸的方向是由 角弧給出。[2]

改變剛體的定向與旋轉其附體參考系相同。

對於內部有一點被約束固定不動的剛體(或原點固定不動的參考系),歐拉旋轉定理將其任意位移等價為繞著某固定軸的一個旋轉。這允許使用旋轉來表達定向的改變。因此,變換矩陣 可以視為三維旋轉旋轉矩陣,將附體參考系B或剛體做任意環繞著固定點的旋轉,從 旋轉成為 。參考軸 之間的關係為

沙勒定理

當剛體移動時,它的位置與定向都可能會隨著時間演進而改變。在運動學裏,這可分為平移運動與旋轉運動。根據沙勒定理Chasles' theorem),歐拉旋轉定律的推論,剛體的最廣義位移等價於一個平移加上一個旋轉。[2]剛體的現在位置與現在定向可以視為是從某個初始位置與初始定向經過平移與旋轉而成,這平移與旋轉可能不依循物體的運動軌道。挑選剛體內部一點G來代表整個剛體,設定附體參考系B的原點於點G,則從空間參考系S觀測,在剛體內部任意一點P的位置

其中, 分別是點G的位置、從點G到點P的位移。

從附體參考系B觀測,剛體內部每一點的位置都固定不變。剛體從時間 到時間 的運動,可以分為點G從 的平移運動,與位移 從時間 到時間 的旋轉運動。

平移速度與角速度

從不同的參考系觀測剛體運動,可能會獲得不同的平移速度和不同的角速度。為了確保測量結果具有實際物理意義,必需先給定參考系。

剛體的平移速度是向量,是其位置向量的時間變化率,是附著於剛體的附體參考系的原點的速度。對於純平移運動(沒有任何旋轉運動),剛體內部所有點的移動速度相同。假設涉及旋轉運動,則通常剛體內部任意兩點的瞬時速度不相等;只有當它們恰巧處於同一直軸,而這直軸平行於轉動瞬軸instantaneous axis of rotation),則它們的瞬時速度相等。

角速度也是向量,描述剛體定向改變的角速率,以及剛體旋轉時的瞬時轉軸(歐拉旋轉定理保證瞬時轉軸的存在)。在任意時間,剛體內部每一個粒子的角速度相同。

向量的時間變化率

假設一剛體呈純旋轉運動,其附體參考系B也會跟著旋轉,因此,對於任意向量 ,從這附體參考系B與從空間參考系S觀測,會得到不同的結果。假設附體參考系B 與空間參考系S 同原點。對於這些參考系,三維含時向量 分解為

對於時間的導數為

單獨計算附體參考軸對於時間的導數:

其中, 是方向餘弦對於時間的導數。

由於 垂直於 只能是其他兩個單位向量的線性組合

其中,列維-奇維塔符號 是係數。

對於任意 , 單位向量 的內積對於時間的導數為

所以, 的下標 多餘無用,可以刪除,變為

思考 方程式最右邊項目 ,對換傀標 ,可以得到

向量 是由三個係數 組成,對應於附體參考系的三個參考軸 ,係數數值可以從歐拉角計算求得:

試想對應於歐拉角 的三個旋轉軸分別為 ,三個角速度分別為

這三個角速度的向量和,對於附體參考系B的分量分別為

注意到附體參考系B的 就是歐拉角的 ,所以,向量 是附體參考系B旋轉的角速度。

總結,向量 對於時間的導數為

設定 分別為從空間參考系S、附體參考系B觀測到的向量 對於時間的導數,上述方程式可以表達為

由於向量 是任意向量,可以將 當作算符,這樣,對應的算符方程式的形式為:

這算符方程式可以作用於任意含時向量。

动力学

主项目:刚体动力学

描述刚体的平移运动(位置、速度、加速度等等),刚体内任意一点皆可被用为参考点(L参考系的原点)。但是,根据实际需要,一个适当的选择是:

  • 整个刚体的质心
  • 平移运动是零,或可以简易研算的点。例如:在铰链上、在万向接头的中心等等。

当质心被选为参考点时:

刚体的转动

  • 刚体的转动定理:,其中为刚体所受合外力的力矩,为刚体转动惯量,为刚体角加速度。
  • 刚体的转动动能定理,其中表示合外力的力矩在角位移上所作的功,为刚体的转动惯量,为刚体角速度。
  • 刚体的转动和平动可以合成为刚体的平面运动,由柯尼希定理,其动能为,其中为刚体质心对参考系的速度。

参阅

參考文獻

  1. ^ Lorenzo Sciavicco, Bruno Siciliano. §2.4.2 Roll-pitch-yaw angles. Modelling and control of robot manipulators 2nd. Springer. 2000: 32. ISBN 1852332212. 
  2. ^ 2.0 2.1 Whittaker, Edmund. A treatise on the analytical dynamics of particles and rigid bodies; with an introduction to the problem of three bodies. Cambridge University Press. 1917: 2–5. 
  • J.L. Meriam, L.G. Kraige, "Engineering Mechanics: Dynamics,"第三版,ISBN 0471592730