竇 唯，褚寶鑫

（北京航天動(dòng)力研究所，北京 100076）

0 引言

低溫液體火箭發(fā)動(dòng)機(jī)渦輪泵轉(zhuǎn)子通常在高轉(zhuǎn)速下工作，轉(zhuǎn)子的振動(dòng)特性直接影響火箭發(fā)動(dòng)機(jī)的性能。軸承為旋轉(zhuǎn)軸系提供剛度，其軸向位置對(duì)轉(zhuǎn)子系統(tǒng)動(dòng)態(tài)特性有著重要的影響，包括轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速、振型、不平衡響應(yīng)以及穩(wěn)定性特性等[1-15]。因此，開(kāi)展軸承軸向位置對(duì)渦輪泵轉(zhuǎn)子系統(tǒng)穩(wěn)定性的影響研究具有重大的理論和工程實(shí)際意義。

本文針對(duì)低溫液體火箭發(fā)動(dòng)機(jī)渦輪泵轉(zhuǎn)子非線性系統(tǒng)開(kāi)展了軸承軸向位置對(duì)動(dòng)力穩(wěn)定性的影響研究。建立了渦輪泵轉(zhuǎn)子系統(tǒng)的非線性動(dòng)力學(xué)模型，分別研究了理想安裝情況下和非理想安裝情況下泵端和渦輪端軸承軸向位置變化對(duì)轉(zhuǎn)子系統(tǒng)穩(wěn)定性的影響，給出了失穩(wěn)轉(zhuǎn)速隨軸向位置的變化規(guī)律，為液體火箭發(fā)動(dòng)機(jī)渦輪泵轉(zhuǎn)子系統(tǒng)結(jié)構(gòu)設(shè)計(jì)、故障診斷與安裝維護(hù)提供理論依據(jù)。

1 液體火箭發(fā)動(dòng)機(jī)渦輪泵轉(zhuǎn)子系統(tǒng)數(shù)學(xué)模型

渦輪泵轉(zhuǎn)子軸系主要包括轉(zhuǎn)子、輪盤、軸承和密封等，本文采用有限元方法形成渦輪泵轉(zhuǎn)子軸系的剛度、質(zhì)量和陀螺矩陣；進(jìn)而建立渦輪泵轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)模型。

1.1 非線性密封力模型

渦輪泵轉(zhuǎn)子系統(tǒng)在實(shí)際運(yùn)行時(shí)，其密封結(jié)構(gòu)不僅起到減弱流體逸出的作用，同時(shí)近似為滑動(dòng)軸承結(jié)構(gòu)，滑動(dòng)軸承的交叉剛度項(xiàng)可引起軸系的失穩(wěn)。本研究在計(jì)算時(shí)，將密封結(jié)構(gòu)做滑動(dòng)軸承處理，并采用Capone非線性密封力模型。該模型由Capone于1986年推導(dǎo)得出，是基于短圓柱瓦軸承假設(shè)得到的[5-6]，即忽略了密封壓力沿周向變化產(chǎn)生的影響。

徑向滑動(dòng)軸承的幾何關(guān)系如下圖1所示。

軸承長(zhǎng)度即垂直于紙面方向?yàn)閦軸；c為軸承的半徑間隙；圖中O為軸承的幾何中心，OJ為軸頸的幾何中心，R為軸頸半徑；W為軸頸所承受的載荷；Ω為轉(zhuǎn)子的轉(zhuǎn)速；e為偏心距；θ為轉(zhuǎn)角；φ為偏位角。

由于Capone非線性密封力模型是基于短軸承理論假設(shè)提出的，因此可以忽略Reynold方程中壓力梯度沿周向的變化，即可得到流體密封中的壓力的微分方程的無(wú)量綱形式，可以簡(jiǎn)寫為下式

可以獲得無(wú)量綱流體密封壓力表達(dá)式為

如圖1所示，當(dāng)轉(zhuǎn)角θ在區(qū)間[β，β+π]內(nèi)時(shí)，流體密封所形成的壓力區(qū)域?yàn)檎龎簠^(qū)，其中，角度β的值可由下式得到

由自由邊界條件可知，流體密封在空穴區(qū)和軸承兩端的壓力為零，故在軸承邊界上流體密封壓力為0。

對(duì)式 (2)沿著軸頸表面密封流體作用的弧段積分，就可得到非線性流體密封力的無(wú)量綱表達(dá)式如式 (4)所示

對(duì)上式積分，就可得到無(wú)量綱非線性密封力fx和fy的最后的表達(dá)式

將根據(jù)短圓柱瓦軸承理論得到的無(wú)量綱的Capone非線性密封力fx和fy轉(zhuǎn)換成有量綱的非線性密封力Fx和Fy，得到的表達(dá)式如式(9)所示

式中：σ為軸承特性參數(shù)；μ為軸承密封流體粘度；Ω為轉(zhuǎn)子轉(zhuǎn)速；R為軸頸半徑；L為軸承有效長(zhǎng)度；D為軸承直徑。

1.2 轉(zhuǎn)子系統(tǒng)數(shù)學(xué)模型

將渦輪泵轉(zhuǎn)子劃分為若干個(gè)軸段，每個(gè)軸段單元采用Timoshenko梁?jiǎn)卧Ｐ?，這樣就可建立渦輪泵轉(zhuǎn)子系統(tǒng)的無(wú)阻尼自由運(yùn)動(dòng)微分方程。

彈性軸的無(wú)阻尼自由運(yùn)動(dòng)微分方程如式 (11)所示

式中：q為彈性軸節(jié)點(diǎn)的位移向量；Ms為彈性軸的一致 (協(xié)調(diào))質(zhì)量矩陣；Ks為彈性軸的剛度矩陣；Gs為彈性轉(zhuǎn)軸的陀螺矩陣。

由動(dòng)能定理及根據(jù)Lagrange方程可以得到輪盤的運(yùn)動(dòng)方程如式 (12)所示

式中：Q(t)為作用在輪盤上的外力，N。

寫成矩陣形式為

式中：Md為輪盤的質(zhì)量矩陣；Gd為輪盤的陀螺矩陣。

由于制造加工及裝配等誤差的影響，實(shí)際運(yùn)行的轉(zhuǎn)子輪盤系統(tǒng)不可避免地會(huì)存在不平衡量，且輪盤部分的不平衡量遠(yuǎn)大于彈性軸段部分。本文在計(jì)算轉(zhuǎn)子軸系的穩(wěn)定性時(shí)只考慮輪盤結(jié)構(gòu)(即渦輪、離心輪、誘導(dǎo)輪)的影響。

不平衡激勵(lì)力的表達(dá)式為

式中：mi為輪盤的質(zhì)量；ei為該輪盤的偏心距；αi為該輪盤的初始偏位角。將不平衡力寫為向量的形式為

式中：Fu,i即為輪盤不平衡量在轉(zhuǎn)子節(jié)點(diǎn)i處所產(chǎn)生的不平衡激勵(lì)力，表達(dá)式中各行對(duì)應(yīng)轉(zhuǎn)子節(jié)點(diǎn)i 處的四個(gè)自由度

然后根據(jù)渦輪泵軸系的軸承所在節(jié)點(diǎn)及作用在軸承上的載荷，從而可計(jì)算四個(gè)軸承剛度系數(shù)和四個(gè)支承結(jié)構(gòu)的阻尼系數(shù)（不考慮游隙的影響）。

若軸頸中心的坐標(biāo)為xs和ys，則軸承作用于渦輪機(jī)節(jié)點(diǎn)處的軸承力和分別為

將線性軸承力移到系統(tǒng)的運(yùn)動(dòng)微分方程的左端，對(duì)軸承所在的節(jié)點(diǎn)而言，系統(tǒng)的運(yùn)動(dòng)微分方程中與節(jié)點(diǎn)的振動(dòng)有關(guān)的項(xiàng)可以寫成下式

因此，可以把每個(gè)軸承的剛度系數(shù)分別根據(jù)軸承所在的節(jié)點(diǎn)位置加到彈性軸的剛度矩陣上Ksi(i=1,2)，得到考慮了軸承剛度的轉(zhuǎn)子－軸承系統(tǒng)的剛度矩陣Ksbi(i=1,2)。

把輪盤的質(zhì)量矩陣Md按照輪盤所在的節(jié)點(diǎn)對(duì)應(yīng)的自由度加到彈性轉(zhuǎn)軸的一致質(zhì)量矩陣Ms對(duì)應(yīng)的位置上，形成轉(zhuǎn)子－輪盤系統(tǒng)的質(zhì)量矩陣。同理，把輪盤的陀螺矩陣Gd按照輪盤所在位置加到轉(zhuǎn)子的陀螺矩陣Gs上，形成轉(zhuǎn)子－輪盤系統(tǒng)的陀螺矩陣。

按照節(jié)點(diǎn)的排列順序?qū)⑥D(zhuǎn)子－輪盤系統(tǒng)的質(zhì)量、剛度和陀螺矩陣組合成系統(tǒng)的整體質(zhì)量矩陣M，整體剛度矩陣K和整體陀螺矩陣G。

根據(jù)系統(tǒng)的剛度矩陣、陀螺矩陣和質(zhì)量矩陣，可得到轉(zhuǎn)子－輪盤－軸承系統(tǒng)的無(wú)阻尼自由振動(dòng)微分方程為

考慮不平衡激勵(lì)力、密封非線性液膜力及轉(zhuǎn)子-軸承系統(tǒng)的阻尼作用時(shí)，可得到轉(zhuǎn)子－輪盤－軸承系統(tǒng)的有阻尼強(qiáng)迫振動(dòng)微分方程為

1.3 渦輪泵轉(zhuǎn)子模型簡(jiǎn)化

渦輪泵轉(zhuǎn)子系統(tǒng)主要由轉(zhuǎn)子、離心輪和誘導(dǎo)輪三部分組成。工作時(shí)，三部分結(jié)構(gòu)一起旋轉(zhuǎn)。建立有限元模型時(shí)，須對(duì)各部分連接區(qū)域進(jìn)行適當(dāng)簡(jiǎn)化。該轉(zhuǎn)子系統(tǒng)主要通過(guò)兩個(gè)滾動(dòng)軸承支承，其中密封結(jié)構(gòu)也能起到一定的支承剛度作用，在穩(wěn)定性分析中均將其作滑動(dòng)軸承近似處理。文中所提支承剛度均指軸承、鼠籠及阻尼器的總剛度。綜上，渦輪泵轉(zhuǎn)子系統(tǒng)的有限元模型如圖2所示，軸系總長(zhǎng)420 mm，共劃分為41個(gè)梁?jiǎn)卧?個(gè)輪盤單元，由兩個(gè)滾動(dòng)軸承及兩個(gè)迷宮密封結(jié)構(gòu)共同支承。

模型主要進(jìn)行以下簡(jiǎn)化：將渦輪轉(zhuǎn)子的軸套結(jié)構(gòu)等效為轉(zhuǎn)子的一部分；輪盤、誘導(dǎo)輪和離心輪簡(jiǎn)化為規(guī)則結(jié)構(gòu)，并附加部分集中質(zhì)量和轉(zhuǎn)動(dòng)慣量；動(dòng)環(huán)結(jié)構(gòu)以集中質(zhì)量、集中轉(zhuǎn)動(dòng)慣量近似處理。

2 液體火箭發(fā)動(dòng)機(jī)渦輪泵轉(zhuǎn)子系統(tǒng)非線性穩(wěn)定性分析

軸系穩(wěn)定性分析主要計(jì)算軸系由于低頻渦動(dòng)產(chǎn)生液膜振蕩的頻率大小及低頻振蕩出現(xiàn)時(shí)的轉(zhuǎn)速，即失穩(wěn)轉(zhuǎn)速。因?yàn)榱黧w密封力是非線性因素，并且系統(tǒng)的動(dòng)力學(xué)模型是高維方程，本文采用直接積分法來(lái)求解系統(tǒng)的動(dòng)態(tài)響應(yīng)。因?yàn)镹ewmark方法的2個(gè)參數(shù)滿足一定值時(shí)，該方法是無(wú)條件穩(wěn)定的，所以本文采用Newmark方法逐步積分，計(jì)算一系列不同轉(zhuǎn)速下轉(zhuǎn)子的動(dòng)力學(xué)響應(yīng)，確定轉(zhuǎn)子的低頻成分，其中轉(zhuǎn)速間隔取為200 r/min。數(shù)值仿真無(wú)量綱時(shí)間步長(zhǎng)取π/500，Newmark方法的收斂偏差標(biāo)準(zhǔn)取2×10-6。仿真500個(gè)周期，舍棄非穩(wěn)態(tài)數(shù)據(jù)，取后50個(gè)周期穩(wěn)態(tài)數(shù)據(jù)進(jìn)行分析。

軸系的穩(wěn)定性計(jì)算采用梁?jiǎn)卧Ｐ途幊虒?shí)現(xiàn)，計(jì)算時(shí)簡(jiǎn)化了轉(zhuǎn)子結(jié)構(gòu)，去掉輪盤結(jié)構(gòu)，將輪盤的質(zhì)量、轉(zhuǎn)動(dòng)慣量附加到相應(yīng)質(zhì)心處；另一方面，穩(wěn)定性計(jì)算時(shí)密封結(jié)構(gòu)采用上文的非線性流體密封力模型。

理想情況下，軸系在安裝時(shí)各軸段軸線及密封、軸承的中心軸線等完全重合。實(shí)際工作中，軸系安裝時(shí)總會(huì)存在偏心，此時(shí)密封結(jié)構(gòu)的液膜厚度會(huì)發(fā)生變化，并影響轉(zhuǎn)子的動(dòng)力特性。為此，分析中引入安裝偏心影響因素?？紤]到安裝偏心小于密封間隙，因此給定為10 μm。

2.1 無(wú)安裝偏心時(shí)泵端軸承位置變化時(shí)穩(wěn)定性分析

將軸承在模型原始位置處軸向調(diào)整±5 mm，計(jì)算軸承位置對(duì)軸系穩(wěn)定性的影響，由圖2可知，離心輪和渦輪分別支承軸承的外側(cè)，其中泵端軸承軸向位置為負(fù)值表示靠近離心輪，渦輪端軸承軸向位置為負(fù)值表示遠(yuǎn)離渦輪輪盤。分別取泵端軸承軸向位置為-5 mm，0 mm，5 mm三個(gè)值，計(jì)算不同泵端軸承軸向位置時(shí)失穩(wěn)轉(zhuǎn)速變化規(guī)律，不同泵端軸承軸向位置下失穩(wěn)轉(zhuǎn)速結(jié)果如表1所示，失穩(wěn)轉(zhuǎn)速隨泵端軸承軸向位置的變化曲線見(jiàn)圖3。

表1 泵端軸承位置對(duì)失穩(wěn)轉(zhuǎn)速影響Tab.1 Influence of position of bearing near pump on instability rotation speed

從表1中可以看出，當(dāng)泵端軸承軸向位置在-5～5 mm范圍內(nèi)調(diào)整時(shí)，失穩(wěn)轉(zhuǎn)速變化明顯，軸向位置為-5 mm時(shí)失穩(wěn)轉(zhuǎn)速值已經(jīng)超出了計(jì)算范圍，且隨著泵端軸承遠(yuǎn)離密封結(jié)構(gòu)軸系失穩(wěn)轉(zhuǎn)速減小，這是因?yàn)殡S著泵端軸承逐漸遠(yuǎn)離密封結(jié)構(gòu)，在轉(zhuǎn)子靜撓曲及泵端軸承處轉(zhuǎn)子靜變形相對(duì)不變的條件下，轉(zhuǎn)子在密封處?kù)o變形加大，密封激勵(lì)力增大，密封非線性激勵(lì)力越強(qiáng)，引發(fā)軸系不穩(wěn)定因素增多，從而導(dǎo)致軸系更容易失穩(wěn)。

為了分析低頻成分的變化規(guī)律，圖4中給出了無(wú)安裝偏心時(shí)不同泵端軸承軸向位置下離心輪后凸肩密封節(jié)點(diǎn)的三維譜圖。

當(dāng)安裝偏心為0時(shí)，不同泵端軸承軸向位置下軸系也只有一階渦動(dòng)，且一階渦動(dòng)頻率均在12 600 r/min附近。從三維譜圖的基頻曲線也可以看出，泵端軸承軸向位置在小幅度內(nèi)變化時(shí)，軸系臨界轉(zhuǎn)速的變化也很小，一階臨界轉(zhuǎn)速約為9 000 r/min，二階臨界轉(zhuǎn)速約為24 000 r/min。

綜上，泵端軸承軸向位置對(duì)軸系的失穩(wěn)轉(zhuǎn)速影響較大，但對(duì)一階渦動(dòng)的幅值及頻率影響較小。隨著泵端軸承逐漸遠(yuǎn)離密封結(jié)構(gòu)，失穩(wěn)轉(zhuǎn)速逐漸減??；一階渦動(dòng)的頻率約為12 600 r/min。

2.2 無(wú)安裝偏心時(shí)渦輪端軸承位置變化時(shí)穩(wěn)定性分析

分別取渦輪端軸承軸向位置為-5 mm，0 mm，5 mm三個(gè)值研究不同渦輪端軸承軸向位置時(shí)失穩(wěn)轉(zhuǎn)速變化規(guī)律，不同渦輪端軸承軸向位置下失穩(wěn)轉(zhuǎn)速結(jié)果見(jiàn)表2所示，失穩(wěn)轉(zhuǎn)速隨渦輪端軸承軸向位置的變化曲線見(jiàn)圖5。

表2 渦輪端軸承軸向位置對(duì)失穩(wěn)轉(zhuǎn)速影響Tab.2 Influence of radial position of bearing near turbine on instability rotation speed

從表2中可以看出，渦輪端軸承軸向位置變化時(shí)，失穩(wěn)轉(zhuǎn)速也變化顯著，但變化幅度小于泵端軸承軸向位置時(shí)的結(jié)果，且隨著渦輪端軸承靠近渦輪輪盤，軸系失穩(wěn)轉(zhuǎn)速增大，這是因?yàn)殡S著渦輪端軸承逐漸遠(yuǎn)離密封結(jié)構(gòu) (即靠近渦輪輪盤)，相當(dāng)于泵端軸承逐漸靠近密封，故穩(wěn)定性結(jié)果與改變泵端軸承軸向位置時(shí)的結(jié)果相反。但由于該作用為間接的，因而沒(méi)有直接改變泵端軸承軸向位置時(shí)效果顯著。

為分析低頻成分的變化規(guī)律，圖6中給出無(wú)安裝偏心時(shí)不同渦輪端軸承軸向位置下離心輪后凸肩密封節(jié)點(diǎn)的三維譜圖。其中，渦輪端軸承軸向位置為0 m時(shí)的結(jié)果在圖4（b）中已經(jīng)給出，此處予以省略。

同樣，當(dāng)安裝偏心為0時(shí)，不同渦輪端軸承軸向位置下軸系也只有一階渦動(dòng)，且一階渦動(dòng)頻率均在12 600 r/min附近；雖然失穩(wěn)轉(zhuǎn)速變化顯著，但一階渦動(dòng)的幅值幾乎不隨軸向位置的變化而變化。從三維譜圖的基頻曲線也可看出，渦輪端軸承軸向位置在小幅度內(nèi)變化時(shí)，軸系臨界轉(zhuǎn)速的變化也很小，一階臨界轉(zhuǎn)速約為9 000 r/min，二階臨界轉(zhuǎn)速約為24 000 r/min。

綜上，渦輪端軸承軸向位置對(duì)軸系的失穩(wěn)轉(zhuǎn)速影響仍較大，但小于泵端軸承軸向位置的影響效果，對(duì)一階渦動(dòng)的幅值及頻率影響很小。隨著渦輪端軸承逐漸靠近渦輪輪盤結(jié)構(gòu)，失穩(wěn)轉(zhuǎn)速逐漸增大；一階渦動(dòng)的頻率約為12 600 r/min。

2.3 有安裝偏心時(shí)泵端軸承軸向位置變化時(shí)穩(wěn)定性分析

與安裝偏心等于0時(shí)的取值相同，計(jì)算安裝偏心為10 μm時(shí)不同泵端軸承軸向位置下的失穩(wěn)轉(zhuǎn)速，其結(jié)果見(jiàn)表3。此時(shí)失穩(wěn)轉(zhuǎn)速隨軸向位置的變化曲線見(jiàn)圖7（與安裝偏心等于0時(shí)的結(jié)果進(jìn)行對(duì)比）。

表3 泵端軸承軸向位置對(duì)失穩(wěn)轉(zhuǎn)速影響Tab.3 Influence of radial position of bearing near pump on instability rotation speed

從圖7中可以看出，增大安裝偏心后，泵端軸承軸向位置對(duì)失穩(wěn)轉(zhuǎn)速的影響仍然較大，此時(shí)失穩(wěn)轉(zhuǎn)速變化規(guī)律與安裝偏心為0時(shí)的結(jié)果基本相同，變化幅度略有增大，但數(shù)值較安裝偏心為0時(shí)增大較大，軸向位置為-5 mm時(shí)失穩(wěn)轉(zhuǎn)速值也超出了計(jì)算范圍。

為分析低頻成分的變化規(guī)律圖8中給出有安裝偏心時(shí)不同泵端軸承軸向位置下離心輪后凸肩密封節(jié)點(diǎn)的三維譜圖。

安裝偏心增大后，軸系出現(xiàn)二階、三階渦動(dòng)；隨階數(shù)增加，渦動(dòng)幅值迅速降低，但各階渦動(dòng)幅值基本不隨泵端軸承軸向位置的改變而改變；一階、二階渦動(dòng)頻率近似為12 600 r/min，24 000 r/min。

綜上所述，增大安裝偏心后，泵端軸承軸向位置對(duì)軸系的失穩(wěn)轉(zhuǎn)速影響仍然較大，但對(duì)各階渦動(dòng)的幅值及頻率大小影響很??；同時(shí)，安裝偏心增大會(huì)導(dǎo)致二階、三階渦動(dòng)的出現(xiàn)。

2.4 有安裝偏心時(shí)渦輪端軸承軸向位置變化時(shí)穩(wěn)定性分析

與安裝偏心等于0時(shí)的取值相同，計(jì)算安裝偏心為10 μm時(shí)不同渦輪端軸承軸向位置下的失穩(wěn)轉(zhuǎn)速，其結(jié)果見(jiàn)表4。此時(shí)失穩(wěn)轉(zhuǎn)速隨軸向位置的變化曲線見(jiàn)圖9（與安裝偏心等于0時(shí)的結(jié)果進(jìn)行對(duì)比）。

表4 有安裝偏心時(shí)渦輪端軸承位置對(duì)失穩(wěn)轉(zhuǎn)速影響Tab.4 Influence of position of bearing near turbine on instability rotation speed when installation eccentricity exists

從圖9中可以看出，增大安裝偏心后，渦輪端軸承軸向位置對(duì)失穩(wěn)轉(zhuǎn)速的影響仍然較大，此時(shí)失穩(wěn)轉(zhuǎn)速變化規(guī)律與安裝偏心為0時(shí)的結(jié)果基本相同，變化幅度略有增大，但數(shù)值較安裝偏心為0時(shí)增大較大，同時(shí)軸向位置為5 mm時(shí)失穩(wěn)轉(zhuǎn)速值超出了計(jì)算范圍。

為分析低頻成分的變化規(guī)律，圖10中給出了有安裝偏心時(shí)不同渦輪端軸承軸向位置下離心輪后凸肩密封節(jié)點(diǎn)的三維譜圖。其中，安裝偏心10 μm，渦輪端軸承軸向位置為0時(shí)的結(jié)果在圖8中已經(jīng)給出，此處予以省略。

與泵端軸承軸向位置變化時(shí)的結(jié)果類似，安裝偏心增大后，軸系開(kāi)始出現(xiàn)二階和三階渦動(dòng)；隨著階數(shù)增加，渦動(dòng)幅值迅速降低，三階渦動(dòng)的幅值幾乎為零，但各階渦動(dòng)幅值基本不隨渦輪端軸承軸向位置的改變而改變；一階、二階渦動(dòng)的頻率近似為12 600 r/min和24 000 r/min。

綜上所述，增大安裝偏心后，渦輪端軸承軸向位置對(duì)軸系的失穩(wěn)轉(zhuǎn)速影響仍然較大，但對(duì)各階渦動(dòng)的幅值及頻率大小影響很小。

3 結(jié)論

通過(guò)研究軸承軸向位置對(duì)液體火箭發(fā)動(dòng)機(jī)轉(zhuǎn)子系統(tǒng)穩(wěn)定性的影響研究，得出以下結(jié)論：

1)安裝偏心為0時(shí)，不論軸承軸向位置如何變化，軸系均只存在一階渦動(dòng)；增大安裝偏心后，軸系開(kāi)始出現(xiàn)二階、三階渦動(dòng)。

2)各階渦動(dòng)幅值和頻率基本不受軸承軸向位置參數(shù)變化影響，隨階數(shù)增加，各階渦動(dòng)幅值迅速降低，一階、二階渦動(dòng)頻率近似為1 260 r/min，24 000 r/min；隨著安裝偏心增大，低頻渦動(dòng)幅值基本呈現(xiàn)先增大后減小的趨勢(shì)，但一階渦動(dòng)幅值隨安裝偏心增大逐漸減小。

3)渦輪端軸承軸向位置對(duì)軸系的失穩(wěn)轉(zhuǎn)速影響較大；泵端、渦輪端軸承軸向位置變化時(shí)失穩(wěn)轉(zhuǎn)速變化規(guī)律相反，泵端軸承遠(yuǎn)離密封結(jié)構(gòu)時(shí)失穩(wěn)轉(zhuǎn)速減小，而渦輪端軸承遠(yuǎn)離密封結(jié)構(gòu)時(shí)失穩(wěn)轉(zhuǎn)速增大。

[1]WYSSMANN H,JENNY R,PHAM T.Prediction of stiffness and damping coefficients for centrifugal compressor labyrinth seals[J].ASME Journal of Engineering for Gas Turbines and Power,1984,106(4):920-926.

[2]NOSAKA M.Self-lubricating performance and durability of ball bearings for the LE-7 liquid oxygen rocketturbopump[J].Lubrication Engineering,1993,49(9):677-688.

[3]TIWARI M,GUPTA K,PRAKASH O.Effect of radial internal clearance of a ball bearing on the dynamics of a balanced horizontal rotor[J].Journal of Sound and Vibration,2000,238:723-756.

[4]CHIRLDS D W.The space shuttle main engine high pressure fuel turbopump rotordynamic instability problem[J].Journal of Eng.Gas Turbines Power,1978,100(1):48-57

[5]ARGHIR M,FRENE J.Rotordynamic coefficients of circumferentially grooved liquid seals using the averaged Navier-Stokes equations[J].ASME Journal of Tribology,1997,119:556-567.

[6]MUSZYNSKA A,BENTLY D E.Frequency-swept rotating input perturbation techniques and identification of the fluid force models in rotor/bearing/seal systems and fluid handling machines[J].Journal of Sound and Vibration,1990,143(1):103-124.

[7]BROMMUNDT E,OSTERMEYER G P.Stability of a rotor partially filled with liquid and attached to an anisotropically mounted shaft[J].Ingenieur-Archiv,1986,56(3):379-388.

[8]HOLM-CHRISTENSEN O,TRAGER K.A note of instability caused by liquid motions[J].Journal of Applied Mech.,1991,58(4):801-811.

[9]ZHANG W.Dynamic stability of a rotor filled or partially filled with liquid[J].Journal of Applied Mech.,1996,64(1):101-105.

[10]沈紫樂(lè),李凌均,雷新望.轉(zhuǎn)子系統(tǒng)中軸承位置對(duì)聯(lián)軸器不對(duì)中量影響的研究 [J].礦山機(jī)械，2011,39(9):113-116.

[11]牛忠俊,花春.軸系橫向振動(dòng)邊界參數(shù)的確定及軸承位置優(yōu)化 [J].2011泰州職業(yè)技術(shù)學(xué)院學(xué)報(bào),2011,11(3):110-114.

[12]白長(zhǎng)青,許慶余,張小龍.滾動(dòng)軸承-火箭發(fā)動(dòng)機(jī)液氫渦輪泵轉(zhuǎn)子系統(tǒng)的動(dòng)力特性分析[J].航空學(xué)報(bào),2006,27(2):258-261.

[13]應(yīng)桂爐.氫渦輪泵試驗(yàn)振動(dòng)分析 [J].導(dǎo)彈與航天運(yùn)載技術(shù),1995(5):40-46.

[14]祝長(zhǎng)生.支承剛度各向異性部分充液轉(zhuǎn)子系統(tǒng)的穩(wěn)定性[J].航空動(dòng)力學(xué)報(bào),2000,15(4):431-434.

[15]張大義,母國(guó)新,洪杰.航空發(fā)動(dòng)機(jī)轉(zhuǎn)子支承系統(tǒng)剛度計(jì)算中的幾個(gè)問(wèn)題[J].戰(zhàn)術(shù)導(dǎo)彈技術(shù),2005(2):20-23.

[16]竇唯,劉占生.液體火箭發(fā)動(dòng)機(jī)渦輪泵轉(zhuǎn)子彎扭耦合振動(dòng)研究[J].火箭推進(jìn),2012,38(4):17-25.

[17]黃智勇,胡鐘兵,李惠敏.大功率、高轉(zhuǎn)速、高揚(yáng)程渦輪泵振動(dòng)分析與減振研究[J].火箭推進(jìn),2005,31(6):1-6.