An Exact Solution to the Quadratic Damping Strong Nonlinearity Duffing Oscillator

Salas, Alvaro H.; El-Tantawy, S. A.; Aljahdaly, Noufe H.

doi:https://doi.org/10.1155/2021/8875589

Mathematical Problems in Engineering

On this page

Abstract Introduction Conclusion Appendix Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 8875589 | https://doi.org/10.1155/2021/8875589

An Exact Solution to the Quadratic Damping Strong Nonlinearity Duffing Oscillator

Alvaro H. Salas,¹S. A. El-Tantawy,^2,3and Noufe H. Aljahdaly⁴

Academic Editor: Kai Zhang

Received31 Aug 2020

Revised28 Dec 2020

Accepted31 Dec 2020

Published19 Jan 2021

Abstract

The nonlinear equations of motion such as the Duffing oscillator equation and its family are seldom addressed in intermediate instruction in classical dynamics; this one is problematic because it cannot be solved in terms of elementary functions before. Thus, in this work, the stability analysis of quadratic damping higher-order nonlinearity Duffing oscillator is investigated. Hereinafter, some new analytical solutions to the undamped higher-order nonlinearity Duffing oscillator in the form of Weierstrass elliptic function are obtained. Posteriorly, a novel exact analytical solution to the quadratic damping higher-order nonlinearity Duffing equation under a certain condition (not arbitrary initial conditions) and in the form of Weierstrass elliptic function is derived in detail for the first time. Furthermore, the obtained solutions are camped to the Runge–Kutta fourth-order (RK4) numerical solution.

1. Introduction

In the last few decades, considerable work has been invested in developing new methods for analytical and numerical solutions for several differential equations [1–4] such as strongly nonlinear oscillators, but it is still difficult to obtain convergent results in cases of strong nonlinearity. The Duffing oscillator model is one of the most important and most popular models in dynamic systems due to its importance in explaining many nonlinear phenomena in science and engineering. For instance, the nonlinear Duffing equation was devoted for studying the dynamical behavior of oscillations in plasma physics [5], plasma physics and rigid rotator [6], magneto-elastic mechanical systems [7], nonlinear vibration of beams and plates [8], fluid flow induced vibration [9], large amplitude oscillation of centrifugal governor systems [10], etc. Since the early century so far, a lot of researchers have been devoted carefully to study the Duffing oscillator model in several physical situations [11]. So far the Duffing oscillator model has been an effective and applicable example to nonlinear oscillations in dynamic systems [12]. In addition to its mentioned applications, it also has various applications in the mechanical engineering [9, 10] and electrical engineering [13, 14]. It has been solved numerically and analytically via using various methods in order to recognize its physical characters. Some methods which consider the conservative Duffing osciullators are the harmonic balance technique [15], the frequency-amplitude formulation [16], the homotopy perturbation method [10, 17], the coupled homotopy-variational approach [18], and the differential transform [19]. It is known that if the damping effect is taken into account in the conservative Duffing oscillators, the model becomes a nonconservative system and the oscillation amplitude reduces over time. There are many methods that have been devoted to solve the nonconservative Duffing oscillators, whether they are semianalytical or numerical solutions [20–26]. However, most of the analytical methods are unable to deal with the Duffing equation in the presence of a damping effect. Therefore, all published solutions are either approximate numerical solutions, semianalytical or approximate analytical solutions, or analytical solutions under certain restrictions [20–31].

The most popular formula for conservative Duffing oscillator for arbitrary initial conditions is given by [32, 33]

This equation is completely integrable, and its exact analytical solutions in the Jacobi elliptic function forms have been obtained [32, 33]. However, equation (1) without damping effect is not realistic physically, and if the damping effect is taken into account, the following nonconservative Duffing oscillator for arbitrary initial conditions may be obtained as [5, 6, 34–37]

In generally, equation (2) is not integrable differential equation due to the presence damping term , and many methods have been devoted to solve this equation, as was previously mentioned in all the upper references such as [5, 6, 34–37].

The generalized form of the Duffing oscillator equation which has been modelled in an arbitrary form of the nonlinear differential equation may be written as follows:

In equation (3), the second term refers to the nonlinear damping term which is proportional to the th power of and is a positive or negative damping parameter where its sign depends on the type of the physical application. The third term in equation (3), i.e., involves the linear term and the nonlinear terms of which is known as the restoring force of the dynamical system. As we mentioned earlier, many authors studied the linear damping for linear and cubic restoring forces and [5, 6, 34–37]. On the contrary, there are a few attempts to solve and investigate the damping Duffing oscillator for quadratic restoring force [38–40]. The elliptic-harmonic balance method has been applied to find an approximate analytical solution in the Jacobi elliptic function form to a quadratic damping system with linear restoring force: [38]. In [40], some solutions to the quadratic damping Helmholtz–Duffing equation in the form of classical Jacobi elliptic functions have been discussed in details. Given the difficulty of solving quadratic damping Duffing equation with higher-order nonlinearity, this problem remains one of the most important hot topics that require a more careful study to obtain more accurate solutions.

In this paper, the quadratic damping Duffing oscillator with linear and cubic restoring forces is considered under the initial conditions and . Thus, the mathematical model of this physical situation may be written as the initial value problem (IVP):

The IVP (4) is considered a nonintegrable differential equation but may admit exact analytical solution under a certain condition. In this study, we will derive this condition in detail in order to find an exact solution to the IVP (4) in terms of the Weierstrass elliptic function for the first time. Moreover, before finding the exact solution to the IVP (4), we will study its stability in order to analyze its characteristics to recognize its dynamical behavior.

The rest of the paper is presented in the following points. The stability analysis of the quadratic damping higher-order nonlinearity Duffing equation will be investigated in Section 2. Some possible analytical solutions to the Duffing equation and its family in the form of Jacobi elliptic function and Weierstrass elliptic function will be discussed in detail. Firstly, we will discuss some possible solutions to the IVP (2) in terms of Jacobi elliptic function depending on the discriminant of the IVP (2). Thereafter, the exact analytical solutions to both the IVP (2) and IVP (4) in terms of Weierstrass elliptic function will be derived and obtained in details. The summary of our investigations will be introduced in Section 4.

2. Stability Analysis to the Quadratic Damping Duffing Equation

The first equation in the IVP (4) may be written as a dynamical system as follows:

The first integral reads

System (5) becomes Hamilton if the following condition is fulfilledwith

It is noted that condition (7) is fulfilled; thus, system (5) becomes Hamilton with the following Hamiltonian:

Now, if , then system (5) has three equilibrium points:and, for , there is only one equilibrium point:

The Jacobian matrix for system (5) at equilibrium point reads

Thus, the determinant and trace of the Jacobian matrix may be estimated as

The behavior of the dynamical system near the equilibria is described by the Poincaré Diagram as given in [41].

Let us analyze the equilibrium points. For , we obtain The origin becomes a central point if and a saddle if , as shown in Figures 1(a) and 1(b), respectively. For and or , we obtain It is clear that if , the equilibrium point becomes a central point and a saddle point is covered if , as illustrated in Figures 2(a) and 2(b), respectively.

(a)

(b)

(a)

(b)

3. Analytical Solutions to the Quadratic Damping Duffing Equation

Before embarking on solving the IVP (4), we will discuss some possible solutions to both the IVP (1) and the IVP (2) in terms of Weierstrass elliptic function.

3.1. The Analytical Solution of Undamping Duffing Equation in Terms of Weierstrass Elliptic Function

Let us solve the following undamped IVP:where , and the primes represent the differentiation with respect to the independent variable say “.” Here, it is assumed that , , , and are real numbers. To find an analytical solution to the IVP (16) in the form of Weierstrass elliptic function , the following ansatz is assumed:where , , , , , and are some constants to be determined later and the condition must be fulfilled. Note that the Weierstrass elliptic function satisfies the following relations:

Substituting both equations (17) and (18) into equation (16) and after some simplifications, we finally obtainwhere the values of the coefficients are given in Appendix A. Equating the coefficients to zero, we get a system of algebraic equations, and by solving this system, we obtain

Thus, the general exact solution to the undamping strong nonlinear Duffing equation in terms of Weierstrass elliptic function is given by

Inserting the initial conditions and into solution (21), we finally get the value of as follows:and the value of may be determined from the conditionand this condition is equivalent to

The roots of equation (24) give the value of as follows:where .

Remark 1. Solution (21) is periodic with period:where is the greatest real root to the cubic

3.2. An Analytical Solution of Quadratic Damping Duffing Equation in Terms of Weierstrass Elliptic Function

Here, to find an exact solution to the IVP (4) in terms of Weierstrass elliptic function, the following ansatz is introduced:where .

Inserting ansatz (28) and equation (18) into , we get , withwhere the values of the coefficients are given in Appendix B. Equating the coefficients to zero gives an algebraic equation system, and by solving this system, we obtainwith

The value of may be determined from the following cubic equation:

We choose the first root of equation (35) as the default value to or any real root. Now, we must determine the value of so that

Let

Making use of the addition formula for the Weierstrass function giveswhere the values of , , and are defined in equations (30), (31), and (35).Using the initial conditions given in equation (36) will give us

By solving system (39), we can get the values of and as

Note that the value of must also satisfy the condition

The resultant of the polynomials given in equations (35) and (41) with respect to reads

This resultant must vanish in order to get exact solutions to the problem under consideration.

We may consider the condition

From this condition, we can calculate the value of as

Thus, the IVP (4) has been solved analytically under restriction (43) or (44).

A comparison between the analytical solution (38) and the approximate numerical RK4 solution is introduced in Figure 3. It is observed that our solution is in agreement with the numerical simulation solution.

In the limiting case, i.e., for , the exact solution to the undamping Duffing problem (1) for arbitrary initial conditions may be coveredwhere

Note that both solutions (21) and (45) are equivalent.

4. Conclusion

The equation of motion with higher-order nonlinear restoring force and higher-order nonlinear damping force in which the damping is proportional to the square of the velocity has been carefully studied. As an example of this type of equation of motion, the stability analysis of the quadratic damping higher-order nonlinear Duffing oscillator equation using the Bifurcation method has been investigated. Then comes the main goal, which is to find some exact solutions to the Duffing oscillator equation and its family. In the form of Weierstrass elliptic function, the exact solution of the undamping higher-order nonlinear Duffing equation has been discussed. As for the second and most important goal of this study, it is to obtain an exact solution to the quadratic damping higher-order nonlinear Duffing oscillator equation. In general, the damping higher-order nonlinear Duffing oscillator equation is nonintegrable differential equation but under a certain condition it becomes integrable. Accordingly, the condition that makes this equation integrable has been deduced, and then, this solution has been obtained in detail in terms of the Weierstrass elliptic function. Furthermore, a comparison between the obtained solutions and the numerical RK4 solution to the Duffing equation and its family has been carried out. Also, in the limiting case, i.e., for , the exact solution of the quadratic damping higher-order nonlinear Duffing oscillator equation has been reduced to the exact solution of the undamping higher-order nonlinear Duffing equation.

In future work, the Weierstrass elliptic functions will be devoted to find some exact and approximate analytical solutions to the forced higher-order damping cubic-quintic Duffing oscillator with external constant and periodic forces.

Appendix

A. Equation (A.1)

B. Equation (B.1)

Data Availability

The datasets generated for this paper are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge the Deanship of Scientific Research (DSR) in King Abdulaziz University for its support.

References

A. M. Wazwaz, Partial Differential Equations and Solitary Waves Theory, Higher Education Press, Beijing, China, 2009.
A. M. Wazwaz, Partial Differential Equations: Methods and Applications, Balkema, cop, Lisse, Netherlands, 2002.
M. B. Arain, M. M. Bhatti, A. Zeeshan, T. Saeed, and A. Hobiny, “Analysis of arrhenius kinetics on multiphase flow between a pair of rotating circular plates,” Mathematical Problems in Engineering, vol. 2020, Article ID 2749105, 17 pages, 2020.
View at: Publisher Site | Google Scholar
M. M. Bhatti and I. Sara, “Abdelsalam, thermodynamic entropy of a magnetized Ree-Eyring particle-fluid motion with irreversibility process: a mathematical paradigm,” Zeitschrift für Angewandte Mathematik und Mechanik, Article ID e202000186, 2020.
View at: Publisher Site | Google Scholar
A. H. Salas and S. A. El-Tantawy, “On the approximate solutions to a damped harmonic oscillator with higher-order nonlinearities and its application to plasma physics: semi-analytical solution and moving boundary method,” The European Physical Journal Plus, vol. 135, pp. 833–917, 2020.
View at: Publisher Site | Google Scholar
M. Abu Hammad, A. H. Salas, and S. A. El-Tantawy, “New method for solving strong conservative odd parity nonlinear oscillators: applications to plasma physics and rigid rotator,” AIP Advances, vol. 10, Article ID 085001, 11 pages, 2020.
View at: Publisher Site | Google Scholar
J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems and Bifurcations of Vector Fields, Springer-Verlag, Heidelberg, Germany, 1983.
F. Bakhtiari-Nejad and M. Nazari, “Nonlinear vibration analysis of isotropic cantilever plate with viscoelastic laminate,” Nonlinear Dynamics, vol. 56, no. 4, pp. 325–356, 2009.
View at: Publisher Site | Google Scholar
N. Srinil and H. Zanganeh, “Modelling of coupled cross-flow/in-line vortex-induced vibrations using double Duffing and van der pol oscillators,” Ocean Engineering, vol. 53, pp. 83–97, 2012.
View at: Publisher Site | Google Scholar
D. Younesian, H. Askari, Z. Saadatnia, and M. K. Yazdi, “Periodic solutions for nonlinear oscillation of a centrifugal governor system using the he’s frequency–amplitude formulation and he’s energy balance method,” Nonlinear Science Letters A, vol. 2, pp. 143–148, 2011.
View at: Google Scholar
I. Kovacic and M. J. Brennan, The Duffing Equation: Nonlinear Oscillators and Their Behaviour, John Wiley & Sons, Ltd., Hoboken, NJ, USA, 1st edition, 2011.
P. S. Landa, Nonlinear Oscillations and Waves in Dynamical Systems, Springer, Heidelberg, Germany, 1996.
U. R. Singh and G. Purohit, “Synchronization of quadratically damped mathieu– duffing based chaotic circuits with applications to secure communications,” International Journal of Science and Research, vol. 5, pp. 80–85, 2016.
View at: Google Scholar
U. R. Singh, G. Purohit, A. Sharma, and V. Patidar, “An analogue circuit to study the forced and quadratically damped Duffing oscillator,” International Journal of Nonlinear Dynamics and Control, vol. 1, no. 1, pp. 87–96, 2017.
View at: Publisher Site | Google Scholar
R. E. Mickens, “Mathematical and numerical study of the duffing-harmonic oscillator,” Journal of Sound and Vibration, vol. 244, no. 3, pp. 563–567, 2001.
View at: Publisher Site | Google Scholar
J. H. He, “Comment on he’s frequency formulation for nonlinear oscillators,” European Journal of Physics, vol. 29, pp. 19–22, 2008.
View at: Publisher Site | Google Scholar
J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.
View at: Publisher Site | Google Scholar
Y. Khan, M. Akbarzade, and A. Kargar, “Coupling of homotopy and the variational approach for a conservative oscillator with strong odd-nonlinearity,” Scientia Iranica, vol. 19, no. 3, pp. 417–422, 2012.
View at: Publisher Site | Google Scholar
M. El-Shahed, “Application of differential transform method to non-linear oscillatory systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, pp. 1714–1720, 2008.
View at: Google Scholar
B. S. Wu and W. P. Sun, “Construction of approximate analytical solutions to strongly nonlinear damped oscillators,” Archive of Applied Mechanics, vol. 81, no. 8, pp. 1017–1030, 2011.
View at: Publisher Site | Google Scholar
L. Cveticanin, “Oscillators with nonlinear elastic and damping forces,” Computers & Mathematics with Applications, vol. 62, no. 4, pp. 1745–1757, 2011.
View at: Publisher Site | Google Scholar
A. Y. T. Leung, Z. Guo, and H. X. Yang, “Residue harmonic balance analysis for the damped duffing resonator driven by a van der pol oscillator,” International Journal of Mechanical Sciences, vol. 63, no. 1, pp. 59–65, 2012.
View at: Publisher Site | Google Scholar
M. Turkyilmazoglu, “An effective approach for approximate analytical solutions of the damped duffing equation,” Physica Scripta, vol. 86, no. 6, Article ID 015301, 2012.
View at: Publisher Site | Google Scholar
Y. Khan and Q. Wu, “Homotopy perturbation transform method for nonlinear equations using he’s polynomials,” Computers & Mathematics with Applications, vol. 61, pp. 1963–1967, 2011.
View at: Google Scholar
Y. Khan and F. Austin, “Application of the laplace decomposition method to nonlinear homogeneous and non-homogenous advection equations,” Zeitschrift für Naturforschung A, vol. 65, no. 10, pp. 849–853, 2010.
View at: Publisher Site | Google Scholar
S. Nourazar and A. Mirzabeigy, “Approximate solution for nonlinear duffing oscillator with damping effect using the modified differential transform method,” Scientia Iranica, vol. 20, pp. 364–368, 2013.
View at: Google Scholar
J.-H. He, “Preliminary report on the energy balance for nonlinear oscillations,” Mechanics Research Communications, vol. 29, no. 2-3, pp. 107–111, 2002.
View at: Publisher Site | Google Scholar
S. Durmaz, S. A. Demirbag, and M. O. Kaya, “High order he’s energy balance method based on collocation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 11, no. 1-5, 2010.
View at: Publisher Site | Google Scholar
M. K. Yazdi, H. Ahmadian, A. Mirzabeigy, and A. Yildirim, “Dynamic analysis of vibrating systems with nonlinearities,” Communications in Theoretical Physics, vol. 57, no. 2, pp. 183–187, 2012.
View at: Publisher Site | Google Scholar
A. Elias-Zungia, “Analytical solution of the damped helmholtz–duffing equation,” Applied Mathematics Letters, vol. 25, pp. 2349–2353, 2012.
View at: Google Scholar
Y. Geng, “Exact solutions for the quadratic mixed-parity helmholtz–duffing oscillator by bifurcation theory of dynamical systems,” Chaos, Solitons and Fractals, vol. 81, pp. 68–77, 2018.
View at: Google Scholar
E. T. Whittaker and G. N. Watson, A Course of Modern Analysis, Cambridge University Press, Cambridge, UK, 4th edition, 1980.
M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, Dover, New York NY, USA, 9th edition, 1980.
K. Johannesen, “The duffing oscillator with damping,” European Journal of Physics, vol. 36, no. 13, Article ID 065020, 2015.
View at: Google Scholar
K. Johannesen, “The duffing oscillator with damping for a softening potential,” International Journal of Applied and Computational Mathematics, vol. 3, pp. 3805–3816, 2017.
View at: Google Scholar
J.-w. Zhu, “A new exact solution of a damped quadratic non-linear oscillator,” Applied Mathematical Modelling, vol. 38, no. 24, pp. 5986–5993, 2014.
View at: Publisher Site | Google Scholar
B. Baumann, J. Schwieger, M. Wolff, F. Manders, and J. Suijker, “Nonlinear behavior in high-intensity discharge lamps,” Journal of Physics D Applied Physics, vol. 49, Article ID 255201, 2015.
View at: Google Scholar
L. Cveticanin, “Oscillator with strong quadratic damping force,” Publications de l’Institut Math?matique (Belgrade), vol. 85, no. 99, pp. 119–130, 2009.
View at: Publisher Site | Google Scholar
R. Porwal and N. S. Vyas, “Damped quadratic and mixed-parity oscillator response using krylov-bogoliubov method and energy balance,” Journal of Sound and Vibration, vol. 309, no. 3-5, pp. 877–886, 2008.
View at: Publisher Site | Google Scholar
S. K. Lai and K. W. Chow, “Exact solutions for oscillators with quadratic damping and mixed-parity nonlinearity,” Physica Scripta, vol. 85, Article ID 045006, 2012.
View at: Publisher Site | Google Scholar
B. R. Smith, “The quadratically damped oscillator: a case study of a non-linear equation of motion,” American Journal of Physics, vol. 80, no. 9, p. 816, 2012.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Alvaro H. Salas et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1423

Downloads

907

Citations