輻射噪聲線譜多渦卷混沌同步化方法研究

劉樹勇, 位秀雷, 方　煊, 王　基

(海軍工程大學 動力工程學院，武漢　430033)

?

輻射噪聲線譜多渦卷混沌同步化方法研究

劉樹勇, 位秀雷, 方煊, 王基

(海軍工程大學 動力工程學院，武漢430033)

摘要：針對非線性隔振系統(tǒng)不容易產(chǎn)生小幅值混沌振動的難題，選擇合適的多渦卷混沌系統(tǒng)及其參數(shù)，實現(xiàn)了非線性隔振Duffing系統(tǒng)與渦卷混沌系統(tǒng)同步，重構(gòu)了隔振系統(tǒng)響應(yīng)信號的頻譜，得到小幅值的寬頻動力學行為。研究了不同參數(shù)條件下渦卷混沌行為以及對應(yīng)的頻譜特征，通過計算系統(tǒng)響應(yīng)的最大Lyapunov指數(shù)，證實該系統(tǒng)處于混沌狀態(tài)。通過參數(shù)驅(qū)動同步策略，發(fā)現(xiàn)非線性隔振系統(tǒng)在周期狀態(tài)下，能夠被有效地同步到混沌狀態(tài)。而且從系統(tǒng)的分岔圖以及吸引子相圖可以看出，同步后系統(tǒng)的混沌響應(yīng)具有明顯的小幅值特性；同時驅(qū)動系統(tǒng)的渦卷參數(shù)可調(diào)，可以滿足不同系統(tǒng)的同步需求。

關(guān)鍵詞：多渦卷；混沌；同步；Duffing系統(tǒng)；小幅值

混沌是確定性系統(tǒng)產(chǎn)生的貌似隨機的運動,混沌振動的一些優(yōu)良特性可用于船舶動力機械的振動噪聲控制[1]。一方面由于混沌能夠改變其輸出響應(yīng)的頻譜結(jié)構(gòu)，從而改變潛艇輻射水聲的頻譜特征；另一方面，混沌使單頻輸入產(chǎn)生寬頻輸出，將集中于某一頻率的能量分散到較寬的頻帶上，從而降低輻射噪聲中的線譜成分[2-3]。

目前，基于混沌理論的線譜控制方法，在理論分析[4]、混沌區(qū)域預測[5]、隔振性能評估[6]以及混沌識別[7]等方面，取得了諸多成果。但實現(xiàn)其工程應(yīng)用等方面還存在若干難點，其中，如何保持持續(xù)的小幅值混沌振動是線譜混沌化在工程實際應(yīng)用的一個關(guān)鍵問題[8]。劉樹勇等[9]提出一種機械式混沌同步系統(tǒng)，在主系統(tǒng)為大振幅混沌狀態(tài)時，實現(xiàn)了子系統(tǒng)的小幅值混沌態(tài)。楊慶超等[10]設(shè)計了投影同步的試驗臺架，實現(xiàn)了隔振系統(tǒng)在變工況和小振幅條件下的混沌響應(yīng)。上述方法大多集中在混沌的單渦卷和雙渦卷形式，振動的能量還不能夠充分地分散到其他頻帶上，阻礙了混沌化目標的實現(xiàn)。

近年來，多渦卷混沌吸引子的研究逐漸成為混沌研究領(lǐng)域的一個新方向[11]，相比于傳統(tǒng)的單渦卷和雙渦卷混沌系統(tǒng)，多渦卷混沌系統(tǒng)呈現(xiàn)出更為復雜的結(jié)構(gòu)和動力學行為[12-13]。本文研究了不同參數(shù)條件下多渦卷混沌行為以及對應(yīng)的頻譜特征，采用多渦卷混沌同步方法將系統(tǒng)由周期運動狀態(tài)誘導進入混沌狀態(tài)，通過參數(shù)驅(qū)動同步策略，實現(xiàn)了非線性隔振Duffing系統(tǒng)的小幅值混沌行為。

1非線性隔振系統(tǒng)無量綱模型

單層非線性隔振系統(tǒng)模型如圖1所示，得系統(tǒng)的運動微分方程如下：

(1)

式中：M為質(zhì)量；C為阻尼器的黏滯阻尼系數(shù)；K1X+K3X3為非線性彈性恢復力；F0為作用在質(zhì)塊M上的激勵力。其中，K3為負時，系統(tǒng)呈現(xiàn)軟特性；K3為正時，系統(tǒng)呈現(xiàn)硬特性。

圖1　Duffing隔振系統(tǒng)Fig.1 Duffing vibration isolation system

(2)

則有

(3)

將式(3)代入式(1)，得

(4)

2混沌同步原理

考慮如下參數(shù)驅(qū)動單向耦合系統(tǒng)[8]：

(5)

式中：x={x1,x2,…,xn}為n維空間Rn中的驅(qū)動狀態(tài)向量，y={y1,y2,…,ym}為m維空間Rm中的響應(yīng)狀態(tài)向量。F和G都為向量場，F(xiàn)：Rn→Rn，G：Rm→Rm。P表示系統(tǒng)之間的耦合狀態(tài)向量，P：Rn→Rm，kt為表示耦合強度系數(shù)。

要判斷式(5)是否存在廣義同步，通常采用的方法是構(gòu)造一個輔助系統(tǒng)y′)，它與式(5)第二式的形式一致：

(6)

式(5)第一式和第二式實現(xiàn)廣義同步的充要條件是式(5)第二式和式(6)實現(xiàn)精確同步，由于式(5)第二式和式(6)實現(xiàn)精確同步與式(5)第二式漸近穩(wěn)定等價，因此式(5)第一式與第二式實現(xiàn)廣義同步當且僅當該系統(tǒng)第二式漸近穩(wěn)定。

3非線性隔振系統(tǒng)的多渦卷混沌同步

輻射噪聲線譜的混沌化，在本質(zhì)上是要使得系統(tǒng)將振動能量從線譜轉(zhuǎn)移到其他頻率成分，從而削減特征頻譜的能量。傳統(tǒng)的混沌系統(tǒng)雖然具有寬頻響應(yīng)，但大多集中在單渦卷和雙渦卷形式，振動的能量還不能夠充分地分散到其他頻帶上，阻礙了混沌化目標的實現(xiàn)。構(gòu)造多渦卷的混沌系統(tǒng)并與非線性隔振系統(tǒng)混沌同步，將為解決該問題提供新的思路。

選取具有鞍焦點的系統(tǒng)：

(7)

對應(yīng)的Jacobi矩陣為：

將式(5)具體化為：

(8)

根據(jù)系統(tǒng)的Jacobi矩陣，當參數(shù)為u=-8,v=0.4,w=-6時，系統(tǒng)的特征值為λ1，2=0.2±2.821i,λ3=-6，這表明系統(tǒng)的行為將限制在某一區(qū)域，并具有類似蝴蝶吸引子的纏繞軌道。為了重構(gòu)多渦卷混沌系統(tǒng)，取非線性函數(shù)為：

sgn[s-(2j-1)]}

(9)

聯(lián)合式(8)和式(9)，可進一步得到具有非線性項的動力學系統(tǒng)：

(10)

當系統(tǒng)的參數(shù)在特定范圍時，系統(tǒng)將呈現(xiàn)混沌行為。將式(10)作為驅(qū)動系統(tǒng)，式(4)作為被驅(qū)動系統(tǒng)，在式(4)中引入剛度系數(shù)后寫成狀態(tài)方程形式：

(11)

同步過程中采用參數(shù)驅(qū)動模式，將驅(qū)動信號p2作用于式(11)的線性剛度k1，因此得到系統(tǒng)的k1=kp2，其中，k為驅(qū)動強度，因此，得到：

(12)

4仿真研究

4.1驅(qū)動系統(tǒng)的動力學特性

對式(10)而言，取u=-8，v=0.4，w=-6，Q=1時，系統(tǒng)的響應(yīng)如圖2所示，此時系統(tǒng)的吸引子是由三個螺旋上升的渦卷構(gòu)成，頻率譜呈現(xiàn)連續(xù)譜狀態(tài)，最大Lyapunov指數(shù)為0. 221。當Q的值進行調(diào)整時，系統(tǒng)渦卷數(shù)不斷改變，Q為2時，生成5個渦卷，Q=3時生成7個渦卷(見圖3,圖4)。由此可見，調(diào)整該值能夠得到不同的工作狀態(tài)，為隔振系統(tǒng)的混沌同步奠定了良好的基礎(chǔ)。

圖2　u=-8;v=0.4;w=-6，Q=1時，系統(tǒng)的響應(yīng)Fig.2 System response when u=-8;v=0.4;w=-6，Q=1

圖3　u=-8;v=0.4;w=-6，Q=2時，系統(tǒng)的響應(yīng)Fig.3 System response when u=-8;v=0.4;w=-6，Q=2

圖4　u=-8;v=0.4;w=-6，Q=3時，系統(tǒng)的響應(yīng)Fig.4 System response when u=-8;v=0.4;w=-6，Q=3

4.2被驅(qū)動系統(tǒng)的動力學特性及同步響應(yīng)

(1) 動力學特性分析

以線性剛度為負的硬彈簧Duffing系統(tǒng)為研究對象，對式(11)取不同的參數(shù)進行分析。首先采用全局分岔方法來分析系統(tǒng)的整體動力學性態(tài)。當系統(tǒng)的參數(shù)為k1=-1，k3=1，ξ=0.3，ω=1.2，以激勵力幅值f為變化參數(shù)，得到系統(tǒng)的分岔圖如圖5所示。從圖中可以得知系統(tǒng)在初始階段出現(xiàn)周期1的解，隨后出現(xiàn)多周期解，經(jīng)過倒分岔進入周期1解，再次經(jīng)過倍周期分岔后，進入混沌狀態(tài)。但是激勵力的幅值和系統(tǒng)響應(yīng)的幅值已經(jīng)明顯增大，更確切地說，這種大幅值混沌振動狀態(tài)和隔振系統(tǒng)小幅值混沌狀態(tài)要求相矛盾，必須采取新的方式誘導系統(tǒng)進入有益的混沌行為中。

圖5　系統(tǒng)的分岔圖Fig.5 Bifurcation diagram of the system

根據(jù)全局分岔圖，當系統(tǒng)的其他參數(shù)不變，激勵力幅值大小f=11.5時，系統(tǒng)的胞映射圖如圖6所示，從中可以看出，系統(tǒng)的對初值具有很強的敏感性，混沌吸引域相互纏繞，具有明顯的分形結(jié)構(gòu)。系統(tǒng)相平面的吸引子如圖7所示，功率譜圖如圖8所示。吸引子不同于不動點、極限環(huán)和環(huán)面等平庸吸引子，而是具有自相似結(jié)構(gòu)的奇怪吸引子，功率譜呈現(xiàn)譜特征。系統(tǒng)的行為具有混沌特性，但混沌響應(yīng)處于較大振幅狀態(tài)。

從系統(tǒng)的分岔規(guī)律可知，當參數(shù)選擇為f=2.3時，系統(tǒng)出現(xiàn)周期1解，通過PMUCR胞映射法得到系統(tǒng)的胞映射圖如圖9所示，在該圖中可以看出，當系統(tǒng)的初始條件取不同的值時，都將呈現(xiàn)周期1響應(yīng)。圖10是系統(tǒng)在該范德玻爾平面上的解，所有的曲線都收斂在一個吸引子上，這驗證了系統(tǒng)的行為具有唯一性。

圖6　系統(tǒng)的胞映射圖Fig.6 Cellular mapping of the system

圖7 相平面吸引子圖Fig.7Phaseplaneattractor圖8 響應(yīng)功率譜圖Fig.8Responsepowerspectrum

圖9 周期1解的胞映射圖Fig.9Cellmappingofperiodic1solutions圖10 在范德玻爾平面上的解Fig.10SolutionintheplaneofvanderBohr

對系統(tǒng)求解，到了系統(tǒng)的相平面圖和Poincare截面圖(見圖11)。系統(tǒng)的相圖是一個橢圓，Poincare截面圖只有一個點，因此當沒有使用驅(qū)動系統(tǒng)時，此時系統(tǒng)不會因為初始條件的不同而產(chǎn)生混沌運動，其響應(yīng)的位移功率譜如圖12所示。

圖11 系統(tǒng)響應(yīng)的相平面圖Fig.11Phaseplanofsystemresponse圖12 系統(tǒng)響應(yīng)的功率譜圖Fig.12Powerspectrumofsystemresponse

(2) 混沌同步結(jié)果

當系統(tǒng)的參數(shù)為f=2.3，k1=-1，k3=1，ξ=0.3，ω=1.2時，系統(tǒng)處于周期1狀態(tài)。采用3渦卷同步時，驅(qū)動強度k=1，能夠?qū)⒃撝芷?狀態(tài)驅(qū)動進入混沌狀態(tài)，而最大幅值只有2.54，如圖13所示，響應(yīng)的最大Lyapunov指數(shù)為0.118，顯然是一種獲得小幅值混沌的有效手段。而不采用同步方式時，系統(tǒng)的參數(shù)為k1=-1，k3=1，ξ=0.3，ω=1.2，但激勵的幅值需要達到f=11.5該系統(tǒng)才能夠出現(xiàn)混沌，而響應(yīng)的幅值達到7.78，兩者的比較如圖14所示。

圖13 同步混沌響應(yīng)相圖Fig.13Phasediagramofchaoticsynchronizedresponse圖14 非同步混沌響應(yīng)相圖Fig.14Phasediagramofchaoticnon-synchronizedresponse

5結(jié)論

通過對非線性隔振Duffing系統(tǒng)的全局分岔分析，得到了系統(tǒng)隨激勵力幅值變化的整體性態(tài)圖，觀察到了系統(tǒng)在在不同參數(shù)時的周期1、周期2解以及由倍周期分岔導致的混沌行為。應(yīng)用PMUCR胞映射方法得到了系統(tǒng)在不同參數(shù)區(qū)域時初始條件對系統(tǒng)解的影響，在混沌參數(shù)范圍時，胞映射圖具有明顯的分形結(jié)構(gòu)，反映了混沌系統(tǒng)內(nèi)部復雜但有規(guī)律的動力學特征。在周期1解情況下，初始條件對系統(tǒng)的影響不明顯，應(yīng)用范德玻爾平面證實了解的唯一性。由于系統(tǒng)本身產(chǎn)生混沌時，處于大幅度振動狀態(tài)，不利于輻射噪聲的控制，因此采用多渦卷混沌同步方法將系統(tǒng)由周期運動狀態(tài)誘導進入混沌狀態(tài)，實現(xiàn)了非線性隔振Duffing系統(tǒng)的小幅值混沌行為，從而為噪聲線譜混沌控制提供有益參考。

參 考 文 獻

[1] 劉樹勇，朱石堅，俞翔. 準周期激勵非線性隔振系統(tǒng)的混沌研究[J]. 船舶力學，2010，14(1)：141-147.

LIU Shu-yong,ZHU Shi-jian,YU Xiang. Study on the chaos of nonlinear vibration isolation system under quasi-periodic excitation[J]. Journal of Ship Mechanics,2010,14(1):141-147.

[2] Lou J J，Zhu S J，He L，et al. Application of chaos method to line spectra reduction[J].Journal of Sound and Vibration, 2005，286(3):645-652.

[3] Zhu S J，Zheng Y F，F(xiàn)u Y M. Analysis of non-linear dynamics of a two-degree-of-freedom vibration system with non-linear damping and non-linear spring[J].Journal of Sound and Vibration,2004，271(1/2):15-24.

[4] Yu X，Zhu S J，Liu S Y．Bifurcation and chaos in multi-degree-of-freedom nonlinear vibration isolation system[J]. Chaos，Solitons，and Fractals, 2008，38 ( 5 ): 1498-1504.

[5] 張振海，朱石堅，樓京俊. 基于跟蹤混沌化方法的線譜控制技術(shù)研究[J]. 振動與沖擊，2011，30(7):40-44.

ZHANG Zhen-hai,ZHU Shi-jian,LOU Jing-jun. Line spectra reduction of a vibration isolation system based on tracking chaotification method[J]. Journal of Vibration and Shock, 2011，30(7):40-44.

[6] 浣石，陶為俊，朱石堅，等. 硬特性非線性隔振裝置混沌動力學特性研究[J]. 振動與沖擊，2011，30(11):245-248.

HUAN Shi, TAO Wei-jun, ZHU Shi-jian, et al. Chaotic dynamics of hardening nonlinear isolation device[J]. Journal of Vibration and Shock, 2011，30(11):245-248.

[7] 劉樹勇，楊慶超，位秀雷，等. 鄰近點快速搜索算法在混沌識別中的應(yīng)用[J]. 華中科技大學學報:自然科學版，2012，40(11):89-92.

LIU Shu-yong, YANG Qing-chao, WEI Xiu-lei, et al.The application of fast searching nearest points method to chaos identification [J]. Journal of Huazhong University of Science and Technology:Natural Science Edition, 2012，40(11):89-92.

[8] Wen G L，Lu Y Z，Zhang Z Y，et al．Line spectra reduction and vibration isolation via modified projective synchronization for acoustic stealth of submarines[J]. Journal of Sound and Vibration， 2009，324: 954-961.

[9] 劉樹勇，楊愛波，楊慶超，等. 機械式混沌同步系統(tǒng)動力學特性研究[J]. 華中科技大學學報:自然科學版，2012，40(8):108-111.

LIU Shu-yong, YANG Ai-bo, YANG Qing-chao,et al. Dynamic characteristics of mechanical chaos synchronization systems [J]. Journal of Huazhong University of Science and Technology:Natural Science Edition, 2012，40(8):108-111.

[10] 楊慶超，朱石堅，劉樹勇，等. 投影混沌同步試驗研究[J]. 振動與沖擊，2012，31(20):126-130.

YANG Qing-chao,ZHU Shi-jian,LIU Shu-yong, et al. Experimental study on projective synchronization to a vibration isolation system[J]. Journal of Vibration and Shock,2012，31(20):126-130.

[11] 諶龍，彭海軍，王德石. 一類多渦卷混沌系統(tǒng)構(gòu)造方法研究[J]. 物理學報，2008，57(6):3337-3341.

CHEN Long, PENG Hai-Jun, WANG De-Shi. Study on the construction of a family of multi-scroll chaotic systems[J]. Acta Physica Sinica, 2008，57(6):3337-3341.

[12] Zhang C X，Yu S M. Generation of grid multi-scroll chaotic attractors via switching piecewise linear controller[J]. Physics Letters A，2010，374(30):3029-3037.

[13] 張朝霞，禹思敏. 用時滯和階躍序列組合生成網(wǎng)格多渦卷蔡式混沌吸引子[J]. 物理學報，2009，58(1):120-130.

ZHANG Chao-xia, YU Si-min. Generation of grid multi-scroll Chua’s chaotic attractors with combination of hysteresis and step series[J]. Acta Physica Sinica, 2009，58(1):120-130.

基金項目：國家自然科學基金(51179197);國家自然科學基金面上項目(51579242);海洋工程國家重點實驗室(上海交通大學)開放課題(1009);國家留學基金委資助

收稿日期：2015-02-11修改稿收到日期：2015-08-10

通信作者位秀雷 男，博士生，1988年10月生

中圖分類號:O322

文獻標志碼：A

DOI:10.13465/j.cnki.jvs.2016.14.022

Multi-scroll chaotic synchronization method for radiated noise line spectrum

LIU Shu-yong, WEI Xiu-lei, FANG Xuan, WANG Ji

(College of Power Engineering, Naval University of Engineering, Wuhan 430033, China)

Abstract:The nonlinear vibration isolation system is hard to produce chaotic vibration with small amplitude. In order to solve the problem, an appropriate multi-scroll chaotic system and its parameters were selected and the synchronization between the multi-scroll chaotic system and Duffing system were achieved. The frequency and spectrum of the isolation system was reconstructed, and then the dynamics behavior with broadband spectrum and small amplitude was realized. The chaotic dynamic characteristics under different parameters of the multi-scroll system and the corresponding frequency spectrum characteristics were studied. The largest Lyapunov exponent of the system response was calculated to confirm the chaotic state of system. By using the parametric driven synchronization strategy, it is found that the synchronization between the periodic state and chaotic state of the nonlinear isolation system can be achieved exactly, and the chaotic response after synchronization has obvious characteristics of small amplitude which can be presented in the bifurcation and attractor diagram. In addition, the scroll number of drive system is adjustable, so the synchronizations of different systems can be implemented.

Key words:multi-scroll; chaos; synchronization; Duffing system; small amplitude

第一作者 劉樹勇 男，副教授，1975年11月生