唐亮，王建梅，康建峰，馬立新，薛亞文

(太原科技大學(xué) 機(jī)械工程學(xué)院，太原 030024)

油膜軸承是一種以潤(rùn)滑油為介質(zhì)的滑動(dòng)軸承，作為重要的支承元件，廣泛應(yīng)用于中厚板軋機(jī)、冷軋機(jī)、高速線材軋機(jī)及核電發(fā)電機(jī)組等設(shè)備中[1]。動(dòng)壓油膜軸承工作時(shí)根據(jù)流體動(dòng)壓潤(rùn)滑理論，當(dāng)安裝于軸頸的錐套以一定的速度旋轉(zhuǎn)時(shí)，錐套和襯套之間被一層油膜完全隔離開，形成純液體摩擦狀態(tài)，能有效減小摩擦和磨損。因此，潤(rùn)滑油膜的工作特性決定著軸承的使用壽命[2]。

油膜軸承性能計(jì)算的可視化界面設(shè)計(jì)是對(duì)油膜軸承潤(rùn)滑理論研究的完善與擴(kuò)展。本例通過Fortran和Visual Basic(簡(jiǎn)稱VB)的混合編程得到了油膜壓力在周向和軸向的分布規(guī)律，給出了不同工況下沿軸向和周向油膜壓力變化曲線的可視化顯示。

1 可視化計(jì)算的思路

圖1所示為油膜軸承可視化計(jì)算系統(tǒng)流程?？梢暬?jì)算界面通過計(jì)算軸承的工作性能參數(shù)模擬油膜軸承的實(shí)際工況?；赩B界面輸入或更改計(jì)算參數(shù)，其中輸入?yún)?shù)包括相對(duì)間隙、軸承直徑、寬徑比和偏心率等；輸出參數(shù)包括最大承載力、最小油膜厚度、偏位角和軋制壓力，計(jì)算結(jié)果以報(bào)表形式輸出[3]。該系統(tǒng)基于Fortran源程序生成動(dòng)態(tài)鏈接庫DLL，通過VB調(diào)用動(dòng)態(tài)鏈接庫進(jìn)行混合編程，整個(gè)計(jì)算過程在后臺(tái)運(yùn)行，將計(jì)算結(jié)果返回VB主程序界面，同時(shí)繪制油膜壓力曲線。該計(jì)算系統(tǒng)關(guān)鍵在于動(dòng)態(tài)鏈接庫的建立與子程序的輸出[4]。

圖1 油膜軸承可視化計(jì)算系統(tǒng)流程

1.1 動(dòng)態(tài)鏈接庫的建立

Fortran PowerStation為建立動(dòng)態(tài)鏈接庫提供了全面的支持。建立一個(gè)動(dòng)態(tài)鏈接庫包括：生成動(dòng)態(tài)鏈接庫，輸出動(dòng)態(tài)鏈接庫中的變量或過程，以供其他程序使用，使其他程序順利地使用動(dòng)態(tài)鏈接庫[5]。

利用Fortran PowerStation建立動(dòng)態(tài)鏈接庫的步驟如下：

(1)在Microsoft FortranStation4.0的環(huán)境下，新建一個(gè)類型為Dynamic Link Library的工程并命名；

(2)在該工程中編寫Fortran源程序；

(3)編譯、建立此源程序，生成動(dòng)態(tài)鏈接庫文件。

1.2 從動(dòng)態(tài)鏈接庫中輸出程序

將Fortran編寫的計(jì)算程序?qū)雱?dòng)態(tài)鏈接庫，首先要在Fortran子程序中聲明該子程序?yàn)橥獠孔映绦颍蛊淠鼙籇LL從外部訪問[5]。Fortran程序名必須從DLL中輸出且與VB的名字完全匹配。

輸入?yún)?shù)類型，INTENT(IN)：：輸入?yún)?shù)表；

輸出參數(shù)類型，INTENT(OUT)：：輸出參數(shù)表。

1.3 動(dòng)態(tài)鏈接庫的調(diào)用

1.3.1 動(dòng)態(tài)鏈接庫DLL的聲明

Fortran動(dòng)態(tài)鏈接庫中的子程序或函數(shù)過程相對(duì)VB系統(tǒng)是個(gè)子過程，使用時(shí)須指定過程的位置和調(diào)用參數(shù)，以便編譯器找到該過程，生成正確的調(diào)用接口，聲明可用 Declare 語句提供這類信息。

如果該過程返回一個(gè)值，應(yīng)將其聲明為：Declare Function<動(dòng)態(tài)鏈接庫中的子過程名>Lib，<動(dòng)態(tài)鏈接庫中文件路徑名>[alis“過程別名”](參數(shù)表)。如果過程沒有返回值，可將其聲明為：Declare Sub<動(dòng)態(tài)鏈接庫中的子過程名>Lib，<動(dòng)態(tài)鏈接庫中文件路徑名>[alis“過程別名”](參數(shù)表)。

1.3.2 動(dòng)態(tài)鏈接庫DLL的調(diào)用

在聲明DLL 過程之后，VB將其作為內(nèi)部過程使用。在聲明函數(shù)之后，其調(diào)用方式與標(biāo)準(zhǔn)的 VB函數(shù)相同。需要先定義被調(diào)用子程序的過程參數(shù)類型，然后賦值，最后使用Call語句調(diào)用。調(diào)用格式為：Call子過程或函數(shù)過程名(實(shí)參，…)，其中無參數(shù)時(shí)，省略括號(hào)[6]。

2 可視化界面的開發(fā)

2.1 油膜壓力計(jì)算

油膜壓力計(jì)算關(guān)鍵在于Reynolds方程的求解。Reynolds方程是流體動(dòng)壓潤(rùn)滑基本方程，普遍形式為[7]

(1)

對(duì)于穩(wěn)態(tài)工況的軸承，假設(shè)潤(rùn)滑油為不可壓縮流體，液體潤(rùn)滑劑密度可視為常數(shù)。當(dāng)潤(rùn)滑油膜的熱效應(yīng)不十分明顯時(shí)，可視為等溫狀態(tài)，即流體黏度在整個(gè)潤(rùn)滑膜中保持不變，可將Reynolds方程簡(jiǎn)化為

(2)

為計(jì)算簡(jiǎn)便，對(duì)Reynolds方程進(jìn)行如下無量綱化[6]

其中，x為油膜軸承周向長(zhǎng)度；z為油膜軸承軸向長(zhǎng)度；y為油膜厚度；R為軸承半徑；c為半徑間隙；L為軸承寬度；ε為偏心率；ω為軸頸自轉(zhuǎn)角速度；η0為入口處潤(rùn)滑油的黏度。

(3)

運(yùn)用有限差分法解上述偏微分方程，利用有限差分近似值代替導(dǎo)數(shù)，把二維偏微分方程轉(zhuǎn)換為一維差分方程，將承載區(qū)沿周向展開成矩形區(qū)域并劃分網(wǎng)格，計(jì)算各個(gè)節(jié)點(diǎn)的壓力值[8]。

2.2 油膜壓力分布的可視化

將油膜軸承承載區(qū)沿周向平均分為M個(gè)節(jié)點(diǎn)，沿軸向平均分為N個(gè)節(jié)點(diǎn)，本例設(shè)定M=49，N=37，運(yùn)用有限差分法計(jì)算出各網(wǎng)格節(jié)點(diǎn)的壓力值，然后采用VB繪圖，其編程方法為：

(1)在VB界面中創(chuàng)建命令控件1(commandbutton)，標(biāo)簽設(shè)定為“繪制壓力曲線”；創(chuàng)建命令控件2，標(biāo)簽設(shè)定為“清屏”，代碼中使用Cls方法清屏[7]。

(2)在VB界面中創(chuàng)建圖片框控件(picturebox)，在控件上輸出圖形并用print方法輸出文本，同時(shí)設(shè)定繪圖區(qū)域大小，其代碼在命令控件1中設(shè)置。

(3)用scale方法設(shè)定坐標(biāo)區(qū)域大小，使用line方法繪制兩坐標(biāo)軸，橫坐標(biāo)為周向或軸向節(jié)點(diǎn)，縱坐標(biāo)為無量綱壓力大小，用print方法輸出坐標(biāo)軸信息，其代碼在命令控件1中設(shè)置。

(4)創(chuàng)建文本框用以輸入所繪制的位置，并在命令控件1中編寫代碼給變量賦值，在循環(huán)語句中寫入line語句，從而繪制出壓力分布曲線[7]。

輸入表1中工況1的參數(shù)，點(diǎn)擊VB界面“計(jì)算”按鈕，并點(diǎn)擊“繪制壓力曲線”按鈕，繪制出承載區(qū)周向壓力分布(圖2)和軸向壓力分布(圖3)；輸入工況2的參數(shù)，進(jìn)行與工況1相同的操作，繪制出承載區(qū)周向壓力分布(圖4)和軸向壓力分布(圖5)。結(jié)合工況1和工況2的計(jì)算結(jié)果，觀察可知油膜軸承的受載特征。

表1 油膜軸承輸入?yún)?shù)

圖2 工況1的周向油膜壓力分布

圖3 工況1的軸向油膜壓力分布

圖4 工況2的周向油膜壓力分布

圖5 工況2的軸向油膜壓力分布

(1)選定工況1中周向第19組數(shù)據(jù)進(jìn)行繪圖(圖2)，周向油膜壓力分布具有對(duì)稱性，在0～15節(jié)點(diǎn)和35～48節(jié)點(diǎn)內(nèi)數(shù)值較小,而且變化平穩(wěn)；在15～35節(jié)點(diǎn)內(nèi)數(shù)值較大；在27節(jié)點(diǎn)附近達(dá)到峰值，確定該點(diǎn)區(qū)域?yàn)橹芟虺惺苘堉屏Φ闹饕獏^(qū)域。

(2)選定工況1中軸向第30組數(shù)據(jù)進(jìn)行繪圖(圖3)，軸向油膜壓力分布均勻且對(duì)稱，在19節(jié)點(diǎn)附近達(dá)到峰值，確定該點(diǎn)區(qū)域是軸向承受軋制力的主要區(qū)域。

(3)分別選定周向第19組數(shù)據(jù)和軸向第30組數(shù)據(jù)進(jìn)行繪圖(圖4、圖5)，油膜壓力分布趨勢(shì)基本相同。對(duì)比圖4與圖2、圖5與圖3可知，相同節(jié)點(diǎn)處在工況2下的油膜壓力高于工況1。由此可知，在其他輸入?yún)?shù)不變的情況下，油膜壓力隨著軋輥轉(zhuǎn)速的升高而增大。

3 結(jié)束語

本系統(tǒng)可依據(jù)不同參數(shù)進(jìn)行計(jì)算和繪圖，可以直觀地顯示油膜壓力沿周向和軸向的分布規(guī)律，能方便地分析出不同工況對(duì)油膜壓力的影響。將可視化技術(shù)作為油膜軸承潤(rùn)滑理論研究的新方法，將會(huì)顯著地縮短油膜軸承的設(shè)計(jì)周期，有效降低設(shè)計(jì)成本。