彭 震,賈培強,李 晶

(唐山學院,河北唐山 063000)

Rayleigh-Ritz法在應用中需注意的幾個問題

彭 震,賈培強,李 晶

(唐山學院,河北唐山 063000)

結(jié)合實例提出了Rayleigh-Ritz法在應用中需注意的幾個問題,有助于Rayleigh-Ritz法的準確運用,從而提高利用該方法解決問題的計算精度。

Rayleigh-Ritz法;試函數(shù);邊界條件;勢能泛函;完備性

Rayleigh-Ritz法[1]是用最小勢能原理推導出的一種計算方法,在結(jié)構(gòu)分析中的重要應用之一是利用它求解復雜問題的近似解。當Rayleigh-Ritz法在1908年被提出時,其理論基礎并不完備,只是在實踐中證明它卓有成效。它的理論基礎在經(jīng)過人們多年努力后才開始有了頭緒。由于教材對此問題闡述的不多,在用其求解收斂性問題時又涉及到較深的泛函分析,故本文僅從該方法應用的角度做些探討。

1 Rayleigh-Ritz法的本質(zhì)

Rayleigh-Ritz法解決問題的出發(fā)點是將彈性結(jié)構(gòu)的總勢能 П表達成結(jié)構(gòu)位移的函數(shù)。對于一個處于穩(wěn)定平衡條件下的彈性結(jié)構(gòu),依據(jù)最小勢能原理,使結(jié)構(gòu)總勢能 П取得最小值的位移就是真實位移。位移試函數(shù)可預先假設為無窮級數(shù)、富里葉級數(shù)或任何其它形式的級數(shù)。

例如,將試函數(shù)設為無窮級數(shù)

式中 u0,v0,w0,ui,vi,wi為滿足邊界條件的假定位移函數(shù), Ai,Bi,Ci為任意常數(shù)。

對于實際問題而言,由于結(jié)構(gòu)的復雜性,事先假定的位移函數(shù)只能近似地表示結(jié)構(gòu)的真實變形形狀,基于近似形狀,利用位移函數(shù)的微分,可求反力和應力合力的近似值。由于反力和應力合力也是近似的,它們就不可能與結(jié)構(gòu)上的實際載荷形成靜力平衡?？梢奟ayleigh-Ritz法的本質(zhì)是在維數(shù)降低后的解空間去尋找近似解。

這里需強調(diào)一點,稱 Rayleigh-Ritz法是一種近似解法(一些書籍中的稱謂),和用它求近似解是不同的?？梢钥隙ǖ氖?若能列出包含真實位移在內(nèi)的所有位移,該方法亦可求得精確解。

2 多項式的項數(shù)對試函數(shù)的影響

將試函數(shù)寫成無窮級數(shù)的形式,泛函將由無窮數(shù)列 a1, a2,…,an,…來確定。根據(jù)Weierstrass定理,若試函數(shù)符合連續(xù)函數(shù)的要求,可以靠增加項數(shù)來提高精度。由于一般情況需聯(lián)立求解方程才能得出系數(shù) ai,數(shù)值條件不利于求解,所以,試函數(shù)在形式上是依賴于有限個參數(shù)an的函數(shù)族。

3 兩類邊界條件對試函數(shù)的影響

使用Rayleigh-Ritz法能否獲得成效,很大程度上與假設的位移函數(shù) ——試函數(shù)選取是否得當有密切關(guān)系。構(gòu)造位移試函數(shù)最好滿足兩類邊界條件:一類是試函數(shù)須事先滿足的幾何邊界條件;另一類是由勢能泛函駐值要求導出的條件,變分學稱之為強制邊界條件和自然邊界條件[2]。

3.1 強制邊界條件的影響

強制邊界條件由對實際約束情況的分析給定。如圖1所示的四邊簡支矩形薄板,其上有垂直于板面的任意分布荷載作用。位移試函數(shù)可假設為雙重正弦級數(shù)

滿足w為垂直板面方向的位移的幾何邊界條件是x=0,x= a處,w=0;y=0,y=b處,w=0。

圖1 四邊簡支矩形薄板

進一步分析還可知,在簡支邊處,限制自由翹曲,而允許繞著邊界的緣線自由轉(zhuǎn)動,故x=0,x=a處,繞x軸的彎矩Mx=0;y=0,y=b處,繞y軸的彎矩My=0。由彈性力學可知,將彎矩用位移試函數(shù)表示,上述條件也是滿足的?？梢娫嚭瘮?shù)的完備性越好,就越接近真實的位移函數(shù)。能否將問題考慮得面面俱到,就目前而言還靠經(jīng)驗和深入分析。

又如圖2所示等截面懸臂梁,試函數(shù)可寫成n次多項式

也可寫成三角級數(shù)式

上兩式均滿足位移邊界條件x=0處,w=0;x=l處,w≠0。

圖2 等截面懸臂梁

3.2 自然邊界條件的影響

選擇試函數(shù)滿足自然邊界條件是不容易做到的。現(xiàn)用下例說明自然邊界條件。圖3所示一端固定一端自由的等截面直桿,受軸向載荷q(x)=x作用,設 EA=1,l=1,用試函數(shù)u(x)=a1x+a2x2來表示桿沿軸線方向位移,勢能泛函可寫成

按駐值要求δП=0導出,在 x=0處,有δu=0為強制

圖3 受軸向載荷的一端固定一端自由的等截面直桿

4 完備性對試函數(shù)的影響

Rayleigh-Ritz法是將含有無限多變量的泛函變分問題轉(zhuǎn)換成有限個多變量的函數(shù)極值問題來處理的一種方法。若試函數(shù)滿足邊界條件,且泛函的積分表達式的積分值接近泛函的極小值,則可認定這樣構(gòu)造的函數(shù)序列將收斂于問題的解。

考察完備性是指假設的位移函數(shù)的近似程度是否真實準確地反映了結(jié)構(gòu)的變形。構(gòu)造滿足整個定義域所有邊界條件的試函數(shù)是困難的,通?？刹捎靡恍┈F(xiàn)有的結(jié)論。如泛函本身有二階導數(shù)項,選取的試函數(shù)就應使二階導數(shù)存在,等等?？傊?對事物認識得越透徹越全面,越利于提高解的精度。有限單元法是Rayleigh-Ritz法局部化的解法,其核心是將解域離散,使每個單元具有完備性較好的試函數(shù),然后再在一定條件下進行總裝。這種方法能適應各種復雜形狀,現(xiàn)已成為行之有效的工程分析手段。

[1] 孫訓方,方孝淑,關(guān)來泰.材料力學[M].北京:高等教育出版社,1994:249.

[2] 胡海昌,胡潤莓.變分法[M].北京:中國建筑工業(yè)出版社,1987:14.

(責任編校:夏玉玲)

Questions in the Application of Rayleigh-Ritz Method

PENG Zhen,JIA Pei-qiang,LIJing
(Tangshan College,Tangshan 063000,China)

Combined with examples,several questions w ere pointed out in the apllication of Ray-leigh-Ritz method.Solving these problems can contribute to the right apllication of Rayleigh-Ritz method and also imp rove the calculation accuracy.

Rayleigh-Ritz method;trial function;boundary condition;potential energy fonctio-nelle;completeness

O177

A

1672-349X(2010)06-0008-02

2010-07-14

彭震(1955-),男,教授,主要從事固體力學的教學與研究。