刚体

来自testwiki
跳转到导航 跳转到搜索

Template:NoteTA

本条目中,向量标量分別用粗體斜體顯示。例如,位置向量通常用 𝐫 表示;而其大小則用 r 來表示。
刚体的空间位形决定于其质心位置与其取向(最多有六个自由度)。[1]

物理学裏,理想刚体Template:Lang-en)是一種有限尺寸,可以忽略形变固体。不论是否感受到外力,在刚体內部,質點與質點之间的距离都不会改变。这种理想模型适用条件是,运动过程比固体中的弹性波的传播要缓慢得多。根據相對論,這種物體不可能實際存在,但物體通常可以假定為完美剛體,前提是必須滿足運動速度遠小於光速的條件。

经典力学裡,刚体通常被視為连续质量分佈体;在量子力学裏,刚体被視為一群粒子的聚集。例如,分子(由假定為質點电子核子组成)时常會被视为刚体。

运动学

位置與取向

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

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

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

方向餘弦

單位向量𝐞i^與參考系的三個單位向量𝐱^1𝐱^2𝐱^3之間的夾角分別為θi1θi2θi3

方向餘弦方法可以用來設定附體參考系B的取向,即剛體的取向。假設沿著參考系S的坐標軸的三個單位向量分別為𝐱^1𝐱^2𝐱^3,沿著參考系B的坐標軸的三個單位向量分別為𝐞^1𝐞^2𝐞^3。定義𝐞^i𝐱^j之間的方向餘弦aij

aij =def cos(θij)=𝐞^i𝐱^j

其中,θij𝐞^i𝐱^j之間的夾角。

𝐞^1𝐞^2𝐞^3𝐱^1𝐱^2𝐱^3之間的關係分別為

𝐞^1=cos(θ11)𝐱^1+cos(θ12)𝐱^2+cos(θ13)𝐱^3=a11𝐱^1+a12𝐱^2+a13𝐱^3
𝐞^2=cos(θ21)𝐱^1+cos(θ22)𝐱^2+cos(θ23)𝐱^3=a21𝐱^1+a22𝐱^2+a23𝐱^3
𝐞^3=cos(θ31)𝐱^1+cos(θ32)𝐱^2+cos(θ33)𝐱^3=a31𝐱^1+a32𝐱^2+a33𝐱^3

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

A=[a11a12a13a21a22a23a31a32a33]

採用愛因斯坦求和約定,由於𝐞^i=aij𝐱^j,給定方向餘弦矩陣A,則可設定附體參考系B的取向,也就是剛體的取向。

反過來,經過一番運算,可以得到𝐱^j=aij𝐞^i

給定位置向量𝐫

𝐫=x1𝐱^1+x2𝐱^2+x3𝐱^3=e1𝐞^1+e2𝐞^2+e3𝐞^3

𝐫𝐞^i內積

𝐫𝐞^i=ei=ai1x1+ai2x2+ai3x3=aijxj

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

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

xj=aijei

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

aijaik=δjk

其中,δjk克羅內克函數

注意到θijθji不同,夾角θji𝐞^j與空間參考系S的坐標軸單位向量𝐱^i之間的夾角。變換矩陣A通常不是對稱矩陣

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

對於二維旋轉,變換矩陣A可以視為旋轉矩陣。例如,將附體參考系B或剛體旋轉,從𝐞^1𝐞^2𝐞^3旋轉θ角弧成為𝐞^'1𝐞^'2𝐞^'3;其中,𝐞^3=𝐞^'3。對於這旋轉,旋轉矩陣A

A=[cosθsinθ0sinθcosθ0001]

參考軸𝐞^'i𝐞^j之間的關係為

𝐞^'i=aij𝐞^j

旋轉矩陣A也可以視為將點P的位置向量𝐫=xi𝐱^i旋轉θ角弧成為點P'的位置向量𝐫=x'i𝐱^i

x'i=aijxj

歐拉角

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

Template:主條目 方向餘弦矩陣A足以設定附體參考系B的取向。但是,矩陣A有九個元素,而剛體只能供給三個自由度來設定取向,因為這九個元素不是[[自變量|Template:Tooltip]]。歐拉角的三個自變量可以用來設定剛體的取向。

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

  • α是x-軸與交點線(N)之間的夾角,
  • β是z-軸與Z-軸之間的夾角,
  • γ是交點線(N)與X-軸之間的夾角。

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

Aα=[cosαsinα0sinαcosα0001]
Aβ=[1000cosβsinβ0sinβcosβ]
Aγ=[cosγsinγ0sinγcosγ0001]

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

A=AγAβAα

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

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

經過一番運算,可以得到A矩陣:[2]

A=[cosαcosγcosβsinαsinγsinαcosγ+cosβcosαsinγsinβsinγcosαsinγcosβsinαcosγsinαsinγ+cosβcosαcosγsinβcosγsinβsinαsinβcosαcosβ]

A逆矩陣是:

A1=[cosαcosγcosβsinαsinγcosαsinγcosβsinαcosγsinβsinαsinαcosγ+cosβcosαsinγsinαsinγ+cosβcosαcosγsinβcosαsinβsinγsinβcosγcosβ]

歐拉旋轉定理

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

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

𝐞^'i=aij𝐞^j

沙勒定理

掛在重力式摩天輪邊緣的乘客座艙在做平移運動。
剛體平移運動示意圖。

當剛體移動時,它的位置與取向都可能會隨著時間演進而改變。[[沙勒定理|Template:Tooltip]]是歐拉旋轉定理的一個推論。根據沙勒定理,剛體的最廣義位移等價於一個平移加上一個旋轉。[3]因此,剛體運動可分為平移運動與旋轉運動。剛體的現在位置與現在取向可以視為是從某個初始位置與初始取向經過平移與旋轉而成。

如右圖所示,從時間t1到時間t2,當剛體在做平移運動時,任意內部兩點,點P與點Q的軌跡(以黑色實線表示)相互平行,線段PQ(以黑色虛線表示)的方向保持恆定。

挑選剛體內部一點G來代表整個剛體,設定附體參考系B的原點於點G(稱為「基點」),則從空間參考系S觀測,在剛體內部任意一點P的位置𝐫P

𝐫P=𝐫G+𝐫P/G

其中,𝐫G𝐫P/G分別是基點G的位置、點P對於基點G的相對位置。

從附體參考系B觀測,剛體內部每一點的位置都固定不變,但從空間參考系S觀測,剛體從時間t1到時間t2的運動,可以分為基點G從𝐫G(t1)𝐫G(t2)的平移運動,與位移𝐫P/G從時間t1到時間t2的旋轉運動。

平移速度與角速度

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

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

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

向量的時間變化率

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

𝐅=fi𝐱^i=Fi𝐞^i

𝐅(t)對於時間的導數為

d𝐅dt=dfidt𝐱^i=dFidt𝐞^i+Fid𝐞^idt

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

d𝐞^idt=ddt(aij𝐱^j)=a˙ij𝐱^j=a˙ijakj𝐞^k

其中,a˙ij =def daijdt是方向餘弦對於時間的導數。

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

ddt𝐞^i=ϵijωij𝐞^

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

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

ddt(𝐞^i𝐞^m)=(ddt𝐞^i)𝐞^m+𝐞^i(ddt𝐞^m)=ϵijωij𝐞^𝐞^m𝐞^iϵmjωmj𝐞^=ϵijmωijϵmjiωmj=ϵijm(ωijωmj)=0

所以,ωij的下標i多餘無用,可以刪除,變為ωj

d𝐞^idt=ϵijωj𝐞^

思考d𝐅dt方程式最右邊項目Fid𝐞^idt,對換傀標i ,可以得到

Fid𝐞^idt=ϵijFiωj𝐞^=ϵij𝐞^iωjF=ω×𝐅

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

ω1=a˙2ja3j=α˙sinβsinγ+β˙cosγ
ω2=a˙1ja3j=α˙sinβcosγ+β˙sinγ
ω3=a˙1ja2j=α˙cosβ+γ˙

試想對應於歐拉角αβγ的三個旋轉軸分別為𝐳^𝐍^𝐙^,三個角速度分別為

ωα=α˙𝐳^=α˙sinβsinγ𝐗^+α˙sinβcosγ𝐘^+α˙cosβ𝐙^
ωβ=β˙𝐍^=β˙cosγ𝐗^β˙sinγ𝐘^
ωγ=γ˙𝐙^

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

ωX=α˙sinβsinγ+β˙cosγ=ω1
ωY=α˙sinβcosγ+β˙sinγ=ω2
ωZ=α˙cosβ+γ˙=ω3

注意到附體參考系B的𝐞^1𝐞^2𝐞^3就是歐拉角的𝐗^𝐘^𝐙^,所以,向量ω是附體參考系B旋轉的角速度。

總結,向量𝐅(t)對於時間的導數為

d𝐅dt=dfidt𝐱^i=dFidt𝐞^i+ω×𝐅

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

(d𝐅dt)space=(d𝐅dt)body+ω×𝐅

項目 ω×𝐅可以想像為,從空間參考系S觀測,剛體內部位置向量為𝐅的質點,由於剛體旋轉而產生的速度。

向量𝐅(t)是任意向量,因此可以將(ddt)space(ddt)body當作算符,這樣,對應的算符方程式的形式為:

(ddt)space=(ddt)body+ω×

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

運動學方程式

根據沙勒定理,剛體的最廣義位移等價於一個平移加上一個旋轉。[3]挑選剛體內部一點G來代表整個剛體,設定附體參考系B的原點於基點G,則從空間參考系S觀測,在剛體內部任意一點P的位置𝐫P

𝐫P=𝐫G+𝐫P/G

其中,𝐫G𝐫P/G分別是基點G的位置、點P對於基點G的相對位置。

點P的速度𝐯P=(d𝐫Pdt)space

𝐯P=𝐯G+𝐯P/G

其中,𝐯G=(d𝐫Gdt)space𝐯P/G=(d𝐫P/Gdt)space分別是基點G的速度、點P對於基點G的相對速度。

應用前段推導出的適用於任意含時向量的算符方程式,可以計算出𝐯P/G。由於從附體參考系B觀測,剛體內部每一點的位置都固定不變,項目(d𝐫P/Gdt)body等於零:

𝐯P/G=(d𝐫P/Gdt)body+ω×𝐫P/G=ω×𝐫P/G

其中,ω是剛體的角速度向量。

所以,點P的速度為

𝐯P=𝐯G+ω×𝐫P/G

點P的加速度𝐚P=(d𝐯Pdt)space

𝐚P=𝐚G+𝐚P/G

其中,𝐚G=(d𝐯Gdt)space𝐚P/G=(d𝐯P/Gdt)space分別是基點G的速度、點P對於基點G的相對速度。

再應用前段推導出的算符方程式,可以計算出

𝐚P/G=(dωdt)space×𝐫P/G+ω×𝐯P/G=α×𝐫P/G+ω×(ω×𝐫P/G)

其中,α=(dωdt)space是附體參考系B旋轉的角加速度向量。

動力學

主项目:刚体动力学

當描述刚体的动力运动時,必需先處理妥善剛體的平移运动,即先選擇一個参考点來代表剛體與其附體參考系B。刚体内部任意一点都可以被選为参考点(附體參考系B的原点)。但是,根据实际應用需要,比較适当的选择為:

  • 刚体的质心:對於自由移動於空間的剛體,以其質心為參考點的方法通常會給出最簡單的運動。
  • 平移运动為零或可以简易研算的点:例如,在铰链上、在万向接头的中心等等。

当质心被选为参考点时:

  • 刚体的动量与其旋转运动无关。在任何时间,动量等于刚体的总质量乘以平移速度。
  • 不论刚体的平移运动为何,对于质心的角动量皆等同。所以,在计算角动量时,可以忽略平移运动。在任何时间,角动量等于惯性张量乘以角速度。假若知道刚体绕主轴的角速度,那麼,角动量对于每一主轴的分量,是对应的主慣性矩乘以对应的角速度;力矩是转动惯量乘以角加速度。
  • 在无外力作用下,可能形成的运动为等速直线运动、稳定[[绕定轴转动|Template:Tooltip]]、零力矩进动等等。
  • 作用于刚体的净外力,等于总质量乘以刚体平移运动的加速度(也就是说,不论净外力矩是否为零,或这刚体是否在作旋转运动,牛頓第二運動定律可以正确地应用于刚体平移运动,)。
  • 动能平移动能旋转动能的总和。

刚体的转动

  • 刚体的转动定律:M=Iα,其中M为刚体所受合外力的力矩,I为刚体转动惯量,α为刚体角加速度。
  • 刚体的转动动能定理W=θ1θ2Mdθ=12Iω2212Iω12,其中θ1θ2Mdθ表示合外力的力矩M在角位移θ2θ1上所作的功,I为刚体的转动惯量,ω为刚体角速度。
  • 刚体的转动和平动可以合成为刚体的平面运动,由柯尼希定理,其动能为Ek=12mVc2+12Iω2,其中Vc为刚体质心对参考系的速度。

参阅

參考文獻

Template:Reflist

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

外部連結


Template:經典力學