彭 震,賈培強,李 晶
(唐山學(xué)院,河北唐山 063000)
Rayleigh-Ritz法在應(yīng)用中需注意的幾個問題
彭 震,賈培強,李 晶
(唐山學(xué)院,河北唐山 063000)
結(jié)合實例提出了Rayleigh-Ritz法在應(yīng)用中需注意的幾個問題,有助于Rayleigh-Ritz法的準(zhǔn)確運用,從而提高利用該方法解決問題的計算精度。
Rayleigh-Ritz法;試函數(shù);邊界條件;勢能泛函;完備性
Rayleigh-Ritz法[1]是用最小勢能原理推導(dǎo)出的一種計算方法,在結(jié)構(gòu)分析中的重要應(yīng)用之一是利用它求解復(fù)雜問題的近似解。當(dāng)Rayleigh-Ritz法在1908年被提出時,其理論基礎(chǔ)并不完備,只是在實踐中證明它卓有成效。它的理論基礎(chǔ)在經(jīng)過人們多年努力后才開始有了頭緒。由于教材對此問題闡述的不多,在用其求解收斂性問題時又涉及到較深的泛函分析,故本文僅從該方法應(yīng)用的角度做些探討。
Rayleigh-Ritz法解決問題的出發(fā)點是將彈性結(jié)構(gòu)的總勢能 П表達(dá)成結(jié)構(gòu)位移的函數(shù)。對于一個處于穩(wěn)定平衡條件下的彈性結(jié)構(gòu),依據(jù)最小勢能原理,使結(jié)構(gòu)總勢能 П取得最小值的位移就是真實位移。位移試函數(shù)可預(yù)先假設(shè)為無窮級數(shù)、富里葉級數(shù)或任何其它形式的級數(shù)。
例如,將試函數(shù)設(shè)為無窮級數(shù)
式中 u0,v0,w0,ui,vi,wi為滿足邊界條件的假定位移函數(shù), Ai,Bi,Ci為任意常數(shù)。
對于實際問題而言,由于結(jié)構(gòu)的復(fù)雜性,事先假定的位移函數(shù)只能近似地表示結(jié)構(gòu)的真實變形形狀,基于近似形狀,利用位移函數(shù)的微分,可求反力和應(yīng)力合力的近似值。由于反力和應(yīng)力合力也是近似的,它們就不可能與結(jié)構(gòu)上的實際載荷形成靜力平衡??梢奟ayleigh-Ritz法的本質(zhì)是在維數(shù)降低后的解空間去尋找近似解。
這里需強調(diào)一點,稱 Rayleigh-Ritz法是一種近似解法(一些書籍中的稱謂),和用它求近似解是不同的??梢钥隙ǖ氖?若能列出包含真實位移在內(nèi)的所有位移,該方法亦可求得精確解。
將試函數(shù)寫成無窮級數(shù)的形式,泛函將由無窮數(shù)列 a1, a2,…,an,…來確定。根據(jù)Weierstrass定理,若試函數(shù)符合連續(xù)函數(shù)的要求,可以靠增加項數(shù)來提高精度。由于一般情況需聯(lián)立求解方程才能得出系數(shù) ai,數(shù)值條件不利于求解,所以,試函數(shù)在形式上是依賴于有限個參數(shù)an的函數(shù)族。
使用Rayleigh-Ritz法能否獲得成效,很大程度上與假設(shè)的位移函數(shù) ——試函數(shù)選取是否得當(dāng)有密切關(guān)系。構(gòu)造位移試函數(shù)最好滿足兩類邊界條件:一類是試函數(shù)須事先滿足的幾何邊界條件;另一類是由勢能泛函駐值要求導(dǎo)出的條件,變分學(xué)稱之為強制邊界條件和自然邊界條件[2]。
3.1 強制邊界條件的影響
強制邊界條件由對實際約束情況的分析給定。如圖1所示的四邊簡支矩形薄板,其上有垂直于板面的任意分布荷載作用。位移試函數(shù)可假設(shè)為雙重正弦級數(shù)
滿足w為垂直板面方向的位移的幾何邊界條件是x=0,x= a處,w=0;y=0,y=b處,w=0。
圖1 四邊簡支矩形薄板
進(jìn)一步分析還可知,在簡支邊處,限制自由翹曲,而允許繞著邊界的緣線自由轉(zhuǎn)動,故x=0,x=a處,繞x軸的彎矩Mx=0;y=0,y=b處,繞y軸的彎矩My=0。由彈性力學(xué)可知,將彎矩用位移試函數(shù)表示,上述條件也是滿足的??梢娫嚭瘮?shù)的完備性越好,就越接近真實的位移函數(shù)。能否將問題考慮得面面俱到,就目前而言還靠經(jīng)驗和深入分析。
又如圖2所示等截面懸臂梁,試函數(shù)可寫成n次多項式
也可寫成三角級數(shù)式
上兩式均滿足位移邊界條件x=0處,w=0;x=l處,w≠0。
圖2 等截面懸臂梁
3.2 自然邊界條件的影響
選擇試函數(shù)滿足自然邊界條件是不容易做到的?,F(xiàn)用下例說明自然邊界條件。圖3所示一端固定一端自由的等截面直桿,受軸向載荷q(x)=x作用,設(shè) EA=1,l=1,用試函數(shù)u(x)=a1x+a2x2來表示桿沿軸線方向位移,勢能泛函可寫成
按駐值要求δП=0導(dǎo)出,在 x=0處,有δu=0為強制
圖3 受軸向載荷的一端固定一端自由的等截面直桿
Rayleigh-Ritz法是將含有無限多變量的泛函變分問題轉(zhuǎn)換成有限個多變量的函數(shù)極值問題來處理的一種方法。若試函數(shù)滿足邊界條件,且泛函的積分表達(dá)式的積分值接近泛函的極小值,則可認(rèn)定這樣構(gòu)造的函數(shù)序列將收斂于問題的解。
考察完備性是指假設(shè)的位移函數(shù)的近似程度是否真實準(zhǔn)確地反映了結(jié)構(gòu)的變形。構(gòu)造滿足整個定義域所有邊界條件的試函數(shù)是困難的,通常可采用一些現(xiàn)有的結(jié)論。如泛函本身有二階導(dǎo)數(shù)項,選取的試函數(shù)就應(yīng)使二階導(dǎo)數(shù)存在,等等。總之,對事物認(rèn)識得越透徹越全面,越利于提高解的精度。有限單元法是Rayleigh-Ritz法局部化的解法,其核心是將解域離散,使每個單元具有完備性較好的試函數(shù),然后再在一定條件下進(jìn)行總裝。這種方法能適應(yīng)各種復(fù)雜形狀,現(xiàn)已成為行之有效的工程分析手段。
[1] 孫訓(xùn)方,方孝淑,關(guān)來泰.材料力學(xué)[M].北京:高等教育出版社,1994:249.
[2] 胡海昌,胡潤莓.變分法[M].北京:中國建筑工業(yè)出版社,1987:14.
(責(zé)任編校:夏玉玲)
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-),男,教授,主要從事固體力學(xué)的教學(xué)與研究。