劉卓乾, 曹樹謙, 郭虎倫， 李利青

(1.天津大學 力學系,天津 300072; 2.天津市非線性動力學與混沌控制重點實驗室,天津 300072;3. 天津大學 內(nèi)燃機燃燒學國家重點實驗室，天津 300072)

?

含螺栓連接轉(zhuǎn)子系統(tǒng)非線性振動特性研究

劉卓乾1,2,3, 曹樹謙1,2,3, 郭虎倫1,2,3， 李利青1,2,3

(1.天津大學 力學系,天津 300072; 2.天津市非線性動力學與混沌控制重點實驗室,天津 300072;3. 天津大學 內(nèi)燃機燃燒學國家重點實驗室，天津 300072)

針對航空發(fā)動機轉(zhuǎn)子大量采用的螺栓法蘭連接，建立了連接結(jié)構的力學模型，定性給出了兩段轉(zhuǎn)子在相對彎曲時連接部位彎矩的分段線性參數(shù)化模型。應用狀態(tài)空間理論和數(shù)值計算方法對一個含螺栓連接的轉(zhuǎn)子系統(tǒng)進行了數(shù)值模擬，研究其模態(tài)特性和穩(wěn)態(tài)動力學響應。研究表明，分段線性變化連接剛度對中轉(zhuǎn)速段尤其是二階臨界轉(zhuǎn)速附近的模態(tài)和穩(wěn)態(tài)動力學響應產(chǎn)生較大影響，而轉(zhuǎn)子系統(tǒng)在低速段和高速段保持相對穩(wěn)定的、可線性化的運動。分析了預緊力對轉(zhuǎn)子運動特性的影響，發(fā)現(xiàn)預緊力對轉(zhuǎn)子系統(tǒng)運動穩(wěn)定性影響顯著，在一定范圍內(nèi)增大預緊力會使系統(tǒng)不穩(wěn)定加劇，但會縮短不穩(wěn)定運動的區(qū)間范圍。

轉(zhuǎn)子系統(tǒng)；螺栓連接；穩(wěn)定性；非線性動力學

航空發(fā)動機轉(zhuǎn)子的連接結(jié)構可分為拉桿連接、螺栓連接和無螺栓連接三種形式?，F(xiàn)代航空發(fā)動機連接結(jié)構的選用考證多種因素，除了保證轉(zhuǎn)子段間較好的一致性和振動穩(wěn)定性，為了滿足推重比要求，對于質(zhì)量的要求逐漸苛刻。螺栓連接結(jié)構簡單，連接穩(wěn)定性較好，較為輕便，因此廣泛應用于航空發(fā)動機中[1-2]。

由于螺栓法蘭連接本身的結(jié)構特性，航空發(fā)動機轉(zhuǎn)子系統(tǒng)存在“準剛性”特性。即除去剛性轉(zhuǎn)子應有的俯仰、平移振動外，實驗中在某些特殊的轉(zhuǎn)速下會發(fā)生錯動和彎曲，甚至兩段轉(zhuǎn)子相位上的不一致。而螺栓連接的剛度小于轉(zhuǎn)子的剛度，這使得螺栓連接成為故障發(fā)生的源頭之一。

壓氣機轉(zhuǎn)子典型的螺栓連接結(jié)構如圖1所示，其各方向力學特性均可用分段線性剛度描述。B?SWALD等[3-4]采用實驗與數(shù)值方法研究了螺栓連接軸向和橫向分段線性剛度的定性規(guī)律，提出了采用非線性傳遞函數(shù)方法識別橫向剛度的途徑。GAO等[5-6]用有限元法給出了拉桿轉(zhuǎn)子連接的數(shù)學模型，并將其數(shù)學模型應用于Jeffcott轉(zhuǎn)子，研究其動力學響應特性和模態(tài)特性。QIN等[7]用有限元法計算了螺栓轂筒連接的彎矩模型，給出了含彎矩模型的動力學響應特性。LIU等[8]結(jié)合有限元仿真和實驗，對螺栓法蘭連接的橫向剛度模型進行了對比研究，發(fā)現(xiàn)了剛度模型隨連接結(jié)構特性變化的一系列規(guī)律。

圖1 壓氣機轉(zhuǎn)子螺栓連接結(jié)構Fig.1 Bolt joint structure of aero-engines

國內(nèi)外學者對于航空發(fā)動機螺栓連接結(jié)構的研究，著重于靜態(tài)實驗研究或者有限元仿真計算，動力學數(shù)值仿真分析的研究相對較少。對于螺栓預緊力的研究給出了一些規(guī)律，但對動力學響應的影響研究較少。本文通過建立一個簡單的轉(zhuǎn)子連接模型，突出螺栓法蘭連接對于轉(zhuǎn)子系統(tǒng)的動力學特性影響，研究螺栓連接對轉(zhuǎn)子系統(tǒng)模態(tài)和非線性響應特性的影響。

1 轉(zhuǎn)子螺栓法蘭連接等效力學模型

為使相連的兩段轉(zhuǎn)子具有較好的轉(zhuǎn)動同步性以及保持足夠大的徑向以及軸向剛度，連接結(jié)構在三維平動方向上對兩段轉(zhuǎn)子提供相對約束，還具有抵抗連接的兩段轉(zhuǎn)子彎曲的作用。

1.1 螺栓法蘭連接的軸向力學特性

考慮一個螺栓法蘭連接的簡化模型，其剖面圖如圖2(a)所示。連接結(jié)構軸向力簡化為兩類，即彈性力和阻尼力，并認為其作用在一直條線上。

圖2 螺栓法蘭連接軸向力簡化模型 Fig.2 Qualitative relationship between axial force and the relative displacement z

由于連接結(jié)構螺栓受拉而法蘭盤受壓，因此其軸向的彈性作用力存在著拉壓不同的特性，定性函數(shù)如圖2(b)中虛線所示。工程中需要在螺栓連接部位施加一個軸向預緊力，以使連接部分的力學特性相對合理。預緊力過大會減小結(jié)構壽命，過小會使得結(jié)構相對運動超過標準，因此通常以圖中實線所代表的定性函數(shù)關系作為軸向彈性力，其中拉力的第二段剛度kT2=kc。簡化成力學模型得到分段線性剛度，在圖2(a)中用非線性彈簧表示。

1.2 螺栓法蘭連接的彎矩特性

圖2(a)所示的單個螺栓軸向力簡化模型代表螺栓法蘭連接的單個螺栓剖面。一個有n個性質(zhì)相同螺栓的連接結(jié)構具有n個如圖剖面。假設轉(zhuǎn)子彎曲量較小，則在截面上一周的n個螺栓連接結(jié)構受拉和受壓的連接區(qū)域數(shù)目相等，如圖3(a)所示。圖3(b)給出了由瞬時對稱位置固定坐標系中兩個螺栓組成的轉(zhuǎn)子連接結(jié)構剖面圖。

圖3 連接部位彎矩模型示意圖Fig.3 Model of the joint moment

左盤上下兩點受力方向相反，使連接部位產(chǎn)生彎矩作用。若連接結(jié)構系統(tǒng)受力處于瞬時平衡，則有上下兩點受力關系[9]:

(1)

式中：FTi為受拉端產(chǎn)生的拉力，F(xiàn)Ci為受壓端產(chǎn)生的壓力，d為螺栓法蘭連接的有效直徑。共有n/2個圖3(b)所示截面，產(chǎn)生n/2個彎矩。對于整個連接截面有:

(2)

航空發(fā)動機轉(zhuǎn)子為準剛性轉(zhuǎn)子，其渦動與橫向振動均不明顯。為了研究方便，把轉(zhuǎn)子連接部分受到的合彎矩分解到繞x軸和y軸兩個方向：

M=Mx+My

(3)

考慮繞x軸方向的彎矩示意圖如圖3(b)，則軸向力學特性使得連接結(jié)構拉壓變形并不相等，分別為:

(4)

(5)

式中：kT、kC分別為螺栓連接的等效軸向拉壓剛度。在小變形條件下，連接結(jié)構繞x軸轉(zhuǎn)角與軸向變形的關系可以表示為:

(6)

將式(4)、式(5)代入式(6)，整理得到螺栓連接彎曲角剛度表達式：

(7)

由于預緊力存在，等效拉伸剛度kT呈現(xiàn)分段線性規(guī)律，對應的kθ也呈現(xiàn)分段線性規(guī)律，如圖3(c)所示，與文獻[10]分析定性一致。由于每個螺栓性質(zhì)相同，則繞y軸彎矩模型可由同樣過程獲得，分段的臨界點對應的相對臨界轉(zhuǎn)角:

(8)

式中：Fpre為等效預緊力。

由上述分析可知，轉(zhuǎn)子系統(tǒng)發(fā)生彎曲的瞬時，連接部位的彎矩模型可由繞x、y軸作用的分段線性角剛度描述，沿連接截面一周的螺栓與固定坐標的四個交點為等效預緊力作用點。

1.3 含螺栓法蘭連接的轉(zhuǎn)子系統(tǒng)動力學建模

為了研究螺栓法蘭連接結(jié)構對轉(zhuǎn)子系統(tǒng)動力學的影響，基于下面假設建立轉(zhuǎn)子系統(tǒng)的動力學模型：

(1) 為體現(xiàn)連接部分的剛度對轉(zhuǎn)子系統(tǒng)的影響，將轉(zhuǎn)子軸考慮為剛性軸，支撐采用剛性支撐結(jié)構；

(2) 將高低壓葉盤簡化為兩個不同質(zhì)量的剛性盤。轉(zhuǎn)子系統(tǒng)的軸承支撐簡化為簡支約束。忽略彈性法蘭盤的質(zhì)量，連接結(jié)構的作用力看作兩葉盤的直接相互作用力；

(3) 兩段轉(zhuǎn)子之間的約束彎矩簡化為角剛度作用下，兩盤相互的耦合彎矩作用；

(4) 假設轉(zhuǎn)子橫向振動較小，螺栓約束轉(zhuǎn)子橫向振動的剛度簡化為線性剛度與線性阻尼；

(5) 忽略轉(zhuǎn)子軸向運動。

根據(jù)上述假設得到含螺栓連接假設(4)自由度轉(zhuǎn)子動力學模型如圖4所示，其廣義坐標為:

q=[θx1θy1θx2θy2]T

(9)

式中：θxi、θyi為別為兩段轉(zhuǎn)子剛體繞固定坐標系中x、y坐標軸的轉(zhuǎn)角，i=1,2。

圖4 含螺栓連接的轉(zhuǎn)子系統(tǒng)模型示意圖Fig.4 Model of the rotor system with bolt joints

這里需要說明，根據(jù)轉(zhuǎn)子動力學建模理論，每個剛性轉(zhuǎn)子運動需要用4個自由度表征，即xi、yi、θxi和θyi。單獨討論盤1的y方向，當盤1在y方向上的位移足夠小時，此時剛體繞固定坐標軸x的轉(zhuǎn)角也足夠小,滿足sinθ≈θ,此時y1=l1θx1。同理，x1、x2、y2也可由θy1、θy2、θx2表示。由于假設(1)，每盤運動使用兩個自由度即可表征。

采用拉格朗日能量法建立4個自由度的動力學微分方程：

(10)

式中:K為線性剛度矩陣。根據(jù)假設(5)，線剛度kl取定值。本文建模認為轉(zhuǎn)子系統(tǒng)為對稱系統(tǒng)，即兩段轉(zhuǎn)子x、y平動方向的剛度作用相等，取klx=kly=kl。同理，由clx=cly=cl組成阻尼矩陣C。G為陀螺矩陣，表示陀螺效應力對系統(tǒng)的影響。

角剛度取分段線性模型，繞x軸旋轉(zhuǎn)與繞y軸旋轉(zhuǎn)產(chǎn)生的角剛度相等，kθx=kθy=kθ(Δθ),其規(guī)律滿足如圖6的定性關系。微分方程(10)中的非線性剛度矩陣:

Knl(Δθ)=

(11)

式中：Δθx=θx1-θx2，Δθy=θy1-θy2，即非線性剛度矩陣Knl(Δθ)隨著系統(tǒng)的運動而變化。

2 含螺栓法蘭連接轉(zhuǎn)子系統(tǒng)的模態(tài)分析

2.1 螺栓法蘭連接對轉(zhuǎn)子渦動頻率的影響

由于陀螺效應的存在，隨著轉(zhuǎn)速增高渦動頻率逐漸分離，轉(zhuǎn)子系統(tǒng)發(fā)生渦動。求解動力學微分方程(10)的特征值問題，通過將物理坐標轉(zhuǎn)化為狀態(tài)坐標，可以得到狀態(tài)坐標式:

(12)

設式(12)的特解為q′=ψeλt，ψ為幅值列陣，將q′代入式(12)得到廣義特征值問題:

λPq′+Qq′=0

(13)

變形為:

Hq′=λq′

(14)

式中：

計算矩陣H的特征值與特征向量可得到轉(zhuǎn)子系統(tǒng)的n階渦動頻率[11]。系統(tǒng)的參數(shù)取值如表1所示。

表1 轉(zhuǎn)子系統(tǒng)簡化模型的參數(shù)Tab.1 Parameters of the rotor system

需要指出，由于阻尼力只在轉(zhuǎn)子系統(tǒng)不平衡響應中起到作用，在計算渦動頻率隨轉(zhuǎn)速變化的過程中不考慮。由于狀態(tài)空間理論只適用線性系統(tǒng)，角剛度矩陣Knl中的組成元素先后取角剛度分段的kθ1、kθ2形成兩個線性系統(tǒng)，來分別計算彎矩呈分段線性變化前后系統(tǒng)渦動頻率的值，對應等效預緊力Fpre取足夠大和取0的兩個極端值時系統(tǒng)的渦動頻率。以此線性化過程，得出一個范圍來定性預測渦動頻率的變化區(qū)間，并認為非線性剛度對應系統(tǒng)的實際情況在這個區(qū)間內(nèi)。使轉(zhuǎn)子轉(zhuǎn)速由0增長到1 800 rad/s，計算每個轉(zhuǎn)速下對應H矩陣的特征值問題，將各階特征值取虛部得到對應該轉(zhuǎn)速下的渦動頻率，繪制坎貝爾圖如圖5所示。

圖5 含螺栓連接轉(zhuǎn)子系統(tǒng)的坎貝爾圖Fig.5 Campbell diagram of the rotor system with bolt joints

圖5中斜線與頻率曲線的交點其橫坐標為臨界轉(zhuǎn)速?？梢钥闯?，角剛度取不同值時，第一階渦動頻率差異并不明顯，而第二階產(chǎn)生了較大差異。角剛度取kθ1時，系統(tǒng)的第二階渦動頻率較??；而角剛度取kθ2時系統(tǒng)第二階渦動頻率較大。第一階臨界轉(zhuǎn)速沒有明顯差異，為445 rad/s；由于兩線性系統(tǒng)第二階渦動頻率的差異，臨界轉(zhuǎn)速也產(chǎn)生了較大的差異，分別為910 rad/s與1 320 rad/s?？梢园l(fā)現(xiàn)角剛度取分段線性時，第一階臨界轉(zhuǎn)速不會受到角剛度變化的太大影響，而第二階臨界轉(zhuǎn)速出現(xiàn)較大差異，其取值在兩個線性系統(tǒng)計算的兩個臨界轉(zhuǎn)速之間。可以發(fā)現(xiàn)系統(tǒng)的渦動并不明顯，渦動頻率的分離較小。這是由于兩盤間的約束剛度足夠大，轉(zhuǎn)軸也簡化為剛性所導致；同時剛性盤的質(zhì)量取值較小，因此陀螺效應對于系統(tǒng)固有頻率的擾動較小，反映在圖5中就表現(xiàn)出無論角剛度取值如何，系統(tǒng)的渦動頻率隨轉(zhuǎn)速升高分離不明顯。這與實際發(fā)動機轉(zhuǎn)子系統(tǒng)的短粗結(jié)構和準剛性特性一致。

2.2 含螺栓法蘭連接轉(zhuǎn)子的振型

計算角剛度取kθ1時對應各個轉(zhuǎn)速下H矩陣特征值問題，可以得到各階特征向量，對應各個轉(zhuǎn)速下的轉(zhuǎn)子系統(tǒng)兩盤形心的振型，如圖6所示。

圖6 靜態(tài)與轉(zhuǎn)速為1 000 rad/s的振型Fig.6 Modes of the rotor at static and 1 000 rad/s

對于圖6結(jié)果需要指出，系統(tǒng)在靜態(tài)時未發(fā)生頻率分離，因此只有兩階振型；而在1 000 rad/s時轉(zhuǎn)子系統(tǒng)渦動頻率發(fā)生了分離，存在四階振型，但由圖5可知渦動頻率分離程度并不明顯，因此一二兩階、三四兩階振型十分接近，圖6中僅給出兩階。通過圖6中各盤的形心形成近似的圓形軌跡可以定性看出，轉(zhuǎn)子系統(tǒng)x,y方向的振型基本一致。這是因為兩盤間的約束在x和y以及繞x和y的方向上一致。從幅值來看，第一階振型兩盤幅值接近，第1盤略大于第2盤；第二階時第2盤幅值明顯大于第1盤，兩盤發(fā)生錯動。

3 含螺栓法蘭連接的轉(zhuǎn)子系統(tǒng)穩(wěn)態(tài)不平衡響應

由上節(jié)分析可以看出，當角剛度值發(fā)生變化時，轉(zhuǎn)子系統(tǒng)在第二階臨界轉(zhuǎn)速附近渦動頻率值不能確定，且兩盤發(fā)生錯動。本節(jié)引入轉(zhuǎn)子偏心不平衡量，采用四階Runge-Kutta方法求解方程組，對轉(zhuǎn)子系統(tǒng)進行數(shù)值計算，研究其最大幅值、臨界轉(zhuǎn)速的變化規(guī)律以及穩(wěn)態(tài)不平衡響應。

假設偏心在兩葉盤上，偏心距ε1=ε2=0.000 2 m，且偏心相位一致。加入轉(zhuǎn)子系統(tǒng)的不平衡激勵后，動力學微分方程(10)變?yōu)?

(16)

式中：f(t)為離心力激勵矩陣。

3.1 含螺栓法蘭連接轉(zhuǎn)子系統(tǒng)不平衡響應

取不同等效預緊力，以轉(zhuǎn)速為分岔參數(shù)，從0逐漸升高到1 800 rad/s，取dθx1/dt=0截面作分岔圖如圖7(a)～(e)所示，并作出等效預緊力分別取5 kN、10 kN、20 kN和40 kN時局部轉(zhuǎn)速范圍內(nèi)最大李雅普諾夫指數(shù)，如圖7(b)～(e)中小圖所示。

通過比較一系列分岔圖可以看出，轉(zhuǎn)子系統(tǒng)在第一階臨界轉(zhuǎn)速附近運動沒有明顯變化，但過第二階臨界轉(zhuǎn)速時發(fā)生分岔；隨著Fpre增大，第二階臨界轉(zhuǎn)速峰值逐漸不明顯甚至發(fā)生突跳；臨界轉(zhuǎn)速值也從910 rad/s(Fpre=2 kN)逐漸升高到1 300 rad/s(Fpre=40 kN)左右，變化范圍與圖5給出的范圍基本一致。

圖7 不同等效預緊力Fpre下第1盤y方向位移分岔圖Fig.7Bifurcationdiagramsofdisplacementsinthedirectionofyunderdifferentequivalentpre-tighteningforceFpre

Fpre取不同值時，系統(tǒng)位移響應總會在某一轉(zhuǎn)速下發(fā)生突跳，與文獻[6]的響應研究中在峰值附近發(fā)生的突跳類似；由圖7(b)～(e)中的最大李雅普諾夫指數(shù)可以看出，轉(zhuǎn)子系統(tǒng)在突跳點之后發(fā)生混沌運動，并且隨著轉(zhuǎn)速繼續(xù)升高發(fā)生倒分岔變?yōu)橹芷谶\動。不同的Fpre取值其對應的岔點以及倒分岔點均不同。當Fpre取值增大時，其突跳點向高轉(zhuǎn)速方向移動，而倒分岔點向低轉(zhuǎn)速方向移動，分岔轉(zhuǎn)速區(qū)間逐漸變小，如表2所示。

圖8為當Fpre=10 kN、轉(zhuǎn)速為1 200 rad/時盤1在y方向的穩(wěn)態(tài)響應。從圖7(c)以及圖8(b)中可以看出，盤1在該轉(zhuǎn)速下發(fā)生混沌運動；除兩階模態(tài)峰值外出現(xiàn)了3倍頻，但幅值較小；轉(zhuǎn)子的兩個橫向振動不再對稱，轉(zhuǎn)子此時運動軌跡變復雜。圖9為Fpre=20 kN、轉(zhuǎn)速為1 200 rad/s時盤1在y方向的穩(wěn)態(tài)響應。對比圖8可以看出，轉(zhuǎn)子系統(tǒng)運動更加復雜，頻譜圖中3倍頻幅值更高，且出現(xiàn)了5倍頻；軸心軌跡也更加不規(guī)則。這表明等效預緊力不僅使得混沌運動的轉(zhuǎn)速區(qū)間變小，而且使轉(zhuǎn)子系統(tǒng)運動更加復雜。通過圖7(c)和圖7(d)的比較，可以看出Fpre=20 kN對應的最大李雅普諾夫指數(shù)高于Fpre=10 kN一個量級，其混沌運動更加無序。

表2 不同預緊力下運動特性比較Tab.2 Comparison of motion characteristics under different pre-tightening force Fpre

圖8 Fpre=10 kN、1 200 rad/s時的穩(wěn)態(tài)響應Fig.8 Steady response at 1 200 rad/s, Fpre=10 kN

圖9 Fpre=20 kN、1 200 rad/s時的穩(wěn)態(tài)響應Fig.9 Steady response at 1 200 rad/s, Fpre=20 kN

3.2 兩盤的相對運動描述

轉(zhuǎn)子受到不平衡激勵后發(fā)生兩盤的相對錯動，螺栓法蘭連接參數(shù)化模型中兩盤相對轉(zhuǎn)角存在一個臨界值Δθ0，臨界值兩邊系統(tǒng)的角剛度發(fā)生變化。由第2節(jié)內(nèi)容可知，該臨界轉(zhuǎn)角的大小隨等效預緊力變化。

圖10 等效預緊力分別取2 kN、10 kN和20 kN時 不同轉(zhuǎn)速下的相對轉(zhuǎn)角Fig.10 Relatively rotation angle at diferent rotate speed when Fpre=2 kN、10 kN、20 kN

等效預緊力分別取2 kN、10 kN、20 kN，系統(tǒng)運動穩(wěn)定后兩盤相對轉(zhuǎn)角如圖10所示，圖中雙直線表示該等效預緊力下的相對臨界轉(zhuǎn)角。由于轉(zhuǎn)子橫向平動和彎曲存在關聯(lián)，圖中只給出了角度相對值。對比圖7(a)、圖7(c)和圖7(d)可以發(fā)現(xiàn)，相對轉(zhuǎn)角未超過臨界值的轉(zhuǎn)速下，盤1維持單周期運動；而當相對轉(zhuǎn)角超過臨界時，盤1發(fā)生混沌運動?？梢哉J為分段線性角剛度發(fā)生變化是轉(zhuǎn)子系統(tǒng)發(fā)生混沌運動的原因。比較圖10(b)和圖10(c)，F(xiàn)pre=10 kN時轉(zhuǎn)子系統(tǒng)穩(wěn)定運動時每一個周期相對臨界轉(zhuǎn)角都超過臨界值，而Fpre=20 kN存在某些周期未超過臨界值的情況，這使得轉(zhuǎn)子系統(tǒng)運動更加復雜。

3.3 含螺栓法蘭連接轉(zhuǎn)子系統(tǒng)的響應幅值

Fpre分別取0 kN、2 kN和預緊力足夠大的極限狀態(tài)作比較，在1 800 rad/s的范圍內(nèi)得到了盤1和盤2中心點的y方向幅頻曲線，如圖11所示。

圖11 轉(zhuǎn)子系統(tǒng)在不同F(xiàn)pre時的幅頻特性曲線Fig.11 Amplitude-frequency curve of rotor system under different Fprein direction y

比較Fpre=0、Fpre足夠大兩種極端情況，此時轉(zhuǎn)子系統(tǒng)為線性系統(tǒng)，角剛度取值不變。兩種狀態(tài)下第一階臨界轉(zhuǎn)速值(450 rad/s)重合，F(xiàn)pre=0對應的響應曲線幅值稍大；第二階臨界轉(zhuǎn)速發(fā)生了較大差異，F(xiàn)pre=0時系統(tǒng)的二階臨界轉(zhuǎn)速為910 rad/s，而Fpre取足夠大時為1 320 rad/s，幅值Fpre=0對應的曲線較小，與圖5一致。二階臨界轉(zhuǎn)速變小與文獻[7]得到的定性規(guī)律一致，但由于本文轉(zhuǎn)子系統(tǒng)模型突出連接結(jié)構影響，角剛度變化時二階臨界轉(zhuǎn)速發(fā)生的變化差異較大。

對比圖11(a)、(b)可以發(fā)現(xiàn)，第一階臨界轉(zhuǎn)速時，盤1較盤2的響應幅值大，而第二階臨界轉(zhuǎn)速下盤2比盤1幅值大，與圖6振型分析得到的定性規(guī)律一致。

對應圖10(a)，由圖11小圖可知，F(xiàn)pre=2 kN對應的響應幅值曲線在低轉(zhuǎn)速階段以及第一階臨界轉(zhuǎn)速附近與Fpre足夠大的曲線重合，此時兩轉(zhuǎn)子的相對轉(zhuǎn)角未突破臨界值；當相對轉(zhuǎn)角開始超過臨界值范圍時，在特定的轉(zhuǎn)速(700 rad/s)下幅值發(fā)生波動跳躍，最終與Fpre=0的響應幅值曲線重合。

令盤2偏心距ε2為0，逐步增大盤1的相對偏心距ε1，系統(tǒng)對應最大響應振幅如圖12所示?？梢钥闯鱿到y(tǒng)振動能量集中在第二階臨界轉(zhuǎn)速附近。聯(lián)系非線性動力學理論，能量集中使得第二階臨界轉(zhuǎn)速附近發(fā)生復雜運動；此時由于兩盤較大的相對轉(zhuǎn)角，非線性現(xiàn)象也趨于明顯。

圖12 不同偏心量兩階臨界轉(zhuǎn)速幅值比較Fig.12 Comparison of amplitude at critical rotation speed with different eccentricity quantity

4 結(jié) 論

本文在研究含螺栓法蘭連接轉(zhuǎn)子系統(tǒng)的振動機理過程中得出了以下結(jié)論:

(1)建立了剛性轉(zhuǎn)子螺栓法蘭連接作用彎矩的一個簡單模型，通過推導得到了轉(zhuǎn)子耦合彎矩模型，與之前學者得到的彎矩特性定性一致。

(2)應用狀態(tài)空間理論計算了動力學系統(tǒng)的特征值問題，給出了系統(tǒng)在動靜條件下的振型。發(fā)現(xiàn)轉(zhuǎn)子渦動頻率分離不明顯，且耦合角剛度函數(shù)的不光滑使轉(zhuǎn)子系統(tǒng)渦動頻率變小，二階臨界轉(zhuǎn)速值發(fā)生變化。定性給出了轉(zhuǎn)子兩階振型兩盤振動幅值之間的關系。

(3)研究了不同等效預緊力下含螺栓法蘭連接模型的轉(zhuǎn)子系統(tǒng)不平衡響應。轉(zhuǎn)子相對彎曲較大時，系統(tǒng)振動響應發(fā)生突跳，之后表現(xiàn)為混沌運動；預緊力增大，則分岔點后移，倒分岔點前移，混沌運動對應的轉(zhuǎn)速區(qū)間有變小趨勢，第二階臨界轉(zhuǎn)速由900 rad/s附近逐漸增加到1 300 rad/s；同時轉(zhuǎn)子系統(tǒng)的運動不穩(wěn)定加劇，混沌運動更加明顯。

(4)在系統(tǒng)能量集中在第二階臨界轉(zhuǎn)速附近，兩盤響應幅值大小有明顯差異。幅頻特性曲線規(guī)律隨著分段角剛度變化而發(fā)生變化，表現(xiàn)為低速段和高速段規(guī)律曲線不同的特性。

[1] 劉長福, 鄧明. 航空發(fā)動機結(jié)構分析 [M].西安：西北工業(yè)大學出版社，2006.[2] 楊帆, 趙普揚, 葛長闖. 無螺栓高壓壓氣機轉(zhuǎn)子結(jié)構分析[J]. 航空發(fā)動機，2012,38(3): 42-45. YANG Fang, ZHAO Puyang, GE Changchuang. Structure analysis of no bolted high pressure compressor rotor[J]. Aeroengine, 2012,38(3): 42-45.

[3] BOESWALD M, LINK M. Experimental and analytical investigations of non-linear cylindrical casing joints using base excitation testing[C]// Proceedings of International Modal Analysis Conference.Kissimmee, FL:(IMAC-XXI), 2003.

[4] BOESWALD M, LINK M. Identification of non-linear joint parameters by using frequency response residuals[C]// Proceedings of the 2004 International Conference on Noise and Vibration Engineering. Leuven: ISMA, 2004.

[5] GAO J, YUAN Q, LI P, et al. Effects of bending moments and pretightening forces on the flexural stiffness of contact interfaces in rod-fastened rotors[J]. Journal of Engineering for Gas Turbines and Power,2012, 134(10): 102503-102510.

[6] YUAN Q, GAO J, LI P. Nonlinear dynamics of the rod-fastened jeffcott rotor[J]. Journal of Vibration and Acoustics, Transactions of the ASME,2014, 136: 021011.

[7] QIN Z Y, HAN Q K, CHU F L. Analytical model of bolted disk-drum joints and its application to dynamic analysis of jointed rotor[J]. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 2014, 228(4): 646-663.

[8] LIU S G, MA Y H, ZHANG D Y, et al. Studies on dynamic characteristics of the joint in the aero-engine rotor system[J]. Mechanical Systems and Signal Processing, 2012, 29(5): 120-136.

[9] 欒宇.航天器結(jié)構中螺栓法蘭連接的動力學建模方法研究[D]. 大連：大連理工大學，2012.

[10] 王建民，鄭常良. 螺栓對接結(jié)構的非線性解析建模與分析[J].振動與沖擊,2013,32(20): 5-8. WANG Jianmin, ZHENG Changliang. Nonlinear analytical modeling and analysis for bolted joint structures[J]. Journal of Vibration and Shock, 2013,32(20): 5-8.

[11] 曹樹謙，張文德，蕭龍翔.振動結(jié)構模態(tài)分析-理論實驗與應用[M]. 天津：天津大學出版社，2014.

Vibration characteristics of rotor systems with bolt joints

LIU Zhuoqian1,2,3, CAO Shuqian1,2,3, GUO Hulun1,2,3, LI Liqing1,2,3

(1. Department of Mechanics, Tianjin University, Tianjin 300072, China;2. Tianjin Key Laboratory of Nonlinear Dynamics and Chaos Control, Tianjin 300072, China;3. State Key Laboratory of Engines, Tianjin University, Tianjin 300072, China)

Bolt joint structures are widely used in the aero-engine rotor systems. The mechanical model of a joint structure was set up and a parametric model for the bilinear moment distribution was given qualitatively. The modal characters and steady responses of the rotor system with bolt joint were investigated, by using the state space theory and numerical simulation. It is shown that the bilinear moment character has influence on modal characters and steady responses at medium rotational speed, especially at the vicinity of second order critical rotational speed, while it keeps steady and can be linearized at lower and higher rotational speeds. It is also found that pre-tightening force has a large influence on the stability of the rotor system. The increasing of pre-tightening force can intensify the instability of the system, but shorten the range of unstable motion.

rotor system; bolt joint; stability; nonlinear dynamics

國家973計劃資助項目(2015CB057400)；天津市自然科學基金項目(11JCZDJC25400)

2015-07-02 修改稿收到日期：2015-09-13

劉卓乾 男，碩士生，1990年生

曹樹謙 男，博士,教授,博士生導師，1964年生 E-mail: sqcao@tju.edu.cn

V231

A

10.13465/j.cnki.jvs.2016.22.002