單搖擺力性分析

本文格式為Word版，下載可任意編輯——單搖擺力性分析單擺的非線性動力學(xué)分析

張亞兵

（蘭州交通大學(xué)車輛工程專業(yè)，甘肅蘭州，730070）

摘要:研究單擺的運動，從是否有無阻尼和驅(qū)動力方面來分析它們對單擺運動的影響。對于小角度單擺的運動，從單擺的動力學(xué)方程入手，借助李雅普諾夫一次近似理論，推導(dǎo)出單擺的運動穩(wěn)定性狀況。再借助繪圖工具matlab，對小角度和大角度單擺的運動進行仿真，通過改變參數(shù)，如阻尼大小、驅(qū)動力大小等繪出單擺運動的不同相圖，對相圖進行分析比較，從驗證單擺運動的穩(wěn)定性狀況。關(guān)鍵詞:單擺；振動；阻尼；驅(qū)動力

Abstract:Thevibrationofsimplependulumisstudiedbyanalyzingwhetherornotdampanddriveforceitsinfluenceofthesimplependulum.Forsmallanglependulummotion,pendulumdynamicequationfromthestart,withanapproximateLyapunovtheoryofstabilityofmotionisderivedpendulumsituation.Drawingtoolswithhelpfrommatlab,smallangleandwide-anglependulummotionsimulation,bychangingtheparameters,suchasdampingsize,drivesizedrawsimplependulumofdifferentphasediagram,analysisandcomparisonofthephasediagram,fromtheverificationthestabilityofthesituationpendulummovement.

Keywords:simplependulum;vibration;ddriveforce

1引言

單擺是一種理想的物理模型[1]，單擺作簡諧振動（擺角小于5°）時其運動微分方程為線性方程，可以求出其解析解，而當(dāng)單擺做大幅度擺角運動時，其運動微分方程為非線性方程，我們很難用解析的方法探討其運動，這個時候可以用MATLAB軟件對單擺的運動進行數(shù)值求解，并可以模擬不可憐況下單擺的運動。隨著擺角的減小，擺球的運動速率將越來越大，而加速度將單調(diào)下降,至??0時,加速度取微小值。本文從動力學(xué)的角度詳細考察了這一過程中擺球的非線性運動,

.得出了在運動過程中??t,???的關(guān)系。

圖1單擺模型

1

2單擺的線性狀況

2.1線性單擺的無阻尼振動

如圖（1）所示，一根不會伸縮、長度為l的細線，上端固定（或一根剛性輕桿，上端與無摩擦的鉸鏈相連），下端懸掛一質(zhì)量為m的小球就構(gòu)成一單擺。當(dāng)擺角?很小時，sin??0，單擺作簡諧振動，令s?l?，其運動方程為

s?Acos(?0t??0)（1）

式中，?0?gl稱為單擺的固有圓頻率，其周期為（2）g式(1)對時間t求一階導(dǎo)數(shù)，得小球的速度

v??A?0sin(?0t??0)（3）

由式(1)與式(3)得

?s??v???1（4）????AA????0?上式說明單擺的速度—位移曲線圖(尋常稱為相圖[2])為一橢圓[3]。

22T0?2?l圖2圖無阻尼單擺的相平面軌跡

2.2線性單擺的阻尼振動

當(dāng)小球搖擺的速度較小時，小球受到一個與速度方向相反的阻力,f???ds/dt，?為阻力系數(shù)，它與物體的形狀以及周邊媒質(zhì)的性質(zhì)有關(guān)。根據(jù)牛頓其次定律有

d2sds??02s?0（5）2?2?dtdt稱為阻尼因數(shù)。對于一定振動系統(tǒng)，根據(jù)比值（?:?0）小2m于、等于、大于1，則稱單擺處于欠阻尼、臨界阻尼和過阻尼振動狀態(tài)。設(shè)振動系統(tǒng)的周期T0?2s，初始條件s(0)?0.01m,v(0)?0，取?=0.3。

用Matlab軟件中的函數(shù)，\與“ezplot〞可求出單擺阻尼振動方程并描繪振動曲線及其相圖。

式中，???2

（a）振動曲線（b）振動相圖

圖3線性單擺的欠阻尼振動曲線及其相圖

線性單擺在臨界阻尼，過阻尼狀態(tài)下的振動曲線及振動相圖，只需改變相應(yīng)的?取值即可

3有阻尼和有驅(qū)動力單擺的運動分析

有阻尼和有驅(qū)動力單擺的運動方程為

??2???sin??fcos?t(6)在任意大振幅下，方程(6)的解變得十分繁雜，下面利用計算機模擬，分別探討單擺運動隨初值的變化和其混沌運動。3.1初值不同所產(chǎn)生的??t曲線

為簡單計，設(shè)??0.10,f?1,??2/3，當(dāng)t=0時，兩振動初始條件相差微小，

.???(0)?0,?1(0)?0.01有?1(7).???2(0)??0.01,?1(0)?0取0?t?120s，對(6)式在初始值(7)式下利用MATLAB繪圖，其??t變化曲線如圖4所示(其中實線為?1?t，虛線為?2?t)。

...圖4有阻尼和驅(qū)動力的??t圖

3

由圖4可以看出，當(dāng)0?t?25s時，兩條曲線重合，兩個解?1(t)、?2(t)不能分辯；但當(dāng)t?25s時，兩條曲線不再重合，兩個解?1(t)、?2(t)、完全不一樣，這種混沌運動對初始條件的敏感性稱為蝴蝶效應(yīng)。

3.2振幅不同所產(chǎn)生的相圖???

為簡單計，設(shè)方程(6)中，除驅(qū)動參數(shù)f取變值外，其余參數(shù)不變，即??0.25,??2/3。

對(6)式在初始值式?(0)?0,?(0)?0下利用MATLAB作計算模擬繪圖，相圖

...???如圖6所示。

當(dāng)f?1.06時，振蕩周期?等于外加周期力的周期T，??T?2?/??3?，應(yīng)應(yīng)單周期解，其相圖???如圖5示。

.

圖5

當(dāng)f?1.07時，??3T，對應(yīng)三倍周期解，其相圖???如圖6示。

.

圖6

當(dāng)周期強迫力的振幅達到某一臨界值f??1.684時，??2?T，系統(tǒng)運動出現(xiàn)混沌，其相圖???如圖7示。

.4

圖7

4無阻尼和無驅(qū)動力單擺的運動(β=0,f=0)

無阻尼無驅(qū)動力狀況下，單擺系統(tǒng)為保守系統(tǒng)，系統(tǒng)的運動微分方程

??2???sin??fcos?t

??????????變?yōu)椋??sin??0，狀態(tài)方程????2???sin??fcos?

???????變?yōu)椋?/p>

????????（8）

?????sin?為了將高階微分方程組降階，我們首先做如下處理：

d?（I）首先把高階微分方程改寫成一階微分方程組，令y1??，y2?，則

dt???y??d??（9）

???dt?于是原二階微分方程可改寫成如下一階微分方程組

?dy1??dt??y2??