在分析力學中,質量矩陣是質量到廣義坐標概念上的推廣,它給出了系統廣義坐標q的變化率和系統動能T的關係,即
![{\displaystyle T={\frac {1}{2}}{\dot {q}}^{\mathrm {T} }M{\dot {q}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c11cc9746d902c7f20da99d2bc46ec5f68e9a89d)
其中
是向量
的 轉置 [1]。若粒子質量
,速度v, 那麼單個粒子的動能T為
![{\displaystyle T\;=\;{\frac {1}{2}}m|v|^{2}\;=\;{\frac {1}{2}}v\cdot mv}](https://wikimedia.org/api/rest_v1/media/math/render/svg/28990fe7cf55cc6c2ba239fc86a7895a89c4c3fd)
將系統中每個粒子的位置用q表示,可推導出上述的一般關係式。
例如,在一維中討論兩體粒子系統。這樣一個系統的位置具有2個自由度,每個粒子的位置都可由廣義位置矢量描述:
![{\displaystyle \mathbf {x} =[x_{1}\,x_{2}]^{\mathrm {T} }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/42e7c751c33bfa097782968cd6355b0100a8837a)
假設粒子具有質量
和
,系統的動能為
![{\displaystyle E=\sum _{i=1}^{2}{\frac {1}{2}}m_{i}{\dot {x_{i}}}^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/733d0f53094635ee1d2938941955877f7d20f5b2)
將質量列寫成矩陣
![{\displaystyle M={\begin{bmatrix}m_{1}&0\\0&m_{2}\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/85222657a80f7986b7c88151e9fdcb39fd546c3b)
那麼總動能由下列公式給出:
![{\displaystyle E={\frac {1}{2}}{\dot {\mathbf {x} }}^{\mathrm {T} }M{\dot {\mathbf {x} }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5504e620a3fcc4011c54ec0bd26b10f2f9b75ea4)
在多維情況下,質量矩陣會變得更為複雜。例如,在二維情況下,一個給定的粒子具有兩個自由度,因此,如果第i 個粒子對應自由度j 和j+1,那麼
![{\displaystyle M_{j,j}=M_{j+1,j+1}=m_{i}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8af0e19f58dd5a7511530448652948d6bf5fab95)
在如剛體動力學之類質量是分布式的情況下的應用中,非對角線元素非零的情況也是存在的。
一般地,質量矩陣M依賴於位置矢量q,且隨時間變化。拉格朗日力學中,常微分方程(組)描述了在系統中由粒子的位置所定義的廣義坐標矢量隨時間的變化。方程中的動能公式表示所有粒子的總動能。
示例
二體線性系統
一維空間中的質量系統.
考慮由僅限於直線軌道的兩個點狀物體。這樣一個系統的位置具有2個自由度,每個粒子的位置都可由廣義位置矢量q描述,即兩個粒子沿著軌道的位置:
.
設粒子的質量分別為 m1, m2, 系統的總動能是
![{\displaystyle T=\sum _{i=1}^{2}{\frac {1}{2}}m_{i}{\dot {x}}_{i}{}^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1d5a11cd95925ee71e0b9de50855fe891ffb4256)
也可以寫為
![{\displaystyle T={\frac {1}{2}}{\dot {q}}^{\mathrm {T} }M{\dot {q}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c11cc9746d902c7f20da99d2bc46ec5f68e9a89d)
其中
![{\displaystyle M={\begin{bmatrix}m_{1}&0\\0&m_{2}\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/85222657a80f7986b7c88151e9fdcb39fd546c3b)
N體系統
更一般地,考慮一個N個粒子的系統,記標號分別為i=1, 2, ..., N,其中粒子 i 的位置是Ni個自由直角坐標(其中Ni是1,2,或3)。令q是由坐標形成的列向量。質量矩陣M是對角塊矩陣,每個塊中的對角元素是對應的粒子的質量:[2]
![{\displaystyle M=\mathrm {diag} [m_{1}I_{n_{1}},m_{2}I_{n_{2}},\cdots ,m_{N}I_{n_{N}}]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8a322355d8350720da759df3cb6b845c51aff8fe)
其中 In i 是 ni × ni 單位陣,更具體地:
轉動系統
旋轉的啞鈴型系統.
考慮兩個點狀的物體,質量分別為m1,m2,連接到一個長度為2R的剛性無質量棒的兩端。該組件可以自由旋轉且在一個固定的平面上。該系統的狀態可以由廣義坐標向量描述如下:
![{\displaystyle q=[x,y,\alpha ]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e4db68719a08176ceb06d8bc377524718ee4697d)
其中 x, y 是以棒中點為原點的直角坐標系,角度 α 是棒轉動的角度。兩個粒子的位置和速度是:
![{\displaystyle {\begin{array}{ll}p_{1}=(x,y)+R(\cos \alpha ,\sin \alpha )&v_{1}=({\dot {x}},{\dot {y}})+R{\dot {\alpha }}(-\sin \alpha ,\cos \alpha )\\p_{2}=(x,y)-R(\cos \alpha ,\sin \alpha )&v_{2}=({\dot {x}},{\dot {y}})-R{\dot {\alpha }}(-\sin \alpha ,\cos \alpha )\end{array}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f23eeb4be978fd05f9e81620c6011e04e9e9202d)
它們的總動能是:
![{\displaystyle T=m{\dot {x}}^{2}+m{\dot {y}}^{2}+mR^{2}{\dot {\alpha }}^{2}+2Rd\cos \alpha {\dot {x}}{\dot {\alpha }}+2Rd\sin \alpha {\dot {y}}{\dot {\alpha }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/16ae3a59a54d34bd2d9cf9bf8ef2cbe0592943ce)
其中
,
。 這個公式寫成矩陣形式為:
![{\displaystyle T={\frac {1}{2}}{\dot {q}}^{\mathrm {T} }M{\dot {q}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c11cc9746d902c7f20da99d2bc46ec5f68e9a89d)
其中
![{\displaystyle M={\begin{bmatrix}m&0&Rd\cos \alpha \\0&m&Rd\sin \alpha \\Rd\cos \alpha &Rd\sin \alpha &R^{2}m\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d303c3b1c85446ca6751c96d717048b5db69aba6)
這裡矩陣依賴於棒的當前角度α。
對於連續介質力學的離散近似,如在有限元分析中,有多種方法可以構造質量方程,這取決於所期望的計算和精度性能。例如,利用忽略每一有限元變形的集中質量法可以構建一個對角質量矩陣並且不需要通過變形元來累積質量。
參考文獻
- ^
Mathematical methods for physics and engineering, K.F. Riley, M.P. Hobson, S.J. Bence, Cambridge University Press, 2010, ISBN 978-0-521-86153-3
- ^
Analytical Mechanics, L.N. Hand, J.D. Finch, Cambridge University Press, 2008, ISBN 978 0 521 57572 0
參見