Fractional Derivatives and Dynamical Systems in Material Instability

: Loss of stability is studied extensively in nonlinear investigations, and classiﬁed as generic bifurcations. It requires regularity, being connected with non-locality. Such behavior comes from gradient terms in constitutive equations. Most fractional derivatives are non-local, thus by using them in deﬁning strain, a non-local strain appears. In such a way, various versions of non-localities are obtained by using various types of fractional derivatives. The study aims for constitutive modeling via instability phenomena, that is, by observing the way of loss of stability of material, we can be informed about the proper form of its mathematical model.


Introduction
There are materials where tests justify models with fractional order derivatives. Moreover, all fractional derivatives are non-local. Thus, by using strain calculated by fractional derivation of the displacement field, a non-local quantity appears instead of the conventional (local) strain. Fractional calculus appears to be a powerful tool to deal with forms of non-localities. When the term non-locality is used in its original meaning, the value of some quantity in an internal point of the body is determined by the values of other quantities in a whole region around that location.
In material instability problems, like shear band or neck formation [1], the post-critical behavior plays an essential role. Especially for numerical analysis [2], the form and thickness of instability zones should be mesh-independent for a proper material model. Such property implies a generic method of the loss of stability in the sense of dynamical systems theory [3,4]. Thus, the complete mathematical description of the solid continuum (a set of differential equations) should undergo a generic bifurcation. While that model is more or less fixed, the only element of the set of basic equations to be varied is the constitutive equation. Consequently, one should use a constitutive equation, which, by adding the equations of motion and the kinematic equation, forms a system, having generic bifurcation at loss of stability. Conventionally, such constitutive equations are formed by gradient regularization [5], which can also be interpreted as addition of non-locality.
Non-locality is an essential part of continuum physics from the early period of rational mechanics. Two forms are distinguished, the weak and the strong non-localities [6]. The earliest distinction was introduced by Maugin in 1977 at a conference in Warsaw [7]. Weak non-locality concepts that use gradient models could be better referred to as gradient theories of the n-th order [6]. At strong non-locality, constitutive equations become integral expressions over space, perhaps with a more or less rapid attenuation with the distance of the spatial kernel [6]. In a continuum physicist's point of view, the use of fractional derivatives fits well into the generalization of strong non-locality theory [8].
Weak non-localities can be introduced like Aifantis [5], by adding gradient term(s) to constitutive equations Term ∇ 2 ε is used in Equation (1) as in [5], to express second "gradient dependence" of material [8]. (Obviously, such an operator is the Laplace operator, ∇ 2 ≡ ∆). Gradient effects have several physical interpretation including internal lengths or micro-structural effect theories. Early studies from the 1960s use the notion of polar bodies [9] and couple stress effect and then a gradient of deformation gradient appears, which can be interpreted as a constitutive equation containing a first order strain gradient.
In case of strong non-locality, local stress is determined by the strain in a neighborhood [10]. Similarly to [11], assume that such interval in a uniaxial case is [a, b]. Then, integration is used to encounter that with attenuation function c (ξ), which reads by transforming into the displacement field and assuming attenuation function as In [11], one could see how Formula (2) can easily be generalized to Caputo's fractional derivative. It shows why non-locality is often modeled by using non-local derivatives [12] and then constitutive equations are in the forms of fractional differential equations [13]. This work deals with the consequences of fractional material models in material instability investigations.
The structure of this paper is as follows. Firstly, a short overview of basic fractional calculus is presented, including Euler's gamma function. A wide range of detailed monographs of that topic is available like [14] or [13] and the reference in them. Secondly, dynamical systems theory is introduced and applied to material instability problems to explain the importance of a property called generic in bifurcation analysis. The most important part is the following section, in which conditions of generic bifurcation behavior are studied for various types of constitutive equations using various fractional derivatives.

Fractional Integral Operators and Fractional Derivatives
Fractional derivatives are constructed from a fractional order generalization of Cauchy's integra formula. In such a way, Riemann-Liouville's fractional integral operator is defined, where 0 < α < 1 and Euler's gamma function Γ (α) is an integral of the two most important functions of analysis. These are power function and exponential functions and then It can also be interpreted as a generalization of factorial to non-integers because, for integer α = n, definition (4) results in being factorial, Integration in Formula (3) goes from left to right, (a < x); for that reason, it is called the left integral operator. By taking a derivative of the left Riemann-Liouville integral operator, a fractional derivative can be defined. The usual form of the left Riemann-Liouville derivative is obtained from Formulas (3) and (6) a D α which is a non-local fractional derivative for interval [a, x]. By varying the interval of integral operator in Formula (7) to [x, b], the so-called right Riemann-Liouville derivative can be defined as (6)

By changing operators of derivation and integration in Formula
and then its integro-differential operator formula reads

Symmetric Fractional Derivatives
For most applications in non-local mechanics, fractional derivative ∂ α ∂x α should be a symmetric derivative. For example, such derivative can be constructed from left and right Riemann-Liouville or Caputo derivatives (11) or similarly from other types of left and right fractional derivatives. As a famous example in fractional mechanics, a symmetric derivative [15] is presented. Let us start from equivalent forms of the left and right Caputo derivatives like [16]. From (8), Then, take "asymptotic cases" (a → ∞, b → ∞) and combine them into a symmetric form of Liouville-Caputo derivatives, like which is closely related to the Riesz potential [15].

On a Two-Sided Fractional Derivative as Infinite Series
There is a more efficient way to construct symmetric fractional derivatives. Its origin is at Grünwald-Letnikov formulation. For repeated n th backward differentiation [17] D n Then, instead of integer (n), take fractional order (α) difference as an infinite series Now, for h > 0, the left-and right-sided Grünwald-Letnikov derivatives are and a symmetric derivative can be constructed by adding them together. In a similar way, so-called centered derivatives could be defined and even integral forms are available [18]. As a generalized version, in this study, the two-sided fractional derivatives (TSFD) are used as it was introduced by where β > −1. The solution of the eigenvalue -eigenvector problem for differential operator (12) is and this fact will be used extensively in our stability investigation.

A Dynamical Systems Approach of Material Instability
The classical description of continuum mechanics consists of the set of basic equations of continua. In this section, a dynamical system is constructed from these equations and material instability of a state of material is studied as an instability of a solution of that dynamical system.
In the simplest possible case, a uniaxial problem is studied with small deformations. Then, basic equations are

•
Cauchy's equations of motion • the kinematic equation in rate formε • and the constitutive equation Note that constitutive Equation (16) is used here in its general form, a function of strain field ε, stress field σ, and various derivatives of them. Moreover, in addition to the variables shown in Equation (16), further dependence may and will appear later in this paper. Strain gradient and fractional strain gradient may also be included into the constitutive equation.
From the basic Equations (14)-(16), a dynamical system can be defined where By studying stability and bifurcation, let us assume that is an equilibrium of Equation (17), Then, by introducing homogeneous perturbations, the linearization of Equations (17) and (18) is satisfied for the perturbations. In that sense, homogeneous perturbation for a uniaxial continuum of length L means that In mechanics, the case of adiabatic localization [20] may justify the case of infinite length [15], (L → ∞), especially for material instability problems. Such case assumes that the effect of the boundaries can be eliminated sufficiently far from both ends of a rod, which correlates the experimental observations. Shear bands or necking regions appear generally "in the middle" of a specimen, sufficiently far from both ends.
As it is presented in [21], the eigenvalues and eigenvectors of linearized operator play the key role in stability and bifurcation analysis. For this reason, the characteristic equation of Formula (22) (LF − λI) y = 0 should be solved and stability conditions are formed for its solutions λ i in a usual way: • state y 0 of the material is stable, if Re (λ i ) < 0 for all solutions of Equation (23), • loss-of-stability happens, when at least for one i, Two ways of instabilities can be distinguished, and • the dynamic bifurcation (Re (λ) = 0, Im (λ) = 0) , for non-trivial y c . Nonlinear analysis can be performed, if operator defined by Formula (22) has non-trivial critical eigenspace, thus Equation (23) should have a non-trivial solution y c at Re (λ) = 0. Such case is referred to as generic bifurcation and nonlinearity should be projected to this eigenspace resulting in bifurcation equations [21]. For this reason, the existence and exact determination of y c are of critical importance. The following section deals with that problem in the case of various constitutive equations.
Characteristic Equation (23) is a partial differential equation, while Formula (22) is a differential operator. Its solution, or its solvability is far from obvious even in the case of integer order derivatives. A kind of simplifying approximation is to restrict the problem to some proper functions. Such functions are harmonics e iκx in the so-called periodic perturbation technique [21]. When homogeneous boundary conditions given by Formulas (19)-(21) are taken into consideration, the perturbation frequencies are κ = nπ L , n = 1, 2, . . . .

Constitutive Equation with Non-Local Strain
Firstly, strong non-locality is assumed. In such case, local stress σ is determined by the strain ε in a neighborhood [10]. A generalization of such approach is the use of a non-local strain, where, instead of local strain ε = 1 2 (u∇ + ∇u) , a non-local one is used, defined by a symmetric fractional derivative (see Formulas (10)- (12)). In Formula (26), absolute value points out the symmetric nature of derivative. Then, a static bifurcation is studied at non-local strain as in [15]. Let the constitutive equation be the so-called Malvern-Cristescu equation and let us extend the scope of investigation to dynamic bifurcation too. Now, characteristic Equation (23) has the form Static bifurcation condition in Formula (24) to Equation (28) result in

Malvern-Cristescu Constitutive Equation at TSFD
By using TSFD, static bifurcation condition reads Formula (29) and for all periodic perturbations e iκx the critical tangent stiffness is E = 0. Here, a kind of lack of generic behavior is that all perturbations are critical. This means that no finite dimensional critical eigenspace can be determined. For the dynamic bifurcation (condition in Formula (25)), by substituting λ = iω into Equation (28), is obtained. Assume that κ > 0; then, by substituting Formula (13) into Equation (30), the bifurcation condition is The next step is to get real and imaginary partition of Equation (31). The real part leads to while the imaginary part implies The solution of the "real part" for ω reads Then, the imaginary part results and the solution for the "real part" is Now, ω should be substituted from Formula (35) into Equation (34) and conditions should be formed for constitutive parameters E, D, H for various values of perturbing frequency κ. The next subsection, for the sake of simplicity, deals with the case of Riesz derivative as a special case of TSFD.

Malvern-Cristescu at Riesz Derivative
To study the dynamic bifurcation of Malvern-Cristescu equation, one should return to Equations (32) and (33). In such case, these are and While ω = 0, from Equation (36) the dynamic bifurcation condition then reads Remark that either D = 0 or H = 0 leads to the "classical" coexistent dynamic and static bifurcation case at E = 0.
In Formula (39), either D or E should be negative, but the sign of H has no effect on dynamic bifurcation. By expressing E from Formula (39), its value is plotted in Figure 1, with data D = −1, H = 1, ρ = 0.01, β 1 = −0.2, β 2 = −0.55, β 3 = −0.8. Now, for a critical perturbation frequency κ crit , the critical tangent stiffness can be determined and vice versa.

Classical Visco-Elasto-Plastic Case at Riesz Derivative
For a classical visco-elasto-plastic material, D = 0. While θ = 0, β > −1 represents the case of symmetric Riesz derivative [19], Equation (34) implies and Equation (32) implies Unfortunately, Formula (40) is the same as the static bifurcation condition, thus such constitutive equation results in non-generic coexistent static and dynamic bifurcations, consequently unavailable for a dynamic stability analysis. Additionally, non-generic behavior appears also in the fact that all periodic perturbations are critical ( all values of κ can be selected as a critical perturbing frequency).

Fractional Gradient Material
As an example of strong non-locality, second gradient materials are studied. As it is used by Aifantis [5], a second gradient term should be included into constitutive equations and Formula (1) is obtained. In the Aifantis-Tarasov material model [22], the idea is to use fractional derivative in Formula (1). In uniaxial form, such equation reads However, non-locality is already involved in Equation (42) and the use of local "fractional" derivatives might even be possible in the sense of non-local mechanics ( [23,24]). However, such derivatives are widely criticized [25]; for this reason, TSFD will be used again: In this case, characteristic Equation (23) reads Then, for periodic perturbations e κ x i (selected to be the eigenfunctions of D β θ ), is obtained. Equation (45) can be solved to λ From Formula (46), static bifurcation condition is in case of symmetric Riesz derivative (θ = 0, β > −1) By expressing B from Equation (47) as a function of perturbation frequency κ, its graph is plotted in Figure 2 at (C = −1) . Dynamic bifurcation study requires adding rate dependence to the constitutive Equation (43) Then, the characteristic equation is for eigenfunctions e κ x i Equation (50) can be solved for λ A necessary condition for dynamic bifurcation is Notice that the inequality in Formula (52) prescribes that parameters should be below the curves in Figure 2. For example, the stability chart for β = −0.8 is presented in Figure 3. When D > 0, and the parameters are in the shaded region, the state of the material is stable.

Non-Local Strain Gradient Material
In case of a generalization of the classical second gradient dependent material [5], the constitutive equation reads In material model of Formula (53), both weak (second gradient term) and strong non-localities (fractional strain) are present.
In this case, a critical perturbation can be identified. Then, static bifurcation condition in Formula (24) for periodic perturbations Equation (54) implies that Then, critical perturbation frequency is and the non-trivial eigenfunction is From Formulas (55)-(57), we find no dependence on the type of fractional derivative, and the result is the same as the integer order gradient problem [21].

Conclusions
Stability analysis was performed for non-local materials. Non-locality was modeled by using fractional derivatives. The main aim was to guarantee non-trivial critical eigenfunctions at the loss of stability. Such eigenfunctions play an important role in nonlinear studies because nonlinearities are projected to them to form bifurcation equations. This case is referred to as generic bifurcation. The existence of non-trivial critical eigenfunctions was studied for various non-local constitutive equations.
Two main groups distinguished the strong and the weak types of non-localities. In the case of strongly non-local materials, the introduction of fractional strain via fractional derivative seems to be an obvious generalization of the classical theories. Fractional derivatives here should be non-local and symmetric. The so-called two-sided fractional derivatives by [19] satisfy all these requirements and its special version, the Riesz derivative, was used in all bifurcation investigations.
Constitutive relations were classified as equations with non-local strain and fractional gradient dependent ones. The first type was represented by the Malvern-Cristescu equation and a simplified version of it, called the classical visco-elasto-plastic material. The second group was represented by the fractional gradient material (Aifantis-Tarasov model), its rate dependent version. At the end, the non-local strain gradient material was studied.
The main results are as follows: • For the Malvern-Cristescu equation at static bifurcation in non-generic because no finite dimensional critical eigenspace exists. By studying dynamic bifurcation, generic behavior was detected as it can be observed in (39) and in Figure 1.

•
In case of classical visco-elasto-plastic material, no generic nature can be assumed. Both infinite dimensional critical eigenspace and coexistent static and dynamic bifurcations are found. Thus, such equation cannot be used for any nonlinear stability analysis.

•
For fractional gradient materials, a generic static bifurcation exists with condition (47). In Figure 2, the curve of critical tangent stiffness and critical eigenfunction was presented. Unfortunately, dynamic bifurcation is non-generic, but a stability chart can be drawn, Figure 3.

•
In the case of non-local strain gradient material, both strong and weak non-localities are present in the constitutive equation. The results are the same as the simple second gradient material [21].
Funding: This research received no external funding.

Conflicts of Interest:
The author declares no conflicts of interest.