Next Article in Journal
Aczel–Alsina Power Aggregation Operators for Complex Picture Fuzzy (CPF) Sets with Application in CPF Multi-Attribute Decision Making
Next Article in Special Issue
The Solitary Solutions for the Stochastic Jimbo–Miwa Equation Perturbed by White Noise
Previous Article in Journal
Ptolemy’s Theorem in the Relativistic Model of Analytic Hyperbolic Geometry
Previous Article in Special Issue
Stochastic Solitons in Birefringent Fibers for Biswas–Arshed Equation with Multiplicative White Noise via Itô Calculus by Modified Extended Mapping Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Wave Propagation for a Strain Wave Equation of a Flexible Rod with Finite Deformation

by
Aqilah Aljuaidan
1,
Mamdouh Elbrolosy
1,2 and
Adel Elmandouh
1,3,*
1
Department of Mathematics and Statistics, College of Science, King Faisal University, P.O. Box 400, Al-Ahsa 31982, Saudi Arabia
2
Department of Mathematics, Faculty of Science, Tanta University, Tanta 31527, Egypt
3
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(3), 650; https://doi.org/10.3390/sym15030650
Submission received: 10 February 2023 / Revised: 27 February 2023 / Accepted: 3 March 2023 / Published: 5 March 2023

Abstract

:
The present work is attentive to studying the qualitative analysis for a nonlinear strain wave equation describing the finite deformation elastic rod taking into account transverse inertia, and shearing strain. The strain wave equation is rewritten as a dynamic system by applying a particular transformation. The bifurcation of the solutions is examined, and the phase portrait is depicted. Based on the bifurcation constraints, the integration of the first integral of the dynamic system along specified intervals leads to real wave solutions. We prove the strain wave equation has periodic, solitary wave solutions and does not possess kink (or anti-kink) solutions. In addition, the set of discovered solutions contains Jacobi-elliptic, trigonometric, and hyperbolic functions. This model contains many kinds of solutions, which are always symmetric or anti-symmetric in space. We study how the change in the physical parameters impacts the solutions that are found. Numerically, the behavior of the strain wave for the elastic rod is examined when particular periodic forces act on it, and moreover, we clarify the existence of quasi-periodic motion. To clarify these solutions, we present a 3D representation of them and the corresponding phase orbit.

1. Introduction

Many natural phenomena, including fluid dynamics, water waves, optical fibers, plasma, and nuclear physics are governed by nonlinear partial differential Equations (NLPDEs). It is widely believed that finding exact solutions to these phenomena is one of the most effective ways to comprehend and interpret them. Thus, the construction of solutions for NLPDEs has become a more significant and necessary tool for researchers. Despite this, there is no unified method that can provide all exact solutions of an NLPDE because these equations involve multifarious states and properties, which make it challenging to identify their exact analytical solutions. The most important of these methods in the literature include the extended modified direct algebraic method [1,2], the exponential function [3], the improved auxiliary equation technique [4], ( G / G 2 ) -expansion function methods [5], the bilinear formalism method [6], the direct method of the Hirota and the linear superposition principle [7], the sinh-Gordon expansion method [8], an improved mapping approach [9], and modified extended rational expansion method [10], a modified direct algebraic method [11], the first integral method [12], the Sardar-subequation method (SSM) [13,14,15,16], ϕ 6 model expansion technique [17], the Riccati-Bernoulli sub-ODE and e x p ( G / G ) -expansion method [18], an extended mapping technique [19], the extended exponential function method [20], a modified F-expansion method [21], Lie symmetry analysis [22,23], and for other several methods, see, e.g., [24,25,26,27,28,29,30]. Many of these methods are based on imposing a solution in a specific formula that contains a polynomial satisfying a certain ordinary differential equation. Another approach is a qualitative analysis of the traveling wave system related to the given partial differential equation utilizing the bifurcation theory of the dynamical system. For more details about this approach, see, e.g., [31,32,33,34,35,36,37,38,39]. Moreover, these methods are also applied to solve stochastic partial differential Equations [40,41,42]. The most important advantage of the bifurcation method lies in predicting the solution before calculating it. Recently, Elmandouh and Elborolosy [43,44,45,46,47,48,49] improved the procedures of this approach by presenting the interval of real propagation which prevents the appearance of complex solutions and permits various solutions for the same energy level, as well as presenting the degeneracy analysis of the solutions through the orbits of the phase portrait owing to the change in the initial conditions.
On the other side, numerous technological issues have been addressed using the nonlinear elastic wave. Nonlinearities of solid structures have various sources, e.g., physical and geometrical nonlinearities, kinetic nonlinearity, and boundary constraint. Solitary wave solutions or shock wave solutions are examples of steady traveling wave solutions that may exist as a result of the interaction between nonlinearity and dispersion or dissipation effect. There has been a great deal of interest in analyzing and solving nonlinear wave problems qualitatively using the nonlinear evolution equation since more and more problems involve nonlinearity. This motivates us to study some nonlinear dynamical behaviors of the wave strain equation in the form [50]
2 u t 2 c 0 2 2 u x 2 = 2 x 2 c 0 2 2 ( 3 u 2 + u 3 ) + ν 2 R 2 2 2 u t 2 c 1 2 2 u x 2 ,
where c 0 = E ρ is the wave velocity of the longitudinal and c 1 = μ ρ is the velocity of shear wave while ρ is the rod density per unite volume, μ is the material shear modulus, and E is the elastic modules, u ( x , t ) is the displacement distribution at each space x and time t. Equation (1) is a model for an elastic circular-rod waveguide that describes a double nonlinear wave equation concerning axial displacement gradient. As a consequence of the transverse Poisson effect, the longitudinal wave propagates simultaneously with the shear wave. Equation (1) is formulated by applying the Hamilton principle of the least action [50]. Recently, few methods have been proposed for constructing exact solutions to this problem; see, e.g., [51,52,53]. Only solitary wave solutions and shock wave solutions can be obtained, which are typically not able to yield an exact periodic solution. By utilizing the Jacobi elliptic function expansion method, one can obtain periodic, shock, and the corresponding solitary wave solutions of the derived nonlinear Equation [54].
This work is organized as follows: In Section 2, we study the bifurcation and introduce the phase portrait corresponding to the traveling wave system. Some wave solutions assorted into periodic, super-periodic, and solitary wave solutions are introduced in Section 3 as well as the proof of non-existing kink (anti-kink) wave solutions. Section 4 includes some graphic representations of periodic, super-periodic, and solitary solutions besides the study of the effect of changing the physical parameters on the obtained solutions. In Section 5, we examine numerically the existence of quasi-behaviour after allowing certain periodic forces to act on the rod. Section 6 includes a discussion about the main results while Section 7 is a summary of the obtained results.

2. Bifurcation Analysis and Phase Portraits

To investigate the dynamical analysis for wave solutions of Equation (1), we apply the next wave transformation to Equation (1)
u = u ( ξ ) , ξ = k ( x ω t ) ,
where k and ω are constants denoting the number of waves and speed of the waves, respectively, and ξ is the wave variable, turning Equation (1) into
( ω 2 c 0 2 ) d 2 u d ξ 2 c 0 2 2 d 2 d ξ 2 ( 3 u 2 + u 3 ) ν 2 R 2 k 2 2 d 2 d ξ 2 [ ( ω 2 c 1 2 ) d 2 u d ξ 2 ] = 0 .
Equation (3) is integrated twice with respect to ξ , we obtain
d 2 u d ξ 2 = c u + 3 a u 2 + a u 3 ,
where the integration constants are ignored and the two constants a , c given below are utilized instead of the main constants for simplicity
c = 2 ( ω 2 c 0 2 ) ν 2 R 2 k 2 ( ω 2 c 1 2 ) and a = c 0 2 ν 2 R 2 k 2 ( ω 2 c 1 2 ) .
Equation (4) is rewritten as a 1D-Hamilton system
u = z , z = a u 3 + 3 a u 2 + c u ,
where the derivative with regard to ξ is referred by primes, with Hamilton function
H = 1 2 z 2 + V ( u ) ,
where
V ( u ) = c 2 u 2 a u 3 a 4 u 4 ,
is a potential function. The Hamilton system (6) is conservative because div ( u , z ) = 0 , and based on H ξ = 0 , H is a constant of the motion, i.e., it takes a constant value along any trajectory of the motion [55,56], i.e.,
1 2 z 2 + V ( u ) = h ,
where h is the value of the constant of the motion along a specific trajectory and it is usually determined from the initial conditions. It is noticeable that both problems of establishing wave solutions and the solution of the Hamilton system (6) are equivalents. This equivalence leads to the acquisition of many properties of wave solutions, which will be discussed later. We acquire the differential form by inserting the first expression in Equation (6) and splitting the variables.
d u Q 4 ( u ) = ± d ξ ,
where
Q 4 ( u ) = 2 ( h V ( u ) ) = 2 h + c u 2 + 2 a u 3 + a 2 u 4 .
The needed range of the three parameters h , a , and c is required to integrate both sides of Equation (10). There are two different methods that can be applied to get this range. They are the bifurcation theory [31] and the complete discriminant system of the polynomial Q 4 ( u ) [57]. The bifurcation theory has various advantages over the other approach, such as the ability to assort the type of constructed solutions prior to forming them. Homoclinic orbits, periodic orbits, and heteroclinic orbits, for example, imply the appearance of solitary, periodic, and kink(or anti-kink) solutions. In other words, utilizing the bifurcation analysis, we determine the criteria on the parameters that ensure the occurrence of such solutions besides studying the dependence of the solutions on the initial conditions.
To study the bifurcation and phase portraits for the Hamilton system (6), we first find the equilibrium points which are the critical points for the potential function (8), i.e., the equilibria are ( u 0 , 0 ) , where ( u 0 , 0 ) satisfies
d V ( u 0 ) d u = u 0 ( a u 0 2 + 3 a u 0 + c ) = 0 .
The nature of the equilibrium points ( u 0 , 0 ) can be determined by evaluating the eigenvalues of the linearized system for the Hamilton system (6), that takes the form
λ 1 , 2 = ± d 2 V ( u 0 ) d u 2 = ± 3 a u 0 2 + 6 a u 0 + c .
Thus, the equilibrium point ( u 0 , 0 ) is either center if it is a local minimum for the potential function (8) or saddle if it is local maximum or cusp if V ( u 0 ) = 0 .
Depending on the value of Δ = a ( 9 a 4 c ) , we find the equilibrium points and specify their nature in the next cases:
Case I: If Δ < 0 , the Hamiltonian system (6) has only one equilibrium point O = ( 0 , 0 ) . The eigenvalues (13) calculated at O are λ 1 , 2 = ± c . It is clear that the two parameters a and c have the same sign, i.e., a c > 0 . Hence, O is either a saddle point if c > 0 ( a > 0 ) as outlined by Figure 1a or a center if c < 0 ( a < 0 ) as outlined by Figure 1b.
Case II: If Δ = 0 , then the two points O = ( 0 , 0 ) and E 1 = ( 3 2 , 0 ) are equilibria for system (6), where a 0 . The eigenvalues (13) evaluated at the points O and E 1 are
λ 1 , 2 ( O ) = ± 3 a 2 , λ 1 , 2 ( E 1 ) = 0 .
Taking into account the existence condition Δ = 0 for two equilibrium points, the two parameters a and c have the same sign, i.e., a c > 0 . Hence, they are either saddle and cusp if a > 0 ( c > 0 ) as outlined by Figure 2a or center and cusp, respectively, as clarified by Figure 2b if a < 0 ( c < 0 ).
Case III: If Δ > 0 , system (6) owns three equilibria. They are O = ( 0 , 0 ) and E 2 , 3 = ( 3 a ± Δ 2 a , 0 ) . The eigenvalues (13) calculated at these points are
λ 1 , 2 ( O ) = ± c , λ 1 , 2 ( E 2 ) = ± 2 a ( Δ 3 a Δ ) 2 a , λ 1 , 2 ( E 3 ) = ± 2 a ( Δ + 3 a Δ ) 2 a .
Based on the eigenvalues (15), we sum up the classification of the equilibria, O and E 2 , 3 , in Table 1 so as to avoid confusion.
Table 1. Classification of the equilibria O , E 2 , 3 if Δ > 0 .
Table 1. Classification of the equilibria O , E 2 , 3 if Δ > 0 .
CaseConditionsNatureFigure
a c O E 2 E 3
1.++saddlecentersaddleFigure 3a
2.centercentersaddleFigure 3b
3.+centersaddlesaddleFigure 4a
4.+saddlecentercenterFigure 4b
Figure 3. Phase portrait for the Hamiltonian system (6) in the phase plane ( u , z ) if Δ > 0 and a c > 0 . The black solid circles are the equilibrium points; (a) a = 1 ,   c = 5 / 4 ; (b) a = −1, c = −5/4.
Figure 3. Phase portrait for the Hamiltonian system (6) in the phase plane ( u , z ) if Δ > 0 and a c > 0 . The black solid circles are the equilibrium points; (a) a = 1 ,   c = 5 / 4 ; (b) a = −1, c = −5/4.
Symmetry 15 00650 g003
Figure 4. Phase plane orbits for the system (6) if Δ > 0 and a c < 0 . The equilibria are indicated by the black solid circles; (a) a = 1 / 15   , c = 7 / 28 ; (b) a = −1/15, c = 7 / 28 .
Figure 4. Phase plane orbits for the system (6) if Δ > 0 and a c < 0 . The equilibria are indicated by the black solid circles; (a) a = 1 / 15   , c = 7 / 28 ; (b) a = −1/15, c = 7 / 28 .
Symmetry 15 00650 g004
The values of the parameter h at the equilibrium points are given by
h 0 = V ( 0 ) = 0 , h 1 = V ( 3 2 ) = 27 a 64 , h 2 = V ( 3 a + Δ 2 a ) = ( 3 a Δ ) 2 32 a 2 [ 2 c 3 a + Δ ] , h 3 = V ( 3 a Δ 2 a ) = ( 3 a + Δ ) 2 32 a 2 [ 2 c 3 a Δ ] .

3. Solutions

The aim of this section is to obtain some bounded wave solutions for Equation (1). These solutions are assorted into periodic, solitary, and kink (anti-kink) wave solutions. Besides, the same techniques can also be employed to construct unbounded wave solutions, but we avoid them because they are not physically meaningful. We now present the subsequent lemma, which will be substantial in the next analysis.
Lemma 1.
(see, e.g., [56,58,59]). Let system (6) has a continuous solution u = u ( k ( x c t ) ) = u ( ξ ) for ξ ] , [ and assume lim ξ u ( ξ ) = κ 1 ,   a n d lim ξ u ( ξ ) = κ 2 .
(i)  
If κ 1 = κ 2 , then the solution u ( ξ ) is solitary which corresponds to a homoclinic orbit for the system (6).
(ii)  
If κ 1 κ 2 , then the solution u ( ξ ) is a kink (or anti-kink) wave that corresponds to a heteroclinic orbit for the system (6).
(iii)  
If system (6) possesses a periodic orbit, then its corresponding solution u ( ξ ) is also periodic.
(iv)  
If system (6) has a closed orbit in the phase portrait evolved by at least two centers and one separatrix layer, then its corresponding solution u ( ξ ) is a super periodic wave.
Based on the bifurcation theory and Lemma 1, we prove the next theorem.
Theorem 1.
Equation (1) does not possess any kink or anti-kink solutions.
Proof. 
The bifurcation analysis, outlined in Section 2, shows that system (6) has not any heteroclinic orbit. Taking into account Lemma 1, we then proved that Equation (1) has not any kink or anti-kink solutions. □

3.1. Periodic Solutions

We are interested in this subsection in deriving all periodic solutions of Equation (1). As we can see from the bifurcation analysis, Equation (1) has several periodic and super periodic solutions of different shapes. The following theorems enumerate these cases
Theorem 2.
Assume Equation (1) has a solution in the form of (2).
(i)  
If ( Δ , a , c , h ) ( R × R × R × R + ) ( { 0 } × R × R × ] 0 , h 1 [ ) ( R + × R × R × ] h 2 , 0 [ ) ( R + × R × R + × ] h 2 , h 3 [ ) , then Equation (1) has new periodic solutions in the form
u ( ξ ) = 1 A B A u 1 B u 2 2 A B ( u 1 u 2 ) A + B + ( A B ) cn ( a A B 2 ξ , k 1 ) ,
for u 1 < u < u 2 , where A 2 = ( u 2 Re . u 3 ) 2 + Im 2 u 3 , B 2 = ( u 1 Re . u 3 ) 2 + Im 2 u 3 , k 1 2 = 1 4 A B [ ( u 1 u 2 ) 2 ( A B ) 2 ] and u 1 , u 2 are the two real roots while u 3 , u 3 * are the two complex conjugate roots of Q 4 ( u ) .
(ii)  
If ( Δ , a , c , h ) ( { 0 } × R × R × ] h 1 , [ ) ( R + × R × R × ] h 3 , [ ) ) ( R + × R × R + × ] 0 , [ ) , then Equation (1) has super periodic solutions as in the form (17) but with different arguments.
Proof. 
(i) 
For fixed values of Δ , a , c , h in the given range, system (6) owns different families of phase orbits. An orbit belonging to these families passes through the u -axis twice which proves the existence of two real roots u 1 , u 2 and two complex conjugate complex roots u 3 , u 3 * for Q 4 ( u ) . This enables us to write Q 4 ( u ) = a 2 ( u 1 u ) ( u u 2 ) ( u u 3 ) ( u u 3 * ) and consequently, the interval of real propagation is u ] u 1 , u 2 [ . Therefore, we assume u ( 0 ) = u 1 . Integration of Equation (10) gives
u 1 u d u ( u 1 u ) ( u u 2 ) ( u u 3 ) ( u u 3 * ) = a 2 0 ξ d ξ .
The last equation gives the new periodic solutions (17) with a period of 4 2 a A B K ( k 1 ) [60].
(ii) 
For certain values of Δ , a , c , h in the specific range, system (6) different families of super periodic orbits. An orbit of them includes all equilibrium points of the system and intersects the u -axis in two points as in case (i). consequently, we get solutions as in (17) but with different arguments.
Theorem 3.
Assume Equation (1) has a solution in the form of (2). If ( Δ , a , c , h ) ( R + × R × R × { 0 } ) ( R + × R × R + × { h 3 } ) , then Equation (1) has periodic solutions in the form
u ( ξ ) = u e 2 ( u 4 u e ) ( u 5 u e ) 2 u e u 4 u 5 ( u 5 u 4 ) cos a 8 ( u 4 u e ) ( u 5 u e ) ( u 5 u 4 ) ξ ,
for u 4 < u < u 5 where
u e = 0 , if h = 0 3 a Δ 2 a , if h = h 3
and u 4 < u 5 < u e are the real roots of Q 4 ( u ) .
Proof. 
For selecting values of Δ , a , c , h in the given range, system (6) possesses a single orbit in red as illustrated by Figure 3b and it intersects u -axis in u 4 , u 5 and u e . Thus, Q 4 ( u ) reads as Q 4 ( u ) = a 2 ( u u e ) 2 ( u 4 u ) ( u u 5 ) . The possible interval of real wave propagation ] u 4 , u 5 [ . Assuming u ( 0 ) = u 4 , the integration of (10) yields
u 4 u d u u ( u u 4 ) ( u 5 u ) = a 2 0 ξ d ξ ,
which gives the solution as in the form (19). □
Theorem 4.
Assume Equation (1) has a solution in the form of (2). If ( Δ , a , c , h ) ( R + × R × R × ] 0 , h 3 [ ) ( R + × R × R + × ] h 3 , 0 [ ) , then Equation (1) has periodic solutions in the form
u ( ξ ) = u 9 + ( u 9 u 7 ) ( u 9 u 6 ) ( u 7 u 9 ) ( u 7 u 6 ) sn 2 ( Ω 1 ξ , k 2 ) ,
for u 6 < u < u 7 , and
u ( ξ ) = u 6 ( u 9 u 6 ) ( u 8 u 6 ) u 6 u 8 + ( u 8 u 9 ) sn 2 ( Ω 1 ξ , k 2 ) ,
for u 8 < u < u 9 , where Ω 1 = a 8 ( u 9 u 7 ) ( u 8 u 6 ) , k 2 2 = ( u 9 u 8 ) ( u 7 u 6 ) ( u 9 u 7 ) ( u 8 u 6 ) and u 6 < u 7 < u 8 < u 9 are the four real roots of Q 4 ( u ) .
Proof. 
For selecting values of Δ , a , c , h in the given range, system (6) has two periodic families of orbits as outlined in Figure 3b and Figure 4b in green. An orbit belonging these families passes four times through u -axis which indicates the existence of four real zeros for Q 4 ( u ) , i.e., Q 4 ( u ) = a 2 ( u 6 u ) ( u u 7 ) ( u u 8 ) ( u u 9 ) . There are two possible intervals of real wave propagation. They are u ] u 6 , u 7 [ ] u 8 , u 9 [ . We calculate the periodic wave solution along each interval individually.
  • If u ] u 6 , u 7 [ , we assume u ( 0 ) = u 6 , the integration of Equation (10) yields
    u 6 u d u ( u 6 u ) ( u u 7 ) ( u u 8 ) ( u u 9 ) = a 2 0 ξ d ξ .
    The last equation gives the new periodic solution (22) with a period of 2 Ω 1 K ( k 2 ) .
  • If u ] u 8 , u 9 [ , we postulate u ( 0 ) = u 9 . The integration of Equation (10) implies
    u 9 u d u ( u 6 u ) ( u u 7 ) ( u u 8 ) ( u u 9 ) = a 2 0 ξ d ξ .
    The last equation gives the new periodic solution (23) with a period of 2 Ω 1 K ( k 2 ) .
Theorem 5.
Assume Equation (1) has a solution in the form of (2). If ( Δ , a , c , h ) ( R + × R + × R + × ] h 2 , 0 [ ) ( R + × R + × R × ] 0 , h 2 [ ) , then Equation (1) has periodic solutions in the form
u ( ξ ) = u 11 + ( u 10 u 12 ) ( u 10 u 11 ) u 12 u 10 + ( u 11 u 12 ) sn 2 ( Ω 2 ξ , k 3 ) ,
for u 11 < u < u 12 where Ω 2 = a 8 ( u 13 u 11 ) ( u 12 u 10 ) , k 3 2 = ( u 12 u 11 ) ( u 13 u 10 ) ( u 13 u 11 ) ( u 12 u 10 ) and u 10 < u 11 < u 12 < u 13 are the real roots of Q 4 ( u ) .
Proof. 
For selecting values of Δ , a , c , h in the given range, system (6) has one periodic families of orbits as outlined in Figure 3a and Figure 4a in green. An orbit belonging these family intersects u -axis in u 10 , u 11 , u 12 , u 13 which indicates the existence of four real zeros for Q 4 ( u ) , i.e., Q 4 ( u ) = a 2 ( u 10 u ) ( u 11 u ) ( u 12 u ) ( u 13 u ) . We calculate the periodic wave solution along u ] u 11 , u 12 [ with u ( 0 ) = u 11 to give
u 11 u d u ( u 10 u ) ( u 11 u ) ( u 12 u ) ( u 13 u ) = a 2 0 ξ d ξ ,
which gives the new periodic solution as in the form (26) with a period of 2 Ω 2 K ( k 3 ) .

3.2. Solitary Solution

This subsection aims to construct solitary wave solutions for Equation (1) which correspond to the homoclinic orbits for the traveling wave system (6).
Theorem 6.
Assume Equation (1) has a solution in the form of (2). If ( Δ , a , c , h ) ( R + × R × R × { h 3 } ) ( R + × R × R + × { 0 } ) , then Equation (1) has solitary wave solutions in the form
u ( ξ ) = u e 2 ( s 2 u e ) ( s 1 u e ) ( s 1 s 2 ) cos h Ω 3 ξ + 2 u e s 1 s 2 .
for s 1 < u < u e , and
u ( ξ ) = u e + 2 ( s 2 u e ) ( s 1 u e ) ( s 1 s 2 ) cos h Ω 3 ξ + s 1 + s 2 2 u e ,
for u e < u < s 2 , where Ω 3 = a 2 ( u e s 1 ) ( s 2 u e ) , and s 1 < u e < s 2 are the real roots of Q 4 ( u ) .
Proof. 
For selecting values of Δ , a , c , h in the given range, system (6) has homoclinic orbits in black as outlined as Figure 3b and Figure 4b, and an orbit of them cuts u -axis in three points. Therefore, the polynomial Q 4 ( u ) has two simple roots s 1 , s 2 and the other u e is double. Q 4 ( u ) reads as Q 4 ( u ) = a 2 ( u e u ) 2 ( u s 1 ) ( s 2 u ) and the intervals of real wave propagation is u ] s 1 , u e [ ] u e , s 2 [ . Integrate both sides of Equation (10) for u ] s 1 , u e [ and u ] u e , s 2 [ assuming u ( 0 ) = s 1 and u ( 0 ) = s 2 , respectively, we get the new solitary wave solutions (28) and (29). □

4. Graphical Representation

There are two goals for this section. The first is to graphically illustrate some of the solutions we have found, and the second is to investigate how altering one physical parameter while keeping the others constants, will affect the solution’s attitude. It is noted from Equation (3) that the variation in the parameters ν , R , and k have the same effect on the solutions, so we will only study one of them, say ν .
Let us assume a = 1 , c = 3 and select a suitable initial condition implying to h = 0.2 . System (6) has a periodic orbit around the center point ( 0 , 0 ) as outlined by Figure 5a. Thus, the polynomial (11) has two real roots u 1 = 0.4221201674 ,   u 2 = 0.3283297824 and two complex conjugate roots u 3 , u 3 * = 1.953104807 ± 1.399146643 i . Thanks to theorem 2, Equation (1) has a periodic solution in the form
u ( x , t ) = 13.83235869 4.750306908 + 0.602283150 c n ( 1.665933324 x + 0.3331866648 t , 0.09501209619 ) 3.006356686 .
The 3D graphic representation of the solution (30) shows it is periodic as it is illustrated by Figure 5b.
The influence of varying the parameters on the periodic solutions is clarified in Figure 6, where the wave contracts as both its wavelength and its amplitude decrease with the increase of c 0 and c 1 as in Figure 6a,b, assuming the other parameters are fixed, while it expands with the increase of ν and w as in Figure 6c,d.
For the values a = 1 , c = 1.25 , we select a suitable initial condition implying h = 0.08 that satisfies h > h 2 = 1.953125 . Consequently, system (6) has a super periodic orbit as outlined by Figure 7a. The polynomial Q ( u ) in (11) has two real zeros which are also the intersection points of the orbit with u -axis (see Figure 7a), i.e., u 1 = 3.237 , u 2 = 0.291637 and the two conjugate complex roots u 3 , u 3 * = 0.5274221118 ± 0.2457027386 i . Following Theorem 2, Equation (1) has a super-periodic solution in the form
u ( x , t ) = 8.807117860 3.576302665 1.865383825 cn ( 1.078788901 x + 0.2157577802 t , 0.9818388414 ) + 1.910446655 ,
that is sketched in Figure 7b.
Remark 1.
It is worth mentioning that the previous two cases clarify the difference between the two parts in Theorem 2. This is due to the two solutions (30) and (31) are obtained from the same Equation (17).
Figure 8 illustrates the impacts of changing the parameters on the super periodic solutions. It is noted that the wave shrinks by increasing the values of each c 0 and c 1 as in Figure 8a,b, while it grows by increasing values of ν and w as in Figure 8c,d.
Assuming a = 1 ,   c = 1.25 and selecting a suitable initial condition providing h = 0.046875 , system (6) has two homoclinic orbits each of them connecting the saddle point ( 0 , 0 ) with it self as outlined by Figure 9a and it intersects u -axis in the points ( s 1 , 0 ) = ( 3.232050808 , 0 ) and ( s 2 , 0 ) = ( 0.2320508076 , 0 ) . Following [56], by virtue of Theorem 6, there are two types of solitary wave solutions in the form
u ( x , t ) = 0.5 4 3.464101616 c o s h ( x + 0.2 t ) + 2 , if 3.232050808 < u 0.5 , 0.5 + 4 3.464101616 c o s h ( x + 0.2 t ) + 2 , if 0.5 u < 0.23205080765 .
These solutions are named rarefactive wave solitary and compressive wave solitary, for more details about these solutions see, e.g., [56]. These solutions are clarified graphically by Figure 9. We examine graphically the influence of the included parameters on the the compressive soliton solution by introducing Figure 10. We find that both wavelength and the height decrease with the increase of c 0 as in Figure 10a, the wavelength decreases while the height remains constant with the increase of c 1 as in Figure 10b, the wavelength increases while the height remains constant with the increase of ν as in Figure 10c, and both wavelength and the height enlarge with the increase of w as in Figure 10d.

5. Quasi Periodic Behaviour

In this section, we are concerned with a discussion of the autoresonance behavior for the non-autonomous system, i.e., the oscillator self-adjusts by subjecting the system to a variable periodic force. The perturbed form of problem (1) after inserting a periodic external force d cos ( e ( x c t ) ) will be
2 u t 2 c 0 2 2 u x 2 = 2 x 2 3 E 2 ρ u 2 + E 2 ρ u 3 + ν 2 R 2 2 2 u t 2 c 1 2 2 u x 2 + d cos ( e ( x c t ) ) .
Performing the transformation (2) to Equation (33), the corresponding perturbed dynamical system takes the form
u = z , z = a u 3 + 3 a u 2 + c u + α cos ( θ ) , θ = β ,
where θ = β ξ while α = 2 d v 2 R 2 k 2 ( c c 1 2 ) and β = e k refer, respectively, to the strength and frequency of the external periodic force. We investigate the behavior of the perturbed system (34) for distinct values of α by fixing both the remaining parameters and initial condition.
For fixed values a = c = 2 ,   β = 0.09 , we select a suitable initial condition u ( 0 ) = 0.5 , z ( 0 ) = 1 which implies the existence of a super periodic phase orbit outlined in Figure 3b in blue, for unperturbed system (6). The value of h corresponding to the initial condition is h = 1.0325 satisfying the condition h > h 3 = 1.0325 . In the absence of the external periodic force, i.e., α = 0 , the 2D and 3D phase portraits for the unperturbed system (6) are displayed in Figure 11a,d. To examine the behaviour for growing values of α , we introduce the 2D and 3D phase portraits for perturbed system (34) for α = 2 ,   α = 5 as outlined in Figure 11b,e and Figure 11c,f, respectively. We see that an increase in the strength of the external force reveals an increase in the periodic irregularity of the wave. It is noticeable the existence of a quasi-periodic behavior is due to the existence of two incommensurable frequencies in system (34).

6. Discussion

This section aims to discuss the obtained results. We investigate the dynamical behavior for Equation (1) that is modeled for the elastic circular-rod wave-guide. It describes a double nonlinear wave equation with respect to the axial displacement gradient. As a consequence of the transverse Poisson effect, the longitudinal wave propagates simultaneously with the shear wave. Certain transformation involving wave variable (2) is applied to Equation (1) in order to convert it to second order differential equation that is followed by writing it as a dynamical system which is equivalent to the Hamilton system describing the one dimension motion of a particle. This equivalent has two significant roles in the current study. We summarize them in the next items:
(a)
It enables us to determine the interval of real wave propagation which corresponds to the interval of real motions in Hamilton systems.
(b)
The process of finding the wave solutions for Equation (1) is reduced to quadrature as we see above. This is due to the Hamilton system (1) having one degree of freedom and consequently, it is integrable. All integrable systems are solved by quadrature.
The integration of both sides of Equation (10) requires the range of the parameters a , c , h . This range can be obtained by two distinct methods. They are bifurcation analysis and complete discriminate method of a polynomial Q 4 ( u ) . We apply the bifurcation analysis in this work for two reasons. They are
(a)
It gives us the required range of the parameters a , c , and h. It also enables us to determine the type of the solution before constructing them via the type of the phase plane orbits as it is clarified in Lemma 1. By virtue of these facts, we prove the non-existence of kink or (anti-kink) wave solutions for Equation (1) as a result of system (6) has no heteroclinic phase orbit.
(b)
It manages us to clarify the dependence of the solutions on the initial conditions. The constant h in Equation (9) is determined from the initial conditions. Thus, for distinct values of h, or equivalently, for different initial conditions, there are different wave solutions. Let us illustrate this point by providing an example. If ( Δ , a , c ) R + × R × R , then Equation (1) has either supper-periodic solution if h ] h 3 , [ , see theorem 2 or solitary solution if h = h 3 , see Theorem 6. Hence, for the same conditions on the physical parameters, the type of the solution depends on h which is always calculated from the initial conditions. Or, equivalently, the type of the solution depends on the initial conditions. Thus, we can also conclude that the bifurcation analysis enables us to find all possible wave solutions.
We illustrate graphically some of the obtained solutions by displaying the 3D-graphical representation and the corresponding phase orbit. Moreover, we investigate the influence of the included parameters c 0 , c 1 , ν , and ω on the solutions when one of them varies while the others are kept fixed. Figure 6, Figure 8 and Figure 10, clarify the influence of the parameters c 0 , c 1 , ν , ω on periodic, supper-periodic and solitary solutions, respectively. These effects appear on both the amplitude and the width of the solutions.
Finally, we study the existence of quasi-behavior for the Hamilton system (6) after allowing a periodic perturbed term which is an external periodic force to act on the rod. After adding this term, system (6) which is named here as an unperturbed system is converted into perturbed system (34). For selected values of the parameters a , c , α , β and suitable choice of the initial condition, the 2D and 3D phase portrait is outlined by Figure 11.

7. Conclusions

The study of some qualitative properties for the problem of a finite deformation flexible rod is really the focus of the present work. The governing equations, which are nonlinear partial differential equations, have been derived by applying Hamilton’s variational principle. The governing equation has been converted into one degree of freedom Hamiltonian system describing the one-dimensional motion of a particle utilizing specific wave transformation. In this context, the problem of finding wave solutions and the problem of establishing the solution of the Hamiltonian are equivalent. This equivalence is significant because it enables us to find the real wave solutions which have been limited by determining the interval of real wave propagation, which is also a permitted interval of particle real motions. In other words, it enables us to find bounded real-wave solutions that are desirable in real-world applications instead of complex solutions. Analyses of bifurcation and phase plane description have been carried out. This study has been utilized to prove the governing equation has no kink or anti-kink solutions and has solitary, periodic, and unbounded wave solutions. Based on the bifurcation constraints on the problem’s parameters, the conserved quantity has been integrated along the bounded phase plane orbits. Thus, we have discovered various periodic and solitary wave solutions that are displayed of them by exhibiting a 3D, contour, and the appropriate phase orbit. We have investigated how altering one physical parameter while preserving the other’s constant will affect the solutions we have found. Finally, We have examined numerically the phase plane after allowing periodic forces to act on the rod. We have concluded the perturbed system has a quasi-periodic behavior.
The investigation of chaotic behavior for system (34) and its applications to image encryption will be considered in our upcoming work.

Author Contributions

The authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

The authors extend their appreciation to the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia for funding this research work through project number Grant No. 2257.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported through the Annual Funding track by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [Grant.2257].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arshad, M.; Seadawy, A.R.; Lu, D. Elliptic function and solitary wave solutions of the higher-order nonlinear Schrödinger dynamical equation with fourth-order dispersion and cubic-quintic nonlinearity and its stability. Eur. Phys. J. Plus 2017, 132, 371. [Google Scholar] [CrossRef]
  2. Arshad, M.; Seadawy, A.; Lu, D.; Wang, J. Travelling wave solutions of Drinfel’d–Sokolov–Wilson, Whitham–Broer–Kaup and (2 + 1)-dimensional Broer–Kaup–Kupershmit equations and their applications. Chin. J. Phys. 2017, 55, 780–797. [Google Scholar] [CrossRef]
  3. Seadawy, A.R.; Nuruddeen, R.; Aboodh, K.; Zakariya, Y. On the exponential solutions to three extracts from extended fifth-order KdV equation. J. King Saud Univ.-Sci. 2020, 32, 765–769. [Google Scholar] [CrossRef]
  4. Nasreen, N.; Seadawy, A.R.; Lu, D.; Arshad, M. Solitons and elliptic function solutions of higher-order dispersive and perturbed nonlinear Schrödinger equations with the power-law nonlinearities in non-Kerr medium. Eur. Phys. J. Plus 2019, 134, 485. [Google Scholar] [CrossRef]
  5. Younas, U.; Seadawy, A.R.; Younis, M.; Rizvi, S. Optical solitons and closed form solutions to the (3 + 1)-dimensional resonant Schrödinger dynamical wave equation. Int. J. Mod. Phys. B 2020, 34, 2050291. [Google Scholar] [CrossRef]
  6. Hossen, M.B.; Roshid, H.O.; Ali, M.Z. Multi-soliton, breathers, lumps and interaction solution to the (2 + 1)-dimensional asymmetric Nizhnik-Novikov-Veselov equation. Heliyon 2019, 5, e02548. [Google Scholar] [CrossRef] [Green Version]
  7. Manukure, S.; Chowdhury, A.; Zhou, Y. Complexiton solutions to the asymmetric Nizhnik–Novikov–Veselov equation. Int. J. Mod. Phys. B 2019, 33, 1950098. [Google Scholar] [CrossRef]
  8. Cattani, C.; Sulaiman, T.A.; Baskonus, H.M.; Bulut, H. On the soliton solutions to the Nizhnik-Novikov-Veselov and the Drinfel’d-Sokolov systems. Opt. Quantum Electron. 2018, 50, 138. [Google Scholar] [CrossRef]
  9. Li, Z. Diverse oscillating soliton structures for the (2 + 1)-dimensional Nizhnik–Novikov–Veselov equation. Eur. Phys. J. Plus 2020, 135, 8. [Google Scholar] [CrossRef]
  10. Iqbal, M.; Seadawy, A.R.; Khalil, O.H.; Lu, D. Propagation of long internal waves in density stratified ocean for the (2 + 1)-dimensional nonlinear Nizhnik-Novikov-Vesselov dynamical equation. Results Phys. 2020, 16, 102838. [Google Scholar] [CrossRef]
  11. Seadawy, A.R.; Arshad, M.; Lu, D. The weakly nonlinear wave propagation of the generalized third-order nonlinear Schrödinger equation and its applications. Waves Random Complex Media 2022, 32, 819–831. [Google Scholar] [CrossRef]
  12. Taghizadeh, N.; Mirzazadeh, M.; Farahrooz, F. Exact solutions of the nonlinear Schrödinger equation by the first integral method. J. Math. Anal. Appl. 2011, 374, 549–553. [Google Scholar] [CrossRef] [Green Version]
  13. ur Rehman, H.; Awan, A.U.; Habib, A.; Gamaoun, F.; El Din, E.M.T.; Galal, A.M. Solitary wave solutions for a strain wave equation in a microstructured solid. Results Phys. 2022, 39, 105755. [Google Scholar] [CrossRef]
  14. Rehman, H.U.; Iqbal, I.; Subhi Aiadi, S.; Mlaiki, N.; Saleem, M.S. Soliton solutions of Klein–Fock–Gordon equation using Sardar subequation method. Mathematics 2022, 10, 3377. [Google Scholar] [CrossRef]
  15. Asjad, M.I.; Inc, M.; Iqbal, I. Exact solutions for new coupled Konno–Oono equation via Sardar subequation method. Opt. Quantum Electron. 2022, 54, 798. [Google Scholar]
  16. Asjad, M.I.; Munawar, N.; Muhammad, T.; Hamoud, A.A.; Emadifar, H.; Hamasalh, F.K.; Azizi, H.; Khademi, M. Traveling wave solutions to the Boussinesq equation via Sardar sub-equation technique. AIMS Math. 2022, 7, 11134–11149. [Google Scholar]
  17. Rehman, H.U.; Awan, A.U.; Allahyani, S.A.; Tag-ElDin, E.M.; Binyamin, M.A.; Yasin, S. Exact solution of paraxial wave dynamical model with Kerr Media by using ϕ6 model expansion technique. Results Phys. 2022, 42, 105975. [Google Scholar] [CrossRef]
  18. El-Ganaini, S. Solitons and other solutions to a new coupled nonlinear Schrodinger type equation. J. Egypt. Math. Soc. 2017, 25, 19–27. [Google Scholar] [CrossRef]
  19. Seadawy, A.R.; Arshad, M.; Lu, D. Dispersive optical solitary wave solutions of strain wave equation in micro-structured solids and its applications. Phys. A Stat. Mech. Its Appl. 2020, 540, 123122. [Google Scholar] [CrossRef]
  20. Kumar, S.; Kumar, A.; Wazwaz, A.M. New exact solitary wave solutions of the strain wave equation in microstructured solids via the generalized exponential rational function method. Eur. Phys. J. Plus 2020, 135, 870. [Google Scholar] [CrossRef]
  21. Seadawy, A.R.; Ali, A.; Baleanu, D.; Althobaiti, S.; Alkafafy, M. Dispersive analytical wave solutions of the strain waves equation in microstructured solids and Lax’fifth-order dynamical systems. Phys. Scr. 2021, 96, 105203. [Google Scholar] [CrossRef]
  22. Jadaun, V.; Kumar, S. Lie symmetry analysis and invariant solutions of (3 + 1)(3 + 1)-dimensional Calogero–Bogoyavlenskii–Schiff equation. Nonlinear Dyn. 2018, 93, 349–360. [Google Scholar] [CrossRef] [Green Version]
  23. Tu, J.M.; Tian, S.F.; Xu, M.J.; Zhang, T.T. On Lie symmetries, optimal systems and explicit solutions to the Kudryashov–Sinelshchikov equation. Appl. Math. Comput. 2016, 275, 345–352. [Google Scholar] [CrossRef]
  24. El-Dessoky, M.; Elmandouh, A. Qualitative analysis and wave propagation for Konopelchenko-Dubrovsky equation. Alex. Eng. J. 2023, 67, 525–535. [Google Scholar] [CrossRef]
  25. Adeyemo, O.D.; Khalique, C.M.; Gasimov, Y.S.; Villecco, F. Variational and non-variational approaches with Lie algebra of a generalized (3 + 1)-dimensional nonlinear potential Yu-Toda-Sasa-Fukuyama equation in Engineering and Physics. Alex. Eng. J. 2023, 63, 17–43. [Google Scholar] [CrossRef]
  26. Wang, K. Exact travelling wave solution for the local fractional Camassa-Holm-Kadomtsev-Petviashvili equation. Alex. Eng. J. 2023, 63, 371–376. [Google Scholar] [CrossRef]
  27. Ahmed, M.S.; Zaghrout, A.A.; Ahmed, H.M. Travelling wave solutions for the doubly dispersive equation using improved modified extended tanh-function method. Alex. Eng. J. 2022, 61, 7987–7994. [Google Scholar] [CrossRef]
  28. Ali, K.K.; Osman, M.; Abdel-Aty, M. New optical solitary wave solutions of Fokas-Lenells equation in optical fiber via Sine-Gordon expansion method. Alex. Eng. J. 2020, 59, 1191–1196. [Google Scholar] [CrossRef]
  29. Rehman, H.U.; Awan, A.U.; Tag-ElDin, E.M.; Alhazmi, S.E.; Yassen, M.F.; Haider, R. Extended hyperbolic function method for the (2 + 1)-dimensional nonlinear soliton equation. Results Phys. 2022, 40, 105802. [Google Scholar] [CrossRef]
  30. Rehman, H.; Seadawy, A.R.; Younis, M.; Rizvi, S.; Anwar, I.; Baber, M.; Althobaiti, A. Weakly nonlinear electron-acoustic waves in the fluid ions propagated via a (3 + 1)-dimensional generalized Korteweg–de-Vries–Zakharov–Kuznetsov equation in plasma physics. Results Phys. 2022, 33, 105069. [Google Scholar] [CrossRef]
  31. Nemytskii, V.V. Qualitative Theory of Differential Equations; Princeton University Press: Princeton, NJ, USA, 2015; Volume 2083. [Google Scholar]
  32. Elbrolosy, M.; Elmandouh, A. Bifurcation and new traveling wave solutions for (2 + 1)-dimensional nonlinear Nizhnik–Novikov–Veselov dynamical equation. Eur. Phys. J. Plus 2020, 135, 533. [Google Scholar] [CrossRef]
  33. Elmandouh, A. Bifurcation and new traveling wave solutions for the 2D Ginzburg–Landau equation. Eur. Phys. J. Plus 2020, 135, 648. [Google Scholar] [CrossRef]
  34. Elmandouha, A.; Ibrahim, A. Bifurcation and travelling wave solutions for a (2 + 1)-dimensional KdV equation. J. Taibah Univ. Sci. 2020, 14, 139–147. [Google Scholar] [CrossRef] [Green Version]
  35. Al Nuwairan, M. The exact solutions of the conformable time fractional version of the generalized Pochhammer–Chree equation. Math. Sci. 2022, 1–12. [Google Scholar] [CrossRef]
  36. Karakoc, S.B.G.; Saha, A.; Sucu, D.Y. A collocation algorithm based on septic B-splines and bifurcation of traveling waves for Sawada–Kotera equation. Math. Comput. Simul. 2023, 203, 12–27. [Google Scholar] [CrossRef]
  37. Pradhan, B.; Saha, A. Coexistence of Chaotic, Quasiperiodic and Multiperiodic Features in Quantum Plasma. In Proceedings of the Seventh International Conference on Mathematics and Computing—ICMC, Online, 2–5 March 2021; pp. 903–914. [Google Scholar]
  38. Saha, A.; Pradhan, B.; Natiq, H. Multiperiodic and chaotic wave phenomena of collective ion dynamics under KP-type equation in a magnetised nonextensive plasma. Phys. Scr. 2022, 97, 095604. [Google Scholar] [CrossRef]
  39. Prasad, P.K.; Mandal, U.K.; Das, A.; Saha, A. Phase plane analysis and integrability via Bäcklund transformation of nucleus-acoustic waves in white dwarf. Chin. J. Phys. 2021, 73, 534–545. [Google Scholar] [CrossRef]
  40. Mohammed, W.W.; Al-Askar, F.M.; Cesarano, C.; Botmart, T.; El-Morshedy, M. Wiener process effects on the solutions of the fractional (2 + 1)-dimensional Heisenberg ferromagnetic spin chain equation. Mathematics 2022, 10, 2043. [Google Scholar] [CrossRef]
  41. Mohammed, W.W.; Al-Askar, F.M.; Cesarano, C.; El-Morshedy, M. The Optical Solutions of the Stochastic Fractional Kundu–Mukherjee–Naskar Model by Two Different Methods. Mathematics 2022, 10, 1465. [Google Scholar] [CrossRef]
  42. Mohammed, W.W.; Alshammari, M.; Cesarano, C.; Albadrani, S.; El-Morshedy, M. Brownian Motion Effects on the Stabilization of Stochastic Solutions to Fractional Diffusion Equations with Polynomials. Mathematics 2022, 10, 1458. [Google Scholar] [CrossRef]
  43. Elbrolosy, M.; Elmandouh, A. Dynamical behaviour of nondissipative double dispersive microstrain wave in the microstructured solids. Eur. Phys. J. Plus 2021, 136, 955. [Google Scholar] [CrossRef]
  44. Elmandouh, A. Integrability, qualitative analysis and the dynamics of wave solutions for Biswas–Milovic equation. Eur. Phys. J. Plus 2021, 136, 638. [Google Scholar] [CrossRef]
  45. Elbrolosy, M. Qualitative analysis and new soliton solutions for the coupled nonlinear Schrödinger type equations. Phys. Scr. 2021, 96, 125275. [Google Scholar] [CrossRef]
  46. Elbrolosy, M.; Elmandouh, A.; Elmandouh, A. Construction of new traveling wave solutions for the (2 + 1) dimensional extended Kadomtsev-Petviashvili equation. J. Appl. Anal. Comput. 2022, 12, 533–550. [Google Scholar] [CrossRef]
  47. Elmandouh, A.A.; Elbrolosy, M.E. New traveling wave solutions for Gilson–Pickering equation in plasma via bifurcation analysis and direct method. Math. Methods Appl. Sci. 2022. [Google Scholar] [CrossRef]
  48. Elbrolosy, M.; Elmandouh, A. Dynamical behaviour of conformable time-fractional coupled Konno-Oono equation in magnetic field. Math. Probl. Eng. 2022, 2022, 3157217. [Google Scholar] [CrossRef]
  49. Elmandouh, A.; Elbrolosy, M. Integrability, variational principal, bifurcation and new wave solutions for Ivancevic option pricing model. J. Math. 2022, 2, 3. [Google Scholar]
  50. Liu, Z.F.; Zhang, S.Y. Solitary waves in finite deformation elastic circular rod. Appl. Math. Mech. 2006, 27, 1255–1260. [Google Scholar] [CrossRef]
  51. Liu, S.K.; Liu, S. Nonlinear Equations in Physics; Peking University Press: Beijing, China, 2000. [Google Scholar]
  52. Porubov, A.V.; Velarde, M.G. On nonlinear waves in an elastic solid. Comptes Rendus Acad. Sci.-Ser. IIB-Mech. 2000, 328, 165–170. [Google Scholar] [CrossRef]
  53. Zhang, S.; Guo, J.; Zhang, N. The dynamics behaviors and wave properties of finite deformation elastic rods with viscous or geometrical-dispersive effects. In Proceedings of the Forth International Conference on Nonlinear Mechanics—ICNM-IV, Shanghai, China, 14–17 August 2002; pp. 728–732. [Google Scholar]
  54. Fu, Z.; Liu, S.; Liu, S. New transformations and new approach to find exact solutions to nonlinear equations. Phys. Lett. A 2002, 299, 507–512. [Google Scholar] [CrossRef]
  55. Goldstein, H. Classical Mechanics; Addison-Wesley Series in Physics; Addison-Wesley: Boston, MA, USA, 1980. [Google Scholar]
  56. Saha, A.; Banerjee, S. Dynamical Systems and Nonlinear Waves in Plasmas; CRC Press: Boca Raton, FL, USA, 2021; pp. x+207. [Google Scholar]
  57. Yang, L.; Hou, X.; Zeng, Z. A complete discrimination system for polynomials. Sci. China Ser. E 1996, 39, 628–646. [Google Scholar]
  58. Meng, Q.; He, B.; Long, Y.; Rui, W. Bifurcations of travelling wave solutions for a general Sine–Gordon equation. Chaos Solitons Fractals 2006, 29, 483–489. [Google Scholar] [CrossRef]
  59. Bin, H.; Jibin, L.; Yao, L.; Weiguo, R. Bifurcations of travelling wave solutions for a variant of Camassa–Holm equation. Nonlinear Anal. Real World Appl. 2008, 9, 222–232. [Google Scholar] [CrossRef]
  60. Byrd, P.F.; Friedman, M.D. Handbook of Elliptic Integrals for Engineers and Scientists, 2nd ed., revised; Die Grundlehren der Mathematischen Wissenschaften, Band 67; Springer: New York, NY, USA; Berlin/Heidelberg, Germany, 1971; pp. xvi+358. [Google Scholar]
Figure 1. Phase plane orbits for the system (6) if Δ < 0 . Equilibria are shown by black solid circles; (a) a = 1 ,   c = 3 ; (b) a = −1, c = −3.
Figure 1. Phase plane orbits for the system (6) if Δ < 0 . Equilibria are shown by black solid circles; (a) a = 1 ,   c = 3 ; (b) a = −1, c = −3.
Symmetry 15 00650 g001
Figure 2. Phase plane orbits for the system (6) if Δ = 0 . The equilibria are indicated by the black solid circles; (a) a = 2 ,   c = 2 , (b) a = −2, c = −2.
Figure 2. Phase plane orbits for the system (6) if Δ = 0 . The equilibria are indicated by the black solid circles; (a) a = 2 ,   c = 2 , (b) a = −2, c = −2.
Symmetry 15 00650 g002
Figure 5. Graphic representation for the periodic solution (30) and its corresponding orbit; (a) Periodic orbit; (b) 3D-periodic solution.
Figure 5. Graphic representation for the periodic solution (30) and its corresponding orbit; (a) Periodic orbit; (b) 3D-periodic solution.
Symmetry 15 00650 g005
Figure 6. The influence of parameters c 0 , c 1 , ν and w on the periodic solutions; (a) Influence of c 0 ; (b) Influence of c 1 ; (c) Influence of ν ; (d) Influence of w.
Figure 6. The influence of parameters c 0 , c 1 , ν and w on the periodic solutions; (a) Influence of c 0 ; (b) Influence of c 1 ; (c) Influence of ν ; (d) Influence of w.
Symmetry 15 00650 g006
Figure 7. Graphic representation for the super-periodic solution (31) and its corresponding orbit; (a) Super-periodic orbit; (b) 3D-super-periodic solution.
Figure 7. Graphic representation for the super-periodic solution (31) and its corresponding orbit; (a) Super-periodic orbit; (b) 3D-super-periodic solution.
Symmetry 15 00650 g007
Figure 8. The influence of parameters c 0 , c 1 , ν and w on the super-periodic solutions; (a) Influence of c 0 ; (b) Influence of c 1 ; (c) Influence of ν ; (d) Influence of w.
Figure 8. The influence of parameters c 0 , c 1 , ν and w on the super-periodic solutions; (a) Influence of c 0 ; (b) Influence of c 1 ; (c) Influence of ν ; (d) Influence of w.
Symmetry 15 00650 g008
Figure 9. Graphic representation for the solitary solutions (28) and (29) and their corresponding orbit; (a) Homoclinic orbit; (b) rarefactive solitary; (c) Compressive solitary.
Figure 9. Graphic representation for the solitary solutions (28) and (29) and their corresponding orbit; (a) Homoclinic orbit; (b) rarefactive solitary; (c) Compressive solitary.
Symmetry 15 00650 g009
Figure 10. The influence of parameters c 0 , c 1 , ν and w on compressive wave solitary solution; (a) Influence of c 0 ; (b) Influence of c 1 ; (c) Influence of ν ; (d) Influence of w.
Figure 10. The influence of parameters c 0 , c 1 , ν and w on compressive wave solitary solution; (a) Influence of c 0 ; (b) Influence of c 1 ; (c) Influence of ν ; (d) Influence of w.
Symmetry 15 00650 g010
Figure 11. 2D and 3D phase portrait for the perturbed system (34) when a = c = 2 and initial condition u ( 0 ) = 0.5 ,   z ( 0 ) = 1 for fixed values of the parameters a , c and various values of α ; (a) α = 0 ; (b) α = 2 ; (c) α = 5 ; (d) α = 0 ; (e) α = 2 ; (f) α = 5 .
Figure 11. 2D and 3D phase portrait for the perturbed system (34) when a = c = 2 and initial condition u ( 0 ) = 0.5 ,   z ( 0 ) = 1 for fixed values of the parameters a , c and various values of α ; (a) α = 0 ; (b) α = 2 ; (c) α = 5 ; (d) α = 0 ; (e) α = 2 ; (f) α = 5 .
Symmetry 15 00650 g011
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

Aljuaidan, A.; Elbrolosy, M.; Elmandouh, A. Nonlinear Wave Propagation for a Strain Wave Equation of a Flexible Rod with Finite Deformation. Symmetry 2023, 15, 650. https://doi.org/10.3390/sym15030650

AMA Style

Aljuaidan A, Elbrolosy M, Elmandouh A. Nonlinear Wave Propagation for a Strain Wave Equation of a Flexible Rod with Finite Deformation. Symmetry. 2023; 15(3):650. https://doi.org/10.3390/sym15030650

Chicago/Turabian Style

Aljuaidan, Aqilah, Mamdouh Elbrolosy, and Adel Elmandouh. 2023. "Nonlinear Wave Propagation for a Strain Wave Equation of a Flexible Rod with Finite Deformation" Symmetry 15, no. 3: 650. https://doi.org/10.3390/sym15030650

APA Style

Aljuaidan, A., Elbrolosy, M., & Elmandouh, A. (2023). Nonlinear Wave Propagation for a Strain Wave Equation of a Flexible Rod with Finite Deformation. Symmetry, 15(3), 650. https://doi.org/10.3390/sym15030650

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