吳超，尹雪梅，崔曉康，羅國富

(鄭州輕工業(yè)學(xué)院機電工程學(xué)院，河南鄭州450002)

隨著工業(yè)技術(shù)的發(fā)展，對旋轉(zhuǎn)機械的要求越來越高，高速、高精度、高負(fù)荷一直是旋轉(zhuǎn)機械技術(shù)發(fā)展追求的目標(biāo)，而油膜軸承－轉(zhuǎn)子系統(tǒng)作為大型旋轉(zhuǎn)機械的“心臟”，其性能參數(shù)的計算、設(shè)計至關(guān)重要。經(jīng)典的研究方法是:通過對軸承油膜進行線性化處理，用油膜剛度和阻尼來表達油膜的動力特性[1]。由于油膜軸承的動力特性求解比較復(fù)雜，很難得到其解析解，常利用計算機編程計算得到其離散值[2]。因Fortran 語言具有強大的計算功能，傳統(tǒng)的油膜軸承－轉(zhuǎn)子系統(tǒng)計算設(shè)計多用Fortran 語言編寫程序，已經(jīng)積累了大量成熟的計算滑動軸承－轉(zhuǎn)子系統(tǒng)性能參數(shù)的程序[2－3]。由于Fortran 語言的DOS 字符界面與現(xiàn)有的開發(fā)工具兼容性不好，因而需要對Fortran油膜軸承計算程序加以升級改造以滿足現(xiàn)代程序開發(fā)的要求。一些學(xué)者為使Fortran 油膜軸承性能計算程序能夠適應(yīng)時代的發(fā)展，對其進行了可視化研究。作者對油膜軸承－轉(zhuǎn)子系統(tǒng)動力學(xué)及其可視化進行綜述并展望其發(fā)展。

1 油膜軸承－轉(zhuǎn)子系統(tǒng)動力學(xué)計算綜述

油膜軸承是靠一層薄的油膜來減少軸和軸承之間的摩擦和磨損的，其承載性能和油膜的特性息息相關(guān)。經(jīng)典的計算方法是對Reynolds 方程、油膜厚度方程和熱力學(xué)方程聯(lián)合求解得到壓力場和溫度場，然后進行積分得到軸承的線性動、靜特性系數(shù)。幾十年前，研究者在對擠壓油膜阻尼器[4]、可控油膜軸承[5－7]的研究中發(fā)現(xiàn)了大量的非線性現(xiàn)象，并注意到對油膜進行線性化處理的結(jié)果不能解釋這些現(xiàn)象，二者之間存在著本質(zhì)的差別。針對油膜的非線性問題，國內(nèi)外研究者提出了短軸承理論[8－10]、油膜力數(shù)據(jù)庫[11]等計算模型，能很好地解決軸承油膜的非線性問題。焦映厚等[12]分析了各種非線性油膜力模型后指出，當(dāng)綜合考慮動力學(xué)特性、計算精度和計算速度等方面的因素時，油膜力數(shù)據(jù)庫方法應(yīng)為優(yōu)先采用的模型。

把油膜力模型納入到軸承－轉(zhuǎn)子系統(tǒng)的非線性動力學(xué)計算中，一些學(xué)者發(fā)展了計算油膜軸承－轉(zhuǎn)子系統(tǒng)非線性動力學(xué)的方法。求解油膜軸承－轉(zhuǎn)子系統(tǒng)動力學(xué)方程最直接的方法是數(shù)值積分。但直接積分不便于研究系統(tǒng)動力學(xué)內(nèi)部的分岔及穩(wěn)定性規(guī)律。為方便研究軸承－轉(zhuǎn)子系統(tǒng)行為規(guī)律，提高計算效率，ADILETTA 等[10]、袁 小 陽 等[13]、SUNDARARAJAN等[14]、李立等人[15]提出了改進的算法，并把這些算法運用到對軸承－轉(zhuǎn)子系統(tǒng)的穩(wěn)定性和非線性動力學(xué)特性的研究中。這些算法在進行時間積分獲得軸心軌跡時，必須計算多次非線性油膜力，消耗大量的時間。李志剛等[16]提出利用非線性油膜力數(shù)據(jù)庫方法來求解滑動軸承的油膜力，解決了耗時問題。并在此基礎(chǔ)上提出了一套新的關(guān)于多自由度的轉(zhuǎn)子－軸承系統(tǒng)的非線性動力學(xué)仿真計算方法，詳細考察了撓性轉(zhuǎn)子－軸承系統(tǒng)不平衡質(zhì)量大小和分布對轉(zhuǎn)子系統(tǒng)穩(wěn)定性的影響。文獻[17－19]是在非線性油膜力數(shù)據(jù)庫的基礎(chǔ)上，提出了轉(zhuǎn)子系統(tǒng)動力學(xué)的計算新方法。以上計算方法的出現(xiàn)為研究油膜軸承－轉(zhuǎn)子系統(tǒng)的可視化提供了強有力的技術(shù)支持。

2 油膜軸承－轉(zhuǎn)子系統(tǒng)動力學(xué)可視化研究綜述

在進行油膜軸承工程計算的過程中，積累了大量Fortran 語言程序，由于Fortran 語言與現(xiàn)有的開發(fā)工具兼容性不好，使用很不方便。為使Fortran 油膜軸承計算程序能夠適應(yīng)時代的發(fā)展，發(fā)揮其長期理論指導(dǎo)意義，因而有必要對原有和現(xiàn)有Fortran 計算程序加以升級改造以滿足現(xiàn)代程序開發(fā)的需要，可視化的開發(fā)手段無疑是延續(xù)油膜軸承性能計算程序生命力的唯一途徑。

崔曉康等[20]基于VB 和Fortran 語言混合編程的方法，研究滑動軸承的壓力分布和溫度分布。王建梅等[21]對比了VB 調(diào)用Fortran程序的3種方法，說明了使用Fortran 動態(tài)鏈接庫的優(yōu)點，并使用MDI 窗體優(yōu)化了VB 界面，實現(xiàn)了油膜軸承可視化的系統(tǒng)性能計算的應(yīng)用程序。陽旭等人[22]開發(fā)了基于OpenGL 交互控制的滑動軸承油膜壓力動態(tài)分布三維可視化仿真系統(tǒng)，能實時顯示任意轉(zhuǎn)速、潤滑油黏度、寬徑比等參數(shù)下的油膜壓力動態(tài)分布三維模型。吳剛等人[23]基于短軸承理論，對滑動軸承軸心軌跡的可視化進行了研究。范紅紅等[24]基于OpenGL技術(shù)實現(xiàn)了油膜特性可視化系統(tǒng)的開發(fā)過程，包括油膜厚度、油膜壓力、油膜溫度的分布規(guī)律。郅剛鎖等[25]基于粒子紋理映射的方法，建立了虛擬軸承流場分析系統(tǒng)，可以更深入地研究軸承運行過程中油膜的動態(tài)變化和傳熱機理。王丹丹等[26]基于VC/MFC，采用面向?qū)ο蠹夹g(shù)，開發(fā)了常見結(jié)構(gòu)動靜壓滑動軸承設(shè)計計算軟件包，可以方便地計算不同結(jié)構(gòu)參數(shù)滑動軸承的靜、動特性系數(shù)并進行結(jié)果分析。這些可視化軟件包的開發(fā)，提高了油膜軸承的設(shè)計效率和設(shè)計質(zhì)量，為油膜軸承特性的研究提供了新的手段，同時豐富了可視化技術(shù)的應(yīng)用領(lǐng)域。

3 油膜軸承－轉(zhuǎn)子系統(tǒng)動力學(xué)可視化研究展望

目前，油膜軸承－轉(zhuǎn)子系統(tǒng)的可視化研究以油膜軸承的性能研究居多，通過這些可視化軟件能夠得到軸承的壓力場、溫度場以及軸承的動、靜特性參數(shù)，可以方便地進行軸承結(jié)構(gòu)參數(shù)的設(shè)計，但一直沒有涉及到轉(zhuǎn)子系統(tǒng)的設(shè)計。隨著旋轉(zhuǎn)機械向著高速高精度的方向發(fā)展，要想提高旋轉(zhuǎn)機械的精度，必須考慮油膜軸承的非線性，并把油膜軸承－轉(zhuǎn)子系統(tǒng)作為整體來進行研究，也更符合轉(zhuǎn)子－軸承系統(tǒng)的實際運行情況。因此，油膜軸承－轉(zhuǎn)子系統(tǒng)可視化研究，尤其在工業(yè)現(xiàn)場要想有很好的應(yīng)用，仍需做如下幾個方面的工作:

(1)納入油膜軸承－轉(zhuǎn)子系統(tǒng)的非線性模型，開發(fā)適合于轉(zhuǎn)子－軸承系統(tǒng)設(shè)計計算的三維可視化軟件系統(tǒng)，使軸承特性參數(shù)和轉(zhuǎn)子的振動、運行軌跡等計算結(jié)果在可視化界面顯示，便于現(xiàn)場人員時時了解轉(zhuǎn)子的運行狀態(tài);

(2)基于全局的優(yōu)化設(shè)計，除了讓軸承－轉(zhuǎn)子系統(tǒng)滿足相應(yīng)的機械要求之外，還應(yīng)從系統(tǒng)的角度考慮油膜軸承－轉(zhuǎn)子系統(tǒng)的穩(wěn)定性和可靠性，給出系統(tǒng)的穩(wěn)定裕度;

(3)考慮擾動等外界因素對轉(zhuǎn)子－軸承系統(tǒng)穩(wěn)定性和振動的影響，開發(fā)更人性化的可視化界面，便于操作和交流。

4 結(jié)束語

油膜軸承－轉(zhuǎn)子系統(tǒng)可視化軟件能夠以圖形的方式顯示油膜軸承－轉(zhuǎn)子系統(tǒng)的壓力場、溫度場、軸心軌跡、穩(wěn)定裕度以及其他特性參數(shù)，能指導(dǎo)設(shè)計人員預(yù)先了解軸承－轉(zhuǎn)子系統(tǒng)的運行狀態(tài)，使軸承－轉(zhuǎn)子系統(tǒng)的設(shè)計過程可視化、智能化、程序化，也方便對轉(zhuǎn)子－軸承系統(tǒng)的設(shè)計、改進，從而提高滑動軸承－轉(zhuǎn)子系統(tǒng)以及機械產(chǎn)品的設(shè)計質(zhì)量和效率，具有很強的實用性。同時也可以為專業(yè)和非專業(yè)用戶提供簡便的滑動軸承－轉(zhuǎn)子系統(tǒng)性能分析結(jié)果。

【1】張直明.滑動軸承的流體動力潤滑理論[M].北京:高等教育出版社，1986.

【2】楊沛然.流體潤滑數(shù)值分析[M].北京:國防工業(yè)出版社，1998.

【3】吳超，王文，陳曉陽，等.推力軸承支承方式及間隙影響研究[J].潤滑與密封，2006，31(5):130－131.

【4】RABINOWITZ M D，HAHN E J.Steady State Performance of Squeeze Film Damper Supported Flexible Rotor[J].Trans ASME，Journal of Engineering Power，1977，99(4):552－558.

【5】高小沖，李運華.基于電液可控擠壓油膜阻尼器的轉(zhuǎn)子系統(tǒng)振動的主動控制研究[J].機床與液壓，2002(5):11－13.

【6】SUN L，KRODKIEWSKI J M.Experimental Investigation of Dynamic Properties of an Active Journal Bearing[J].Journal of Sound and Vibration，2000，230(5):1103－1117.

【7】吳超，楊春燕，吳溢華，等.主動可控的徑向油膜軸承研究綜述[J].潤滑與密封，2009，34(8):111－113.

【8】TAYLOR R I.Inclusion of Lubricant Shear Thinning in the Short Bearing Approximation [J].Journal of Engineering Tribology，1999，213(1):35－46.

【9】HATTORI H.Dynamic Analysis of a Rotor-journal Bearing System with Large Dynamic Loads (Stiffness and Damping Coefficient Variations in Bearing Oil Film)[J].JSME International，1993，36(2):251－257.

【10】ADILETTA G，GUIDO A R，ROSSI C.Chaotic Motions of a Rigid Rotor in Short Journal Bearings[J].Nonlinear Dynamics，1996，10(3):251－269.

【11】王文，張直明.油葉型軸承非線性油膜力數(shù)據(jù)庫[J].上海工業(yè)大學(xué)學(xué)報，1993(4):299－305.

【12】焦映厚，李明章，陳照波.不同油膜模型下轉(zhuǎn)子－圓柱軸承系統(tǒng)的動力學(xué)分析[J].哈爾濱大學(xué)學(xué)報，2007，39(1):46－50.

【13】袁小陽，朱均.多圓盤轉(zhuǎn)子滑動軸承系統(tǒng)自激振動的計算與分析[J].西安交通大學(xué)學(xué)報，1997，31(8):80－86.

【14】SUNDARARAJAN P，NOAH S T.Dynamics of Forced Nonlinear Systems Using Shooting/Arc-length Continuation Method-application to Rotor Systems[J].ASME Journal of Vibration and Acoustics，1997，119:9－20.

【15】李立，許慶余，鄭鐵生.非線性轉(zhuǎn)子－軸承系統(tǒng)瞬態(tài)響應(yīng)求解的分塊直接積分法[J].振動工程學(xué)報，1996，9(l):64－68.

【16】李志剛，張直明.多跨轉(zhuǎn)子－滑支軸承系統(tǒng)非線性動力學(xué)仿真[J].自然雜志，l997，l9(A10):76－82.

【17】CHEN Z，JIAO Y，XIA S，et al.An Efficient Calculation Method of Nonlinear Fluid Film Forces in Journal Bearing[J].Tribology Transactions，2002，45(3):324－329.

【18】秦平，沈鉞，徐華，等.考慮進油壓力的滑動軸承非線性油膜力數(shù)據(jù)庫[J].摩擦學(xué)學(xué)報，2004，24(3):258－262.

【19】呂延軍，張永芳，于楊冰，等.固定瓦－可傾瓦組合徑向軸承－柔性轉(zhuǎn)子系統(tǒng)非線性運動分析[J].振動與沖擊，2011，30(5):257－262.

【20】崔曉康，吳超，吳溢華，等.基于VB 和Fortran 語言混合編程的滑動軸承性能分析的可視化研究[J].機床與液壓，2011，39(21):91－93.

【21】王建梅，趙春江，姚建斌，等.油膜軸承可視化性能計算系統(tǒng)的設(shè)計及實現(xiàn)[J].太原科技大學(xué)學(xué)報，2006，27(S1):17－19.

【22】陽旭，王培俊，楊利明，等.滑動軸承三維油膜壓力動態(tài)分布可視化研究[J].機械工程與自動化，2010(2):14－16.

【23】吳剛，李濤生，嚴(yán)新平.滑動徑向軸承軸心軌跡的可視化及其應(yīng)用[J].機械科學(xué)與技術(shù)，2004，23(1):117－119.

【24】范紅紅，張小棟，李單龍.基于OpenGL技術(shù)的滑動軸承油膜特性的可視化研究[J].潤滑與密封，2008，33(11):81－83.

【25】郅剛鎖，朱均.油膜流場可視化的實時性研究[J].計算機輔助設(shè)計與圖形學(xué)學(xué)報，2004，16(1):116－120.

【26】王丹丹，岑少起，郭紅.基于VC/MFC的滑動軸承設(shè)計軟件開發(fā)[J].煤礦機械，2010，31(4):217－219.