Mathematical Modeling of Linear Fractional Oscillators

Abstract

In this work, based on Newton’s second law, taking into account heredity, an equation is derived for a linear hereditary oscillator (LHO). Then, by choosing a power-law memory function, the transition to a model equation with Gerasimov–Caputo fractional derivatives is carried out. For the resulting model equation, local initial conditions are set (the Cauchy problem). Numerical methods for solving the Cauchy problem using an explicit non-local finite-difference scheme (ENFDS) and the Adams–Bashforth–Moulton (ABM) method are considered. An analysis of the errors of the methods is carried out on specific test examples. It is shown that the ABM method is more accurate and converges faster to an exact solution than the ENFDS method. Forced oscillations of linear fractional oscillators (LFO) are investigated. Using the ABM method, the amplitude–frequency characteristics (AFC) were constructed, which were compared with the AFC obtained by the analytical formula. The Q-factor of the LFO is investigated. It is shown that the orders of fractional derivatives are responsible for the intensity of energy dissipation in fractional vibrational systems. Specific mathematical models of LFOs are considered: a fractional analogue of the harmonic oscillator, fractional oscillators of Mathieu and Airy. Oscillograms and phase trajectories were constructed using the ABM method for various values of the parameters included in the model equation. The interpretation of the simulation results is carried out.