變分法&能量原理(上)

本篇也需要用到一些張量內容和變分內容,給兩篇我之前的文章參考。

QuYln:解放雙手——愛因斯坦的張量表示

QuYln:彈性力學基礎(含張量與變分法)

一、應變能:

V_{varepsilon}=int_{V}^{}frac{1}{2}[frac{E}{left( 1+mu 
ight)} varepsilon_{ij} varepsilon_{ij} + frac{Emu}{(1+mu)(1-2mu)} 	heta^{2}] d V

這是張量(分量)記法,但是一般書用的不是張量應變,而是工程應變,修改一下:

varepsilon_{ij} varepsilon_{ij}=varepsilon_{11}^{2} +varepsilon_{22}^{2} + varepsilon_{33}^{2} + 2varepsilon_{12}^{2} + 2varepsilon_{13}^{2}+ 2 varepsilon_{23}^{2}

varepsilon_{11}^{2} +varepsilon_{22}^{2} + varepsilon_{33}^{2} + 2varepsilon_{12}^{2} + 2varepsilon_{13}^{2}+ 2 varepsilon_{23}^{2} =varepsilon_{x}^{2} + varepsilon_{y}^{2} +varepsilon_{z}^{2}

+ frac{1}{2} gamma_{xy}^{2} +frac{1}{2} gamma_{xz}^{2} +frac{1}{2} gamma_{yz}^{2}

因為工程切應變是張量切應變的二倍。

上式經過推導,得:

V_{varepsilon}=frac{1}{2} frac{E}{left ( 1+mu 
ight)} int_{V}^{}[ ( varepsilon_{x}^{2} + varepsilon_{y}^{2} +varepsilon_{z}^{2} )+frac{1}{2}( gamma_{xy}^{2} + gamma_{xz}^{2} + gamma_{yz}^{2} )

 + frac{mu}{(1-2mu)} 	heta^{2}] d V

接下來推導卡斯提安諾定理

frac{partial v_{varepsilon}}{partial varepsilon_{x}} = sigma_{x}

以此類推,應變能密度對於應變的導數,等於應力,其實就是材料力學裡的卡氏定理

推導的過程一定要注意這個步驟:

frac{partial(	heta^{2})}{partial varepsilon_{x}}=2	heta

因為:

	heta=varepsilon_{kk}=varepsilon_{x}+varepsilon_{y}+varepsilon_{z}

另外,我們要得到一個結論,結構整體的應變能,可以由位移場決定: V_{varepsilon}=frac{1}{2} frac{E}{left ( 1+mu 
ight)} int_{V}^{}[ ( varepsilon_{x}^{2} + varepsilon_{y}^{2} +varepsilon_{z}^{2} )+frac{1}{2}( gamma_{xy}^{2} + gamma_{xz}^{2} + gamma_{yz}^{2} )

 + frac{mu}{(1-2mu)} 	heta^{2}] d V

在上式中,只需要用幾何關係,將應變代換,即可得到位移表示的應變能。

二、應變余能:

V_{c}=frac{1}{2E} int_{V}^{}[(1+mu)sigma_{ij} sigma_{ij} - mu Theta^{2} ] dV

將張量的表達展開,如下:

V_{c}=frac{1}{2E} int_{V}^{}[2(1+mu)(	au_{xy}^{2}+ 	au_{xz}^{2} + 	au_{yz}^{2}) -2 mu (sigma_{x} sigma_{y}+sigma_{x} sigma_{z} +sigma_{y} sigma_{z})

+ (sigma_{x}^{2}+sigma_{y}^{2}+sigma_{z}^{2}) ] dV

同樣的,繼續推導卡斯提安諾定理,得到:

frac{partial v_{c}}{partial sigma_{x}} = varepsilon_{x}

通過應變余能對應力求導數,可以得到應變,這其實也是卡氏定理的另一種形式。

三、位移變分方程,虛位移原理與最小勢能原理:

先總體說一下我的理解,位移變分方程,虛位移原理,最小勢能原理,原方程加邊界條件,四者是完全等價的(可以通過變分法,把PDE邊值問題轉化為泛函問題)。

位移變分方程和虛位移原理是一種基於能量守恆的變體。

位移變分方程沒有把應變能的變分具體化,而虛位移原理則具體化了應變能的變分。

最小勢能原理實則表達了把位移變分方程或者虛位移的變分符號整體提出之後,所表達成的某個泛函的一階變分為零,這個泛函就是總體的勢能。

(1)那麼首先介紹位移變分原理

delta V_{varepsilon} = int_{V}^{} f_{i} delta u_{i} dV + int_{S_{sigma}}^{} ar{f_{i}} delta u_{i} dS

由於這裡,位移邊界的位移是給定的,故位移邊界上位移的變分(虛位移)為0,故外力的虛功不計算位移邊界上的部分。

(2)接下來,介紹虛位移原理

delta V_{varepsilon} = int_{V}^{} delta v_{varepsilon} d V

這個積分與變分符號可以調換順序,在前文證明過。

delta V_{varepsilon} = int_{V}^{} delta v_{varepsilon} d V = int_{V}^{} frac{partial v_{varepsilon}}{partial varepsilon_{i} }delta varepsilon_{i} d V (六項求和)

應變能密度的變分,其實類似於全微分,上式是需要求和的。

int_{V}^{} frac{partial v_{varepsilon}}{partial varepsilon_{i} }delta varepsilon_{i} d V = int_{V}^{} sigma_{i} delta varepsilon_{i} d V

此項的推導是基於卡氏定理,物理意義則是應力在虛應變上做的功。

將上表達式帶入位移變分方程,得到虛位移原理:

int_{V}^{} sigma_{i} delta varepsilon_{i} d V = int_{V}^{} f_{i} delta u_{i} dV + int_{S_{sigma}}^{} ar{f_{i}} delta u_{i} dS

物理意義就是,虛位移上外力做的虛功,等於虛應變上,應力做的虛功,整體的角度,就是能量守恆的一種體現。

(3)現在,介紹最小勢能原理

首先,給出一個錯誤的表達:

int_{V}^{} sigma_{i} delta varepsilon_{i} d V =deltaint_{V}^{} sigma_{i} varepsilon_{i} d V

這裡的變分符號,是不能隨便提出來的,這個錯誤原因是:

應力並不是常數,是不可以隨便和變分號交換的。

那麼,基於位移變分方程開始推導:

delta V_{varepsilon} = int_{V}^{} f_{i} delta u_{i} dV + int_{S_{sigma}}^{} ar{f_{i}} delta u_{i} dS

delta [V_{varepsilon} -( int_{V}^{} f_{i} u_{i} dV + int_{S_{sigma}}^{} ar{f_{i}} u_{i} dS)]=0

這裡的變分符號不僅和積分符號換序,又因為外力是常數,所以可以和變分符號換序。

V_{p}= -( int_{V}^{} f_{i} u_{i} dV + int_{S_{sigma}}^{} ar{f_{i}} u_{i} dS)

這一項稱之為外力勢能,定義彈性體未變形時的狀態為零勢能狀態。

規定:從當前狀態到零勢能狀態所做的功稱為勢能。

由於此處的位移表達式為變形位移,故「恢複位移」需要帶負號。

delta (V_{varepsilon} +V_{p}) =0

那麼定義彈性系統的總能量為E,則有:

E=V_{varepsilon} +V_{p}

並且,E是位移的泛函,在平衡狀態的時候,有:

delta E =0

可以證明:

delta^{2} E > 0

即,得到了最小勢能原理。

推薦閱讀:

二次方程, 3次方程, 4次方程的解法
更漂亮的十二,十二進之於生命的合理性
配極理論與調和點列
DeepLearning筆記:多節點神經網路

TAG:有限元分析FEA | 數學 | 力學 |