Optimal Auxiliary Functions Method for a Pendulum Wrapping on Two Cylinders

: In the present work, the nonlinear oscillations of a pendulum wrapping on two cylinders is studied by means of a new analytical technique, namely the Optimal Auxiliary Functions Method (OAFM). The equation of motion is derived from the Lagrange’s equation. Analytical solutions and natural frequency of the system are calculated. Our results obtained through this new procedure are compared with numerical ones and a very good agreement was found, which proves the accuracy of the method. The presented numerical examples show that the proposed approach is simple, easy to implement and very accurate.


Introduction
The study of the simple pendulum has a long history.During the Renaissance, Leonardo da Vinci made some drawings related to the motion of pendulum, without realizing at that time its great importance for timekeeping.Beginning around 1602, Galileo Galilei studied for the first time the properties of pendulum, isochronisms, and found that the period of this system is approximately independent of the amplitude or with the swing.Additionally, he demonstrated that the period is proportional to the square root of the length of the pendulum, but independent on the mass.Forty years later he conceived and dictated to his son a design for a pendulum clock.The pendulum was the first harmonic oscillator used by human being [1].In 1673, Huygens discovered that the period of the pendulum is identical, no matter if it hung from its centre of oscillation or from its pivot [2].In 1818 Henry Kater invented the so-named reversible Kater's pendulum, making very accurate measurements of gravity possible.In 1851 Foucault made his investigations known, and a "pendulum mania" broke out [3].Around 1900, the need for higher precision clocks led to the use of low-thermal-expansion materials for pendulum rods.In 1921, the quartz crystal oscillator was invented, and in 1927 quartz clocks replaced pendulum clocks [4].Pendulum gravimeters were replaced by "free fall" gravimeters in the 1950s [5], but pendulum instruments continued to be considered into the 1970s.In 1721, G. Graham [6] invented the mercury pendulum, whose weight is represented by a container of mercury, in which case the pendulum rod gets longer with rising temperature.In 1726, J. Harrison invented the gridiron pendulum, consisting of alternating rods made of different metals, with totally different thermal expansion properties (steel and zinc or brass, respectively).In 1896, C.E. Guillaume invented the nickel-steel alloy [4].The invar pendulum was used for the first time in the Riefler regulator clock, achieving an excellent accuracy.In 1826 G. Airy proved the smallest disturbing effect of the drive force on the period if given as a short type of pendulum, such as the Repsold-Bessel pendulum [7], Van Sterneck and Mendelhall gravimeters, double pendulum gravimeters, Gulf gravimeter [8], and so on.Dynamic mechanical systems possessing the pendulum arise in many domains of activity and many scientists paid attention to obtaining a governing equation of pendulums.The above-mentioned studies were later extended to other types of pendulum with different conditions along their dynamic behavior.
Hamouda and Pierce [9] analyzed the blades of a helicopter rotor (similar to a simple pendulum) to suppress the root reactions.The general nonlinear equations of motion are linearized.They consider the hingeless rotor blade excited by a harmonic variation of span wise air load distribution.Simple flap and lead-lag pendulum are treated individually.The pendulum mass effectiveness was also investigated.
A comprehensive discussion of the corrections needed to accurately measure the acceleration of gravity using a plane pendulum is provided by Nelson and Olson [10].A simple laboratory experiment was described, in which g was determined to four significant figures of accuracy.The effect of the Coriolis force acting on the bob during station is evaluated, adapting a spring-pendulum system analysis to the nearly stiff limit.In their study, the linear and quadratic damped were used and perturbation expansion of the small dimensionless parameter was developed.
Ge and Ku [11] extended the Melnikov approach (which is traditionally restricted to study ingweak non-linear phenomena including sufficient small harmonic excitation) to a pendulum suspended on a rotating arm described by two-dimensional differential equations.These equations possess strongly odd nonlinear function of the displacement and are subjected to large harmonic excitation.
Nester et al. [12] presented an experimental investigation into the dynamic response of rotor systems fitted with centrifugal pendulum vibration absorbers.Two types of absorbers are considered, which exhibit different types of nonlinear behavior.
The spatial double pendulum, comprising two pendulums that swing in different planes is analyzed in [13] by Bendersky and Sandler.Some Mathlab codes were proposed to solve the nonlinear differential equations.The frequency spectra were obtained using Fourier transformation.Solutions of free vibrations and frequency spectra were employed in dynamic investigations for different initial conditions of motion.
A small ellipticity of the driving, perturbing the classical parametric pendulum, was studied by Horton et al. [14].Warminski and Kecik analyzed the motion of a nonlinear oscillator with attached pendulum, excited by the moment of its suspension point, the oscillator, and the pendulum being strongly coupled by inertial terms [15].In [16], Kecik and Warminski proposed a new suspension composed of a semiactive magnetorheological damper and a nonlinear spring in order to control motions.In this way, unstable areas and the chaotic or rotating motion of the pendulum are reduced.
A variation of the simple pendulum involving square plates was investigated by Rafat et al. [17].The equilibrium configurations and normal modes of oscillations are obtained.The equations of motion were solved numerically to produce Poincare sections.The accurate analytic solution of the nonlinear pendulum differential equation is obtained using homotopy analysis technique by Turkyilmazoglu [18].The obtained explicit analytical expressions for the frequency, period and displacement are compared with numerical ones.
Awrejcewicz [19] studied the mathematical pendulum motion oscillating in a plane rotating with angular velocity.The three-dimensional double pendulum, which is coupled by two universal joints, is investigated in [20].The multiple scales method was used in [21] for recognizing resonances occurring in a parametrically and externally excited nonlinear spring pendulum.Energy balance method was employed in [22] to obtain approximations for achieving the nonlinear frequency for pendulum attached to rolling wheels that is restrained by a spring.The nonlinear oscillations of pendulum wrapping on two cylindrical bases were investigated by Mazaheri et al. [23].To obtain an analytical solution, the multiple scale method is used and there are analyzed effects of amplitude and radius of cylinder.
Boubaker presented in [24] a survey on the inverted pendulum in nonlinear control theory offering an overall picture of historical, current and trend developments.Synchronization of two pendulums mounted on a mutual base is investigated by Alevras et al. [25] and the response of pendulum was obtained when the base was excited by a random sinusoidal force.The influence of an external harmonic excitation on a chain of nonlinear pendulum was explored by Jallouli et al. in [26] in case of simultaneous external and parametric excitations.
In this paper we propose a novel procedure, the Optimal Auxiliary Functions Method (OAFM), to investigate the nonlinear oscillations of a simple pendulum bounded by two cylinders at the point of suspension.The length of this pendulum varies due to wrapping around the cylinders.Such systems of pendulum with such additional conditions along with their dynamic behavior could find applications in aerospace engineering and shipping engineering.
Unlike other solution procedures applied to find approximate analytical solutions to nonlinear dynamical systems, the proposed approach is based upon original construction of the solution using a moderate number of convergence-control parameters, which are basic components of the original auxiliary functions introduced in the present developments.These parameters lead to a high precision, comparing our approximate solutions with exact or numerical ones.
The accuracy of the obtained results is proved by numerical developments, which validate the analytical results.

The Optimal Auxiliary Functions Method
The basics of OAFM can be found in [27,28], where OAFM is applied to solve different problems.In order to develop an application of the OAFM, let us consider the nonlinear differential equation [27][28][29]: where L is the linear operator, N is the nonlinear operator, and g is a known function, x being an independent variable and u(x) an unknown function at this stage.The initial or boundary conditions are: It is well-known that it is often very hard to find an exact solution for strongly nonlinear equations of type (1) and (2) [30].In order to find the approximate solution u(x), we suppose this can be expressed as where the initial and the first approximation will be obtained, as described below.After the substitution of Equation (3) into Equation (1), one obtains where C i , i = 1,2,...,s are the convergence-control parameters, which will be rigorously determined.
The initial approximation u 0 (x) may be determined from the linear equation while the first approximation is obtained from Equations ( 4) and ( 5): The nonlinear term from Equation ( 6) is expanded as In order to avoid the difficulties appearing in solving Equation ( 6) and also to accelerate the convergence of the solution u(x, C i ), instead of the last term, one can suggest another expression, so that this equation may be rewritten as where A 1 and A 2 are auxiliary functions which depend on initial approximation u 0 (x) and some convergence-control parameters C j , and C k , j = 1,2,...,p, k = p + 1, p + 2,...,s, and F(N[u 0 (x)]) are functions which depend on the expressions which appear within the nonlinear term N[u 0 (x)].It should be emphasized that the auxiliary functions A 1 and A 2 (namely optimal auxiliary functions) and F(N[u 0 (x)]) are not unique, but these auxiliary functions are of the same form, similar to u 0 (x).
More precisely, if u 0 (x) is a polynomial function, then A 1 and A 2 are sums of polynomial functions.If u 0 (x) is an exponential function, then A 1 and A 2 are sums of exponential functions.In the case of u 0 (x), which is a trigonometric function, it follows that A 1 and A 2 are sums of trigonometric functions, and so on.
In the case when N[u 0 (x)] = 0, then u 0 (x) is the exact solution of the original equation.
The initially unknown convergence-control parameters C j and C k may be rigorously and optimally determined via various methods, among them being the least square method, Galerkin method, collocation method, Ritz method, but the preferred one should be minimizing the square residual error: where in which the approximate solution u(x, C i ) is given by Equation (3).The unknown parameters C 1 , C 2 ,..., C s can be identified from the conditions Similar results could be obtained by imposing the conditions By using this above presented approach, the approximate solution is completed after the determination of the optimal values of convergence-control parameters C i , i = 1,2,...,s.Hence, our procedure involves the auxiliary functions A 1 and A 2 which provide an effective way to adjust and control the convergence of the final solutions u(x, C i ).It is necessary to remark the importance of carefully choosing the functions A 1 and A 2 involved in the construction of the first-order approximation u 1 (x, C i ).It was already proved that our method is easily applicable to solve nonlinear problems without small or large parameters, including systems with more degrees of freedom [27].

Equation of Motion
In what follows, we present the governing equation of the simple pendulum wrapping around two cylinders at the point of suspension [23].The length of the pendulum is L while the radius of cylinders is r (Figure 1).The motion of the system is described by the generalized coordinate θ, but the string length is changing.The kinetic energy can be expressed in the form where m is the mass of pendulum and the dot denotes differentiation with respect to time.
The potential energy becomes From the Lagrange's equation one can put After some manipulation, one obtains: ..
where a = r/L.The initial conditions for Equation (16) are

Application of OAFM to a Pendulum Wrapping on Two Cylinders
If one inserts the independent variable τ = Ωt and the dependent variable ϕ = θA −1 , then the governing Equations ( 16) and ( 17) become where Ω is the frequency of the system and prime denotes differentiation with respect to τ.
Having in view that cos τ|cos τ| = cos 2 τ(sgn(cos τ), sin(A cos τ) = A − A and substituting Equations ( 25)-( 28) into Equation (24), one can get Taking into account Equations ( 8) and ( 29), we can choose the auxiliary functions in the form where C 1 , C 2 , C 3 , and C 4 are unknown parameters and α, β, γ are obtained from Equation ( 29): We also may choose the auxiliary functions A 1 and A 2 , and the function F as follows or and so on.Substituting Equation (30) into Equation ( 8), the result is In order to avoid secular terms, the following condition should be imposed From Equations (31) and (35) one retrieves where The solution (34) is given by From Equations ( 3), ( 23) and (37), and from the transformations τ = Ωt and ϕ = θA −1 , one obtains the first-order approximate solution of Equation (16) as (cos Ωt − cos 7Ωt) where the coefficients α, β and γ are given in Equation (31) and Ω in Equation (36).

Results and Discussion
In order to emphasize the accuracy of our approach, we consider various sets of values for the parameters a, A, and L. We analyze the solution θ in 10 different cases and we develop comparisons between analytical and numerical integration results.Additionally, we represent a graphical comparison of the phase plane and a comparison between the frequencies Ω given by analytical developments (36) and numerical integration results, respectively.The calculation parameters were chosen as to reflect real cases, which could be encountered in practice.

Case 1
First, we consider A = 0.1, a = 0.2, L = 0.6, and g = 9.8.Using the proposed procedure, by minimizing the residual function, the optimal values of the convergence-control parameters C i and the frequency (36) are The solution given by (38) can be written as follows: In Figures 2 and 3 is plotted the comparison between approximate solution (39) and numerical integration results, and the phase plane in this case, respectively.
First, we consider A = 0.1, a = 0.2, L = 0.6, and g = 9.8.Using the proposed procedure, by minimizing the residual function, the optimal values of the convergence-control parameters Ci and the frequency (36) are The solution given by (38) can be written as follows: In Figures 2 and 3 is plotted the comparison between approximate solution (39) and numerical integration results, and the phase plane in this case, respectively.Analyzing the comparison between the approximate and numerical integration results ented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means ur procedure are almost identical with the results obtained using a numerical integration oach.Moreover, from Table 1 one can be observe that the accuracy of the approximate uency is remarkably good when compared to numerical results.From Figures 2-19, it can be seen that the errors of the approximate solutions increase with ect to increasing values of the parameters a and A. Additionally, for the particular case of ical simple pendulum, the results obtained through our procedure are in very good agreement numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that frequency of the system is increased by increasing the radius of cylinders (parameter a).itionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is ased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and r 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping   Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution.In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping (50 In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results. From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case o classical simple pendulum, the results obtained through our procedure are in very good agreemen with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce tha the frequency of the system is increased by increasing the radius of cylinders (parameter a) Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution (39).

Case 2
For A = 0.1, a = 0.4, L = 0.6 we obtain The solution given by (38) in this case can be written as follows: The comparison between analytical solution (40) and numerical integration results is presented in Figures 4 and 5.
The solution given by (38) in this case can be written as follows: The comparison between analytical solution (40) and numerical integration results is presented in Figures 4 and 5.    Analyzing the comparison between the approximate and numerical integration results sented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means ur procedure are almost identical with the results obtained using a numerical integration roach.Moreover, from Table 1 one can be observe that the accuracy of the approximate uency is remarkably good when compared to numerical results.From Figures 2-19, it can be seen that the errors of the approximate solutions increase with ect to increasing values of the parameters a and A. Additionally, for the particular case of sical simple pendulum, the results obtained through our procedure are in very good agreement h numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that frequency of the system is increased by increasing the radius of cylinders (parameter a).itionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is eased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping numerical A = 0.4, L = 0.6._____ numerical_ _ _ _approximate solution.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution.In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping   Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution (40).Graphical comparisons between analytical and numerical results in this case are presented in Figures 6 and 7.

Case 3
For A = 0.1, a = 0.6, L = 0.  (41) Graphical comparisons between analytical and numerical results in this case are presented in Figures 6 and 7.   (50) In Figure 20, we compared the results obtained through OAFM with numerical integration lts for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a parison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results ented in Figures  In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with approximate solution.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping  Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution (41).

Case 4
For A = 0.2, a = 0.2, L = 0.6, it holds that      (50) In Figure 20, we compared the results obtained through OAFM with numerical integration lts for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a parison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results ented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means ur procedure are almost identical with the results obtained using a numerical integration oach.Moreover, from Table 1 one can be observe that the accuracy of the approximate uency is remarkably good when compared to numerical results.From Figures 2-19, it can be seen that the errors of the approximate solutions increase with ect to increasing values of the parameters a and A. Additionally, for the particular case of ical simple pendulum, the results obtained through our procedure are in very good agreement numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that frequency of the system is increased by increasing the radius of cylinders (parameter a).itionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is ased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and r 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping   (50) In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21  (50) In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with approximate solution (42).

Case 5
For A = 0.2, a = 0.4, L = 0.6, it holds that A comparison between the analytical solution (43) and corresponding numerical integration results is presented in Figures 10 and 11.

Case 5
For A = 0.2, a = 0.4, L = 0.6, it holds that      (50) In Figure 20, we compared the results obtained through OAFM with numerical integration lts for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a parison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results ented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means ur procedure are almost identical with the results obtained using a numerical integration oach.Moreover, from Table 1 one can be observe that the accuracy of the approximate uency is remarkably good when compared to numerical results.From Figures 2-19, it can be seen that the errors of the approximate solutions increase with ect to increasing values of the parameters a and A. Additionally, for the particular case of ical simple pendulum, the results obtained through our procedure are in very good agreement numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that frequency of the system is increased by increasing the radius of cylinders (parameter a).itionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is ased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and r 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping numerical comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution.In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration result presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by mean of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximat frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case o classical simple pendulum, the results obtained through our procedure are in very good agreemen with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce tha the frequency of the system is increased by increasing the radius of cylinders (parameter a) Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system i increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution (43).

Case 6
For A = 0.2, a = 0.6, L = 0.6, we obtain In Figures 12 and 13 is plotted the comparison between approximate solution (44) and numerical integration results in this case.In Figures 12 and 13 is plotted the comparison between approximate solution (44) and numerical integration results in this case.Analyzing the comparison between the approximate and numerical integration results sented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means ur procedure are almost identical with the results obtained using a numerical integration roach.Moreover, from Table 1 one can be observe that the accuracy of the approximate uency is remarkably good when compared to numerical results.From Figures 2-19, it can be seen that the errors of the approximate solutions increase with ect to increasing values of the parameters a and A. Additionally, for the particular case of sical simple pendulum, the results obtained through our procedure are in very good agreement h numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that frequency of the system is increased by increasing the radius of cylinders (parameter a).itionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is eased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping  Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution.In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping   Analyzing the comparison between the approximate and numerical integration result presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by mean of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximat frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case o classical simple pendulum, the results obtained through our procedure are in very good agreemen with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce tha the frequency of the system is increased by increasing the radius of cylinders (parameter a) Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system i increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution (44).

Case 7
In this case, for A = 0.3, a = 0.2, L = 0.6, yields Graphical comparisons between analytical and numerical results are presented for this case in Figures 14 and 15.
Graphical comparisons between analytical and numerical results are presented for this case in Figures 14 and 15    (50) In Figure 20, we compared the results obtained through OAFM with numerical integration lts for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a parison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results sented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping  Analyzing the comparison between the approximate and numerical integration result presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by mean of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case o classical simple pendulum, the results obtained through our procedure are in very good agreemen with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce tha the frequency of the system is increased by increasing the radius of cylinders (parameter a) Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system i increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution (45).

Case 8
Considering A = 0.3, a = 0.4, L = 0.6, it follows that    (50) In Figure 20, we compared the results obtained through OAFM with numerical integration lts for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a parison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results sented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means ur procedure are almost identical with the results obtained using a numerical integration roach.Moreover, from Table 1 one can be observe that the accuracy of the approximate uency is remarkably good when compared to numerical results.From Figures 2-19, it can be seen that the errors of the approximate solutions increase with ect to increasing values of the parameters a and A. Additionally, for the particular case of sical simple pendulum, the results obtained through our procedure are in very good agreement h numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that frequency of the system is increased by increasing the radius of cylinders (parameter a).itionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is eased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping   In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration (50 In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration result presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by mean of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results. From Figures 2-19, it can be seen that the errors of the approximate solutions increase with approximate solution (46).

Case 9
In this case, we consider A = 0.3, a = 0.6, L = 0.6, such that Graphical comparisons between analytical and numerical results in this case are presented in Figures 18 and 19.Moreover Table 1 presents a comparison between the values of the frequency obtained in the above considered cases.Graphical comparisons between analytical and numerical results in this case are presented in Figures 18 and 19.Moreover Table 1 presents a comparison between the values of the frequency obtained in the above considered cases.(50) In Figure 20, we compared the results obtained through OAFM with numerical integration lts for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a parison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results sented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means ur procedure are almost identical with the results obtained using a numerical integration roach.Moreover, from Table 1 one can be observe that the accuracy of the approximate uency is remarkably good when compared to numerical results.From Figures 2-19, it can be seen that the errors of the approximate solutions increase with ect to increasing values of the parameters a and A. Additionally, for the particular case of sical simple pendulum, the results obtained through our procedure are in very good agreement h numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that frequency of the system is increased by increasing the radius of cylinders (parameter a).itionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is eased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping numerical In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution.
In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).
Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration result presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by mean of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximat frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case o classical simple pendulum, the results obtained through our procedure are in very good agreemen with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce tha the frequency of the system is increased by increasing the radius of cylinders (parameter a) Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system i increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping approximate solution (47).

Case 10
The classical simple pendulum is obtained from Equation ( 16) in the case when no cylinder exists.Therefore, for a = r/L = 0 we obtain from (36) the approximate frequency The approximate solution for the simple pendulum is obtained from Equation (38) with the following coefficients given by Equation (31) for this particular case: The optimal values of the control parameters and the approximate frequency in this case are, respectively In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.(50) ure 20, we compared the results obtained through OAFM with numerical integration the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a between the phase planes in this case is presented.zing the comparison between the approximate and numerical integration results n Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means cedure are almost identical with the results obtained using a numerical integration Moreover, from Table 1 one can be observe that the accuracy of the approximate s remarkably good when compared to numerical results.Figures 2-19, it can be seen that the errors of the approximate solutions increase with increasing values of the parameters a and A. Additionally, for the particular case of ple pendulum, the results obtained through our procedure are in very good agreement rical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that cy of the system is increased by increasing the radius of cylinders (parameter a).ly, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approach.Moreover, from Table 1 one can be observe that the accuracy of the approximate frequency is remarkably good when compared to numerical results.
From Figures 2-19, it can be seen that the errors of the approximate solutions increase with respect to increasing values of the parameters a and A. Additionally, for the particular case of classical simple pendulum, the results obtained through our procedure are in very good agreement with numerical integration results.From the cases 5.4-5.6, and 5.7-5.8respectively, we deduce that the frequency of the system is increased by increasing the radius of cylinders (parameter a).Additionally, from the cases 5.1, 5.4 and 5.7 it can be seen that the frequency of the system is increased by increasing the amplitude A. The same conclusion is obtained from the cases 5.2, 5.5 and 5.8 or 5.3, 5.6 and 5.9, respectively.The sources of nonlinear oscillations of the pendulum wrapping   In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.In Figure 20, we compared the results obtained through OAFM with numerical integration results for the particular care of simple pendulum for A = 0.4, L = 0.6, while in Figure 21 a comparison between the phase planes in this case is presented.Analyzing the comparison between the approximate and numerical integration results presented in Figures 2-21 for the cases 5.1-5.10, it can be observed that the results obtained by means of our procedure are almost identical with the results obtained using a numerical integration approximate solution (50).

Figure 1 .
Figure 1.Simple pendulum wrapping around the cylinders.Figure 1. Simple pendulum wrapping around the cylinders.

Figure 1 .
Figure 1.Simple pendulum wrapping around the cylinders.Figure 1. Simple pendulum wrapping around the cylinders.

32 5
Figures 8 and 9 emphasize the comparison of the analytical solution (42) with numerical integration results.

Figures 8
Figures 8 and 9 emphasize the comparison of the analytical solution (42) with numerical integration results.
approximate solution.

32 5
Figures 16 and 17 emphasize the comparison of the analytical solution (46) with numerical integration results.
for the cases 5.1-5.10, it can be observed that the results obtained by means

Table 1 .
Comparison between the numerical solution of the frequency and the approximate frequency (36).