摘 要：建立應(yīng)變能密度積分公式，以彈性力學(xué)的物理方程為基礎(chǔ)，研究在彈性體中應(yīng)變能和應(yīng)變與應(yīng)力之前的關(guān)系。主要通過(guò)廣義格林公式下的積分和全微分形式下的積分。在這兩種方式的求解下，都證明了應(yīng)變能密度和應(yīng)變與應(yīng)力的關(guān)系是一致的。

關(guān)鍵詞：彈性體；應(yīng)變能密度；應(yīng)變

1 引言

在普遍的彈性力學(xué)的教材中[1]，都只給出了

這個(gè)公式，而沒(méi)有直接對(duì)應(yīng)變能密度的積分給出詳細(xì)的求解過(guò)程。這就使得讀者產(chǎn)生疑惑，當(dāng)應(yīng)力展開成與彈性模量的關(guān)系式代入積分方程中，就出現(xiàn)了問(wèn)題。

在 中，我們?cè)诶脧埩抗降恼归_：當(dāng) ？著kl的下標(biāo)k，l的取值和deij的下標(biāo)i，j的下標(biāo)取值相同：k=m=i；l=n=j；我們可以直接積分得出 （注：此公式不帶啞標(biāo)，且？著mn=？著nm），當(dāng)k取值不等與i和j，或者l取值不等于i和j時(shí)，那么？著kl是不是就與？著ij沒(méi)有關(guān)系，而可以直接把eij積出來(lái)。例如

。對(duì)于？著ij與？著kl取值不相同的部分，直接積分得出的結(jié)果不是■倍的關(guān)系而是1倍的關(guān)系。這就與應(yīng)變能密度的公式矛盾，顯然是我們把積分簡(jiǎn)單化了，因此我們要改變積分方式，由于彈性力學(xué)里的物理方程是肯定成立，所以我們從積分路徑考慮。

2 線彈性應(yīng)變

對(duì)圖1所示積分，從面積就知道 成立。但是對(duì)非線彈性應(yīng)變時(shí)，面積不是那么輕易算出，所以我們必須從積分本身出發(fā)。

3 非線彈性應(yīng)變

圖1 圖2

對(duì)應(yīng)變能密度 進(jìn)行證明，其中？滓ij=Cijkl？著kl，？滓ij為應(yīng)力，？著ij為應(yīng)變， Cijkl為彈性模量。

4 數(shù)學(xué)求解

首先對(duì)Cijkl進(jìn)行簡(jiǎn)化，通過(guò)？滓ij=Cijkl？著kl，對(duì)于各向異性彈性體應(yīng)力張量的對(duì)稱性？滓ij=？滓ji，于是Cijkl=Cjikl，再由應(yīng)變張量的對(duì)稱性？著kl=？著lk，得出Cijkl=Cijlk，于是把原來(lái)81個(gè)分量降為36個(gè)分量，整理成：

（1）

對(duì)U（？著）求二階偏導(dǎo)得出：

（2）

又因?yàn)?/p>

將 代入可得出 此式就稱之為廣義格林公式[2]

所以有

（3）

對(duì)U（？著）展開的

（4）

再將（1）式代入（4）式得

對(duì)此式整理

（5）

我們?nèi)∑渲腥我鈳捉M研究：

其中■x和■y是終止?fàn)顟B(tài)的彈性應(yīng)變，在求解

時(shí)，由（3）式得到C12=C21，所以由格林公式可知此積分與路徑無(wú)關(guān)。不妨設(shè)？著y=k？著x，

對(duì)其他部分采用同樣的方法，整理得

（6）

再把（6）式比照（1）式進(jìn)行整理得到（4）式，因此得應(yīng)變能密度公式得證 。

如果從全微分形式 出發(fā)，因?yàn)樵瘮?shù)U（？著）的存在，可以利用牛頓-萊布尼茨公式同樣可以得到結(jié)果： 。

這說(shuō)明積分與路徑無(wú)關(guān)，所以 即

廣義格林公式 是成立的。這就從正反兩面證明了廣義格林公式，從而說(shuō)明積分與路徑無(wú)關(guān)。

5 結(jié)束語(yǔ)

對(duì)于多元函數(shù)的積分應(yīng)同時(shí)聯(lián)系實(shí)際的物理背景，不能簡(jiǎn)單的從一方面考慮，否則就會(huì)出現(xiàn)差錯(cuò)或使問(wèn)題復(fù)雜化，此題也可以從功的方面考慮。高等數(shù)學(xué)中在曲線積分中提到了格林公式，這里需要用到廣義格林公式。這兩種解法都是為了尋找應(yīng)變能密度求解中積分與路徑無(wú)關(guān)這條件。

參考文獻(xiàn)

[1]陸明萬(wàn)，羅學(xué)富.彈性理論基礎(chǔ)[M].北京：清華大學(xué)出版社，1990.

[2]徐芝綸.彈性力學(xué)簡(jiǎn)明教程[M].北京：高等教育出版社，2002.8.