曹振邦，樂 源

(西南交通大學力學與工程學院，成都 610031)

近年來耦合非線性系統(tǒng)在物理、化學、生物、機電設備系統(tǒng)等許多科學領域中都有廣泛的應用，引起了眾多研究者的關注。在這些系統(tǒng)中最具有代表性的，也是研究最深入的要數(shù)Duffing振蕩器系統(tǒng)［1－4］，而具有耦合項是這類系統(tǒng)區(qū)別于其他系統(tǒng)的本質(zhì)特征，因而眾多的國內(nèi)外學者也對其理論推導、數(shù)值模擬等方面進行了大量的研究。Pastor［5］研究了具有對稱耦合項的van der Pol系統(tǒng)，分析了不動點的結構及系統(tǒng)通向混沌的道路。Stagliano［6］研究了一個具有不可約頻率的加性耦合項的兩自由度Duffing系統(tǒng)，觀察到在一定的參數(shù)條件下系統(tǒng)會出現(xiàn)倍周期分岔。進而，Koziowski［7］分析了這類系統(tǒng)的全局分岔現(xiàn)象。目前關于耦合項對系統(tǒng)響應及混沌運動的研究還不夠充分，因此有必要對其做進一步研究。

1 系統(tǒng)模型描述

船舶在海洋中的航行HE航空航天設備的運行［8］可用如下數(shù)學模型來描述:

式中:圓點表示對時間的導數(shù);μ1和μ2為阻尼系數(shù);ω1和ω2為線性固有頻率;k1和k2為非線性剛度系數(shù);α1，α2，β1和 β2表征了兩個單自由度系統(tǒng)的耦合程度;Fcosωt為系統(tǒng)所受的周期性外力。

2 系統(tǒng)穩(wěn)定性分析

無論是對穩(wěn)定性分析或?qū)?shù)值模擬分析來說，把二階微分方程(1)化成一個自治一階微分方程(2)會更方便一些:

由于穩(wěn)定性分析比較簡單，這部分的分析在本文下一節(jié)一并給出結論。

3 系統(tǒng)的分岔及混沌分析

非線性耦合系統(tǒng)的分岔研究對于考察系統(tǒng)與混沌運動的關系具有重要的意義，而分岔圖是研究系統(tǒng)拓撲結構隨參數(shù)變化的良好工具。因此，本文利用系統(tǒng)中較易調(diào)整的參數(shù)，即外加周期力的強度F和外加周期力的頻率ω作為控制參數(shù)研究耦合系統(tǒng)的分岔過程，從而分析其通向混沌的道路。參數(shù)取值為:μ1=0.2;μ2=0.2;ω1=0.2;ω2=0.4;k1=0.3;k2=0.5;α1=0.4;α2=0.2;β1=0.12;β2=0.08;ω =1.0。迭代初值取(－0.3，0.5，0.6，－0.8，0)，控制參數(shù) F取值范圍F∈(0，1.3)，采用四節(jié)龍格庫塔方法進行數(shù)值模擬，結果見圖1～8。

圖1 F變化過程中x的分岔圖

圖2 F變化過程中y的分岔圖

圖3 F變化過程中的最大Lyapunov指數(shù)圖

圖4 F=0.1時的Poincaré映射圖

圖5 F=0.1時的相圖

圖6 F=0.3時的Poincare映射圖

圖7 F=0.3時的相圖

圖8 F=0.85時的Poincaré映射圖

通過比較圖1和圖2可知:在F變化過程中x的分岔圖和y的分岔圖在結構上是相似的，這在定性分析上是必然的。由于最大Lyapunov指數(shù)中μ大于0說明系統(tǒng)處于混沌狀態(tài)，小于0處于非混沌狀態(tài)，故結合圖3的最大Lyapunov指數(shù)圖就可以看出分岔圖是完全正確的。由圖4和圖5可以看出:在很小時(本文選的點F為0.1)系統(tǒng)做的是周期1運動，當F達到0.21時系統(tǒng)突然達到混沌狀態(tài)，而這種混沌狀態(tài)并不能說是由周期1引發(fā)的，只能說在F=0.21點處時周期1運動和混沌是共存的，故該點處的系統(tǒng)的行為是比較復雜的。而在這里通過圖6和圖7可以簡單說明混沌是一種很復雜的混亂狀態(tài)。當達到0.43時系統(tǒng)出現(xiàn)倒分岔現(xiàn)象，又進入到周期1運動繼而進入混沌狀態(tài)，然后系統(tǒng)由周期6通過周期倍化分岔變?yōu)橹芷?2運動(見圖1和圖8)，然后進入混沌區(qū)，最后回歸周期1運動。

4 結束語

通過對周期性驅(qū)動力作用下兩個自由度的Duffing系統(tǒng)分岔及混沌分析可以發(fā)現(xiàn)，由于耦合項的存在使得系統(tǒng)在隨控制參數(shù)變化的過程中其運動形態(tài)非常復雜，可以看出系統(tǒng)經(jīng)由周期倍化分岔可以進入混沌狀態(tài)，而系統(tǒng)由周期1直接突發(fā)進入混沌狀態(tài)這更體現(xiàn)出其復雜性。

［1］冷永剛，賴志慧，范勝波，等.二維Duffing振子的大參數(shù)隨機共振及微弱信號檢測研究［J］.物理學報，2012(23):71－80.

［2］劉海波，吳德偉，金偉，等.Duffing振子微弱信號檢測方法研究［J］.物理學報，2013(5):42－47.

［3］劉海波，吳德偉，戴傳金，等.基于Duffing振子的弱正弦信號檢測方法研究［J］.電子學報，2013(1):8－12.

［4］楊紅英，葉昊，王桂增，等.Duffing振子的Lyapunov指數(shù)與Floquet指數(shù)研究［J］.儀器儀表學報，2008(5):927－932.

［5］Pastor L，Perez-Garcia V M，Encinas-Sanz F，et al.Orderd and chaotic behavior of two coupled Vander Por Oscillators［J］.phys Rev E，1933，48(1):171－182.

［6］Stagliano J J，Wersinger J M，Slamminka E E.Doubling bifurcation of destroyed tori［J］.Physica D，1996，92(3):164－177.

［7］Koziowski J，Pralitz U，Lauterborn W.Bifurcation analysis of two couple periodically driven Duffing oscillator［J］.Phys Rev E，1995，51(3):1861－1867.

［8］陳予恕.非線性振動系統(tǒng)的分岔理論和混沌理論［M］.北京:高等教育出版社，1993:180－198.