浮環(huán)式擠壓油膜阻尼器減振機理研究

陳以彪，羅貴火

(南京航空航天大學能源與動力學院，江蘇南京210016)

摘要：從廣義雷諾方程出發(fā)，推導(dǎo)了浮環(huán)式擠壓油膜阻尼器(FSFD)內(nèi)、外層油膜的穩(wěn)態(tài)雷諾方程。根據(jù)文中的理論模型，基于有限差分法編寫了計算程序，研究分析了結(jié)構(gòu)參數(shù)對FSFD動力特性的影響。研究表明：與傳統(tǒng)SFD比較，F(xiàn)SFD改善了油膜力的非線性；相同條件下，F(xiàn)SFD內(nèi)層油膜力要大于外層油膜力，內(nèi)、外層油膜力都隨油膜寬度的增大而增大，隨油膜間隙的增大而減??；在結(jié)構(gòu)參數(shù)一定的情況下，F(xiàn)SFD抑制突加不平衡的能力要強于傳統(tǒng)SFD。

關(guān)鍵詞:廣義雷諾方程有限差分法浮環(huán)式擠壓油膜阻尼器結(jié)構(gòu)參數(shù)

中圖分類號：V231.92文獻標識碼：A

作者簡介：陳以彪，男，南京航空航天大學碩士研究生，研究方向：結(jié)構(gòu)強度與完整性。

收稿日期：2014-12-20

Mechanism study of the floating-ring squeeze film dampers

CHEN Yibiao, LUO Guihuo

Abstract:In this paper, Reynolds equations of the inner and outer film of floating-ring squeeze film dampers(FSFD)are derived from the general Reynolds equation. A calculation program is written based on finite difference method and according to the theoretical model in the article. The influence of the dynamic characteristics is studied for FSFD about the structural parameters. The results show that Nonlinear is obviously improved by the FSFD compared with traditional SFD, and the inner film force is greater than the outer film force under the same conditions. Furthermore, both inner and outer film forces increases with the increase of oil film width, but decreases with the increase of oil film clearance. Under the condition of fixed structure parameters, the performance of FSFD is better than traditional SFD in preventing sudden unbalance response.

Keywords:general Reynolds equation; finite difference method; floating-ring squeeze film dampers; structure parameters

0引言

擠壓油膜阻尼器(SFD)因結(jié)構(gòu)簡單、占用空間小、成本低、減振效果好等優(yōu)點，廣泛應(yīng)用于高速旋轉(zhuǎn)機械中。但由于SFD油膜力的高度非線性[1]，SFD-轉(zhuǎn)子系統(tǒng)容易出現(xiàn)雙穩(wěn)態(tài)響應(yīng)或非協(xié)調(diào)進動響應(yīng)等現(xiàn)象[2-3]，進而引發(fā)轉(zhuǎn)子系統(tǒng)產(chǎn)生較大振動甚至出現(xiàn)碰磨、疲勞等故障，限制了SFD的使用范圍。為此，國內(nèi)外學者提出了多種新型SFD，試圖抑制或改善SFD由于油膜力高度非線性帶來的危害。張世平、晏礪堂[4-5]研制了多孔擠壓油膜阻尼器，能提高轉(zhuǎn)子系統(tǒng)的不平衡承載能力以及工作的穩(wěn)定性和可靠性，但由于小孔容易被堵塞，故未能在實際中得到應(yīng)用。Hooshang Heshmat 和J.F. walton[6-7]等人研究設(shè)計了螺旋箔片多層擠壓油膜阻尼器(SFMSFD)，能有效改善油膜剛度的高度非線性，但對制造工藝水平要求較高。俄羅斯學者設(shè)計了彈性環(huán)式擠壓油膜阻尼器(ERSFD)，并將它成功地應(yīng)用于多種現(xiàn)役的航空發(fā)動機上[8-9]。

圖1　浮環(huán)式擠壓油膜 阻尼器結(jié)構(gòu)簡圖

此外，還有分段式擠壓油膜阻尼器(Segmented Dampers)、整體擠壓油膜阻尼器(ISFD)、結(jié)構(gòu)參數(shù)可控擠壓油膜阻尼器(VPSFD)、金屬橡膠阻尼器等[10-11]。這些新型SFD在某些方面改善了傳統(tǒng)SFD的不足，但由于結(jié)構(gòu)復(fù)雜或自身存在缺陷等方面的原因，未能在實際中廣泛應(yīng)用。本文對浮環(huán)式擠壓油膜阻尼器(FSFD)進行了研究，F(xiàn)SFD結(jié)構(gòu)簡圖見圖1。

1內(nèi)外層Reynolds方程的推導(dǎo)

取X軸沿軸切線方向，Y軸沿直徑方向，Z軸為軸向，u、v、w，分別表示沿X、Y、Z方向的油膜流動速度，h為油膜厚度。假設(shè)油膜壓力P沿Y方向不發(fā)生變化，油膜流動為層流，符合牛頓粘性定律，參考文獻[1]，則三維可壓縮流連續(xù)方程和運動方程可表示為：

(1)

(2)

邊界條件：

當y=0時：u=U1，v=V1，w=0

當y=h時：u=U2，v=V2，w=0

對方程(2)沿0到h連續(xù)兩次積分，代入方程(1)，結(jié)合邊界條件，整理可得不可壓縮流廣義雷諾方程的極坐標表示：

(3)

1.1內(nèi)層油膜Reynolds方程的推導(dǎo)

圖3　FSFD外層油膜間隙示意圖

軸頸M點速度分量如下：

(4)

浮動環(huán)內(nèi)表面M′點速度分量如下：

(5)

將式(4)和(5)代入式(3)可得內(nèi)層油膜壓力的Reynolds方程

(6)

式中：vd表示內(nèi)外層油膜上下兩表面向外流動的流體凈流速，θ、Z分別表示周向、周向坐標。

1.2外層油膜Reynolds方程的推導(dǎo)

參照圖3，可近似得到外層油膜的油膜間隙為：

h2=Rb-Rfo+e2cosθ=C2+e2cosθ

N點和N′點沿浮動環(huán)外表面的切向速度及徑向速度分量分別為U3、U4和V3、V4，油膜環(huán)靜止，故：

U4=V4=0

(7)

1.3浮動環(huán)自轉(zhuǎn)轉(zhuǎn)速的確定

參考壁面摩擦方程和Hirs整體流動理論[12-13]，則壁面剪應(yīng)力為

τn，w=0.5ρun,wfn,w

式中的un、uw分別表示浮環(huán)內(nèi)、外表面相對于內(nèi)、外層油膜的流速，滿足式：

un=[w2+(Rfiω1)2]1/2

uw=[w2+((1-ξ)Rfoω1)2]1/2

w=e1Ω1

式中，ξ為流速比，取為0.2；fn，w為摩擦參數(shù)，c1、c2、c3為Moody實驗數(shù)值(c1=1.38×10-3，c2=2×104，c3=5×105)；en、ew分別為內(nèi)、外表面粗糙度，en/h1=ew/h2=2×10-3。

穩(wěn)定狀態(tài)下，浮動環(huán)轉(zhuǎn)速滿足力矩平衡條件

∫02πτnRfidθ=∫02πτwRfodθ

(8)

由方程(8)可得出浮動環(huán)的自轉(zhuǎn)轉(zhuǎn)速。

2構(gòu)造參數(shù)對FSFD油膜力特性的影響

根據(jù)前面推導(dǎo)的FSFD雷諾方程，基于有限差分法，構(gòu)造求解油膜壓力場的非線性Newton松弛迭代格式，結(jié)合短軸承半油膜邊界條件，用MATLAB編寫程序，從而求得內(nèi)外層油膜壓力分布。轉(zhuǎn)子穩(wěn)定運轉(zhuǎn)時，軸頸進動速度與浮動環(huán)進動速度相等，為3 000rad/s，軸頸半徑R=14mm，內(nèi)、外層油膜間隙相等。初始條件下，間隙C=0.2mm，油膜寬度L=13mm，偏心率ε=0.1，滑油粘度μ=0.011 67Pa·s，進出口壓力均為一個大氣壓。

圖4為初始參數(shù)下FSFD內(nèi)、外層油膜力分布圖，內(nèi)層油膜力幅值要比外層油膜壓力幅值大。圖5為不同偏心率下，F(xiàn)SFD內(nèi)、外層最大油膜壓力沿周向坐標的變化曲線。圖6和圖7分別為內(nèi)、外層油膜力隨油膜寬度、油膜間隙的變化曲線。由圖可知：

圖4　偏心率ε=0.1時內(nèi)外層油膜壓力分布圖

圖5　不同偏心率下內(nèi)、外層最大油膜壓力沿周向變化曲線

1)隨著偏心率ε的增加，F(xiàn)SFD內(nèi)、外層最大油膜壓力逐漸增大，即使在0.4≤ε≤0.6范圍內(nèi)，油膜壓力增加的線性度仍然較好。而傳統(tǒng)SFD在偏心率ε>0.4時，油膜力增加的非線性會很明顯，故FSFD改善了油膜力高度的非線性，擴大工程適用范圍，具有較好的應(yīng)用前景。

2)在相同條件下，F(xiàn)SFD內(nèi)層油膜徑向力要大于外層油膜徑向力，內(nèi)層油膜切向力要大于外層油膜切向力；內(nèi)、外層油膜的徑向力、切向力都隨油膜寬度的增大而增大，隨油膜間隙的增大而減小。

3雙向激勵試驗

FSFD內(nèi)外層雷諾方程推導(dǎo)過程中作了一系列的假設(shè)，導(dǎo)致與實際工況有所不同，本文設(shè)計了雙向激勵試驗器來模擬FSFD實際工況，采用進口的壓力傳感器測試了油膜壓力，并將試驗結(jié)果與計算結(jié)果進行對比分析。圖8為雙向激勵試驗簡圖，通過激振器對芯棒產(chǎn)生兩個激勵力，激勵力幅值和頻率相同，相角相差90°，模擬轉(zhuǎn)子進動對油膜的擠壓作用。

圖8　雙向激勵試驗裝置簡圖

圖9為油膜力隨油膜間隙的變化曲線，試驗時內(nèi)外油膜間隙相同，從0.12 mm增加到0.22 mm。圖10為油膜力隨油膜寬度的變化曲線，油膜寬度從13 mm增加動33 mm。油膜壓力試驗測試結(jié)果與理論結(jié)果變化趨勢大體一致，即油膜壓力隨油膜間隙的增加而減小，隨油膜寬度的增加而增大。同時試驗結(jié)果與計算結(jié)果的誤差在10%以內(nèi)，表明本文所建計算模型的正確性。

圖9　油膜力與油膜間隙關(guān)系

圖10　油膜力與油膜寬度關(guān)系

4突加不平衡響應(yīng)分析

為比較傳統(tǒng)SFD與FSFD抑制轉(zhuǎn)子系統(tǒng)突加不平衡能力，本文選取一個對稱的Jeffcott轉(zhuǎn)子模型，仿真計算了兩種擠壓油膜阻尼器轉(zhuǎn)子系統(tǒng)在突加不平衡時的響應(yīng)，并分析了轉(zhuǎn)速對FSFD-轉(zhuǎn)子系統(tǒng)突加不平衡響應(yīng)的影響，結(jié)果見圖11-13所示。

圖11　SFD-轉(zhuǎn)子系統(tǒng)突　　圖12　FSFD-轉(zhuǎn)子系統(tǒng) 　　　加不平衡響應(yīng)　　　　　　突加不平衡響應(yīng)

圖13　不同轉(zhuǎn)速下FSFD-轉(zhuǎn)子系統(tǒng)突加不平衡響應(yīng)

由圖11和圖12知，F(xiàn)SFD-轉(zhuǎn)子系統(tǒng)的突加不平衡過程明顯更短，瞬態(tài)幅值明顯小。因此，與傳統(tǒng)SFD相比，F(xiàn)SFD具有更強的抑制突加不平衡的能力。由圖13知，隨著轉(zhuǎn)速的增加，F(xiàn)SFD-轉(zhuǎn)子系統(tǒng)的突加不平衡瞬態(tài)幅值增大且瞬態(tài)過程變長。ω=300 rad/s時，接近轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速，故瞬態(tài)響應(yīng)幅值較大。

5結(jié)束語

在考慮浮動環(huán)自轉(zhuǎn)的情況下，基于廣義雷諾方程，推導(dǎo)了FSFD內(nèi)外層油膜雷諾方程。根據(jù)本文推導(dǎo)的理論模型，基于有限差分法編制了FSFD 油膜特性分析程序，研究分析了結(jié)構(gòu)參數(shù)對FSFD動力特性的影響：

1)隨著偏心率ε的增加，F(xiàn)SFD內(nèi)、外層最大油膜壓力逐漸增大，即使在0.4≤ε≤0.6范圍內(nèi)，油膜壓力增加的線性度仍然較好。而傳統(tǒng)SFD在偏心率ε>0.4時，油膜力增加的非線性會很明顯，故FSFD改善了油膜力高度的非線性，擴大工程適用范圍，具有較好的應(yīng)用前景。

2)相同條件下，F(xiàn)SFD內(nèi)層油膜徑向力要大于外層油膜徑向力，內(nèi)層油膜切向力要大于外層油膜切向力；內(nèi)、外層油膜的徑向力、切向力都隨油膜寬度的增大而增大，隨油膜間隙的增大而減小。

3)設(shè)計了雙向激勵試驗，測量FSFD油膜壓力，試驗結(jié)果與理論結(jié)果變化趨勢是一致的，表明本文推導(dǎo)計算模型的正確性。試驗結(jié)果均小于理論結(jié)果，誤差在10%以內(nèi)，可能是由于滑油在油管中流動損耗掉部分油壓以及試驗裝置未完全密封引起的。

4)在結(jié)構(gòu)參數(shù)一定的情況下，F(xiàn)SFD抑制突加不平衡的能力要強于傳統(tǒng)SFD。隨著轉(zhuǎn)速的增加，F(xiàn)SFD-轉(zhuǎn)子系統(tǒng)的突加不平衡瞬態(tài)幅值增大且瞬態(tài)過程變長。

參考文獻

[1]周明, 李其漢, 晏礪堂. 彈性環(huán)式擠壓油膜阻尼器油膜力特性研究[J]. 燃氣渦輪試驗與研究, 1998, (3): 17-23

[2]馬艷紅, 洪杰, 趙福安. 自適應(yīng)擠壓油膜阻尼器減振機理理論研究[J]. 北京航空航天大學學報, 2004, (1): 5-8

[3]曹磊, 高德平, 江和甫. 彈性環(huán)式擠壓油膜阻尼器設(shè)計因素研究[J]. 燃氣渦輪試驗與研究, 2006, (4): 30-34

[4]張世平, 晏礪堂. 高效油膜阻尼器的研究與開發(fā)[J]. 航空動力學報, 1991,(2): 180-184

[5]晏礪堂, 張世平, 李其漢. 高效多孔環(huán)擠壓油膜阻尼器的減振特性研究[J]. 航空動力學報, 1993, (3): 225-233

[6]H. Heshmat, J. F. Walton. Advanced multi-squeeze film dampers for rotor vibration control [J]. Tribology transactions, 1991, 34(4): 489-496

[7]J. F.Walton, H. Heshmat. Rotor dynamic evaluation of an advanced multi squeeze film damper-imbalance response and blade-loss simulation [J].Journal of Engineering for Gas Turbines and Power, 1993, 115(2): 347-352

[8]M. Leontiev.Damper supports[R].Samsung aerospace. Engine Research & Development Center, 1996

[9]Andres L S, Kim T H. Analysis of gas foil bearings integrating FE top foil models [J]. Tribology International,2009,42:111-120

[10]D.L.Taylor, & V.S.Fehr. Analysis and design of segmented dampers for rotor dynamic control[J].Journal of Lubracation Technology, 1982：104

[11]Zeidan, F.L. San Andrés, and J. Vance. Design and application of squeeze film dampers in rotating machinery [C]. Proceedings of the 25th Turbomachinery Symposium, Turbomachinery Laboratory, Texas A&M University, September, 1996:169-188

[12]孫啟國, 姜培琳, 虞烈. 大間隙環(huán)流壁面摩擦及偏心轉(zhuǎn)子靜特性研究[J]. 摩擦學學報, 1999,19(3):261-265

[13]孫啟國, 虞烈. 有限長大間隙環(huán)流中同心轉(zhuǎn)子動特性系數(shù)研究[J]. 摩擦學學報,2001,21(6): 473-477