The paper presents a model of fractal parametric oscillator. Showing that the solution of such a model exists and is unique. A study of the solution with the aid of diagrams Stratton-Ince. The regions of instability, which can occur parametric resonance. It is suggested that this solution can be any signal, including acoustic.
It is known that the natural medium (geological medium) may have fractal properties. These properties characterize the spatial-temporal nonlocality or “memory” of the medium, which in turn is determined by the power laws.
Usually geological medium with fractal properties described in terms of fractional calculus [
In this paper the nonlinear parametric oscillatory process in the geological medium with fractal properties. Feature of this process is that the displacement of the points geomedium a result of its stress-strain of the state of can occur with increasing amplitude due to changes in the parameters of the medium itself.
Through this process may be described cracks in the avalanche geomedium, which in most cases is preceded by seismic activity in the region (Kamchatka), which can be used in prediction of strong earthquakes.
The geomedium formed loose deposits of rocks. Assume that this medium has fractal properties. Then the problem of displacement geomedium points in time can be stated as follows:
with the initial conditions of the problem
Here shift function geomedium;
—the generalized cosine with parameter. Taking get the usual cosine, i.e., [
—the fractional differential order, and -parameters of medium. Assume, and parameters of geomedium, as they are depend on its fractal dimension.
The Equation (1) is a generalization of parametric resonance the classical Mathieu’s equation in a case.
Note if put, then Equation (1) is known as equation of fractional oscillator, which is investigated in study [
Since Equation (1) is considered first, then call it a fractal equation parametric oscillator.
In study [
The kernel of the Equation (3)
and
.
If put in solution (3), it is the solution of the fractional oscillator
Solve the Equation (3), use the composite trapezoidal quadrature formula. Take a grid with step. Put in (3), obtain:
The integral in expression (5) approximate the sum, considering, obtain:
—are coefficients of the quadrature formula,—the approximation error. Solution (6) can be reduced a system of algebraic equations:
The denominator of (7) must satisfy. This condition can be achieved by changing the step. Trapezoidal quadrature formula on the interval has an error, and the total error in the segment-.
The numerical solution of (7) allows the study of fractal parametric oscillator in particular can make the visualization of calculation results.
Numerical simulation of (1) and (2) was realized using a mathematical software MAPLE. First there was the case when in (1). It’s the classical Mathieu equation. It is known that solution of the Cauchy problem (1)- (2), taking and can be written in terms of the Mathieu function. The MAPLE gives the following result:
The solution (8), will used as a test for the analysis of the numerical solution obtained by the method (7). The simulation results for of the method (7) and formula (8) are shown in
In
According to this diagram, it’s impossible to determine at what values of parameters A and m parametric resonance occurs, for example when A = m = 1 parametric resonance occurs in
In
In this case the solutions haven’t a property of periodic, and have a decaying character. These solutions are characteristic of media with dissipation, in particular for inhomogeneous, fractal mediums.
In
Curves are also damped character at short times behave the same, while at large time intervals is regrouping of curves in the reverse order.
Indeed the values of parameters ξ and δ are entered the so-called zone of instability that can be constructed using of diagrams Strutt-Ince [
In
According to this diagram, it’s impossible to determine at what values of parameters ξ and δ parametric resonance occurs, for example when parametric resonance occurs in
Consider the differential Equation (1) and fractional derivatives for and:
Define the conditions under which there is a parametric resonance in (9). To do this, in the plane must construct diagrams Strutt-Ince. As a rule, there is a region of instability, parametric resonance, which leads to an increase in the amplitude of oscillations Usually in the area of instability exists parametric resonance, which leads to an increase in amplitude. Estimate the parameters δ.
Consider the derivative of fractional order on the left side of (9):
Use the method of harmonic balance for Equation (9), its solution formed of a harmonic series [
For the first resonance take the first harmonic (11), i.e. and substitute (9) in view of the representation (10). After some transformations go to the following result:
If in (12) to put, get the known relation for the first classical parametric resonance Mathieu
In
It can be noted that in (12) imposes constraints on the parameters. The value must satisfy the following inequality:
Spend the visualization of the results of research solutions of (9). According to the above analysis, it was the expression (12). Below is its visualization: