Open access peer-reviewed chapter

Periodic Solution of Nonlinear Conservative Systems

By Akuro Big-Alabo and Chinwuba Victor Ossia

Submitted: December 23rd 2018Reviewed: October 26th 2019Published: January 15th 2020

DOI: 10.5772/intechopen.90282

Downloaded: 346


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. Introduction

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 [1], human eardrum [2], relativistic oscillator [3], planetary orbits around the sun [3], and current-carrying conductor in the electric field of an infinite rod [4]. 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 (T) and potential (V) energies except at critical points where the total energy may be only kinetic (Tmax) or potential (Vmax). Generally, it is expected that T=Tqq̇and V=Vqq̇, where qis the generalized displacement. Naturally, qand q̇are not expected to form a product in the function Tqq̇, 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 Hamiltonian systemsbecause they admit a Hamiltonian function (Hqq̇=Tqq̇+Vqq̇) that is always constant at any point in time. While this definition of Hamiltonian systems is a physical one, a mathematical definition has been discussed by Jordan and Smith [3]. Examples of Hamiltonian systems include mass-spring oscillator, simple pendulum, and a mass attached to the mid-point of an elastic spring. On the other hand, there are conservative systems in which the kinetic energy cannot be expressed as a pure quadratic function of the velocity because the kinetic energy expression contains a product of velocity and displacement. This second group of conservative systems is referred to as nonnatural because their kinetic energy is not a pure quadratic function of velocity. Although the total energy in such systems is conserved, their Hamiltonian function (Hqq̇) is not constant [4]. Hence, the nonnatural conservative systems may be referred to as non-Hamiltonian conservative systems. Examples of this category of conservative systems abound in artificial systems and include slider-crank mechanism [1], particle sliding on a vertical rotating parabola [4], pendulum attached to massless rolling wheel [4], rigid rod rocking on a circular surface without slip [4], and circular sector oscillator [5]. An important quality of the non-Hamiltonian conservative systems is that their vibration equation, which is normally derived by the Lagrangian approach, does not conform to the standard representation of conservative systems that clearly shows the restoring force. Rather, the derived vibration equation has a quadratic velocity term, which represents a coordinate-dependent parameter rather than a dissipative parameter.

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) [6], 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 classical perturbation methods. The perturbation methods are formulated based on the concept that an unknown nonlinear system can be studied by introducing a small disturbance to a known linear system in equilibrium. For this reason, the classical perturbation methods (see Nayfeh [9] for a comprehensive treatment of classical perturbation methods) depend on the assumption of a small parameter. The problem with the small parameter assumption is that it has a small range of validity and only produces reliable solutions for cases of small-amplitude oscillations and weak nonlinearity. Nevertheless, the classical perturbation methods are still very relevant today for introducing and investigating various nonlinear concepts.

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 [10], Homotopy perturbation method [11] and its variants [12, 13, 14, 15, 16, 17], modified Lindstedt-Poincare method [18, 19, 20, 21], book-keeping parameter method [22], iteration perturbation method [23], parameterized perturbation method [24], perturbation incremental method [25], and linearized perturbation method [26]. A review article on some of the recent perturbation methods has been published by He [27]. 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 [28]. 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 [29], Homotopy analysis method [30], Variational iteration method [31], Energy balance method [2] and its modifications [32, 33, 34], He Chengtian’s interpolation method [27] also called max-min approach, amplitude-frequency formulation [35], Hamiltonian approach [36], global error minimization method [37], Harmonic balance method [4] and its modifications [38, 39, 40, 41, 42], cubication methods [43, 44, 45, 46], variational methods [47, 48, 49], differential transform method [50], and continuous piecewise linearization method [8]. Nonperturbation methods also have various limitations. For instance, a study [51] 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 [1] and the bifilar pendulum [52]. Furthermore, it has been observed that higher order estimates do not always improve the solution of the oscillation history [27]. 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 [8].

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 [53], 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 [55]. 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 nequal 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].

Figure 1.

Discretization of the restoring force of a typical nonlinear oscillator.

2.2 Mathematical formulation of the continuous piecewise linearization method

The standard form for representation of a nonlinear conservative oscillator moving in the u-direction is given as:


where fuis 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 sthdiscretization represents a general discretization with start point at rand endpoint at s=r+1. 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 sthdiscretization can be expressed as:


where Krs=fusfur/usuris the linear slope of Frsuand Fr=fur.

Since Frsuis an approximation of fufor the sthdiscretization, then substituting Eq. (2) in (1) gives the approximate equation of motion for each discretization as follows:


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 Krs>0, the solution for the displacement and velocity can be expressed as:


where ωrs=Krs, Crs=urFr/Krs, and Rrs=urCrs2+u̇r/ωrs21/2. The initial conditions and other parameters are determined based on the oscillation stage. For the oscillation stage that moves from +Ato A, the initial conditions for each discretization are ur=ur0=Aruand u̇r=u̇r0=2Aurfudu, where u=A/n, and the other parameters are calculated as:


For the oscillation stage that moves from Ato +A, the initial conditions are ur=ur0=A+ruand u̇r=u̇r0=2Aurfudu; the other parameters are calculated as:


The time at the end of each discretization is ts=tr+t, and the end conditions usand u̇sare calculated by replacing rwith sin the formulae for the initial conditions.

2.2.2 Solution for negative linearized stiffness

If Krs<0, the solution for the displacement and velocity can be expressed as follows:


where ωrs=Krs; Crs=ur+Fr/Krs. Applying the initial conditions to Eqs. (7a) and (7b) gives: Ars=12ur+u̇r/ωrsCrs; Brs=12uru̇r/ωrsCrs. The initial and end conditions are determined in the same way as for Krs>0above. 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 +Ato Aand vice versa. We note that if u̇r=0, then Ars=Brs=12urCrsand




2.2.3 Solution for zero linearized stiffness

In very rare situations, we may have Krs=0for one or two discretization around the turning points or relatively flat regions of the restoring force. This is likely when uis very small, i.e., for very large n, and can be eliminated by increasing or decreasing nslightly. However, if we want to account for Krs=0, then we get [53]:


where Grs=u̇r+Frtrand Hrs=uru̇rtr12Frtr2. Hence, the time interval is derived from Eq. (11) as:


2.3 Remarks on the CPLM algorithm

  1. 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.

  2. The CPLM is inherently stable and does not have convergence issues [8].

  3. For few discretization, say n10, 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.

  4. 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 A. This means that there will be 2ndiscretizations from Ato A.

  5. The CPLM algorithm retains the same simplicity in implementation irrespective of the complexity of the restoring force. Only the Krsconstant and the integral of the restoring force are to be evaluated anew for any oscillator.

  6. 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.

  7. 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.

  8. For few discretization, say n20, 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 n50, there is no need to extract values from any discretization.

  9. 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 [57] 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 uis the angular displacement, Ω0=g/l, lis the length of the pendulum, and g=9.8m/s2. The initial conditions are given as u0=Aand u̇0=0. 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 [58]:


Where snis the Jacobi elliptic sine function, k=sinA/2, and Kk2is the complete elliptic integral of the first kind given as:


From Eq. (13), the restoring force for the pendulum is fu=Ω02sinuand looks like the plot in Figure 1. This means that Krs=Ω02sinussinur/usurand u̇=±Ω02cosucosA. The initial and final velocities for each discretization are determined using the expression for u̇, and all the other constants for the solution of each discretization are determined based on Krs. A plot of the frequency-amplitude response for the simple pendulum when Ω0=1.0is 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 A178°, the maximum error in the CPLM estimate is less than 0.45% for n=50and 0.14% for n=100. Also, the oscillation history for moderate-amplitude (A=45°) and large-amplitude (A=135°) is shown in Figure 3 and an excellent agreement between the CPLM solution for n=100and 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., n=50100) are required to get an accuracy that is within 1.0% of the exact solution during large-amplitude oscillations (90°<A<180°) of the simple pendulum.

Figure 2.

(a) Frequency-amplitude response of the simple pendulum for0°<A<180°. (b) CPLM error analysis.

Figure 3.

Oscillation history of the simple pendulum for (a)A=45°and (b)A=135°. CPLM estimate—Lines; exact solution—Markers.

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 2dand the coordinate of the satellite motion is u, then the equation of motion of the satellite is given as [3]:


where Mis the mass of a star and the restoring force is fu=2Mu/d2+u23/2. 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 [3] and Arnold [59], but here, the periodic solution was investigated.

Figure 4.

Restoring force forEq. (18):M=105kg;d=100m;A=500.

The main CPLM constant is Krs=2Md2+us23/2d2+ur23/2and the velocity was derived as: u̇=±2Md2+u21/2d2+A21/2. 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 M=105kgand d=1000m. 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 n=10and 0.20% for n=20. Significantly higher accuracies can be achieved by increasing n, but the results show that n=10gives reasonably accurate estimates.

Figure 5.

(a) Frequency-amplitude response for satellite. (b) CPLM error analysis.

Figure 6.

Oscillation history of satellite for (a)A=50and (b)A=1500. CPLM estimate—Lines; exact solution—Markers.

On the other hand, Figure 6 shows the oscillation history of the satellite during small-amplitude (A=50) and large-amplitude (A=1500) oscillations. The former gives a simple harmonic response with a natural frequency that is independent of the amplitude and approximately equal to 2M/d3, while the latter exhibits an anharmonicresponse with a natural frequency that depends strongly on the amplitude.

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 [65]:


where the restoring force, fu=u1/2m+1, has a fractional index for all mϵ+. The main CPLM constant is evaluated as Krs=us1/2m+1ur1/2m+1/usurand the velocity as u̇=±2m+1/m+11/2A2m+2/2m+1u2m+2/2m+1. The periodic solution for the case of m=1, i.e., u1/3oscillator, is shown in Figures 7 and 8. The exact frequency-amplitude response used for the verification of the CPLM solution is [63]:


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.

(a) Frequency-amplitude response foru1/3oscillator. (b) CPLM error analysis.

Figure 8.

Oscillation history ofu1/3oscillator for (a)A=0.01and (b)A=10. CPLM estimate—lines; exact solution—markers.

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 n=10and 0.076% for n=20. The error in the CPLM estimate is approximately constant for all amplitudes, and the maximum error is well below 1.0% for n=10. The results also reveal that the frequency approaches zero as A. In Figure 8, the oscillation history for small-amplitude (A=0.01) and large-amplitude (A=10.0) oscillations is shown to exhibit similar anharmonic response, which is an indication of strong nonlinearity. Therefore, it can be concluded that the u1/3oscillator 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 [1].

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 y=qu2for q>0and rotating at a constant speed, Ω, about the y-axis (Figure 9a). The u-axis represents the perpendicular displacement of the mass from the y-axis. The kinetic and potential energies of the system are given as [4]:


Hence, the Lagrangian is:


Figure 9.

(a) Schematic of particle on a rotating parabola. (b) Restoring force whenq=1.0;Λ=10.0;A=2.0.

Next, we substitute Eq. (21) into the Euler–Lagrange equation to derive the equation of motion. The Euler–Lagrange equation can be written as:


Therefore, using Eqs. (21) and (22), the motion of a particle on a rotating parabola is governed by:


where Λ=2gqΩ2and the initial conditions are: u0=Aand u̇0=0.

To solve Eq. (23) using the CPLM, it must be recast in the form of Eq. (1). The conservation of energy for Eq. (23) is given as [4]:


where his a constant representing the total energy in the system. Eq. (24) confirms that the Hamiltonian, H=T+V=12m1+4q2u2u̇2+Λ+2Ω2u2=12mh+2Ω2u2,is not constant. Applying the initial conditions, we get h=ΛA2so that u̇2=ΛA2u2/1+4q2u2. Substituting this expression for u̇2in Eq. (24) and after algebraic simplification, we get:


Therefore, the restoring force is fu=Λ1+4q2A2u/1+4q2u22. Figure 9b shows that fuis linear at small amplitudes and strongly nonlinear at large amplitudes. The main CPLM constant was calculated as Krs=Λ1+4q2A21/1+4q2us221/1+4q2ur22, and the velocity was evaluated as u̇=±ΛA2u2/1+4q2u2. The exact time period for this oscillator can be derived in terms of elliptic function as follows [4]:


where Ek2=0π/21k2sin2ϕ1/2is the complete elliptic integral of the second kind and k2=4q2A2/1+4q2A2. Then, the exact frequency was computed as ωex=2π/Tex, 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 A=0.50and A=2.0are shown in Figure 11. As demonstrated in [4], periodic solutions for this system exist only for Λ>0. Hence, the simulations in Figures 10 and 11 were conducted for Λ=10and q=1.0. An excellent agreement is observed between the CPLM estimates and the exact results. For 0<A20, the maximum error in the CPLM estimate of the frequency-amplitude response is 0.642% for n=10and 0.101% for n=20, 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.

Figure 10.

(a) Frequency-amplitude response for particle on a rotating parabola. (b) CPLM error analysis.

Figure 11.

Oscillation history of particle on a rotating parabola for (a)A=0.50and (b)A=2.0. CPLM estimate—Lines; exact solution—Markers.

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 n20, 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.

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Akuro Big-Alabo and Chinwuba Victor Ossia (January 15th 2020). Periodic Solution of Nonlinear Conservative Systems, Progress in Relativity, Calin Gheorghe Buzea, Maricel Agop and Leo Butler, IntechOpen, DOI: 10.5772/intechopen.90282. Available from:

chapter statistics

346total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes

By Joás Venâncio and Carlos Batista

Related Book

First chapter

Complementarity in Quantum Mechanics and Classical Statistical Mechanics

By Luisberis Velazquez Abad and Sergio Curilef Huichalaf

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us