On a credit oscillatory system with the inclusion of stick-slip

The work was a mathematical model that describes the effect of the sliding attachment (stick-slip), taking into account hereditarity. This model can be regarded as a mechanical model of earthquake preparation. For such a model was proposed explicit finitedifference scheme, on which were built the waveform and phase trajectories hereditarity effect of stick-slip.


Introduction
The effect of stick-slip occurs into tribology problems.For example, we study movement of cargo on a spring along the surface.Due to the adhesion of the goods adheres to the surface, and due to the spring tension breaks off and slides along it at the same place his hesitation [1][2][3].Also stick-slip effect may be laid in the foundation of the mechanical model of the earthquake in the subduction zone of lithospheric plates [4].
In this paper we will explore hereditarity nonlinear oscillator which characterizes the stick-slip effect [5,6].Recall that in [5,6], we considered the model equation stick-slip effect, c considering hereditarity with derivatives of fractional order permanent.
In this paper we consider the oscillator hereditarity stick-slip effect with fractional derivatives variables governmental orders, which are bounded functions of.And the orders of fractional derivatives may depend not only on the transfer variable t, and shift function x (t).Introduction thus fractional control parameters in the model equation generalizes considered into the third chapter hereditarity nonlinear oscillator based on stick-slip effects that can lead to new properties or effects as a result of mathematical modeling.
Next, we construct an explicit finite-difference scheme for the numerical calculation of the solution of the generalized model of stick-slip effect with variable fractional order.Next, look at some examples of the application of this scheme with different values of the control parameters, we plot and investigate phase trajectories.
Note.Gerasimov-Caputo Operators into the model equation ( 1) can be written as: Therefore, considering the representations (3) are placed two Cauchy problems (1) and ( 2).Note.The Cauchy problem ( 1) and ( 2) when the fractional parameters β and γ are constants, becomes the Cauchy problem in [5], if these constants are β = 2 and γ = 1, we arrive at the Cauchy problem for the classical stick-slip effect [3].
The Cauchy problem (1) and ( 2) generally has no clear decision, so we will solve this problem numerically using the theory of finite difference schemes [7,8].
We construct an explicit finite-difference scheme, because it is easiest to implement on a computer.We divide the time interval [0, T ] for N equal parts increments τ = T/N.Then the function solutions x (t) of the differential of the Cauchy problem (1) and (2) pass grid function in x t j , where t j = jτ, j = 0, 1, • • • , N − 1, which will be solution of the corresponding Cauchy difference problem.

Solar-Terrestrial Relations and Physics of Earthquakes Precursors
Note.Explicit finite difference scheme as are usually conditional stable, ie, there is a limitation on the step τ of difference grid.The scheme ( 6) are no exceptions.The step of the grid can be controlled by the method of double recalculation or Runge rule.In this paper we carry out a numerical experiment and see how the maximum absolute value error between the exact and numerical solutions, depending on the increase in the number of N nodes of the computational grid.
Note.The explicit finite difference scheme (6) in the case where the fractional parameters are constants and have values of β = 2 and γ = 1, become to the local settlement explicit finite-difference scheme for the classical optical effect stick-slip [3].

Results of numerical simulation
In [3] it was said that it is enough to take the first seven coefficients a k in the expansion of the function f (x(t), t) to obtain reliable solutions (1).The values of these coefficients take from [3] a 1 = 0.436, a 2 = 0.344, a 3 = 0.164, a 4 = 0.058, a 4 = 0.021, a 6 = 0.004, a 7 = 0.003.
Numerical modeling hereditarity stick-slip effect on the formulas (6) and obtained in [5] depending on the control parameters.First, we consider the dependence of the fractional control parameters only on the time t, ie, β = β (t) and γ = γ (t).Then we look at the more general case, when the fractional parameters also depend on x (t) of displacement.
In Fig. 1 shows the calculated curves of displacement, velocity displacement and the phase trajectory.In Fig. 1a shows an oscillogramma for example 1.It can be seen that in the separation of the load fluctuates, the rate of these oscillations in the potential well decays slowly enough (Fig. 1b).This effect is ereditarnostyu process.In the phase trajectory (Fig. 1c) shows that potential wells are stable focuses.
Example 2. The control parameters: In Fig. 2 shows the calculated curves obtained by numerical formulas (6) and [5].Unlike the examples discussed above, in this example the fractional parameters depend not only on the time t, but on the x (t) displacement.
In Fig. 2a shows the oscillogrammas obtained by the formulas (6) and [5], which characterize the effect of stick-slip given hereditarity.You may notice that the goods slide along the surface, testing damped oscillations of small amplitude.It can be seen that the amplitude of these oscillations increases slightly for the subsequent potential wells (Fig. 2b).Phase trajectories to this case are given in Fig. 2c.