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

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

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

4、（1-1） （1-2） 以P為潤滑膜壓力，雷諾方程的二維二階偏微分方程的標準形式為： （1-3）其中A,B,C,D和E都為已知量。然后將上述方程應用到各個節點，根據中差分公式（1-1）和（1-2）用差商代替偏導數，即可求得各個節點的變量于相鄰各個節點變量的關系。這種關系可以寫成： （1-4）其中 （1-5）式（1-4）中各系數值隨節點位置而改變。方程（1-4）是有限差分法的計算方程，對于每個節點都可以寫出一個方程，而在邊界上的節點變量應滿足邊界條件，它們的數值是已知量。這樣，就可以求得一組線性代數方程。方程與未知量數目相一致，所以可以求解。采用消去法或者迭代法求解代數方程組，并使計算結果滿足一

5、定的收斂精度，最終求得整個求解域上各節點的變量值。求解代數方程使用迭代法求解。1、 雷諾方程的無量綱化定常雷諾方程 (2-1)將軸承表面沿平面展開，如圖1-1所示，并代入得等式兩邊同時乘以則雷諾方程變為 (2-2)若令代入后得化簡得將 代入得 （2-3）由得代入（2-3）式，得 再次化簡得無量綱雷諾方程 （2-4）R為軸承半徑，L為軸承長度，為偏心率，為偏心距，為半徑間隙，采用有限元差分法進行迭代計算。式（1-4）為標準形式，參考標準式（1-3）可求得標準式中A,B,C,D,E的值。將以上各值代入式（1-5）求得將已知值代入式（1-4）得 （2-5）將代入式（2-5）得迭代方程:將代入上式中，得 （2-6）上式為最終迭代方程。邊界問題：將軸承表面沿平面展開，如圖2-1圖.2-1對于徑向軸承，方程（2-4）中兩個自變量的變化范圍是：在軸承中間斷面上Y=0：在邊緣上 Y=1。而在之間變化，這一問題的邊