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

研究了含分数阶Caputo导数的非线性振动系统响应的数值计算方法。首先,由Caputo分数阶导数算子的叠加关系,得到含分数阶导数项非线性振动系统状态方程的标准形式。其次,基于Caputo导数与Riemann-Liouville导数和Grunwald-Letnikov导数间的关系,推导计算了Caputo导数的一般数值迭代格式。本文方法不要求状态方程中各分数阶导数阶数相等,弱化了已有算法中对分数阶导数阶数的限制,并可推广到多自由度的情形。随后,选择若干有解析解的算例验证了本文方法的正确性。最后,以多吸引子共存的分数阶Duffing振子系统为例,比较Caputo和GL两种算法所得结果,说明了用GL算法求解存在的问题。     

三种分形和分数阶导数阻尼振动模型的比较研究

标准的整数阶导数方程不能准确描述粘弹性材料的记忆性参考文献[1]和阻尼的分数次幂频率依赖[2],因此分形导数、分数阶导数及正定分数阶导数被用于描述粘弹性介质中的阻尼振动.该文通过分析模型和数值模拟,比较了三种模型描述的振动过程.结果显示,当p小于约O.75或大于约1.9时(p为非整数阶导数的阶数),分形导数模型衰减最快;当P大于约0.75且小于约1.9时,正定分数阶导数模型衰减最快,衰减最慢的分别为分数阶导数模型(p<1)和分形导数模型(p>1).且正定分数阶导数模型衰减快于分数阶导数模型,当p接近2时,两种模型较为相近.     

分数阶导数粘弹性模型的有限元法

在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数值解。并以二维沥青路面结构为例进行了路面动态粘弹性响应分析。算例分析表明,该方法能够正确有效地进行路面动态粘弹性分析。     

分数阶微分方程的Haar小波算法研究

在BPFs的Caputo分数阶微分算子矩阵的基础上,建立了Haar小波的分数阶微分算子矩阵,提出了一种有效的求解分数阶微分方程的Haar小波数值方法,并将该方法应用于线性和非线性分数阶常微分方程求解中.数值算例表明,该算法简单,数值精确度高,是一种高效的数值求解方法.     

具有分数导数本构关系的非线性粘弹性Timoshenko梁动力学行为分析

本文利用分数导数型本构关系建立了在有限变形情况下Timoshenko梁的控制方程并利用Galerkin方法进行简化。然后利用一种存储部分历史数据的分数积分的计算方法对梁的控制方程进行求解。考察了载荷参数和分数导数参数对梁振动的影响,并采用非线性动力学中的各种数值方法,如时程曲线、功率谱、相图、Poincare截面等,揭示了非线性粘弹性Timoshenko梁丰富的动力学行为。     

分数阶粘弹性地基上薄板波动和振动特性分析

本文研究了粘弹性地基上薄板的波动和振动问题．主要讨论了基于分数导数理论的粘弹性地基模型上 薄板弯曲波的传播特性以及固有频率对地基的依赖特性．推导了三种经典粘弹性地基模型的复模量．并利用分 数导数的性质得到分数阶粘弹性地基上 Kirchhoff板中弯曲波的传播速度、衰减系数以及自由振动的复固有频 率．数值算例表明粘弹性地基对弯曲波传播特性存在显著影响,不同粘弹性模型所对应的色散和衰减特性也存 在较大差别．分数阶导数可以实现相邻整数阶导数之间的光滑过渡．利用分数导数的本构关系可以更加真实地 描述粘弹性地基的历史依赖行为,更准确地表现出粘弹性地基板中弯曲波的色散和衰减特性．     

分数阶导数阻尼下非线性随机振动结构响应的功率谱密度估计

对于受到由分数阶导数模拟的粘弹性阻尼的非线性随机振动结构,本文给出了一种计算响应的功率谱密度方法。借助标准的随机平均法,首先得到了振动结构随机响应振幅的稳态概率密度。对于原振动结构的非线性项,运用改进的统计线性化方法得到了均方意义下的等价线性振动结构,并求得了其响应的依赖于振幅的条件功率谱密度。综合以上的结果,针对随机振动响应的功率谱密度的估计,通过与数值模拟结果进行验证,从而证明了所提方法的有效性和准确性。     

基于非局部理论的分数阶黏弹性场地地震放大效应分析

基于非局部理论和分数阶导数理论,研究上覆黏弹性场地土的地震放大效应。利用Eringen非局部理论考虑土体颗粒尺度等非局部效应的影响,通过分数阶黏弹性本构模型刻画场地土的应力应变本构关系,建立基于非局部理论的分数阶黏弹性场地土的振动微分方程;考虑分数阶导数的性质和黏弹性场地土的边界条件,得到了简谐地震波作用下黏弹性场地土的位移和剪切应力的解析解,并在频率域内给出了位移放大系数和应力放大系数的表达式;最后通过数值算例分析了非局部效应、分数阶导数的阶数和土体黏性参数等对黏弹性场地地震放大效应的影响。数值分析结果表明,在低频时位移放大系数和应力放大系数随频率变化曲线存在波动,高频时逐渐趋于稳定;非局部效应对场地土位移放大系数的影响与频率有关,对应力放大系数的影响较大,在研究场地土振动效应时有必要考虑土体非局部效应的影响;分数阶导数的阶数越小,位移放大系数和应力放大系数随频率变化曲线波动越大;场地土的力学性质对场地土的振动效应的影响较大;上覆场地土的黏性对位移放大系数的影响与频率有关,高频时,土体黏性越大,位移放大系数越大;越接近基岩,土体的应力放大系数越大,且土体深度对应力放大系数的影响越大。     

上覆分数导数粘弹性场地土地震放大效应

基于一维波动模型和分数导数粘弹性本构关系,分析了在竖直方向上传播的剪切地震波作用下,基岩上分数导数粘弹性模型描述的场地土的横向振动问题,用直接刚度矩阵法求得了场地土的地震放大效应系数,并用数值算例讨论了相关参量对分数导数粘弹性场地土地震放大效应系数的影响。研究结果表明:在简谐剪切地震波作用下,分数导数粘弹性场地土存在共振现象;分数导数的阶数、模型参数和基岩与上覆场地土层底部之间的阻抗比对场地土的地震放大效应系数有较大的影响。     

Caputo导数下分数阶Birkhoff系统的准对称性与分数阶Noether定理

应用分数阶模型可以更准确地描述和研究复杂系统的动力学行为和物理过程,同时Birkhoff力学是Hamilton力学的推广,因此研究分数阶Birkhoff系统动力学具有重要意义.分数阶Noether定理揭示了Noether对称变换与分数阶守恒量之间的内在联系,但是当变换拓展为Noether准对称变换时,该定理的推广遇到了很大的困难.本文基于时间重新参数化方法提出并研究Caputo导数下分数阶Birkhoff系统的Noether准对称性与守恒量.首先,将时间重新参数化方法应用于经典Birkhoff系统的Noether准对称性与守恒量研究,建立了相应的Noether定理;其次,基于分数阶Pfaff作用量分别在时间不变的和一般单参数无限小变换群下的不变性给出分数阶Birkhoff系统的Noether准对称变换的定义和判据,基于Frederico和Torres提出的分数阶守恒量定义,利用时间重新参数化方法建立了分数阶Birkhoff系统的Noether定理,从而揭示了分数阶Birkhoff系统的Noether准对称性与分数阶守恒量之间的内在联系.分数阶Birkhoff系统的Noether对称性定理和经典Birkhoff系统的Noether定理是其特例.最后以分数阶Hojman-Urrutia问题为例说明结果的应用.     

On the definition of fractional derivatives in rheology

During the last two decades fractional calculus has been increasingly applied to physics, especially to rheology. It is well known that there are obivious differences between Riemann-Liouville (R-L) definition and Caputo definition, which are the two most commonly used definitions of fractional derivatives. The multiple definitions of fractional derivatives have hindered the application of fractional calculus in rheology. In this paper, we clarify that the R-L definition and Caputo definition are both rheologically unreasonable with the help of the mechanical analogues of the fractional element model. We also find that to make them more reasonable rheologically, the lower terminals of both definitions should be put to ?∞. We further prove that the R-L definition with lower terminal ?∞ and the Caputo definition with lower terminal ?∞ are equivalent in the differentiation of functions that are smooth enough and functions that have finite number of singular points. Thus we can define the fractional derivatives in rheology as the R-L derivatives with lower terminal ?∞ (or, equivalently, the Caputo derivatives with lower terminal ?∞ ) not only for steady-state processes, but also for transient processes.     

Stable fractional Chebyshev differentiation matrix for the numerical solution of multi-order fractional differential equations

This paper presents an algorithm to obtain numerically stable differentiation matrices for approximating the left- and right-sided Caputo-fractional derivatives. The proposed differentiation matrices named fractional Chebyshev differentiation matrices are obtained using stable recurrence relations at the Chebyshev–Gauss–Lobatto points. These stable recurrence relations overcome previous limitations of the conventional methods such as the size of fractional differentiation matrices due to the exponential growth of round-off errors. Fractional Chebyshev collocation method as a framework for solving fractional differential equations with multi-order Caputo derivatives is also presented. The numerical stability of spectral methods for linear fractional-order differential equations (FDEs) is studied by using the proposed framework. Furthermore, the proposed fractional Chebyshev differentiation matrices obtain the fractional-order derivative of a function with spectral convergence. Therefore, they can be used in various spectral collocation methods to solve a system of linear or nonlinear multi-ordered FDEs. To illustrate the true advantages of the proposed fractional Chebyshev differentiation matrices, the numerical solutions of a linear FDE with a highly oscillatory solution, a stiff nonlinear FDE, and a fractional chaotic system are given. In the first, second, and forth examples, a comparison is made with the solution obtained by the proposed method and the one obtained by the Adams–Bashforth–Moulton method. It is shown the proposed fractional differentiation matrices are highly efficient in solving all the aforementioned examples.     

An optimization method based on the generalized polynomials for nonlinear variable-order time fractional diffusion-wave equation

In this paper, an optimization method based on a new class of basis functions, namely generalized polynomials (GPs), is proposed for nonlinear variable-order time fractional diffusion-wave equation. Variable-order time fractional derivative is expressed in the Caputo sense. In the proposed method, solution of the problem under consideration is expanded in terms of GPs with unknown free coefficients and control parameters. In this way, some new operational matrices of the ordinary and fractional derivatives are derived for these basis functions. The residual function and its 2-norm are employed for converting the problem under study to an optimization one and then choosing the unknown free coefficients and control parameters optimally. As a useful result, the necessary conditions of optimality are derived as a system of nonlinear algebraic equations with unknown free coefficients and control parameters. The validity and effectiveness of the method are demonstrated by solving some numerical examples. The results demonstrate that the proposed method is a powerful algorithm with good accuracy for solving such kind of problems.     

Analytical solution for the time-fractional heat conduction equation in spherical coordinate system by the method of variable separation

In this paper,using the fractional Fourier law,we obtain the fractional heat conduction equation with a time-fractional derivative in the spherical coordinate system.The method of variable separation is used to solve the timefractional heat conduction equation.The Caputo fractional derivative of the order 0 < α≤ 1 is used.The solution is presented in terms of the Mittag-Leffler functions.Numerical results are illustrated graphically for various values of fractional derivative.     

Discussion of the published manuscript lazopoulos K.A., karaoulanis D., lazopoulos Α.Κ. (2016) «on fractional modelling of viscoelastic mechanical systems», mechanics research communications, 78, pp. 1–5

The paper [20] attempts to solve a one dimensional model of the Zener type viscoelastic model (see Fig. 1 of the paper) replacing the classical time derivatives with a new fractional time derivative that was presented earlier by two of the authors (see ref. [1], [9] and [10] cited in the work). In addition, the authors present the classical solution of what is called in the paper “conventional viscoelastic model” and the solution due to Caputo fractional derivative, which is the solution of what the authors call “existing fractional order derivative viscoelastic models”. The paper seems to have several problems, typos and mistakes that will be discussed as follows     

On fallacies in the decision between the Caputo and Riemann–Liouville fractional derivatives for the analysis of the dynamic response of a nonlinear viscoelastic oscillator

Recently Dal [Dal, F., 2011. Multiple time scale solution of an equation with quadratic and cubic nonlinearities having fractional-order derivative. Mathematical and Computational Applications 16 (1), 301–308] presented ‘a new analytical scheme’ to calculate the dynamic response of a fractionally damped nonlinear oscillator possessing both quadratic and cubic nonlinearities via the method of multiple time scales. It has been claimed that damping features are modeled via the Caputo fractional derivative. In the present paper, it is shown that both the scheme and the object of investigation are not new, and moreover, the above mentioned author's statement is inconsistent, since under the assumptions made in the paper under consideration these two fractional-order derivatives coincide. Besides, the utilized procedure was inconsequential. It has been proved that the investigation of the dynamic response of a nonlinear viscoelastic oscillator presents the case that, with some minimal restrictions, the Riemann–Liouville and Caputo definitions produce completely equivalent mathematical models of the nonlinear viscoelastic phenomenon.     

Solution for a Fractional Diffusion-Wave Equation Defined in a Bounded Domain

A general solution is given for a fractional diffusion-wave equation defined in a bounded space domain. The fractional time derivative is described in the Caputo sense. The finite sine transform technique is used to convert a fractional differential equation from a space domain to a wavenumber domain. Laplace transform is used to reduce the resulting equation to an ordinary algebraic equation. Inverse Laplace and inverse finite sine transforms are used to obtain the desired solutions. The response expressions are written in terms of the Mittag–Leffler functions. For the first and the second derivative terms, these expressions reduce to the ordinary diffusion and wave solutions. Two examples are presented to show the application of the present technique. Results show that for fractional time derivatives of order 1/2 and 3/2, the system exhibits, respectively, slow diffusion and mixed diffusion-wave behaviors.     

Galerkin Projections and Finite Elements for Fractional Order Derivatives

Ordinary differential equations (ODEs) with fractional order derivatives are infinite dimensional systems and nonlocal in time: the history of the state variable is needed to calculate the instantaneous rate of change. This nonlocal nature leads to expensive long-time computations (O(t ²) computations for solution up to time t). A finite dimensional approximation of the fractional order derivative can alleviate this problem. We present one such approximation using a Galerkin projection. The original infinite dimensional system is replaced with an equivalent infinite dimensional system involving a partial differential equation (PDE). The Galerkin projection reduces the PDE to a finite system of ODEs. These ODEs can be solved cheaply (O(t) computations). The shape functions used for the Galerkin projection are important, and given attention. The approximation obtained is specific to the fractional order of the derivative; but can be used in any system with a derivative of that order. Calculations with both global shape functions as well as finite elements are presented. The discretization strategy is improved in a few steps until, finally, very good performance is obtained over a user-specifiable frequency range (not including zero). In particular, numerical examples are presented showing good performance for frequencies varying over more than 7 orders of magnitude. For any discretization held fixed, however, errors will be significant at sufficiently low or high frequencies. We discuss why such asymptotics may not significantly impact the engineering utility of the method.     

Steady state response analysis for fractional dynamic systems based on memory-free principle and harmonic balancing

A semi-analytic approach is proposed to analyze steady state responses of dynamic systems containing fractional derivatives. A major purpose is to efficiently combine the harmonic balancing (HB) technique and Yuan–Agrawal (YA) memory-free principle. As steady solutions being expressed by truncated Fourier series, a simple yet efficient way is suggested based on the YA principle to explicitly separate the Caputo fractional derivative as periodic and decaying non-periodic parts. Neglecting the decaying terms and applying HB procedures result into a set of algebraic equations in the Fourier coefficients. The linear algebraic equations are solved exactly for linear systems, and the non-linear ones are solved by Newton–Raphson plus arc-length continuation algorithm for non-linear problems. Both periodic and triple-periodic solutions obtained by the presented method are in excellent agreement with those by either predictor–corrector (PC) or YA method. Importantly, the presented method is capable of detecting both stable and unstable periodic solutions, whereas time-stepping integration techniques such as YA and PC can only track stable ones. Together with the Floquet theory, therefore, the presented method allows us to address the bifurcations in detail of the steady responses of fractional Duffing oscillator. Symmetry breakings and cyclic-fold bifurcations are found and discussed for both periodic and triple-periodic solutions.