劉宗民，梁立孚，樊濤

(哈爾濱工程大學(xué)航天與建筑工程學(xué)院，黑龍江哈爾濱150001)

基于基面力的非保守系統(tǒng)有限變形問題是一個(gè)新興的課題?；媪ψ鳛橐环N描述應(yīng)力狀態(tài)的新方法，較傳統(tǒng)的應(yīng)力張量表示方法簡單[1]。利用基面力，各種應(yīng)力張量都可以由表征位移的協(xié)變矢基的并矢來表示，可以得到彈性力學(xué)基本方程(平衡方程、邊界條件、本構(gòu)關(guān)系)的簡潔表達(dá)式?；媪哂袀鹘y(tǒng)的二階應(yīng)力張量無法比擬的優(yōu)越性，尤其是將基面力應(yīng)用于有限變形問題，這種特點(diǎn)就顯得十分明顯。有限變形的勢能原理可以和小變形情況類似的建立起來。但是，有限變形的余能原理經(jīng)歷了近百年，卻未給出適當(dāng)?shù)谋磉_(dá)形式。高玉臣利用基面力作為基本未知量，給出了保守系統(tǒng)有限變形余能原理的適當(dāng)表達(dá)形式[2]。

對(duì)于非保守系統(tǒng)，國外以Leipholz為代表，提出廣義自共軛的概念，建立了廣義的Hamilton原理，給出了著名的Leipholz桿模型[3]。我國學(xué)者也作了大量工作，提出并完善了各類擬變分原理。劉殿魁提出了彈性非保守系統(tǒng)的一般擬變分原理[4]。梁立孚、宋海燕研究了非保守系統(tǒng)彈性靜(動(dòng))力學(xué)兩類變量的廣義擬變分原理[5-6]。

對(duì)于有限變形問題，系統(tǒng)也常常是非保守的。文獻(xiàn)[7]對(duì)非線性非保守系統(tǒng)彈性力學(xué)擬變分原理作了系統(tǒng)的研究。但是，由于有限變形中“裂隙函數(shù)”[8]的存在，使得對(duì)于有限變形非保守系統(tǒng)擬變分原理的研究變得困難。本文以基面力和位移梯度為基本未知量，“裂隙函數(shù)”自動(dòng)消失了，采用變積法建立了類似于小變形情況的彈性動(dòng)力學(xué)時(shí)域邊值問題的有限變形非保守系統(tǒng)三類變量的廣義擬Hamilton原理。本文研究的非保守系統(tǒng)專指“伴生力”系統(tǒng)，即作用于系統(tǒng)的非保守力隨物體的變形而變化。

1 基本方程

按照文獻(xiàn)[9-11]，基于基面力的有限變形時(shí)域邊值問題的基本方程為平衡方程:

幾何方程:

本構(gòu)關(guān)系:

時(shí)域邊值條件:

式中:Ti為基面力，gi為位移梯度，u為位移，ū為給定位移，mi為變形前面元單位法矢量分量，τ-0為變形前應(yīng)力邊界上的面力，f為單位質(zhì)量物體所受的體力，R為動(dòng)量矢量，V為速度矢量，ρ0為變形前的密度，VP為變形前的基容，A為單位質(zhì)量的應(yīng)變能，B為單位質(zhì)量的余能，Sσ為給定面力的邊界，Su為給定位移的邊界，ū0為t0時(shí)刻位移，ū1為t1時(shí)刻位移,為對(duì)時(shí)間的導(dǎo)數(shù)。

2 三類變量的廣義擬Hamilton原理

應(yīng)用變積方法[12-13]，根據(jù)廣義力和廣義位移之間的對(duì)應(yīng)關(guān)系，將式(1)～(6)和式(8)乘上相應(yīng)的虛量，然后積分并代數(shù)相加，再對(duì)時(shí)間域積分，可得

應(yīng)用Green定理:

及分部積分公式:

將式(13)和(14)代入式(12)，并按慣例在時(shí)域邊界處取δu=0，經(jīng)整理后可將式(12)變換為

上式可進(jìn)一步表示為

將式(16)簡記為

式中:

式(17)稱為基于基面力的非保守系統(tǒng)有限變形三類變量廣義擬Hamilton原理。

如果應(yīng)用Green定理:

及分部積分公式:

將式(20)和(21)代入式(12)，并將式中各項(xiàng)改變符號(hào)，并按慣例在時(shí)域邊界處取δR=0，經(jīng)整理后可將式(12)變換為

上式可進(jìn)一步表示為

將式(23)簡記為

式中:

式(24)稱為基于基面力的非保守系統(tǒng)有限變形三類變量廣義擬余Hamilton原理。

由以上的推導(dǎo)過程可見，泛函 δΓH3+δQH+δPH=0是 δΠH3- δQH- δPH=0的對(duì)偶形式[14]。

3 正確性的驗(yàn)證

按照基本變量之間的對(duì)應(yīng)關(guān)系，將基于基面力的有限變形廣義擬Hamilton原理轉(zhuǎn)化為以第二類P-K應(yīng)力張量和Green應(yīng)變張量為基本變量的有限變形廣義擬Hamilton原理。進(jìn)而證明本文建立的基于基面力的有限變形廣義擬Hamilton原理的正確性。

基面力Ti與第二類P-K應(yīng)力張量Σ，以及位移梯度gi與Green應(yīng)變張量εG之間的對(duì)應(yīng)關(guān)系可以表示為[9-10]

圖1表示變形前與變形后的位置關(guān)系。如圖1所示，位移u和位移梯度gi可以表示為式中:F為變形梯度，P和Q分別代表一點(diǎn)變形前和變形后的矢徑，Pi和Qi分別為變形前和變形后的協(xié)變基矢，Pi和Qi分別為變形前和變形后的逆變基矢。

圖1 位移和變形Fig.1 Displacement and deformation

將位移梯度gi表示成在標(biāo)架Pj中的分量形式，可得

式中:tij為基面力Ti在標(biāo)架Qj上的分量.

將式(29)～(31)代入式(32)，可得

考慮式(27)、(28)、(31)和(33)，式(25)中等號(hào)右端第4項(xiàng)的變分式、以及等號(hào)右端第5項(xiàng)可以分別表示為[11，15-16]

式中:Σij為第二類P-K應(yīng)力張量Σ的分量，Eij為Green應(yīng)變張量εG的分量.

考慮式(27)，(28)，(33)，(34)和(35)，式(17)和(24)可以分別改寫為

式中:

式(36)為以第二類P-K應(yīng)力張量和Green應(yīng)變張量為基本變量非保守系統(tǒng)有限變形三類變量廣義擬Hamilton原理。

式中:

式(37)為以第二類P-K應(yīng)力張量和Green應(yīng)變張量為基本變量非保守系統(tǒng)有限變形三類變量廣義擬余Hamilton原理。式中，vk，rk，F(xiàn)-k和P-k分別表示為速度分量，動(dòng)量分量，體力分量和面力分量。

式(36)和(37)的表達(dá)式與文獻(xiàn)[7]中的結(jié)果是一致的，從而證明本文建立的基于基面力的有限變形廣義擬Hamilton原理的正確性。

4 結(jié)束語

以基面力和位移梯度為基本變量，應(yīng)用變積法建立了有限變形非保守系統(tǒng)三類變量的廣義擬Hamilton原理。當(dāng)把基面力和位移梯度作為基本變量時(shí)，“裂隙函數(shù)”自動(dòng)消失了，繁冗的公式得以簡化，建立有限變形的擬變分原理就變得比較容易。本文建立的有限變形的廣義擬Hamilton原理，具有和小變形情況一樣的結(jié)構(gòu)形式。本文的工作表明，應(yīng)用基面力理論解決有限變形問題，接近于彈性小變形問題的理論和方法，這不僅使得有限變形理論和應(yīng)用變得較為容易，而且對(duì)借鑒較為成熟的彈性小變形問題的理論和方法來解決有限變形的理論和應(yīng)用問題提供了方便。

[1]GAO Yuchen.A new description of the stress state at a point with applications[J].Archive of Applied Mechanics，2003，73:171-183.

[2]GAO Yuchen.Complementary energy principle for large elastic deformation[J].Science in China:Series G，2006，49(3):341-356.

[3]LEIPHOLZ H.Direct variational methods and eigenvalue problems in engineering[M].Leyden:Noordhoff International Publishing，1977:207-217.

[4]劉殿魁，張其浩.彈性理論中非保守問題的一般變分原理[J].力學(xué)學(xué)報(bào)，1981(6):562-570.LIU Diankui，ZHANG Qihao.Some general variational principles for non-conservative problems in theory of elasticity[J].Acta Mechanica Sinica，1981(6):562-570.

[5]LIANG Lifu，LIU Diankui，SONG Haiyan.The generalized quasi-variational principles of non-conservative systems with two kinds of variables[J].Science in China:Series G，2005，48(5):600-613.

[6]宋海燕，梁立孚，劉殿魁.非保守彈性動(dòng)力系統(tǒng)兩類變量的廣義擬變分原理及其應(yīng)用[J].地震工程與工程振動(dòng)，2005，25(4):24-30.SONG Haiyan，LIANG Lifu，LIU Diankui.On generalized quasi-variational principles of two-kinds of variables in nonconservative elastodynamic systems and their application[J].Journal of Earthquake Engineering and Engineering Vibration，2005，25(4):24-30.

[7]樊濤.非線性非保守系統(tǒng)彈性力學(xué)擬變分原理研究[D].哈爾濱:哈爾濱工程大學(xué)，2007:57-91.FAN Tao.Research on the quasi-variational principle in nonlinear non-conservative elasticity[D].Harbin:Harbin Engineering University，2007:57-91.

[8]GAO D Y，STRANG G.Geometric nonlinearity:potential energy，complementary energy and the gap function[J].Quarterly of Applied Mathematics，1989，47(3):487-504.

[9]高玉臣.固體力學(xué)基礎(chǔ)[M].北京:中國鐵道出版社，1999:14-60.

[10]金明.非線性連續(xù)介質(zhì)力學(xué)教程[M].北京:清華大學(xué)出版社，北京交通大學(xué)出版社聯(lián)合出版，2005:80-146.

[11]劉宗民.基于基面力的彈性大變形擬變分原理[D].哈爾濱:哈爾濱工程大學(xué)，2008:21-152.LIU Zongmin.Quasi-variational principles for large elastic deformation based on base forces[D].Harbin:Harbin Engineering University，2008:21-152.

[12]梁立孚.變分原理及其應(yīng)用[M].哈爾濱:哈爾濱工程大學(xué)出版社等五社聯(lián)合出版，2005:97-115.

[13]LIANG Lifu，HU Haichang.Generalized variational principles of three kinds of variables in general mechanics[J].Science in China:Series A，2001，44(6):770-776.

[14]宋海燕，梁立孚，周軼雷.一般力學(xué)廣義變分原理的對(duì)偶形式[J].哈爾濱工程大學(xué)學(xué)報(bào)，2004，25(6):740-743，760.SONG Haiyan，LIANG Lifu，ZHOU Yilei.Dual form of generalized variational principles in general mechanics[J].Journal of Harbin Engineering University，2004，25(6):740-743，760.

[15]LIU Zongmin，LIANG Lifu，F(xiàn)AN Tao.Quasi-variational principles of large elastic deformation in non-conservative systems based on base forces theory and its application[C]//Proceedings of 2008 ASME International Mechanical Engineering Congress and Exposition.Boston，USA，2008，12:487-493.

[16]LIU Zongmin，WU Guohui，LIANG Lifu，et al.Geometric nonlinear elasto-dynamics problems solved by Variational method[C]//Proceedings of 2010 ASME International Mechanical Engineering Congress and Exposition.Vancouver，Canada，2010，8:299-305.