Next Article in Journal
Topology Optimization of Elastoplastic Behavior Conditions by Selectively Suppressing Plastic Work
Next Article in Special Issue
Hereditary Oscillator Associated with the Model of a Large-Scale αω-Dynamo
Previous Article in Journal
On the Affine Image of a Rational Surface of Revolution
Previous Article in Special Issue
Mathematical Modeling of Linear Fractional Oscillators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Model of Fractional Duffing Oscillator with Variable Memory

by
Valentine Kim
1,2,† and
Roman Parovik
1,2,3,*,†
1
Department of Mathematics and Physics, Vitus Bering Kamchatka State University, Pogranichnaya, 4, 683032 Petropavlovsk-Kamchatskiy City, Russia
2
Department of Control Systems, Kamchatka State Technical University, Kluchevskaya, 35, 683003 Petropavlovsk-Kamchatskiy City, Russia
3
Institute of Cosmophysical Research and Radio Wave Propagation, Far East Branch, Russian Academy of Sciences, Mirnaya, 7, 684034 Paratunka, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2020, 8(11), 2063; https://doi.org/10.3390/math8112063
Submission received: 30 October 2020 / Revised: 15 November 2020 / Accepted: 16 November 2020 / Published: 19 November 2020
(This article belongs to the Special Issue Mathematical Modeling of Hereditarity Oscillatory Systems)

Abstract

:
The article investigates a mathematical model of the Duffing oscillator with a variable fractional order derivative of the Riemann–Liouville type. The study of the model is carried out using a numerical scheme based on the approximation of the fractional derivative of the Riemann–Liouville type by a discrete analog—the fractional derivative of Grunwald–Letnikov. The adequacy of the numerical scheme is verified using specific examples. Using a numerical algorithm, oscillograms and phase trajectories are constructed depending on the values of the model parameters. Chaotic regimes of the Duffing fractional oscillator are investigated using the Wolf–Bennetin algorithm. The forced oscillations of the Duffing fractional oscillator are investigated using the harmonic balance method. Analytical formulas for the amplitude-frequency, phase-frequency characteristics, and also the quality factor are obtained. It is shown that the fractional Duffing oscillator possesses different modes: regular, chaotic, multi-periodic. The relationship between the order of the fractional derivative and the quality factor of the oscillatory system is established.

1. Introduction

At present, one of the scientific directions has received wide development—hereditary dynamics, the founder of which is the Italian mathematician Vito Volterra [1]. Hereditarity is a property of a dynamic system in which its current state depends on previous states. Hereditary dynamical systems find their application in hereditary mechanics [2], biology [3], economics [4] and other fields of knowledge.
From the point of view of mathematics, a hereditary dynamical system can be described using integro-differential equations with difference kernels, which are called memory functions [1]. If the memory functions are chosen as exponential, then we turn to equations with derivatives of fractional derivatives, which are studied in the framework of fractional calculus [5,6,7]. Fractional derivatives have many definitions and unique properties, but all of them to some extent describe the memory effect that characterizes information about previous states of the system. This effect determines additional degrees of freedom—orders of fractional derivatives. The orders of fractional derivatives can be constant or functions of time, which in the latter case significantly complicates the study of a dynamical system.
It should be noted that in article [8] the Duffing oscillator with a constant order Riemann–Liouville derivative was considered. Chaotic and regular regimes are investigated using a numerical scheme and various tests. Similar studies were carried out in article [9], but with the Caputo derivative of fractional constant order. In article [10], on the basis of Melnikov’s method, the necessary condition for a chaotic regime of a Duffing oscillator with a fractional order derivative under harmonic excitation is investigated. In article [11], forced oscillations of the Duffing oscillator with a constant-order Caputo derivative were investigated. The amplitude-frequency, phase-frequency characteristics, as well as the bistable mode are investigated.
Article [12] proposes an analytical method for solving the Duffing equation with fractional Caputo derivatives based on a combination of the perturbation and Laplace methods. Paper [13] uses He’s fractional derivative to describe the fractional Duffing equation, which is determined using a variational iterative algorithm. To solve such an equation, a complex fractional transformation was used. In article [14], a solution of the Duffing oscillator with fractional derivatives was obtained using the method of homotopy analysis.
This article proposes a generalization of the Duffing oscillator, a nonlinear oscillatory system with chaotic dynamics. This generalization is based on the introduction of a variable order derivative of the Riemann–Liouville type to describe the energy dissipation of an oscillatory system. This article analyzes the chaotic and regular modes, frequency characteristics of the Duffing oscillator with variable memory, which is a generalization of the research from article [8].
The work is organized as follows. Section 2 gives some definitions from the theory of fractional disappearance; Section 3 presents a statement of the Cauchy problem for the Duffing equation with a fractional variable-order derivative of the Riemann–Liouville type in the dissipative term. It also shows a numerical scheme obtained on the basis of the discrete fractional Grunwald–Letnikov derivative. It is shown on a specific example that the circuit has the first order of accuracy. Section 4 contains an algorithm for calculating the maximum Lyapunov exponents. It is shown with a specific example that chaotic regimes are possible under certain conditions. In Section 5, we studied forced oscillations of the Duffing fractional oscillator. Using the harmonic balance method, formulas are derived for calculating the amplitude-frequency and phase-frequency characteristics, as well as the quality factor. Surfaces of the corresponding characteristics are constructed. The conclusion is given in Section 6.

2. Preliminary Material on Fractional Calculus

Definition 1.
The fractional derivative of variable order 0 < q t < 1 of the Riemann–Liouville type has the form:
D 0 t q t x τ = 1 Γ 1 q t d d t 0 t x τ d τ t τ q t , t > 0 ,
where Γ 1 q t is Euler’s gamma function.
Definition 2.
The fractional derivative of variable order 0 < q t < 1 of the Grunwald–Letnikov type has the form [15]:
D 0 t q t x τ = lim h 0 , k h = t τ q t j = 0 k 1 1 j q t j x t j h ,
where q t j —binomial coefficient.
Remark 1.
The weighting factor in Formula (2) can be determined as follows:
c j = 1 j q t j = 1 q t + 1 j c j 1 , j = 1 , 2 , 3 , , c 0 = 1 .

3. Problem Statement and Numerical Solution

Consider the following Cauchy problem with offset [16]:
x ¨ t + λ D 0 t q t x τ + x t + x 3 t = δ cos ω t , x 0 = x 0 , x ˙ 0 = y 0 ,
where x t —displacement, solution function, x 0 and y 0 —given constants, λ is the coefficient of friction, δ and ω are the amplitude and frequency of the external periodic impact, 0 < q t < 1 , t 0 , T , T > 0 —simulation time.
Remark 2.
The variable order q t of the fractional derivative determines the intensity of energy dissipation in the oscillatory system. If this order is constant and equal to one, then the Cauchy problem (3) becomes the Cauchy problem for the classical Duffing oscillator [4].
Let x t C 3 0 , T to achieve the required smoothness. Divide the time interval 0 , T into N equal parts with a constant step h = T T N N . The solution function x t will go to the grid function x t k = x k , where t k + 1 = k + 1 h , k = 1 , , N 1 .
Remark 3.
The following relation is valid
D 0 t q k x τ = h q k j = 0 k 1 c j x k j + O h ,
The Cauchy problem (3) may be represented as a system:
x ˙ t = y t , D 0 t q ( t ) x τ = w t , y ˙ t = δ cos ω t λ w t x 3 t x t , x 0 = x 0 , x ˙ 0 = y 0 ,
Let us rewrite the Cauchy problem (5) in the difference setting taking into account the relation (4)
x k = h y k 1 + x k 1 , w k = h q k i = 0 k 1 c i x k i , y k = y k 1 + h ( δ cos ( ω h ( k 1 ) ) x k 1 x k 1 3 λ w k ) ,
Remark 4.
Note that relation (6) is an explicit Euler approximation, which is conditionally stable, and taking into account relation (4), the scheme has the first order of accuracy.
Let us show the last remark with a specific example.
Example 1.
Consider the case:
q ( t ) = sin 2 ( ω t ) , N = 2500 , T = 500 , λ = 0.05 , δ = 0.12 , ω = 0.075 , x 0 = 0.3 , y 0 = 0.5 .
Consider the change of the absolute error ε and computational accuracy α = ln ε / ln h / 2 of scheme (6) when changing the step h. To calculate the absolute error, we use the Runge rule: ε = max i x i x 2 i —absolute magnitude error, x 2 i is the numerical solution at the step h h 2 2 .
From Table 1, we see that with an increase in the number of nodes in the computational grid, the error e decreases, while the computational accuracy tends to unity.
Figure 1 shows the oscillogram (Figure 1a) and the phase trajectory (Figure 1b).

4. Lyapunov Exponents

In the study of nonlinear systems, one of the important tasks is to determine the type of oscillations-periodic, quasi-periodic, random, chaotic. It is especially difficult to distinguish quasiperiodic oscillations from chaotic and random ones, since quasiperiodic oscillations often have a very complex shape that is visually indistinguishable from “random” [17].
A feature of chaotic oscillations is their high sensitivity to small changes in initial conditions. Therefore, one of the most reliable ways to detect chaos is to determine the speed of trajectory run-up, which is estimated using a spectrum of Lyapunov exponents [18].
The geometric meaning of the Lyapunov exponents is that two solutions whose initial values are located in a certain neighborhood of radius ϵ will diverge in time T into an n-dimensional ellipsoid along n main semi-axes, and at time t the radii will be determined by the values
ε e Λ i t , i = 1 , , n ,
where Λ i —Lyapunov exponents.
If two trajectories remain close over time, then their change obeys the law
u j + 1 = J u j u j , J u j = f 1 x j f 1 y j f 2 x j f 2 y j .
where J is Jacobi matrix, u j = ( x j , y j ) T is column vector, f 1 , f 2 —right sides of the system (5). The Formula (7) defines the perturbed equations.
Remark 5.
The presence of at least one positive Lyapunov exponent in the spectrum means the presence of a chaotic regime (asymptotic instability) of the phase trajectory under consideration. A negative reading indicates a limit cycle, while a zero reading indicates that analysis is not enough to draw a conclusion about stability and instability according to Lyapunov. The vectors of the initial conditions for Equation (7) must be orthogonal.
Consider the Benettin–Wolf algorithm for calculating Lyapunov exponents and implement it in the MAPLE environment [19].
  • Select the starting point of the vector x 0 and together with it we will track the K disturbed trajectories. Let K = 3 .
  • We solve the original equation numerically together with three sets of perturbed equations or equations in variations (7). As the initial vectors for equations in variations, you must select a set of vectors x ˜ 0 0 , y ˜ 0 0 , z ˜ 0 0 , that are orthogonal and normalized by one.
  • After time T, the trajectory will move to the point vector x 1 , the perturbation vectors x ˜ 1 , y ˜ 1 , z ˜ 1 , which we renormalize using the Gram–Schmidt method.
    x ˜ 1 0 = x ˜ 1 x ˜ 1 ,
    y ˜ 1 = y ˜ 1 y ˜ 1 , x ˜ 1 0 x ˜ 1 0 , y ˜ 1 0 = y ˜ 1 y ˜ 1 ,
    z ˜ 1 = z ˜ 1 z ˜ 1 , x ˜ 1 0 x ˜ 1 0 z ˜ 1 , y ˜ 1 0 y ˜ 1 0 , z ˜ 1 0 = z ˜ 1 z ˜ 1 ,
  • Next, we continue counting from the point x 1 and the perturbation vectors x ˜ 1 0 , y ˜ 1 0 , z ˜ 1 0 . After a regular time interval T, we get a new set of perturbation vectors x ˜ 2 , y ˜ 2 , z ˜ 2 , which is subjected to orthogonalization and renormalization.
  • Steps 2–4 are repeated M times, and the sums of the calculations are calculated: S 1 = i = 1 ln x ˜ i , S 2 = i = 1 ln y ˜ i , S 3 = i = 1 ln z ˜ i , which include perturbation vectors before renormalization, but after orthogonalization.
  • The estimation of the first three Lyapunov exponents is calculated using the formula:
    Λ i = S i M T , i = 1 , 2 , 3 .
For a more simplified calculation, take K = 2 . As a result, we get the following equations:
x k = h y k 1 + x k 1 , y k = y k 1 + h ( x k 1 x k 1 3 + δ cos ( ω h ( k 1 ) ) λ h q k i = 0 k 1 c i x k i ) ,
x x k = h y y k 1 + x x k 1 , y y k = y y k 1 + h ( ( 1 3 x k 1 2 ) x x k 1 λ h q k i = 0 k 1 c i x x k i ) , x x 0 = 1 , y y 0 = 0 .
X X k = h Y Y k 1 + X X k 1 , Y Y k = Y Y k 1 + h ( ( 1 3 x k 1 2 ) X X k 1 λ h q k i = 0 k 1 c i X X k i ) , X X 0 = 0 , Y Y 0 = 1 .
where x x k , y y k , X X k , Y Y k solutions of perturbed Equations (9) and (10) according to Formula (7).
Applying the Benettin–Wolf algorithm on Equations (8)–(10), namely steps 3–6, we get the values of the senior Lyapunov exponent. Let us consider the work of the Bennetin–Wolf algorithm using specific examples.
Example 2.
Task parameter values: N = 2000 , h = 0.05 , δ = 0.3 , ω = 1 , q t = 1 / exp ( cos ( ω t ) ) , x 0 = 0.2 , y 0 = 0.3 .
Figure 2 shows the spectrum of Lyapunov exponents relative to the coefficient of dissipation λ (Figure 2a) and phase trajectories at different values λ (Figure 2b–d). The figures show that for positive values of Lyapunov exponents, chaotic modes are observed, and the lower the value λ , the more closely the trajectory approaches that of the classical Duffing oscillator. For negative values of Lyapunov exponents, the phase trajectories are limit cycles, and the greater the modulus of the negative Lyapunov exponent, the more stable the limit cycle will be.
Example 3.
Task parameter values: N = 2000 , h = 0.025 , δ = 0.3 , ω = 1 , λ = 0.15 , x 0 = 0.2 , y 0 = 0.3 .
Figure 3 shows the spectrum of Lyapunov exponents relative to the fractional derivative q (Figure 3a) and phase trajectories at different values of q (Figure 3b,c). The figures show that for positive values of Lyapunov exponents, chaotic modes are observed, and the closer the value of q is to one, the more closely the trajectory approaches that of the classical Duffing oscillator. With these values of control parameters, the spectrum does not show negative values.

5. Forced Oscillations

As shown by the results of articles [20,21], the orders of fractional derivatives in oscillatory systems with memory are associated with the quality factor and are responsible for the intensity of energy dissipation. Therefore, let us investigate the forced oscillations of the Duffing fractional oscillator (3).
Consider the dependence of the amplitude A and phase φ of steady-state oscillations for operator (1) on the frequency of the external force. To calculate the frequency response and frequency response, we use the harmonic linearization method [21,22,23].
Consider the Cauchy problem (3). We look for a solution in the form:
x t = A cos ( ω t + φ ) = A cos ( U ) .
D 0 t q ( t ) x ( τ ) = 1 Γ ( 1 q ( t ) ) d d t 0 t x ( τ ) d τ ( t τ ) q ( t ) = v = t τ d v = d τ = 1 Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) x ( t v ) d v = A Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) cos ( ω ( t v ) + φ ) d v
= A Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) [ cos ( ω t ω v ) cos ( φ ) sin ( ω t ω v ) sin ( φ ) ] d v = A cos ( φ ) Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) cos ( ω t ω v ) d v + A sin ( φ ) Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) sin ( ω t ω v ) d v
= A cos ( φ ) Γ ( 1 q ( t ) ) d d t cos ( ω t ) 0 t v q ( t ) cos ( ω v ) d v A cos ( φ ) Γ ( 1 q ( t ) ) d d t sin ( ω t ) 0 t v q ( t ) sin ( ω v ) d v + A sin ( φ ) Γ ( 1 q ( t ) ) d d t sin ( ω t ) 0 t v q ( t ) cos ( ω v ) d v A sin ( φ ) Γ ( 1 q ( t ) ) d d t cos ( ω t ) 0 t v q ( t ) sin ( ω v ) d v
= A cos ( φ ) Γ ( 1 q ( t ) ) ω sin ( ω t ) 0 t v q ( t ) cos ( ω v ) d v + cos ( ω t ) d d t 0 t v q ( t ) cos ( ω v ) d v A cos ( φ ) Γ ( 1 q ( t ) ) ω cos ( ω t ) 0 t v q ( t ) sin ( ω v ) d v + sin ( ω t ) d d t 0 t v q ( t ) sin ( ω v ) d v
+ A sin ( φ ) Γ ( 1 q ( t ) ) ω cos ( ω t ) 0 t v q ( t ) cos ( ω v ) d v + sin ( ω t ) d d t 0 t v q ( t ) cos ( ω v ) d v A sin ( φ ) Γ ( 1 q ( t ) ) ω sin ( ω t ) 0 t v q ( t ) sin ( ω v ) d v + cos ( ω t ) d d t 0 t v q ( t ) sin ( ω v ) d v
= A ω cos ( φ ) sin ( ω t ) Γ ( 1 q ( t ) ) 0 t v q ( t ) cos ( ω v ) d v A cos ( φ ) cos ( ω t ) Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) cos ( ω v ) d v A ω cos ( φ ) cos ( ω t ) Γ ( 1 q ( t ) ) 0 t v q ( t ) sin ( ω v ) d v A cos ( φ ) sin ( ω t ) Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) sin ( ω v ) d v
+ A ω sin ( φ ) cos ( ω t ) Γ ( 1 q ( t ) ) 0 t v q ( t ) cos ( ω v ) d v + A sin ( φ ) sin ( ω t ) Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) cos ( ω v ) d v ) + A ω sin ( φ ) sin ( ω t ) Γ ( 1 q ( t ) ) 0 t v q ( t ) sin ( ω v ) d v A sin ( φ ) cos ( ω t ) Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) sin ( ω v ) d v )
We group expressions by integrals from sine and cosine and differentials from integrals from sine and cosine [21].
= A ω Γ ( 1 q ( t ) ) [ cos ( φ ) sin ( ω t ) + sin ( φ ) cos ( ω t ) ] 0 t v q ( t ) cos ( ω v ) d v + A ω Γ ( 1 q ( t ) ) [ cos ( φ ) cos ( ω t ) + sin ( φ ) sin ( ω t ) ] 0 t v q ( t ) sin ( ω v ) d v
+ A Γ ( 1 q ( t ) ) [ cos ( φ ) cos ( ω t ) + sin ( φ ) sin ( ω t ) ] d d t 0 t v q ( t ) cos ( ω v ) d v + A Γ ( 1 q ( t ) ) [ cos ( φ ) sin ( ω t ) sin ( φ ) cos ( ω t ) ] d d t 0 t v q ( t ) sin ( ω v ) d v
= A ω sin ( ω t + φ ) Γ ( 1 q ( t ) ) 0 t v q ( t ) cos ( ω v ) d v A ω cos ( ω t + φ ) Γ ( 1 q ( t ) ) 0 t v q ( t ) sin ( ω v ) d v A cos ( ω t + φ ) Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) cos ( ω v ) d v A sin ( ω t + φ ) Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) sin ( ω v ) d v
Thus, the Riemann–Liouville operator for solution (11) has the form:
D 0 t q ( t ) x ( τ ) = 1 Γ ( 1 q ( t ) ) d d t 0 t x ( τ ) d τ ( t τ ) q ( t ) = v = t τ d v = d τ = A ω sin ( ω t + φ ) Γ ( 1 q ( t ) ) 0 t v q ( t ) cos ( ω v ) d v A ω cos ( ω t + φ ) Γ ( 1 q ( t ) ) 0 t v q ( t ) sin ( ω v ) d v A cos ( ω t + φ ) Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) cos ( ω v ) d v A sin ( ω t + φ ) Γ ( 1 q ( t ) ) d d t 0 t v q ( t ) sin ( ω v ) d v
In the case of steady oscillations at large times t ( t ), the integrals (12) can be written in the form [23]:
I sin = 0 v q ( t ) sin ( ω v ) d v = Γ ( 1 q ( t ) ) ω 1 q ( t ) cos q ( t ) π 2
I cos = 0 v q ( t ) cos ( ω v ) d v = Γ ( 1 q ( t ) ) ω 1 q ( t ) sin q ( t ) π 2
Substituting (13) and (14) in (12) we’ll get
D 0 t q ( t ) x ( τ ) = A ω q ( t ) sin ( ω t + φ ) sin q ( t ) π 2 A ω q ( t ) cos ( ω t + φ ) cos q ( t ) π 2 A d q d t ω q ( t ) 1 cos ( ω t + φ ) ( ln ( ω ) ψ ( 1 q ( t ) ) ) sin q ( t ) π 2 + π cos q ( t ) π 2 2 A d q d t ω q ( t ) 1 sin ( ω t + φ ) ( ln ( ω ) ψ ( 1 q ( t ) ) ) cos q ( t ) π 2 π sin q ( t ) π 2 2
Substituting (11) and (15) in Equation (3), given that:
x ( t ) = A cos ( ω t + φ ) , x ( t ) = A ω sin ( ω t + φ ) , x ( t ) = A ω 2 cos ( ω t + φ ) ,
cos 3 ( ω t + φ ) = 3 4 cos ( ω t + φ ) + 1 4 cos ( 3 ω t + 3 φ ) ,
We get the equation [23]:
x ¨ t + p x ˙ t + s 2 x t = δ cos ω t , x 0 = x 0 , x ˙ 0 = y 0 ,
where
p = 2 λ ω q ( t ) 1 sin q ( t ) π 2 + 2 λ d q d t ω q ( t ) 2 ( ln ( ω ) ψ ( 1 q ( t ) ) ) cos q ( t ) π 2 π sin q ( t ) π 2 2 ,
s 2 = ω 0 2 2 λ ω q ( t ) cos q ( t ) π 2 + 3 A 2 b 4 2 λ d q d t ω q ( t ) 1 ( ln ( ω ) ψ ( 1 q ( t ) ) ) sin q ( t ) π 2 + π cos q ( t ) π 2 2 .
There ψ ( 1 q ( t ) ) = γ + n = 1 1 n 1 n + 1 q ( t ) —digamma function. γ 1 2 ( 10 3 1 ) —the constant of Euler. In practice it is enough to take n = 50 .
Equation (16) is a classical linear oscillator for which the formulas for the Amplitude-frequency and Phase-frequency characteristics (AFC and PFC) are derived using the harmonic linearization method:
A = δ s 2 ω 2 2 + p 2 ω 2
φ = arctan p ω s 2 ω 2 .
The Q-factor for a fractional Duffing oscillator can be determined from Equation (16) using the formula [21,23]:
Q = s p ,
where s and p are defined by the Equalities (17) and (18).
Example 4.
Task parameter values: N = 3000 , T = 100 , λ = 0.1 , δ = 1 , ω = 1.5 , x 0 = 0.3 , y 0 = 0.5 .
Figure 4 shows that when the q parameter is reduced, the quality-factor increases (21). The maximum amplitude corresponds to the maximum Q-factor, and when the frequency decreases, the Q-factor decreases. However, most of all, Q-factor depends on the parameter q ( t ) .
Figure 5 shows AFC and PFC versus amplitude and phase versus frequency and the exponent of the derivative q ( t ) . In Figure 5a, we see separate areas from low frequencies to high frequencies.

6. Conclusions

The article proposes a mathematical model of a fractional Duffing oscillator with a variable fractional order derivative of the Riemann–Liouville type in the dissipative term. Using the discrete Grunwald–Letnikov derivative, a first-order numerical scheme is constructed. Using Runge’s method, it is shown that with an increase in the number of nodes in the computational grid, the order of computational accuracy tends to unity. A numerical algorithm is proposed for calculating the spectra of the maximum Lyapunov exponents. It is shown that the spectrum of the maximum Lyapunov exponents can contain positive values, which indicates the existence of chaotic regimes. Knowing the conditions for the existence of chaotic regimes, we can investigate regular regimes, for example, regimes of steady oscillations under the action of an external harmonic force. The study of forced vibrations using analytical formulas for the frequency characteristics and quality factor of the oscillatory system showed that the variable order of the fractional derivative determines the intensity of energy dissipation in the system. With an increase in the order, the quality factor of the oscillatory system increases, and with a decrease, it increases. This fact suggests that oscillatory modes can be controlled using the order of the fractional derivative. This fact is very important for the development of applied engineering research. For example, in machine learning and intelligent communications [24].
It is necessary to study the fractional Duffing oscillator with a variable fractional derivative of the Riemann–Liouville type, taking into account different definitions and types of its function [15]. It is of interest to study forced oscillations, resonance curves as in the Figure 5a. Furthermore, also the presence of bistable behavior [11].
It is of interest to study the fractional Duffing oscillator with higher order nonlinearities, for example, by analogy with the article [25].

Author Contributions

Conceptualization, R.P.; methodology, R.P.; software, V.K.; validation, V.K.; formal analysis, V.K.; investigation, V.K.; writing—original draft preparation, V.K.; writing—review and editing, R.P.; visualization, V.K. All authors have read and agreed to the published version of the manuscript.

Funding

The work was carried out according to the Subject AAAA-A17-117080110043-4 «Dynamics of physical processes in the active zones of near space and geospheres» IKIR FEB RAS. The study was carried out with the support of the Russian Foundation for Basic Research within the framework of scientific project No. 19-31-50033. Research topic of Vitus Bering KamSU “Mathematical model of the Duffing oscillator with a derivative of variable fractional order” No. AAAA-A20-120021190004-8.

Acknowledgments

I would like to express my gratitude to Arsen Pskhu for valuable advice, which served to better understand the research results.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AFCAmplitude–Frequency Characteristics
PFCPhase-Frequency Characteristics
Q-factorQuality Factor

References

  1. Volterra, V. Sur les equations integro-differentielles et leurs applications. Acta Math. 1912, 35, 295–356. [Google Scholar] [CrossRef]
  2. Vlasov, V.; Rautian, N. A Study of Operator Models Arising in Problems of Hereditary Mechanics. J. Math Sci. 2020, 244, 170–182. [Google Scholar] [CrossRef]
  3. Li, Z.; Yang, Q. Systems and synthetic biology approaches in understanding biological oscillators. Quant Biol. 2018, 6, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Danylchuk, H.; Kibalnyk, L.; Serdiuk, O. Study of critical phenomena in economic systems using a model of damped oscillations. SHS Web Conf. 2019, 65, 06008. [Google Scholar] [CrossRef]
  5. Oldham, K.; Spanier, J. The Fractional Calculus. Theory and Applications of Differentiation and Integration to Arbitrary Order; Academic Press: London, UK, 1974; p. 240. [Google Scholar]
  6. Miller, K.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differntial Equations; A Wiley-Interscience Publication: New York, NY, USA, 1993; p. 384. [Google Scholar]
  7. Kilbas, A.; Srivastava, H.; Trujillo, J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204, p. 523. [Google Scholar]
  8. Syta, A.; Litak, G.; Lenci, S.; Scheffler, M. Chaotic vibrations of the duffing system with fractional damping. Chaos Interdiscip. J. Nonlinear Sci. 2014, 24, 013107. [Google Scholar] [CrossRef]
  9. Jiménez, S.; Zufiria, P. Characterizing chaos in a type of fractional Duffing’s equation. Conf. Publ. Am. Inst. Math. Sci. 2015, 2015, 660–669. [Google Scholar] [CrossRef]
  10. Xing, W.; Chen, E.; Chang, Y.; Wang, M. Threshold for chaos of a duffing oscillator with fractional-order derivative. Shock Vib. 2019, 2019, 1230194. [Google Scholar] [CrossRef]
  11. Liu, Q.X.; Liu, J.K.; Chen, Y.M. An analytical criterion for jump phenomena in fractional Duffing oscillators. Chaos Solitons Fractals 2017, 98, 216–219. [Google Scholar] [CrossRef]
  12. Eze, S.C. Analysis of fractional Duffing oscillator. Rev. Mex. Fis. 2019, 66, 187–191. [Google Scholar] [CrossRef]
  13. Wang, Y.; An, J.Y. Amplitude–frequency relationship to a fractional Duffing oscillator arising in microphysics and tsunami motion. J. Low Freq. Noise Vib. Act. Control 2019, 38, 1008–1012. [Google Scholar] [CrossRef]
  14. Ejikeme, C.L.; Oyesanya, M.O.; Agbebaku, D.F.; Okofu, M.B. Solution to nonlinear Duffing oscillator with fractional derivatives using homotopy analysis method (HAM). Glob. J. Pure Appl. Math. 2018, 14, 1363–1388. [Google Scholar]
  15. Atanackovic, T.; Pilipovic, S. Hamilton’s principle with variable order fractional derivatives. Fract. Calc. Appl. Anal. 2011, 14, 94–109. [Google Scholar] [CrossRef] [Green Version]
  16. Kim, V.A. Duffing oscillator with external harmonic action and variable fractional Riemann–Liouville derivative characterizing viscous friction. Bull. KRASEC Phys. Math. Sci. 2016, 13, 46–49. [Google Scholar]
  17. Oseledets, V.I. A Multiplicative ergodic theorem. Characteristic Lyapunov exponents of dynamical systems. Proc. Mosc. Math. Soc. 1968, 19, 197–231. [Google Scholar]
  18. Wolf, A.; Swift, J.B.; Swinney, H.L.; Vastano, J.A. Determining Lyapunov exponents from a time series. Physica 1985, 16, 285–317. [Google Scholar] [CrossRef] [Green Version]
  19. Kim, V.A.; Parovik, R.I. Calculation of maximal Lyapunov exponents for an oscillatory Duffing system with power memory. Bull. KRASEC Phys. Math. Sci. 2018, 23, 98–105. [Google Scholar] [CrossRef]
  20. Pskhu, A.; Rekhviashvili, S. Analysis of Forced Oscillations of a Fractional Oscillator. Tech. Phys. Lett. 2018, 44, 1218–1221. [Google Scholar] [CrossRef]
  21. Parovik, R.I. Amplitude-frequency and phase-frequency performances of forced oscillations of a nonlinear fractional oscillator. Tech. Phys. Lett. 2019, 45, 660–663. [Google Scholar] [CrossRef]
  22. Parovik, R. Mathematical Modeling of Linear Fractional Oscillators. Mathematics 2020, 8, 1879. [Google Scholar] [CrossRef]
  23. Kim, V.A.; Parovik, R.I. Investigation of forced oscillations of a Duffing oscillator with a derivative of a variable fractional Riemann–Liouville order. Izv. Kabard. Balkar Sci. Cent. Russ. Acad. Sci. 2020, 93, 46–56. [Google Scholar] [CrossRef]
  24. Huang, X.L.; Ma, X.; Hu, F. Machine learning and intelligent communications. Mob. Netw. Appl. 2018, 23, 68–70. [Google Scholar] [CrossRef] [Green Version]
  25. Pirmohabbati, P.; Sheikhani, A.R.; Najafi, H.S.; Ziabari, A.A. Numerical solution of full fractional Duffing equations with Cubic-Quintic-Heptic nonlinearities. AIMS Math. 2020, 5, 1621–1641. [Google Scholar] [CrossRef]
Figure 1. Oscillogram (a) and phase trajectory (b).
Figure 1. Oscillogram (a) and phase trajectory (b).
Mathematics 08 02063 g001
Figure 2. Range of Lyapunov exponents from λ (a) and phase trajectories at (b) λ = 0.47 , (c) λ = 0.64 and (d) λ = 0.15 .
Figure 2. Range of Lyapunov exponents from λ (a) and phase trajectories at (b) λ = 0.47 , (c) λ = 0.64 and (d) λ = 0.15 .
Mathematics 08 02063 g002
Figure 3. (a) Range of Lyapunov exponents from q and (b) phase trajectories at q = 0.6 and (c) q = 0.75 .
Figure 3. (a) Range of Lyapunov exponents from q and (b) phase trajectories at q = 0.6 and (c) q = 0.75 .
Mathematics 08 02063 g003
Figure 4. Q-factor: (a) A = 1.5 , t 96.4 , 97.5 ; (b) q t = sin 2 ω t , t = 97.4 .
Figure 4. Q-factor: (a) A = 1.5 , t 96.4 , 97.5 ; (b) q t = sin 2 ω t , t = 97.4 .
Mathematics 08 02063 g004
Figure 5. (a) AFC and (b) PFC with A = 1.5 , t 96.4 , 97.5 , q t = sin 2 ω t .
Figure 5. (a) AFC and (b) PFC with A = 1.5 , t 96.4 , 97.5 , q t = sin 2 ω t .
Mathematics 08 02063 g005
Table 1. Results of the study of the computational accuracy of the method (6).
Table 1. Results of the study of the computational accuracy of the method (6).
Nh ε α
100.10.0046870831.790194986
200.050.0022354791.654513226
400.0250.0010913171.556442049
800.01250.0005414371.481975798
1600.006250.000273771.422116237
3200.0031250.0001444431.368516573
6400.00156250.00008526321.309611649
12800.000781250.00006411111.230273342
25600.0003906250.0000671121.125073091
51200.00019531250.0000906361.008079025
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, V.; Parovik, R. Mathematical Model of Fractional Duffing Oscillator with Variable Memory. Mathematics 2020, 8, 2063. https://doi.org/10.3390/math8112063

AMA Style

Kim V, Parovik R. Mathematical Model of Fractional Duffing Oscillator with Variable Memory. Mathematics. 2020; 8(11):2063. https://doi.org/10.3390/math8112063

Chicago/Turabian Style

Kim, Valentine, and Roman Parovik. 2020. "Mathematical Model of Fractional Duffing Oscillator with Variable Memory" Mathematics 8, no. 11: 2063. https://doi.org/10.3390/math8112063

APA Style

Kim, V., & Parovik, R. (2020). Mathematical Model of Fractional Duffing Oscillator with Variable Memory. Mathematics, 8(11), 2063. https://doi.org/10.3390/math8112063

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop