搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

经由脉冲式爆炸连接的复合式张弛振荡

宋锦 魏梦可 姜文安 张晓芳 韩修静 毕勤胜

经由脉冲式爆炸连接的复合式张弛振荡

宋锦, 魏梦可, 姜文安, 张晓芳, 韩修静, 毕勤胜
PDF
HTML
导出引用
导出核心图
  • 张弛振荡现象普遍存在于自然科学以及工程技术的各个领域, 探索张弛振荡的可能路径是张弛振荡研究的重要问题之一. 最近, 一种名为“脉冲式爆炸”(pulse-shaped explosion, PSE)的可以诱发张弛振荡的新机制被相继报道. PSE意味着平衡点和极限环表现出了与参数变化相关的脉冲式急剧量变, 这导致系统出现急剧转迁现象, 进而诱发张弛振荡. 本文以多频激励Mathieu-van der Pol-Duffing系统为例, 探讨了复合式的张弛振荡现象. 当参数激励和外部激励存在相位差时, 快子系统包含了两个不同的向量场部分, 由此得到了系统的双稳定特性. 特别地, 在狭小的参数范围内, 分岔会随着PSE的产生而产生, 这使得PSE更具复杂性. 基于此, 揭示了两种复合式的张弛振荡, 其特征是每一周期的演化过程包含了由PSE连接的两个张弛振荡簇. 我们的研究深化了对PSE及张弛振荡复杂动力学行为的理解.
      通信作者: 韩修静, xjhan@mail.ujs.edu.cn
    • 基金项目: 国家级-基于延迟分岔的混合模式振动的动力学机制与分类(11572141)
    [1]

    Belhaq M, Ghouli Z, Hamdi M 2018 Nonlinear Dyn. 94 2537

    [2]

    Mani A K, Narayanan M D, Sen M 2018 Nonlinear Dyn. 93 945

    [3]

    Miwadinou C H, Hinvi L A, Monwanou A V, Orou J B C 2017 Nonlinear Dyn. 88 97

    [4]

    Miwadinou C H, Monwanou A V, Yovogan J, Hinvi L A, Tuwa P R N, Orou J B C 2018 Chin. J. Phys. 56 1089

    [5]

    于文婷, 张娟, 唐军 2017 物理学报 66 7

    Yu W T, Zhang J, Tang J 2017 Acta Phys. Sin. 66 7

    [6]

    Ferraro M, D'annessa I, Moroni E, Morra G, Paladino A, Rinaldi S, Compostella, F, Colombo G 2019 J. Med. Chem. 62 60

    [7]

    谢献忠, 沈伟成, 彭剑, 龙昊 2016 防灾减灾工程学报 36 972

    XieX Z, Sheng W C, Peng J, Long H 2016 Journal of Disaster Prevention and Mitigation Engineering 36 972

    [8]

    Tuttle T D, Seering W P 1996 IEEE Trans. Rob. Autom. 12 368

    [9]

    梁晓冰, 刘希顺, 刘安芝, 王博亮 2009 物理学报 58 5065

    Liang X B, Liu X S, Liu A Z, Wang B L 2009 Acta Phys. Sin. 58 5065

    [10]

    Li Y X, Rinzel J 1994 J. Theor. Biol. 166 461

    [11]

    Newby J M, Bressloff P C, Keener J P 2013 Phys. Rev. Lett. 111 128101

    [12]

    Zhang J J, Chen D Y, Zhang H, Xu B B, Li H H, Aggidis G A, Chatterton S 2019 J. Vib. Control 25 2863

    [13]

    He J H, Ji F Y 2019 Therm. Sci. 23 2131

    [14]

    Deng B 2004 Chaos 14 1083

    [15]

    Roberts A 2016 SIAM J. Appl.Dyn. Syst. 15 609

    [16]

    Krupa M, Touboul J D 2016 J. Dyn. Differ. Equ. 28 471

    [17]

    Golmakani A, Homburg A J 2011 Dynam Syst. 26 61

    [18]

    Cheng W, Xiang Z 2019 J. Differ. Equ. 267 3397

    [19]

    Nikolay D D, Stepan A T 2019 Appl. Math. Lett. 96 208

    [20]

    LI JZ, Tomsovic S 2019 Phys. Rev. E 99 032212

    [21]

    Han X J, Xia F B, Ji P, Bi Q S, Kurths J 2016 Commun. Nonlinear Sci. Numer. Simul. 36 517

    [22]

    Cao J D, Guerrini L, Cheng Z S 2019 Appl. Math. Comput. 343 21

    [23]

    Li L, Wang Z, Li Y X, Shen H, Lu J W 2018 Appl. Math. Comput. 330 152

    [24]

    Han M A, Zhang W N 2010 J. Differ. Equ. 248 2399

    [25]

    Jiang H P, Zhang T H, Song Y L 2015 Int. J. Bifurcation Chaos 25 1550058

    [26]

    Monica C, PItchaimani M 2016 Nonlinear Anal. Real World Appl. 27 55

    [27]

    Wang Z, Wang X H, Li Y X, Huang X 2017 Int. J. Bifurcation Chaos 27 1750209

    [28]

    魏梦可, 韩修静, 张晓芳, 毕勤胜 2019 力学学报 51 904

    Wei M k, Han X J, Zhang X F, Bi Q S 2019 Chinese Journal of Theoretical and Applied Mechanics 51 904

    [29]

    Han X J, Bi Q S, Kurths J 2018 Phys. Rev. E 98 010201

    [30]

    Mohamad M A, Sapsis T P 2016 Ocean Eng. 120 289

    [31]

    Lakrad F, Azouani A, Abouhazim N, Belhaq M 2005 Chaos, Solitons Fractals 24 813

    [32]

    Wiggers V, Rech P C 2018 Eur. Phys. J. B 91 144

    [33]

    Shukla A K, Ramamohan T R, Srinivas S 2014 Phys. Scr. 89 075202

    [34]

    Kumar P, Kumar A, Erlicher S 2017 Physica D 358 1

    [35]

    Tamba V K, Kingni S T, Kuiate G F, Fotsin H B, Talla P K 2018 Pramana-J Phys. 91 12

    [36]

    He J H 2005 Chaos, Solitons Fractals 26 695

    [37]

    He J H 1999 Int. J. NonLinear Mech. 34 699

    [38]

    He J H, Abdou M A 2007 Chaos, Solitons Fractals 34 1421

    [39]

    He J H 2019 Results in Physics 15 102546

    [40]

    Li Y J, Wu Z Q, Wang F, Zhang G Q, Wang Y C 2019 J. Low Freq. Noise Vib. Act. Control UNSP 1461348419878534

    [41]

    He J H 2014 Indian J. Phys. 88 193

    [42]

    Han X J, Bi Q.S, Ji P, Kurths J 2015 Phys. Rev. E 92 012911

    [43]

    Sprott J C, Li C B 2014 Phys. Rev. E 89 066901

    [44]

    Arumugam R, Sarkar S, Banerjee T, Sinha S, Dutta P S 2019 Phys. Rev. E 99 012911

  • 图 1  系统(1)中典型的复合式张弛振荡 (a) $\alpha = 1.5$; (b) $\alpha = 0.4$; (c) $\alpha = 0.2$. 其他参数固定在$\gamma = 4$, $\delta = 1.00$, ${\beta _1} = 0.99$, ${\beta _2} = 1$, $\omega = 0.01$$\theta = - {{\text{π}}}/{2}$

    Fig. 1.  Typical compound relaxation oscillations in system (1): (a) $\alpha = 1.5$; (b) $\alpha = 0.4$; (c) $\alpha = 0.2$. Other parameters are fixed at $\gamma = 4$, $\delta = 1$, ${\beta _1} = 0.99$, ${\beta _2} = 1$, $\omega = 0.01$ and $\theta = - {{\text{π}}}/{2}$.

    图 2  (a)子系统(2a)和(b)子系统(2b)在参数平面$(w, \alpha )$上的分岔集. 其中GH为广义Hopf分岔SubH为亚临界Hopf分岔, SupH为超临界Hopf分岔, LPC为极限环的分岔. 系统参数的取值与图1相同

    Fig. 2.  Bifurcation sets of the subsystem (2a) (a) and (2b) (b) in the parameter plane $(w, \alpha )$. Here GH represent the generalized Hopf bifurcation, SubH represent the subcritical Hopf bifurcation, SupH represent the supercritical Hopf bifurcation, LPC represent the limit point cycle bifurcation. The values of system parameters are the same as those in Fig. 1.

    图 3  快子系统(2)在A, B, C各区域中典型的稳定性和分岔行为 (a1), (a2) $\alpha = 1.5$; (b1), (b2) $\alpha = 0.4$; (c1), (c2) $\alpha = 0.2$. 其他参数的取值与图1相同

    Fig. 3.  Typical stability and bifurcation behaviors of the fast subsystem (2) in the areas A, B and C: (a1), (a2) $\alpha = 1.5$; (b1), (b2) $\alpha = 0.4$; (c1), (c2) $\alpha = 0.2$. The values of other parameters are the same as those in Fig. 1.

    图 4  子系统(2b)在参数平面$\left( {w, {\beta _{\rm{1}}}} \right)$上的分岔集. 其他参数的取值与图1(a)相同

    Fig. 4.  Bifurcation sets of the subsystem (2b) in the parameter plane $\left( {w, {\beta _{\rm{1}}}} \right)$. The values of other parameters are the same as those in Fig. 1(a).

    图 5  为子系统(2b)的分岔图 (a) ${\beta _1} = 1.1$; (b) ${\beta _1} = 1$. 其他参数的取值与图1(a)相同

    Fig. 5.  Bifurcation diagrams of the subsystem (2b): (a) ${\beta _1} = $1.1; (b) ${\beta _1} = 1$. The values of other parameters are the same as those in Fig. 1(a).

    图 6  图1(a)中的张弛振荡的快慢分析 (a)张弛振荡的转换相图与图3(a1)中的分岔图的叠加(与子系统(2a)相关); (b)张弛振荡的转换相图与图3(a2)中分岔图的叠加(与子系统(2b)相关); (c)一个完整周期下的张弛振荡. 这里$\alpha = 1.5$, 而其他参数与图1相同

    Fig. 6.  Fast-slow analysis of the relaxation oscillations in Fig. 1(a): (a) Overlay of the transformed phase diagram of the relaxation oscillations and the bifurcation diagram in Fig. 3(a1) (related to the subsystem (2a)); (b) overlay of the transformed phase diagram of the relaxation oscillations and the bifurcation diagram in Fig. 3(a2) (related to the subsystem (2b)); (c) a whole period of the relaxation oscillations. Here $\alpha = 1.5$and other parameters are the same as those in Fig. 1.

    图 7  图1(b)中的张弛振荡的快慢分析

    Fig. 7.  Fast-slow analysis of the relaxation oscillations in Fig. 1(b).

    图 8  图1(c)中的张弛振荡的快慢分析

    Fig. 8.  Fast-slow analysis of the relaxation oscillations in Fig. 1(c).

  • [1]

    Belhaq M, Ghouli Z, Hamdi M 2018 Nonlinear Dyn. 94 2537

    [2]

    Mani A K, Narayanan M D, Sen M 2018 Nonlinear Dyn. 93 945

    [3]

    Miwadinou C H, Hinvi L A, Monwanou A V, Orou J B C 2017 Nonlinear Dyn. 88 97

    [4]

    Miwadinou C H, Monwanou A V, Yovogan J, Hinvi L A, Tuwa P R N, Orou J B C 2018 Chin. J. Phys. 56 1089

    [5]

    于文婷, 张娟, 唐军 2017 物理学报 66 7

    Yu W T, Zhang J, Tang J 2017 Acta Phys. Sin. 66 7

    [6]

    Ferraro M, D'annessa I, Moroni E, Morra G, Paladino A, Rinaldi S, Compostella, F, Colombo G 2019 J. Med. Chem. 62 60

    [7]

    谢献忠, 沈伟成, 彭剑, 龙昊 2016 防灾减灾工程学报 36 972

    XieX Z, Sheng W C, Peng J, Long H 2016 Journal of Disaster Prevention and Mitigation Engineering 36 972

    [8]

    Tuttle T D, Seering W P 1996 IEEE Trans. Rob. Autom. 12 368

    [9]

    梁晓冰, 刘希顺, 刘安芝, 王博亮 2009 物理学报 58 5065

    Liang X B, Liu X S, Liu A Z, Wang B L 2009 Acta Phys. Sin. 58 5065

    [10]

    Li Y X, Rinzel J 1994 J. Theor. Biol. 166 461

    [11]

    Newby J M, Bressloff P C, Keener J P 2013 Phys. Rev. Lett. 111 128101

    [12]

    Zhang J J, Chen D Y, Zhang H, Xu B B, Li H H, Aggidis G A, Chatterton S 2019 J. Vib. Control 25 2863

    [13]

    He J H, Ji F Y 2019 Therm. Sci. 23 2131

    [14]

    Deng B 2004 Chaos 14 1083

    [15]

    Roberts A 2016 SIAM J. Appl.Dyn. Syst. 15 609

    [16]

    Krupa M, Touboul J D 2016 J. Dyn. Differ. Equ. 28 471

    [17]

    Golmakani A, Homburg A J 2011 Dynam Syst. 26 61

    [18]

    Cheng W, Xiang Z 2019 J. Differ. Equ. 267 3397

    [19]

    Nikolay D D, Stepan A T 2019 Appl. Math. Lett. 96 208

    [20]

    LI JZ, Tomsovic S 2019 Phys. Rev. E 99 032212

    [21]

    Han X J, Xia F B, Ji P, Bi Q S, Kurths J 2016 Commun. Nonlinear Sci. Numer. Simul. 36 517

    [22]

    Cao J D, Guerrini L, Cheng Z S 2019 Appl. Math. Comput. 343 21

    [23]

    Li L, Wang Z, Li Y X, Shen H, Lu J W 2018 Appl. Math. Comput. 330 152

    [24]

    Han M A, Zhang W N 2010 J. Differ. Equ. 248 2399

    [25]

    Jiang H P, Zhang T H, Song Y L 2015 Int. J. Bifurcation Chaos 25 1550058

    [26]

    Monica C, PItchaimani M 2016 Nonlinear Anal. Real World Appl. 27 55

    [27]

    Wang Z, Wang X H, Li Y X, Huang X 2017 Int. J. Bifurcation Chaos 27 1750209

    [28]

    魏梦可, 韩修静, 张晓芳, 毕勤胜 2019 力学学报 51 904

    Wei M k, Han X J, Zhang X F, Bi Q S 2019 Chinese Journal of Theoretical and Applied Mechanics 51 904

    [29]

    Han X J, Bi Q S, Kurths J 2018 Phys. Rev. E 98 010201

    [30]

    Mohamad M A, Sapsis T P 2016 Ocean Eng. 120 289

    [31]

    Lakrad F, Azouani A, Abouhazim N, Belhaq M 2005 Chaos, Solitons Fractals 24 813

    [32]

    Wiggers V, Rech P C 2018 Eur. Phys. J. B 91 144

    [33]

    Shukla A K, Ramamohan T R, Srinivas S 2014 Phys. Scr. 89 075202

    [34]

    Kumar P, Kumar A, Erlicher S 2017 Physica D 358 1

    [35]

    Tamba V K, Kingni S T, Kuiate G F, Fotsin H B, Talla P K 2018 Pramana-J Phys. 91 12

    [36]

    He J H 2005 Chaos, Solitons Fractals 26 695

    [37]

    He J H 1999 Int. J. NonLinear Mech. 34 699

    [38]

    He J H, Abdou M A 2007 Chaos, Solitons Fractals 34 1421

    [39]

    He J H 2019 Results in Physics 15 102546

    [40]

    Li Y J, Wu Z Q, Wang F, Zhang G Q, Wang Y C 2019 J. Low Freq. Noise Vib. Act. Control UNSP 1461348419878534

    [41]

    He J H 2014 Indian J. Phys. 88 193

    [42]

    Han X J, Bi Q.S, Ji P, Kurths J 2015 Phys. Rev. E 92 012911

    [43]

    Sprott J C, Li C B 2014 Phys. Rev. E 89 066901

    [44]

    Arumugam R, Sarkar S, Banerjee T, Sinha S, Dutta P S 2019 Phys. Rev. E 99 012911

  • [1] 邢雅清, 陈小可, 张正娣, 毕勤胜. 多平衡态下簇发振荡产生机理及吸引子结构分析. 物理学报, 2016, 65(9): 090501. doi: 10.7498/aps.65.090501
    [2] 姜海波, 李涛, 曾小亮, 张丽萍. 周期脉冲作用下Logistic映射的复杂动力学行为及其分岔分析. 物理学报, 2013, 62(12): 120508. doi: 10.7498/aps.62.120508
    [3] 吴天一, 陈小可, 张正娣, 张晓芳, 毕勤胜. 非对称型簇发振荡吸引子结构及其机理分析. 物理学报, 2017, 66(11): 110501. doi: 10.7498/aps.66.110501
    [4] 张正娣, 刘亚楠, 李静, 毕勤胜. 分段Filippov系统的簇发振荡及擦边运动机理. 物理学报, 2018, 67(11): 110501. doi: 10.7498/aps.67.20172421
    [5] 巩华荣, 唐昌建, 王文祥, 魏彦玉, 黄民智, 宫玉彬. 微波管中离子张弛振荡的混沌现象. 物理学报, 2005, 54(1): 159-163. doi: 10.7498/aps.54.159
    [6] 丁晓玲, 王健, 王旭明, 何大韧. 一类不可逆保守系统中的混沌类吸引子. 物理学报, 2002, 51(3): 482-486. doi: 10.7498/aps.51.482
    [7] 何希哲, 邹立寿. 张弛限累振荡电路分析及其阻抗检验. 物理学报, 1975, 143(3): 223-230. doi: 10.7498/aps.24.223
    [8] 高超, 毕勤胜, 张正娣. 一个跃变电路切换系统的振荡行为及分岔机理分析. 物理学报, 2013, 62(2): 020504. doi: 10.7498/aps.62.020504
    [9] 唐驾时, 梁翠香, 萧寒. 单频外激励弹簧摆的鞍结分岔控制. 物理学报, 2009, 58(5): 2989-2995. doi: 10.7498/aps.58.2989
    [10] 吴天一, 张正娣, 毕勤胜. 切换电路系统的振荡行为及其非光滑分岔机理. 物理学报, 2012, 61(7): 070502. doi: 10.7498/aps.61.070502
    [11] 李绍龙, 张正娣, 吴天一, 毕勤胜. 广义BVP电路系统的振荡行为及其非光滑分岔机理. 物理学报, 2012, 61(6): 060504. doi: 10.7498/aps.61.060504
    [12] 张正娣, 刘杨, 张苏珍, 毕勤胜. 余维-1非光滑分岔下的簇发振荡及其机理. 物理学报, 2017, 66(2): 020501. doi: 10.7498/aps.66.020501
    [13] 季颖, 毕勤胜. 高维广义蔡氏电路中的快慢动力学行为及其分岔分析. 物理学报, 2012, 61(1): 010202. doi: 10.7498/aps.61.010202
    [14] 葛烨, 胡以华, 舒嵘, 洪光烈. 一种新型的用于差分吸收激光雷达中脉冲式光学参量振荡器的种子激光器的频率稳定方法. 物理学报, 2015, 64(2): 020702. doi: 10.7498/aps.64.020702
    [15] 于永吉, 陈薪羽, 王超, 吴春婷, 董渊, 李述涛, 金光勇. 基于MgO:APLN的多光参量振荡器实验研究及其逆转换过程演化分析. 物理学报, 2015, 64(4): 044203. doi: 10.7498/aps.64.044203
    [16] 余跃, 张春, 韩修静, 姜海波, 毕勤胜. 周期切换下Chen系统的振荡行为与非光滑分岔分析. 物理学报, 2013, 62(2): 020508. doi: 10.7498/aps.62.020508
    [17] 肖金标, 王登峰. 硅基槽式纳米线多模干涉型模阶数转换器全矢量分析. 物理学报, 2017, 66(7): 074203. doi: 10.7498/aps.66.074203
    [18] 姜海波, 张丽萍, 陈章耀, 毕勤胜. 脉冲作用下Chen系统的非光滑分岔分析. 物理学报, 2012, 61(8): 080505. doi: 10.7498/aps.61.080505
    [19] 季颖, 毕勤胜. 参外联合激励复合非线性振子的分岔分析. 物理学报, 2009, 58(7): 4431-4438. doi: 10.7498/aps.58.4431
    [20] 夏志林, 郭培涛, 薛亦渝, 黄才华, 李展望. 短脉冲激光诱导薄膜损伤的等离子体爆炸过程分析. 物理学报, 2010, 59(5): 3523-3530. doi: 10.7498/aps.59.3523
  • 引用本文:
    Citation:
计量
  • 文章访问数:  151
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-11-29
  • 修回日期:  2020-01-17

经由脉冲式爆炸连接的复合式张弛振荡

  • 江苏大学土木工程与力学学院, 镇江 212013
  • 通信作者: 韩修静, xjhan@mail.ujs.edu.cn
    基金项目: 国家级-基于延迟分岔的混合模式振动的动力学机制与分类(11572141)

摘要: 张弛振荡现象普遍存在于自然科学以及工程技术的各个领域, 探索张弛振荡的可能路径是张弛振荡研究的重要问题之一. 最近, 一种名为“脉冲式爆炸”(pulse-shaped explosion, PSE)的可以诱发张弛振荡的新机制被相继报道. PSE意味着平衡点和极限环表现出了与参数变化相关的脉冲式急剧量变, 这导致系统出现急剧转迁现象, 进而诱发张弛振荡. 本文以多频激励Mathieu-van der Pol-Duffing系统为例, 探讨了复合式的张弛振荡现象. 当参数激励和外部激励存在相位差时, 快子系统包含了两个不同的向量场部分, 由此得到了系统的双稳定特性. 特别地, 在狭小的参数范围内, 分岔会随着PSE的产生而产生, 这使得PSE更具复杂性. 基于此, 揭示了两种复合式的张弛振荡, 其特征是每一周期的演化过程包含了由PSE连接的两个张弛振荡簇. 我们的研究深化了对PSE及张弛振荡复杂动力学行为的理解.

English Abstract

    • 从生物到化学, 从物理到大气科学, 多尺度耦合效应问题[1-4]普遍存在, 例如神经元系统的信息传递[5]、生物代谢过程中的变构效应[6]、输电塔与塔线之间的耦合振动问题[7]以及减速器系统的复杂振动问题[8]等. 一般地, 多尺度耦合作用下的非线性系统往往能够表现出复杂的动力学特性, 如张弛振荡. 张弛振荡是一类复杂的振荡模式, 其特征是在每一周期的演化过程中可以观测到大幅振荡和小幅振荡的相互交替.

      张弛振荡的动力学机理问题是张弛振荡研究的重要问题之一. 从1963年诺贝尔奖获得者Hodgkin和Huxley[9,10]在研究神经元放电过程建立的数学模型开始, 张弛振荡就逐渐受到了学者们的关注, 而后Rinzel提出的快慢分析法[11,12]为研究张弛振荡提供了理论框架. 快慢分析法在使用时需将多尺度系统分解为维数较低的快、慢子系统. 然而, 在系统降维分解的过程中, 通常会涉及到信息损失, 损失的信息可采用双尺度数学[13]加以分析. 基于快慢分析法, 诸如Canard现象[14-16], Shilnikov同宿轨的失稳[17-20], Hopf分岔[21-24]以及经由延迟分岔的慢通道效应[25-27]等多种路径先后被揭示与张弛振荡的产生有关.

      “脉冲式爆炸”(pulse-shaped explosion, PSE)是最近报道的一种可以诱发张弛振荡的新机制, 其特征是在平衡点/极限环的解支上出现了与参数的变化相关的脉冲状的急剧量变[28]. 现有的研究表明, PSE的产生通常与不同激励的频率关系有关; 特别地, 激励频率比增加能够导致PSE数量的增加[29]. 本文考虑一个由Mathieu[30,31]和van der Pol-Duffing振荡器[32-35]耦合的非线性方程, 即Mathieu-van der Pol-Duffing方程(MVD):

      $\ddot x + [1 + {\beta _1}\cos (\omega t)](x + \delta {x^3}) - \gamma \dot x(\alpha - {x^2}) \!=\! {\beta _2}\cos (\omega t + \theta ),$

      其中${\beta _1}\cos (\omega t)$${\beta _2}\cos(\omega t + \theta )$为系统的参数激励和外部激励, 且激励频率$\omega $远小于系统的固有频率${\omega _0}$. 系统(1)是一类典型的非线性振动方程, 可采用同伦摄动法[36]、变分迭代法[37]以及指数函数法[38]等多种解析方法加以分析. 特别地, 幅频关系问题是非线性振动研究的重要问题, 系统(1)的振动频率可以采用文献[39]提出的最简方法加以估算. 此外, 注意到系统(1)含有多个参数, 参数的变化往往会引起系统行为的定性变化[40], 而不同参数之间的耦合作用关系可以考虑采用双参数同伦摄动法[41]加以探讨.

      本文焦注于系统(1)的PSE及其诱导的张弛振荡. 我们发现, 在一定的条件下系统会产生复合式的张弛振荡(见图1). 我们分析系统在相位差作用下的双稳定性和PSE现象, 由此揭示复合型张弛振荡的产生机制. 研究表明, 系统在每个周期的演化过程中, 先后产生了由PSE连接的两个张弛振荡簇, 由此形成了所谓的复合式的张弛振荡.

      图  1  系统(1)中典型的复合式张弛振荡 (a) $\alpha = 1.5$; (b) $\alpha = 0.4$; (c) $\alpha = 0.2$. 其他参数固定在$\gamma = 4$, $\delta = 1.00$, ${\beta _1} = 0.99$, ${\beta _2} = 1$, $\omega = 0.01$$\theta = - {{\text{π}}}/{2}$

      Figure 1.  Typical compound relaxation oscillations in system (1): (a) $\alpha = 1.5$; (b) $\alpha = 0.4$; (c) $\alpha = 0.2$. Other parameters are fixed at $\gamma = 4$, $\delta = 1$, ${\beta _1} = 0.99$, ${\beta _2} = 1$, $\omega = 0.01$ and $\theta = - {{\text{π}}}/{2}$.

    • 注意到本文所考虑的激励频率远远小于系统的固有频率. 因此, 系统的激励项$\cos (\omega t)$$\cos (\omega t + \theta )$在较慢的尺度上演化, 而原系统在较快的尺度上演化, 即系统(1)是一个典型的含有两个慢变量的快慢系统. 根据文献[42]提出的分析方法, 即频率转换快慢分析法, 可将系统(1)转化为仅含一个慢变量的快慢系统. 然后, 通过分析转化后的快慢系统, 进而可以揭示原系统中张弛振荡的动力学机制.

      下面, 采用频率转换快慢分析法对系统(1)进行分析. 令外部激励$\cos (\omega t + \theta ) = w$为基准慢变量. 对于固定的相位差$\theta = -{{\text{π}}}/{2}$, 可得$\cos (\omega t + \theta ) =$$ \sin (\omega t) = w $. 另一方面, 注意到$\cos (\omega t) = \pm \sqrt {1 - {w^2}} $, 故转化后的快慢系统的快子系统可以表示为:

      $\left\{ \begin{aligned} &\ddot x + (1 + {\beta _1}\sqrt {1 - {w^2}} )(x + \delta {x^3}) \\ &\;\; - \gamma \dot x(\alpha - {x^2}) = {\beta _2}w,\quad \quad \quad\quad\quad \quad \quad \quad\;\;\;\;\;\left( {\rm 2a} \right)\\ &\ddot x + (1 - {\beta _1}\sqrt {1 - {w^2}} )(x + \delta {x^3}) \\ & \;\; - \gamma \dot x(\alpha - {x^2}) = {\beta _2}w.\quad \quad\quad\quad\quad\quad \quad \quad\;\;\;\;\;\left( {\rm 2b} \right) \end{aligned} \right.$

      对于(2a)式来说, $\omega t \in \left( {2 n{\text{π}} - \dfrac{{\text{π}}}{2}, 2 n\pi + \dfrac{{\text{π}}}{{\rm{2}}}} \right)$, 其中n为整数; 而对于(2b)式来说, $\omega t \in \left(2 n{\text{π}} + \dfrac{{\text{π}}}{2}\right. $, $2 n{\text{π}} + \left.\dfrac{{3{\text{π}}}}{2} \right) $, 其中n为整数.

      显然, 快子系统(2)的平衡点可以表示为E(x, 0), 其中x由方程

      $\left\{ \begin{aligned} & (1 + {\beta _1}\sqrt {1 - {w^2}} )(x + \delta {x^3}) - {\beta _2}w = 0\;\;\;\;\;\;\;\;\;\;\;\left( {\rm 3a} \right)\\ & (1 - {\beta _1}\sqrt {1 - {w^2}} )(x + \delta {x^3}) - {\beta _2}w = 0\;\;\;\;\;\;\;\;\;\;\;\left( {\rm 3b} \right) \end{aligned} \right.$

      的实根决定. 在本文的研究中, 取定参数$\delta > 0$${\beta _1} > 0$. 在方程(3a)中, 注意到关系式

      $\varDelta = \frac{{\beta _2^2{w^2}}}{{4{\delta ^2}{{(1 + {\beta _1}\sqrt {1 - {w^2}} )}^2}}} + \frac{1}{{27{\delta ^3}}} > 0$

      恒成立, 因此子系统(2a)始终存在一个平衡点. 而对于方程(3b), 当$\left( {1 - {\beta _1}\sqrt {1 - {w^2}} } \right) \ne 0$时,

      $\varDelta = \frac{{\beta _2^2{w^2}}}{{4{\delta ^2}{{(1 - {\beta _1}\sqrt {1 - {w^2}} )}^2}}} + \frac{1}{{27{\delta ^3}}} > 0$

      恒成立. 此时, 子系统(2b)存在唯一的平衡点. 特别地, 当$\left( {1 - {\beta _1}\sqrt {1 - {w^2}} } \right) = 0$, 即$w = \pm \sqrt {\dfrac{{\beta _1^2 - 1}}{{\beta _1^2}}} $时, 快子系统(2b)不存在平衡点.

    • 当系统仅含一个吸引子时, 系统表现出所谓的单稳定性; 而当两个或两个以上吸引子共存时, 便得到了双稳定性或多稳定性[43,44]. 对于本文所考虑的系统来说, 在相位差的作用下, 其快子系统被分解为由方程(2a)和(2b)决定的两个部分. 而每个向量场部分都可能会产生吸引子, 因此相位差的存在可能会诱发快子系统的双稳定性.

      为了揭示相位差下的双稳定性, 作出快子系统(2)关于$(w, \alpha )$参数平面的分岔集(见图2), 系统参数的取值与图1相同. 此外, 为了清晰地展示各向量场部分的稳定性和分岔行为, 将分岔集切分为对应于子系统(2a)和(2b)的两个部分, 并分别绘制于图2(a)图2(b)中.

      图  2  (a)子系统(2a)和(b)子系统(2b)在参数平面$(w, \alpha )$上的分岔集. 其中GH为广义Hopf分岔SubH为亚临界Hopf分岔, SupH为超临界Hopf分岔, LPC为极限环的分岔. 系统参数的取值与图1相同

      Figure 2.  Bifurcation sets of the subsystem (2a) (a) and (2b) (b) in the parameter plane $(w, \alpha )$. Here GH represent the generalized Hopf bifurcation, SubH represent the subcritical Hopf bifurcation, SupH represent the supercritical Hopf bifurcation, LPC represent the limit point cycle bifurcation. The values of system parameters are the same as those in Fig. 1.

      图2所示, 参数平面$\left( {w, \alpha } \right)$被直线$\alpha = 0.466$$\alpha = 0.333$划分为A, B, C三个区域. 首先, 考虑参数$\alpha $属于区域A的情形(即$\alpha > 0.466$). 以$\alpha = 1.5$时的情形为例, 当w从–1开始逐渐增大时(对应图2(a)), 系统仅含唯一的一个吸引子, 即极限环吸引子(见图3(a1)), 且该吸引子无分岔行为的发生. 当w增大到1之后, 便逐渐减小. 此时, 系统从向量场(2a)切换到(2b). 随着w的逐渐减小, w将依次穿越8条分岔曲线, 发生8次分岔(见图3(a2)).

      图  3  快子系统(2)在A, B, C各区域中典型的稳定性和分岔行为 (a1), (a2) $\alpha = 1.5$; (b1), (b2) $\alpha = 0.4$; (c1), (c2) $\alpha = 0.2$. 其他参数的取值与图1相同

      Figure 3.  Typical stability and bifurcation behaviors of the fast subsystem (2) in the areas A, B and C: (a1), (a2) $\alpha = 1.5$; (b1), (b2) $\alpha = 0.4$; (c1), (c2) $\alpha = 0.2$. The values of other parameters are the same as those in Fig. 1.

      然后, 以$\alpha = 0.4$为例, 考虑参数$\alpha $属于区域B时(即$0.333 < \alpha < 0.466$)的情形. 当w从–1逐渐增大时, 系统将先后产生4次分岔, 即极限环的fold分岔, 亚临界Hopf分岔, 亚临界Hopf分岔, 极限环的fold分岔(见图3(b1)). 随后, w开始减小, 系统切换到了向量场(2b). 接着, 与区域B中的4条分岔曲线相关的4种分岔行为依次发生(见图3(b2)).

      最后, 考虑参数$\alpha $属于区域C时(即$\alpha < 0.333$)的情形. 以$\alpha = 0.2$为例, 当系统由向量场(2a)支配时, w将两次穿越超临界Hopf分岔值(见图3(c1)); 而当系统由向量场(2b)支配时, w也将两次穿越超临界Hopf分岔值(见图3(c2)).

      综上, 相位差的存在使得快子系统被切割为(2a)和(2b)两部分. 图3给出了子系统(2a)和(2b)在各参数区域中典型的稳定性和分岔行为. 可以发现, 不论是子系统(2a), 还是(2b), 它们均表现出了单稳态的动力学特性. 另一方面, 注意到子系统(2a)和(2b)联合构成了快子系统(2), 因此快子系统(2)呈现出了因相位差的存在而诱发的双稳定性.

    • 考虑到系统(2b)中, ${\beta _1}$在1附近时存在临界峰值. 为了更加深入地揭示系统的PSE行为, 固定参数$\alpha = 1.5$, 其他系统参数同图1(a). 对转换后的快子系统(2b)作关于$\left( {w, {\beta _1}} \right)$参数平面的分岔集(如图4所示), 的确, 当${\beta _1}$在1附近时系统会产生不同的分岔行为. 为了定性分析快子系统(2b)的分岔行为, 固定其他参数, 仅改变${\beta _1}$, 分别取定${\beta _1}$分别为1.1, 1, 0.99作关于$\left( {w, x} \right)$相平面的分岔图, 如图5(a),图5(b)图3(a2)所示.

      图  4  子系统(2b)在参数平面$\left( {w, {\beta _{\rm{1}}}} \right)$上的分岔集. 其他参数的取值与图1(a)相同

      Figure 4.  Bifurcation sets of the subsystem (2b) in the parameter plane $\left( {w, {\beta _{\rm{1}}}} \right)$. The values of other parameters are the same as those in Fig. 1(a).

      图  5  为子系统(2b)的分岔图 (a) ${\beta _1} = 1.1$; (b) ${\beta _1} = 1$. 其他参数的取值与图1(a)相同

      Figure 5.  Bifurcation diagrams of the subsystem (2b): (a) ${\beta _1} = $1.1; (b) ${\beta _1} = 1$. The values of other parameters are the same as those in Fig. 1(a).

      由于临界峰值的存在, 且临界值为$w = \pm \sqrt {\dfrac{{\beta _1^2 - 1}}{{\beta _1^2}}} $, 代入可知, 当${\beta _1} = 1.1$时, 临界值为$w = \pm 0.4166$, 当稳定的平衡线出现在临界值附近时, 平衡点发生急剧的转迁, 在极短的时间, 极短的空间内迅速趋于无穷, 这导致出现急剧陡峭的平衡点迅速离开原曲线, 也称吸引子的急速逃逸. 而当${\beta _1} = 1$时, 临界值为$w = 0$, 此时系统存在一个发散点, 当稳定的平衡线靠近临界值0时, 平衡点发生急剧的转迁, 平衡点迅速趋于无穷, 且临界线左右方向相反. 当${\beta _1} = 0.99$时, 发散区域消失, 平衡线合二为一. 即在给定系统参数情况下, 仅改变参数${\beta _1}$从1.1到1再到0.99时, 穿越三个不同区域, 在临界峰值处, 左右两条平衡线逐渐合二为一, 产生了正负双向脉冲式爆炸, 即正负双向PSE现象.

    • 已经分析了相位差下系统的双稳定性和PSE现象. 本部分探讨与此相关的张弛振荡的产生. 两类复合式的张弛振荡模式, 即“subHopf/fold-cycle”型和“subHopf/subHopf”型将被揭示.

    • 由于快子系统被划分为对应于不同动力学行为的三个参数区域, 因此当参数取在不同的区域时可能会产生不同的张弛振荡模式. 本部分探讨复合式subHopf/fold-cycle型张弛振荡, 它与参数$\alpha $取在区域A和B有关. 首先考虑$\alpha $取在区域A的情形, 即情形A.

      情形A: 为了便于分析, 固定$\gamma = 4$, $\delta = 1$, $\omega = 0.01$, ${\beta _1} = 0.99$, ${\beta _2} = 1$, 图2给出了$\left( {w, \alpha } \right)$参数平面上的分岔集. 如图2所示, 当$\alpha > 0.466$时, 在所考虑的参数间隔内不同的$\alpha $不会产生定性的变化. 因此取定$\alpha = 1.5$为情形A. 通过数值模拟可得到时间历程图, 如图1(a)所示, 在每个周期内, 此时系统表现为: 在两个大幅振荡簇之间存在一个正负双向PSE.

      为了更好地揭示该系统的动力学行为, 对系统(1)进行快慢分析并引入转换相图, 令$\cos (\omega t + \theta ) = $$\sin (\omega t) = w$, 由于在该多频激励系统中, 慢变参数均可以用关于w的代数式表示. 从而原系统又可以表示为快子系统(2). 将$\sin (\omega t) = w$作为分岔参数, 在$(w, x)$相平面上作分岔图与慢子系统的转换相图的叠加图, 由于在该多频激励系统中还存在着参数激励$\cos (\omega t)$, 在用含有w的参数表示时, 原系统的分岔图及转换相图应由两部分组成, 如图6所示. 在系统(2)中, $\alpha = 1.5$为对应的情形A.

      图  6  图1(a)中的张弛振荡的快慢分析 (a)张弛振荡的转换相图与图3(a1)中的分岔图的叠加(与子系统(2a)相关); (b)张弛振荡的转换相图与图3(a2)中分岔图的叠加(与子系统(2b)相关); (c)一个完整周期下的张弛振荡. 这里$\alpha = 1.5$, 而其他参数与图1相同

      Figure 6.  Fast-slow analysis of the relaxation oscillations in Fig. 1(a): (a) Overlay of the transformed phase diagram of the relaxation oscillations and the bifurcation diagram in Fig. 3(a1) (related to the subsystem (2a)); (b) overlay of the transformed phase diagram of the relaxation oscillations and the bifurcation diagram in Fig. 3(a2) (related to the subsystem (2b)); (c) a whole period of the relaxation oscillations. Here $\alpha = 1.5$and other parameters are the same as those in Fig. 1.

      结合单周期时间历程图6(c)可知, 在系统(2a)中, 当慢变量$\sin (0.01 t)$从–1到1逐渐增大时, 对应时间历程图6(c)左侧, 轨线始终被稳定的极限环吸引产生大幅振荡(见图6(a)). 如图6(b)所示, 当慢变量$\sin (0.01 t)$达到其最大值1时, 慢变量开始减小, 对应时间历程图5(c)右侧. 结合图6(b)所示, 即在图6(b)中最右侧, 轨线仍旧被稳定的极限环吸引产生大幅振荡. 当慢变量减小到LPC1分岔点时, 极限环失稳脱离原来的轨道, 随着轨线继续运动, 运动到SubH1分岔点, 轨线由极限环吸引子转迁到平衡点吸引子. 极限环吸引子消失, 从而导致了大幅振荡的消失和小幅振荡的开始, 即系统由激发态向沉寂态转迁. 特别地, 分析结果表明, 此时平衡点类型为结点, 故系统轨线向极限环收敛的速度非常快, 之后轨线沿着稳定的平衡点运动, 根据数值模拟的分岔情况可知, 轨线应该在SubH1分岔点沿着稳定的平衡线向左运动后应该接着运动到SubH2分岔点时, 开始沿着不稳定的极限环起振直至运动到LPC2分岔点产生极限环的fold分岔. 吸引子变成稳定的极限环, 继续向左运动. 当慢变继续减小轨线运动到LPC3分岔点, 经由LPC3分岔点产生不稳定的极限环, 随后跳向SubH3分岔点, 系统由激发态向沉寂态转迁. 轨线应脱离极限环继续向左运动. 然而由于SubH2分岔点及SubH3分岔点间隔时间非常短, 轨线还未来得及跳上极限环, 轨线就已经穿过该区域, 又迅速转迁到平衡线上. 即系统产生了正负双向PSE. 平衡线急剧转迁, 直接越过该区域到达SubH3分岔点, 并运动到稳定的平衡线上. 随着慢变量的继续减小, 轨线运动到SubH4处开始起振, 经由LPC4点后, 产生极限环的fold分岔, 轨线跳到稳定的极限环上, 但由于“慢通道效应”, 轨线并没有立刻起振, 而是继续运动一段时间后起振到极限环上, 到$\sin (0.01 t)$最小值时, 完成一次循环, 开始下个周期. 也可以说正负双向脉冲式爆炸连接了subHopf/fold cycle型脉冲式爆炸.

      情形B: 当参数$\alpha $取在区域B时, 即其他参数同情形A仅改变$\alpha = 0.4$时, 在参数变化范围内, 改变$\alpha $的取值不会产生定性的变化. 通过数值模拟得到了系统分岔图与转换相图的叠加, 如图7所示.

      图  7  图1(b)中的张弛振荡的快慢分析

      Figure 7.  Fast-slow analysis of the relaxation oscillations in Fig. 1(b).

      情形B与情形A相似但不同, 相似点表现在, 情形B定性分析与情形A均为subHopf/fold-cycle型复合式张弛振荡. 不同点具体表现为以下两个方面: 一方面, 当$\sin (0.01 t)$从–1到1逐渐增大的过程中, 从图7(a)可知, 轨线沿着稳定的平衡线运动, 经由“慢通道效应”导致的延迟后, 就已经发生并完成了SubHopf/fold-cycle分岔. 不同的是, 在情形A图6(a)中, 轨线仅仅沿着稳定的极限环运动, 系统的稳定性并没有发生任何改变, 直到$\sin (0.01 t)$从1到–1逐渐减小的过程中才产生分岔, 然后发生了PSE现象; 另一方面, 从图2就可看出, 当w从–1到1增大的过程中, 在情形A中并没有产生任何分岔, 而情形B中已经产生了4次分岔行为, 从单个周期时间历程图6(c)图7(c)中也可以很明显发现, 在每个周期内, 情形B中经由PSE连接的大幅振荡簇的起振范围及起振区域比情形A都更小.

    • 图2可知, 当参数$\alpha $落在区域C(即$\alpha < 0.333$)时, 慢变量可以穿越不同的分岔曲线, 由此诱发不同的分岔模式和复杂动力学行为. 本部分探讨与此相关的张弛振荡. 当参数$\alpha $落在区域C时, 不同的$\alpha $不会导致定性不同的分岔行为. 因此, 不失一般性, 在下面的讨论中固定参数$\alpha = 0.2$.

      $\alpha = 0.2$时, 作关于$(w, x)$相平面的分岔图与转换相图叠加图(见图8(a)图8(b)). 在图8(a)中, 当慢变量$\sin (0.01 t)$逐渐增大时, 从时间历程图8(c)可知, 轨线从最左侧沿着稳定的平衡线运动, 然后运动到SupH1分岔点, 由于“慢通道效应”, 轨线并未沿着稳定的平衡线运动, 而是接着沿着平衡线运动很长一段时间后才开始起振, 形成稳定的极限环, 然后运动到SubH2分岔点后, 再经过一段时间后从稳定的极限环上下来, 运动到平衡线上. 然后到最右侧当慢变量$\sin (0.01 t)$达到其最大值1时, 结合时间历程图8(c)可知, 慢变量开始减小, 即轨线从图7(b2)最右端沿着稳定的平衡线向左运动. 一直向左运动到SubH3分岔点和SubH4分岔点, 与区域A和区域B相同, 由于两点之间间隔非常短, 轨线还未来得及起振就已经跳到SupH4分岔点左侧的平衡线处, 沿着稳定的平衡线向左运动, 当$\sin (0.01 t)$达到其最小值–1时, 完成一次循环. 即正负双向脉冲式爆炸连接了subHopf/subHopf型张弛振荡.

      图  8  图1(c)中的张弛振荡的快慢分析

      Figure 8.  Fast-slow analysis of the relaxation oscillations in Fig. 1(c).

    • 探讨张弛振荡的动力学机制是张弛振荡研究的重要问题之一. PSE作为一种诱发张弛振荡的新路径, 吸引着许多科研工作者. 在以往的研究中, 不同激励的频率关系被认为是诱发PSE的重要因素, 例如频率关系的改变会导致系统向量场的变化, 这使得极速逃逸现象的数量增多, 进而导致PSE数量的增加, 并由此诱发复杂的张弛振荡. 本文的研究表明, 当频率相等的两个激励存在相位差时, 会导致快子系统由两个不同的向量场部分组成, 这诱发了系统的双稳定性和PSE现象. 特别地, 其中一个向量场部分展现了关于PSE对称的两组分岔行为, 且每组分岔行为可以诱发一个张弛振荡簇. 基于此, 得到了经由PSE连接的复合式的张弛振荡模式. 本文的研究基于一个特定的非线性系统. 然而, 由本文的分析可知, 当两个慢变激励存在相位差时, 必然会导致快子系统包含两个不同的向量场部分. 因此, 对于其他系统来说, 有可能也会诱发双稳定性和PSE现象. 综上所述, 本文报道的相位差下的双稳定性、PSE现象以及基于此而产生的复杂的张弛振荡模式也可能在其他系统中被观测到.

参考文献 (44)

目录

    /

    返回文章
    返回