參考系的形變【彈性應(yīng)變在彎曲時空下的推廣】
在彈性力學(xué)中,應(yīng)變是一個非常重要的幾何概念,我們通過胡克定律可以直接將其與應(yīng)力建立聯(lián)系,從而發(fā)展連續(xù)介質(zhì)的彈性理論。而一般彈性力學(xué)是建立在三維歐氏空間中的,因此這種應(yīng)變也是歐氏空間中的幾何概念,自然我們希望能將其推廣到一般的流形上去。這里介紹在廣義相對論背景下的應(yīng)變。
參考資料:《張量分析與彈性力學(xué)》- 申文斌,張朝玉
《The large?scale structure of space-time》- Stephen Hawking, George Ellis
*注:此處不涉及Fermi導(dǎo)數(shù)的概念,只需要了解基本的張量分析和黎曼幾何即可。
一:引入 Introduction
在四維時空中,取一個向量場,它給出了一個全等測地線集,稱為Congurence of geodesics,這些測地線的切矢場均為
,仿射參數(shù)為間隔
。從物理角度,測地線是在引力場中做自由運動的粒子的世界線,因此這些測地線給出的四維速度也均為
,在局部標(biāo)架下分量記為
。
同時,這些測地線可以用一個連續(xù)參數(shù)來區(qū)分,因此生成了一個單參變化群,它對應(yīng)了一個向量場
,具體參數(shù)值由其積分曲線與測地線的交點確定。

利用Levi-Civita聯(lián)絡(luò),在局部標(biāo)架下可以得到著名的Jacobi方程或測地偏離方程:
,其中
是沿測地線的二階協(xié)變導(dǎo)數(shù),
為黎曼曲率張量。
并且由定義,向量場在向量場
的單參變化群下保持不變,也就是
,利用無撓性定義得
,在局部標(biāo)架下為無撓性方程?:
。
現(xiàn)在考慮單個自由粒子在測地線集中的狀態(tài)。從曲線系的角度出發(fā),在時空中的每一點,該粒子所觀察到周圍的測地線都是不同的,時而遠(yuǎn)離時而扭曲,這些都可以反映到向量場在各點的值上。而且很容易觀察到,沿著測地線切矢,也就是在速度的方向,是觀察不到這些變化的,因為向量場
在向量場
上同樣是不變的。
因此,總結(jié)可知,向量場的變化發(fā)生在垂直于四維速度
的超曲面上。由此可以選擇這樣的一個四維正交標(biāo)架,以
為其中一個基,和超曲面
的交點
處的,垂直于
的三個基為初始基向量,以它們沿曲線的平移為另三個基,記作
,稱為無自轉(zhuǎn)參考系。由度規(guī)的平行性易知
,滿足正交要求,是一個正交子空間。于是,向量場
在垂直面上的分量為
,通過上文可見其分量會發(fā)生變化,滿足無撓性方程。下文中以
代表垂直分量,
代表四維分量。

根據(jù)Cauchy-Lipchitz定理,這樣的一階系統(tǒng)在確定初始條件后存在唯一解,寫作,其中
是一個3x3的單位矩陣滿足方程
。
可見向量場的變化完全被矢量的變化
控制,解出這個量就等于解出了相應(yīng)的變化。
二:三維歐氏空間分析?Analysis in 3D Euclidien Space
現(xiàn)在我們從彈性力學(xué)中應(yīng)變(strain)的角度來分析一下三維歐氏空間中的表達(dá)式。
首先在三維歐氏空間中,聯(lián)絡(luò)系數(shù)為0,則,并且不區(qū)分逆/協(xié)變指標(biāo),但啞標(biāo)體系和愛因斯坦求和約定依舊成立。
其次,我們將位移矢量記為,這里僅代表位置的移動,并不表示速度。當(dāng)然也可以表示速度,這樣能更好的適配流體力學(xué)中的粘性力項。
最后,三維歐氏空間下的坐標(biāo)系記為。
現(xiàn)取歐氏空間中的體積微元,
在平面上,稱沿軸方向的應(yīng)變?yōu)?strong>線應(yīng)變/正應(yīng)變,記為
,稱角度變化的應(yīng)變?yōu)?strong>角應(yīng)變/剪應(yīng)變,記為
,則對形變量近似到一階有:
,具體計算可查看參考資料。

引入應(yīng)變張量(strain tensor):,
容易證明在歐式空間下這是一個對稱的雙線性型,同時稱這六個獨立分量與三個位移量的關(guān)系方程為Cauchy方程。
同時,在應(yīng)變前后體積也會發(fā)生變化,算上所有的線應(yīng)變和角應(yīng)變,展開到一階后可以得到體積膨脹系數(shù)/體積應(yīng)變:,可以看到它是應(yīng)變張量的跡,因此是一個應(yīng)變不變量。具體計算是以混合積表示體積展開即可。
對于一個對稱雙線性型,從矩陣的角度可以表示為有跡與無跡部分之和,以張量實體表示為,其中
稱為偏斜張量/應(yīng)變偏量,是應(yīng)變張量的無跡部分,其分量形式為
。
最后,考慮到在歐氏空間下偏導(dǎo)數(shù)是一個二階張量的分量,應(yīng)變張量是其對稱部分,定義其反對稱部分為轉(zhuǎn)動張量,于是根據(jù)二階張量分解定理有
。

還可定義相應(yīng)的角速度,為轉(zhuǎn)動張量的對偶矢量,星號為Hodge星算子。
至此三維歐氏空間中的位移變化就解完了,其余內(nèi)容在書中都有介紹,這里就不寫了。
三:四維廣義黎曼流形分析?Analysis in?4D Riemannian Manifold
現(xiàn)在考慮引力場下的情形,我們將二中的內(nèi)容推廣至任意廣義黎曼流形上。
首先,在彎曲時空中聯(lián)絡(luò)系數(shù)不為0,因此偏導(dǎo)數(shù)不再代表一個二階張量,同時逆/協(xié)變指標(biāo)也需要認(rèn)真的區(qū)分。
其次,根據(jù)一中的四維正交標(biāo)架,這里我們主要研究在子空間
上的分量變化,而其上的度規(guī)是對應(yīng)的超曲面
上的誘導(dǎo)度規(guī),因此下文中均默認(rèn)使用誘導(dǎo)度規(guī)做英文字母指標(biāo)的升降。
最后,對二階張量的張量指標(biāo)簡化記號,以圓括號表示對稱化:,以對易子表示反對稱化:
。
引入投影映射張量,其中
是超曲面的單位法向量,利用該張量,垂直部分可以寫作
,這里本質(zhì)上是施密特正交化定理,
是考慮到超曲面的性質(zhì),類時取正號。考慮到
在標(biāo)架內(nèi)的指標(biāo)只可以是0,因此上述投影只是將四維指標(biāo)換為三維值。
由二中的名稱和定義,給出在流形中的情形:
應(yīng)變張量 Expansion Tensor :;?
轉(zhuǎn)動張量 Vorticity Tensor :,角速度
;
體積系數(shù) Volume expansion :;
偏斜張量 Shear Tensor :。
于是在測地線下可得到。
若定義測地線簇之間的相對速度為,則根據(jù)無撓性方程易知
。
現(xiàn)在考慮這三項的物理意義,我們記為
上的單位向量,則:
,說明此時測地線上各點的相對速度在相對間隔上的投影為一個定值,且無其他分量,代表了體積的膨脹效應(yīng);
,說明此時測地線上各點的相對速度與相對間隔垂直,根據(jù)圓周運動模型,此時代表了旋轉(zhuǎn)效應(yīng);
,說明此時測地線上各點的相對速度在相對間隔上的投影分量不唯一,但保持體積不變,代表了剪切效應(yīng)。
或者也可以取法坐標(biāo)系,或者測地坐標(biāo)系,結(jié)合二中的內(nèi)容也能得出一樣的結(jié)論。
四:進(jìn)一步的計算 Further Calculation
現(xiàn)在我們來考慮用來表示上述的所有張量,記號依舊不變,記
的逆為
。
應(yīng)變張量 Expansion Tensor :;?
轉(zhuǎn)動張量 Vorticity Tensor :;
體積系數(shù) Volume expansion :。
偏斜張量 Shear Tensor :。
利用Jacobi方程可以求出,同時由這個方程做進(jìn)一步的變化可以得到在測地線下的情況:
?;?
?;?
以及Raychaudhuri-Landau方程,其中
。
于是在知道了曲率張量之后,我們就可以通過上述方程計算具體張量在場中的傳播速度,這對與研究引力波對時空的擾動效應(yīng)由重要的啟示。
最后提一下,這里的內(nèi)容對于理想牛頓粘性流體也一樣有效,考慮到平面牛頓粘性力方程,這里切應(yīng)力與速度的梯度有關(guān),在三維中是
,
自然會想到,對于引力場中的粘性流體,其粘性張量項自然為。并且,對于速度梯度的分解也和流體力學(xué)中關(guān)于流場旋度的描述一致,是微分幾何在流體力學(xué)中的應(yīng)用的另一個有力的證明。
此處我們完全是從測地線的角度出發(fā)的,對于任意的積分曲線簇,還應(yīng)加上對應(yīng)的加速度項,即,這樣就可以得到Jacobi方程,以及各張量的傳播速度了。