非線性能量阱系統的強調制響應研究

柴凱 李爽 樓京俊 朱石堅

摘要：為了得到諧波激勵下本質非線性能量阱系統出現強調制響應的充要條件，開展非線性能量阱系統的強調制響應研究.利用復變量平均法推導主共振下系統的慢變動力流方程；結合多尺度法實現系統快、慢變流形的分離，得到不同參數下系統的慢不變流形以及全局分岔特性；通過構建快變流形的一維映射函數，揭示非線性能量阱系統呈現強調制響應狀態的充要條件；搭建非線性能量阱系統的模擬電路，進行強調制響應檢測電路試驗.仿真和試驗結果表明，強調制響應是由耦合系統慢變動力流中極限環的鞍結分岔引起的且真實存在，系統出現穩定的強調制響應必須滿足兩個條件：一是非線性能量阱系統響應超越慢不變流形上的極值點，又不吸引至慢不變流形某一分支；二是形成連續跳躍環路而不陷入局部循環時.

關鍵詞：非線性能量阱；強調制響應；非線性動力系統；復變量平均法；振動控制

中圖分類號：O322 文獻標志碼：A

Investigation on Strongly Modulated Response of Nonlinear Energy Sink System

CHAI Kai，LI Shuang，LOU Jingjun，ZHU Shijian

（College of Naval Architecture and Ocean，Naval University of Engineering，Wuhan 430033，China）

Abstract：In order to obtain the necessary and sufficient conditions for the strongly modulated response （SMR）of a nonlinear energy sink （NES）system under harmonic excitation，the SMR study of the NES system is carried out. Firstly，the equation of systematic slow variation flow is derived by using the complex-averaging method. Secondly，the multi-scale method is used to realize the separation of the rapid variable manifold and the slow variable manifold. And the systematic slow invariant manifold under the different parameters and the global bifurcation property are obtained. Then，by building the one-dimensional mapping function of the rapid variable manifold，the sufficient and necessary conditions for the NES system in the SMR state are revealed. Finally，the simulation circuit of the NES system is constructed，and the response detection circuit test is carried out. The simulation and experimental results show that the SMR is caused by the saddle-node bifurcation of limit cycles in the slowly varying power flow of the coupling system，and it is an actual phenomenon in engineering. The system in which SMR state can appear must satisfy the following two conditions：the response of the NES system exceeds the amplitude of the extremum point on the slow invariant manifold，but it does not attract to a branch of the slow invariant manifold，and a continuous jump loop without falling into a local cycle is formed.

Key words：nonlinear energy sink（NES）；strongly modulated response（SMR）；nonlinear dynamical systems；complex-averaging method；vibration control

非線性能量阱（Nonlinear Energy Sink，NES）是指剛度近似為立方剛度用于振動控制的非線性振子，它在一定條件下會出現靶能量傳遞現象.與傳統線性吸振器相比，NES能有效增加吸振帶寬，大幅提升減振效率，目前在航空航天[1]、房屋橋梁抗震[2]、能量采集[3]以及振動噪聲控制[4]等領域得到了廣泛的應用.在受外加激勵作用的耦合NES系統中，由于非線性參數的作用，系統會呈現幾種具有明顯差異的響應形式，即周期、弱調制、強調制和混沌響應.廣義上前3種類型均屬于穩態響應，但其對應相軌跡具有明顯差異，使得NES的振動抑制效果也具有很大的不同[5].

NES因其寬頻控制的特性引起了相關學者的廣泛關注.Jiang等[6]研究了正弦激勵作用下，NES在較寬頻帶從線性振子吸收能量，且不論向前還是向后掃頻，均能實現能量定向傳遞.Zhang等[7]研究了主結構受正弦周期力、多約束下碰撞非線性系統的能量傳遞問題，指出當激勵幅值達到一定閾值時，能量會向NES振子聚積，而且根據激勵幅值不同，系統會呈現周期、弱調制、強調制等多種具有明顯差異的響應形式.Starosvetsky等[8]通過相軌跡法給出了系統響應類型與平衡點之間的關系，指出系統慢變方程對應周期解是否分岔是NES振子實現靶能量傳遞的關鍵.Gourdon等[9]研究了NES系統的強調制響應（Strongly Modulated Response，SMR），給出了系統產生強調制響應的條件，通過對比說明了NES的減振效果在非周期響應時要比穩態周期響應更好.李爽等[10]采用柔性鉸鏈結構提出了一種NES構造方法，分析了簡諧激勵下耦合系統的局部分岔特性.張也弛等[11]通過數值方法研究了兩自由度NES系統在簡諧激勵下的力學特性與抑振效果.由于強調制響應下NES系統的振動抑制效率比穩態振動更優越，因此有必要對系統產生強調制響應的充要條件作進一步研究.

本文從近似解析計算角度出發，重點研究NES系統的強調制響應.利用復變量平均法推導主共振下系統響應幅值的慢變動力流方程，深入探討NES 系統產生強調制響應的充要條件，通過電路仿真和試驗驗證理論分析的正確性和強調制響應的真實性，從而間接證明NES振子在簡諧激勵能否實現靶能量傳遞.

1NES系統的慢不變流形

1.1動力學建模

如圖1所示，建立機械設備耦合非線性能量阱的兩自由度非線性吸振系統動力學模型（以下簡稱NES系統）.機械設備主系統中m₁為待減振機械設備，通過線性剛度彈簧k₁、阻尼λ₁與剛性基座連接；NES子系統安裝于機械設備上層并與其耦合連接，其組成元素包括質量m₂、剛度k₂以及阻尼λ₂；f_b=F cos ΩT為作用于機械設備上的外界激勵信號，F為激勵力幅值，Ω為激勵力頻率，機械設備和NES產生的垂向位移分別為z₁和z₂.本文主要研究具有本質非線性的立方剛度型NES，故其非線性彈簧回復力為f_NES=k₂（z₂-z₁）³.

相應的動力學方程為：

采用復變量平均法推導系統的慢變動力流方程.以質量比ε為小參量，引入長度量綱l₀=m₁g/k₁，l₀對應線性彈簧k₁在重量G₁=m₁g作用下的靜態形變量，并作如下無量綱變換：

將式（2）代入式（1）可得：

考慮激勵頻率接近主系統固有頻率時的1：1：1主共振響應，令ω=1+εσ，其中σ為調諧參數，用來描述內共振頻率之間的接近程度.將系統振動響應分解為質心運動以及兩個振子之間的相對運動，對系統再次簡化，引入新變量：u=x₁+εx₂，v=x₁- x₂，在不考慮主系統阻尼（即ξ₁=0）的情形下，將u和v代入式（3）可得：

1.2系統慢不變流形

引入新的時間尺度τ_k=ε_kt，k=0，1，…，并令ρ₂=ρ₂（τ₀，τ₁，…），采用多尺度法展開，并忽略高階項，可得式（6）關于ε的前兩階方程：

考慮慢變幅值ρ₂關于快變時間尺度τ₀的近似解，對式（7）第一式進行積分，可得：

式中：R為關于τ₁尺度的任意函數.令?ρ₂/?τ₀=0，可得系統平衡點方程：

由式（9）可知，平衡點Φ（τ₁）只與τ₁有關，將復變量用模和相角表示，即Φ（τ₁）=N（τ₁）exp[jθ（τ₁）]，（N，θ）∈（R⁺×S¹），代入式（9）求模后可得：

Z（τ₁）=N²（τ₁），對應相位如下：

由圖2可知，由于流形的不變性，系統慢變幅值只能沿著曲線運動，當達到折疊線N₁、N₂時，可能從一支穩定分支跳躍至另一支穩定分支，從而出現跳躍現象.此外，在折疊線N₁、N₂處可能存在鞍結分岔過程，從而形成N₁→N_u→N₂→N_d→N₁的連續跳躍環路并使系統響應慢變幅值出現周期性變化，類似于目標能量轉移過程中“泵能”與“放能”現象.文獻[13]將其稱為強調制響應，該響應下NES系統出現靶能量傳遞現象，在其主共振附近會出現準周期響應，且幅值變化較為劇烈，能使NES具有優異的振動抑制效果.

2強調制響應產生的必要條件

2.1折奇點產生條件和慢不變流形相軌跡特征

式中：Φ（τ₁）=N（τ₁）exp[jθ（τ₁）].對式（12）兩邊取復數共軛，通過化簡得到：

1）當g（N）≠0時為普通平衡點，滿足如下表達式：

2）當g（N）=0時為折奇點，滿足如下表達式：

式中：

選取系統參數C=4/3，ξ₂=0.2，σ=0，普通平衡點分岔情況如圖3所示.由圖3可知，當f變化時，系統只存在一個周期解，不穩定解位于折疊線N₁= 0.595 0、N₂= 0.989 6之間.另外，通過聯合求解式（14）可得折疊線對應的激勵力分岔值為f_b1= 0.241 9、f_b2= 0.989 2.

進一步求解f_ck可得折奇點對應分岔值為f_c1=0.1760、f_c2=0.9844，滿足f_c1＜f_b1＜f_c2＜f_b2系統的全局分岔如圖4所示，圖中rp為周期吸引子，fs1、fs2分別為上、下折奇點，bp1、bp2為周期吸引子分岔點，普通平衡點與折奇點穩定性通過慢變系統線性化擾動方程的雅可比矩陣特征值判斷.同時，結合Runge-Kutta數值方法與Matlab Streamline命令繪制不同激勵力幅值下系統的慢不變流形，如圖5所示. 圖5中縱坐標為慢變幅值N，橫坐標為相位θ，

由圖4和圖5可知，當f=0時，從N₂上方出發的相軌跡都可以回到N₂，但從N₁出發的相軌跡不能返回至N₁.與圖2對比可知，相軌跡從上穩定分支出發時能跳躍至下穩定分支，而從下穩定分支出發的相軌跡不能跳躍至上穩定分支，這與實際系統也是相符合的.當不存在外界激勵時，由于阻尼存在，系統能量會逐漸被耗散，直至趨于穩定.當f = 0.1時，對應f ＜f_c1，存在一個穩定的周期平衡點，不存在折奇點，所有的相軌跡都流入該吸引子.當f =f_c1時，普通平衡點出現亞臨界分岔，在折疊線N₁處產生一對不穩定的下折奇點.當f = 0.18時，對應f略大于f_c1，存在一個周期平衡點與一對下折奇點（左側下折奇點為結點，右側下折奇點為鞍點），從鞍點右邊出發的相軌跡都被吸引至周期平衡點，結點左邊以及折疊線N₁右邊部分出發的相軌跡吸引至結點，而結點與鞍點之間的相軌跡則有可能返回至折疊線N₁，這意味著強調制響應可能發生.當f_c1＜f＜f_b1時，存在一對不穩定的下折奇點和一個穩定的周期吸引子.當f = f_b1時，周期吸引子退化為不穩定解，同時下折奇點由不穩定變為穩定.當f = 0.5時，對應f_b1＜f ＜f_c2，系統性態發生了非常顯著的變化：首先，和f = 0.18 相比，普通平衡點消失，同時鞍點與結點沿著下折疊線往兩側移動，其中結點往左側移動，而鞍點往右側移動，鞍結點之間的距離擴大，表明最終能回到N₁的相軌跡區域變大，同時出現強調制響應的可能性也將增加；另外，由于在繪制相軌跡過程中，相位只取一個周期（0，2π），當結點運動超過最左側時，又重新在右端出現，此時折奇點性質已由一對鞍結點通過碰撞演化成了一對穩定的焦點.當f =f_c2時，平衡點再次出現亞臨界分岔，此時演變出一對穩定的上折奇點.當f = 0.987時，對應f略大于f_c2，在上折疊線N₂出現一對鞍結點，而下折疊線N₁的一對穩定焦點依然存在，仍有可能出現強調制響應.當f=f_b2時，上折奇點由穩定退化為不穩定，同時周期吸引子再次發生Hopf分岔，由不穩定變為穩定.當f = 1時，對應f ＞f_b2，折疊線N₂上鞍結點演化成了一對不穩定的鞍點，同時在折疊線N₂附近出現了穩定的普通平衡點，此時強調制響應依然有可能出現.

2.2系統響應慢變流形

仍考慮時間尺度的前兩階，令ρ₁（t）=a₁（t）+jb₁（t）、ρ₂（t）=a₂（t）+jb₂（t），對式（5）利用多尺度法直接展開，取ε⁰階分離實部與虛部，可得四階常微分方程組：

對比圖7中的系統響應時間歷程也可驗證以上分析結論，圖中實線為式（5）數值計算得到的實際響應，虛線為式（16）數值計算得到的慢變幅值.由圖7 可知，通過慢變方程計算得到的幅值響應與真實幅值存在一定的誤差，這主要是由設定質量比為小參量條件引起的，但總體來說，預測值還是可靠的.

圖8給出了C = 4/3、ξ₂= 0.2、f = 0.5、σ=0，而質量比分別為0.05、0.01和0.001時的系統響應慢變流形.質量比并不影響慢不變流形的形狀，由圖8可知，3種情形均能出現強調制響應，且質量比參數越小，吻合程度越高，但總體而言，通過觀察慢變流形變化趨勢就可判斷系統是否出現了強調制響應.

總結而言，f>f_c1并不能說明系統一定會出現強調制響應，當激勵頻率變化時，系統響應有可能被吸引至慢不變流形的某一穩定分支或陷入局部循環，從而導致無法形成N₁→N_u→N₂→N_d→N₁的連續跳躍環路，也就不能產生強調制響應.因此，系統出現折奇點只是出現強調制響應的必要條件，仍需進一步探求系統出現強調制響應的充分條件.

3強調制響應產生的充分條件

3.1一維映射函數構建

仍考慮慢變系統，由式（15）可得N₁折疊線上折奇點對應跳躍的邊界條件為Θ₁和Θ₂，令相軌跡可能從折疊線N₁跳躍至上穩定分支的相位區間為R=[Θ₁，Θ₂].考慮R→R的一維映射，若從區間R出發的相軌跡經多次N₁→N_u→N₂→N_d→N₁4個階段后仍可返回至R，則表明慢變系統出現了穩定的極限環，強調制響應必然會發生.而系統能夠產生穩定極限環的相軌跡區域則對應出現強調制響應的初始條件，這顯然也是產生強調制響應的充分條件.

2）第2階段：N_u→N₂.由于慢不變流形的確定性，從（N_u，θ_u）出發的相軌跡必然只能沿著上穩定分支運動，對式（13）積分即可得到映射終點（N₂，θ₂）.

3）第3階段：N₂→N_d.類似于N₁→N_u的映射過程，映射表達式如下：

4）第4階段：N_d→N₁.類似于N_u→N₂過程，從（N_d，θ_d）出發的相軌跡必然只能沿著下穩定分支運動，同樣對式（13）積分即可得到映射終點（N₁，θ₁₀）.因此，可得到一維映射關系：R→R：（N₁，θ₁）→（N₁，θ₁₀）.

顯然，當系統存在普通平衡點時，從R=[Θ₁，Θ₂]出發的相軌跡，有可能被吸引至慢不變流形的上下穩定分支，只有通過多次映射后映射終點仍落于R區間時，才能出現穩定的極限環，從而產生強調制響應.

3.2數值討論與分析

當參數ε= 0.05、C = 4/3、ξ₂= 0.2、f = 0.5、σ=0時，對應Θ₁=-0.9110，Θ₂= 1.5114，各階段映射過程如圖9所示.圖9（a）對應第1階段映射過程，相軌跡從R=［Θ₁，Θ₂］出發跳躍至上穩定分支N處；圖9（b）對應第2階段映射過程，相軌跡通過上穩定分支慢變至N₂折疊線處；圖9（c）對應第3階段映射過程，相軌跡從N₂折疊線跳躍至下穩定分支；圖9（d）對應第4階段映射過程，相軌跡通過下穩定分支慢變至N₁折疊線處.由圖9可知，經過一次完整的映射過程后，從R=［Θ₁，Θ₂］出發的相軌跡最終都會落在該區間內，且新得到的相位區間向內收縮，表明無論經過多少次映射，終點都會位于R=［Θ₁，Θ₂］內.因此，在該參數條件下，系統能在θ₁≈0處出現穩定的極限環.

進一步考慮在相同參數條件下，通過局部分岔得到的幅頻特性曲線，如圖10所示.圖10中，HP_i（i=1，2，3）表示第i個Hopf分岔點；SN_j（j=1，2，3，4）表示j個鞍結分岔點.由圖10可知，當σ=0時，系統不存在穩定的周期平衡點，從R=［Θ₁，Θ₂］出發的相軌跡無法被普通周期平衡點吸引，都能產生穩定的強調制響應，這也驗證了一維映射的分析結果.

當σ= 2.5時，系統存在一個幅值較小的穩定周期平衡點，與σ= 0相比，其一維映射過程也發生了明顯的變化，如圖11所示.由圖11可知，從R出發的相軌跡在第4階段映射過程中有一部分被下穩定分支上的周期平衡點吸引，從而不能返回至R，在經過一次或多次映射過程后能夠返回至R區間的映射終點較映射起始點而言整體向右偏移，表明在經過一次或多次映射過程后，系統響應將從強調制響應狀態逃逸；從R區間出發的相軌跡都會逐漸吸引至穩定周期吸引子，從而不能產生穩定極限環，強調制響應也不會持續發生.

由圖12可知，當超過該頻率區域時，不能產生強調制響應，同時根據一維映射圖可以大致判斷產生強調制響應的相角范圍，該相角范圍代表觸發強調制響應的系統初始條件；若在頻率邊界附近產生強調制響應的初始條件范圍大大縮小，則意味著出現強調制響應的條件將會更加苛刻.

綜上所述，NES系統能產生穩定的強調制響應必須滿足以下條件：

2）必須形成N₁→N_u→N₂→N_d→N₁的連續跳躍環路，且軌線不被周期平衡點吸引.

顯然以上兩個條件對系統參數、外界激勵以及初始條件都有嚴格的要求，特別是在多解共存區域，需要特定的初始條件才能激發穩定的強調制響應. 同時，在某些激勵頻率區域，系統只存在穩定的強調制響應.仍以σ=0為例，強調制響應的穩定極限環與慢變系統數值計算結果如圖13所示.由圖13可知：一是穩定極限環對應的初始相位為θ₁≈-0.05；二是一維映射結果與實際系統響應吻合良好，可清晰看到調制響應的兩個慢變過程與兩個快變過程.

4強調制響應檢測電路設計與試驗

采用Tina-Ti軟件中模擬電路開發模塊對特定參數下的NES系統進行仿真計算，并開展強調制響應檢測電路試驗對理論分析結論進行驗證.

4.1電路原理圖設計

考慮式（3）所示的動力學系統，將其改寫為四維狀態方程：

依據微分狀態方程與電路狀態方程的等價關系，強調制響應檢測電路主要由積分電路、反相比例電路、乘法器等幾個關鍵模塊組成，其中運算放大器選用UA741封裝芯片，而乘法器選用具有8引腳的AD633JN封裝芯片，通過2個AD633JN乘法器串聯來實現狀態方程中的立方項.Tina-Ti仿真軟件中模擬電路的原理如圖14所示，各元器件工作電壓為-15～15 V，輸入余弦電壓信號幅值為25 mV，頻率為159 mHz，輸出電壓測點布置為VF1與VF2，分別與式（19）中的x₁、x₂相對應.

4.2電路仿真與試驗結果分析

利用Tina-Ti電路仿真軟件中的瞬態響應求解器得到兩個測點的電壓響應，并與式（19）系統參數分別取ε= 0.05、C= 2、ξ₂= 0.2、f = 0.5、σ=0，初始條件均為0時的數值計算結果進行對比，如圖15所示.由圖15可知，VF1測點與VF2測點均呈現穩定的強調制響應，這與數值計算結果中位移響應x₁、x₂的變化趨勢保持一致，存在的略微差異是由于輸入頻率只能近似取整為159 mHz，導致相位上存在一定延遲.

依據電路原理圖制作電路板，并開展電路試驗. 其中信號發生器用于提供外界余弦激勵信號，穩壓電源設置為±14.9V，用于提供元器件工作電壓，示波器用于采集測點VF1和VF2的電壓信號，相互連接關系以及電路試驗結果分別如圖16和圖17所示. 由圖可知，試驗測試結果與電路仿真結果以及Runge-Kutta數值方法結果吻合良好，表明強調制響應檢測電路是有效的，也驗證了強調制響應在NES系統中是真實存在的.

5結論

利用復變量平均法推導了主共振下NES系統的慢變動力流方程，通過多尺度法研究快變與慢變兩個時間尺度上系統平衡點的特性，從慢不變流形相軌跡、慢變流形特征以及松弛振子一維映射等多視角探究了NES系統產生強調制響應的充要條件，通過電路仿真與試驗驗證了理論分析的正確性.得出如下結論：

到的平衡點.

3）強調制響應檢測電路仿真和試驗結果與理論分析結果保持一致，驗證了NES系統中強調制響應存在的真實性和靶能量傳遞的可行性.

參考文獻

[1] ZHANG Y W，LU Y N，ZHANG W，et al. Nonlinear energy sink with inerter [J]. Mechanical Systems and Signal Processing，2019，125：52-64.

[2]吳長青，張志田，張偉峰.考慮幾何非線性的橋梁后顫振極限環特性[J].湖南大學學報（自然科學版），2018，45（5）：1-10.

WU C Q，ZHANG Z T，ZHANG W F. Limit cycle oscillation of bridge post-flutter with geometric nonlinearities [J]. Journal of Hunan University （Natural Sciences），2018，45（5）：1-10. （In Chinese ）

[3]劉良坤，譚平，陳洋洋，等.非線性能量阱減振系統受基底簡諧激勵的強調制反應分析[J].北京工業大學學報，2019，45（2）：177-185.

LIU L K，TAN P，CHEN Y Y，et al. Strongly modulated response analysis of nonlinear energy sink absorption system under ground harmonic excitation [J]. Journal of Beijing University of Technology，2019，45（2）：177-185.（In Chinese）

[4] XIONG L Y，TANG L H，LIU K F，et al. Broadband piezoelectric vibration energy harvesting using a nonlinear energy sink[J]. Journal of Physics D：Applied Physics，2018，51（18）：185502.

[5] WANG J J，WANG B，WIERSCHEM N E，et al. Dynamic analysis of track nonlinear energy sinks subjected to simple and stochastice excitations [J]. Earthquake Engineering & Structural Dynamics，2020，49（9）：863-883.

[6]JIANG X A，MCFARLAND D M，BERGMAN L A，et al. Steady state passive nonlinear energy pumping in coupled oscillators：theoretical and experimental results[J]. Nonlinear Dynamics，2003，33：87-102.

[7]ZHANG Y W，ZHANG H，HOU S，et al. Vibration suppression of composite laminated plate with nonlinear energy sink[J]. Acta Astronautica，2016，123：109-115.

[8]STAROSVETSKY Y，GENDELMAN O V. Bifurcations of attractors in forced system with nonlinear energy sink：the effect of mass asymmetry[J].Nonlinear Dynamics，2010，59（4）：711-731.

[9]GOURDON E，LAMARQUE C H. Nonlinear energy sink with uncertain parameters[J]. Journal of Computational and Nonlinear Dynamics，2006，1（3）：187-195.

[10]李爽，樓京俊，柴凱，等.激勵力幅值對非線性能量阱系統全局分岔特性的影響研究[J].船舶力學，2021，25（4）：479-488.

LI S，LOU J J，CHAI K，et al. Research on the influence of excitation amplitude on global bifurcation characteristics of nonlinear energy sink system[J]. Journal of Ship Mechanics，2021，25（4）：479-488.（In Chinese）

[11]張也弛，孔憲仁，張紅亮.非線性耦合振子間的靶能量傳遞研究：保守系統中的完全能量傳遞[J].振動與沖擊，2012，31 （1）：150-155.

ZHANG Y C，KONG X R，ZHANG H L. Targeted energy transfer among coupled nonlinear oscillators：complete energy exchange in a conservative system[J]. Journal of Vibration and Shock，2012，31（1）：150-155.（In Chinese）

[12] XIE L，BAGUET S，PRABEL B，et al. Bifurcation tracking by Harmonic Balance Method for performance tuning of nonlinear dynamical systems[J].Mechanical Systems and Signal Processing，2017，88：445-461.

[13] CHEN J E，HE W，ZHANG W，et al. Vibration suppression and higher branch responses of beam with parallel nonlinear energy sinks[J].Nonlinear Dynamics，2018，91（2）：885-904.