外圈故障滾動(dòng)軸承周期運(yùn)動(dòng)Neimark-Sacker分岔研究

王 強(qiáng)，劉永葆，徐慧東，賀 星，劉樹勇

（1.海軍工程大學(xué) 動(dòng)力工程學(xué)院，武漢430033；2.湖南大學(xué) 機(jī)載與運(yùn)載工程學(xué)院，長(zhǎng)沙 410082）

外圈故障滾動(dòng)軸承周期運(yùn)動(dòng)Neimark-Sacker分岔研究

王 強(qiáng)1，劉永葆1，徐慧東2，賀 星1，劉樹勇1

（1.海軍工程大學(xué) 動(dòng)力工程學(xué)院，武漢430033；2.湖南大學(xué) 機(jī)載與運(yùn)載工程學(xué)院，長(zhǎng)沙 410082）

文章基于非線性理論，建立了三自由度非光滑系統(tǒng)軸承外圈故障模型，研究了該情況下，系統(tǒng)周期運(yùn)動(dòng)的Neimark-Sacker分岔現(xiàn)象和混沌等非線性行為。求出系統(tǒng)的切換矩陣，將得到的切換矩陣結(jié)合Floquet理論確定了該非光滑系統(tǒng)周期運(yùn)動(dòng)發(fā)生Neimark-Sacker分岔的條件。通過在碰撞面處建立Poincaré映射，用數(shù)值方法進(jìn)一步揭示軸承系統(tǒng)的周期運(yùn)動(dòng)經(jīng)Neimark-Sacker分岔通向混沌的現(xiàn)象。發(fā)現(xiàn)當(dāng)旋轉(zhuǎn)頻率接近臨界分岔點(diǎn)時(shí)，系統(tǒng)有一對(duì)Floquet特征乘子的模接近1，其余特征乘子模都小于1，系統(tǒng)發(fā)生Neimark-Sacker分岔，隨著旋轉(zhuǎn)頻率的增加，系統(tǒng)經(jīng)歷了典型的Neimark-Sacker分岔通向混沌的非線性行為。同時(shí)研究了不同的阻尼系數(shù)對(duì)系統(tǒng)分岔的影響，發(fā)現(xiàn)阻尼可以有效地延遲系統(tǒng)的分岔點(diǎn)。對(duì)該故障軸承系統(tǒng)分岔和混沌的研究，可為實(shí)際裝備安全運(yùn)行及故障診斷提供依據(jù)，同時(shí)為設(shè)計(jì)提供理論指導(dǎo)和技術(shù)支持。

軸承；Floquet理論；倍化分岔；混沌

0 引 言

隨著船舶事業(yè)的大發(fā)展，燃?xì)廨啓C(jī)作為艦船的主動(dòng)力正在逐步地發(fā)展，而目前燃?xì)廨啓C(jī)在運(yùn)行過程中極易出現(xiàn)故障，幾年前某船在海上航行過程中，出現(xiàn)了多級(jí)葉片損壞，軸承嚴(yán)重失效的故障，造成了重大損失。在燃機(jī)用來發(fā)電的設(shè)備使用中，國(guó)內(nèi)外也出現(xiàn)了多起軸承損壞造成嚴(yán)重事故的情況。據(jù)統(tǒng)計(jì)軸承作為一個(gè)旋轉(zhuǎn)核心部件發(fā)生故障的概率占所有設(shè)備故障的30%，因此，本文針對(duì)這種情況，從最根本的軸承碰撞出發(fā)，研究了軸承在運(yùn)轉(zhuǎn)過程中出現(xiàn)的分岔混沌等非線性振動(dòng)行為，為設(shè)備的軸承設(shè)計(jì)、健康管理和預(yù)測(cè)診斷提供重要的理論依據(jù)。而碰撞問題屬于強(qiáng)非線性、非光滑系統(tǒng)的一種，常用的線性化假設(shè)已經(jīng)不能滿足需要，必須考慮非線性、非光滑因素的的影響。在光滑的非線性領(lǐng)域已有了很多成熟的理論，Lyapunov、Bendix、Dulac和Birkoff等人討論的系統(tǒng)的穩(wěn)定性、極限環(huán)以及不變流形等問題，而Poincare、Andronov和Bogolubov等人的如小參數(shù)法、坐標(biāo)變形法、多尺度法、慢變參數(shù)法、KBM法、諧波平衡法、等效線性化方法等來逼近非線性系統(tǒng)的精確解[1]。隨著非線性理論的發(fā)展，出現(xiàn)了如混沌的結(jié)構(gòu)和普適性、分岔、分形、突變、奇異性、奇怪吸引子、混純同步等等[2-3]，形成了非線性動(dòng)力學(xué)中的主要概念以及理論基礎(chǔ)。而中心流行理論、分岔理論、奇異性理論、攝動(dòng)理論等光滑系統(tǒng)的非線性理論在非光滑系統(tǒng)中，這些理論并不能完全適用或做很大修改，而且出現(xiàn)了許多常規(guī)光滑系統(tǒng)沒有的特殊非線性現(xiàn)象，如角點(diǎn)碰撞、粘滯—滑移分岔、C—型奇怪吸引子等。

目前非線性理論在各個(gè)領(lǐng)域有了突破和應(yīng)用。國(guó)外，Shaw和Holmes[4]對(duì)一類在簡(jiǎn)諧激振力作用下有單側(cè)約束的單自由度振子做了研究，用中心流形定理分析了周期運(yùn)動(dòng)的局部分岔，并通過同宿相截條件討論了混沌運(yùn)動(dòng)。Peterka[5]研究了具有粘滯阻尼的碰撞振子中的擦邊分岔、周期倍化分岔和鞍結(jié)分岔之間的轉(zhuǎn)遷現(xiàn)象。Kleczka等人[6]通過對(duì)不穩(wěn)定周期解延續(xù)及胞映射方法發(fā)現(xiàn)并討論了含間隙振子中混沌的激變現(xiàn)象。Leine等人[7]對(duì)非光滑系統(tǒng)周期解的不連續(xù)分岔作了進(jìn)一步的研究，分析了伴隨基解矩陣的跳躍而發(fā)生的各種不連續(xù)分岔現(xiàn)象。Luo[8]研究了一個(gè)分段線性周期激勵(lì)系統(tǒng)，通過建立相應(yīng)映射，研究了各類穩(wěn)定和不穩(wěn)定的周期運(yùn)動(dòng)。國(guó)內(nèi)，在醫(yī)學(xué)方面，陳方方，洪靈等人[9]根據(jù)一類具有時(shí)滯及非線性特性發(fā)生率的SIRS傳染病模型，將非線性發(fā)生率的引入其中，運(yùn)用中心流形定理和規(guī)范型理論進(jìn)行研究分析，給出了分岔方向及分岔周期解穩(wěn)定性的計(jì)算公式，利用特征值理論分析了模型的Neimark-Sacker分岔行為及平衡點(diǎn)的局部穩(wěn)定性。夏小飛，王俊松[10]基于分岔理論，通過神經(jīng)群模型興奮性和抑制性突觸增益的余維一分岔分析，分別給出了神經(jīng)群運(yùn)行于單穩(wěn)、雙穩(wěn)、正常和異常極限環(huán)振蕩狀態(tài)的興奮性和抑制性突觸增益的單參數(shù)區(qū)間，進(jìn)而通過興奮性和抑制性突觸增益的余維二分岔分析給出了神經(jīng)群運(yùn)行于上述多種狀態(tài)的雙參數(shù)區(qū)域。在電力系統(tǒng)方面，王曉東，陳予恕[11]運(yùn)用多尺度法和C-L方法研究了單機(jī)無窮大電力系統(tǒng)在外部周期性負(fù)荷擾動(dòng)作用下的解析解及其穩(wěn)定性和主共振響應(yīng)在不同系統(tǒng)參數(shù)下的不同分岔模式，發(fā)現(xiàn)了系統(tǒng)具有倍周期分岔、混沌、增幅振蕩失步等豐富動(dòng)力學(xué)行為，可為電力系統(tǒng)中同步發(fā)電機(jī)使用提供理論指導(dǎo)。陳章耀等人[12]研究了自治與非自治電路系統(tǒng)在周期切換連接下的動(dòng)力學(xué)行為及機(jī)理，采用系統(tǒng)周期切換的方法分析了Rayleigh振子的不同參數(shù)下，兩子系統(tǒng)的周期振蕩行為，及切換系統(tǒng)隨參數(shù)變化的最大李雅普諾夫指數(shù)及對(duì)應(yīng)的分岔圖，研究了其復(fù)雜的非線性行為及振蕩機(jī)理。而研究最多的主要集中在非線性理論的研究和機(jī)械系統(tǒng)中存在各種非線性，鄭小武，謝建華等人[13]采用Floquet-Lyapunov理論將常系數(shù)系統(tǒng)的控制分岔行為的方法應(yīng)用于一類具有周期系數(shù)的力學(xué)微分系統(tǒng)，研究了其控制平衡點(diǎn)分岔行為的有效性，設(shè)計(jì)了相應(yīng)的控制器，將的Flip分岔行為和Neimark-Sacker分岔行為最后控制在周期一點(diǎn)。侯東曉等人[14]應(yīng)用多尺度法和奇異性理論分析一類具有三勢(shì)阱Mathieu-Duffing振子系統(tǒng)在非自治情況下的余維三分岔特性，獲得系統(tǒng)在Smale馬蹄意義下混沌的閾值，研究了系統(tǒng)的混沌行為和安全盆分岔。王學(xué)弟等人[15]基于非線性動(dòng)力系統(tǒng)的分岔控制理論，研究了一個(gè)新的非線性動(dòng)力系統(tǒng)的Neimark-Sacker分岔極限環(huán)幅值的控制問題，對(duì)應(yīng)本系統(tǒng)應(yīng)用具體的非線性控制器實(shí)現(xiàn)了Neimark-Sacker分岔極限環(huán)幅值的反饋控制，時(shí)得到了計(jì)算極限環(huán)幅值近似值的計(jì)算公式。吳鴻濤，張艷龍[16]研究了在隨機(jī)干擾條件下對(duì)兩自由度碰撞振動(dòng)系統(tǒng)倍化分岔的影響。柴林，吳曉明[17]建立了單自由度碰撞振動(dòng)機(jī)械的通用動(dòng)力學(xué)模型，研究了碰撞間隙、阻尼、剛度、激振頻率等參數(shù)對(duì)系統(tǒng)運(yùn)動(dòng)出現(xiàn)的分岔和混沌現(xiàn)象的影響。于海，陳予恕[18]針對(duì)航空發(fā)動(dòng)機(jī)低壓轉(zhuǎn)子系統(tǒng)，采用降維方法將26維裂紋故障的高維非線性動(dòng)力學(xué)模型降為2自由度的含有特征的低維非線性系統(tǒng)，得到了各種不同的分岔模式，反映了裂紋轉(zhuǎn)子的動(dòng)力學(xué)特征。茍向鋒等人[19]建立了三自由度單級(jí)直齒輪副傳動(dòng)系統(tǒng)的扭轉(zhuǎn)振動(dòng)非線性動(dòng)力學(xué)模型，結(jié)合系統(tǒng)相圖、Poincaré映射圖及FFT頻譜圖，分析了系統(tǒng)在激勵(lì)頻率變化時(shí)的動(dòng)力學(xué)特性。李曉貞，朱如鵬[20]基于集中參數(shù)理論，建立了正交面齒輪多自由度耦合振動(dòng)模型，采用龍格庫(kù)塔數(shù)值積分法對(duì)系統(tǒng)的動(dòng)力學(xué)方程求解，得到隨摩擦系統(tǒng)變換的系統(tǒng)動(dòng)態(tài)響應(yīng)分岔特性。而目前單獨(dú)針對(duì)軸承本身的非光滑、非線性振動(dòng)的研究還沒有。

軸承作為旋轉(zhuǎn)機(jī)械最核心的關(guān)鍵部件之一，軸承本身輕微異常振動(dòng)反應(yīng)在外圍的旋轉(zhuǎn)機(jī)械部件會(huì)有一個(gè)放大，尤其在高速旋轉(zhuǎn)的情況下，可能造成意想不到的后果，如燃機(jī)輪機(jī)等葉片的碰磨，振動(dòng)破壞，機(jī)械設(shè)備的損壞，甚至造成人員的傷亡。軸承一旦發(fā)生破損等故障時(shí)會(huì)造成更嚴(yán)重的后果，因此本文針對(duì)軸承本身，建立軸承三自由度故障模型，研究軸承非光滑、非線性振動(dòng)行為。

1 故障滾動(dòng)軸承建模分析及運(yùn)動(dòng)過程

圖1 滾動(dòng)軸承簡(jiǎn)化模型Fig.1 The simplified rolling bearing model

根據(jù)軸承的滾動(dòng)模型圖1（a）建立外圈存在故障的情況下軸承單個(gè)滾動(dòng)體的簡(jiǎn)化模型如圖1（b）所示。M1，M2，M3分別為內(nèi)圈與軸的等效質(zhì)量、滾動(dòng)體質(zhì)量、外圈及機(jī)座的質(zhì)量，X1，X2，X3分別為M1，M2，M3的運(yùn)動(dòng)的位移；C1，C2，C3別為M1與M2的阻尼，M2與M3之間阻尼，M3與固定端阻尼，K1，K2，K3分別為三個(gè)物體之間的剛度，K4為發(fā)生故障時(shí)，外圈與滾動(dòng)體之間的接觸剛度。F1sin（ Ω1T），F(xiàn)2sin（ Ω2T ）分別為作用在M1，M2的等效作用力，D為故障深度，Ω1與Ω2分別為轉(zhuǎn)軸和保持架的頻率，發(fā)生故障時(shí)Ω1與Ω2的關(guān)系為：

圖2 軸承模型二維相平面圖Fig.2 The two-dimensional phase plane of bearing model

內(nèi)圈轉(zhuǎn)動(dòng)頻率為fr=N/60（N為轉(zhuǎn)軸轉(zhuǎn)速）；外圈轉(zhuǎn)動(dòng)頻率fo=0；滾動(dòng)體個(gè)數(shù)為Nb；接觸角為α；軸承節(jié)圓直徑；滾動(dòng)體直徑d；外圈半徑R，內(nèi)圈半徑r，Ω1=fr，Ω2=fr/Nb。

為了描述該軸承系統(tǒng)的運(yùn)動(dòng)過程，引入一個(gè)分界面。首先定義邊界函數(shù)，E=X2-X3-D，分界面可表示如下：

表示物塊M2與M3剛接觸或分離，這樣狀態(tài)空間被分界面分成兩部分，如圖2所示。

根據(jù)上面的分析，可建立系統(tǒng)的運(yùn)動(dòng)學(xué)方程如下：

2 系統(tǒng)周期運(yùn)動(dòng)的Floquet特征乘子分析

將系統(tǒng)（3）和（4）寫為如下的規(guī)范式：

設(shè)系統(tǒng)（5）的一個(gè)解x（t）從區(qū)域v-出發(fā)，即x（ t0）∈v-。在t=tp時(shí)刻到達(dá)分界面Σ。系統(tǒng)在區(qū)間是連續(xù)的，相應(yīng)的基解矩陣也是連續(xù)的。然而由于向量場(chǎng)f（ t， x（t））在分界面處的非光滑性使得相應(yīng)的Jacobian矩陣在分界面處通常是不連續(xù)的，這將引起系統(tǒng)整個(gè)基解矩陣不連續(xù)，因此在不連續(xù)處需要求出相應(yīng)的切換矩陣。

下面來求分界面處的切換矩陣。

（1）從區(qū)域v-進(jìn)入?yún)^(qū)域v+時(shí)，對(duì)超平面Σ：e=x2-x3-d=0，有法向量n=［0，0，1，0，-1，0 ］T，設(shè)一周期解x（t）到達(dá)分界面Σ的時(shí)間為t1并交于點(diǎn)。在t1時(shí)刻計(jì)算切換矩陣如下

（2）從區(qū)域v+進(jìn)入?yún)^(qū)域v-，設(shè)周期解x（t）到達(dá)分界面Σ的時(shí)間為t2并交于點(diǎn)，在時(shí)刻t2有切換矩陣

下面求各光滑區(qū)域的基解矩陣。在區(qū)域v-系統(tǒng)的運(yùn)動(dòng)方程為

因?yàn)槲覀冄芯康氖窍到y(tǒng)（5）穿越了分界面的周期解，將切換矩陣（8）和（9）結(jié)合各子空間相應(yīng)的基解矩陣（13）和（14）經(jīng)過合成可得全局的單值矩陣：

（16）式中，T=2nπ/ω1為系統(tǒng)周期的整數(shù)倍，n為自然數(shù)；t0表示初始時(shí)刻，t1和t2分別表示軌線到達(dá)分界面的兩段時(shí)間。

于是系統(tǒng)的Floquet特征乘子即為單值矩陣（16）的特征值。對(duì)于系統(tǒng)（5）這樣的非光滑系統(tǒng)，由于分界面是光滑的，系統(tǒng)的Floquet特征乘子是連續(xù)穿越單位圓周的。當(dāng)一對(duì)Floquet特征乘子的模等于1，其它特征乘子仍位于單位圓內(nèi)時(shí)，系統(tǒng)（5）穩(wěn)定的周期解將發(fā)生Neimark-Sacker分岔。

3 數(shù)值分析

為了通過數(shù)值仿真進(jìn)一步揭示滾動(dòng)軸承系統(tǒng)（3）和（4）的倍化分岔通向混沌的現(xiàn)象，在分界面Σ處取Poincare截面如下：

其中：θ=ω1t；S=R（ mod 2π）為1個(gè)實(shí)數(shù)對(duì)2π取余數(shù)。

選取系統(tǒng)（3）和（4）的一組無量綱化參數(shù)：

以旋轉(zhuǎn)頻率ω1為分岔參數(shù)。

當(dāng)ω1=2.831時(shí)（ωs為軸承臨界分岔的旋轉(zhuǎn)頻率），軸承故障系統(tǒng)對(duì)應(yīng)的其中一對(duì)Floquet特征乘子為，接近單位圓周上的（1，0）點(diǎn)，其余特征值的模都小于1，位于單位圓內(nèi)，所有的特征值為：

由此可見軸承外圈故障系統(tǒng)在ω1=2.831時(shí)發(fā)生了Neimark-Sacker分岔。系統(tǒng)隨ω1變化的分岔圖如圖3所示。系統(tǒng)起初處于穩(wěn)定的單周期運(yùn)動(dòng)，如ω1=2.75時(shí)圖4所示的相圖和龐相萊截面圖。當(dāng)ω1=2.95時(shí)，系統(tǒng)經(jīng)過了分岔點(diǎn)ω1=2.831發(fā)生了Neimark-Sacker分岔，如圖5所示的相圖和龐加萊截面圖。研究了當(dāng)外圈m3發(fā)生Neimark-Sacker分岔時(shí)內(nèi)圈m1和滾動(dòng)體m2的振動(dòng)行為，其運(yùn)動(dòng)行為的龐相萊截面圖如圖6所示。同時(shí)研究了碰撞面處不同的阻尼系數(shù)對(duì)系統(tǒng)分岔點(diǎn)及相圖和龐相萊截面的影響，可以明顯看出，阻尼系數(shù)增加對(duì)系統(tǒng)的振動(dòng)有很大的衰減作用，隨著阻尼系數(shù)的增大，系統(tǒng)的分岔點(diǎn)會(huì)逐漸增大甚至消失如表1所示。

表1 不同的阻尼系數(shù)ζ3對(duì)系統(tǒng)分岔的影響Tab.1 The influence of different damping coefficients ζ3on the bifurcation of system

圖3 系統(tǒng)的分岔圖Fig.3 The bifurcation diagram of the system

圖4 ω1=2.75時(shí)單周期運(yùn)動(dòng)的相圖和龐加萊截面圖Fig.4 The phase diagram and poincare section diagram of single periodic motion at ω1=2.75

圖5 ω1=2.95時(shí)運(yùn)動(dòng)的相圖和龐加萊截面圖Fig.5 The phase diagram and poincare section diagram at ω1=2.95

圖6 ω1=2.95時(shí)內(nèi)圈m1和滾動(dòng)體m2的龐加萊截面圖Fig.6 The poincare section diagrams of inner ring m1and cage ball m2at ω1=2.95

當(dāng)選取系統(tǒng)（3）和（4）的無量綱化參數(shù)：

d=0.000 1；f1=15；f2=0；m2=1；m3=2；ζ1=0.1；ζ2=0.2；k4=5.5；k2=2；k3=2；基本參數(shù)保持不變，當(dāng)ζ3減小到ζ3=0.02時(shí)，以旋轉(zhuǎn)頻率ω1為分岔參數(shù)，依據(jù)Floquet理論及前面的公式推導(dǎo)，確定此時(shí)系統(tǒng)的分岔點(diǎn)。

當(dāng)ω1=2.772 4時(shí)（ωs為軸承臨界分岔的旋轉(zhuǎn)頻率），軸承故障系統(tǒng)對(duì)應(yīng)的其中一對(duì)Floquet特征乘子為，接近單位圓周上的（1，0）點(diǎn)，其余特征值的模都小于1，位于單位圓內(nèi)，所有的特征值為：

圖7 系統(tǒng)的分岔圖ζ3=0.02Fig.7 The bifurcation diagram of the system ζ3=0.02

由此可見軸承內(nèi)圈故障系統(tǒng)在ω1=2.772 4時(shí)發(fā)生了Neimark-Sacker分岔。系統(tǒng)隨ω1變化的分岔圖如圖7所示。

圖8 ω1=2.75時(shí)運(yùn)動(dòng)的相圖和龐加萊截面圖Fig.8 The phase diagram and poincare section diagram at ω1=2.75

圖9 ω1=2.9時(shí)運(yùn)動(dòng)的相圖和龐加萊截面圖Fig.9 The phase diagram and poincare section diagram at ω1=2.9

圖10 系統(tǒng)的分岔圖ζ3=0.04Fig.10 The bifurcation diagram of the system ζ3=0.04

同時(shí)選取系統(tǒng)（3）和（4）的無量綱化參數(shù)，基本參數(shù)保持不變，當(dāng)ζ3增大到ζ3=0.04時(shí)，以旋轉(zhuǎn)頻率ω1為分岔參數(shù)，依據(jù)Floquet理論及前面的公式推導(dǎo)，確定此時(shí)系統(tǒng)的分岔點(diǎn)。

當(dāng)ω1=2.879時(shí)（ωs為軸承臨界分岔的旋轉(zhuǎn)頻率），軸承故障系統(tǒng)對(duì)應(yīng)的其中一對(duì) Floquet特征乘子為，接近單位圓周上的（1，0）點(diǎn)，其余特征值的模都小于1，位于單位圓內(nèi)，該所有的特征值為：

圖11 ω1=2.85時(shí)運(yùn)動(dòng)的相圖和龐加萊截面圖Fig.11 The phase diagram and poincare section diagram at ω1=2.85

圖12 ω1=2.98時(shí)運(yùn)動(dòng)的相圖和龐加萊截面圖Fig.12 The phase diagram and poincare section diagram at ω1=2.98

由此可見軸承內(nèi)圈故障系統(tǒng)在ω1= 2.879時(shí)時(shí)發(fā)生了Neimark-Sacker分岔。系統(tǒng)隨ω1變化的分岔圖如圖10所示。

同時(shí)選取系統(tǒng)（3）和（4）的無量綱化參數(shù)，基本參數(shù)保持不變，當(dāng)ζ3增大到ζ3=0.1時(shí)，系統(tǒng)的分岔圖如圖13所示。

圖13 系統(tǒng)的分岔圖ζ3=0.1Fig.13 The bifurcation diagram of the system ζ3=0.1

4 結(jié) 論

本文針對(duì)軸承內(nèi)圈破損故障，建立了軸承三自由度分段非光滑模型，應(yīng)用F-loquet理論分析了該系統(tǒng)周期運(yùn)動(dòng)發(fā)生倍化分岔的條件。結(jié)果表明系統(tǒng)有1對(duì)Floquet特征乘子的模接近1，其余特征乘子的模小于1，系統(tǒng)發(fā)生了Neimark-Sacker分岔，數(shù)值仿真進(jìn)一步調(diào)查了系統(tǒng)由單周期經(jīng)過分岔點(diǎn)通向Neimark-Sacker分岔的過程；研究了當(dāng)外圈發(fā)生Neimark-Sacker分岔時(shí)對(duì)內(nèi)圈和滾動(dòng)體的影響發(fā)生了類似Neimark-Sacker分岔行為；同時(shí)通過研究不同的阻尼系數(shù)對(duì)系統(tǒng)分岔的影響，發(fā)現(xiàn)阻尼可以有效地延長(zhǎng)系統(tǒng)的分岔點(diǎn)，減少系統(tǒng)的分岔混沌等非線性行為。

致謝

非常感謝國(guó)家自然科學(xué)基金資助項(xiàng)目（51179197），國(guó)家海洋工程重點(diǎn)實(shí)驗(yàn)室（上海交通大學(xué)）基金資助項(xiàng)目（1009）。

[1]周紀(jì)卿,朱因遠(yuǎn).非線性振動(dòng)[M].西安:西安交通大學(xué)出版社,1997.

[2]Golubitsky,Schaeffer.Singularities and groups in bifurcation theory[M].Springer-Verlag,New York,1985.

[3]Kurths J,mayer-Kress.Nonlinear analysis of physiological data[M].Springer-Verlag,Berlin,1998.

[4]Shaw S W,Holmes P J.A periodically forced piecewise linear oscillator[J].Journal of Sound and Vibration,1983,90(1): 129-155.

[5]Peterka E.Bifurcation and transition phenomena in an impact oscillator[J].Chaos,Solitons and Fractals,1996,7(10): 1635-1647.

[6]Kleczka M,et a1.Nonlinear dynamics in engineering systems[M].Ed.By Schiehlen,Springer-Verlag,1990:141-148.

[7]Leine R I,Nijmeijer H.Dynamics and bifurcation of non-smooth mechanical systems[M].Berlin,Springer,2004:101-118.

[8]Luo A C J,Che L.Periodic motions and grazing in a harmonically forced,piecewise,linear oscillator with impacts[J]. Chaos,Solitons and Fractals,2005,24:567-578.

[9]陳方方,洪 靈.一類具有時(shí)滯和非線性發(fā)生率的SIRS傳染病模型穩(wěn)定性與Neimark-Sacker分岔分析[J].動(dòng)力學(xué)與控制學(xué)報(bào),2014,12(1):79-85. Chen Fangfang,Hong Ling.Stability and neimark-sacker bifurcation analysis of a delayed sirs epidemic model with nonlinear saturation incidence[J].Journal of Dynamic and Control,2014,12(1):79-85.

[10]夏小飛,王俊松.基于分岔理論的突觸可塑性對(duì)神經(jīng)群動(dòng)力學(xué)特性調(diào)控規(guī)律研究[J].物理學(xué)報(bào),2014,63(14): 140503-1:10. Xia Xiaofei,Wang Junsong.Influence of synaptic plasticity on dynamics of neural mass model:A bifurcation study[J]. Acta Phys.Sin.2014,63(14):140503-1:10.

[11]王曉東,陳予恕.一類電力系統(tǒng)的分岔和奇異性分析[J].振動(dòng)與沖擊,2014,33(4):1-6. Wang Xiaodong,Chen Yushu.Bifurcation and singularity analysis for a class of power systems[J].Journal of Vibration and Shock,2014,33(4):1-6.

[12]陳章耀,雪增紅,張 春,季 穎,畢勤勝.周期切換下Rayleigh振子的振蕩行為及機(jī)理[J]物理學(xué)報(bào),2014,63(1): 010504-1:8. Chen Zhangyao,Xue Zenghong,Zhang Chun,Ji Ying,Bi Qinsheng.Oscillation behaviors and mechanism of Rayleigh oscillator with periodic switches[J].Acta Phys.Sin,2014,63(1):010504-1:8.

[13]鄭小武,謝建華.一類周期系數(shù)力學(xué)系統(tǒng)分岔控制[J].西南交通大學(xué)學(xué)報(bào),2014,49(4):741-745. Zheng Xiaowu,Xie Jianhua.Bifurcation control of mechanical system with periodic coefficients[J].Journal of Southwest Jiaotong University,2014,49(4):741-745.

[14]劉 彬,趙紅旭,侯東曉.一類含三勢(shì)阱Mathieu-Duffing振子的相對(duì)轉(zhuǎn)動(dòng)系統(tǒng)的分岔和混沌[J].物理學(xué)報(bào),2014,63 (17):174502-1:9. Liu Bin,Zhao Hongxu,Hou Dongxiao.Bifurcation and chaos of some relative rotation system with triple-well Mathieu-Duffing oscillator[J].Acta Phys.Sin.,2014,63(17):174502-1:9.

[15]王學(xué)弟,彭 淼,楊天宇.一個(gè)新系統(tǒng)Neimark-Sacker分岔極限環(huán)幅值控制[J].工程數(shù)學(xué)學(xué)報(bào),2014,31(4):557-566. Wang Xuedi,Peng Miao,Yang Tianyu.Amplitude control of limit cycle from neimark-sacker bifurcation of a new system [J].Chinese Journal of Engineering Mathematics,2014,31(4):557-566.

[16]吳鴻濤,張艷龍.隨機(jī)干擾對(duì)兩自由度碰撞振動(dòng)系統(tǒng)的倍化分岔的影響[J].蘭州交通大學(xué)學(xué)報(bào),2014,33(3):59-63. Wu Hongtao,Zhang Yanlong.The effect of stochastic disturbance on doubling bifurcation of a TWO-DOF Vibro-impact system[J].Journal of Lanzhou Jiaotong University,2014,33(3):59-63.

[17]柴 林,吳曉明.機(jī)械碰撞振動(dòng)系統(tǒng)分岔與混沌的參數(shù)演化[J].廈門大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,53(4):508-513. Chai Lin,Wu Xiaoming.Evolution of bifurcation and chaos in mechanical vibro-impact system with parameters[J].Journal of Xiamen University(Natural Science),2014,53(4):508-513.

[18]于 海,陳予恕,曹慶杰.多自由度裂紋轉(zhuǎn)子系統(tǒng)非線性動(dòng)力學(xué)特性分析[J].振動(dòng)與沖擊,2014,33(7):92-98. Yu Hai,Chen Yushu.Cao Qingjie.Nonlinear dynamic behavior analysis for a cracked multi-DOF rotor system[J].Journal of Vibration and Shock,2014,33(7):92-98.

[19]茍向鋒,呂小紅,陳代林.單級(jí)齒輪傳動(dòng)系統(tǒng)的Hopf分岔與混沌研究[J].中國(guó)機(jī)械工程,2014,25(5):679-683,691. Gou Xiangfeng,Lu Xiaohong,Chen Dailin.Research on Neimark-Sacker Bifurcation and chaos of the single-stage gear transmission system[J].China Mechanical Engineering,2014,25(5):679-683,691.

[20]李曉貞,朱如鵬,李政民卿,靳廣虎.齒面摩擦對(duì)面齒輪傳動(dòng)系統(tǒng)振動(dòng)特性的影響分析[J].振動(dòng)工程學(xué)報(bào),2014,27 (4):583-588. Li Xiaozhen,Zhu Rupeng,Li Zhengminqing,Jin Guanghu.Influences of frictional coefficient on vibration characteristic of face-gear transmission system[J].Journal of Vibration Engineering,2014,27(4):583-588.

Neimark-Sacker bifurcation of rolling bearing system with fault in outer ring

WANG Qiang1,LIU Yong-bao1,XU Hui-dong2,HE Xing1,LIU Shu-yong1
(1.College of Power Engineering,Naval University of Engineering,Wuhan 430033,China;2.College of Mechanical and Vehicle Engineering,Hunan University,Changsha 410082,China)

Piecewise non-smooth model of three-degree-of-freedom rolling bearing system with fault in outer ring is established by the method of the nonlinear theory.The bifurcations and chaos of bearing system are studied.The switching matrixes of system are obtained at the switching boundaries,and the Neim-ark-Sacker bifurcation of non-smooth bearing system is analyzed by combining the switching matrixes with the Floquet theory for smooth systems.The numerical method is used to further reveal the bifurcations and chaos of bearing system through estabilshing the Poincare mapping on the collision plane.Results show that when the rotating frequency is decreased to a critical bifurcation point,a pair of complex conjugate Floquet multipliers is on the unit circle and others into a unit circle,and the Neimark-Sacker bifurcation appears.With the increase of rotating frequency,the system has experienced the nonlinear dynamical behaviors of classical Neimark-Sacker bifurcations to chaos.Also the influence of different damping coefficients on the bifurcation of system is analyzed and it is found that the damping of system can effectivelyreduce the nonlinear behaviors of bifurcation and chaos.The study of bifurcation and chaos of the fault bearing system provides reliable basis for the design and fault diagnosis and provides theoretical guidance and technical support for the actual design in the safe and stable operation of large high-speed rotating machinery.

bearing;floquet theory;period-doubling bifurcation;chaos

TH212 TH213.3

：Adoi:10.3969/j.issn.1007-7294.2017.05.013

1007-7294（2017）05-0621-12

2016-09-11

國(guó)家自然科學(xué)基金資助項(xiàng)目（51179197）；國(guó)家海洋工程重點(diǎn)實(shí)驗(yàn)室（上海交通大學(xué)）基金資助項(xiàng)目（1009）

王 強(qiáng)（1985-），男，博士研究生，E-mail：wangqiang13000306@163.com。