Next Article in Journal
A New Efficient Method for the Numerical Solution of Linear Time-Dependent Partial Differential Equations
Next Article in Special Issue
New Bell–Sheffer Polynomial Sets
Previous Article in Journal
Characterizations of the Total Space (Indefinite Trans-Sasakian Manifolds) Admitting a Semi-Symmetric Metric Connection
Previous Article in Special Issue
Solutions to Abel’s Integral Equations in Distributions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Periodically Forced Nonlinear Oscillatory Acoustic Vacuum

by
Makrina Agaoglou
1,2,
Michal Fečkan
1,3,*,
Michal Pospíšil
1,3,
Vassilis M. Rothos
2 and
Alexander F. Vakakis
4
1
Mathematical Institute of Slovak Academy of Sciences, Štefánikova 49, 814 73 Bratislava, Slovakia
2
Department of Mechanical Engineering, Faculty of Engineering, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
3
Department of Mathematical Analysis and Numerical Mathematics, Comenius University in Bratislava, Mlynská dolina, 842 48 Bratislava, Slovakia
4
Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
*
Author to whom correspondence should be addressed.
Axioms 2018, 7(4), 69; https://doi.org/10.3390/axioms7040069
Submission received: 31 July 2018 / Revised: 14 September 2018 / Accepted: 19 September 2018 / Published: 22 September 2018
(This article belongs to the Special Issue Mathematical Analysis and Applications)

Abstract

:
In this work, we study the in-plane oscillations of a finite lattice of particles coupled by linear springs under distributed harmonic excitation. Melnikov-type analysis is applied for the persistence of periodic oscillations of a reduced system.

1. Introduction

We analytically study the persistence of periodic oscillations for certain three-dimensional systems of ordinary differential equations (ODEs) with periodic perturbations and a slowly-varying variable. The considered ODEs are derived from a model of a finite lattice of particles coupled by linear springs under distributed harmonic excitation, which is described in detail in Section 2. This model presents a low-energy nonlinear acoustic vacuum. We refer the reader for more motivations, further details and applications to [1,2]. Following the computations of [1], we consider just two modes in Section 3, and we postpone higher modes investigation to our future paper, since another approach will be used. Melnikov analysis is demonstrated in Section 4 for finding conditions for the existence of periodic solutions for the perturbed ODEs corresponding to two modes. More precisely, following [1], we derive a three-dimensional periodically-perturbed system of ODEs with a slowly-varying variable. Then, we analyze an unperturbed autonomous system of ODEs to compute its family of periodic solutions by revising the results of [1] in more detail. Since we are interested in the persistence of periodic solutions for the perturbed ODEs, we compute the corresponding Melnikov functions by [3]. Due to the difficulty of finding simple roots of these Melnikov functions explicitly, we outline an asymptotic approach for the location of some of them. Note that the simple roots of Melnikov functions predict the persistence and location of periodic solutions for perturbed ODEs. This is a novelty and a contribution of our paper. Section 5 outlines possible future research along with summarizing our achievements in this paper.

2. The Model

We consider a lattice consisting of N identical particles coupled by identical linear springs (they are un-stretched when the lattice is in the horizontal position) and executing in-plane oscillations (see Figure 1). Fixed boundary conditions and dissipative terms are imposed, and the transverse harmonic forces are also applied. The equations of motion can be expressed as follows,
m d 2 u i d t 2 + T i ξ d ϵ i d t cos ϕ i T i + 1 ξ d ϵ i + 1 d t cos ϕ i + 1 = 0 m d 2 v i d t 2 + T i ξ d ϵ i d t sin ϕ i T i + 1 ξ d ϵ i + 1 d t sin ϕ i + 1 F i = 0
with u i , v i being the longitudinal and transversal displacements of the i-th particle, respectively, ϕ i the angle between the i-th spring and the horizontal direction, ξ the damping coefficient, ϵ i = l i l i the deformation of the i-th spring, F i the exciting transverse force and m the mass of each particle of the lattice. The tensile forces are proportional to the deformations of the springs, and considering the geometry of the deformed state of the lattice (see Figure 1), one may write:
T i = k ( l i l i ) , ϵ i = l i l i = [ ( v i v i 1 ) 2 + ( l i + u i u i 1 ) 2 ] 1 / 2 l i
with l i being the equilibrium length of the i-th spring (each spring has the same length) and k the linear stiffness coefficient of each coupling spring. Introducing δ i = ϵ i / l i , s i = u i / l i , w i = v i / l i , c = ξ / ( k m ) 1 / 2 , where s i and w i are the normalized axial and transverse displacements, and the “slow” time scale τ = ( k m ) 1 / 2 t , Equation (1) can be rewritten in normalized form:
d 2 s i d τ 2 = δ i cos ϕ i δ i 1 cos ϕ i 1 + c δ i cos ϕ i c δ i 1 cos ϕ i 1 d 2 w i d τ 2 = δ i sin ϕ i δ i 1 sin ϕ i 1 + c δ i sin ϕ i c δ i 1 sin ϕ i 1 + f i ,
where:
δ i = [ ( w i + 1 w i ) 2 + ( 1 + s i + 1 s i ) ] 1 / 2 1 , δ i = ( w i + 1 w i ) ( w i + 1 w i ) + ( 1 + s i + 1 s i ) ( s i + 1 s i ) [ ( w i + 1 w i ) 2 + ( 1 + s i + 1 s i ) 2 ] 1 / 2
and:
cos ϕ i = 1 + s i + 1 s i [ ( w i + 1 w i ) 2 + ( 1 + s i + 1 s i ) 2 ] 1 / 2 , sin ϕ i = w i + 1 w i [ ( w i + 1 w i ) 2 + ( 1 + s i + 1 s i ) 2 ] 1 / 2 , f i = F i k .
The normalized system (2) is referred to as the “exact lattice” in the following sections.
According to the previous research [1], when we introduce this system (2) without extra transverse force and damping terms, an interesting feature is that in the low energy limit and under the assumption that the axial displacements are assumed to be an order of magnitude smaller compared to the transverse ones, it was shown that, correct to the leading order of approximation, the transverse oscillations decouple from the axial ones and are governed by the following reduced system of equations for predominantly transverse oscillations of the particles:
d 2 w i d τ 2 + 2 1 ( N + 1 ) 1 q = 0 N ( w q + 1 w q ) 2 ( 2 w i w i + 1 w i 1 ) = 0 , i = 1 , , N , w 0 ( 0 ) = w N + 1 ( 0 ) = 0 .
Then, the nearly linear axial oscillations are driven by the transverse responses (see [1,2] for more details). Therefore, we focus our analysis like in [1] just on Equation (3), which presents a low-energy nonlinear acoustic vacuum, because in the absence of linear terms, it possesses zero speed of sound in the context of classical linear acoustics. What is more, it is notable that the existence of the strongly nonlocal multiplier 2 1 ( N + 1 ) 1 q = 0 N ( w q + 1 w q ) 2 indicates that the response of each particle is dependent (and hence, it is coupled) on the responses of all other particles. Equation (3) admits N exact nonlinear standing waves, or nonlinear normal modes (NNMs), in the form:
w i ( τ ) = A p ( τ ) sin π p i N + 1 , i = 1 , , N
for the p-th NNM, 1 p N , where A p ( τ ) denotes the p-th modal amplitude. These, by construction, are mutually orthogonal, and there are no other NNMs in this system, nor any NNM bifurcations [1]. Substituting this NNM ansatz into Equation (3) yields a set of N uncoupled nonlinear equations governing the time-dependent amplitudes of the NNMs:
A p ( τ ) + ( 1 / 4 ) ω p 4 A p 3 ( τ ) = 0
with:
ω p = 2 sin π p 2 ( N + 1 ) ,
which is the p-th natural frequency of the corresponding linear system Equation (3) and the prime denoting differentiation with respect to τ . The exciting force, which is applied on each particle in the transverse direction, is expressed as:
f i = F p cos ω p τ sin p i π n + 1
where i = 1 , , N , for the p-th NNM, 1 p N , and this exciting force includes NNMs in the form sin p i π n + 1 , i = 1 , , N , for the p-th NNM, 1 p N and the p-th natural linear frequency ω p .
The frequency of the p-th NNM is tunable with the force and energy, and it also paves the way for nonlinear resonances between NNMs widely separated in the nonlinear spectrum, given that their energies tune their frequencies to satisfy certain rational relationships.
Summarizing, we can write (3) as:
d 2 w d τ 2 + 2 1 ( N + 1 ) 1 M w , w M w = 0 ,
where w = [ w 1 , , w N ] R N , M is a symmetric matrix given by
M = 2 1 0 0 1 2 1 0 0 1 2 1 0 0 1 2
and · , · is the standard scalar product on R N . The eigenvectors of M are ϕ ̲ p = [ sin p π N + 1 , , sin p N π N + 1 ] with the corresponding eigenvalues ω p 2 , 1 p N . Moreover, it holds (see [4], p. 37):
ϕ ̲ p , ϕ ̲ p = i = 1 N sin 2 p i π N + 1 = N + 1 2 , ϕ ̲ p , ϕ ̲ k = i = 1 N sin p i π N + 1 sin k i π N + 1 = 1 2 i = 1 N cos ( p k ) i π N + 1 cos ( p + k ) i π N + 1 = 0 , p k .
The forced (3) has the form:
d 2 w d τ 2 + 2 1 ( N + 1 ) 1 M w , w M w = p = 1 N F p cos ω p τ ϕ ̲ p .
Therefore, considering the basis { ϕ ̲ p } p = 1 N of R N , we take w ( τ ) = p = 1 N C p ( τ ) ϕ ̲ p in (4) to get:
C p ( τ ) + 1 4 i = 1 N C i 2 ( τ ) ω i 2 ω p 2 C p ( τ ) = F p cos ω p τ , 1 p N .
Next, applying the coordinate transformation to (5):
A p ( τ ) = ω p 2 C p ( τ ) , 1 p N ,
we get:
A p ( τ ) + i = 1 N A i 2 ( τ ) ω p 2 A p ( τ ) = F p ω p 2 cos ω p τ , 1 p N .

3. Two-Mode System

In this paper, we consider just two modes: k and p, so we study the system:
A k ( τ ) + [ A k 2 ( τ ) + A p 2 ( τ ) ] ω k 2 A k ( τ ) + ϵ μ 1 cos ( ω k τ ) = 0 A p ( τ ) + [ A k 2 ( τ ) + A p 2 ( τ ) ] ω p 2 A p ( τ ) + ϵ μ 2 cos ( ω p τ ) = 0 ,
for ϵ 0 small and parameters μ 1 , μ 2 . Using the transformation:
ψ 1 ( τ ) = A k ( τ ) + j Ω A k ( τ ) ζ 1 ( τ ) e j Ω τ ψ 2 ( τ ) = A p ( τ ) + j Ω A p ( τ ) ζ 2 ( τ ) e j Ω τ ,
and performing an averaging approach like in [1], Equation (6) is modified to the form:
ζ 1 + j Ω 2 ζ 1 j ω k 2 8 Ω 3 ( ζ 2 2 ζ 1 + 2 | ζ 2 | 2 ζ 1 + 3 | ζ 1 | 2 ζ 1 ) + ϵ μ 1 cos ( ω k τ ) e j Ω τ = 0 ζ 2 + j Ω 2 ζ 2 j ω p 2 8 Ω 3 ( ζ 1 2 ζ 2 + 2 | ζ 1 | 2 ζ 2 + 3 | ζ 2 | 2 ζ 2 ) + ϵ μ 2 cos ( ω p τ ) e j Ω τ = 0 .
Introducing ζ i = a i e j β i and Δ = β 2 β 1 , we get:
a 1 + ω k 2 8 Ω 3 a 2 2 a 1 sin 2 Δ + ϵ μ 1 cos ( ω k τ ) cos ( Ω τ + β 1 ) = 0 a 2 ω p 2 8 Ω 3 a 1 2 a 2 sin 2 Δ + ϵ μ 2 cos ( ω p τ ) cos ( Ω τ + β 1 + Δ ) = 0 Δ ω p 2 8 Ω 3 ( 3 a 2 2 + a 1 2 cos 2 Δ + 2 a 1 2 ) + ω k 2 8 Ω 3 ( 3 a 1 2 + a 2 2 cos 2 Δ + 2 a 2 2 ) ϵ μ 2 a 2 cos ( ω p τ ) sin ( Ω τ + β 1 + Δ ) + ϵ μ 1 a 1 cos ( ω k τ ) sin ( Ω τ + β 1 ) = 0
where we consider β 1 as a constant parameter. Now, by introducing the coordinate transformations a 1 = ( ρ ω p ) sin θ and a 2 = ( ρ ω k ) cos θ into Equation (7), we get:
ρ = ϵ μ 1 ω p sin θ cos ( ω k τ ) cos ( Ω τ + β 1 ) ϵ μ 2 ω k cos θ cos ( ω p τ ) cos ( Ω τ + Δ + β 1 ) θ + ρ 2 16 Ω 3 sin 2 θ sin 2 Δ + ϵ μ 1 ω p cos θ ρ cos ( ω k τ ) cos ( Ω τ + β 1 ) ϵ μ 2 ω k sin θ ρ cos ( ω p τ ) cos ( Ω τ + Δ + β 1 ) = 0 Δ ρ 2 8 Ω 3 3 ω p 2 ω k 2 cos 2 θ 3 ω k 2 ω p 2 sin 2 θ cos 2 θ ( 2 + cos 2 Δ ) ϵ μ 2 ω k ρ cos θ cos ( ω p τ ) sin ( Ω τ + Δ + β 1 ) + ϵ μ 1 ω p ρ sin θ cos ( ω k τ ) sin ( Ω τ + β 1 ) = 0 .
In the rest of the paper, we assume ω p = ω k = P and Ω = k P for a natural number k, so we study the periodically-perturbed system:
ρ = ϵ μ 1 sin θ cos ( P τ ) cos ( k P τ + β 1 ) ϵ μ 2 cos θ cos ( P τ ) cos ( k P τ + Δ + β 1 ) θ + ρ 2 16 k 3 P 3 sin 2 θ sin 2 Δ + ϵ μ 1 cos θ ρ cos ( P τ ) cos ( k P τ + β 1 ) ϵ μ 2 sin θ ρ cos ( P τ ) cos ( k P τ + Δ + β 1 ) = 0 Δ ρ 2 4 k 3 P 3 cos 2 θ sin 2 Δ + ϵ μ 1 ρ sin θ cos ( P τ ) sin ( k P τ + β 1 ) ϵ μ 2 ρ cos θ cos ( P τ ) sin ( k P τ + Δ + β 1 ) = 0
where we scaled P μ i μ i , i = 1 , 2 . We may suppose:
μ 1 2 + μ 2 2 = 1 .
First, consider the unperturbed case where ϵ = 0 , so the system:
ρ = 0 θ + ρ 2 16 k 3 P 3 sin 2 θ sin 2 Δ = 0 Δ ρ 2 4 k 3 P 3 cos 2 θ sin 2 Δ = 0 .
By introducing the temporal variable τ 2 = ρ 2 8 k 3 P 3 τ in Equation (9), we get:
d θ d τ 2 = 1 2 sin 2 θ sin 2 Δ d Δ d τ 2 = 2 cos 2 θ sin 2 Δ
which is fully integrable and gives us the first integral I = sin 2 θ sin Δ = K = c o n s t . of the degenerate slow flow. If we consider the following initial conditions: θ ( 0 ) = θ 0 , where 0 < θ 0 < π 4 and Δ ( 0 ) = π / 2 , then we get K = sin 2 θ 0 ( 0 , 1 ) . To find exact solutions of (10), we first derive:
d τ 2 d θ = 2 sin 2 θ sin 2 Δ = 1 K cos Δ = 1 K 1 sin 2 Δ = 1 K 1 K 2 sin 2 2 θ = 1 K sin 2 θ sin 2 2 θ K 2 ,
for τ 2 > 0 small, since (10) gives 0 < θ < π 2 , 0 < Δ < π (see Figure 2) and d Δ d τ 2 ( 0 ) = 2 cos 2 θ 0 > 0 , so for τ 2 > 0 small, we have Δ ( τ 2 ) > π 2 .
By using the formula in [4] ((2.599.4) p. 205), we obtain:
τ 2 = θ 0 θ 1 K sin 2 θ sin 2 2 θ K 2 d θ = 1 2 K π 2 sin 1 cos 2 θ 1 K 2 ,
which gives:
θ ( τ 2 ) = 1 2 cos 1 1 K 2 cos ( 2 K τ 2 ) ,
recalling 0 < θ < π 2 . Of course, Formula (11) holds for all τ 2 , not just for small positive ones. Next, using (10) and (11), we have:
d Δ d τ 2 = 2 1 K 2 cos ( 2 K τ 2 ) sin 2 Δ ,
which can be easily solved to arrive at:
Δ ( τ 2 ) = π cot 1 1 K 2 K sin ( 2 K τ 2 ) .
Hence, the exact solution of the system (10) is given as:
θ ( τ 2 ) = 1 2 cos 1 1 K 2 cos ( 2 K τ 2 ) Δ ( τ 2 ) = π cot 1 1 K 2 K sin ( 2 K τ 2 )
where the period is:
T ( θ 0 ) = 2 K θ 0 π 2 θ 0 d θ 1 sin 2 2 θ 0 sin 2 2 θ 1 / 2 = π K .
It is easy to verify the following symmetry property (see Figure 3):
Δ τ 2 + T 2 = π Δ ( τ 2 ) , θ τ 2 + T 2 = π 2 θ ( τ 2 ) .
Summarizing, the exact solution of the unperturbed (9) is the following:
ρ ( τ ) = c o n s t . θ ( τ ) = 1 2 cos 1 1 K 2 cos K ρ 2 4 k 3 P 3 τ Δ ( τ ) = π cot 1 1 K 2 K sin K ρ 2 4 k 3 P 3 τ
with the period:
T = 8 k 3 P 3 π K ρ 2 .
Consequently, for any:
ρ > 2 k P 2 k ,
taking:
K ( ρ ) = 4 k 3 P 4 ρ 2 ( 0 , 1 ) ,
Equation (9) has the T = 2 π / P -periodic solution:
θ ( ρ , τ ) = 1 2 cos 1 ρ 4 16 k 6 P 8 ρ 2 cos ( P τ ) Δ ( ρ , τ ) = π cot 1 ρ 4 16 k 6 P 8 4 k 3 P 4 sin ( P τ ) .

4. Melnikov Analysis for Periodic Oscillations

Writing (8) as:
θ + ρ 2 16 k 3 P 3 sin 2 θ sin 2 Δ = ϵ ρ 1 g 1 ( ρ , θ , Δ , τ ) Δ ρ 2 4 k 3 P 3 cos 2 θ sin 2 Δ = ϵ ρ 1 g 2 ( ρ , θ , Δ , τ ) ρ = ϵ g 3 ( ρ , θ , Δ , τ )
and using [5], (3.5.11), p. 111, with α = 0 , and [6], Lemma 2.5, p. 283, we compute the Melnikov function:
M ( β 1 , ρ ) = M 1 ( β 1 , ρ ) , M 2 ( β 1 , ρ )
as:
M 1 ( β 1 , ρ ) = 0 T I θ ρ 1 g 1 + I Δ ρ 1 g 2 I θ θ ρ g 3 I Δ Δ ρ g 3 d τ = 0 T I θ ρ 1 g 1 + I Δ ρ 1 g 2 d τ d K d ρ 0 T g 3 d τ
since differentiating by ρ the identity:
I ( θ ( ρ , τ ) , Δ ( ρ , τ ) ) = K ( ρ ) ,
we get:
I θ θ ρ + I Δ Δ ρ = d K d ρ ,
which is independent of τ , and:
M 2 ( β 1 , ρ ) = 0 T g 3 d τ .
Formulas (13) and (14) are similar to [3], (2.7). We are looking for a simple zero of M, which is equivalent to considering:
M ¯ ( β 1 , ρ ) = M ¯ 1 ( β 1 , ρ ) , M ¯ 2 ( β 1 , ρ ) M ¯ 1 ( β 1 , ρ ) = 0 T I θ g 1 + I Δ g 2 d τ , M ¯ 2 ( β 1 , ρ ) = 0 T g 3 d τ .
Since:
g i = μ 1 g i 1 + μ 2 g i 2 , i = 1 , 2 , 3
for:
g 11 ( ρ , θ , Δ , τ ) = cos θ cos ( P τ ) cos ( k P τ + β 1 ) g 12 ( ρ , θ , Δ , τ ) = sin θ cos ( P τ ) cos ( k P τ + Δ + β 1 ) g 21 ( ρ , θ , Δ , τ ) = cos ( P τ ) sin θ sin ( k P τ + β 1 ) g 22 ( ρ , θ , Δ , τ ) = cos ( P τ ) cos θ sin ( k P τ + Δ + β 1 ) g 31 ( ρ , θ , Δ , τ ) = sin θ cos ( P τ ) cos ( k P τ + β 1 ) g 32 ( ρ , θ , Δ , τ ) = cos θ cos ( P τ ) cos ( k P τ + Δ + β 1 )
we get:
M ¯ i ( β 1 , ρ ) = μ 1 M ¯ i 1 ( β 1 , ρ ) + μ 2 M ¯ i 2 ( β 1 , ρ ) , i = 1 , 2 ,
for
M ¯ 1 j ( β 1 , ρ ) = 0 T I θ g 1 j + I Δ g 2 j d τ , M ¯ 2 j ( β 1 , ρ ) = 0 T g 3 j d τ , j = 1 , 2 .
To solve:
μ 1 M ¯ i 1 ( β 1 , ρ ) + μ 2 M ¯ i 2 ( β 1 , ρ ) = 0 , i = 1 , 2 ,
we first solve the scalar equation:
M ˜ ( β 1 , ρ ) = M ¯ 11 ( β 1 , ρ ) M ¯ 22 ( β 1 , ρ ) M ¯ 12 ( β 1 , ρ ) M ¯ 21 ( β 1 , ρ ) = 0
to get its root β 1 , 0 and ρ 0 . Then, we look for μ 1 , 0 and μ 2 , 0 with μ 1 , 0 2 + μ 2 , 0 2 = 1 such that:
μ 1 , 0 M ¯ i 1 ( β 1 , 0 , ρ 0 ) + μ 2 , 0 M ¯ i 2 ( β 1 , 0 , ρ 0 ) = 0 , i = 1 , 2 det μ 1 , 0 M ¯ 11 ( β 1 , 0 , ρ 0 ) + μ 2 , 0 M ¯ 12 ( β 1 , 0 , ρ 0 ) μ 1 , 0 M ¯ 21 ( β 1 , 0 , ρ 0 ) + μ 2 , 0 M ¯ 22 ( β 1 , 0 , ρ 0 ) 0 .
Summarizing, we have the following result.
Theorem 1.
If there are β 1 , 0 [ 0 , 2 π ) , ρ 0 satisfying (15), μ 1 , 0 and μ 2 , 0 with μ 1 , 0 2 + μ 2 , 0 2 = 1 solving (16), then for any μ 1 near μ 1 , 0 and μ 2 near μ 2 , 0 with μ 1 2 + μ 2 2 = 1 and ϵ 0 small, there are β 1 ( ϵ ) near β 1 , 0 and ρ ( ϵ ) near ρ 0 such that (8) with β 1 = β 1 ( ϵ ) and ρ = ρ ( ϵ ) has a T = 2 π / P -periodic solution near (12) with ρ = ρ ( ϵ ) .
Note:
M ¯ 11 ( β 1 , ρ ) = 0 T ( 2 cos ( 2 θ ( ρ , τ ) ) sin Δ ( ρ , τ ) cos θ ( ρ , τ ) cos ( P τ ) cos ( k P τ + β 1 ) + 2 cos θ ( ρ , τ ) cos Δ ( ρ , τ ) cos ( P τ ) sin ( k P τ + β 1 ) ) ) d τ M ¯ 12 ( β 1 , ρ ) = 0 T ( 2 cos ( 2 θ ( ρ , τ ) ) sin Δ ( ρ , τ ) sin θ ( ρ , τ ) cos ( P τ ) cos ( k P τ + Δ ( ρ , τ ) + β 1 ) + 2 sin θ ( ρ , τ ) cos Δ ( ρ , τ ) cos ( P τ ) sin ( k P τ + Δ ( ρ , τ ) + β 1 ) ) d τ M ¯ 21 ( β 1 , ρ ) = 0 T sin θ ( ρ , τ ) cos ( P τ ) cos ( k P τ + β 1 ) d τ M ¯ 22 ( β 1 , ρ ) = 0 T cos θ ( ρ , τ ) cos ( P τ ) cos ( k P τ + Δ ( ρ , τ ) + β 1 ) d τ .
Next, taking ρ in (12), we obtain:
θ ( , τ ) = P τ 2 τ 0 , π P , π P τ 2 τ π P , 2 π P , Δ ( , τ ) = π τ 0 , π P , 0 τ π P , 2 π P , π 2 τ 0 , π P , 2 π P , θ ρ ( , τ ) = 0 , Δ ρ ( , τ ) = 0 .
Hence:
M ¯ 11 ( β 1 , ) = 0 2 π P 2 cos P τ 2 cos ( P τ ) sin ( k P τ + β 1 ) d τ = 16 k ( 4 k 2 5 ) cos β 1 ( 16 k 4 40 k 2 + 9 ) P M ¯ 12 ( β 1 , ) = 0 2 π P 2 sin P τ 2 cos ( P τ ) sin ( k P τ + β 1 ) ) d τ = 8 ( 4 k 2 + 3 ) sin β 1 ( 16 k 4 40 k 2 + 9 ) P M ¯ 21 ( β 1 , ) = 0 2 π P sin P τ 2 cos ( P τ ) cos ( k P τ + β 1 ) d τ = 4 ( 4 k 2 + 3 ) cos β 1 ( 16 k 4 40 k 2 + 9 ) P M ¯ 22 ( β 1 , ) = 0 2 π P cos P τ 2 cos ( P τ ) cos ( k P τ + β 1 ) d τ = 8 k ( 4 k 2 5 ) sin β 1 ( 16 k 4 40 k 2 + 9 ) P .
Then, (15) gives as ρ ,
M ˜ ( β 1 , ) = 16 sin 2 β 1 ( 4 k 2 9 ) P 2
with asymptotic solutions β 1 , 0 satisfying either sin β 1 , 0 = 0 or cos β 1 , 0 = 0 . The asymptotic equation of (16) is as follows:
0 = μ 1 , 0 M ¯ 11 ( β 1 , 0 , ) + μ 2 , 0 M ¯ 12 ( β 1 , 0 , ) = μ 1 , 0 16 k ( 4 k 2 5 ) cos β 1 , 0 ( 16 k 4 40 k 2 + 9 ) P μ 2 , 0 8 ( 4 k 2 + 3 ) sin β 1 , 0 ( 16 k 4 40 k 2 + 9 ) P ,
which has solutions: μ 1 , 0 = 0 and μ 2 , 0 = ± 1 when sin β 1 , 0 = 0 or μ 1 , 0 = ± 1 and μ 2 , 0 = 0 when cos β 1 , 0 = 0 . However, M ¯ i j ρ ( β 1 , ) = 0 for i , j { 1 , 2 } , so Theorem 1 cannot be applied directly. Note we consider just the first asymptotic equation of (16), since the second one is a scalar multiple of the first one.
On the other hand, following the method of [5], p. 111, we get the following result.
Corollary 1.
For any ρ > 0 sufficiently large and ϵ 0 sufficiently small, there is μ 1 ( ρ , ϵ ) , μ 2 ( ρ , ϵ ) with μ 1 2 ( ρ , ϵ ) + μ 2 2 ( ρ , ϵ ) = 1 and β 1 ( ρ , ϵ ) such that (8) has a T-periodic solution near (17) for ϵ 0 small and either μ 1 ( , 0 ) = 0 , μ 2 ( , 0 ) = ± 1 , sin β 1 ( , 0 ) = 0 or μ 1 ( , 0 ) = ± 1 , μ 2 ( , 0 ) = 0 , cos β 1 ( , 0 ) = 0 .
Proof. 
The bifurcation equation has the form (see [5], p. 111):
μ 1 M ¯ i 1 ( β 1 , ρ ) + μ 2 M ¯ i 2 ( β 1 , ρ ) = O ( ϵ ) , i = 1 , 2 , μ 1 2 + μ 2 2 = 1
i.e.,
sin Γ M ¯ i 1 ( β 1 , ρ ) + cos Γ M ¯ i 2 ( β 1 , ρ ) = O ( ϵ ) , i = 1 , 2
for μ 1 = sin Γ and μ 2 = cos Γ . For ρ = and ϵ = 0 , (18) takes the form:
16 k ( 4 k 2 5 ) ( 16 k 4 40 k 2 + 9 ) P sin Γ cos β 1 8 ( 4 k 2 + 3 ) ( 16 k 4 40 k 2 + 9 ) P cos Γ sin β 1 = 0 4 ( 4 k 2 + 3 ) cos β 1 ( 16 k 4 40 k 2 + 9 ) P sin Γ cos β 1 8 k ( 4 k 2 5 ) ( 16 k 4 40 k 2 + 9 ) P cos Γ sin β 1 = 0 .
The determinant of (19) is 32 ( 4 k 2 9 ) P 2 0 , so (19) has the only solutions:
sin Γ = sin β 1 , 0 = 0 , cos Γ = cos β 1 , 0 = 0 .
The determinants of Jacobians of (19) at these zeros are as follows:
det 16 k ( 4 k 2 5 ) ( 16 k 4 40 k 2 + 9 ) P cos Γ cos β 1 , 0 8 ( 4 k 2 + 3 ) ( 16 k 4 40 k 2 + 9 ) P cos Γ cos β 1 , 0 4 ( 4 k 2 + 3 ) cos β 1 , 0 ( 16 k 4 40 k 2 + 9 ) P cos Γ cos β 1 , 0 8 k ( 4 k 2 5 ) ( 16 k 4 40 k 2 + 9 ) P cos Γ cos β 1 , 0 = 32 cos 2 Γ cos 2 β 1 , 0 ( 4 k 2 9 ) P 2 0
and:
det 8 ( 4 k 2 + 3 ) ( 16 k 4 40 k 2 + 9 ) P sin Γ sin β 1 , 0 16 k ( 4 k 2 5 ) ( 16 k 4 40 k 2 + 9 ) P sin Γ sin β 1 , 0 8 k ( 4 k 2 5 ) ( 16 k 4 40 k 2 + 9 ) P sin Γ sin β 1 , 0 4 ( 4 k 2 + 3 ) cos β 1 , 0 ( 16 k 4 40 k 2 + 9 ) P sin Γ sin β 1 , 0 = 32 sin 2 Γ sin 2 β 1 , 0 ( 4 k 2 9 ) P 2 0 ,
respectively. Hence, the zeroes Γ and β 1 , 0 are simple, so we can apply the implicit function theorem to get the result. The proof is complete. ☐

5. Discussion

Melnikov analysis is applied for the persistence of periodic oscillations for periodically-perturbed systems of ODEs with a slowly-varying variable. The ODEs are obtained from a model of a finite lattice of particles coupled by linear springs under distributed harmonic excitation, which presents a low-energy nonlinear acoustic vacuum (see [1]), but we consider just two modes. We extend the study of [1] to a problem with small exciting harmonic forces. We apply an analytical method based on derivation of Melnikov functions and then on the location of their simple roots. Melnikov functions are derived by using the approach of [3,5] developed for slowly-varying ODEs. Since the Melnikov functions are rather complicated, we follow an asymptotic way for solving the corresponding Melnikov equations. It would be nice to solve these Melnikov equations numerically for finding other simple roots, which is postponed to our next research. These roots determine and locate periodic solutions of the periodically-perturbed systems of ODEs derived from the two-mode low-energy nonlinear acoustic vacuum system. Moreover, our next investigation will be also to consider higher numbers of modes represented by a system of ODEs in (5). The method used will be different from that in this paper, since it will be based on the results of Section 3.3 of [5]. Note that higher modes of (2) were numerically studied in [2], which is another challenge for our study.

Author Contributions

The authors contributed equally to this work.

Funding

M.A. is supported by the National Scholarship Programme of the Slovak Republic for the Support of Mobility of Students, Ph.D. Students, University Teachers, Researchers and Artists. M.F. and M.P. are supported by the Slovak Research and Development Agency (Grant Number APVV-14-0378) and the Slovak Grant Agency VEGA(Grant Numbers 2/0153/16 and 1/0078/17).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Manevich, L.I.; Vakakis, A.F. Nonlinear oscillatory acoustic vacuum. SIAM J. Appl. Math. 2014, 74, 1742–1762. [Google Scholar] [CrossRef]
  2. Zhang, Z.; Manevitch, L.I.; Smirnov, V.; Bergman, L.; Vakakis, A.F. Extreme nonlinear energy exchanges in a geometrically nonlinear lattice oscillating in the plane. J. Mech. Phys. Solids 2018, 110, 1–20. [Google Scholar] [CrossRef]
  3. Wiggins, S.; Holmes, P. Periodic orbits in slowly varying oscillators. SIAM J. Math. Anal. 1987, 18, 592–611. [Google Scholar] [CrossRef]
  4. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 7th ed.; Academic Press: Cambridge, MA, USA, 2007. [Google Scholar]
  5. Fečkan, M. Topological Degree Approach to Bifurcation Problems; Springer: Berlin, Germany, 2008. [Google Scholar]
  6. Hale, J.K. Ordinary Differential Equations, 2nd ed.; Robert E. Krieger Publishing Company, Inc.: Malabar, FL, USA, 1980. [Google Scholar]
Figure 1. Forced and damped lattice oscillating in the plane (see [2]).
Figure 1. Forced and damped lattice oscillating in the plane (see [2]).
Axioms 07 00069 g001
Figure 2. Top panel: Graph for θ for different initial values of θ 0 . Bottom panel: Graph for Δ for different initial values of θ 0 .
Figure 2. Top panel: Graph for θ for different initial values of θ 0 . Bottom panel: Graph for Δ for different initial values of θ 0 .
Axioms 07 00069 g002
Figure 3. Orbits in the phase portrait of (10), where ( θ , Δ ) [ 0 , π 2 ] × [ 0 , π ] .
Figure 3. Orbits in the phase portrait of (10), where ( θ , Δ ) [ 0 , π 2 ] × [ 0 , π ] .
Axioms 07 00069 g003

Share and Cite

MDPI and ACS Style

Agaoglou, M.; Fečkan, M.; Pospíšil, M.; Rothos, V.M.; Vakakis, A.F. Periodically Forced Nonlinear Oscillatory Acoustic Vacuum. Axioms 2018, 7, 69. https://doi.org/10.3390/axioms7040069

AMA Style

Agaoglou M, Fečkan M, Pospíšil M, Rothos VM, Vakakis AF. Periodically Forced Nonlinear Oscillatory Acoustic Vacuum. Axioms. 2018; 7(4):69. https://doi.org/10.3390/axioms7040069

Chicago/Turabian Style

Agaoglou, Makrina, Michal Fečkan, Michal Pospíšil, Vassilis M. Rothos, and Alexander F. Vakakis. 2018. "Periodically Forced Nonlinear Oscillatory Acoustic Vacuum" Axioms 7, no. 4: 69. https://doi.org/10.3390/axioms7040069

APA Style

Agaoglou, M., Fečkan, M., Pospíšil, M., Rothos, V. M., & Vakakis, A. F. (2018). Periodically Forced Nonlinear Oscillatory Acoustic Vacuum. Axioms, 7(4), 69. https://doi.org/10.3390/axioms7040069

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