杜芬振子分岔及通往混沌道路的分析.doc

基于增量谐波平衡法的马休-杜芬振子分岔及通往混沌道路的分析 陈树辉， 沈建和 (中山大学应用力学与工程系， 广州， 510275) 摘要 本文提出一种分析Mathieu-Duffing振子动力学行为: 包括稳定性、分岔及通往混沌道路的方法。在该方法中，基于增量谐波平衡法, 求得特定参数状态下的周期解; 基于Floquet理论，考察求得的周期解的稳定性. 以系统的控制参数作为主动增量, 在每一步增量过程中，计算出周期解及周期解的Floquet乘子. 以Floquet乘子随着控制参数变化的演化作为判据,了解周期解的分岔类型及分岔值. 若周期解出现周期倍化分岔, 则通过引入适当的谐波项, 把分岔点处的周期解写成倍周期的形式, 然后把此时得到的形式解作为下一级周期倍化轨道的初始值. 再次利用增量谐波平衡法, 计算出周期倍化轨道及其Floquet乘子. 进一步追踪Floquet乘子随控制参数的变化的演化, 得到下一周期倍化分岔点. 重复上述过程, 则可半解析地证实马休-杜芬振子经连续的周期倍化分岔进入混沌的道路. 同时, 近似地获到周期倍化分岔的系列分岔点及混沌产生的阈值. 给出的分岔图及不同参数状态下的相图也定性地说明了该振子经连续的周期倍化分岔进入混沌的道路。 关键词 分岔；混沌产生的道路；Mathieu-Duffing振子；增量谐波平衡法；Floquet理论 引 言对非线性系统的动力学行为稳定性、分岔及通往混沌道路的研究引人们的兴趣这不仅动力系统在特定参数条件下的动力学行为了解从有序向转化的机制也是分岔控制及混沌控制的需要至今已有大量的文 E-mail: stscsh@mail.sysu.edu.cn 献动力系统经的分岔进入混沌的问题发现几种混沌典型途径. 但是,所采用的研究方法基本上集中于数值研究与实验研究或者两者的结合, 解析即使是半解析的研究手段也很少见. 本文基于增量谐波平衡法, 结合Floquet理论, 半解析地考察马休-杜芬振子经多次的分岔进入混沌的道路问题. 增量谐波平衡法(IHB法)由Lau为处理弹性系统的结构周期振动而提出. 之后, 该方法被广泛地运用于研究动力系统的多种动力学现象, 并由E. H. Dowell首次推荐用于求解Lorenz系统的混沌运动. 之后, Leung应用IHB法, 构造了Duffing振子的分岔边界及混沌区域. Ge等利用IHB法, 求得具有参数激励的单轴速率陀螺的周期解; 并基于多变量Floquet理论, 考察求得的周期解的稳定性; 最后, 利用数值积分, 计算周期解各类分岔的分岔点以及系统混沌产生的阈值. Xu等基于IHB法求得具有分段刚度与分段阻尼的单自由度谐激振子的周期解, 并通过不同周期响应的振幅峰值与控制参数的关系, 获得系统的周期倍化分岔点. A. Raghothama等也利用IHB法, 分别研究了转子轴承系统、具有分段非线性刚度的连接载荷平台模型及逃逸模型的周期解、周期解的分岔以及系统混沌运动的形成问题。但是，纵观这些文献，有几个不足之处：第一，这些文献都没有给出连续求解周期倍化轨道的定量方法, 因此, 对于高周期的轨道, 一般难以求得; 第二, 分岔点及混沌产生的阈值的确定依靠数值技巧, 例如are-length法等. 因此存在几个不足: 例如过程较为复杂, 精确度较低, 只能获得系统的前几个分岔点, 没有办法求得系统的不稳定解等; 第三，给出的计算周期轨道monodromy矩阵的公式有误。 本文基于IHB法与Floquet理论的结合，研究Mathieu-Duffing振子经连续的分岔进入混沌的道路问题。以IHB法定量地求出周期轨道及高周期轨道的表达式；以Floquet理论考察周期轨道与高周期轨道的稳定性；以周期轨道的Floquet乘子随系统控制参数变化的演化作为判据，定性地分析周期轨道的分岔类型及对应的分岔值。通过以上方法，可得到一系列的周期倍化分岔点。这些分岔点之间的距离越来越小，趋于某个极限值，该极限值代表系统混沌产生的阈值；同时发现，这一系列周期倍化分岔点近似地服从Feigenbaum于Logistic模型中所发现的普适规律. 1 增量谐波平衡计算公式 考虑如下Mathieu-Duffing方程: ， (1) 方程(1)可用于描述带有变速度的轴向运动梁的单模横向振动等力学与工程问题。引入新的无量纲时间尺度: , (2) 则方程(1)变为：(3) 设x0(?)为对应于初始振幅?0、频率 ?0 及阻尼?0的周期解, 则x, ?, ?及? 可写成 , (4) 这里为增量符号. 把(4)式代入(3)式展开,由于等带有增量符号的量是小量, 故可略去增量的高阶