Next Article in Journal
Oscillation and Asymptotic Properties of Differential Equations of Third-Order
Next Article in Special Issue
An Improved Tikhonov-Regularized Variable Projection Algorithm for Separable Nonlinear Least Squares
Previous Article in Journal
Advances in the Theory of Compact Groups and Pro-Lie Groups in the Last Quarter Century
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Periodic Property and Instability of a Rotating Pendulum System

1
School of Science, Xi’an University of Architecture & Technology, Xi’an 710055, China
2
School of Mathematics and Information Science, Henan Polytechnic University, Jiaozuo 454000, China
3
National Engineering Laboratory for Modern Silk, College of Textile and Clothing Engineering, Soochow University, 199 Ren-Ai Road, Suzhou 215000, China
4
Mathematics Department, Faculty of Science, Tanta University, Tanta 31527, Egypt
5
Kafrelsheikh Higher Institute of Engineering and Technology, Kafrelsheikh 33511, Egypt
6
Department of Engineering Physics and Mathematics, Faculty of Engineering, Tanta University, Tanta 31734, Egypt
*
Author to whom correspondence should be addressed.
Axioms 2021, 10(3), 191; https://doi.org/10.3390/axioms10030191
Submission received: 12 July 2021 / Revised: 14 August 2021 / Accepted: 16 August 2021 / Published: 18 August 2021

Abstract

:
The current paper investigates the dynamical property of a pendulum attached to a rotating rigid frame with a constant angular velocity about the vertical axis passing to the pivot point of the pendulum. He’s homotopy perturbation method is used to obtain the analytic solution of the governing nonlinear differential equation of motion. The fourth-order Runge-Kutta method (RKM) and He’s frequency formulation are used to verify the high accuracy of the obtained solution. The stability condition of the motion is examined and discussed. Some plots of the time histories of the gained solutions are portrayed graphically to reveal the impact of the distinct parameters on the dynamical motion.

1. Introduction

It is known that many engineering problems can be formulated by nonlinear ordinary or partial differential equations. With an exception of few problems, their exact solutions seem to be extremely complex and sometimes unreachable. Therefore, asymptotic solutions have shed the interest of many scientists to deal with various nonlinear equations, such as the averaging method and the small parameter method for some weak nonlinear problems [1,2,3,4]. In addition, the multiple scales (MS) method and the Lindstedt-Poincaré (LP) method have great advantages in obtaining the solutions of vibratory systems [5,6]. However, these methods depend on a small parameter, and improper selection of this parameter leads to wrong solutions.
On the other side, the homotopy perturbation method (HPM), that goes back to Ji-Huan He [7], doesn’t depend on a small parameter and it can transform a non-linear problem to a limited number of linear ones which are easy to be solved analytically. In [7], the famous Lighthill equation and Duffing equations were solved. In [8], the method was found to be powerful to fractional differential equations. In [9], this method was applied to solve the nonlinear damped equation of Mathieu with periodic coefficients, and the behavior of stability at both cases of resonance and non-resonance were studied. In [10], the solution of the dynamical motion of a vibrating system was obtained using HPM. This system consists of two masses, one of them attached with a fixed spring and it moves horizontally. The second one relates to the first mass with a massless string and moves vertically. The stability of the motion was examined and discussed. In [11], the motion of a rocking uniform rigid rod on a circular surface was investigated. The approximate solution was obtained applying the HPM and Laplace transform in which the stability conditions were also obtained. This problem was studied previously in [12] using the method of variational approach, the comparison shows that HMP is very accurate, and it is easy to use. The motion of a strong nonlinear system was investigated in [13]. The author obtained the approximate solution using a combination of multiple scales method and the homotopy perturbation method. The stability of an excited delayed Mathieu equation using the He-multiple-scales perturbation method [14] was investigated in [15]. HPM was utilized in [16] to deal with a cubic nonlinearity problem of a conservative couple mass-spring system dynamical system in which the periodic solutions are obtained. This problem was examined also in [17] through the development of an iteration technique based on the method of Mickens iteration to get the asymptotic angular frequencies. In [18,19] HMP was applied to study of the pull-in instability of N/MEMS systems. In [20], HPM-based dynamic analysis was proposed. In [21] the method was extended to solve fractional evolution equation. So far there were many modifications of the homotopy perturbation method, see for examples, Refs. [22,23,24,25,26,27,28,29,30], and have many advantages over others in open literature [31,32,33,34,35].
In this work, the HPM is applied to obtain an asymptotic solution of the controlling equation of motion of a simple pendulum fixed in a rotating rigid frame with constant angular velocity. The numerical solution of this equation is obtained applying RKM from fourth-order. A comparison between them, through some tables and their corresponding figures, emphasizes the accuracy of the HPM. The stability condition of the motion is obtained and discussed.

2. Problem’s Description

The aim of this section is to derive the equation of motion of a simple pendulum of an arm r fixed from one end in a rotating axis represented by a rigid rod and the other end is attached with a mass m . To visualize the motion, we consider two Cartesian systems of coordinates in which the first one O x y z is fixed in space and the other O x y z is fixed in the body and rotates with it, see Figure 1. The rod is connected with a rotating rigid frame with constant angular velocity Ω = φ ˙ about the vertical axes O y and O y . Here, φ is the angle of rotation of the rotary frame O x y z about the vertical axis O y .
Therefore, we can write
z = z cos φ x sin φ , x = x cos φ + z sin φ ,
where
x = r sin θ , y = y = l r cos θ , z = 0 ,
where θ is the angle of inclination of the T-shaped with the vertical axis and l is the distance from the rotating axis (rigid rod) of the pendulum to z -axis. Based on the above, one obtains easily
x = r sin θ cos φ , y = l r cos θ , z = r sin θ sin φ .
Making use of the above projections of the point p on the coordinates system O x y z to write the kinetic and the potential energies in the form
T = 1 2 m r 2 ( θ ˙ 2 + Ω 2 sin 2 θ ) , V = m g ( l r cos θ ) ,
where dots denote the derivative with respect to time t and g is the gravitational acceleration.
According to the variational theory [36,37,38], the Lagrange’s equation for conservative dynamical systems is
d d t ( L θ ˙ ) L θ = 0 ,
where L = T V is the Lagrangian, the equation of motion (EOM) has therefore the form
r g θ ¨ + ( 1 Λ cos θ ) sin θ = 0 ,
where Λ = Ω 2 r / g .
Now, let us introduce a new independent variable τ in the form
τ = g r t .
According to (7), we can rewrite the EOM (6) in the form
θ + ( 1 Λ cos θ ) sin θ = 0 ,
where prims denote the derivative with respect to time τ .

3. The Homotopy Perturbation Method

In this section, we outline on the HPM through consideration of the following nonlinear equation
K ( θ ) f ( r ) = 0 , r Ω * ,
with the boundary condition
B ( θ , θ n ) = 0 , r Γ .
Here K , B , f and Γ are a general differential operator, a boundary one, analytical function, and the boundary of a domain Ω * . Moreover, θ n represent the differential along the normal drawn outwards from Ω * .
According to HPM, we can separate the operator K into linear and nonlinear parts L and N respectively. Consequently, one can rewrite Equation (9) in the form
L ( θ ) + N ( θ ) f ( r ) = 0
An inspection of Equation (10) shows that we can formulate a homotopy of (10) to satisfy θ ( r , ρ ) : Ω * × [ 0 , 1 ] R ,
H ( θ , ρ ) = L ( θ ) L ( U ) + ρ L ( U ) + ρ [ N ( θ ) f ( r ) ] = 0 , ρ [ 0 , 1 ]
where ρ is an embedding parameter and U is an initial approximation guess of Equation (9), in which the boundary conditions are fulfilled.
Based on the HPM, the solution of (12) can be expanded into a power series of ρ as
θ = θ 0 + ρ θ 1 + ρ 2 θ 2 + .
When ρ 1 , will find that Equation (12) matches with Equation (9) and therefore, we can express the asymptotic solution of Equation (9) in the form
v = lim ρ 1 θ = θ 0 + θ 1 + θ 2 + .
It is noteworthy that, in many cases, the series (14) is convergent. For the convergence of this sequence, some certain conditions are proposed [7].

4. Method of Solution

The purpose of this section is to employ the HPM to obtain the asymptotic solution of the EOM. Substituting the following expansion of the trigonometric function sin θ into the EOM (8)
sin θ = θ θ 3 3 ! + ,
to get its approximation form
θ + ( 1 Λ ) θ + ( 4 Λ 1 3 ! ) θ 3 + = 0 .
Let’s consider the initial conditions of the above equation in the form
θ ( 0 ) = 1 , θ ( 0 ) = 0 .
According to (11), we can express the linear and nonlinear parts L ( θ ) and N ( θ ) of Equation (16) as follows
L ( θ ) = θ + ( 1 Λ ) θ ,
N ( θ ) = ( 4 Λ 1 3 ! ) θ 3 ,
where
f ( t ) = 0 .
It is possible to rewrite Equation (12) in the form
H ( θ , ρ ) = ( 1 ρ ) [ L ( θ ) L ( U ) ] + ρ [ L ( θ ) + N ( θ ) f ( t ) ] = 0 , ρ [ 0 , 1 ]
The substituting of (18)–(20) into (21) yields
H ( θ , ρ ) = ( 1 ρ ) [ θ + ( 1 Λ ) θ U ( 1 Λ ) U ] + ρ [ θ + ( 1 Λ ) θ + ( 4 Λ 1 3 ! ) θ 3 ] = 0 ,
where U is the initial approximation guess. It does not depend on the boundary conditions. So that we may choose a specific value or function to make the solution easier and shorter.
For simplicity we set U = 0 , Equation (22) has the following form
H ( θ , ρ ) = θ + ( 1 Λ ) θ + ρ ( 4 Λ 1 3 ! ) θ 3 = 0 .
Substituting the power series (13) into (23), and equating the coefficients of different powers of ρ with zero, we get
Coefficients of ρ 0 :
θ 0 + ( 1 Λ ) θ 0 = 0 .
Coefficients of ρ 1 :
θ 1 + ( 1 Λ ) θ 1 + ( 4 Λ 1 6 ) θ 0 3 = 0 .
Coefficients of ρ 2 :
θ 2 + ( 1 Λ ) θ 2 + ( 4 Λ 1 2 ) θ 1 θ 0 2 = 0 .
The inspection of the above Equations (24)–(26) shows that we can solve them sequentially with the use of the following initial conditions
θ 0 ( 0 ) = 1 , θ 0 ( 0 ) = 0 , θ 1 ( 0 ) = 0 , θ 1 ( 0 ) = 0 , θ 2 ( 0 ) = 0 , θ 2 ( 0 ) = 0 ,
to obtain
θ 0 = cos ( τ 1 Λ ) ,
θ 1 = [ 4 ( 1 Λ ) 3 ] 96 ( 1 Λ ) [ 6 1 Λ τ + sin ( 2 1 Λ τ ) ] sin ( 1 Λ τ ) ,
θ 2 = [ 3 4 ( 1 Λ ) ] 2 36864 ( 1 Λ ) 2 { [ 72 ( 1 Λ ) τ 2 23 ] cos ( 1 Λ τ ) + 24 cos ( 3 1 Λ τ ) cos ( 5 1 Λ τ ) + 12 τ 1 Λ [ 3 sin ( 3 1 Λ τ ) 8 sin ( 1 Λ τ ) ] }
Substituting (28)–(30) into (13) to obtain
θ = cos ( 1 Λ τ ) + ρ [ [ 4 ( 1 Λ ) 3 ] 96 ( 1 Λ ) [ 6 1 Λ τ + sin ( 2 1 Λ τ ) ] sin ( 1 Λ τ ) ] + ρ 2 { [ 3 4 ( 1 Λ ) ] 2 36864 ( 1 Λ ) 2 { [ 72 τ 2 ( 1 Λ ) 23 ] cos ( 1 Λ τ ) + 24 cos ( 3 1 Λ τ ) cos ( 5 1 Λ τ ) + 12 τ 1 Λ [ 3 sin ( 3 1 Λ τ ) 8 sin ( 1 Λ τ ) ] } } .
According to (14), the above approximate solution has the form
v = lim ρ 1 θ = cos ( 1 Λ τ ) + [ 4 ( 1 Λ ) 3 ] 96 ( 1 Λ ) [ 6 1 Λ τ + sin ( 2 1 Λ τ ) ] sin ( 1 Λ τ ) [ 3 4 ( 1 Λ ) ] 2 36864 ( 1 Λ ) 2 { [ 72 τ 2 ( 1 Λ ) 23 ] cos ( 1 Λ τ ) + 24 cos ( 3 1 Λ τ ) cos ( 5 1 Λ τ ) + 12 τ 1 Λ [ 3 sin ( 3 1 Λ τ ) 8 sin ( 1 Λ τ ) } ] .
Numerically, the successive approximations by the HPM (with a finite series) are guaranteed to converge to the exact solution over some intervals. Therefore, v ( τ ) converges locally uniform.

5. Stability Analysis

The main objective of this section is to study the stability of the considered dynamical model by examining its EOM (16), in which it isn’t useful to study stability through the obtained solution (32). Therefore, we are going to consider the linear and nonlinear parts of Equation (16) in which they are represented by Equations (18) and (19) respectively. It is worthy to mention that, the stability of the linear part depends on the frequency term ( 1 Λ ) which is always positive. Therefore, this term can be expanded in a power series of ρ as follows
( 1 Λ ) = ω 2 + ρ ω 1 + ρ 2 ω 2 +
Substituting from Equation (13) about θ and from (33) about the expanded frequency into Equation (16) and equating the coefficients of equal powers of ρ for both sides of the resulted equation to get
Coefficient of ρ 0 :
θ 0 + ω 2 θ 0 = 0 ,
Coefficient of ρ :
θ 1 + ω 2 θ 1 + ω 1 θ 0 + 1 6 ( 4 Λ 1 ) θ 0 3 = 0 ,
Taking into account the previous initial conditions (27), we can write the solution of the homogenous Equation (34) in the form
θ 0 = cos ( τ ω ) .
Therefore Equation (35) becomes
θ 1 + ω 2 θ 1 = ω 1 cos ( τ ω ) 1 6 ( 4 Λ 1 ) cos 3 ( τ ω ) .
Elimination of the secular terms demands that
ω 1 = 1 8 ( 4 Λ 1 ) .
Considering the initial conditions θ 1 ( 0 ) = θ 1 ( 0 ) = 0 , we can write the solution of (37) after elimination of the secular term in the form
θ 1 = ( 4 Λ 1 ) 24 ( 7 4 Λ ) [ cos ( 3 τ ω ) cos ( τ ω ) ] .
Making use of (38) into (33) and considering ρ = 1 , one can write the obtained frequency in the form
ω 2 = 1 8 ( 7 4 Λ ) .
In order to keep the system stable, we must consider the following stability condition
Λ < 7 4 .

6. He’s Frequency Formulation

In order to verify our above results, this section introduces briefly He’s frequency formulation [38,39]. Considering a nonlinear oscillator in the form
θ + g ( θ ) = 0
with initial conditions
θ ( 0 ) = A , θ ( 0 ) = 0
where g is a smooth function. Equation (42) has periodic solution when
d 2 d θ 2 g ( θ ) > 0
He’s frequency formulation is [38,39]
ω 2 = d d θ g ( θ ) | θ = A / 2
In our study,
g ( θ ) = ( 1 Λ ) θ + ( 4 Λ 1 3 ! ) θ 3
It is easy to find that
d d θ g ( θ ) = 1 Λ + ( 4 Λ 1 2 ) θ 2
In our study A = 1, according to Equation (45), we obtain
ω 2 = 1 Λ + 4 Λ 1 8 = 7 4 Λ 8
This is exactly same as that given in Equation (40). Applications of He’s frequency formulation are referred to refs [40,41,42,43,44,45].

7. Results and Discussion

This section sheds light on the great accuracy of the obtained results that are achieved by using HPM through the comparison of these results with the numerical ones that are gained by utilizing the fourth-order RKM [46,47,48].
A beneficial way for a good comparison between the attained asymptotic results by HPM and the numerical ones obtained by RKM, is to look reviews them through the Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6 in addition to the error between them. The results included in Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6 correspond to the curves of Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 respectively for the same corresponding values of r and Ω .
The curves displayed Figure 2, Figure 3 and Figure 4 are calculated when r = 0.6 m with the distinct values of Ω = ( 1.5 , 2 , 2.5 ) rad · s 1 while Figure 6 and Figure 7 are plotted when Ω = 2 at r = 0.6 . The objective of these curves is to reveal the comparison between the approximate analytical solutions that are represented in Equation (32) (with red color) and the numerical solutions of the governing EOM (8) (with blue color). These drawings indicate that the comparison between both results reaches a peak of congruence at r = 0.6 when Ω = 1.5 and Ω = 2 as seen in Figure 2 and Figure 3 respectively. On the other side, this comparison is not completely consistent when r = 0.6 and Ω = 2.5 for the attained solutions after the elapse of half period time, as shown in Figure 4. A closer look at these figures shows that the plotted curves have aperiodic behavior, which confirms the stability of the obtained solutions.
The purpose of the graphically generated results Figure 8 and Figure 9 is to investigate the impact of different values of Ω = ( 1.2 , 2 , 2.5 ) rad · s 1 and r ( = 0.3 , 0.6 , 0.9 ) m respectively, with the constancy of r = 0.6 in Figure 8 and Ω = 2 in Figure 9 on the behavior of the considered dynamical model. It is clear that when Ω and r increase, the amplitudes of the waves increase to some extent besides the constancy of the oscillations number and wavelengths.
The phase plane diagrams that assert the stability of the attained solution, at different values of Ω and r , are represented graphically in Figure 10 and Figure 11 respectively.

8. Conclusions

The asymptotic periodic solution of the EOM of a simple pendulum fixed in a rotating rigid frame is obtained using HPM. The numerical solution of the governing EOM is achieved utilizing the fourth-order RKM. A comparison between the attained solutions, whether analytical or numerical, showed a clear match between them which emphasizes the accuracy of the used HPM. These solutions are performed through computer codes to represent the time histories of the motion graphically at the distinct values of the physical parameters of the studied model. The stability condition of the motion is obtained.

Author Contributions

Conceptualization, J.-H.H., A.A.G.; methodology, J.-H.H., A.A.G.; software, T.S.A., S.E.; validation, J.-H.H., A.A.G.; formal analysis, J.-H.H.; investigation, J.-H.H., T.S.A., S.E., A.A.G.; writing—original draft preparation, J.-H.H., A.A.G.; writing—review and editing, J.-H.H., T.S.A., S.E., A.A.G.; visualization, T.S.A., S.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Ji, W.M.; Wang, H.; Liu, M. Dynamics analysis of an impulsive stochastic model for spruce budworm growth. Appl. Comput. Math. 2021, 19, 336–359. [Google Scholar]
  2. Janevski, G.; Kozic, P.; Pavlovic, R.; Posavljak, S. Moment Lyapunov exponents and stochastic stability of a thin-walled beam subjected to axial loads and end moments. Facta Univ. Ser. Mech. Eng. 2021, 19, 209–228. [Google Scholar] [CrossRef]
  3. Pavlovic, I.R.; Pavlovic, R.; Janevski, G.; Despenić, N.; Pajković, V. Dynamic behavior of two elastically connected nanobeams under a white noise process. Facta Univ. Ser. Mech. Eng. 2020, 18, 219–227. [Google Scholar] [CrossRef]
  4. Zuo, Y.-T. A gecko-like fractal receptor of a three-dimensional printing technology: A fractal oscillator. J. Math. Chem. 2021, 59, 735–744. [Google Scholar] [CrossRef]
  5. Yeasmin, I.A.; Rahman, M.S.; Alam, M.S. The modified Lindstedt-Poincare method for solving quadratic nonlinear oscillators. J. Low Freq. Noise Vib. Act. Control 2020, 1461348420979758. [Google Scholar] [CrossRef]
  6. El-Sabaa, F.M.; Amer, T.S.; Gad, H.M.; Bek, M.A. On the motion of a damped rigid body near resonances under the influence of harmonically external force and moments. Results Phys. 2020, 19, 103352. [Google Scholar] [CrossRef]
  7. He, J.H. Homotopy perturbation technique. Comput. Meth. Appl. Mech. Eng. 1999, 178, 257–262. [Google Scholar] [CrossRef]
  8. Anjum, N.; Ain, Q.T. Application of He’s fractional derivative and fractional complex transform for time fractional Camassa-Holm equation. Therm. Sci. 2020, 24, 3023–3030. [Google Scholar] [CrossRef]
  9. El-Dib, Y.O. Homotopy perturbation for excited nonlinear equations. Sci. Eng. Appl. 2017, 2, 96–108. [Google Scholar] [CrossRef]
  10. Amer, T.S.; Galal, A.A.; Elnaggar, S. The vibrational motion of a dynamical system using homotopy perturbation technique. Appl. Math. 2020, 11, 1081–1099. [Google Scholar] [CrossRef]
  11. El-Dib, Y.O.; Moatimid, G.M. Stability configuration of a rocking rigid rod over a circular surface using the homotopy perturbation method and Laplace transform. Arab. J. Sci. Eng. 2019, 44, 6581–6591. [Google Scholar] [CrossRef]
  12. Ganji, S.S.; Ganji, D.D.; Babazadeh, H.; Karimpour, S. Variational approach method for nonlinear oscillations of the motion of a rigid rocking back and forth and cubic-quintic Duffing oscillators. Prog. Electromagn. Res. M 2008, 4, 23–32. [Google Scholar] [CrossRef] [Green Version]
  13. El-Dib, Y.O. Periodic solution and stability behavior for nonlinear oscillator having a cubic nonlinearity time-delayed. Int. Annu. Sci. 2018, 5, 12–25. [Google Scholar] [CrossRef] [Green Version]
  14. El-Dib, Y.O. Periodic solution of the cubic nonlinear Klein-Gordon equation and the stability criteria via the He-multiple-scales method. Pramana J. Phys. 2019, 92, 7. [Google Scholar] [CrossRef]
  15. El-Dib, Y.O. Stability approach for periodic delay Mathieu equation by the He-multiple-scales method. Alex. Eng. J. 2018, 57, 4009–4020. [Google Scholar] [CrossRef]
  16. Hosen, M.A.; Chowdhury, M.S.H. Accurate approximations of the nonlinear vibration of couple-mass-springs systems with linear and nonlinear stiffnesses. J. Low Freq. Noise Vib. Act. Control 2021, 40, 1072–1090. [Google Scholar] [CrossRef] [Green Version]
  17. Hosen, M.A. Analysis of nonlinear vibration of couple-mass–spring systems using iteration technique. Multidiscip. Modeling Mater. Struct. 2020, 16, 1539–1558. [Google Scholar] [CrossRef]
  18. Anjum, N.; He, J.H. Homotopy perturbation method for N/MEMS oscillators. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  19. Tian, D.; Ain, Q.-T.; Anjum, N.; He, C.-H.; Cheng, B. Fractal N/MEMS: From pull-in instability to pull-in stability. Fractals 2021, 29, 2150030. [Google Scholar] [CrossRef]
  20. Meshki, H.; Rezaei, A.; Sadeghi, A. Homotopy Perturbation-Based Dynamic Analysis of Structural Systems. J. Eng. Mech. 2020, 146, 04020136. [Google Scholar] [CrossRef]
  21. Wang, K.L.; Yao, S.W. He’s fractiona; derivative for the evolution equation. Therm. Sci. 2020, 24, 2507–2513. [Google Scholar] [CrossRef]
  22. He, J.H.; El-Dib, Y.O. The enhanced homotopy perturbation method for axial vibration of strings. Facta Univ. Ser. Mech. Eng. 2021. [Google Scholar] [CrossRef]
  23. He, C.H.; Tian, D.; Moatimid, G.M.; Salman, H.F.; Zekry, M.H. Hybrid Rayleigh-Van der Pol-Duffing Oscillator (HRVD): Stability Analysis and Controller. J. Low Freq. Noise Vib. Act. Control 2021. [Google Scholar] [CrossRef]
  24. He, J.-H.; Galal, M.M.; Mostapha, D.R. Nonlinear Instability of Two Streaming-Superposed Magnetic Reiner-Rivlin Fluids by He-Laplace Method. J. Electroanal. Chem. 2021, 895, 115388. [Google Scholar] [CrossRef]
  25. Anjum, N.; He, J.H.; Ain, Q.T.; Tian, D. Li-He’s modified homotopy perturbation method for doubly-clamped electrically actuated microbeams-based microelectromechanical system. Facta Univ. Ser. Mech. Eng. 2021. [Google Scholar] [CrossRef]
  26. El-Dib, Y.O.; Matoog, R.T. The Rank Upgrading Technique for a Harmonic Restoring Force of Nonlinear Oscillators. J. Appl. Comput. Mech. 2021, 7, 782–789. [Google Scholar] [CrossRef]
  27. Elgazery, N.S. A Periodic Solution of the Newell-Whitehead-Segel (NWS) Wave Equation via Fractional Calculus. J. Appl. Comput. Mech. 2020, 6, 1293–1300. [Google Scholar]
  28. Koochi, A.; Goharimanesh, M. Nonlinear Oscillations of CNT Nano-resonator Based on Nonlocal Elasticity: The Energy Balance Method. Rep. Mech. Eng. 2021, 2, 41–50. [Google Scholar] [CrossRef]
  29. Filobello-Nino, U.; Vazquez-Leal, H.; Herrera-May, A.; Ambrosio, R.; Jimenez-Fernandez, V.M.; Sandoval-Hernandez, M.A.; Alvarez-Gasca, O.; Palma-Grayeb, B.E. The study of heat transfer phenomena by using modified homotopy perturbation method coupled by Laplace transform. Therm. Sci. 2020, 24, 1105–1115. [Google Scholar] [CrossRef] [Green Version]
  30. Li, X.X.; He, C.H. Homotopy perturbation method coupled with the enhanced perturbation method. J. Low Freq. Noise Vib. Act. Control 2019, 38, 1399–1403. [Google Scholar] [CrossRef]
  31. Mahmudov, N.I.; Huseynov, I.T.; Aliev, N.A.; Aliev, F.A. Analytical approach to a class of Bagley-Torvik equations. TWMS J. Pure Appl. Math. 2020, 11, 238–258. [Google Scholar]
  32. Qalandarov, A.A.; Khaldjigitov, A.A. Mathematical and Numerical modeling of the coupled dynamic thermoelastic problems for isotropic bodies. TWMS J. Pure Appl. Math. 2020, 11, 119–126. [Google Scholar]
  33. Fikret, A.A.; Aliev, N.A.; Mutallimov, M.M.; Namazov, A.A. Algorithm for solving the identification problem for determining the fractional-order derivative of an oscillatory system. Appl. Comput. Math. 2020, 19, 435–442. [Google Scholar]
  34. Zhang, J.J.; Shen, Y.; He, J.H. Some analytical methods for singular boundary value problem in a fractal space: A review. Appl. Comput. Math. 2019, 18, 225–235. [Google Scholar]
  35. Song, H.Y. A thermodynamic model for a packing dynamic system. Therm. Sci. 2020, 24, 2331–2335. [Google Scholar] [CrossRef]
  36. Yao, S.W. Variational principle for non-linear fractional wave equation in a fractal space. Therm. Sci. 2021, 25, 1243–1247. [Google Scholar] [CrossRef]
  37. Liu, H.Y.; Li, Z.M.; Yao, S.W.; Yao, Y.J.; Liu, J. A variational principle for the photocatalytic Nox abatement. Therm. Sci. 2020, 24, 2515–2518. [Google Scholar] [CrossRef]
  38. He, J.H. The simpler, the better: Analytical methods for nonlinear oscillators and fractional oscillators. J. Low Freq. Noise Vib. Act. Control 2019, 38, 1252–1260. [Google Scholar] [CrossRef] [Green Version]
  39. Qie, N.; Hou, W.F.; He, J.H. The fastest insight into the large amplitude vibration of a string. Rep. Mech. Eng. 2020, 2, 1–5. [Google Scholar] [CrossRef]
  40. Liu, C.X. Periodic solution of fractal Phi-4 equation. Therm. Sci. 2021, 25, 1345–1350. [Google Scholar] [CrossRef]
  41. Liu, C.X. A short remark on He’s frequency formulation. J. Low Freq. Noise Vib. Act. Control 2021, 40, 672–674. [Google Scholar] [CrossRef]
  42. Feng, G.Q.; Niu, J.Y. He’s frequency formulation for nonlinear vibration of a porous foundation with fractal derivative. Int. J. Geomath. 2021, 12, 14. [Google Scholar] [CrossRef]
  43. Feng, G.Q. He’s frequency formula to fractal undamped Duffing equation. J. Low Freq. Noise Vib. Act. Control 2021, 1461348421992608. [Google Scholar] [CrossRef]
  44. Elías-Zúñiga, A.; Palacios-Pineda, L.M.; Jiménez-Cedeño, I.H.; Martínez-Romero, O.; Trejo, D.O. He’s frequency-amplitude formulation for nonlinear oscillators using Jacobi elliptic functions. J. Low Freq. Noise Vib. Act. Control 2020, 39, 1216–1223. [Google Scholar] [CrossRef]
  45. Wu, Y.; Liu, Y.P. Residual calculation in He’s frequency-amplitude formulation. J. Low Freq. Noise Vib. Act. Control 2021, 40, 1040–1047. [Google Scholar] [CrossRef] [Green Version]
  46. Gilat, A. Numerical Methods for Engineers and Scientists; Wiley: Hoboken, NJ, USA, 2013. [Google Scholar]
  47. Simos, T.E.; Tsitouras, C. 6th order Runge-Kutta pairs for scalar autonomous IVP. Appl. Comput. Math. 2020, 19, 412–421. [Google Scholar]
  48. Hussanan, A.; Khan, I.; Khan, W.A.; Chen, Z.M. Micropolar mixed convective flow Cattaneo-Christove heat flux: Non-Fourier Heat Conduction Analysis. Therm. Sci. 2020, 24, 1345–1356. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The dynamical model.
Figure 1. The dynamical model.
Axioms 10 00191 g001
Figure 2. Illustrates the time history of the approximate solutions (red color) and numerical solution (blue color) at r = 0.6 and Ω = 1.5 .
Figure 2. Illustrates the time history of the approximate solutions (red color) and numerical solution (blue color) at r = 0.6 and Ω = 1.5 .
Axioms 10 00191 g002
Figure 3. Shows the time history of the approximate solutions (red color) and numerical solution (blue color) at r = 0.6 and Ω = 2 .
Figure 3. Shows the time history of the approximate solutions (red color) and numerical solution (blue color) at r = 0.6 and Ω = 2 .
Axioms 10 00191 g003
Figure 4. Reveals a comparison between the analytical solution obtained by HPM (red color) and the numerical one obtained by RKM (blue color) r = 0.6 and Ω = 2.5 .
Figure 4. Reveals a comparison between the analytical solution obtained by HPM (red color) and the numerical one obtained by RKM (blue color) r = 0.6 and Ω = 2.5 .
Axioms 10 00191 g004
Figure 5. Shows a comparison between the time histories of the approximate solution (red color) and the numerical one (blue color) Ω = 2 and r = 0.3 .
Figure 5. Shows a comparison between the time histories of the approximate solution (red color) and the numerical one (blue color) Ω = 2 and r = 0.3 .
Axioms 10 00191 g005
Figure 6. Portrays the comparison between the analytic solution (red color) and the numerical one (blue color) Ω = 2 and r = 0.9 .
Figure 6. Portrays the comparison between the analytic solution (red color) and the numerical one (blue color) Ω = 2 and r = 0.9 .
Axioms 10 00191 g006
Figure 7. Shows a comparison between the homotopy solution (red color) and the numerical one (blue color) Ω = 2 and r = 1.1 .
Figure 7. Shows a comparison between the homotopy solution (red color) and the numerical one (blue color) Ω = 2 and r = 1.1 .
Axioms 10 00191 g007
Figure 8. Describes the impact of distinct values of Ω ( = 1.5 , 2 , 2.5 ) on the solution θ at r = 0.6 .
Figure 8. Describes the impact of distinct values of Ω ( = 1.5 , 2 , 2.5 ) on the solution θ at r = 0.6 .
Axioms 10 00191 g008
Figure 9. Explores the effect of the different values of r ( = 0.3 , 0.6 , 0.9 ) on the solution θ at Ω = 2 .
Figure 9. Explores the effect of the different values of r ( = 0.3 , 0.6 , 0.9 ) on the solution θ at Ω = 2 .
Axioms 10 00191 g009
Figure 10. Describes the phase plane of the solution at Ω ( = 1.5 , 2 , 2.5 ) when r = 0.6 .
Figure 10. Describes the phase plane of the solution at Ω ( = 1.5 , 2 , 2.5 ) when r = 0.6 .
Axioms 10 00191 g010
Figure 11. Describes the phase plane of the solution at r ( = 0.3 , 0.6 , 0.9 ) when Ω = 2 .
Figure 11. Describes the phase plane of the solution at r ( = 0.3 , 0.6 , 0.9 ) when Ω = 2 .
Axioms 10 00191 g011
Table 1. Error percentage of HPM for r = 0.6   m , Ω = 1.5   rad · s 1 .
Table 1. Error percentage of HPM for r = 0.6   m , Ω = 1.5   rad · s 1 .
TimeNumerical Results (NR)HPM Results (HPMR) | H P M R N R N R |
0110
10.6308690.6277760.0049026
20.21956−0.2253840.0265278
3−0.899482−0.9031690.00409939
4−0.908171−0.9030950.00558896
5−0.239734−0.224690.0627527
60.6147910.6287340.02268
70.9997961.000710.000911276
80.6466720.6301530.0255458
9−0.199286−0.2250650.129353
10−0.890417−0.9035410.0147397
Table 2. Error percentage of HPM for r = 0.6 m, Ω = 2 rad s−1.
Table 2. Error percentage of HPM for r = 0.6 m, Ω = 2 rad s−1.
TimeNumerical Results (NR)HPM Results (HPMR) | H P M R N R N R |
0110
10.6547640.646950.0119336
2−0.148467−0.163550.101592
3−0.847042−0.8582980.013288
4−0.956805−0.9466170.0106482
5−0.404543−0.3663630.0943759
60.4321580.4731180.0947819
70.9651350.9779070.0132337
80.8305830.7921240.046304
90.1183010.04652040.606763
10−0.677385−0.732010.0806414
Table 3. Error percentage of HPM for r = 0.6 m, Ω = 2.5 rad s−1.
Table 3. Error percentage of HPM for r = 0.6 m, Ω = 2.5 rad s−1.
TimeNumerical Results (NR)HPM Results (HPMR) | H P M R N R N R |
0110
10.6854040.6713870.0204506
20.0527615−0.08087870.532912
3−0.7589940.7838830.0327925
4−0.99427−0.9855330.00878713
5−0.6041640.5384740.108728
60.1577240.2425820.538009
70.8240480.8758860.0629072
80.9771470.9421310.0358349
90.516240.3852230.253791
10−0.261008−0.4041210.548309
Table 4. Error percentage of HPM for r = 0.3 m, Ω = 2 rad s−1.
Table 4. Error percentage of HPM for r = 0.3 m, Ω = 2 rad s−1.
TimeNumerical Results (NR)HPM Results (HPMR) | H P M R N R N R |
0110
10.6274520.6250230.00387119
2−0.229458−0.234060.0200579
3−0.905921−0.9087520.00312526
40.900042−0.8962090.00425873
50.215822−0.2042510.0536163
60.6381760.6488820.0167751
70.9999071.000720.000812474
80.6166030.6052410.0184257
9−0.2430470.2623630.0794747
100.911628−0.9212370.0105405
Table 5. Error percentage of HPM for r = 0.9 m, Ω = 2 rad s−1.
Table 5. Error percentage of HPM for r = 0.9 m, Ω = 2 rad s−1.
TimeNumerical Results (NR)HPM Results (HPMR) | H P M R N R N R |
0110
10.6820040.6686810.0195354
2−0.0636153−0.09027690.419106
30.7700030.7932140.0301438
40.9917210.9823830.0094154
5−0.5829430.5204990.107117
60.1898560.2691910.417868
70.8454180.8917070.0547537
80.9670250.9296640.038635
90.4745080.351220.259822
10−0.313132−0.4431940.415355
Table 6. Error percentage of HPM for r = 1.1 m, Ω = 2 rad s−1.
Table 6. Error percentage of HPM for r = 1.1 m, Ω = 2 rad s−1.
TimeNumerical Results (NR)HPM Results (HPMR) | H P M R N R N R |
0110
10.7001160.683090.0243185
2−0.00515091−0.03935586.64055
3−0.707562−0.7397030.0454255
4−0.9999440.9952520.00469164
5−0.6925950.6111150.117644
60.01545230.1296927.39304
70.7149330.7939760.11056
80.9997750.9838710.0159081
90.6850010.5187030.24277
10−0.0257521−0.2526578.81111
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

He, J.-H.; Amer, T.S.; Elnaggar, S.; Galal, A.A. Periodic Property and Instability of a Rotating Pendulum System. Axioms 2021, 10, 191. https://doi.org/10.3390/axioms10030191

AMA Style

He J-H, Amer TS, Elnaggar S, Galal AA. Periodic Property and Instability of a Rotating Pendulum System. Axioms. 2021; 10(3):191. https://doi.org/10.3390/axioms10030191

Chicago/Turabian Style

He, Ji-Huan, Tarek S. Amer, Shimaa Elnaggar, and Abdallah A. Galal. 2021. "Periodic Property and Instability of a Rotating Pendulum System" Axioms 10, no. 3: 191. https://doi.org/10.3390/axioms10030191

APA Style

He, J. -H., Amer, T. S., Elnaggar, S., & Galal, A. A. (2021). Periodic Property and Instability of a Rotating Pendulum System. Axioms, 10(3), 191. https://doi.org/10.3390/axioms10030191

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