謝秀峰，李俊林，劉 迪

(1.太原科技大學(xué) 應(yīng)用科學(xué)學(xué)院，太原 030024；2.山西大學(xué) 數(shù)學(xué)科學(xué)學(xué)院，太原 030006)

隨著新材料、新結(jié)構(gòu)的應(yīng)用發(fā)展，粘彈性材料[1]被廣泛用于航空航天、造船、汽車、鐵路、建筑、紡織等行業(yè)。粘彈性系統(tǒng)的動力學(xué)行為的研究具有重要的理論意義和工程意義，引起了越來越多學(xué)者的重視。BISHOP et al[2]討論了具有粘彈阻尼船舶模型的響應(yīng)。ADHIKARI et al[3]研究了粘彈阻尼的線性振子的特征值和定性動力學(xué)特征。噪聲廣泛存在于實際工程中，其對系統(tǒng)的影響不容忽視，所以研究隨機激勵下粘彈性系統(tǒng)的相關(guān)問題受到很多學(xué)者的關(guān)注[4-6]。ARIARATNAM[7]應(yīng)用隨機平均方法研究了線性粘彈性系統(tǒng)的隨機穩(wěn)定性。XIE[8]研究了有界噪聲激勵下二維粘彈性系統(tǒng)的矩Lyapunov穩(wěn)定性。

由碰撞、干摩擦等非光滑因素增加了動力系統(tǒng)的復(fù)雜性[9]，光滑系統(tǒng)的很多理論與成果不能直接應(yīng)用到非光滑系統(tǒng)的分析當(dāng)中。DIMENTBERG et al[10]利用Dirac delta函數(shù)和符號函數(shù)對非光滑系統(tǒng)進行光滑化處理，將碰撞系統(tǒng)轉(zhuǎn)化為不含碰撞的動力系統(tǒng)，然后用能量平均法，分析了碰撞振動系統(tǒng)的隨機響應(yīng)問題。FENG et al[11]借助平均Poincare映射研究了隨機激勵下的線性碰撞振動系統(tǒng)。ZHAO et al[12]研究了隨機激勵下粘彈性碰撞系統(tǒng)的穩(wěn)態(tài)響應(yīng)。

綜上所述，對于粘彈性碰撞系統(tǒng)的研究較少，因此有必要研究隨機噪聲作用下粘彈性碰撞系統(tǒng)的響應(yīng)。本文研究了高斯白噪聲激勵下粘彈性碰撞系統(tǒng)的矩穩(wěn)定性，用恢復(fù)系數(shù)來描述粘彈性碰撞前后的能量損失，結(jié)合Zhuravlev變換和隨機平均法給出了系統(tǒng)任意階矩Lyapunov指數(shù)。并討論了恢復(fù)系數(shù)和粘彈性系數(shù)對系統(tǒng)矩穩(wěn)定性的影響。

1 模型分析

考慮受高斯白噪聲激勵的粘彈性碰撞系統(tǒng)可表述為

(1)

式(1)中的粘彈性效應(yīng)可表示為：

(2)

(3)

2 隨機平均法

系統(tǒng)(3)的響應(yīng)可近似表示為:

x=A(t)cosΦ(t) .

(4)

(5)

其中，Φ(t)=ωt+φ(t)，且A(t),Φ(t),φ(t)都是隨機過程。系統(tǒng)(3)可寫為關(guān)于幅值A(chǔ)和相位φ的隨機微分關(guān)系如下:

(6)

(7)

由文獻[13]可知，對于一個隨機系統(tǒng)的狀態(tài)變量X(t)，其矩穩(wěn)定性可以用矩Lyapunov指數(shù)Λ(p)表示

(8)

式中：E[·]表示數(shù)學(xué)期望，‖·‖2表示二階范數(shù)，系統(tǒng)解p階矩漸進穩(wěn)定的充分必要條件是Λ(p)<0，并且Λ'(0)等于系統(tǒng)的最大Lyapunov指數(shù)

(9)

而系統(tǒng)的幾乎必然穩(wěn)定的充分必要條件是λ<0.

為了求出系統(tǒng)的p階矩Lyapunov指數(shù)Λ(p)，作變換P(t)=Ap(t),并解式(6)和式(7)，得P與φ的微分關(guān)系

(10)

(11)

(12)

其中W(t)是標(biāo)準(zhǔn)Wiener過程，并且

推導(dǎo)過程見附錄A.

選取Maxwell型粘彈性核函數(shù)

(13)

經(jīng)過式(A3)和式(A4)的變換，得

(14)

(15)

(16)

(17)

得到系統(tǒng)的p階矩Lyapunov指數(shù)為

(18)

最大Lyapunov指數(shù)為

(19)

3 穩(wěn)定性分析

根據(jù)式(18)的解析式給出相應(yīng)的數(shù)值結(jié)果，可分析恢復(fù)因子、噪聲強度和粘彈性參數(shù)對系統(tǒng)(1)隨機穩(wěn)定性的影響。圖1給出系統(tǒng)響應(yīng)的矩Lyapunov指數(shù)Λ(p)作為階數(shù)p的函數(shù)，隨不同的恢復(fù)因子r和噪聲強度σ變化的函數(shù)曲線，其余系統(tǒng)參數(shù)為：γ1=κ1=1,γ2=κ2=0.5,β=0.5,ω=1.0,ε=0.1.實線表示由式(18)得到的解析解，虛線由Monte Carlo模擬的數(shù)值結(jié)果。該圖表明，Λ(p)是p的非線性函數(shù)，在Λ(p)<0的區(qū)域內(nèi)，系統(tǒng)響應(yīng)的p階矩漸近穩(wěn)定。Λ(p)隨恢復(fù)因子r的增大而增大，即r越大系統(tǒng)的矩穩(wěn)定性越弱。Λ(p)隨噪聲強度σ的增大而增大，即噪聲強度越大系統(tǒng)的矩穩(wěn)定性越弱。

圖1 系統(tǒng)(1)的p階矩Lyapunov指數(shù)變化曲線Fig.1 The pth moment Lyapunov exponent of system (1)

圖2、3分別表示矩Lyapunov指數(shù)Λ(p)隨粘彈性特征參數(shù)γ,κ和階數(shù)p的變化三維圖。對于p>0，隨粘彈性參數(shù)γ增大穩(wěn)定區(qū)域變大，說明粘彈性越強有助于系統(tǒng)的穩(wěn)定；而當(dāng)松弛時間減小，穩(wěn)定區(qū)域變窄，表明長的松弛時間有助于系統(tǒng)的穩(wěn)定。

圖2 系統(tǒng)響應(yīng)的矩Lyapunov指數(shù)與粘彈性參數(shù)γ的變化關(guān)系Fig.2 Moment Lyapunov exponent of system with viscoelastic parameter γ

圖3 系統(tǒng)響應(yīng)的矩Lyapunov指數(shù)與粘彈性參數(shù)κ的變化關(guān)系Fig.3 Moment Lyapunov exponent of system with viscoelastic parameter κ

圖4給出了矩Lyapunov指數(shù)Λ(p)與噪聲強度σ和階數(shù)p之間的函數(shù)關(guān)系。該圖表明，隨著噪聲強度σ的增大，矩指數(shù)增大，因此，噪聲強度增大使得系統(tǒng)響應(yīng)的矩穩(wěn)定性減弱。

圖4 系統(tǒng)響應(yīng)的矩Lyapunov指數(shù)與噪聲強度σ的變化關(guān)系Fig.4 Moment Lyapunov exponent of system with noise intensity of σ

4 結(jié)論

本文研究了受高斯白噪聲激勵下粘彈性碰撞系統(tǒng)的隨機穩(wěn)定性。首先，應(yīng)用Zhuravlev變換將粘彈性碰撞系統(tǒng)轉(zhuǎn)換為非碰撞系統(tǒng)，然后用隨機平均法得到系統(tǒng)的隨機微分方程，求得了矩Lyapunov指數(shù)，最后通過解析結(jié)果與Monte Carlo模擬數(shù)值結(jié)果的分析，得到恢復(fù)因子、噪聲強度、粘彈性參數(shù)對系統(tǒng)穩(wěn)定性的影響。研究表明，恢復(fù)因子越大系統(tǒng)穩(wěn)定性越弱；系統(tǒng)的穩(wěn)定性隨噪聲強度的增大而減弱；隨粘彈性參數(shù)的增大，穩(wěn)定區(qū)域增大，隨松弛時間增大，穩(wěn)定區(qū)域變窄。

附錄A

(A1)

其中

R(τ)=E[ξ(τ)ξ(t+τ)]為高斯白噪聲ξ(t)的相關(guān)函數(shù).

(A2)

應(yīng)用變量代換s=t-τ和改變積分次序，得

(A3)

同理可得

(A4)

：

[1] GURTIN M E,STERNBERG E.On the linear theory of viscoelasticity[J].Arch Ration Mech Anal,1962,11:291-356.

[2] BISHOP R E D,PRICE W G.An investigation into the linear theory of ship response to waves[J].Journal of Sound and Vibration,1979,62:353-363.

[3] ADHIKARI S,PASCUAL B.Eigenvalues of linearviscoelastic systems[J].J Sound Vib,2009,325:1000-1011.

[5] HAN X,WANG M.General decay of energy for a viscoelastic equation with nonlinear damping[J].Franklin Inst,2010,347:806-817.

[6] MESSAOUDI S A.On the control of solutions of a viscoelastic equation[J].Franklin Inst,2007,344:765-776.

[7] ARIARATNAM S T.Stochastic stability of linear viscoelastic systems[J].Probab Eng Mech,1993,8:153-155.

[8] XIE W C.Moment Lyapunov exponents of a two dimensional viscoelastic system under bounded noise excitation[J].ASME J Appl Mech,2002,69:346-357.

[9] 金棟平，胡海巖.碰撞振動與控制[M].北京：科學(xué)出版社，2005.

[10] DIMENTBERG M F,IOURTCHENKO D V.Random vibrations with impacts:a review[J].Nonlinear Dyn,2004,36:229-254.

[11] FENG J Q,XU W,RONG H W,et al.Stochastic responses of Duffing-van der Pol vibro-impact system under additive and multiplicative random excitations[J].Int J Non-linear Mech,2009,44:51-57.

[12] ZHAO X R,XU W,GU X D,et al.Stochastic stationary responses of a viscoelastic system with impacts under additive Gaussian white noise excitation[J].Physica A,2015,431:128-139.

[13] ARNOLD L.A formula connection sample and moment stability of linear stochastic systems[J].SIAM Journal of Applied Mathematics，1984，44：793-802.

[14] 朱位秋.隨機振動[M].北京：科學(xué)出版社，1992.