分数阶Duffing振子的动力学研究-动力学与控制学报.PDF

第６卷第２期２００８年６月 动力 学 与控 制 学报 Ｖｏｌ．６Ｎｏ．２ １６７２６５５３／２００８／０６ ／１２２４ Ｊｕｎ．２００８ ⑵ ＪＯＵＲＮＡＬＯＦＤＹＮＡＭＩＣＳＡＮＤＣＯＮＴＲＯＬ 分数阶Ｄｕｆｆｉｎｇ振子的动力学研究 廖少锴　张卫 （暨南大学理工学院力学与土木工程系，暨南大学应用力学研究所，广州　５１０６３２） 摘要　在经典Ｄｕｆｆｉｎｇ振子中引入分数微分型阻尼项，推导了高效率的数值计算格式，对其表现出来的特有 的非线性现象进行讨论．研究表明：分数微分型阻尼的分数阶值较小时，振子将出现倍周期分岔并导致混 沌．在不同的外激励频率下，分数微分型Ｄｕｆｆｉｎｇ振子会呈现对称性破缺、分岔、混沌等强烈的非线性现象；在 一定参数范围内，分数微分型Ｄｕｆｆｉｎｇ振子较经典Ｄｕｆｆｉｎｇ振子，在较小的激励下即可进入混沌． 关键词　分数微分，　Ｄｕｆｆｉｎｇ振子，　Ｎｅｗｍａｒｋ法 ２ ｑ ３ Ｄｘ＋ｃＤｘ－ｘ＋ｘ＝Ｆｃｏｓ（ｔ） ω 引言 {０＜ｑ＜１ （２） 经典Ｄｕｆｆｉｎｇ振子可看作是含粘性阻尼的单摆 上式分数微分为Ｒｉｅｍａｎｎ－Ｌｉｏｕｖｉｌｅ定义的分数微 作大振幅振荡时非线性振动的描述，其阻尼项是基 分： 于整数阶微积分描述的牛顿粘性阻尼．近年来越来 ｔ １ ｄ ｆ（） ｑ τ Ｄｆ（ｔ）＝ ｄ， τ ∫ ｑ ０ Γ（１－ｑ）ｄｔ （ｔ－ ） 越多的研究表明：基于整数阶微积分算子描述的经 τ 典的粘弹性本构关系对于相当一大类高分子阻尼 　０＜ｑ＜１ （３） 材料并不能很好地描述其力学行为．越来越多的实 考虑到分数微分型Ｄｕｆｆｉｎｇ振子方程解析解很 验研究都表明：分数微积分型本构关系能在相当宽 难得到，本文将采用无条件收敛的Ｎｅｗｍａｒｋ法构 的频率域范围内仅用少量参数就可以对相当一大 造其数值求解单步算法，由文献［２－５］可知，分数 类高分子阻尼材料的力学行为进行很好地描 微分的Ｎｅｗｍａｒｋ型的Ｚｈａｎｇ－Ｓｈｉｍｉｚｕ数值算法为： 述［１－２］．