旋轉(zhuǎn)機械徑向碰摩模型分析

何 青，李 展，李 紅

（華北電力大學(xué) 能源動力與機械工程學(xué)院，北京 102206）

導(dǎo)致旋轉(zhuǎn)機械轉(zhuǎn)靜子碰摩發(fā)生的原因一般分為制造、安裝、運行三個方面。如在制造中由于技術(shù)不精而造成的不對稱或者偏心，在安裝過程中轉(zhuǎn)靜子之間的間隙偏小，在運行中升降負(fù)荷過快、暖機時間過短都會使轉(zhuǎn)靜子發(fā)生碰摩。而碰摩一般分為徑向、軸向和組合碰摩三類。從目前的文獻(xiàn)中可以看出徑向碰摩發(fā)生的概率最大。所謂徑向碰摩就是指轉(zhuǎn)子的外緣與靜止部位發(fā)生接觸而引起的碰摩。它一般發(fā)生在轉(zhuǎn)子的動葉圍帶汽封上、軸端汽封、隔板汽封、油擋、軸承、擋汽片合密封瓦等部位。同時徑向碰摩又可分成局部碰摩和全周碰摩。當(dāng)轉(zhuǎn)子在轉(zhuǎn)動中始終與靜子接觸是全周碰摩，否則為局部碰摩。相對于徑向碰摩來說由于軸向間隙一般比徑向間隙大，所以軸向碰摩和組合碰摩的發(fā)生率較低。

對于轉(zhuǎn)子與定子碰摩的研究主要包括理論實驗研究和現(xiàn)場故障排除研究兩類。其中理論實驗研究一般又包括碰摩模型的建立和基于模型基礎(chǔ)上的系統(tǒng)動力學(xué)行為分析兩方面。文獻(xiàn)［1，2］就研究了周期碰摩的渦動性和穩(wěn)定性。文獻(xiàn)［3］又研究了次諧周期碰摩運動的響應(yīng)特性。文獻(xiàn)［4］首先應(yīng)用動量守恒定理對碰摩的機理進(jìn)行了詳細(xì)闡述，然后又利用恢復(fù)系數(shù)，忽略轉(zhuǎn)靜子系統(tǒng)的彈簧阻尼力和系統(tǒng)碰摩后的位移提出碰摩的數(shù)學(xué)模型。文獻(xiàn)［5］也提出了適用于徑向局部碰摩的模型。文獻(xiàn)［6］等提出了單盤轉(zhuǎn)子系統(tǒng)的非線性動力學(xué)模型。文獻(xiàn)［7］建立了非線性彈性碰摩力模型。文獻(xiàn)［8］提出了非對稱轉(zhuǎn)子系統(tǒng)的碰摩力模型等等。目前碰摩的模型種類很多，但實際效果并不太理想。因為，一方面旋轉(zhuǎn)機械的種類繁多、結(jié)構(gòu)差異較大，且在不同部位（密封、軸承、葉片等）、以不同方式（徑向或軸向、單點或多點等）、不同力學(xué)特性的靜、動件（材料特性異同、支承特性差異等）的碰摩的模型不盡相同。另一方面系統(tǒng)碰摩模型中參數(shù)很多，所以受環(huán)境、溫度等的影響較大，以致造成模型的精度較低。模型的精度又將直接影響動力學(xué)行為的分析，造成碰摩故障無法識別和排除，給實際生產(chǎn)造成重大損失。

本文基于轉(zhuǎn)子與定子徑向碰摩發(fā)生時的系統(tǒng)結(jié)構(gòu)特征、發(fā)生部位和所含參數(shù)等不同，主要總結(jié)和分析了各類碰摩模型的優(yōu)缺點和使用范圍，并對將來模型的發(fā)展和研究提出了合理的建議。

1 徑向碰摩模型分類

由于大型轉(zhuǎn)子系統(tǒng)本身是精密且復(fù)雜的非線性系統(tǒng)，再加上轉(zhuǎn)定子的碰摩一般為其他故障的間接結(jié)果，所以碰摩這種故障的分析方法必定也是多樣的。對不同的機構(gòu)形式，從不同角度出發(fā)，可以建立不同的數(shù)學(xué)和物理模型。這些模型的使用對象和精度不盡相同。通過對目前使用的模型進(jìn)行總結(jié)可以得出如圖1所示各類模型。

圖1 徑向碰摩模型的分類

2 物理模型

物理模型是建立在分析現(xiàn)象與機理認(rèn)識基礎(chǔ)上的模型。描述的是對象系統(tǒng)“如何做”、“如何實現(xiàn)”的系統(tǒng)的物理過程。一個理想的物理模型的建立需要舍棄研究對象中的次要因素，抓住導(dǎo)致此問題的主要因素，以便用最簡單的方法解決問題。而轉(zhuǎn)靜子的系統(tǒng)結(jié)構(gòu)形式不同勢必會造成它的物理模型的差異。一般根據(jù)轉(zhuǎn)子的支承是否是剛性的而把系統(tǒng)分為剛性支承的物理模型［9，10］和彈性支承的物理模型［11］。

（1）基于剛性支承的物理模型

剛性支承的物理模型如圖2和圖3所示，也是研究碰摩問題的最原始的模型。此模型大大簡化了實際的轉(zhuǎn)子系統(tǒng)，所以分析起來比較簡單。

從圖3中可知當(dāng)轉(zhuǎn)子運動時的幾何中心偏移量O1O2的值大于轉(zhuǎn)靜子之間的間隙時碰摩就會產(chǎn)生。剛性支承結(jié)構(gòu)簡單，分析方便，但是若一味的加固支承將會得到相反的效果。且現(xiàn)代研究表明應(yīng)用彈性支承的物理模型更能反映現(xiàn)實的系統(tǒng)結(jié)構(gòu)。因此，剛性支承相對于彈性支承來說，適合精度要求不高的初步分析。

圖2 剛性支承的結(jié)構(gòu)圖

圖3 剛性支承的碰摩力圖

（2）基于彈性支承的物理模型

彈性支承的物理模型如圖4和圖5，定子的支承為彈性的，此時定子可以隨軸一起振動，它可以減少振幅，降低支承力，因此采用彈性支承是使超臨界機組運行安全可靠的有效設(shè)計途徑。從圖5中可以看出OO1和OO2分別為定子和轉(zhuǎn)子在運動時所產(chǎn)生的偏移量，當(dāng)兩個偏移量之差大于轉(zhuǎn)靜子之間的間隙時碰摩現(xiàn)象就會出現(xiàn)。

圖4 彈性支承的結(jié)構(gòu)圖

圖5 彈性支承的碰摩力圖

隨著振動力學(xué)的發(fā)展，文獻(xiàn)［8］等進(jìn)一步考慮了轉(zhuǎn)子和定子之間的阻尼力作用即讓彈性支承的物理模型包含了轉(zhuǎn)靜子在運動中的能量耗散作用，且從實驗結(jié)果來看這時的模型更符合實際。它的物理結(jié)構(gòu)模型圖如圖6所示。

圖6 彈性支承的阻尼結(jié)構(gòu)圖

3 力學(xué)模型

簡單而準(zhǔn)確有效的力學(xué)模型是對碰摩進(jìn)行理論分析的基礎(chǔ)。徑向碰摩的力學(xué)模型一般包括法向力模型和切向力模型［12-14］。

3.1 法向力模型

根據(jù)求碰摩法向力時所應(yīng)用的原理不同，法向力模型分為古典碰撞理論模型和動態(tài)接觸理論模型兩類。

3.1.1 古典碰撞理論模型

古典碰撞理論模型，就是運用動量守恒定理求出轉(zhuǎn)靜子之間的相互作用力。然后利用廣義沖量-動量方程和恢復(fù)系數(shù)方程來獲得碰撞后系統(tǒng)的運動狀態(tài)。但對多體系統(tǒng)碰撞問題釆用沖量-動量法時需要一些假設(shè)條件，且在非對心斜碰撞問題中，由于切向摩擦力的作用，僅采用恢復(fù)系數(shù)來表達(dá)系統(tǒng)前后速度的變化是不精確的。所以古典碰撞理論模型的使用有局限性。

3.1.2 動態(tài)接觸理論模型

動態(tài)接觸理論把碰撞過程微觀化，充分考慮碰撞過程中兩碰撞體之間的接觸時間、接觸變形和接觸力［12］。在此基礎(chǔ)上建立法向力模型時，不但要考慮是否完全表現(xiàn)了所要研究碰摩過程的主要特征，還要考慮能否有效分析轉(zhuǎn)靜子碰摩系統(tǒng)。動態(tài)接觸理論模型可以分為線性彈簧模型、基于Hertz理論模型、非線性彈簧阻尼模型。

（1）線性彈簧模型

當(dāng)轉(zhuǎn)軸或圓盤與靜子發(fā)生碰摩時，如果我們采用分段線性彈簧來模擬法向上的相互作用力則徑向法向力［15］為

式中e——轉(zhuǎn)子的徑向偏移量，mm；δ——碰撞面的接觸變形量，mm；k——碰摩面的剛度系數(shù)，N／m。

此模型是上述假設(shè)分段線性彈簧模型的近似描述方程，此方程雖然簡單但準(zhǔn)確性較差。且沒有考慮模型中線性彈簧的阻尼和剛度系數(shù)的選用等問題。文獻(xiàn)［16］提出了一種簡單常用的線性阻尼模型為

式中c——線性彈簧阻尼，N·s／m。

且在線形彈簧阻尼模型中要把碰撞過程看成一個單自由度振動模型。

（2）基于Hertz理論模型

轉(zhuǎn)子與定子碰摩系統(tǒng)是一類多參數(shù)、高維非光滑且強非線性系統(tǒng)，其動力學(xué)行為非常復(fù)雜。如果只在系統(tǒng)為線性的假設(shè)上建立模型，則難免產(chǎn)生誤差。假設(shè)碰摩系統(tǒng)滿足Hertz接觸理論條件，那么在此基礎(chǔ)上可以建立Hertz非線性碰摩力模型［17］。

文獻(xiàn)［18］提出當(dāng)在轉(zhuǎn)靜子間相對速度小于0.5m／s的情況下，徑向法向力可以表示為

式中δ——轉(zhuǎn)靜子的徑向嵌入深度，mm；E1，E2——轉(zhuǎn)靜子的楊氏模量，Pa；R1——轉(zhuǎn)子的半徑，mm；v1，v2——轉(zhuǎn)靜子的泊松比；h——轉(zhuǎn)靜子之間的間隙，mm。

文獻(xiàn)［17］經(jīng)過實驗驗證得出如圖7所示的Hertz非線性模型和線性剛度模型中徑向力Fn隨轉(zhuǎn)—定子嵌入深度之間的關(guān)系圖，在O-A段，Hertz非線性模型的法向力小于線性剛度法向力，在A點附近，此兩種模型的差別不太大，但隨著碰摩程度的加深，兩種模型的區(qū)別也增大。由于Hertz模型是以彈性接觸理論為基礎(chǔ)，通過引入Hertz接觸原理得出，因此在對嚴(yán)重碰摩的研究中使用Hertz碰摩力模型能減小誤差。

圖7 Hertz模型與線性模型的關(guān)系

（3）非線性彈簧阻尼模型

線性剛度模型和Hertz模型屬于建立在準(zhǔn)靜態(tài)接觸理論基礎(chǔ)上的模型，它需要建立在特殊的物理幾何條件下且無法計算系統(tǒng)能量損失，所以適用范圍較窄，不具備通用性。鑒于此有學(xué)者提出了一種改進(jìn)的Hertz理論模型即非線性彈簧阻尼模型，以求適應(yīng)更寬的范圍。

文獻(xiàn)［19］首先提出了同時考慮法向力和法向阻尼的法向力模型為

文獻(xiàn)［20］在式（4）的基礎(chǔ)上提出具有非線性法向力的模型為

根據(jù)Hertz理論，q一般取3／2。

非線性彈簧阻尼模型采用一個非線性彈簧-阻尼來模擬系統(tǒng)，彈簧用來表示接觸過程中的彈性變形，阻尼用來模擬能量耗散。由于非線性彈簧-阻尼模型即考慮轉(zhuǎn)、定子接觸過程中的塑性變形又提供了一種合理的確定法向接觸過程中能量消耗的方法，因此被很多學(xué)者廣泛應(yīng)用。

除此之外還有分段塑變模型和葉片與機箱碰撞力模型。分段塑變模型的特點是碰撞膨脹階段會產(chǎn)生不恢復(fù)的塑性變形，一般用于研究兩碰撞體之間有脆性材料墊存在時所用。這類模型的特點是碰撞膨脹階段與壓縮階段剛度存在一定的滯回，但不會產(chǎn)生塑性變形。葉片與機箱碰撞力模型主要用于飛機的雙轉(zhuǎn)子發(fā)動機的碰摩系統(tǒng)中。由于發(fā)動機中葉片尖端處有很高的線速度，碰撞能量較大，對轉(zhuǎn)子動力學(xué)行為的影響也大。所以需要建立一種專門的模型。以上兩種模型表達(dá)式都較為復(fù)雜，且應(yīng)用范圍較窄。

3.2 切向力模型

切向力模型就是碰摩產(chǎn)生的摩擦力的計算方法，摩擦可能阻礙轉(zhuǎn)子系統(tǒng)的運動，產(chǎn)生能耗。碰摩的切向力模型目前一般分為線性庫侖力、鬃毛模型、非線性Hertz模型、黏性模型等。

（1）線性庫侖力模型

線性庫侖摩擦力模型是計算切向摩擦力時最常用到的模型。因為它的形式簡單，且可以取得不錯的近似效果。在具體計算時如果轉(zhuǎn)靜子系統(tǒng)處于滑移狀態(tài)則為

如果處于黏合狀態(tài)則為

式中μm，μn——靜摩擦系數(shù)和動摩擦系數(shù)；Fn——法向力，N；vr——相對轉(zhuǎn)靜子的相對速度，m／s。

式（7）中動摩擦系數(shù)一般可以根據(jù)不同的因素來進(jìn)行估計。如果假設(shè)轉(zhuǎn)靜子的碰摩產(chǎn)生的熱全部轉(zhuǎn)化為熱能的話，文獻(xiàn)［21］提出了一種計算μn方法。文獻(xiàn)［22］也通過反復(fù)試驗得出一個較為統(tǒng)一的μn隨滑動速度變化的關(guān)系式。

（2）鬃毛模型

文獻(xiàn)［23，24］提出了一種計算切向力的鬃毛模型如圖8所示。此模型是用一組鬃毛來模擬轉(zhuǎn)靜子接觸表面上眾多的隨機分布的凸觸點。定義其摩擦力為

式中k——鬢毛剛度，N／m；S——鬢毛位移量，mm。

鬢毛模型是一種基于連續(xù)接觸觸動力學(xué)和大量的隨機行為模型，能很好地計算出摩擦力隨時間的變化，最重要的是切向摩擦力在接觸點黏合階段具有明確的物理意義和確定性。

文獻(xiàn)［25，26］在鬃毛模型的基礎(chǔ)上提出了LuGre模型如圖9所示。此模型采納了鬃毛模型建模的思想，但其與鬃毛模型不同，LuGre模型是基于鬃毛的平均變形來建模，而不是大量的隨機行為。LuGre模型的計算式為

式中k——剛性系數(shù)，N／m；z——鬢毛的均勻變形量，mm；c——阻尼系數(shù)，N（m／s）；μ——粘性摩擦力系數(shù)；v——兩表面之間的相對速度，m／s。

圖8 鬃毛模型

圖9 LuGre模型

LuGre模型能很好地描述轉(zhuǎn)靜子在尚未滑動時和靜摩擦與滑動摩擦轉(zhuǎn)換時的摩擦特點，因為它是利用鬃毛的平均變形來反映碰摩力的特征的，這一點比單純的鬢毛模型更有優(yōu)勢。除此它能完善的描述各種碰摩特性，如預(yù)滑動位移，可變的最大靜摩擦力、摩擦滯后現(xiàn)象等等。

（3）非線性Hertz模型

由于碰摩現(xiàn)象是一種強非線性故障，而庫倫模型等畢竟是一種對碰摩的簡化描述。則這樣的模型一定會給計算結(jié)果帶來一些影響。

如果不考慮轉(zhuǎn)靜子系統(tǒng)碰摩的熱效應(yīng)，再假設(shè)轉(zhuǎn)子與碰摩體間的局部碰撞變形為彈性變形，建立非線性 Hertz模型［27］為

式中μ——轉(zhuǎn)靜子的摩擦系數(shù)；Fn——法向力，N；vr——碰 摩 體 間 的 相 對 滑 動 速 度，m／s；B1，B2——相對速度影響系數(shù)，其值按系統(tǒng)的相關(guān)條件進(jìn)行假設(shè)。

非線性Hertz模型充分考慮了碰摩的強非線性，高維性及其它相關(guān)因素。所以說此模型可以剛好的和碰摩現(xiàn)象向匹配，精確性更高。

除此之外，還有兩類比較特殊的切向力模型。那就是葉片與機匣的模型［12］和渦輪機摩耗模型［13］，前一種模型在法向力模型中也做了簡要的介紹。由于葉片的高速運轉(zhuǎn)，當(dāng)它和定子碰撞時會產(chǎn)生高溫，以使葉片發(fā)生溶化進(jìn)而熔化金屬（黏性流體）就會產(chǎn)生相應(yīng)的剪切力。后一種模型是由于為了減少更換葉片的成本，而將可摩耗的材料作為機殼的內(nèi)壁材料以使葉片保持完好。但這樣就會產(chǎn)生能量耗散，并由此推出切向力模型。但這兩種模型使用對象較為專一。

3.3 系統(tǒng)動力學(xué)模型

由于研究的側(cè)重點以及基于上面各模型的不同，系統(tǒng)碰摩的動力學(xué)模型一般分為約束系統(tǒng)模型和分段光滑系統(tǒng)模型［12，13］。

3.3.1 約束系統(tǒng)模型

假設(shè)碰摩過程是瞬間完成的，且定子固定而又忽略轉(zhuǎn)靜子的彈性變形量，僅通過恢復(fù)系數(shù)去描述碰撞后的速度和能耗變化，此系統(tǒng)為約束微分系統(tǒng)，其相空間中的流在碰摩處是不連續(xù)的，其對應(yīng)的模型就是約束系統(tǒng)模型［28］。由于此模型不能直接反映轉(zhuǎn)靜子本身在碰摩面的物理參數(shù)變化，更適用于定性地研究碰摩系統(tǒng)可能出現(xiàn)的動力學(xué)現(xiàn)象，一次約束系統(tǒng)模型有局限性［29］。

3.3.2 分段光滑系統(tǒng)模型

如果認(rèn)為碰摩的發(fā)生持續(xù)一定的時間，且碰摩力連續(xù)但非光滑。考慮轉(zhuǎn)靜子碰摩時的能耗和變形，即當(dāng)轉(zhuǎn)子與定子接觸時，剛度增大，反之減小。則此時其在相空間中的流是連續(xù)但分段光滑的，與之相對應(yīng)的模型稱為分段光滑系統(tǒng)模型。為了對碰摩系統(tǒng)動力學(xué)模型進(jìn)行更詳細(xì)的說明，再按支承是剛性和彈性兩種類型對分段光滑系統(tǒng)模型進(jìn)行分類［13，14］。

（1）剛性支承系統(tǒng)模型

剛性碰摩的支承系統(tǒng)模型認(rèn)為定子是非慣性的，根據(jù)碰摩面是否具有嵌入型形變而把剛性支承模型分為變形和非變形碰摩模型兩類，文獻(xiàn)表明 后 者 的 應(yīng) 用 較 為 普 遍［30-32］。且 Wiercigroch［33，34］和文獻(xiàn)［12］都對此模型進(jìn)行了相關(guān)的研究與總結(jié)工作，此模型為

式中m——轉(zhuǎn)子質(zhì)量，kg；x，y——轉(zhuǎn)子水平和豎直方向

的位移，m；ks，kb——轉(zhuǎn)軸剛度和碰摩面（定子）的剛度，N／m；r0——轉(zhuǎn)子與定子之間的間隙，mm；r——轉(zhuǎn)子幾何中心的位移，mm；rd——碰撞處轉(zhuǎn)子的半徑，m；e——轉(zhuǎn)子的偏心量，mm；ω——轉(zhuǎn)子的轉(zhuǎn)速，r／min；ωw——轉(zhuǎn)子的渦動速度r／s；vr——碰摩處的相對線速度，m／s。當(dāng)發(fā)生碰摩時Θ取1。

此模型比起彈性模型的形式更為簡單，所以使用范圍更廣。此外還可以根據(jù)所研究的對象包含的具體的因素如是否考慮重力、裂紋、陀螺效應(yīng)、干摩擦失穩(wěn)等進(jìn)行適當(dāng)?shù)母淖儭?/p>

（2）彈性支承系統(tǒng)模型

彈性支承模型的實質(zhì)是考慮定子具有慣性，更進(jìn)一步可以根據(jù)碰摩面是否有嵌入型形變把彈性支承模型分為變形碰摩和非變形碰摩模型。由于慣性的存在系統(tǒng)的自由度數(shù)增加，在運算時非常復(fù)雜，所以在實際計算時很少選擇此模型。

4 討論

在徑向碰摩的物理模型中，彈性支承的物理模型考慮的因素即包含的參數(shù)較多，則它更符合實際的轉(zhuǎn)子系統(tǒng)。但剛性支承的物理模型所含參數(shù)較少，但計算較為簡單，所以更符合初步的計算和分析。在法向力模型中基于動態(tài)接觸的理論模型充分考慮了碰撞過程中兩碰撞體之間的接觸時間、接觸變形和接觸力等因素，所以它較古典碰撞理論模型更符合實際的系統(tǒng)，應(yīng)用也較廣泛。而在動態(tài)接觸的理論模型中非線性阻尼彈簧模型的應(yīng)用范圍和精度也較其他兩種模型更好。在切向力模型中線性庫侖力模型最簡單，可以用在精度不高的計算中，而鬢毛模型和非線性Hertz模型所含參數(shù)較多，相應(yīng)的精度更高，可以用在精確的分析中。在系統(tǒng)動力學(xué)模型中分段光滑系統(tǒng)模型較約束系統(tǒng)模型來說更符合系統(tǒng)的實際運行狀況。而分段光滑系統(tǒng)模型中的基于彈性支承的系統(tǒng)模型雖然在精度上比剛性支承的系統(tǒng)更接近實際情況，但由于其自由度數(shù)太大，計算繁瑣，所以剛性支承的系統(tǒng)模型的使用率更高。

（1）對于法向力碰摩模型，主要從經(jīng)典碰撞動力學(xué)理論出發(fā)，充分考慮動靜件的局部塑性變形、局部熱效應(yīng)和摩損等因素，其中包含有準(zhǔn)確給出各種法向碰摩力的正確形式以及選擇合理的符合現(xiàn)實的影響參數(shù)。把握主要影響參數(shù)去除次要的，從而構(gòu)建適合特定碰摩條件下的碰摩模型。

（2）對于切向力碰摩模型，主要以經(jīng)典庫侖摩擦理論為出發(fā)點，通過進(jìn)一步的細(xì)化分析摩擦?xí)r動靜件的狀態(tài)變化，改進(jìn)庫侖摩擦模型，研究在不同條件下的摩擦系數(shù)的精確選取，以求能準(zhǔn)確反應(yīng)摩擦過程中動靜件的不同運動狀態(tài)下（未接觸、黏合和滑移）的位移，速度以及受力的大小［35］。

（3）對于系統(tǒng)動力學(xué)模型，努力建立能充分反應(yīng)轉(zhuǎn)子和定子系統(tǒng)的各種的動力學(xué)特征的碰摩系統(tǒng)模型。這其中包含有合理的轉(zhuǎn)子和定子參與響應(yīng)時的模態(tài)變化，能夠精確反映系統(tǒng)發(fā)生局部單點或多點碰摩等。

（4）要充分考慮在碰摩過程系統(tǒng)周的接觸介質(zhì)的作用力變化（如油膜力、激振力、氣尖力等），除此還應(yīng)該考慮碰摩時產(chǎn)生的熱效應(yīng)和瞬時沖擊力大小等等。

［1］ H F Black.Interaction of a whirling rotor with a vibrating stator across a clearance annulus［J］.International Journal of Mechanical Engineering Science，1968，3（10）：1-12.

［2］ E F Enrich.The dynamics stability of rotor／stator radial rubs in rotating machinery ［J］.Journal of Engineering Industry，1969，4（4）：1025-1028.

［3］ D E Bently.Forced subrotative speed dynamic action of rotating machinery［J］.Wastepaper 1974，74-PET-16.

［4］ A.Muszynska.Rotor-to-stationary element rub-related vibration phenomena in rotating machinery literature survey［J］.Shock and Vibration Digest，1989，6（21）：3-11.

［5］ D W Childs.Fractional-frequency rotor motion due to nonsymmetrical clearance effects［C］.ASME 81-GT-145.

［6］ F L Chu，Z S Zhang.Bifurcation and chaos in a rub impact Jeffcott rotor system ［J］.Journal of Sound and Vibration，1998，210（1）：1-18.

［7］ 何韶君.非線性彈性碰摩轉(zhuǎn)子周期運動穩(wěn)定性分析［J］.機械科學(xué)與技術(shù)，2007，10（15）：13-17.HE Shao-jun，ZHANG Song-he.Stability analysis of periodic motion of nonlinear elastics rotor system with rubimpact fault［J］.Mechanical Science and Technology，2007，26（10）：1317-1320.

［8］ 吳敬東.轉(zhuǎn)子系統(tǒng)碰摩的若干非線性動力學(xué)問題研究［D］.沈陽：東北大學(xué)，2006.

［9］ W D Marscher.Test simulation of turbo machinery rotor／stator interactions［J］.Lubricate on engineering journal of the American society of lubrication Engineers，1983，39（9）：577.

［10］ 施維新.汽輪機發(fā)電機組振動及事故［M］.北京：中國電力出版社，2008.

［11］ 陶海亮，潘波.滾動軸承-轉(zhuǎn)子系統(tǒng)非線性動力響應(yīng)分析［J］.燃?xì)廨啓C技術(shù)，2013，26（1）：15-20.TAO Hai-liang，PAN Bo，GAO Qing，et al.Nonlinear dynamic response analysisof a rotor-ball bearing system［J］.Gas Turbine Technology，2013，26（1）：15-20.

［12］ 江俊，陳艷華.轉(zhuǎn)子與定子碰摩的非線性動力學(xué)研究［J］.力學(xué)進(jìn)展，2013，43（1）：132-148.JIANG Jun，CHEN Yan-hua.Advances in the research on nonlinear phenomona in rotor／stator rubbing systems［J］.Advances in Mechanics，2013，43（1）：132-148.

［13］ 董富祥，洪嘉振.多體系統(tǒng)動力學(xué)碰撞問題研究綜述［J］.力學(xué)進(jìn)展，2009，39（3）：352-358.DONG Fu-xiang，HONG Jia-zhen.Review of impact problem for dynamics of multibody system［J］.Advances in Mechanics，2009，39（3）：352-359.

［14］ 馬輝，太興宇，李煥軍，等.旋轉(zhuǎn)葉片-機匣碰摩模型及試驗研究綜述［J］.航空動力學(xué)報，2013，28（9）：2055-2069.MA Hui，TAI Xing-yu，LI Huan-jun，et al.Literature survey of rub-impact model and experiment of rotating blade and casing［J］.Journal of Aerospace Power，2013，28（9）：2055-2069.

［15］ 李永強，劉杰.非線性剛度不平衡轉(zhuǎn)子徑向碰摩動力學(xué)研究［J］.應(yīng)用力學(xué)學(xué)報，2005，22（3）：475-478.LI Yong-qiang，LIU Jie.Dynamics behavior of nonlinear rigid unbalanced rotor system supported on oil film bearing［J］.Chinese Journal of Applied Mechanics，2005，22（3）：475-478.

［16］ B K Mishra.A review of computer simulation of tumbling mills by the discrete element method：Part I-contact mechanics［J］.International Journal of Mineral Processing，2003，71，73-93.

［17］ 王濤，傅行軍.基于Hertz接觸模型的碰摩轉(zhuǎn)子動力學(xué)響應(yīng)研究［J］.汽輪機技術(shù)，2009，51（1）：25-27.WANG Tao，F(xiàn)U Xing-jun.Dynamic response of rub-impact rotor based on hertz contact model［J］.Turbine Technology，2009，51（1）.

［18］ S Timoshenko，J N Goodier.Theory of Elasticity［M］.McGraw Book Company，Inc，1951：377-379.

［19］ K H Hunt，F(xiàn) R Crossley.Coefficient of restitution interpreted as damping in vibrant impact［J］.Journal of Applied Mechanics，1975，42（4）：440-445.

［20］ M Fumagalli，G Schweitzer.Measurements on a rotor contacting its housing［C］.IMachE Conference Transactions 6th International Conference on Vibration in Rotating Machinery.1996，779-788.

［21］ Ettles.The thermal control of fiction at high sliding speeds［J］.Journal of Tribology，1986，108（1）：98-104.

［22］ 溫詩鑄，黃平.第3版.摩擦學(xué)原理［M］.北京：清華大學(xué)出版社，2008.

［23］ D A Haessig，F(xiàn)riesland.On the modeling and simulation of friction［J］.Journal of Dynamic systems，Measurement and control，1991，113：354-362.

［24］ Mao.Contact dynamics modeling for the simulation of the space station manipulators handing payloads［C］.IEEE International Conference on Robotics and Automation，1995：1252-1258.

［25］ W C Canudas，H Olsson，P Lischinsky.A new model for control of systems with friction［J］，IEEE Transactions on Automatic Control，1995，40（3）：419-425.

［26］ 劉占生，張敏，張廣輝，等.基于LuGre摩擦模型的葉片碰撞摩擦特性研究［J］.振動與沖擊，2012，31（12）：172-178.LIU Zhan-sheng，ZHANG Min，ZHANG Guang-hui，et al.Characteristics of impact-contact and friction between tips of blades based on LuGre model［J］.Journal of Vibration and Shock，2012，31（12）：172-178.

［27］ 王正浩，袁惠群，聞邦椿.相對速度對轉(zhuǎn)子碰摩運動特性的影響［J］.機械強度，2009，31（2）：175-181.WANG Zheng-hao，YUAN Hui-qun，WEN Bang-chun.Influence of relative velocity on motion properties of the rub-impact rotor system ［J］.Journal of Mechanical Strength，2009，31（2）：175-181.

［28］ 王濤.基于Hertz接觸模型的碰摩轉(zhuǎn)子動力學(xué)響應(yīng)研究［J］.汽輪機技術(shù)，2009，51（1）：13-16.WANG Tao，F(xiàn)U Xing-jun.Dynamic response of rub-impact rotor based on hertz contact model［J］.Turbine Technology，2009，51（1）：13-16.

［29］ 楊樹華.碰摩轉(zhuǎn)子系統(tǒng)的動力學(xué)行為研究［J］.風(fēng)機技術(shù)，2003，2（5）：12-16.YANG Shu-hua.Research on dynamics behaviour of impactive ＆frictive rotor system［J］.Compressor，Blower＆ Fan Technology，2003，2（5）：12-16.

［30］ 張思進(jìn)，陸啟韶.碰摩轉(zhuǎn)子系統(tǒng)的非光滑分析［J］.力學(xué)學(xué)報，2000，32（1）：59-69.ZHANG Si-jin，LU Qi-shao.A non-smooth analysis to the rubimpacting rotor system［J］.chinese journal of theoretical and applied mechanics，2000，32（1）：59-69.

［31］ E V Karpenko，M Wiercigroch，M P Cartmell.Piecewise approximate analytical solutions for a Jeffcott rotor with a snobbery ring［J］.International Journal of Mechanical Science，2002，44（5）：475-488.

［32］ E V Karpenko，M Wiercigroch，M P Cartmell.Regular and chaotic dynamics of a discontinuously nonlinear rotor system［J］.Chaos，Solitons ＆Fractals，2002，13：1231-1242.

［33］ E V Karpenko，M Wiercigroch，E E Pavlovskaia.Bifurcation analysis of a preloaded Jeffcott rotor.Chaos［J］.Solitons ＆Fractals，2003，15：407-416.

［34］ M Wiercigroch，E V Karpenko，M P Cartmell.Piecewise approximate analytical solutions for a Jeffcott rotor with a snubber ring［J］.International Journal of Mechanical Science，2002，44：475-488.

［35］ 秦志英，趙月靜，侯書軍.物料沖擊破碎過程中的一種非線性力模型［J］.振動與沖擊，2006，25（2）：35-37.QIN Zhi-ying，ZHAO Yue-jing，HOU Shu-jun.Nonlinear force model during material′s impact crushing［J］.Journal of Vibration and Shock，2006，25（2）：35-37.