數(shù)值求解雷諾方程推導(dǎo)

《數(shù)值求解雷諾方程推導(dǎo)》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、主要參數(shù)R=20mm,L=40mm,n=1000rpm,=0.3,c=2mm.各種流體潤滑問題都涉及在狹小間隙中的流體粘性流動，描寫這種物理現(xiàn)象的基本方程為雷諾方程，他的普遍形式是這個橢圓形的偏微分方程僅僅對于特殊的間隙形狀才可能求得解析解，而對于復(fù)雜的幾何形狀或者工況條件下的問題，無法用解析方法求得精確解。隨著迅速發(fā)展的點算技術(shù)，數(shù)值算法成為求解潤滑問題的有效途徑。數(shù)值法師講偏微分方程轉(zhuǎn)化為代數(shù)方程組的變換方法。它的一般原則是：首先將求解域劃分成有限個數(shù)的單元，并使每一個單元充分的微小。以至于可以認(rèn)為在各單元內(nèi)的未知量（本人畢業(yè)設(shè)計中設(shè)油膜壓力為P）相等或者依照線性變化

2、，而不會造成很大的誤差。然后，通過物理分析或數(shù)學(xué)變換方法，將求解的偏微分方程寫成離散形式，即使將它轉(zhuǎn)化成一組線性代數(shù)方程。該代數(shù)方程組表示了各個單元的待求未知量于周圍各單元未知量的關(guān)系。最后根據(jù)消去法或者迭代法求解代數(shù)方程組，從而求得整個求解域上的未知量。用來求解雷諾方程的數(shù)值方法很多，最常用的是有限元差分方法、有限元法和邊界元法，這些方法都是將求解域劃分成許多個單元，但是處理方法各不相同。在有限差分法和有限元法中，代替基本方程的函數(shù)在求解域內(nèi)是近似的，但完全滿足邊界條件。而邊界元法所用的函數(shù)在求解域內(nèi)完全滿足基本方程，但是在邊界上則近似的滿足邊界條件。一、雷諾方程的數(shù)值

3、解法根據(jù)邊界條件求解雷諾方程，這在數(shù)學(xué)上稱為邊值問題。首先將所求解的偏微分方程無量綱化。這樣做的目的是減少自變量和因變量的數(shù)目，同時用無量綱參數(shù)表示的解具有通用性。然后，將求解域劃分成等距的或者不等距的網(wǎng)格，如圖1-1為等距網(wǎng)格。圖1-1沿軸向?qū)劃分為8個等距區(qū)間，沿周向從劃分為12個等距區(qū)間。這樣在Y方向有13個節(jié)點，方向有9個節(jié)點，總計個節(jié)點。則。有限差分法如果用P代表所求的未知量例如油膜壓力，則變量P在整個域中的分布可以用各節(jié)點的P值來表示。根據(jù)差分原理，任意節(jié)點O(i,j)的一階和二階偏導(dǎo)數(shù)都可以由其周圍的節(jié)點變量值來表示。如圖1-2所示，如果采用中差分公式，則

4、變量P在O(i,j)點的偏導(dǎo)數(shù)為圖.1-2（1-1）（1-2）以P為潤滑膜壓力，雷諾方程的二維二階偏微分方程的標(biāo)準(zhǔn)形式為：（1-3）其中A,B,C,D和E都為已知量。然后將上述方程應(yīng)用到各個節(jié)點，根據(jù)中差分公式（1-1）和（1-2）用差商代替偏導(dǎo)數(shù)，即可求得各個節(jié)點的變量于相鄰各個節(jié)點變量的關(guān)系。這種關(guān)系可以寫成：（1-4）其中（1-5）式（1-4）中各系數(shù)值隨節(jié)點位置而改變。方程（1-4）是有限差分法的計算方程，對于每個節(jié)點都可以寫出一個方程，而在邊界上的節(jié)點變量應(yīng)滿足邊界條件，它們的數(shù)值是已知量。這樣，就可以求得一組線性代數(shù)方程。方程與未知量數(shù)目相一致，所以可以求解。

5、采用消去法或者迭代法求解代數(shù)方程組，并使計算結(jié)果滿足一定的收斂精度，最終求得整個求解域上各節(jié)點的變量值。求解代數(shù)方程使用迭代法求解。1、雷諾方程的無量綱化定常雷諾方程(2-1)將軸承表面沿平面展開，如圖1-1所示，并代入得等式兩邊同時乘以則雷諾方程變?yōu)?2-2)若令代入后得化簡得將代入得（2-3）由得代入（2-3）式，得再次化簡得無量綱雷諾方程（2-4）R為軸承半徑，L為軸承長度，為偏心率，為偏心距，為半徑間隙，采用有限元差分法進(jìn)行迭代計算。式（1-4）為標(biāo)準(zhǔn)形式，參考標(biāo)準(zhǔn)式（1-3）可求得標(biāo)準(zhǔn)式中A,B,C,D,E的值。將以上各值代入式（1-5）求得將已知值代入式（1-

6、4）得（2-5）將代入式（2-5）得迭代方程:將代入上式中，得（2-6）上式為最終迭代方程。邊界問題：將軸承表面沿平面展開，如圖2-1圖.2-1對于徑向軸承，方程（2-4）中兩個自變量的變化范圍是：在軸承中間斷面上Y=0：在邊緣上Y=1。而在之間變化，這一問題的邊界條件為：（1）軸向方向在邊緣Y=1處，P=0；在中間斷面Y=0上，.（2）周向方向按雷諾邊界條件：油膜起點在處，取P=0；油膜終點在發(fā)散區(qū)間內(nèi)符合P=0及的地方。