屈鳴鶴， 吳少培， 俞力洋， 丁旺才， 李國(guó)芳， 黃 然

(蘭州交通大學(xué) 機(jī)電工程學(xué)院，蘭州 730070)

黏彈性隔振元件被廣泛應(yīng)用于機(jī)械設(shè)備中，如航空器的輔助動(dòng)力裝置隔振器[1]、軌道車輛二系懸掛的橡膠堆、汽車底盤中的懸架緩沖塊[2]、機(jī)床的減震墊等。其中橡膠只有在變形較小時(shí)才可近似看作線性彈性材料，超出這個(gè)范圍表現(xiàn)為非線性彈性。由于其特殊的結(jié)構(gòu)可以同時(shí)表現(xiàn)出黏彈性，能集緩沖、隔振、降噪等功能于一身，并且可根據(jù)實(shí)際需要來(lái)設(shè)計(jì)外形、剛度及阻尼[3]。為了更好地描述黏彈性隔振系統(tǒng)的力學(xué)特性，有學(xué)者提出三元件固體整數(shù)階模型，可以更好地反映黏彈性材料的松弛和蠕變特性[4]。學(xué)者們將可以等效橡膠材料的三參數(shù)力學(xué)模型稱為Zener模型和Ruzicka模型[5-7]。

分?jǐn)?shù)階微積分由于缺乏實(shí)際應(yīng)用背景使其長(zhǎng)期以來(lái)沒(méi)有得到研究和發(fā)展，近幾年分?jǐn)?shù)階微積分的定義、特性和計(jì)算才得以在工程領(lǐng)域應(yīng)用[8]。雖然傳統(tǒng)的整數(shù)階模型可以描述橡膠材料的力學(xué)特性，但不足以描述黏彈材料的頻率相關(guān)性。分?jǐn)?shù)階力學(xué)模型可以對(duì)黏彈性加以修正，進(jìn)而將整數(shù)階模型優(yōu)化為分?jǐn)?shù)階模型描述材料的本構(gòu)關(guān)系[9]。羅文波等[10]為準(zhǔn)確描述瀝青混合料的動(dòng)態(tài)黏彈性力學(xué)行為，在分?jǐn)?shù)階Zener模型的基礎(chǔ)上提出了改進(jìn)的分?jǐn)?shù)階Zener模型。Martin等[11]對(duì)黏彈性納米梁的分?jǐn)?shù)階動(dòng)力學(xué)行為進(jìn)行研究，借助拉普拉斯變換等方法對(duì)黏彈性特性進(jìn)行了研究分析。Bratu等[12]通過(guò)分?jǐn)?shù)階Zener模型描述了道路的復(fù)合土結(jié)構(gòu)，并分析了系統(tǒng)的動(dòng)態(tài)響應(yīng)。Lewandowski等[13]用分?jǐn)?shù)階Zener模型來(lái)反應(yīng)材料的流變特性。Ciniello等[14]分析了溫度對(duì)黏彈性材料的影響。

針對(duì)隔振系統(tǒng)的非線性特性，常宇健等[15]提出一種含有分?jǐn)?shù)階微分的金屬橡膠黏彈性本構(gòu)模型，在此基礎(chǔ)上建立了非線性動(dòng)力學(xué)模型，并結(jié)合實(shí)驗(yàn)驗(yàn)證了模型的準(zhǔn)確性。秦浩等對(duì)Caputo定義下的分?jǐn)?shù)階Duffing振子解析解和數(shù)值解的比較，驗(yàn)證了分?jǐn)?shù)階項(xiàng)化簡(jiǎn)為一階三角函數(shù)形式系統(tǒng)產(chǎn)生分岔和混沌的必要條件。孔凡等[16]通過(guò)諧波平衡法研究了簡(jiǎn)諧激勵(lì)下同時(shí)具有滯回特性和分?jǐn)?shù)階阻尼單元系統(tǒng)的穩(wěn)態(tài)響應(yīng)，采用不同方法求解系統(tǒng)的遲滯回線，發(fā)現(xiàn)諧波平衡法與逐步積分法得到的相關(guān)結(jié)果吻合較好。Zhen等[17]通過(guò)高次諧波平衡法求解了含負(fù)剛度幾何非線性系統(tǒng)的動(dòng)力學(xué)響應(yīng)，并對(duì)其隔振性能進(jìn)行了分析。余慧杰等[18]通過(guò)三次非線性函數(shù)描述金屬橡膠的非線性特性，用分?jǐn)?shù)階模型描述其黏彈性，所建立的分?jǐn)?shù)階非線性模型可以更加準(zhǔn)確描述橡膠動(dòng)態(tài)特性。

由于恢復(fù)力包含三次方非線性項(xiàng)，系統(tǒng)的動(dòng)力學(xué)響應(yīng)更加復(fù)雜，低頻區(qū)甚至存在超諧波共振。零部件的疲勞破壞與不同頻率的幅值相關(guān)，忽略低頻的幅值會(huì)對(duì)零部件的疲勞壽命設(shè)計(jì)帶來(lái)誤差[19]。對(duì)于隔振系統(tǒng)而言，準(zhǔn)確預(yù)測(cè)每個(gè)頻帶區(qū)間的動(dòng)力學(xué)響應(yīng)是有必要的。精密儀器設(shè)備經(jīng)常處于低頻微幅的振動(dòng)環(huán)境中，金屬橡膠隔振器是一種廣泛應(yīng)用于微振動(dòng)隔離的裝置[20-22]。隨著激勵(lì)頻率改變，系統(tǒng)的動(dòng)力學(xué)響應(yīng)也會(huì)發(fā)生變化，甚至出現(xiàn)分岔和混沌，進(jìn)而影響了系統(tǒng)的穩(wěn)定性，對(duì)系統(tǒng)的隔振性能以及疲勞壽命預(yù)測(cè)帶來(lái)影響。非線性系統(tǒng)還存在周期運(yùn)動(dòng)多樣性，運(yùn)動(dòng)轉(zhuǎn)遷過(guò)程更加復(fù)雜，甚至存在多態(tài)共存現(xiàn)象，這影響了低頻區(qū)隔振系統(tǒng)動(dòng)力學(xué)響應(yīng)的準(zhǔn)確預(yù)測(cè)。針對(duì)整數(shù)階非線性系統(tǒng)的動(dòng)態(tài)特性研究已相對(duì)成熟，而分?jǐn)?shù)階非線性系統(tǒng)的周期運(yùn)動(dòng)多樣性及轉(zhuǎn)遷規(guī)律有待揭示[23-26]。

為了揭示分?jǐn)?shù)階非線性隔振系統(tǒng)的動(dòng)力學(xué)響應(yīng)，本文對(duì)系統(tǒng)的超諧波共振和周期運(yùn)動(dòng)多樣性研究分析。首先對(duì)分?jǐn)?shù)階項(xiàng)進(jìn)行化簡(jiǎn)處理，其次采用高階諧波平衡法求解系統(tǒng)的穩(wěn)態(tài)響應(yīng)，數(shù)值仿真系統(tǒng)的動(dòng)力學(xué)性能并對(duì)近似解析結(jié)果進(jìn)行比較。接著建立兩個(gè)Poincaré映射描述系統(tǒng)周期運(yùn)動(dòng)的多樣性及周期運(yùn)動(dòng)的轉(zhuǎn)遷規(guī)律，采用Floquet理論對(duì)分岔類型加以判定。最后對(duì)超諧波共振和周期運(yùn)動(dòng)多樣性之間的關(guān)系分析研究，進(jìn)一步揭示系統(tǒng)參數(shù)對(duì)幅頻特性、分岔及混沌的影響。

1 模型建立及響應(yīng)求解

1.1 分?jǐn)?shù)階非線性Zener模型建立及分?jǐn)?shù)階項(xiàng)處理

采用分?jǐn)?shù)階非線性Zener模型描述橡膠的黏彈性如圖1所示，其中M為系統(tǒng)的質(zhì)量；F為外激勵(lì)幅值；X1為質(zhì)量塊的位移；X2為節(jié)點(diǎn)的位移；Ω為外激勵(lì)頻率；K1為線性彈性恢復(fù)力剛度系數(shù)；C為黏性阻尼系數(shù)；FK為三次非線性彈性恢復(fù)力；D為分?jǐn)?shù)階微分項(xiàng)；K為分?jǐn)?shù)階項(xiàng)的系數(shù)。引入外激勵(lì)幅值F的對(duì)照參數(shù)進(jìn)Fs，進(jìn)行如下變量代換

圖1 分?jǐn)?shù)階非線性Zener模型Fig.1 Fractional nonlinear Zener model

分?jǐn)?shù)階非線性Zener模型無(wú)量綱微分方程為

(1)

(2)

Γ(y+1)=yΓ(y)

(3)

式中,Γ(n)為Gamma函數(shù)，滿足式(3)。設(shè)質(zhì)量塊的穩(wěn)態(tài)響應(yīng)為

(4)

將式(4)代入式(2)并引入文獻(xiàn)[24]中的公式

(5a)

(5b)

取分?jǐn)?shù)階項(xiàng)的一階近似化簡(jiǎn)可得

(6)

式(6)所得分?jǐn)?shù)階項(xiàng)不僅具有阻尼作用也具有剛度作用與很多文獻(xiàn)結(jié)果相同，其中文獻(xiàn)[16]通過(guò)參數(shù)識(shí)別驗(yàn)證了分?jǐn)?shù)階項(xiàng)采用一階三角函數(shù)表示在實(shí)際工程應(yīng)用的合理性。本文主要研究隔振系統(tǒng)高次諧波幅值相比基波幅值是小量的情況，進(jìn)而分?jǐn)?shù)階化簡(jiǎn)忽略高階項(xiàng)，只取一階近似。

1.2 等效線性剛度系數(shù)和等效阻尼系數(shù)分析

將分?jǐn)?shù)階項(xiàng)化簡(jiǎn)所得一階三角函數(shù)式(6)代入到式(1)中，系統(tǒng)的無(wú)量綱微分方程變?yōu)?/p>

βx1(t)]=fcos(ωt)

(7)

系統(tǒng)等效線性剛度系數(shù)Keq和等效阻尼系數(shù)Ceq為

(8)

(9)

由式(8)和式(9)可知，分?jǐn)?shù)階的系數(shù)λ、分?jǐn)?shù)階的階數(shù)p對(duì)等效線性剛度系數(shù)Keq和等效阻尼系數(shù)Ceq都有一定影響。下面對(duì)系統(tǒng)的等效線性剛度系數(shù)Keq和等效阻尼系數(shù)Ceq進(jìn)行研究。選取系統(tǒng)結(jié)構(gòu)參數(shù)ξ=0.1，μk=2，ε=0.2，f=5，取分?jǐn)?shù)階的階數(shù)p=0.5。由圖2(a)可以看出分?jǐn)?shù)階系數(shù)λ越大，等效線性剛度系數(shù)Keq越大，且隨頻率ω的增大而增大。分?jǐn)?shù)階系數(shù)λ越小，等效阻尼系數(shù)Ceq越小，且隨頻率ω的增大在低頻區(qū)快速減小，離開(kāi)低頻區(qū)等效阻尼系數(shù)Ceq無(wú)明顯變化。其次取分?jǐn)?shù)階系數(shù)λ=0.01，由圖2(b)可以看出頻率較小時(shí)，分?jǐn)?shù)階項(xiàng)的階數(shù)p越小，等效線性剛度系數(shù)Keq越大；而頻率ω較大時(shí)，分?jǐn)?shù)階的階數(shù)p越小，等效線性剛度系數(shù)Keq越小；等效線性剛度系數(shù)Keq隨著頻率ω增大而增大。在超低頻區(qū)，分?jǐn)?shù)階的階數(shù)p越小，等效阻尼系數(shù)Ceq越大；離開(kāi)超低頻區(qū)，分?jǐn)?shù)階的階數(shù)p越大，等效阻尼系數(shù)Ceq越大。在低頻區(qū)等效阻尼系數(shù)Ceq隨著頻率ω的增大迅速減小，隨著頻率ω繼續(xù)增大離開(kāi)低頻區(qū)，等效阻尼系數(shù)Ceq無(wú)明顯變化。

(a) 系數(shù)λ與等效線性剛度系數(shù)和等效阻尼系數(shù)的關(guān)系

(b) 階數(shù)p與等效線性剛度系數(shù)和等效阻尼系數(shù)的關(guān)系圖2 分?jǐn)?shù)階項(xiàng)與等效線性剛度系數(shù)和阻尼系數(shù)的關(guān)系Fig.2 Relationship between fractional order term and equivalent linear stiffness and damping

不同分?jǐn)?shù)階的系數(shù)λ和分?jǐn)?shù)階的階數(shù)p會(huì)改變系統(tǒng)的等效阻尼系數(shù)Ceq和等效線性剛度系數(shù)Keq，進(jìn)而會(huì)改變系統(tǒng)的固有頻率。

1.3 諧波平衡法求解系統(tǒng)穩(wěn)態(tài)響應(yīng)

設(shè)節(jié)點(diǎn)的穩(wěn)態(tài)響應(yīng)為

(10)

其中(2n-1)表示高階諧波的階次。當(dāng)n=1時(shí)根據(jù)式(4)和式(10)可得質(zhì)量塊和節(jié)點(diǎn)的一階諧波響應(yīng)為

x1=A1cos(ωt)+B1sin(ωt)=A11sin(ωt+θ1)

(11a)

x2=a1cos(ωt)+b1sin(ωt)=a11sin(ωt+φ1)

(11b)

將式(11)代入到式(1)中進(jìn)行一階諧波平衡，可得質(zhì)量塊與節(jié)點(diǎn)幅值和相位的關(guān)系

(12)

將式(12)代入到式(11a)中可得到質(zhì)量塊的一階穩(wěn)態(tài)響應(yīng)

(13)

將式(13)和式(11b)代入式(7)可得

(14)

(15)

進(jìn)而可以求解出節(jié)點(diǎn)的一階幅值與相位

根據(jù)以上方法獲取的水體信息如圖1所示，同時(shí)根據(jù)Landsat OLI成像的準(zhǔn)確時(shí)間，借助于磨盤山水庫(kù)同一時(shí)間的水文觀測(cè)資料、水庫(kù)水位和流量等水文數(shù)據(jù)，基于庫(kù)容曲線計(jì)算出當(dāng)時(shí)的庫(kù)區(qū)水面面積值為22.012 km2，并將不同方法所獲取的水庫(kù)面積與實(shí)測(cè)值進(jìn)行對(duì)比，如表1所示。

(16)

(17)

系統(tǒng)的高次超諧波響應(yīng)需先求解節(jié)點(diǎn)處的各階諧波幅值系數(shù)，再選取時(shí)間歷程圖中節(jié)點(diǎn)的最大位移作為幅值，進(jìn)而得到質(zhì)量塊的高次響應(yīng)。當(dāng)n=2時(shí)采用三階諧波平衡法設(shè)節(jié)點(diǎn)處的響應(yīng)為

(18)

對(duì)式(18)進(jìn)行三階諧波平衡可得如下方程組。

(19a)

(19b)

(19c)

36εω3a1b1b3ξ3-30εω2a1a3b1ξ2-3εωa1b1b3ξ+

(19d)

系統(tǒng)參數(shù)確定的條件下借助計(jì)算機(jī)求解并對(duì)結(jié)果進(jìn)行數(shù)據(jù)處理，即可得到節(jié)點(diǎn)幅值系數(shù)a1、b1、a3、b3的值[17]，再將節(jié)點(diǎn)幅值系數(shù)代入式(1)可得質(zhì)量塊的幅值系數(shù)

(20)

進(jìn)而得到質(zhì)量塊的三次穩(wěn)態(tài)響應(yīng)。

(21)

根據(jù)三階諧波平衡法的求解過(guò)程以此類推可得系統(tǒng)的高次諧波響應(yīng)。

1.4 數(shù)值方法對(duì)比

隨機(jī)選取一組系統(tǒng)參數(shù)：μk=0.5，f=5，ξ=0.1，ε=0.1，λ=0.1。采用四階Runge-Kutta法對(duì)系統(tǒng)數(shù)值求解，通過(guò)動(dòng)力學(xué)仿真軟件UM對(duì)隔振系統(tǒng)進(jìn)行虛擬實(shí)驗(yàn)仿真，不同方法求解系統(tǒng)動(dòng)力學(xué)響應(yīng)結(jié)果接近，如圖3所示。

(a) 質(zhì)量塊瞬態(tài)響應(yīng)

(b) 質(zhì)量塊穩(wěn)態(tài)響應(yīng)圖3 質(zhì)量塊時(shí)間歷程圖Fig.3 Time history diagram of mass block

一階諧波平衡法只能求解系統(tǒng)的主共振，無(wú)法求解系統(tǒng)的超諧波共振，為進(jìn)一步研究系統(tǒng)在低頻區(qū)的動(dòng)力學(xué)響應(yīng)，采用三階諧波平衡法對(duì)系統(tǒng)進(jìn)行求解進(jìn)而得到超諧波共振如圖4所示。

圖4 不同方法求解系統(tǒng)的幅頻響應(yīng)曲線Fig.4 Amplitude-frequency response curve by different methods

2 系統(tǒng)運(yùn)動(dòng)狀態(tài)與分岔類型

3 超諧波共振下的周期運(yùn)動(dòng)多樣性

為方便描述系統(tǒng)運(yùn)動(dòng)規(guī)律，定義符號(hào)T-N來(lái)描述系統(tǒng)的運(yùn)動(dòng)狀態(tài)，其中T為系統(tǒng)運(yùn)動(dòng)的周期數(shù)，N為周期運(yùn)動(dòng)內(nèi)簡(jiǎn)諧振動(dòng)的次數(shù)。為避免混沌出現(xiàn)，選取阻尼系數(shù)相對(duì)較大的基準(zhǔn)參數(shù)：μk=0.15，ε=0.6，p=0.5，λ=0.1，f=5，ξ=0.5。選取n=5采用高階諧波平衡法可得系統(tǒng)質(zhì)量塊的幅頻曲線如圖5所示。

(a) 主共振及超諧波共振

(b) 超諧波共振局部放大圖圖5 幅頻特性曲線圖Fig.5 Amplitude-frequency response curve

從圖5中可見(jiàn)高次超諧波共振幅值相比較主共振幅值非常小，但是在掃頻的過(guò)程中高次超諧波也存在跳躍現(xiàn)象如圖5(b)所示，系統(tǒng)在低頻區(qū)也存在多態(tài)共存現(xiàn)象且動(dòng)力學(xué)行為更加復(fù)雜。隨著激勵(lì)頻率改變，系統(tǒng)的動(dòng)力學(xué)響應(yīng)出現(xiàn)分岔和混沌，這影響了隔振系統(tǒng)的穩(wěn)定性，對(duì)系統(tǒng)的動(dòng)力學(xué)響應(yīng)預(yù)測(cè)帶來(lái)影響。隨著頻率ω的降低，數(shù)值模擬系統(tǒng)的運(yùn)動(dòng)狀態(tài)規(guī)律：1-1→1-2→1-1→1-3→1-5→1-7→1-9，如圖6所示。

(a) 定相位面

(b) 定極大位移面圖6 系統(tǒng)質(zhì)量塊的運(yùn)動(dòng)狀態(tài)Fig.6 Motion state of system mass block

對(duì)比超諧波共振的階次可見(jiàn)超諧波共振的階次越高，周期運(yùn)動(dòng)內(nèi)的簡(jiǎn)諧振動(dòng)次數(shù)越多，周期運(yùn)動(dòng)內(nèi)簡(jiǎn)諧振動(dòng)的次數(shù)與超諧波共振的階次近似一致。隨著超諧波次數(shù)的增加，幅頻曲線跳躍的距離減小且跳躍方向與分岔圖中SNB突變的方向一致。由于設(shè)解的形式均為奇數(shù)次諧波，無(wú)論諧波平衡法設(shè)其解的階次為多大，該方法都具有一定局限性。從圖7相圖可以看到在ω<1的低頻區(qū)，隨著頻率ω減小系統(tǒng)始終保持周期一運(yùn)動(dòng)，但是周期內(nèi)的簡(jiǎn)諧振動(dòng)次數(shù)增加。

(a) ω=1的Poincaré截面圖和相圖(1-1運(yùn)動(dòng))

(b) ω=0.696 1的Poincaré截面圖和相圖(1-3運(yùn)動(dòng))

(c) ω=0.440 6的Poincaré截面圖和相圖(1-5運(yùn)動(dòng))

(d) ω=0.328 6的Poincaré截面圖和相圖(1-7運(yùn)動(dòng))

(e) ω=0.230 6的Poincaré截面圖和相圖(1-9運(yùn)動(dòng))圖7 系統(tǒng)的Poincaré截面圖和相圖Fig.7 Poincaré section diagram and phase diagram of the system

由于Duffing系統(tǒng)的恢復(fù)力為奇函數(shù)，這也導(dǎo)致系統(tǒng)的運(yùn)動(dòng)狀態(tài)在相圖中也呈現(xiàn)為反對(duì)稱的運(yùn)動(dòng)軌跡，運(yùn)動(dòng)狀態(tài)及運(yùn)動(dòng)轉(zhuǎn)遷過(guò)程更加復(fù)雜。

4 分?jǐn)?shù)階項(xiàng)對(duì)超諧共振及周期運(yùn)動(dòng)轉(zhuǎn)遷的影響

下面主要研究分?jǐn)?shù)階微分項(xiàng)的系數(shù)λ和階數(shù)p對(duì)系統(tǒng)超諧波共振和周期運(yùn)動(dòng)轉(zhuǎn)遷的影響，選取n=5采用高階諧波平衡法求解系統(tǒng)響應(yīng)，數(shù)值仿真系統(tǒng)動(dòng)力學(xué)響應(yīng)。

首先分析分?jǐn)?shù)階系數(shù)λ對(duì)超諧波共振的影響，選取基準(zhǔn)參數(shù)為：μk=0.15，ε=0.6，f=5，ξ=0.01，取分?jǐn)?shù)階的階數(shù)p=0.5，分?jǐn)?shù)階系數(shù)分別取λ=0.1、λ=0.2、λ=0.3，系統(tǒng)在低頻區(qū)的超諧波共振如圖8(a)所示。由于分?jǐn)?shù)階項(xiàng)不僅具有阻尼作用也具有剛度作用，進(jìn)而隨著分?jǐn)?shù)階系數(shù)λ增大，超諧波共振的幅值與彎曲程度減小。當(dāng)分?jǐn)?shù)階系數(shù)繼續(xù)增大到λ=0.3時(shí)，9次超諧波共振消失，可見(jiàn)較大的分?jǐn)?shù)階系數(shù)λ不僅可以有效降低超諧波共振的幅值，還可以抑制高次超諧波共振的出現(xiàn)。

(a) 不同分?jǐn)?shù)階的系數(shù)λ

(b) 不同分?jǐn)?shù)階的階數(shù)p圖8 超諧波共振幅頻曲線圖Fig.8 Super harmonic resonance amplitude frequency characteristic curve

其次分析分?jǐn)?shù)階的不同階數(shù)p對(duì)系統(tǒng)超諧波共振的影響。選取基本參數(shù)為μk=0.15，ε=0.6，f=5，ξ=0.01，分?jǐn)?shù)階系數(shù)取λ=0.1，分?jǐn)?shù)階的階數(shù)分別取p=0.2、p=0.5、p=0.8，系統(tǒng)在低頻區(qū)的超諧波共振如圖8(b)所示。從圖8(b)可見(jiàn)在頻率ω<0.3的局部放大圖中，分?jǐn)?shù)階的階數(shù)p越小，7、9次超諧波共振的峰值越大；當(dāng)頻率ω>0.3時(shí)，3、5次超諧波共振的峰值隨階數(shù)p的增加而增高。

接著研究分?jǐn)?shù)階系數(shù)λ對(duì)系統(tǒng)周期運(yùn)動(dòng)轉(zhuǎn)遷的影響。數(shù)值仿真得到不同分?jǐn)?shù)階系數(shù)λ所對(duì)應(yīng)系統(tǒng)質(zhì)量塊的分岔圖如圖9所示，虛線表示系統(tǒng)發(fā)生鞍結(jié)分岔SNB。從圖9可見(jiàn)分?jǐn)?shù)階系數(shù)λ由0.1增加到0.2時(shí)，分岔圖中的SNB的個(gè)數(shù)由4減小為1并且3次超諧波和5次超諧波轉(zhuǎn)遷過(guò)程中的混沌消失；隨著分?jǐn)?shù)階系數(shù)λ繼續(xù)增加為0.3時(shí)，鞍結(jié)分岔SNB也消失?？梢?jiàn)分?jǐn)?shù)階系數(shù)λ的增大還可有效避免系統(tǒng)出現(xiàn)分岔和混沌。

(a) λ=0.1的頻率正掃分岔圖和頻率反掃分岔圖

(b) λ=0.2的頻率正掃分岔圖和頻率反掃分岔圖

(c) λ=0.3的頻率正掃分岔圖和頻率反掃分岔圖圖9 不同分?jǐn)?shù)階系數(shù)對(duì)應(yīng)質(zhì)量塊的分岔圖Fig.9 Bifurcation diagram of mass block corresponding to different fractional order coefficients

其次研究分?jǐn)?shù)階的不同階數(shù)p對(duì)系統(tǒng)周期運(yùn)動(dòng)轉(zhuǎn)遷規(guī)律的影響如圖10所示，虛線表示系統(tǒng)發(fā)生鞍結(jié)分岔SNB。當(dāng)分?jǐn)?shù)階的階數(shù)p由0.2增加為0.5時(shí)，SNB的個(gè)數(shù)由5減小為4，在頻率ω=0.234 2的SNB消失，在頻率ω=0.31處的混沌消失；隨著分?jǐn)?shù)階的階數(shù)p繼續(xù)增加為0.8時(shí)，SNB的個(gè)數(shù)由4又增加為5，在頻率ω=0.234 2處再次出現(xiàn)SNB，可見(jiàn)在本節(jié)基準(zhǔn)參數(shù)下為有效抑制在超低頻區(qū)出現(xiàn)分岔和混沌，分?jǐn)?shù)階的階數(shù)p應(yīng)取值適中。

(a) p=0.2的頻率正掃分岔圖和頻率反掃分岔圖

(b) p=0.5的頻率正掃分岔圖和頻率反掃分岔圖

對(duì)比幅頻曲線和分岔圖可見(jiàn)在相鄰次數(shù)超諧波共振的轉(zhuǎn)遷過(guò)渡區(qū)存在混沌，并且存在多態(tài)共存現(xiàn)象。為進(jìn)一步確定多態(tài)共存下的運(yùn)動(dòng)狀態(tài)及轉(zhuǎn)遷規(guī)律，在本節(jié)選取基準(zhǔn)參數(shù)下，選取分?jǐn)?shù)階的系數(shù)λ=0.1和分?jǐn)?shù)階的階數(shù)p=0.5，系統(tǒng)多態(tài)共存區(qū)域主要分布在主共振和3次超諧波共振轉(zhuǎn)遷過(guò)渡區(qū)，以及3次超諧波共振和5次超諧波共振轉(zhuǎn)遷過(guò)渡區(qū)，系統(tǒng)運(yùn)動(dòng)狀態(tài)的多樣性如圖11所示。

圖11 多態(tài)共存及其相鄰區(qū)域分岔圖Fig.11 Multi state coexistence and its adjacent region bifurcation diagram

由于系統(tǒng)存在SNB、倍周期分岔(PDB)、叉式分岔(PFB)以及邊界激變(BC)等多種分岔，在這些分岔的誘導(dǎo)下甚至?xí)霈F(xiàn)混沌，并且該過(guò)程不可逆，具體轉(zhuǎn)遷規(guī)律如圖12所示。

圖12 多態(tài)共存及相鄰區(qū)域運(yùn)動(dòng)轉(zhuǎn)遷規(guī)律Fig.12 Multi state coexistence and movement transition law of adjacent regions

在橡膠工業(yè)中，可根據(jù)上述規(guī)律選擇恰當(dāng)?shù)南到y(tǒng)參數(shù)，進(jìn)而避免隔振系統(tǒng)動(dòng)力學(xué)響應(yīng)出現(xiàn)非線性跳躍現(xiàn)象以及分岔和混沌，使橡膠隔振系統(tǒng)的動(dòng)力學(xué)性能更佳。

5 結(jié) 論

本文采用分?jǐn)?shù)階非線性Zener模型描述黏彈性隔振系統(tǒng)，采用高階諧波平衡法求解了系統(tǒng)的穩(wěn)態(tài)響應(yīng)并結(jié)合多種方法對(duì)結(jié)果進(jìn)行比較，闡述了低頻區(qū)系統(tǒng)幅頻特性與周期運(yùn)動(dòng)多樣性之間的關(guān)系，研究了分?jǐn)?shù)階項(xiàng)對(duì)中低頻范圍內(nèi)隔振系統(tǒng)的分岔、混沌和多態(tài)共存等復(fù)雜非線性動(dòng)力學(xué)行為的影響，得出以下結(jié)論：

(1) 針對(duì)系統(tǒng)在低頻區(qū)域的超諧波共振，分岔圖中出現(xiàn)鞍結(jié)分岔SNB對(duì)應(yīng)的頻率與幅頻曲線發(fā)生跳躍的頻率一致，鞍結(jié)分岔SNB突變的距離隨諧波次數(shù)增加而減小。

(2) 分?jǐn)?shù)階系數(shù)λ的增大不僅可有效降低超諧波共振的峰值，還可有效避免系統(tǒng)發(fā)生分岔和混沌。

(3) 相鄰次數(shù)超諧波共振轉(zhuǎn)遷過(guò)程中存在多態(tài)共存的現(xiàn)象，轉(zhuǎn)遷過(guò)程中相繼出現(xiàn)PFB、SNB、PDF和BC等多種分岔類型，在分岔誘導(dǎo)下甚至?xí)霈F(xiàn)混沌且該過(guò)程不可逆。

上述研究結(jié)果與方法可為高次諧波為小量的隔振系統(tǒng)在微振動(dòng)下的動(dòng)力學(xué)響應(yīng)預(yù)測(cè)提供一定理論依據(jù)。