搜索

x

留言板

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

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

形状记忆合金薄板系统全局激变现象分析

岳晓乐 向以琳 张莹

形状记忆合金薄板系统全局激变现象分析

岳晓乐, 向以琳, 张莹
PDF
HTML
导出引用
导出核心图
  • 形状记忆合金在工程应用中的难点主要来自于系统在温度和外载荷作用下产生的复杂全局动力学行为. 本文以形状记忆合金薄板动力系统为研究对象, 分析在温度和激励振幅两个控制参数作用下系统的全局动力学. 通过全局分岔图, 可以观测到系统会发生复杂的激变现象, 然后利用复合胞坐标系方法, 获取系统的吸引子、吸引域、鞍和域边界等信息, 展现系统的全局演变过程. 研究发现, 系统随着振幅和温度变化会呈现复杂的全局结构, 并发生一系列的边界激变、合并激变现象, 同时多次发生分形-Wada, Wada-Wada, Wada-分形等域边界突变. 通过对指定区域细化, 可以清晰地显示域边界的分形特征. 研究结果对于如何通过调控温度与外载荷强度, 使形状记忆合金薄板在系统中发挥最佳性能具有理论指导意义.
      通信作者: 岳晓乐, xiaoleyue@nwpu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11672230, 11672232)和陕西省自然科学基础研究计划(批准号: 2017JM1029)资助的课题.
    [1]

    Yuan B, Zhu M, Chung C 2018 Materials 11 1716

    [2]

    Hartl D J, Lagoudas D C 2007 Proc. Inst. Mech. Eng. Part G: J. Aerosp. Eng. 221 535

    [3]

    Lee J, Jin M, Ahn K K 2013 Mechatronics 23 310

    [4]

    Jani J M, Leary M, Subic A, Gibson M A 2014 Mater. Des. 56 1078

    [5]

    Song G, Ma N, Li H N 2006 Eng. Struct. 28 1266

    [6]

    Bernardini D, Rega G 2011 Int. J. Bifurcation Chaos 21 2769

    [7]

    Paula A S, Savi M A, Lagoudas D C 2012 J. Braz. Soc. Mech. Sci. Eng. 34 401

    [8]

    Sado D, Pietrzakowski M 2010 Int. J. Non-Linear Mech. 45 859

    [9]

    Hashemi S M T, Khadem S E 2006 Int. J. Mech. Sci. 48 44

    [10]

    Savi M A 2015 Int. J. Non-Linear Mech. 70 2

    [11]

    Han Q, Xu W, Yue X 2014 Int. J. Bifurcation Chaos 24 1450051

    [12]

    Grebogi C, Ott E, Yorke J A 1982 Phys. Rev. Lett. 48 1507

    [13]

    Grebogi C, Ott E, Yorke J A 1983 Physica D 7 181

    [14]

    Chian A C L, Borotto F A, Rempel E L, Rogers C 2005 Chaos Solitons Fractals 24 869

    [15]

    Yue X, Xu W, Zhang Y 2012 Nonlinear Dyn. 69 437

    [16]

    刘莉, 徐伟, 岳晓乐, 韩群 2013 物理学报 62 200501

    Liu L, Xu W, Yue X L, Han Q 2013 Acta Phys. Sin. 62 200501

    [17]

    刘晓君, 洪灵, 江俊 2016 物理学报 65 180502

    Liu X J, Hong L, Jiang J 2016 Acta Phys. Sin. 65 180502

    [18]

    Yue X, Xu W, Wang L 2013 Commun. Nonlinear Sci. Numer. Simul. 18 3567

    [19]

    Zhang Y 2013 Phys. Lett. A 377 1269

    [20]

    Hsu C S 1992 Int. J. Bifurcation Chaos 2 727

    [21]

    Hsu C S 1995 Int. J. Bifurcation Chaos 5 1085

    [22]

    Tongue B H 1987 Physica D 28 401

    [23]

    Falk F 1980 Acta Metall. Sin. 28 1773

    [24]

    Machado L G, Savi M A, Pacheco P M C L 2004 Shock Vib. 11 67

    [25]

    黄志华, 刘平, 杜长城, 李映辉 2009 力学季刊 30 71

    Huang Z H, Liu P, Du C C, Li Y H 2009 Chin. Quarterly Mech. 30 71

  • 图 1  SMA薄板的横向振动

    Fig. 1.  Transverse vibration of SMA thin plate.

    图 2  系统(7)随振幅g变化的多值分岔图

    Fig. 2.  Multivalued bifurcation diagram of the system (7) with the variation of amplitude g.

    图 3  系统(7)的全局图 (a) g = 0.0471; (b) g = 0.0472

    Fig. 3.  Global diagram of the system (7): (a) g = 0.0471; (b) g = 0.0472.

    图 4  系统(7)的全局图 (a) g = 0.0482; (b) g = 0.0483

    Fig. 4.  Global diagram of the system (7): (a) g = 0.0482; (b) g = 0.0483.

    图 5  系统(7)的全局图 (a) g = 0.0491; (b) g = 0.0492

    Fig. 5.  Global diagram of the system (7): (a) g = 0.0491; (b) g = 0.0492.

    图 6  系统(7)的全局图 (a) g = 0.0490; (b) g = 0.0491

    Fig. 6.  Global diagram of the system (7): (a) g = 0.0490; (b) g = 0.0491.

    图 7  系统(7)的全局图 (a) g = 0.0596; (b) g = 0.0597; (c) g = 0.0726; (d) g = 0.0727

    Fig. 7.  Global diagram of the system (7): (a) g = 0.0596; (b) g = 0.0597; (c) g = 0.0726; (d) g = 0.0727.

    图 8  系统(7)的全局图 (a) g = 0.0460; (b) g = 0.0461

    Fig. 8.  Global diagram of the system (7): (a) g = 0.0460; (b) g = 0.0461.

    图 9  系统(7)的全局图 (a) g = 0.0478; (b) g = 0.0479; (c), (d) 分别对应于(a), (b)图的区域细化

    Fig. 9.  Global diagram of the system (7): (a) g = 0.0478; (b) g = 0.0479; (c), (d) the region refinement of panels (a) and (b).

    图 10  系统(7)的全局图 (a) g = 0.0533; (b) g = 0.0534

    Fig. 10.  Global diagram of the system (7): (a) g = 0.0533; (b) g = 0.0534.

    图 11  系统(7)随温度θ变化的多值分岔图

    Fig. 11.  Multivalued bifurcation diagram of the system (7) with the variation of temperature θ.

    图 12  系统(7)的全局图 (a) θ = 0.8379; (b) θ = 0.8380

    Fig. 12.  Global diagram of the system (7): (a) θ = 0.8379; (b) θ = 0.08380.

    图 13  系统(7)的全局图 (a) θ = 0.4088; (b) θ = 0.4089; (c) θ = 0.7182; (d) θ = 0.7183

    Fig. 13.  Global diagram of the system (7): (a) θ = 0.4088; (b) θ = 0.4089; (c) θ = 0.7182; (d) θ = 0.7183.

    图 14  系统(7)的全局图 (a) θ = 0.0950; (b) θ = 0.0951

    Fig. 14.  Global diagram of the system (7): (a) θ = 0.0950; (b) θ = 0.0951.

  • [1]

    Yuan B, Zhu M, Chung C 2018 Materials 11 1716

    [2]

    Hartl D J, Lagoudas D C 2007 Proc. Inst. Mech. Eng. Part G: J. Aerosp. Eng. 221 535

    [3]

    Lee J, Jin M, Ahn K K 2013 Mechatronics 23 310

    [4]

    Jani J M, Leary M, Subic A, Gibson M A 2014 Mater. Des. 56 1078

    [5]

    Song G, Ma N, Li H N 2006 Eng. Struct. 28 1266

    [6]

    Bernardini D, Rega G 2011 Int. J. Bifurcation Chaos 21 2769

    [7]

    Paula A S, Savi M A, Lagoudas D C 2012 J. Braz. Soc. Mech. Sci. Eng. 34 401

    [8]

    Sado D, Pietrzakowski M 2010 Int. J. Non-Linear Mech. 45 859

    [9]

    Hashemi S M T, Khadem S E 2006 Int. J. Mech. Sci. 48 44

    [10]

    Savi M A 2015 Int. J. Non-Linear Mech. 70 2

    [11]

    Han Q, Xu W, Yue X 2014 Int. J. Bifurcation Chaos 24 1450051

    [12]

    Grebogi C, Ott E, Yorke J A 1982 Phys. Rev. Lett. 48 1507

    [13]

    Grebogi C, Ott E, Yorke J A 1983 Physica D 7 181

    [14]

    Chian A C L, Borotto F A, Rempel E L, Rogers C 2005 Chaos Solitons Fractals 24 869

    [15]

    Yue X, Xu W, Zhang Y 2012 Nonlinear Dyn. 69 437

    [16]

    刘莉, 徐伟, 岳晓乐, 韩群 2013 物理学报 62 200501

    Liu L, Xu W, Yue X L, Han Q 2013 Acta Phys. Sin. 62 200501

    [17]

    刘晓君, 洪灵, 江俊 2016 物理学报 65 180502

    Liu X J, Hong L, Jiang J 2016 Acta Phys. Sin. 65 180502

    [18]

    Yue X, Xu W, Wang L 2013 Commun. Nonlinear Sci. Numer. Simul. 18 3567

    [19]

    Zhang Y 2013 Phys. Lett. A 377 1269

    [20]

    Hsu C S 1992 Int. J. Bifurcation Chaos 2 727

    [21]

    Hsu C S 1995 Int. J. Bifurcation Chaos 5 1085

    [22]

    Tongue B H 1987 Physica D 28 401

    [23]

    Falk F 1980 Acta Metall. Sin. 28 1773

    [24]

    Machado L G, Savi M A, Pacheco P M C L 2004 Shock Vib. 11 67

    [25]

    黄志华, 刘平, 杜长城, 李映辉 2009 力学季刊 30 71

    Huang Z H, Liu P, Du C C, Li Y H 2009 Chin. Quarterly Mech. 30 71

  • [1] 洪灵, 徐健学. 一类新的边界激变现象:混沌的边界激变. 物理学报, 2001, 50(4): 612-618. doi: 10.7498/aps.50.612
    [2] 孙光爱, 王虹, 汪小琳, 陈波, 常丽丽, 刘耀光, 盛六四, Woo W, Kang MY. 原位中子衍射研究两相NiTi合金的微力学相互作用和相变机理. 物理学报, 2012, 61(22): 226102. doi: 10.7498/aps.61.226102
    [3] 马明全, 王文秀, 何大韧. 一个二维分段光滑映象中的边界激变. 物理学报, 2000, 49(9): 1679-1682. doi: 10.7498/aps.49.1679
    [4] 王清周, 陆东梅, 崔春翔, 韩福生. 利用内耗研究淬火空位对Cu-11.9Al-2.5Mn(wt%)形状记忆合金逆马氏体相变温度的影响. 物理学报, 2008, 57(11): 7083-7087. doi: 10.7498/aps.57.7083
    [5] 张莹, 雷佑铭, 方同. 混沌吸引子的对称破缺激变. 物理学报, 2009, 58(6): 3799-3805. doi: 10.7498/aps.58.3799
    [6] 何 阅, 姜玉梅, 申 影, 何大韧. 一个分段连续系统中的胖奇异集激变. 物理学报, 2005, 54(3): 1071-1080. doi: 10.7498/aps.54.1071
    [7] 巢小刚, 戴 俊, 王文秀, 何大韧. 一个强、弱耗散功能可分隔系统中的特征激变. 物理学报, 2006, 55(1): 47-53. doi: 10.7498/aps.55.47
    [8] 刘莉, 徐伟, 岳晓乐, 韩群. 一类含非黏滞阻尼的Duffing单边碰撞系统的激变研究. 物理学报, 2013, 62(20): 200501. doi: 10.7498/aps.62.200501
    [9] 洪 灵, 徐健学. 常微分方程系统中内部激变现象的研究. 物理学报, 2000, 49(7): 1228-1234. doi: 10.7498/aps.49.1228
    [10] 刘晓君, 洪灵, 江俊. 非自治分数阶Duffing系统的激变现象. 物理学报, 2016, 65(18): 180502. doi: 10.7498/aps.65.180502
    [11] 陆云清, 王文秀, 何大韧. 一个电张弛振子中的瞬态激变. 物理学报, 2003, 52(5): 1079-1084. doi: 10.7498/aps.52.1079
    [12] 屈世显, 孙建刚, 官山, 纪西萍, 何大韧, 汪秉宏. 动力学不连续性引起的新型激变. 物理学报, 1996, 45(12): 1970-1977. doi: 10.7498/aps.45.1970
    [13] 夏彬凯, 李剑锋, 李卫华, 张红东, 邱枫. 基于离散变分原理的耗散动力学模拟方法:模拟三维囊泡形状. 物理学报, 2013, 62(24): 248701. doi: 10.7498/aps.62.248701
    [14] 刘洪涛, 孙光爱, 王沿东, 陈波, 汪小琳. 冲击诱发NiTi形状记忆合金相变行为研究. 物理学报, 2013, 62(1): 018103. doi: 10.7498/aps.62.018103
    [15] 罗诗裕, 邵明珠, 韦洛霞, 刘曾荣. 位错动力学与系统的全局分叉. 物理学报, 2004, 53(6): 1940-1945. doi: 10.7498/aps.53.1940
    [16] 刘俊明, 刘治国, 吴状春. CuCo合金失稳分解动力学. 物理学报, 1995, 44(2): 319-327. doi: 10.7498/aps.44.319
    [17] 刘耀民, 刘中良, 黄玲艳. 分形理论结合相变动力学的冷表面结霜过程模拟. 物理学报, 2010, 59(11): 7991-7997. doi: 10.7498/aps.59.7991
    [18] 张永伟, 唐刚, 韩奎, 寻之朋, 谢裕颖, 李炎. 分形基底上刻蚀模型动力学标度行为的数值模拟研究. 物理学报, 2012, 61(2): 020511. doi: 10.7498/aps.61.020511
    [19] 杨毅, 唐刚, 宋丽建, 寻之朋, 夏辉, 郝大鹏. 分形基底上受限固-固模型动力学性质的数值模拟研究. 物理学报, 2014, 63(15): 150501. doi: 10.7498/aps.63.150501
    [20] 张程宾, 程启坤, 陈永平. 分形结构纳米复合材料热导率的分子动力学模拟研究. 物理学报, 2014, 63(23): 236601. doi: 10.7498/aps.63.236601
  • 引用本文:
    Citation:
计量
  • 文章访问数:  799
  • PDF下载量:  10
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-27
  • 修回日期:  2019-06-06
  • 上网日期:  2019-09-01
  • 刊出日期:  2019-09-20

形状记忆合金薄板系统全局激变现象分析

    基金项目: 国家自然科学基金(批准号: 11672230, 11672232)和陕西省自然科学基础研究计划(批准号: 2017JM1029)资助的课题.

摘要: 形状记忆合金在工程应用中的难点主要来自于系统在温度和外载荷作用下产生的复杂全局动力学行为. 本文以形状记忆合金薄板动力系统为研究对象, 分析在温度和激励振幅两个控制参数作用下系统的全局动力学. 通过全局分岔图, 可以观测到系统会发生复杂的激变现象, 然后利用复合胞坐标系方法, 获取系统的吸引子、吸引域、鞍和域边界等信息, 展现系统的全局演变过程. 研究发现, 系统随着振幅和温度变化会呈现复杂的全局结构, 并发生一系列的边界激变、合并激变现象, 同时多次发生分形-Wada, Wada-Wada, Wada-分形等域边界突变. 通过对指定区域细化, 可以清晰地显示域边界的分形特征. 研究结果对于如何通过调控温度与外载荷强度, 使形状记忆合金薄板在系统中发挥最佳性能具有理论指导意义.

English Abstract

    • 形状记忆合金(shape memory alloy, SMA)表现出的独特力学性能来源于马氏体相变及其逆变化, 相变的驱动力可以由温度变化和机械载荷提供, 温度诱发的相变引起形状记忆效应, 外加应力诱导的相变产生超弹性. 利用这两个性质, SMA被制成各种参数可控的智能元件, 在机翼结构的变形、精密机械系统的高速驱动及机器人仿生等方面得到广泛应用[15]. 目前关于SMA元件的力学性能研究主要集中在温度和载荷变化下SMA梁、振子及支架系统的平衡点稳定性、分岔行为及混沌现象等非线性动力学特征方面[610], 对于SMA薄板系统的全局动力学研究较少, 相较于局部分析, 全局分析能够揭示更多的动力学信息[11], 有利于从力学理论角度突破SMA薄板系统在机械驱动和振动控制等领域的局限性.

      激变[12,13]作为混沌系统中较为常见的全局分岔现象, 主要刻画了混沌吸引子和混沌鞍的不连续变化. 常见的激变包括边界激变、内部激变和合并激变. 边界激变是指混沌吸引子与吸引域边界上的不稳定周期轨道碰撞, 导致吸引子的突然消失. 内部激变是混沌吸引子与其所在吸引域内的不稳定周期轨道发生碰撞引起的. 合并激变是指两个及以上的吸引子同时与吸引域边界上的不稳定周期轨道发生碰撞, 合并形成新的吸引子. 尽管通过分岔图可以观测到激变现象, 但具体类型并不清楚, 需进一步绘制全局图进行判断. 近年来, 不少系统的激变现象均有所研究, 例如分数阶和单边碰撞系统[1417].

      当系统中存在多个吸引子时, 吸引域边界可以是光滑的, 也可以呈现出“你中有我, 我中有你”的分形结构. 如果吸引子的个数为三个及以上, 且边界上任一点的任意小领域, 覆盖三个及以上吸引域, 则称此域边界具有Wada特性. 随着系统参数变化, 域边界结构会突然发生改变, 称为域边界突变, 包括光滑-Wada域边界突变、分形-Wada域边界突变和Wada-Wada域边界突变等. 域边界突变作为混沌动力学的一个研究重点, 对于确定系统的整体结构具有重要的物理意义, 在非光滑系统和准周期强迫系统[18,19]中都有所研究, 但分形性使得动力系统的力学行为很难预测, 随着胞映射技术的发展, 域边界突变和激变现象的研究得到了突破. 其中经典的广义胞映射、图胞映射、插值胞映射等方法均可获得系统的全局信息[2022], 但在刻画域边界时有一定的局限性. 复合胞坐标系方法通过对连续相空间的多级分割构造一个复合胞空间, 基于点映射的原理, 不仅可以获取动力系统的吸引子、吸引域和鞍等信息, 还能细化任意小区域, 在全局图中优化对域边界的刻画.

      本文以SMA薄板动力系统为研究对象, 温度和外部激励的振幅为分岔参数, 在全局分岔图的基础上, 通过复合胞坐标系方法进一步分析系统在演变过程中出现的激变类型和域边界突变现象, 并通过对指定区域的细化, 展示域边界的分形特征. 研究结果在工程领域中对控制系统的动态响应有重要意义.

    • 在SMA薄板中, 相对于奥氏体, 马氏体更具延展性, 当外加应力低于马氏体的屈服强度时, SMA薄板可表现出超弹性. 马氏体的转变与温度相关, 低温相称为马氏体, 高温相称为奥氏体, 降温时SMA薄板发生马氏体相变, 呈现形状记忆效应. 定义TM为马氏体相变临界温度, 该温度之下, 马氏体稳定; TA为奥氏体相变临界温度, 该温度之上, 奥氏体稳定, a, b, c均为材料常数, 满足:

      ${T_{\rm{A}}} = {T_{\rm{M}}} + \frac{{{b^2}}}{{4ac}}.$

      根据Falk [23]的研究, SMA薄板的本构模型可以由一个五次多项式表示, 该模型给出了应变(ε)-应力(σ)的关系:

      $\sigma = a(T - {T_{\rm{M}}})\varepsilon - b{\varepsilon ^3} + c{\varepsilon ^5}.$

      Savi等[10,24]给出了很多关于SMA弹簧振子的数学模型, 文献[25]在此基础上演变出SMA薄板横向振动的非线性动力学模型. 本文考虑密度ρ, 厚度h, 长宽分别为l1, l2的SMA薄板, 在环境黏性阻尼e下, 受到横向简谐激振力F0 = FsinΩt作用, 如图1所示.

      图  1  SMA薄板的横向振动

      Figure 1.  Transverse vibration of SMA thin plate.

      薄板的横向位移为w (x, y, t), 单位长度上的内力矩分别为Ux(x, y, t), Uy(x, y, t), Uxy(x, y, t), 其横向振动方程为[25]

      $\frac{{{\partial ^2}{U_x}}}{{\partial {x^2}}} + 2\frac{{{\partial ^2}{U_{xy}}}}{{\partial x\partial y}} + \frac{{{\partial ^2}{U_y}}}{{\partial {y^2}}} + {F_0} - e\dot w - \rho h\frac{{{\partial ^2}w}}{{\partial {t^2}}} = 0.$

      薄板面内的3个应力分别为σx, σy, τxy, 内力矩为[25]

      $\begin{array}{l} {U_x} = \displaystyle\int_{{{ - h} / 2}}^{{h / 2}} {{\sigma _x}z{\rm{d}}z},\; {U_y} = \displaystyle\int_{{{ - h} / 2}}^{{h / 2}} {{\sigma _y}z{\rm{d}}z},\\ {U_{xy}} = \displaystyle\int_{{{ - h} / 2}}^{{h / 2}} {{\tau _{xy}}z{\rm{d}}z} . \end{array}$

      在满足四周简支边界条件下取位移模式[25]:

      $w(x,y,t) = W(t)\sin \frac{{{\text{π}}x}}{{{l_1}}}\sin \frac{{{\text{π}}y}}{{{l_2}}}.$

      通过Galerkin原理和应变-应力本构关系, 并对系统参数量纲归一化, 得到SMA薄板的横向振动非线性动力学方程[25]:

      $\ddot W + \xi \dot W + \alpha (\theta - 1)W - \beta {W^3} + \gamma {W^5} = g{\rm{sin}}\varOmega t.$

      方程(6)可化为以下形式的一阶微分方程组:

      $\left\{ \begin{array}{l} \dot x = y, \\ \dot y = - \xi y - \alpha (\theta - 1)x + \beta {x^3} - \gamma {x^5} + g{\rm sin}\varOmega t, \\ \end{array} \right.$

      其中α, ξ, β, γ均为正常数; θ, g分别是无量纲化的温度和简谐激振力振幅. 参数选取同文献[24]所示, ξ = 0.2, α = 1, β = 1300, γ = 470000, Ω = 1. 接下来基于复合胞坐标系方法, 分析温度θ和激励振幅g两个参数变化时, SMA薄板系统的全局动力学特性变化过程.

    • 当温度θ = 0.7, 此时SMA薄板位于马氏体临界温度下, 含有两个稳定的马氏体相. 选取不同的初值点(–0.0553, –0.09), (–0.031, 0.045), (0.083, –0.032), (0.047, 0.059)构造Poincare映射, 保留稳态映射点, 得到对应的多值分岔图, 如图2所示. 随着幅值g值的变化, 系统呈现明显的多吸引子共存特征, 并伴随有丰富的激变和域边界突变等现象.

      图  2  系统(7)随振幅g变化的多值分岔图

      Figure 2.  Multivalued bifurcation diagram of the system (7) with the variation of amplitude g.

      在区域D = {(x, y)|–0.1 ≤ x ≤ 0.1, –0.12 ≤ y ≤ 0.12}上, 均匀划分1000 × 1200个胞, 利用复合胞坐标系方法获取系统的吸引子、吸引域和边界鞍等全局信息, 阐述g在0.0460到0.0727的范围内激变现象出现的机理. 为更直观理解全局图, 对以下图形做以下说明: 用△表示吸引子A, 不同颜色代表不同的吸引子; 吸引域用B表示, 不同吸引域颜色不同; S表示边界鞍; IS表示吸引域内部鞍; 下标表示吸引子、吸引域和内部鞍的个数, 例如A1, A2表示系统两个吸引子共存, B1, B2为其对应的吸引域.

    • g从0.0471增大到0.0472时, 发生两次逆的边界激变, 如图3所示. A5A6表示周期为2的吸引子, A5与嵌在吸引域B1, B5边界的鞍S碰撞, A6和嵌入吸引域B2, B6边界的鞍S碰撞, 使得A5, A6和部分边界鞍S突然消失, 成为内部周期鞍IS1, IS2, 同时吸引域B5, B6消失, 吸引域B1, B2变大.

      图  3  系统(7)的全局图 (a) g = 0.0471; (b) g = 0.0472

      Figure 3.  Global diagram of the system (7): (a) g = 0.0471; (b) g = 0.0472.

      g从0.0482变化到0.0483时, 系统再次发生两次边界激变, 两个周期1吸引子A1A2与边界上的鞍S发生碰撞, 变为吸引域B5, B6的内部鞍IS5, IS6, 同时吸引域B1B2消失, 如图4所示.

      图  4  系统(7)的全局图 (a) g = 0.0482; (b) g = 0.0483

      Figure 4.  Global diagram of the system (7): (a) g = 0.0482; (b) g = 0.0483.

      g从0.0491增大到0.0492时, 混沌吸引子A2与吸引域边界上的混沌鞍S发生碰撞, 发生边界激变, 变成新的更大的边界混沌鞍, 吸引域B2随之消失, 状态空间中为四个周期吸引子和混沌边界鞍共存, 如图5所示.

      图  5  系统(7)的全局图 (a) g = 0.0491; (b) g = 0.0492

      Figure 5.  Global diagram of the system (7): (a) g = 0.0491; (b) g = 0.0492.

    • g从0.0490变化到0.0491时, 混沌吸引子A2A4不断接近吸引域边界上的混沌鞍S, 发生合并激变, 成为新的混沌吸引子A2, 与此同时吸引域B2, B4合并为新的吸引域B2, 如图6所示.

      图  6  系统(7)的全局图 (a) g = 0.0490; (b) g = 0.0491

      Figure 6.  Global diagram of the system (7): (a) g = 0.0490; (b) g = 0.0491.

      g从0.0596变化到0.0597时, 混沌吸引子A1, A2与边界上的混沌鞍S碰撞, 发生合并激变, 变为新的混沌吸引子A1, 如图7(a)图7(b)所示. 当g从0.0726增大到0.0727时, 发生逆合并激变, 混沌吸引子A1消失, 出现两个新的周期1吸引子A1, A2及混沌边界鞍S, 如图7(c)图7(d)所示.

      图  7  系统(7)的全局图 (a) g = 0.0596; (b) g = 0.0597; (c) g = 0.0726; (d) g = 0.0727

      Figure 7.  Global diagram of the system (7): (a) g = 0.0596; (b) g = 0.0597; (c) g = 0.0726; (d) g = 0.0727.

    • g = 0.0460时, 系统存在三个周期1吸引子A1, A2, A3和嵌在吸引域B1, B2, B3边界上的混沌鞍S, 此时的域边界呈现出Wada特性. 当g增大到0.0461时, 系统新出现一个周期3吸引子A4, 此时的域边界由4个吸引域构成, 仍具有Wada特性, 系统发生Wada-Wada域边界突变, 如图8所示.

      图  8  系统(7)的全局图 (a) g = 0.0460; (b) g = 0.0461

      Figure 8.  Global diagram of the system (7): (a) g = 0.0460; (b) g = 0.0461.

      g从0.0478增大到0.0479时, 吸引子的个数从6个变为7个, 为判断吸引子变化过程, 对区域v ={(x, y)|–0.075 ≤ x ≤ –0.055, 0.068 ≤ y ≤ 0.092}进行细化. 可以发现, 原周期3吸引子A4消失, 并在其附近出现两个新的周期3吸引子A4A7, 如图9(c)图9(d). 原吸引域B4分裂为新吸引域B4B7, 且参数变化前后域边界均呈现Wada特性, 域边界结构更加复杂, 系统再次发生Wada-Wada域边界突变.

      图  9  系统(7)的全局图 (a) g = 0.0478; (b) g = 0.0479; (c), (d) 分别对应于(a), (b)图的区域细化

      Figure 9.  Global diagram of the system (7): (a) g = 0.0478; (b) g = 0.0479; (c), (d) the region refinement of panels (a) and (b).

      g = 0.0533时, 状态空间中有4个周期1吸引子共存, 此时域边界仍具有Wada特性. 当g增大到0.0534时, 周期吸引子A3A4消失, 状态空间中仅剩2个周期1吸引子A1, A2, 此时域边界的Wada特性消失, 只呈现分形特性, 系统发生Wada-分形域边界突变, 如图10所示.

      图  10  系统(7)的全局图 (a) g = 0.0533; (b) g = 0.0534

      Figure 10.  Global diagram of the system (7): (a) g = 0.0533; (b) g = 0.0534.

    • 温度作为SMA薄板的一个可控参数, 利用温度改变SMA薄板的力学特性在工程领域中有重要运用. 取g为0.06, 分析温度θ变化对系统的影响, 选取不同的初值点(–0.0553, –0.09), (–0.031, 0.045), (0.083, –0.032), (0.047, 0.059)构造Poincare映射, 得到对应的分岔图, 如图11所示. 可以发现, 随温度θ的变化, 系统中吸引子个数、类型及大小会发生改变, 并出现激变现象. 为从全局角度分析温度变化对系统激变的影响, 本节基于复合胞坐标系方法, 将区域D = {(x, y)| –0.1 ≤ x ≤ 0.1, –0.12 ≤ y ≤ 0.12}, 均匀划分为1000 × 1200个胞, 获得系统的吸引子、吸引域、鞍和域边界等全局特性.

      图  11  系统(7)随温度θ变化的多值分岔图

      Figure 11.  Multivalued bifurcation diagram of the system (7) with the variation of temperature θ.

    • θ = 0.8379时, 状态空间中两个周期1吸引子A1A2共存, S为嵌入在分形域边界上的混沌鞍, IS1, IS2是吸引域B1, B2内部的混沌鞍. 当θ增大为0.8380时, 系统出现两个新的吸引子A3A4, 内部鞍IS1, IS2消失, 系统发生两次逆边界激变, 如图12所示.

      图  12  系统(7)的全局图 (a) θ = 0.8379; (b) θ = 0.8380

      Figure 12.  Global diagram of the system (7): (a) θ = 0.8379; (b) θ = 0.08380.

    • θ从0.4088变化到0.4089时, 周期2吸引子A1A2同时与域边界上的混沌鞍S发生碰撞, 合并为一个新的混沌吸引子A1, 吸引域B1B2合并成为新吸引域B1, 系统发生合并激变, 如图13(a)图13(b)所示. 当θ从0.7182增大到0.7183时, 混沌吸引子A1消失, 出现两个新的混沌吸引子A1A2, 以及嵌在域边界上的混沌鞍S, 系统发生逆的合并激变, 如图13(c)图13(d)所示.

      图  13  系统(7)的全局图 (a) θ = 0.4088; (b) θ = 0.4089; (c) θ = 0.7182; (d) θ = 0.7183

      Figure 13.  Global diagram of the system (7): (a) θ = 0.4088; (b) θ = 0.4089; (c) θ = 0.7182; (d) θ = 0.7183.

    • g = 0.0950时, 状态空间中两个周期2吸引子A1, A2, 和两个周期1吸引子A3, A4共存, 吸引域边界呈现Wada特性. 当g增大为0.0951时, 两个周期1吸引子A3A4消失, 吸引域边界变为由B1B2构成的分形边界, 系统发生Wada-分形域边界突变, 如图14所示.

      图  14  系统(7)的全局图 (a) θ = 0.0950; (b) θ = 0.0951

      Figure 14.  Global diagram of the system (7): (a) θ = 0.0950; (b) θ = 0.0951.

    • 考虑实际工程中温度和应力对于SMA的力学特性的影响, 本文选取SMA薄板为研究对象, 以温度θ和激振力振幅g作为分岔参数, 采用复合胞坐标系方法分析其全局分岔特性, 探究在参数的连续变化下, 系统激变现象及域边界突变的演化过程.

      在一定的参数变化范围内, 系统呈现丰富的激变现象, 如周期或混沌吸引子与域边界上的周期鞍或混沌鞍发生碰撞的边界激变, 周期吸引子或混沌吸引子同时与边界上的混沌鞍发生碰撞的合并激变等. 当多吸引子共存时, 域边界会呈现分形结构, 并随着参数的变化, 发生Wada-Wada, Wada-分形和分形-Wada等域边界突变现象. 本文的研究结果对于通过控制温度和激励强度等参数, 调控SMA薄板系统的动态响应, 优化机械设备的变形及振动控制等问题上提供有效的分析手段.

参考文献 (25)

目录

    /

    返回文章
    返回