基于MATLAB的圓柱軸承動力系數(shù)計算軟件的設(shè)計

周秀文 孫長年

摘 要：通過對圓柱軸承動壓油膜流場雷諾方程的推導，在MATLAB軟件中對雷諾方程采用差分法求解偏微分方程，獲取動壓油膜的相關(guān)動力學系數(shù)。

關(guān)鍵詞：圓柱軸承;動壓油膜;動力學系數(shù);雷諾方程;差分法;MATLAB

中圖分類號：O347.6 ?文獻標識碼：A ?文章編號：1671-7988（2020）16-91-05

Abstract： The Reynolds equation of dynamic oil film flow field of the journal bearing was derived， and the relevant dynamic coefficients of the oil film were obtained using finite difference method based on MATLAB.

Keywords： Journal bearing; Dynamic oil film; Dynamic Coefficient; Reynolds equation; Difference method; MATLAB

CLC NO.： O347.6 ?Document Code： A ?Article ID： 1671-7988（2020）16-91-05

引言

圓柱徑向動壓軸承廣泛應(yīng)用于汽輪機、軌道車輛等高速重載的軸系中。由于轉(zhuǎn)子系統(tǒng)的剛度受軸承油膜剛度的影響十分顯著，且轉(zhuǎn)子系統(tǒng)的阻尼很大程度上取決于軸承油膜阻尼，滑動軸承的動力特性是高速軸系轉(zhuǎn)子動力學計算的基礎(chǔ)。

本文在給定工況、軸承結(jié)構(gòu)參數(shù)、潤滑油特性等條件下，根據(jù)軸承內(nèi)油膜各點的雷諾數(shù)的大小，自動選擇層流和紊流狀態(tài)，通過聯(lián)立雷諾方程、能量方程求解擾動雷諾方程，并考慮潤滑流場的粘溫效應(yīng)，先確定靜平衡位置和靜平衡時的油膜壓力場、溫度場，并由此求得油膜力、油膜動力系數(shù)（剛度系數(shù)、阻尼系數(shù)）、摩擦力及進、出油流量和端泄流量。計算程序基于MATLAB 2010b和VC++ 6.0混合編程。

1 理論模型的建立

同理，可由Simpson數(shù)值積分求得油膜反力在x，y方向上的分力Fx，F(xiàn)y，無量綱偏位角的收斂判據(jù)為 ?，在每個迭代循環(huán)做如下修正來修正偏位角θ的值， ，其中θ的初值按經(jīng)驗公式 求得。

同樣地，摩擦力 ，進油流量 ，出油流量 ?，端泄流量 均可由復(fù)化Simpson積分求得。

4 計算程序的實現(xiàn)

為使本文的成果更便捷更好的服務(wù)于工程實際，計算軟件基于MATLAB和VC++6.0混合編程[3]，主程序由MATLAB編寫并通過mbuilder mcc打包成VC++可調(diào)用的動態(tài)鏈接庫文件（.dll），利用VC++編寫界面程序。

5 算例

某徑向動壓軸承在轉(zhuǎn)速3000RPM、軸系載重107910N、N46#透平潤滑下的計算結(jié)果（油膜壓力、溫度、黏度、密度、厚度分布以及動力學系數(shù)）如下：

溫度、密度、厚度分布以及動力學系數(shù)：

6 結(jié)語

在MATLAB中采用差分法計算了徑向動壓軸承的油膜流場偏微分方程。獲得了高速軸系轉(zhuǎn)子動力學計算必需的動壓軸承的動力學參數(shù)。此程序可快速服務(wù)于轉(zhuǎn)子動力學仿真指導高速軸系設(shè)計。

參考文獻

[1] 鐘一諤.轉(zhuǎn)子動力學[M].北京：清華大學出版社，1987.

[2] 張祖立.液體動壓徑向滑動軸承靜態(tài)特性參數(shù)計算機程序設(shè)計與計算[J].沈陽農(nóng)業(yè)大學學報，1999.

[3] 劉維.精通MATLAB與C/C++混合程序設(shè)計[M].北京：北京航空航天大學出版社，2005.