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

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

陳以彪，羅貴火

(南京航空航天大學(xué)能源與動(dòng)力學(xué)院，江蘇南京210016)

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

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

中圖分類號(hào)：V231.92文獻(xiàn)標(biāo)識(shí)碼：A

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

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

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

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

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

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

(1)

(2)

邊界條件：

當(dāng)y=0時(shí)：u=U1，v=V1，w=0

當(dāng)y=h時(shí)：u=U2，v=V2，w=0

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

(3)

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

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

軸頸M點(diǎn)速度分量如下：

(4)

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

(5)

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

(6)

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

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

參照?qǐng)D3，可近似得到外層油膜的油膜間隙為：

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

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

U4=V4=0

(7)

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

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

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

式中的un、uw分別表示浮環(huán)內(nèi)、外表面相對(duì)于內(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í)驗(yàn)數(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)下，浮動(dòng)環(huán)轉(zhuǎn)速滿足力矩平衡條件

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

(8)

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

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

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

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

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

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

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

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

3雙向激勵(lì)試驗(yàn)

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

圖8　雙向激勵(lì)試驗(yàn)裝置簡圖

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

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

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

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

為比較傳統(tǒng)SFD與FSFD抑制轉(zhuǎn)子系統(tǒng)突加不平衡能力，本文選取一個(gè)對(duì)稱的Jeffcott轉(zhuǎn)子模型，仿真計(jì)算了兩種擠壓油膜阻尼器轉(zhuǎn)子系統(tǒng)在突加不平衡時(shí)的響應(yīng)，并分析了轉(zhuǎn)速對(duì)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具有更強(qiáng)的抑制突加不平衡的能力。由圖13知，隨著轉(zhuǎn)速的增加，F(xiàn)SFD-轉(zhuǎn)子系統(tǒng)的突加不平衡瞬態(tài)幅值增大且瞬態(tài)過程變長。ω=300 rad/s時(shí)，接近轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速，故瞬態(tài)響應(yīng)幅值較大。

5結(jié)束語

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

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

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

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

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

參考文獻(xiàn)

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

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

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

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

[5]晏礪堂, 張世平, 李其漢. 高效多孔環(huán)擠壓油膜阻尼器的減振特性研究[J]. 航空動(dòng)力學(xué)報(bào), 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]. 摩擦學(xué)學(xué)報(bào), 1999,19(3):261-265

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