Next Article in Journal
Performance Analysis and Simulation of IRS-Aided Wireless Networks Communication
Previous Article in Journal
Analysis of B(s)0μ+μ Decays at the Large Hadron Collider
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Dynamics of a New Class of Micro-Electromechanical Oscillators—Open Problems

1
Faculty of Mathematics and Informatics, University of Plovdiv Paisii Hilendarski, 24, Tzar Asen Str., 4000 Plovdiv, Bulgaria
2
Institute of Mathematics and Informatics, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Bl. 8, 1113 Sofia, Bulgaria
3
Faculty of Mathematics and Informatics, Sofia University “St. Kliment Ohridski”, 5, James Bourchier Blvd., 1164 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(2), 253; https://doi.org/10.3390/sym16020253
Submission received: 12 January 2024 / Revised: 8 February 2024 / Accepted: 15 February 2024 / Published: 19 February 2024

Abstract

:
In this paper, we propose a new class of micro-electromechanical oscillators. Some investigations based on Melnikov’s approach are applied for identifying some chaotic possibilities. We demonstrate also some specialized modules for investigating the dynamics of these oscillators. This will be included as an integral part of a planned much more general Web-based application for scientific computing. It turns out that the theoretical apparatus for studying the circuit implementation (design, fabricating, etc.) of the considered differential model for large parameter values is extremely complex and requires a serious investigation. This is the reason to offer this model to the attention of specialists working in this scientific direction. Some open problems related to the use of existing computer algebraic systems for the study of this class of oscillators for large values of n , m and N are also posed. In general, the entire article is subordinated to this frank conversation with the readers with the sole purpose being the professional upgrading of the specialized modules provided for this purpose in subsequent licensed versions of CAS.

1. Introduction

The micro-electromechanical oscillators are important apparatus that have a large usage in many different real life fields which need very fine clock synchronization. We can mention the telecommunications, video transfers, GPS systems, the work of all microprocessors and many others. This wide applicability leads to a necessity of further investigations of their theoretical foundations as well as the possible implementation in the practice.
A variety of micro-electromechanical oscillators is governed by a version of the Mathieu equation that harbors both linear and cubic nonlinear time-varying stiffness terms—the subharmonic resonance is investigated in [1]; the possible chaotic behavior generated by such models is studied in [2], etc. An application of the oscillator dynamics to the filtering problem is provided in [3]. The time-varying capacitors are included in the study of [4] to investigate the nonlinear dynamics of a micro-electromechanical system. A comprehensive study devoted on the topic can be found in [5].
Various results related to hetero/homoclinic bifurcation, optimal control in the nonlinear dynamics, higher-order Melnikov analysis of homo/hetero clinic bifurcation in mechanical oscillators, and bifurcation and chaos in a single d.o.f. mechanical system can be found in [6,7,8,9,10,11,12,13].
In [14], the authors study the dynamics of a class of hypothetical micro-electromechanical oscillator of the form
d x d t = y d y d t = b x i = 0 n 2 1 b i x n 2 i ϵ α y + ( 1 + cos ( ω t ) ) a x + i = 0 m 2 1 a i x m 2 i .
In the present article, we suggest a new trigonometric class of micro-electromechanical oscillators and study their dynamics. The first important task we investigate is the possible chaotic behavior of our model. Here, we use the method proposed by Melnikov in [15]. It gives us an analytic tool for establishing the existence of transfer homoclinic points of the Poincare map for a periodic orbit of a perturbed dynamical system. Some analytical conditions are obtained. Also, the effects of all free parameters in the planar differential system are analyzed.
Our second task is to consider the performance these new models produce. Several simulations are composed, and the generated results are discussed in the light of the theoretical and practical significance they lead. We demonstrate also some specialized modules for investigating the dynamics of these hypothetical oscillators. The derived results can be used as an integral part of a much more general application for scientific computing—for some details, see [14,16,17,18]. Other dynamic models can be found in [19,20,21,22,23,24,25,26].
Last but not least, some open problems related to the use of existing computer algebraic systems for the study of this class of oscillators for large values of n , m and N are also posed.
In the proposed article (for the first time), “friendly fire” is directed at the thousand-strong army of collaborators (full-time and voluntary) of companies creating CAS and Web-based applications (with paid or free access) for scientific computing.
The plan of the paper is as follows. We state our model in Section 2. Our results in the light of Melnikov’s approach can be found in Section 3. Some open problems are discussed in Section 4. Some simulations are presented in Section 5. An alternative to control the chaotic behavior under probability distributed parameters a i , i = 0 , 1 , , m 2 , b i , i = 0 , 1 , , n 2 , g j , j = 1 , 2 , , N is given in Section 5.1. A special case is considered in Section 6. In Section 7, we will focus on some difficulties the user encounters when using Computer Algebraic Systems for scientific computing. We conclude in Section 8.

2. The Model

We consider the following new class of micro-electromechanical oscillators
d x d t = y d y d t = b x i = 0 n 2 1 b i x n 2 i ϵ α y + j = 1 N g j sin ( j ω t ) a x + i = 0 m 2 1 a i x m 2 i ,
where 0 ϵ < 1 , n and m are odd, N is an integer, and the coefficients a, b, α , a i , i = 0 , 1 , , m 2 , b i , i = 0 , 1 , , n 2 , g j , j = 1 , 2 , , N are real numbers.
Let us consider, for example, the case n = m = 3 . We can fix without community restriction b and b 0 at b = b 0 = 1 . For ϵ = 0 , the resulting Hamiltonian of the system (2) is
H ( x , y ) = 1 2 y 2 1 2 x 2 + 1 4 x 4 .
There is a saddle point at the origin, centers at ( ± 1 , 0 ) , and a double homoclinic orbit given by (see Figure 1)
x 0 ( t ) = ± 2 s e c h t y 0 ( t ) = 2 s e c h t t a n h t .
We refer to [14,27,28,29,30] for more details. The potential energy V ( x ) = 1 2 x 2 + 1 4 x 4 is shown in Figure 2.

3. Results in the Light of the Melnikov’s Approach

The Melnikov function gives a measure of the leading order distance between the stable and unstable manifolds when ϵ 0 and can be used to tell where the stable and unstable manifolds intersect transversely. By definition, the homoclinic integral of Melnikov is given by
M ( t 0 ) = y 0 ( t ) α y 0 ( t ) + j = 1 N g i sin ( j ω ( t + t 0 ) ) a x 0 ( t ) + i = 0 m 2 1 a i x 0 m 2 ( t ) d t ,
where the functions x 0 ( t ) and y 0 ( t ) are defined by Equation (3). From a numerical point of view, the task of finding a root of M ( t 0 ) is more interesting given that the parameters appearing in the proposed differential model are subject to a number of restrictions of a physical and practical nature.

3.1. The Case n = m = 3 and N = 1

Suppose that n = m = 3 and N = 1 . We compute the integral of Melnikov replacing formulas (3) in Equation (4):
M ( t 0 ) = y 0 ( t ) α y 0 ( t ) + g 1 sin ( ω ( t + t 0 ) ) a x 0 ( t ) + a 0 x 0 3 ( t ) d t = ( 2 s e c h t t a n h t ) α ( 2 s e c h t t a n h t ) + + g 1 sin ( ω ( t + t 0 ) ) a 2 s e c h t ) + a 0 ( 2 s e c h t ) 3 d t = 4 3 α 1 6 g 1 π ω 2 ( 6 a + a 0 ( 4 + ω 2 ) ) c s c h π ω 2 cos ( t 0 ω ) ,
Thus, we obtain the following statement.
Proposition 1. 
If m = n = 3 and N = 1 , then the Melnikov function M ( t 0 ) has a root which is calculated from the nonlinear equation
4 3 α 1 6 g 1 π ω 2 ( 6 a + a 0 ( 4 + ω 2 ) ) c s c h π ω 2 cos ( t 0 ω ) = 0 .
For example, for α = 1.39627 , a 0 = 1.2 , a = 1 , ω = 2.2 , and g 1 = 0.7 , the root of the nonlinear Equation (6) is t 0 0.0003209 with multiplicity two (see Figure 3).
Remark 1. 
If M ( t 0 ) = 0 and d M ( t 0 ) d t 0 0 for some t 0 and some sets of parameters, then chaos occurs. This Melnikov criterion for the appearance of the intersection between the perturbed and unperturbed separatrixes leads to the following condition for chaotic behavior
g 1 > P Q ,
where
P = 1 6 π ω 2 ( 6 a + a 0 ( 4 + ω 2 ) ) c s c h π ω 2 Q = 4 3 α .
It is important to note that a large region of parameter space will satisfy this inequality.

3.2. The Case n = m = 3 and N = 2

Similarly to Proposition 1, we can prove the following proposition when n = m = 3 and N = 2 .
Proposition 2. 
If m = n = 3 and N = 2 , then the Melnikov function M ( t 0 ) has a root which is a solution of the nonlinear equation
M ( t 0 ) = 1 6 e 2 i t 0 ω c s c h ( π ω ) e i t 0 ω ( 1 + e 2 i t 0 ω ) g 1 π ω 2 ( 6 a + a 0 ( 4 + ω 2 ) ) c o s h π ω 2 + 4 ( ( 1 + e 4 i t 0 ω ) g 2 π ω 2 ( 3 a + 2 a 0 ( 1 + ω 2 ) ) 2 α e 2 i t 0 ω s i n h ( π ω ) ) = 0 .

3.3. The Case n = m = 3 and N = 3

The same approach leads to the following statement when n = m = 3 and N = 3 .
Proposition 3. 
If m = n = 3 and N = 3 , then the root of Melnikov function (4) can be obtained via the nonlinear equation
M ( t 0 ) = 1 24 ( 1 + 2 c o s h ( π ω ) ) e 3 i t 0 ω c s c h π ω 4 s e c h π ω 4 s e c h π ω 2 × ( 1 + e 2 i t 0 ω ) π ω 2 6 a e 2 i t 0 ω ( 2 g 1 9 g 3 ) + 9 g 3 + 9 e 4 i t 0 ω g 3 + + a 0 9 g 3 ( 4 + 9 ω 2 ) + 9 e 4 i t 0 ω g 3 ( 4 + 9 ω 2 ) + + e 2 i t 0 ω 2 g 1 ( 4 + ω 2 ) 9 g 3 ( 4 + 9 ω 2 ) ) c o s h π ω 2 + + e i t 0 ω 8 ( 1 + e 4 i t 0 ω ) g 2 π ω 2 ( 3 a + 2 a 0 ( 1 + ω 2 ) ) c o s h ( π ω ) + + e i t 0 ω ( 1 + e 2 i t 0 ω ) g 1 π ω 2 ( 6 a + a 0 ( 4 + ω 2 ) ) c o s h 3 π ω 2 + + 4 ( ( 1 + e 4 i t 0 ω ) g 2 π ω 2 ( 3 a + 2 a 0 ( 1 + ω 2 ) ) 2 α e 2 i t 0 ω s i n h ( π ω ) 2 α e 2 i t 0 ω s i n h ( 2 π ω ) = 0 .
For example, for α = 1.7 , a 0 = 0.8 , a = 1 , ω = 0.75 , g 1 = 0.7 , g 2 = 0.3 , and g 3 = 0.1 , the root of the nonlinear Equation (8) is t 0 0.144 (see Figure 4).

4. Open Problems

Problem 1. 
The reader can consider the corresponding approximation problem for arbitrarily chosen n , m and N.
For n > 3 , the study of the critical levels of the Hamiltonian H ( x , y ) is very complicated. In this regard, we recommend the study of Gavrilov and Iliev [31]. Note that even with the weakened choice n = 3 and m , N sufficiently large numbers, some difficulties (of a user nature) are encountered when calculating the Melnikov integrals using known Computer Algebraic Systems.
Problem 2. 
The Melnikov function M ( t 0 ) can be evaluated via the method of residuals.
Modern computer algebraic systems for scientific calculations provide this opportunity for the users. The upgrade of the Web application planned by us foresees the use of an algorithm (hidden to the user) to define the limit of the type | I m ( ω ) | C o n s t . For example, see Proposition 3, we calculate M ( t 0 ) but with a pre-set limit by us: | I m ( ω ) | 2 3 .

5. Some Simulations

Here, we will focus on some interesting simulations.
Example 1. For the given n = m = 3 , N = 3 , a = 1 , α = 2.7 , a 0 = 1.2 , b = b 0 = 1 , ϵ = 0.01 , ω = 2.2 , g 1 = 0.7 , g 2 = 0.8 , and g 3 = 0.9 the simulations on the system (2) for x 0 = 0.9 ; y 0 = 0.7 are depicted in Figure 5.
Example 2. For a given n = m = 5 , N = 4 , a = 1 , α = 2.9 , a 0 = 3.35 , a 1 = 2.2 , b = b 0 = b 1 = 1 , ϵ = 0.01 , ω = 2.1 , g 1 = 0.9 . g 2 = 1.1 , g 3 = 1.6 , and g 4 = 2.5 , the simulations on the system (2) for x 0 = 0.7 ; y 0 = 0.6 are depicted in Figure 6.
Example 3. For a given n = m = 7 , N = 4 , a = 1 , α = 2.9 , a 0 = 3.35 , a 1 = 3.35 ,   a 2 = 1.5 ,   b = 0.1 ,   b 0 = 0.3 , b 1 = 0.01 , b 2 = 0.2 , ϵ = 0.01 , ω = 2.1 , g 1 = 0.9 . g 2 = 1.1 ,   g 3 = 1.6 , and g 4 = 2.5 , the simulations on the system (2) for x 0 = 0.7 ; y 0 = 0.6 are depicted in Figure 7.
Problem 3. 
The task of researching the schematic implementation (design, fabrication, etc.) of the considered differential model can be considered open.
Note that one of the possibilities for conducting qualitative experimental studies for fabricated micro-electromechanical oscillator systems of type (2) and verifying their chaotic behavior is the use of some classical methods from the field of approximation theory, particularly from the direction approximations with constraints.
The characteristic restrictions for this class of oscillators, more precisely on the y-component of the differential model solution, are visualized in Figure 8.
Let us present one more example. With the fixed values of the parameters from Example 3, but with the special choice of the parameters g i of the form: g i = 0.9 0.2 ( i 1 ) ;   i = 1 , 2 , 3 , 4 (or a similar uniform network), a relatively good hit in the restriction fork is achieved (see Figure 9a).

5.1. Parameters Generated by a Probability Distribution

An alternative to control the mentioned chaotic behavior is to study this class of oscillators under probability distributed parameters a i , i = 0 , 1 , , m 2 , b i , i = 0 , 1 , , n 2 , g j , j = 1 , 2 , , N . Of course, many distributions can be considered, but we focus mainly on two examples—binomial and truncated geometric. The first one is proposed since its asymptotic behavior is the normal distribution, which may be the most used in all practical fields. In addition, the truncated geometric distribution is suggested, because the original geometric distribution is the unique discrete one which exhibits the property memoryless—a necessary feature for many real-life fields. This way, the values of a i , b i , and g i can be calculated through the probability mass functions of the mentioned above distributions—we denote them by p b i n and p T G . If we assume that their support is 0 , 1 , 2 , , M , then they are
p b i n k = M k p k 1 p M k , k 0 , 1 , 2 , , M , 0 < p < 1 p T G k = p 1 p k 1 1 p M + 1 , k 0 , 1 , 2 , , M , 0 < p < 1 .
Thus, the necessary coefficients a i , b i , and g i can be set to be equal to the corresponding values p b i n k or p T G k .

6. An Example with a Relatively Large Number of Terms: n = 5, m = 5

We suppose now that n = m = 5 , and consider the model (2) for b < 0 , b 1 < 0 , b 0 > 0 . For ϵ = 0 , the resulting Hamiltonian of the system (2) is
H ( x , y ) = 1 2 y 2 1 2 b x 2 + 1 4 b 1 x 4 + 1 6 b 0 x 6 .
The potential energy
V ( x ) = 1 2 b x 2 + 1 4 b 1 x 4 + 1 6 b 0 x 6
is shown in Figure 10. In [32], the authors study the resonant oscillation and the homoclinic bifurcation in a ϕ 6 —Van der Pol oscillator. Following the ideas given in [32], we obtain the Hamiltonian system for ϵ = 0 for our model (2) through a pair of heteroclinic orbits defined as (Figure 11b)
x h e t = ± 2 x 1 s i n h γ 2 t ξ + c o s h ( γ t ) y h e t = ± 2 γ x 1 ( 1 ξ ) c o s h γ 2 t 2 ( ξ + c o s h ( γ t ) ) 3 2
and a pair of symmetric homoclinic trajectories connected each unstable point to itself given by (Figure 11a)
x h o m = ± 2 x 1 c o s h γ 2 t ξ + c o s h ( γ t ) ; y h o m = ± 2 γ x 1 ( ξ 1 ) s i n h γ 2 t 2 ( ξ + c o s h ( γ t ) ) 3 2 ,
where
δ = b 1 2 + 4 b b 0 ρ = b 1 δ b 1 + δ ξ = 5 3 ρ 2 2 ρ 2 1 x 1 = 1 2 b 0 ( b 1 + δ ) γ = x 1 2 2 b 0 ( ρ 2 1 ) .
These results are presented in Figure 11. We refer to [32,33,34,35] for more details. Next, we present a few more illustrative examples.
Example 4. For a given n = m = 5 , N = 4 , a = 1 , α = 0.1 , a 1 = 0.1 ,   a 0 = 0.05 ,   b = 0.4 ,   b 0 = 0.1 , b 1 = 0.4 , ϵ = 0.01 , ω = 1.1 , g 1 = 0.9 , g 2 = 0.6 , g 3 = 0.3 , and g 4 = 0.01 , the simulations on system (2) for x 0 = 0.6 ; y 0 = 0.2 are depicted in Figure 12.
Example 5. With the fixed values of the parameters from Example 5, but with the special choice of the parameters g i of the form: g 1 = 1.7 , g 2 = 1.5 , g 3 = 1.3 ,   g 4 = 1 , and g 5 = 0.7 , the simulations on system (2) for x 0 = 0.7 ; y 0 = 0.2 are depicted in Figure 13.
We have to note that the calculation of homoclinic and heteroclinic Melnikov integrals as well as the corresponding criterion for the chaos occurrence in the dynamical system (for n = m = 5 , N 4 ) are problematic tasks for users. For other problems of a computational nature, see [35].

7. Some Difficulties That the User Encounters When Using Computer Algebraic Systems for Scientific Calculation

Here, we will focus only on those difficulties that the user encounters when calculating homoclinic Melnikov integrals.
First of all, we will note that the user (who is not necessarily to be a professional mathematician) can hardly find out what option to choose for the mandatory field in the existing specialized module (for example, in CAS Mathematica) for calculating the above-mentioned integrals. Usually, this requirement is related to the imposed constraint (for the class of dynamical models) of type | I m ( ω ) | C o n s t (see Section 4). The upgrade of the Web application planned by us foresees the use of an algorithm (hidden for the user) to define the limit of the above-mentioned constraint.
In a number of cases, a comment message appears (after a long pause) of the type: “No more memory available”, or ”SystemException[MemoryAllocationFailure, …]”, which is not particularly useful information for the nonprofessional user.
We note that when calculating the mentioned integrals for large values of the parameters in the considered differential model, the following comment message often appears: “A very large output was generated. Here is a simple of it:”. In this situation, the user’s hands are tied again.
We will give one more example (our unpublished result) related to studying the dynamics of an extended escape model of the form:
d x d t = y d y d t = x + x 2 ϵ A p = 1 n y | y | p 1 + j = 1 N g j sin ( j ω t )
where 0 ϵ < 1 , A is the damping level, p 1 is the damping exponent, and N is the integer. In particular, we consider the following model
d x d t = y d y d t = x + x 2 ϵ A y | y | p 1 + j = 1 N g j sin ( j ω t ) ,
The following statement is valid.
Proposition 4. 
If p = 2 and N = 3 , then the roots of Melnikov function M ( t 0 ) are given as solutions of the equation
M ( t 0 ) = 1 60 e 3 i t 0 ω 32 2 A e 3 i t 0 ω + 15 2 e 2 i t 0 ω g 1 ω P o l y G a m m a [ 0 , 1 4 i ω 4 ] + 15 2 e 4 i t 0 ω g 1 ω P o l y G a m m a [ 0 , 1 4 i ω 4 ] 15 2 e 2 i t 0 ω g 1 ω P o l y G a m m a [ 0 , 3 4 i ω 4 ] 15 2 e 4 i t 0 ω g 1 ω P o l y G a m m a [ 0 , 3 4 i ω 4 ] + 15 2 e 2 i t 0 ω g 1 ω P o l y G a m m a [ 0 , 1 4 + i ω 4 ] + 15 2 e 4 i t 0 ω g 1 ω P o l y G a m m a [ 0 , 1 4 + i ω 4 ] 15 2 e 2 i t 0 ω g 1 ω P o l y G a m m a [ 0 , 3 4 + i ω 4 ] 15 2 e 4 i t 0 ω g 1 ω P o l y G a m m a [ 0 , 3 4 + i ω 4 ] + 30 2 e i t 0 ω g 2 ω P o l y G a m m a [ 0 , 1 4 i ω 2 ] + 30 2 e 5 i t 0 ω g 2 ω P o l y G a m m a [ 0 , 1 4 i ω 2 ] 30 2 e i t 0 ω g 2 ω P o l y G a m m a [ 0 , 3 4 i ω 2 ] 30 2 e 5 i t 0 ω g 2 ω P o l y G a m m a [ 0 , 3 4 i ω 2 ] + 30 2 e i t 0 ω g 2 ω P o l y G a m m a [ 0 , 1 4 + i ω 2 ] + 30 2 e 5 i t 0 ω g 2 ω P o l y G a m m a [ 0 , 1 4 + i ω 2 ] 30 2 e i t 0 ω g 2 ω P o l y G a m m a [ 0 , 3 4 + i ω 2 ] 30 2 e 5 i t 0 ω g 2 ω P o l y G a m m a [ 0 , 3 4 + i ω 2 ] + 45 2 g 3 ω P o l y G a m m a [ 0 , 1 4 3 i ω 4 ] + 45 2 e 6 i t 0 ω g 3 ω P o l y G a m m a [ 0 , 1 4 3 i ω 4 ] 45 2 g 3 ω P o l y G a m m a [ 0 , 3 4 3 i ω 4 ] 45 2 e 6 i t 0 ω g 3 ω P o l y G a m m a [ 0 , 3 4 3 i ω 4 ] + 45 2 g 3 ω P o l y G a m m a [ 0 , 1 4 + 3 i ω 4 ] + 45 2 e 6 i t 0 ω g 3 ω P o l y G a m m a [ 0 , 1 4 + 3 i ω 4 ] 45 2 g 3 ω P o l y G a m m a [ 0 , 3 4 + 3 i ω 4 ] 45 2 e 6 i t 0 ω g 3 ω P o l y G a m m a [ 0 , 3 4 + 3 i ω 4 ] = 0 .
where P o l y G a m m a [ n , z ] gives the n t h derivative of the digamma function ψ ( n ) ( z ) , i.e., ψ ( 0 ) ( z ) = Γ ( z ) Γ ( z ) and Γ ( z ) = 0 t z 1 e t d t (see, for example, [36]).
For example, the equation M ( t 0 ) = 0 (for p = 2 ; N = 3 ; A = 1.01 ; ω = 0.35 ;   g 1 = 0.5 ;   g 2 = 0.95 ; g 3 = 0.8 ) is depicted in Figure 14.
Of course, the user can make some simplifications in Equation (12) himself. But naturally, the question arises as to why, for example, the existing HarmonicNumber[n] operator (in the specified programming environment) was not used (automatically) to obtain a compact record. The examples we can cite are many and varied. This necessitates the professional upgrading of the existing specialized modules.

8. Concluding Remarks

We have investigated in this paper a generalized differential model for micro-electromechanical oscillators in the light of Melnikov’s approach. The derived results can be used to estimate the associated total energy potential of the considered differential system. Nonstandard numerical methods connected to deriving the roots of the nonlinear equation M ( t 0 ) = 0 can be found in [37,38,39]. In some cases (when generating chaos in dynamical systems via x | x | ), M ( t ) is of polynomial type. Numerical methods for solving polynomial equations can be found in [40,41,42,43,44]. Some dynamic modules implemented in CAS Wolfram Mathematica for investigating the dynamics of some new hypothetical oscillators have been demonstrated. A cloud version that requires only a browser and internet connection is offered for some of them. They will be an integral part of the above-mentioned Web-based application for scientific computing. Let us note that the theoretical apparatus for studying the circuit implementation (design, fabricating, etc.) of the considered differential model is extremely complex and requires a serious investigation before being adapted for its possible inclusion in our planned Web platform. Other interesting problems are posed in the article [45]. Questions related to the computer and hardware modeling of periodically forced oscillators discussed in Section 6, especially for large values of N, can be considered open.
In a number of cases, Melnikov’s function M ( t ) can be used to approximate rectangular-type signals (see Figure 15) as well as to synthesize digital filters. Melnikov’s function M ( t ) as a typical emitting diagram is depicted in Figure 16. For some details, see [46]. Such reviews will also be the subject of our planned Web application.
Our planned platform implements the following: the Melnikov functions of higher kind for analysis of homo/heteroclinic bifurcation in mechanical systems; bifurcation diagrams; and chaos in a nonautonomous system. Some of these algorithms are based on known classical and more recent research [6,7,8,9,10,11,12,13,20]. These modules are improving on similar ones realized in computer algebraic systems designed for scientific calculations.
We fully understand that the construction of such an ambitious Web-based platform for scientific computing can only be realized with the active participation of specialists from various branches of scientific knowledge.
Future investigations. We have plans to explore modified escape oscillators of the form:
d x d t = y d y d t = P o l y ( x ) ϵ A p = 1 n y | y | p 1 + j = 1 N g j sin ( j ω t ) ,
where 0 ϵ < 1 , A is the damping level, p 1 is the damping exponent, N is the integer and P o l y ( x ) = { P ( x ) , H ( x ) , T ( x ) , U ( x ) } where P ( x ) , H ( x ) , T ( x ) , U ( x ) are Legendre, Hermite, Chebyshev (first kind) and Chebyshev (second kind) polynomials, respectively.
We envisage the study of such interesting dynamic models to be implemented in the Web application cited above.

Author Contributions

Conceptualization, N.K. and T.Z.; methodology, N.K. and T.Z.; software, V.K., A.I. and A.R.; validation, N.K., A.R., T.Z. and A.I.; formal analysis, N.K. and A.I.; investigation, T.Z., V.K., A.R. and N.K.; resources, A.I., T.Z., V.K. and N.K.; data curation, N.K. and V.K.; writing—original draft preparation, V.K., N.K. and A.I.; writing—review and editing, A.R.; visualization, V.K.; supervision, N.K.; project administration, N.K.; and funding acquisition, A.R., T.Z. and A.I. All authors have read and agreed to the published version of the manuscript.

Funding

The first, third, fourth and fifth authors are supported by the European Union-NextGenerationEU through the National Plan for Recovery and Resilience of the Republic Bulgaria, project No BG-RRP-2.004-0001-C01. The second author was financed by the European Union-NextGenerationEU through the National Recovery and Resilience Plan of the Republic of Bulgaria, project No BG-RRP-2.004-0008.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zounes, R.; Rand, R. Subharmonic resonance in the nonlinear Mathieu equation. Int. J. Non-Linear Mech. 2002, 37, 43–73. [Google Scholar] [CrossRef]
  2. DeMartini, B.; Butterfield, H.; Moehlis, J.; Turner, K. Chaos for a microelectromechanical oscillator governed by the nonlinear Mathieu equation. J. Microelectromech. Syst. 2007, 16, 1314–1323. [Google Scholar] [CrossRef]
  3. DeMartini, B.; Moehlis, J.; Turner, K.; Rhoads, J.; Shaw, S.; Zhang, W. Modelling of parametrically excited microelectromechanical oscillator dynamics with application to filtering. In Proceedings of the IEEE Sensors Conference, Irvine, CA, USA, 30 October–3 November 2005; pp. 345–348. [Google Scholar]
  4. Luo, A.; Wang, F. Nonlinear dynamics of a micro-electro-mechanical system with time-varying capacitors. J. Vib. Acoust. 2004, 126, 77–83. [Google Scholar] [CrossRef]
  5. Wiggins, S. Introduction to Applied Nonlinear Dynamical Systems and Chaos, 2nd ed.; Springer: New York, NY, USA, 2003. [Google Scholar]
  6. Lenci, S.; Lupini, R. Homoclinic and heteroclinic solutions for a class of two-dimensional Hamiltonian systems. Z. Angew. Math. Phys. 1996, 47, 97–111. [Google Scholar] [CrossRef]
  7. Lenci, S.; Rega, G. Controlling nonlinear dynamics in a two-well impact system. I. Attractors and bifurcation scenario under symmetric excitations. Int. J. Bifurc. Chaos 1998, 8, 2387–2407. [Google Scholar] [CrossRef]
  8. Lenci, S.; Rega, G. Controlling nonlinear dynamics in a two-well impact system. II. Attractors and bifurcation scenario under unsymmetric optimal excitation. Int. J. Bifurc. Chaos 1998, 8, 2409–2424. [Google Scholar] [CrossRef]
  9. Lenci, S.; Rega, G. Higher-order Melnikov analysis of homo/heteroclinic bifurcations in mechanical oscillators. In Proceedings of the 16th AIMETA Congress of Theoretical and Applied Mechanics (Ferrara), Ferrara, Italy, 9–12 September 2003. CD-Rom. [Google Scholar]
  10. Lenci, S.; Rega, G. Optimal control of homoclinic bifurcation: Theoretical treatment and practical reduction of safe basin erosion in the Helmholtz oscillator. J. Vib. Control 2003, 9, 281–315. [Google Scholar] [CrossRef]
  11. Lenci, S.; Rega, G. Bifurcations and chaos in single-d.o.f. mechanical systems: Exploiting nonlinear dynamics properties for their control. In Recent Research Developments in Structural Dynamics; Luongo, A., Ed.; Research Signpost: Kerala, India, 2003; pp. 331–369. [Google Scholar]
  12. Soto-Trevino, C.; Kaper, T. Higher-order Melnikov theory for adiabatic systems. J. Math. Phys. 1996, 37, 6220–6249. [Google Scholar] [CrossRef]
  13. Thompson, M.; Stewart, H. Nonlinear Dynamics and Chaos, 2nd ed.; John Wiley & Sons: Chichester, UK, 2002. [Google Scholar]
  14. Kyurkchiev, N.; Zaevski, T. On a hypothetical oscillator: Investigations in the light of Melnikov’s approach, some simulations. Int. J. Differ. Equations Appl. 2023, 22, 67–79. [Google Scholar]
  15. Melnikov, V.K. On the stability of a center for time–periodic perturbation. Trans. Mosc. Math. Soc. 1963, 12, 1–57. [Google Scholar]
  16. Kyurkchiev, V.; Iliev, A.; Rahnev, A.; Kyurkchiev, N. On a class of orthogonal polynomials as corrections in Lienard differential system. Applications. Algorithms 2023, 16, 297. [Google Scholar] [CrossRef]
  17. Kyurkchiev, N.; Iliev, A. On the hypothetical oscillator model with second kind Chebyshev’s polynomial–correction: Number and type of limit cycles, simulations and possible applications. Algorithms 2022, 15, 462. [Google Scholar] [CrossRef]
  18. Vasileva, M.; Kyurkchiev, V.; Iliev, A.; Rahnev, A.; Kyurkchiev, N. Applications of Some Orthogonal Polynomials and Associated Polynomials of Higher Kind; Plovdiv University Press: Plovdiv, Bulgaria, 2023; ISBN 978-619-7663-72-3. [Google Scholar]
  19. Blows, T.; Perko, L. Bifurcation of Limit Cycles from Centers and Separatrix Cycles of Planar Analytic Systems. SIAM Rev. 1994, 36, 341–376. [Google Scholar] [CrossRef]
  20. Gavrilov, L.; Iliev, I. The limit cycles in a generalized Rayleigh–Lienard oscillator. Discret. Contin. Dyn. Syst. 2023, 43, 2381–2400. [Google Scholar] [CrossRef]
  21. Li, X.; Shen, Y.; Sun, J.; Yang, S. New periodic attractors is slow-fast Duffing system with periodic parametric excitation. Sci. Rep. 2019, 9, 11185. [Google Scholar] [CrossRef]
  22. Holmes, P.J. Chaos in Duffing’s equation: Strange Attractors? In Encyclopedia of Vibration; Braun, S.G., Ewins, D.J., Rao, S.S., Eds.; Academic Press: London, UK; San Diego, CA, USA, 2001; pp. 227–236. [Google Scholar]
  23. Sanchez, N.; Nayfeh, A. Prediction of bifurcations in a parametrically excited Duffing oscillator. Int. J. Nonlinear Mech. 1990, 25, 163–176. [Google Scholar] [CrossRef]
  24. Michon, G.; Manin, L.; Parker, R.; Dufour, R. Duffing oscillator with parametric excitation: Analytical and experimental investigation on a Belt–Pulley systems. ASME 2008, 3, 031001. [Google Scholar] [CrossRef]
  25. Yue, X.; Lv, G.; Zhang, Y. Rare and hidden attractors in a periodically forced Duffing system with absolute nonlinearity. Chaos Solitons Fractals 2021, 150, 111108. [Google Scholar] [CrossRef]
  26. Tang, K.; Man, K.; Zhong, G.; Chen, G. Generating chaos via x|x|. IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 2001, 48, 636–641. [Google Scholar] [CrossRef]
  27. Guckenheimer, J.; Holmes, P. Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields; Sprinder: New York, NY, USA, 1983. [Google Scholar]
  28. Perko, L. Differential Equations and Dynamical Systems; Springer: New York, NY, USA, 1991. [Google Scholar]
  29. Holmes, P.; Marsden, J. Horseshoes in perturbation of Hamiltonian systems with two degrees of freedom. Comm. Math. Phys. 1982, 82, 523–544. [Google Scholar] [CrossRef]
  30. Holmes, P.; Marsden, J. A partial differential equation with infinitely many periodic orbits: Chaotic oscillations. Arch. Ration. Mech. Anal. 1981, 76, 135–165. [Google Scholar] [CrossRef]
  31. Gavrilov, L.; Iliev, I. Complete hyperelliptic integrals of the first kind and their non–oscillation. Trans. Am. Math. Soc. 2003, 356, 1185–1207. [Google Scholar] [CrossRef]
  32. Siewe, M.; Kakmeni, F.; Tchawoua, C. Resonant oscillation and homoclinic bifurcation in ϕ6—Van der Pol oscillator. Chaos Solut. Fractals 2004, 21, 841–853. [Google Scholar] [CrossRef]
  33. Adelakun, A.; Njah, A.; Olusola, O.; Wara, S. Computer and hardware modeling of periodically forced ϕ6—Van der Pol oscillator. Act. Passiv. Electron. Components 2016, 7, 3426713. [Google Scholar]
  34. Yu, J.; Li, J. Investigation on dynamics of the extended Duffing–Van der Pol system. Zetschrift Fur Naturforschung 2014, 64, 341–346. [Google Scholar] [CrossRef]
  35. Siewe, M.; Tchawoua, C.; Rajasekar, S. Homoclinic bifurcation and chaos in ϕ6 Rayleigh oscillator with three wells driven by an amplitude modulated force. Int. J. Bifurc. Chaos 2011, 21, 1583–1593. [Google Scholar] [CrossRef]
  36. Abramowicz, A.; Stegun, I. Handbook of Mathematical Functions; Dover: New York, NY, USA, 1970. [Google Scholar]
  37. Thangkhenpau, G.; Panday, S.; Bolunduet, L.C.; Jantschi, L. Efficient Families of Multi-Point Iterative Methods and Their Self-Acceleration with Memory for Solving Nonlinear Equations. Symmetry 2023, 15, 1546. [Google Scholar] [CrossRef]
  38. Akram, S.; Khalid, M.; Junjua, M.-U.-D.; Altaf, S.; Kumar, S. Extension of King’s Iterative Scheme by Means of Memory for Nonlinear Equations. Symmetry 2023, 15, 1116. [Google Scholar] [CrossRef]
  39. Cordero, A.; Janjua, M.; Torregrosa, J.R.; Yasmin, N.; Zafar, F. Efficient four parametric with and without-memory iterative methods possessing high efficiency indices. Math. Probl. Eng. 2018, 1–12. [Google Scholar] [CrossRef]
  40. Proinov, P.; Vasileva, M. Local and Semilocal Convergence of Nourein’s Iterative Method for Finding All Zeros of a Polynomial Simultaneously. Symmetry 2020, 12, 1801. [Google Scholar] [CrossRef]
  41. Proinov, P.; Ivanov, S. On the convergence of Halley’s method for simultaneous computation of polynomial zeros. J. Numer. Math. 2015, 23, 379–394. [Google Scholar] [CrossRef]
  42. Proinov, P.; Vasileva, M. On the convergence of high-order Gargantini-Farmer-Loizou type iterative methods for simultaneous approximation of polynomial zeros. Appl. Math. Comput. 2019, 361, 202–214. [Google Scholar] [CrossRef]
  43. Ivanov, S. Families of high-order simultaneous methods with several corrections. Numer. Algor. 2024. [Google Scholar] [CrossRef]
  44. Kyncheva, V.; Yotov, V.; Ivanov, S. Convergence of Newton, Halley and Chebyshev iterative methods as methods for simultaneous determination of multiple zeros. Appl. Numer. Math. 2017, 112, 146–154. [Google Scholar] [CrossRef]
  45. Stojanovic, V.; Petkovic, M.; Deng, J. Instability of vehicle systems moving along an infinite beam on a viscoelastic foundation. Eur. J. Mech.-A/Solids 2018, 69, 238–254. [Google Scholar] [CrossRef]
  46. Kyurkchiev, N.; Andreev, A. Approximation and Antenna and Filter Synthesis: Some Moduli in Programming Environment MATHEMATICA; LAP Lambert Academic Publishing: Saarbrucken, Germany, 2014; ISBN 978-3-859-53322-8. [Google Scholar]
Figure 1. Double homoclinic orbit [27].
Figure 1. Double homoclinic orbit [27].
Symmetry 16 00253 g001
Figure 2. Potential energy V ( x ) .
Figure 2. Potential energy V ( x ) .
Symmetry 16 00253 g002
Figure 3. The nonlinear Equation (6): (a) for α = 1.39627 , a 0 = 1.2 , a = 1 , ω = 2.2 , g 1 = 0.7 , M ( t 0 ) has root t 0 0.0003209 with multiplicity two; (b) for α = 1.39627 , a 0 = 1.15 , a = 1 , ω = 2.5 , q 1 = 0.6 , M ( t 0 ) has no roots.
Figure 3. The nonlinear Equation (6): (a) for α = 1.39627 , a 0 = 1.2 , a = 1 , ω = 2.2 , g 1 = 0.7 , M ( t 0 ) has root t 0 0.0003209 with multiplicity two; (b) for α = 1.39627 , a 0 = 1.15 , a = 1 , ω = 2.5 , q 1 = 0.6 , M ( t 0 ) has no roots.
Symmetry 16 00253 g003
Figure 4. The nonlinear Equation (8): (a) for α = 1.7 , a 0 = 0.8 , a = 1 , ω = 0.75 , g 1 = 0.7 , g 2 = 0.3 , g 3 = 0.1 , M ( t 0 ) has root t 0 0.144 ; (b) for α = 1.95 , a 0 = 0.7 , a = 1 , ω = 0.6 , g 1 = 0.8 , g 2 = 0.4 , g 3 = 0.15 , M ( t 0 ) has no roots.
Figure 4. The nonlinear Equation (8): (a) for α = 1.7 , a 0 = 0.8 , a = 1 , ω = 0.75 , g 1 = 0.7 , g 2 = 0.3 , g 3 = 0.1 , M ( t 0 ) has root t 0 0.144 ; (b) for α = 1.95 , a 0 = 0.7 , a = 1 , ω = 0.6 , g 1 = 0.8 , g 2 = 0.4 , g 3 = 0.15 , M ( t 0 ) has no roots.
Symmetry 16 00253 g004
Figure 5. (a) x-time series; (b) y-time series; (c) phase space; (Example 1).
Figure 5. (a) x-time series; (b) y-time series; (c) phase space; (Example 1).
Symmetry 16 00253 g005
Figure 6. (a) x-time series; (b) y-time series; (c) phase space; (Example 2).
Figure 6. (a) x-time series; (b) y-time series; (c) phase space; (Example 2).
Symmetry 16 00253 g006
Figure 7. (a) x-time series; (b) y-time series; (c) phase space; (Example 3).
Figure 7. (a) x-time series; (b) y-time series; (c) phase space; (Example 3).
Symmetry 16 00253 g007
Figure 8. Possible restrictions.
Figure 8. Possible restrictions.
Symmetry 16 00253 g008
Figure 9. (a) y-time series (with restriction); (b) phase space; (Example 3).
Figure 9. (a) y-time series (with restriction); (b) phase space; (Example 3).
Symmetry 16 00253 g009
Figure 10. Potential energy V ( x ) for b = 1 , b 1 = 0.7 , b 0 = 0.1 .
Figure 10. Potential energy V ( x ) for b = 1 , b 1 = 0.7 , b 0 = 0.1 .
Symmetry 16 00253 g010
Figure 11. For b = 0.4 , b 1 = 0.7 , b 0 = 0.1 (a) homoclinic orbits; (b) heteroclinic orbits.
Figure 11. For b = 0.4 , b 1 = 0.7 , b 0 = 0.1 (a) homoclinic orbits; (b) heteroclinic orbits.
Symmetry 16 00253 g011
Figure 12. (a) x-time series; (b) y-time series; (c) phase space; (Example 4).
Figure 12. (a) x-time series; (b) y-time series; (c) phase space; (Example 4).
Symmetry 16 00253 g012
Figure 13. (a) x-time series; (b) y-time series; (c) phase space; (Example 5).
Figure 13. (a) x-time series; (b) y-time series; (c) phase space; (Example 5).
Symmetry 16 00253 g013
Figure 14. The equation M ( t 0 ) = 0 (Proposition 5.) for fixed parameters: (a) p = 2 , N = 3 , A = 1.01 , ω = 0.35 , g 1 = 0.5 , g 2 = 0.95 , g 3 = 0.8 ; (b) A = 0.1 , ω = 0.9 , g 1 = 0.01 , g 2 = 0.02 , g 3 = 0.015 , M ( t 0 ) has no roots.
Figure 14. The equation M ( t 0 ) = 0 (Proposition 5.) for fixed parameters: (a) p = 2 , N = 3 , A = 1.01 , ω = 0.35 , g 1 = 0.5 , g 2 = 0.95 , g 3 = 0.8 ; (b) A = 0.1 , ω = 0.9 , g 1 = 0.01 , g 2 = 0.02 , g 3 = 0.015 , M ( t 0 ) has no roots.
Symmetry 16 00253 g014
Figure 15. Melnikov function M ( t ) (thick) for α = 1.95 , a 0 = 0.76 , a = 1 , ω = 0.6 , g 1 = 0.8 , g 2 = 0.4 , g 3 = 0.15 in interval [0,10] (see Proposition 3).
Figure 15. Melnikov function M ( t ) (thick) for α = 1.95 , a 0 = 0.76 , a = 1 , ω = 0.6 , g 1 = 0.8 , g 2 = 0.4 , g 3 = 0.15 in interval [0,10] (see Proposition 3).
Symmetry 16 00253 g015
Figure 16. Melnikov function M ( t ) for A = 0.011 , ω = 0.9 , g 1 = 0.01 , g 2 = 0.02 , g 3 = 0.03 in interval [ 5 , 5 ] (see Proposition 4).
Figure 16. Melnikov function M ( t ) for A = 0.011 , ω = 0.9 , g 1 = 0.01 , g 2 = 0.02 , g 3 = 0.03 in interval [ 5 , 5 ] (see Proposition 4).
Symmetry 16 00253 g016
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kyurkchiev, N.; Zaevski, T.; Iliev, A.; Kyurkchiev, V.; Rahnev, A. Nonlinear Dynamics of a New Class of Micro-Electromechanical Oscillators—Open Problems. Symmetry 2024, 16, 253. https://doi.org/10.3390/sym16020253

AMA Style

Kyurkchiev N, Zaevski T, Iliev A, Kyurkchiev V, Rahnev A. Nonlinear Dynamics of a New Class of Micro-Electromechanical Oscillators—Open Problems. Symmetry. 2024; 16(2):253. https://doi.org/10.3390/sym16020253

Chicago/Turabian Style

Kyurkchiev, Nikolay, Tsvetelin Zaevski, Anton Iliev, Vesselin Kyurkchiev, and Asen Rahnev. 2024. "Nonlinear Dynamics of a New Class of Micro-Electromechanical Oscillators—Open Problems" Symmetry 16, no. 2: 253. https://doi.org/10.3390/sym16020253

APA Style

Kyurkchiev, N., Zaevski, T., Iliev, A., Kyurkchiev, V., & Rahnev, A. (2024). Nonlinear Dynamics of a New Class of Micro-Electromechanical Oscillators—Open Problems. Symmetry, 16(2), 253. https://doi.org/10.3390/sym16020253

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