Conservative systems represent a large number of naturally occurring and artificially designed scientific and engineering systems. A key consideration in the theory and application of nonlinear conservative systems is the solution of the governing nonlinear ordinary differential equation. This chapter surveys the recent approximate analytical schemes for the periodic solution of nonlinear conservative systems and presents a recently proposed approximate analytical algorithm called continuous piecewise linearization method (CPLM). The advantage of the CPLM over other analytical schemes is that it combines simplicity and accuracy for strong nonlinear and large-amplitude oscillations irrespective of the complexity of the nonlinear restoring force. Hence, CPLM solutions for typical nonlinear Hamiltonian systems are presented and discussed. Also, the CPLM solution for an example of a non-Hamiltonian conservative oscillator was presented. The chapter is aimed at showcasing the potential and benefits of the CPLM as a reliable and easily implementable scheme for the periodic solution of conservative systems.
- Hamiltonian system
- conservative system
- nonlinear vibration
- continuous piecewise linearization method
- periodic solution
- nonnatural system
- perturbation method
1.1 Hamiltonian and non-Hamiltonian conservative systems
Conservative systems can be defined as oscillating or vibrating systems in which the total energy content of the system remains constant. In order words, the total energy in the system is conserved. Ideally, such a system will continue to be in periodic oscillatory motion ad infinitum because the energy content of the system does not diminish due to the absence of dissipative force or increase due to additional energy input. However, for real cases where dissipative mechanisms such as friction or viscous damping cannot be completely eliminated, a conservative system can be thought of as one in which the energy dissipated is negligible during the time range under consideration. For example, the first few seconds of the oscillation of a simple pendulum may be considered conservative since the effect of air friction is negligible, but in the long run, the initial energy content is gradually dissipated until the pendulum comes to a halt. Other examples of practical conservative systems include mass-spring oscillator, structural elements (i.e., beams, plates, and shells), slider-crank mechanism , human eardrum , relativistic oscillator , planetary orbits around the sun , and current-carrying conductor in the electric field of an infinite rod . Hence, a large number of oscillating physical systems can be studied as conservative systems.
At any point in time, the energy of a conservative system is composed of kinetic () and potential () energies except at critical points where the total energy may be only kinetic () or potential (). Generally, it is expected that and , where is the generalized displacement. Naturally, and are not expected to form a product in the function , but in some cases, they do Therefore, two types of conservative systems are distinguished namely: natural and nonnatural conservative systems. The natural conservative systems are those in which the kinetic energy can be expressed as a pure quadratic function of velocity, i.e., does not contain a product of the velocity and displacement. They are also known as
1.2 Recent advances in solution schemes for nonlinear conservative oscillators
Exact analytical solutions for the nonlinear vibration models of conservative systems can be derived only in very few situations, and the solutions are usually derived in terms of special functions. Alternatively, highly accurate numerical solutions can be obtained for the nonlinear vibration model of any conservative system. However, as it is well recognized among the nonlinear science community, numerical solutions often have the limitations of lack of physical insight and convergence issues. Furthermore, there is the possibility of obtaining inaccurate convergent solutions for a nonlinear ordinary differential equation (ODE) , thus necessitating the independent verification of the convergent numerical solution by another numerical or analytical method. These limitations have driven the search for approximate analytical schemes capable of providing periodic solutions to nonlinear conservative oscillators. It can be rightly concluded that this search has been very fruitful considering the many approximate analytical schemes that now appear in the nonlinear science literature. The purpose of this section is to provide a brief survey of some of the notable achievements in the development of approximate analytical schemes for the periodic solution of nonlinear conservative oscillators. It should be noted that an approximate analytic method for nonlinear oscillators is considered adequate if it gives accurate predictions for the frequency-amplitude response and the oscillation history as well [7, 8].
Approximate analytical techniques to solve the nonlinear ODE governing the oscillations of a conservative system have been formulated for at least 100 years and can be classified as perturbation and nonperturbation methods. The first attempts were based on perturbation theory and are referred to as
More recently, in the last four decades, a number of approximate analytical schemes have been proposed. Most of these recent schemes are nonperturbation methods, but some recent perturbation methods that attempt to improve on their classical counterparts have been formulated too. The recent perturbation methods include -method , Homotopy perturbation method  and its variants [12, 13, 14, 15, 16, 17], modified Lindstedt-Poincare method [18, 19, 20, 21], book-keeping parameter method , iteration perturbation method , parameterized perturbation method , perturbation incremental method , and linearized perturbation method . A review article on some of the recent perturbation methods has been published by He . The main point of the recent perturbation methods is to deal with the issue of the small parameter in order to formulate solutions that are applicable to small- and large-amplitude oscillations and also weak and strong nonlinear oscillations. Although the higher order approximations of the recent perturbation methods have been very successful in producing accurate estimates of the frequency-amplitude response, the same cannot be said of their estimation of the oscillation history. Studies [7, 28] have shown that the higher order approximations of the recent perturbation methods produce large unbounded errors in the oscillation history during large-amplitude oscillations and are, therefore, not better than the classical perturbation methods in this regard. A plausible explanation for this observation is that it occurs because perturbation methods are based on asymptotic series that are inherently divergent for amplitudes greater than unity . Therefore, it may not be possible to formulate perturbation schemes that would correctly predict the oscillation history of large-amplitude vibrations.
In contrast to the perturbation methods, the nonperturbation methods do not use any small or artificial parameter. Examples include Adomian decomposition method , Homotopy analysis method , Variational iteration method , Energy balance method  and its modifications [32, 33, 34], He Chengtian’s interpolation method  also called max-min approach, amplitude-frequency formulation , Hamiltonian approach , global error minimization method , Harmonic balance method  and its modifications [38, 39, 40, 41, 42], cubication methods [43, 44, 45, 46], variational methods [47, 48, 49], differential transform method , and continuous piecewise linearization method . Nonperturbation methods also have various limitations. For instance, a study  showed that the Adomian decomposition method does not converge to the correct solution in some cases, and the study proposed an optimal convergence acceleration parameter to deal with this issue. Also, methods that rely on a simple harmonic approximation of the oscillation history, such as the energy balance method, amplitude-frequency formulation, Hamiltonian approach, max-min approach, and variational methods, can only give reliable estimate of the frequency-amplitude response. Sometimes these methods perform poorly in predicting the oscillation history during large-amplitude and/or strong nonlinear vibrations. Other methods that usually require high-order approximations, such as Adomian decomposition method, harmonic balance method, and variational iteration method, present algebraic complexities in their determination of higher order solutions and may be impractical for oscillators with highly complex nonlinearities such as the slider-crank mechanism  and the bifilar pendulum . Furthermore, it has been observed that higher order estimates do not always improve the solution of the oscillation history . Finally, some nonperturbation methods are heuristic in nature (e.g., energy balance method and variational methods) and require experience to choose the initial trial function and the condition for error minimization .
The continuous piecewise linearization method (CPLM) is an iterative analytic algorithm that was formulated to overcome most of the above challenges by providing simple and accurate solutions for the oscillation history and frequency-amplitude response of Duffing-type oscillators. In another study , the CPLM was modified in order to generalize it so that it can handle more complicated nonlinear conservative oscillators. Interestingly, the CPLM does not require higher order approximations or any small, artificial, or embedded parameter. Also, the algorithm is inherently stable, straightforward, and based on closed-form analytical approximations. This chapter is aimed at presenting the generalized CPLM algorithm as a veritable approach for accurate periodic solution of Hamiltonian and non-Hamiltonian conservative oscillators with complex nonlinearity. As is shown later, the CPLM retains the same simplicity in its implementation irrespective of the complexity of the nonlinear conservative system.
2. Continuous piecewise linearization method
2.1 Concept of the continuous piecewise linearization method
The main idea of the CPLM is based on the piecewise linearization of the nonlinear restoring force of a conservative oscillator. The linearization technique used by the CPLM was first applied in another algorithm [6, 54] for the solution of half-space impact models called force indentation linearization method (FILM). The FILM has been applied to formulate theoretical solutions for rigid body motions and local compliance response during nonlinear elastoplastic impact of dissimilar spheres . However, because the FILM is limited to impact excitations that are nonoscillatory, it cannot be applied to nonlinear conservative oscillators. Hence, the CPLM applies the piecewise linearization technique of the FILM to provide a periodic solution for nonlinear conservative oscillators.
Essentially, the linearization technique of the CPLM involves equal discretization of the nonlinear restoring force with respect to displacement (Figure 1) and formulating a linear restoring force for each discretization. Therefore, a linear ODE can be derived for each discretization. The solution of the linear ODE approximates the solution of the original nonlinear oscillator for a time-range that is automatically determined by the CPLM and updated continuously from one discretization to the next. Details on the discretization and linearization technique of the CPLM can be found in the following references [6, 8, 54], while the applications of the CPLM to nonlinear conservative systems are presented in [8, 56].
2.2 Mathematical formulation of the continuous piecewise linearization method
The standard form for representation of a nonlinear conservative oscillator moving in the -direction is given as:
where is the nonlinear restoring force as shown in Figure 1. In Figure 1, the numbering on the horizontal axis represents the boundary points of the discretization. The discretization represents a general discretization with start point at and endpoint at . Figure 1 shows that for each discretization, the slope of the linear approximation of the restoring force can either be positive or negative. To account for this possibility, the linearized force for the discretization can be expressed as:
where is the linear slope of and .
Eq. (10) is a nonhomogeneous linear ODE and its solution depends on whether the sign is positive or negative.
2.2.1 Solution for positive linearized stiffness
If , the solution for the displacement and velocity can be expressed as:
where , , and . The initial conditions and other parameters are determined based on the oscillation stage. For the oscillation stage that moves from to , the initial conditions for each discretization are and , where , and the other parameters are calculated as:
For the oscillation stage that moves from to , the initial conditions are and ; the other parameters are calculated as:
The time at the end of each discretization is , and the end conditions and are calculated by replacing with in the formulae for the initial conditions.
2.2.2 Solution for negative linearized stiffness
If , the solution for the displacement and velocity can be expressed as follows:
where ; . Applying the initial conditions to Eqs. (7a) and (7b) gives: ; . The initial and end conditions are determined in the same way as for above. The end conditions are applied in Eq. (7a) to get the time interval for each discretization as:
The sign before the square root in Eq. (8) is negative for the oscillation stage that moves from to and vice versa. We note that if , then and
2.2.3 Solution for zero linearized stiffness
In very rare situations, we may have for one or two discretization around the turning points or relatively flat regions of the restoring force. This is likely when is very small, i.e., for very large , and can be eliminated by increasing or decreasing slightly. However, if we want to account for , then we get :
where and . Hence, the time interval is derived from Eq. (11) as:
2.3 Remarks on the CPLM algorithm
From the above presentation of the CPLM formulation, it is obvious that the CPLM algorithm is simple and can be implemented by undergrads without difficulty.
The CPLM is inherently stable and does not have convergence issues .
For few discretization, say , the CPLM algorithm can be implemented with reasonable accuracy using a pocket calculator. However, the CPLM is better executed using a simple code in any programming language such as MATLAB and Mathematica or using a customized MS Excel spreadsheet.
When dealing with conservative oscillators with odd nonlinearity, which are symmetrical about the origin, discretization of the restoring force is only required for 0 to . This means that there will be discretizations from to .
The CPLM algorithm retains the same simplicity in implementation irrespective of the complexity of the restoring force. Only the constant and the integral of the restoring force are to be evaluated anew for any oscillator.
The CPLM relies on the explicit expression of restoring force, which means that the model for the oscillator must be expressed in the form of Eq. (1). For Hamiltonian systems, the oscillator model is formulated naturally in the form of Eq. (1). For non-Hamiltonian conservative systems, the oscillator model is not formulated naturally in the form of Eq. (1). Therefore, there is a need to transform the model of non-Hamiltonian conservative systems into the form of Eq. (1) before applying the CPLM algorithm. Fortunately, this transformation only requires a simple algebraic manipulation as demonstrated in Section 4.
The phase equation gives the relationship between the state variables (displacement and velocity) and can be derived exactly in closed-form for all conservative systems. The CPLM extends this bilateral relationship into a tripartite one by finding the value of the corresponding independent variable (i.e., time) that matches the values of the state variables in each discretization.
For few discretization, say , it would be necessary to extract values within each discretization in order to obtain a smooth plot of the oscillation history. The values can be extracted using the approximate solution of the displacement. However, for many discretizations, say , there is no need to extract values from any discretization.
The usual initial conditions investigated for nonlinear conservative oscillators are nonzero displacement and zero velocity. However, the CPLM algorithm can comfortably handle nonzero initial conditions for displacement and velocity.
3. Periodic solution of typical Hamiltonian systems
3.1 Nonlinear simple pendulum
The simple pendulum is arguably the most investigated physical system and provides very interesting insights into nonlinear phenomena. Butikov  calls it “an antique but evergreen physical model.” The undamped oscillation of a simple pendulum is a Hamiltonian system governed by the well-known nonlinear ODE as shown:
where is the angular displacement, , is the length of the pendulum, and . The initial conditions are given as and . The same initial conditions are applicable to all other oscillators discussed subsequently. The exact solution to Eq. (13) is expressed in terms of elliptic functions. The displacement and natural frequency are given as :
Where is the Jacobi elliptic sine function, , and is the complete elliptic integral of the first kind given as:
From Eq. (13), the restoring force for the pendulum is and looks like the plot in Figure 1. This means that and . The initial and final velocities for each discretization are determined using the expression for , and all the other constants for the solution of each discretization are determined based on . A plot of the frequency-amplitude response for the simple pendulum when is given in Figure 2a and the corresponding error of the CPLM solution in comparison with the exact solution (Eq. (15)) is shown in Figure 2b. We see that for , the maximum error in the CPLM estimate is less than 0.45% for and 0.14% for . Also, the oscillation history for moderate-amplitude () and large-amplitude () is shown in Figure 3 and an excellent agreement between the CPLM solution for and the exact solution (Eq. (14)) is observed. We noted that trigonometric nonlinearity is usually difficult to deal with and that is why the CPLM shows a relatively slow convergence to accurate results. Hence, many discretizations (e.g., ) are required to get an accuracy that is within 1.0% of the exact solution during large-amplitude oscillations (90) of the simple pendulum.
3.2 Motion of satellite equidistant from twin stars
Consider the motion of a satellite along a path that is equidistant from two identical massive stars with mutually interacting gravitational fields. If the distance between the two stars is and the coordinate of the satellite motion is , then the equation of motion of the satellite is given as :
where is the mass of a star and the restoring force is . Eq. (17) shows that the satellite-star interaction results in a conservative oscillation of the satellite. Figure 4 shows the nonlinear restoring force, which is an irrational force because of the bottom square root. The restoring force spikes on both sides of the vertical axis close to the origin. The spikes indicate the point when the satellite is most influenced by the mutual gravitational field of the stars. Away from the origin, the restoring force decreases gradually and approaches the horizontal axis asymptotically. This means that the satellite is far away from the stars and experiences a much smaller gravitational force. This problem was discussed qualitatively by Jordan and Smith  and Arnold , but here, the periodic solution was investigated.
The main CPLM constant is and the velocity was derived as: . The periodic solutions obtained by the CPLM and exact numerical solution are shown in Figures 5 and 6. The numerical solution was obtained by solving Eq. (17) using the NDSolve function in Mathematica™. The NDSolve function is a Mathematica subroutine for solving ordinary, partial, and algebraic differential equations numerically. In its basic form, it automatically selects the numerical method to use from a list of standard methods such as implicit Runge–Kutta, explicit Runge–Kutta, symplectic partitioned Runge–Kutta, predictor–corrector Adams, and backward difference methods. In some cases, the NDSolve function can combine two or more methods to obtain the required solution. This basic form is preferable because the NDSolve function uses the method(s) that best solves the differential equation considering accuracy and solution time. Hence, the NDSolve function was used in its basic form for all numerical solutions obtained in this chapter.
The input values used for simulation are and . In contrast to the simple pendulum, the oscillation of the satellite requires less discretization for accurate results because there is no trigonometric nonlinearity. The maximum error of the CPLM estimate for the frequency-amplitude response is less than 0.55% for and 0.20% for . Significantly higher accuracies can be achieved by increasing , but the results show that gives reasonably accurate estimates.
On the other hand, Figure 6 shows the oscillation history of the satellite during small-amplitude () and large-amplitude () oscillations. The former gives a simple harmonic response with a natural frequency that is independent of the amplitude and approximately equal to , while the latter exhibits an
3.3 Mass-spring oscillator with fractional nonlinearity
An interesting oscillator that has been the subject of several studies [60, 61, 62, 63, 64, 65] is the Hamiltonian oscillator with odd fractional nonlinearity. For the purpose of the present investigation, we consider an oscillator that is characterized by a general fractional nonlinearity as follows :
where the restoring force, , has a fractional index for all . The main CPLM constant is evaluated as and the velocity as . The periodic solution for the case of , i.e., oscillator, is shown in Figures 7 and 8. The exact frequency-amplitude response used for the verification of the CPLM solution is :
while the exact oscillation history was obtained by the numerical solution of Eq. (18) using the NDSolve function in Mathematica™. The CPLM solution demonstrates an excellent agreement with the exact solution.
Figure 7 compares the CPLM estimates of the frequency-amplitude response with Eq. (19), and the maximum error of the CPLM solution is 0.22% for and 0.076% for . The error in the CPLM estimate is approximately constant for all amplitudes, and the maximum error is well below 1.0% for . The results also reveal that the frequency approaches zero as . In Figure 8, the oscillation history for small-amplitude () and large-amplitude () oscillations is shown to exhibit similar anharmonic response, which is an indication of strong nonlinearity. Therefore, it can be concluded that the oscillator is highly nonlinear for all amplitudes. This quality of possessing strong nonlinearity for all amplitudes is in contrast to most Hamiltonian oscillators that are linear for small amplitudes, e.g., the oscillators considered in Sections 3.1 and 3.2 above. Another Hamiltonian oscillator that possesses strong nonlinearity for all amplitudes is the geometrically nonlinear crank .
4. Periodic solution of non-Hamiltonian conservative systems
The non-Hamiltonian conservative systems are generally more complex and difficult to solve compared with the Hamiltonian systems. In order to demonstrate the application of the CPLM algorithm to non-Hamiltonian conservative systems, we consider the motion of a particle on a rotating parabola. This system consists of a frictionless mass sliding on a vertical parabolic wire described by for and rotating at a constant speed, , about the -axis (Figure 9a). The -axis represents the perpendicular displacement of the mass from the -axis. The kinetic and potential energies of the system are given as :
Hence, the Lagrangian is:
Next, we substitute Eq. (21) into the Euler–Lagrange equation to derive the equation of motion. The Euler–Lagrange equation can be written as:
where and the initial conditions are: and .
where is a constant representing the total energy in the system. Eq. (24) confirms that the Hamiltonian, is not constant. Applying the initial conditions, we get so that . Substituting this expression for in Eq. (24) and after algebraic simplification, we get:
Therefore, the restoring force is . Figure 9b shows that is linear at small amplitudes and strongly nonlinear at large amplitudes. The main CPLM constant was calculated as , and the velocity was evaluated as . The exact time period for this oscillator can be derived in terms of elliptic function as follows :
where is the complete elliptic integral of the second kind and . Then, the exact frequency was computed as , while the exact oscillation history was obtained by solving Eq. (25) numerically using the NDSolve function in Mathematica™.
A comparison of CPLM frequency estimate and the exact frequency is shown in Figure 10, while the oscillation histories for and are shown in Figure 11. As demonstrated in , periodic solutions for this system exist only for . Hence, the simulations in Figures 10 and 11 were conducted for and . An excellent agreement is observed between the CPLM estimates and the exact results. For , the maximum error in the CPLM estimate of the frequency-amplitude response is 0.642% for and 0.101% for , both of which are well below 1.0%. Also, the CPLM solution gives an accurate prediction of the strong anharmonic response in the oscillation history as shown in Figure 11.
5. Concluding remarks
Conservative oscillators generally exhibit nonlinear response, and they form a large class of natural and artificially vibrating systems. Hence, the study of the dynamic response of nonlinear conservative systems is important for understanding many physical phenomena and the design of systems. The main challenge in the theoretical analysis of nonlinear conservative systems is that exact solutions are normally not available except for a few special cases where exact solutions are derived in terms of special functions.
To date, many approximate analytical methods have been formulated for the periodic solution of nonlinear conservative oscillators. This chapter provides a brief survey of the recent advances in the formulation of approximate analytical schemes and then introduced a recent approximate analytical algorithm called the continuous piecewise linearization method. The CPLM has been shown to overcome the challenges of solution accuracy and simplicity usually encountered in using most of the existing approximate analytical methods. The CPLM combines major desirable features of solution schemes such as inherent stability, accuracy, and simplicity. It is simple enough to be introduced at the undergraduate level and is capable of handling conservative oscillators with very complex nonlinearity. Conservative systems of broad interest were used to demonstrate the wide applicability of the CPLM algorithm. As demonstrated above, an accuracy of less than 1.0% relative error can be achieved for most oscillators using few discretizations, say , except for oscillators with trigonometric nonlinearity where such accuracy is achieved with many discretizations. This chapter has been designed to stimulate interest in the use of CPLM for analyzing various types of conservative systems, especially those with complex nonlinearity.
The authors are grateful to the Vice-Chancellor of the University of Port Harcourt, Prof. Ndowa E.S. Lale and the Dean of Engineering, Prof. Ogbonna Joel, for approving the publication grant for this chapter. The OAPF was funded by the African Centre of Excellence in Oilfield Chemicals Research (ACE-FOR), University of Port Harcourt.
Fidlin A. Nonlinear Oscillations in Mechanical Engineering. New York: Springer; 2005
He JH. Preliminary report on the energy balance for nonlinear oscillations. Mechanics Research Communications. 2002; 29:107-111
Jordan DW, Smith P. Nonlinear Ordinary Differential Equations: Problems and Solutions. Oxford: Oxford University Press; 2007
Nayfeh AH, Mook DT. Nonlinear Oscillations. New York: John Wiley & Sons; 1995
Shaban M, Ganji DD, Alipour MM. Nonlinear fluctuation, frequency and stability analyses in free vibration of circular sector oscillation systems. Current Applied Physics. 2010; 10:1267-1285
Big-Alabo A, Cartmell MP, Harrison P. On the solution of asymptotic impact problems with significant localised indentation. Journal of Mechanical Engineering Sciences. 2017; 231(5):807-822
Sanchez NE. A view to the new perturbation technique valid for large parameters. Journal of Sound and Vibration. 2005; 282:1309-1316
Big-Alabo A. Periodic solutions of Duffing-type oscillators using continuous piecewise linearization method. Mechanical Engineering Research. 2018; 8(1):41-52
Nayfeh AH. Perturbation Methods. New York: John Wiley & Sons; 1973
Bender CM, Milton KA, Pinsky SS, Simmons LM Jr. A new perturbative approach to nonlinear problems. Journal of Mathematical Physics. 1989; 30(7):1447-1455
He JH. Homotopy perturbation technique. Computer Methods in Applied Mechanics and Engineering. 1999; 178(3/4):257-262
Odibat ZM. A new modification of the homotopy perturbation method for linear and nonlinear operators. Applied Mathematics and Computation. 2007; 189(1):746-753
Khan Y, Wu Q. Homotopy perturbation transform method for nonlinear equations using He’s polynomials. Computer and Mathematics with Applications. 2019; 61(8):1963-1967
El-Dib YO, Moatimid GM. Stability configuration of a rocking rigid rod over a circular surface using the homotopy perturbation method and Laplace transform. Arabian Journal for Science and Engineering. 2019; 44(7):6581-6591
Kuang W, Wang J, Huang C, Lu L, Gao D, Wang Z, et al. Homotopy perturbation method with an auxiliary term for optimal design of a tangent nonlinear packaging system. Journal of Low Frequency, Noise, Vibration and Active Control. 2019. DOI: 10.1177/1461348418821204
Adamu MY, Ogenyi P, Tahir AG. Analytical solutions of nonlinear oscillator with coordinate-dependent mass and Euler-Lagrange equation using the parameterized homotopy perturbation method. Journal of Low Frequency, Noise, Vibration and Active Control. 2019. DOI: 10.1177/1461348418821222
Song H. A modification of homotopy perturbation method for a hyperbolic tangent oscillator arising in nonlinear packaging system. Journal of Low Frequency, Noise, Vibration and Active Control. 2019. DOI: 10.1177/1461348418822135
Cheung YK, Chen SH, Lau SL. A modified Lindstedt-Poincare method for certain strongly non-linear oscillators. International Journal of Non-Linear Mechanics. 1991; 26(3/4):367-378
He JH. Modified Lindstedt-Poincare methods for some strongly nonlinear oscillations part III: Double series expansion. International Journal of Nonlinear Science and Numerical Simulation. 2001; 2:317-320
He JH. Modified Lindstedt-Poincare methods for some strongly nonlinear oscillations part I: Expansion of a constant. International Journal of Non-Linear Mechanics. 2002; 37:309-314
He JH. Modified Lindstedt-Poincare methods for some strongly nonlinear oscillations part II: A new transformation. International Journal of Non-Linear Mechanics. 2002; 37:315-320
He JH. Bookkeeping parameter in perturbation methods. International Journal of Nonlinear Science and Numerical Simulation. 2001; 2:257-264
He JH. Iteration perturbation method for strongly nonlinear oscillations. Journal of Vibration and Control. 2001; 7(5):631-642
He JH. Some new approaches to Duffing equation with strongly and high order nonlinearity (II) parameterized perturbation technique. Communications in Nonlinear Science and Numerical Simulation. 1999; 4(1):81-83
Chan HSY, Chung KW, Xu Z. A perturbation incremental method for strongly nonlinear oscillators. International Journal of Non-Linear Mechanics. 1996; 31(1):59-72
He JH. Linearized perturbation technique and its applications to strongly nonlinear oscillators. Computers & Mathematics with Applications. 2003; 45(1-3):1-8
He JH. Some asymptotic methods for strongly nonlinear equations. International Journal of Modern Physics B. 2006; 20(10):1141-1199
He JH. Author’s reply to ‘A view to the new perturbation technique valid for large parameters’ by N.E. Sanchez. Journal of Sound and Vibration. 2005; 282:1317-1320
Adomian GA. Review of the decomposition method in applied mathematics. Journal of Mathematical Analysis and Applications. 1988; 135:501-544
Laio SK. On the homotopy analysis method for nonlinear problems. Applied Mathematics and Computation. 1994; 147(2):499-513
He JH. Variational iteration method: A kind of nonlinear analytical technique: Some examples. International Journal of Nonlinear Mechanics. 1999; 34(4):699-708
Durmaz S, Kaya MO. High-order energy balance method to nonlinear oscillators. Journal of Applied Mathematics. 2012; 2012. Article ID 518684, 7p
Khan Y, Mirzabeigy A. Improved accuracy of He’s energy balance method for analysis of conservative nonlinear oscillator. Neural Computing and Application. 2014; 25(3-4):889-895
Molla MHU, Razzak MA, Alam MS. An analytical technique for solving quadratic nonlinear oscillator. Multidiscipline Modeling in Materials and Structures. 2017; 13(3):424-433
He JH. An improved amplitude-frequency formulation for nonlinear oscillators. International Journal of Nonlinear Science and Numerical Simulation. 2008; 9(2):211-212
He JH. Hamiltonian approach to nonlinear oscillators. Physics Letters A. 2010; 374(23):2312-2314
Farzaneh Y, Tootoonchi AA. Global error minimization method for solving strongly nonlinear oscillator differential equations. Computers and Mathematics with Applications. 2010; 59:288-295
Hosen MA, Chowdhury MSH, Ali MY, Ismail AF. A new analytical technique for solving nonlinear non-smooth oscillators based on the rational harmonic balance method. In: Saian R, Abbas MA, editors. Proceedings of the Second International Conference on the Future of ASEAN (ICoFA) 2017. Vol. 2. 2018
Qian YH, Pan JL, Chen SP, Yao MH. The spreading residue harmonic balance method for strongly nonlinear vibrations of a restrained cantilever beam. Advances in Mathematical Physics. 2017; 2017. Article ID 5214616, 8p
Lai SK, Lim CW. Nonlinear vibration of a two-mass system with nonlinear stiffnesses. Nonlinear Dynamics. 2007; 49:233-249
Gimeno E, Belendez A. Rational-harmonic balancing approach to nonlinear phenomena governed by pendulum-like differential equations. Zeitschrift für Naturforschung. 2009; 64a:819-826
Lim CW, Lai SK, Wu BS. Accurate higher-order analytical approximate solution to large-amplitude oscillating systems with a general non-rational restoring force. Nonlinear Dynamics. 2005; 42:267-281
Bravo YS, Sanchez MA. A weighted mean-square method of “cubication” for nonlinear oscillators. Journal of Sound and Vibration. 1989; 134(3):423-433
Belendez A, Alvarez ML, Fernandez E, Pascual I. Cubication of conservative nonlinear oscillators. European Journal of Physics. 2009; 30(5):973-981
Elıas-Zuniga A, Martınez-Romero O. Accurate solutions of conservative nonlinear oscillators by the enhanced cubication method. Mathematical Problems in Engineering. 2013; 2013:9. Article ID 842423
Big-Alabo A. A simple cubication method for approximate solution of nonlinear Hamiltonian oscillators. International Journal of Mechanical Engineering Education. 2019. DOI: 10.1177/0306419018822489
He JH. Variational approach for nonlinear oscillators. Chaos, Solitons and Fractals. 2007; 34:1430
Ismail GM. An analytical coupled homotopy-variational approach for solving strongly nonlinear differential equation. Journal of the Egyptian Mathematical Society. 2017; 25(4):434-437
Yazdi MK. Approximate solutions to nonlinear oscillations via an improved He's variational approach. Karbala International Journal of Modern Science. 2016; 2:289-297
Ebaid AE, Ali E. On a new aftertreatment technique for differential transformation method and its application to non-linear oscillatory systems. International Journal of Nonlinear Science. 2009; 8(4):488-497
Zhang X, Zou L, Liang S, Lui C. A novel analytic approximate method with a convergence acceleration parameter for solving nonlinear problems. Communications in Nonlinear Science and Numerical Simulation. 2018; 56:354-364
Uhler HS. Period of the bifilar pendulum for large amplitudes. Journal of the Optical Society of America. 1923; 7(3):263-269
Big-Alabo A. Approximate periodic solution to the large-amplitude oscillations of a simple pendulum. International Journal of Mechanical Engineering Education. 2019. DOI: 10.1177/0306419019842298
Big-Alabo A, Harrison P, Cartmell MP. Algorithm for the solution of elastoplastic half-space impact: Force-indentation linearisation method. Journal of Mechanical Engineering Sciences. 2015; 229(5):850-858
Big-Alabo A. Equivalent impact system approach for elastoplastic impact analysis of dissimilar spheres. International Journal of Impact Engineering. 2018; 113:168-179
Big-Alabo A. Continuous piecewise linearization method for approximate periodic solution of the relativistic oscillator. International Journal of Mechanical Engineering Education. 2018. DOI: 10.1177/0306419018812861
Butikov EI. The rigid pendulum–An antique but evergreen physical model. European Journal of Physics. 1999; 20:429-441
Lima FMS. Simple but accurate periodic solutions for the nonlinear pendulum equation. Revista Brasileira de Ensino de Física. 2019; 41(1):e20180202-1-6
Arnold VI. Mathematical Methods of Classical Mechanics. 2nd ed. New York: Springer; 1989. pp. 261-264
Mickens RE. Oscillations in an x4/3 potential. Journal of Sound and Vibration. 2001; 246:375-378
Gottlieb H. Frequencies of oscillators with fractional-power nonlinearities. Journal of Sound and Vibration. 2003; 261(3):557-566
Hu H, Xiong ZG. Oscillations in an x(2m+2)/(2n+1) potential. Journal of Sound and Vibration. 2003; 259:977-980
Belendez A. Homotopy perturbation method for a conservative x1/3 force nonlinear oscillator. Computers and Mathematics with Applications. 2009; 58:2267-2273
Cveticanin L, Kovacic I, Rakaric Z. Asymptotic methods for vibrations of the pure non-integer order oscillator. Computers and Mathematics with Applications. 2010; 60:2616-2628
Cveticanin L. Strong Nonlinear Oscillators: Analytical Solutions. 2nd ed. Switzerland: Springer; 2018