基于Caputo导数的分数阶非线性振动系统响应计算  

Calculation of the response of nonlinear vibration system based on the Caputo fractional derivative

作　　者：李亚杰 吴志强[1,2] 章国齐 LI Ya-jie;WU Zhi-qiang;ZHANG Guo-qi(Department of Mechanics,School of Mechanical Engineering,Tianjin University,Tianjin 300072,China;Tianjin Key Laboratory of Nonlinear Dynamics and Chaos Control,Tianjin University,Tianjin 300072,China)

机构地区：[1]天津大学机械工程学院力学系,天津300072 [2]天津大学天津市非线性动力学与混沌控制重点实验室,天津300072

出　　处：《计算力学学报》2018年第4期466-472,共7页Chinese Journal of Computational Mechanics

基　　金：国家自然科学基金(11672349;11372211);国家973计划(2014CB046805)资助项目

摘　　要：研究了含分数阶Caputo导数的非线性振动系统响应的数值计算方法。首先,由Caputo分数阶导数算子的叠加关系,得到含分数阶导数项非线性振动系统状态方程的标准形式。其次,基于Caputo导数与Riemann-Liouville导数和Grunwald-Letnikov导数间的关系,推导计算了Caputo导数的一般数值迭代格式。本文方法不要求状态方程中各分数阶导数阶数相等,弱化了已有算法中对分数阶导数阶数的限制,并可推广到多自由度的情形。随后,选择若干有解析解的算例验证了本文方法的正确性。最后,以多吸引子共存的分数阶Duffing振子系统为例,比较Caputo和GL两种算法所得结果,说明了用GL算法求解存在的问题。In this paper,the method of calculating the response of a nonlinear vibration system with fractional order Caputo derivative is studied.Firstly,by the superposition relation of the fractional operator of Caputo derivative,we obtain the standard form of the state equation of nonlinear vibration system with fractional derivative.Secondly,the general numerical iterative scheme of Caputo derivative is derived based on the relationship between Caputo derivative and Riemann-Liouville derivative and Grunwald-Letnikov derivative.This method does not requires the orders of the various fractional order derivatives in the state equations to be equal,which not only can reduce the restriction on the orders of the fractional derivatives in existing algorithms,but also can be extended to multi-degree-of-freedom systems.Then,we select some examples which have the analytical solutions to verify the correctness of the algorithm presented in this paper.Finally,a fractional order Duffing oscillator with multiple basins of attraction is taken as an example,and then comparison of the results obtained by the Caputo and GL algorithms respectively shows the existing problems to solve the Caputo derivative with GL algorithm.

关 键 词：分数阶导数 算法 数值迭代格式 导数阶数 系统初值 

分 类 号：O322[理学—一般力学与力学基础]