Sound & Vibration |

DOI: 10.32604/sv.2022.014166

ARTICLE

Damped Mathieu Equation with a Modulation Property of the Homotopy Perturbation Method

Department of Mathematics, Faculty of Education, Ain Shams University, Roxy, Cairo, Egypt

*Corresponding Author: Nasser S. Elgazery. Email: nasser522000@gmail.com

Received: 06 September 2020; Accepted: 25 June 2021

Abstract: In this article, the main objective is to employ the homotopy perturbation method (HPM) as an alternative to classical perturbation methods for solving nonlinear equations having periodic coefficients. As a simple example, the nonlinear damping Mathieu equation has been investigated. In this investigation, two nonlinear solvability conditions are imposed. One of them was imposed in the first-order homotopy perturbation and used to study the stability behavior at resonance and non-resonance cases. The next level of the perturbation approaches another solvability condition and is applied to obtain the unknowns become clear in the solution for the first-order solvability condition. The approach assumed here is so significant for solving many parametric nonlinear equations that arise within the engineering and nonlinear science.

Keywords: Damped Mathieu Equation; parametric nonlinear oscillator; resonance instability; homotopy perturbation method (HPM)

In wide engineering and physical applications, the nonlinear oscillators exist. Also, the parametric excitation takes place when a modifying physical parameter, such as a moment of stiffness or inertia, acts in a forcing model. This excitement yields a variable time coefficient, commonly an oscillation, in the governing system of motion. On the other hand, an external excitation outcome acting as an inhomogeneous part in the model of motion. Furthermore, minor parametric excitement produces a major response when the frequency of the excitement is far from the fundamental resonance, as shown in [1–5].

A classical example of parametric excitation is the swinging pendulum with oscillating support. The equation of motion describing the model is the well-known Mathieu equation. In 1868 Mathieu studied the vibration of elliptical membranes [6]. Consequently, he introduced the Mathieu equation that is an example of a linear differential equation (LDE) with parametric excitation. The Mathieu equation has application to the dynamics of passive towed arrays in submarines, as well as serving as a useful model for many interesting problems in physics, biology, applied mathematics, and engineering mechanics fields [7]. For some of the non-linear variations of the Mathieu, the equation has been presented in [8,9]. Moreover, the oscillations of the mechanical systems under the action of an oscillatory external force may reveal a Duffing problem, for instance, see references [10–14]. Recently, Moatimid [15] attempted to study the stability analysis of a parametric Duffing oscillator. In this investigation Moatimid showed that the cubic stiffness parameter and the damped parameter have a destabilizing influence, however, the parametric and natural frequencies are of stabilizing influences.

The main target in the present work is how to achieve accurate approximate solutions of the nonlinear oscillators with highly strong nonlinearity. In recent centuries, many analytical approaches were developed to work out the periodic motion of nonlinear oscillators, such as the averaging method, perturbation methods, harmonic balance method, and the generalized harmonic method. The classical perturbation procedure depends on small parameters and chooses unsuitable small parameters that can lead to wrong solutions [16]. Therefore, a new perturbation technique was first proposed by He et al. [17–20]. This technique is named as the HPM, which represents a combination of the Homotopy analysis and classical perturbation methods. It has a full promise of the traditional perturbation techniques. The major property of the HPM is in its ability and flexibility to deal with many types of linear and nonlinear differential equations conveniently and accurately. Further, the HPM provides us with an appropriate direction to calculate an approximate or an analytic solution to several models arising in different fields. He [21] was built the most two considerable steps in the criteria of the HPM with a suitable initial guess and suggested an alternative approach to the construction of the Homotopy equation. Hence, He applied HPM to solve the Lighthill equation [17], Duffing equation [22], and Blasius equation [23], then the idea goes through and has been applied to solve nonlinear wave equations [17], boundary value problems [20]. Babolian et al. [24] applied the homotopy perturbation method to solve the Burgers, the modified Korteweg-de Vries, and regularized long-wave equations.

On the other hand, the HPM has more improved and developed by many engineers and scientists, for instance, a couple of the Laplace transforms and Homotopy perturbation method was implemented by El-Dib et al. [25]. The HPM with two expanding parameters that efficient for some partial nonlinear equations was suggested by He [26] and El-Dib [27]. Also, El-Dib [28] introduced a modified version of the HPM via the multiple scales technique. This new modification works particularly well for the nonlinear oscillators. Furthermore, away from the traditional HPM, Ren et al. [29] made another couple of the HPM and multiple time scales to become a powerful mathematical tool for many nonlinear equations. They displayed that the present procedure may be further afflicted by incorporating several known technologies. It provides solutions to nonlinear equations, whilst the classical perturbation technique became unsuccessful. Moreover, Rabbani [30] introduced a new homotopy perturbation approach for solving main non-linear models through the projection method. A new homotopy perturbation technique for solving linear and nonlinear Schrödinger equations has been addressed by Ayati et al. [31]. Further, by utilizing the HPM, a novel approach in examining the nonlinear Rayleigh-Taylor instability is conducted by El-Dib et al. [32]. Recently, a periodic solution of the cubic nonlinear Klein–Gordon equation using the He-multiple-scales method has been investigated by El-Dib [33]. Also, El-Dib et al. [34] investigated the impact of fractional derivative properties on the periodic analytic solution of the nonlinear oscillations using the HPM. Moreover, He’s-multiple-scale scale to analyze the cubic-quintic Duffing equation has been analyzed by El-Dib et al. [35]. El-Dib [36] presented a stability approach of a fractional-delayed Duffing oscillator. A Nonlinear Instability of a Cylindrical Interface between two MHD Darcian flows has been studied by Moatimid et al. [37]. Further, El-Dib [38] introduced a modified multiple scale technique for the stability of the fractional delayed nonlinear oscillator. Besides, a periodic solution of the fractional sine-Gordon equation has been studied by Shen et al. [39]. Elgazery [40] applied the HPM to give a periodic solution of the Newell-Whitehead-Segel model. Further, for more very useful modification of the homotopy perturbation approach, Yu et al. [41] introduced HPM with an auxiliary parameter for nonlinear oscillators. Also, HPM for Fangzhu oscillator has been used by He et al. [42]. Finally, for more very useful modification of the homotopy perturbation approach, see [43–45].

Motivated by possibility applications in engineering, biology, and physics, which is based on studying the solution of the damped Mathieu equation. Hence, in the present work, our objective is to apply the HPM to linear or nonlinear equations having periodic coefficients such Mathieu equation which has been of great importance among researchers. For the presentation of this article; the rest of the manuscript is systematized as follows: Section 2 is introducing the HPM for the mathematical formulation. The modulation procedure, in detail, is displayed in Section 3. The non-resonance case, stability analysis of the non-resonance case, the resonance case of

To explain the proposed technique, consider the following parametric pendulum equation as an illustrative example:

where

The homotopy perturbation method can be considered as a combination of the classical perturbation technique and the homotopy (whose origin is in the topology [46]), but not restricted to the limitations of traditional perturbation methods. For example, this method does require neither small parameter nor linearization and only requires little iteration to obtain accurate solutions [17] and [18].

We define the two parts of Eq. (2) as

where,

Construct the homotopy statement as

As in He’s a homotopy perturbation method, it is obvious that when

According to linear differential equations theory, the general solution of Eq. (5) is expressed in terms of two linearly independent solutions, say,

where

It can be noticed that the homotopy function (7) is essentially the same as (4), except for the function

where

Eq. (9) can be satisfied by

Substituting (12) into Eq. (10) gets the form

Before analyzing the first-order problem, we must distinguish between the two cases. The case of the frequency

For arbitrary frequency

To obtain uniform expansions for problems of this kind, the expansion (8) needs to be modified. If we modulate the initial solution (6) so that the constant

Then Eq. (12) in the modulate case becomes

where

Consequently, the homotopy state, Eq. (7), in the modulated form becomes

It is convenient to choose the modulated function

The function

where

where cc. indicates to the complex conjugate for the preceding terms and

It is noted that:

and

where dots indicate differentiation concerning the time

Eq. (24) remains to obey the same homotopy concept because it’s become the same harmonic Eq. (5) as

In the light of Eq. (19), the modulate homotopy Eq. (24) will be expanded as a power series in

It is noted that Eq. (25) has been satisfied by Eq. (16) and the zero-order solution for Eq. (17) as approved in Eq. (15). Substituting Eq. (16) into Eq. (26) becomes

This equation contains secular terms at the non-resonance case and other secular terms when the applied frequency

The analysis in this case concerned with the arbitrary chosen for the applied frequency

with its complex conjugate one. This leads to obtaining the valid function

Consequently, the solution of the first-order problem is formulated as

Substituting Eqs. (16) and (30) into Eq. (27), using Eq. (29), yields

The valid solution requires to be removed the terms that producing unbounded solution. These terms imply the following nonlinear solvability condition:

The second-order solution is found to be

If the accuracy to the second-order perturbation is enough, then the approximate solution at the non-resonance case is formulated by substituting Eqs. (15), (16), (31) and (34) into Eq. (20), and setting

5 Stability Analysis for the Non-Resonance Case

The stability criteria in the non-resonance case can be obtained from solving Eq. (29). One may use the following polar form [16]:

with real the unknown functions

where,

6 The Resonance Case Ω is Near

Return to the first-order problem Eq. (28) and re-analyzed it because of the nearness of

Accordingly, we have

Elimination of secular terms from Eq. (28), because of Eqs. (38) and (39) yields

The first-order solution in this case is

Using Eq. (41) with Eq. (27), we obtain the uniform solution for the second-order problem, and the following solvability is presented:

with its complex conjugate. The valid function

The approximate solution up to the second-order is formulated by substituting from Eqs. (15), (16), (41) and (43) into Eq. (16) gets

7 Stability Analysis of the Linear Mathieu Equation

In the limiting case as

The first-order solvability condition (45) can be used to find the stability picture at the resonance case. The second-order solvability condition (46) can be used to find the value of the detuning parameter

It is easy to show that the Eq. (45) can be satisfied by the form

where, the parameter

The parameter

The use of the first value of

The stability criteria require that the right-hand-side of Eq. (50) be positive, which implies that

Stability condition (51) can be rearranged in powers of the applied frequency

and

The transition curves separating stable state from an unstable state corresponding to

and

8 Stability Analysis for the Nonlinear Case

The first-order solvability condition (40) can be used to find the stability picture at the resonance case. The second-order solvability condition (42) can be used to find the value of the detuning parameter

To relax the periodic term into the Eq. (40) we let

with real functions

In order to solve the above coupled nonlinear Eqs. (57) and (58), we may discuss the behavior at the steady-state response. This case is corresponding to the case of

where

Besides, the constants

Squaring both equations in (62) and adding we get

Combing Eq. (61) with Eq. (63) yields

In order to find a constrain for a bounded solution we may modulate the functions

where the functions

The above system is two coupled linear differential equations of first-order in the two functions

where

where relations (62) are used. This characteristic equation depends on the two related parameters

With the help of the second-order solvability condition (42) one can find an expression for both the unknowns

into the second-order solvability condition (42). Separating the real and imaginary parts, produces the following relations, between the parameters

where relations (62) are used. Removing the parameter

Replacing

This equation gives two values

The stabilization for the problem requires that the right-hand side of Eq. (70) be positive provided that the exponential in Eqs. (68) and (69) has positive values. It is noted that the stability reveals as the coefficient of the periodic term in Eq. (2) tends to zero. The instability arrived as the parameter

Removing the parameter

The transition curves separating the stable state of unstable one are corresponding to

and

Using the definition (38) the above transition curves can be sought within the parameter

To obtain the transition curves, independent of the parameter

It is noted that the instability state lies between the above transition curves. The constant coefficients

The homotopy perturbation method (HPM) is one of a easy, powerful, efficient, and accurate approach for evalueting solutions of a large class of nonlinear equations without the need of a discretization or linearization process. HPM is a combination of the homotopy and perturbation methods. That can take the advantages of the conventional perturbation method and eliminating its restrictions. It yields a rapid convergence of the solution series with a few iterations leading to accurate solutions, and the round-off errors are avoided. In general, this method has been successfully used to solve different kinds of linear and nonlinear problems in engineering and science. So, in the present work, we propose a variation of the homotopy perturbation approach via a modulation method that allows finding analytic solutions for ordinary differential models with periodic coefficients. This article is prepared to analyze a parametrically excited oscillator in the presence of strong cubic nonlinearity. The simplest model of this kind is the Mathieu equation that contains a small parameter [47,48]. The present analysis that employs the homotopy perturbation approach [17], has no dependence on models having a small parameter. Due to the present modulation approach, at each level of perturbation, a solvability condition is enjoined. By solving these solvability conditions drives to examining the stability behavior. In each resonance/non-resonant cases stability conditions, are obtained.

Author Contributions: The first author proposed and developed the mathematical modeling of the problem and examines the theory validation; the second author analyzed the cubic damping Mathieu equation. The manuscript was written through the contribution of all authors. All authors discussed the outcomes, reviewed, and approved the final version of the manuscript.

Funding Statement: The authors received no financial support for this research, authorship and publication of this article.

Conflicts of Interest: The authors declare that there are no competing interests regarding the publication of the present paper.

References

- 1. Rong, H., Xu, W., Fang, T. (1998). Principle response of Duffing oscillator to combined deterministic and narrow-band random parametric excitation. Journal of Sound and Vibration, 210(4), 483–515. DOI 10.1006/jsvi.1997.1325.
- 2. Morrison, T. M. (2016). Three problems in nonlinear dynamics with 2:1 parametric excitation (Ph.D. thesis). Department of Applied Mechanics, University of Glasgow, Scotland, UK.
- 3. Dutta, T. K., Prajapati, P. K. (2016). Some dynamical properties of the Duffing equation. International Journal of Engineering Research and Technology, 5(12), 500–503. DOI 10.17577/IJERTV5IS120339.
- 4. Al-Jawary, M. A., Abd-Al-Razaq, S. G. (2016). Analytic and numerical solution for Duffing equations. International Journal of Basic and Applied Sciences, 5(2), 115–119. DOI 10.14419/ijbas.v5i2.5838.
- 5. Sunday, J. (2017). The Duffing oscillator: Applications and computational simulations. Asian Research Journal of Mathematics, 2(3), 1–13. DOI 10.9734/ARJOM/2017/32011.
- 6. Mclachlan, N. W. (1951). Theory and application of Mathieu functions. London, UK: Oxford University Press.
- 7. Ramani, D. V., Keith, W. L., Rand, R. H. (2004). Perturbation solution for secondary bifurcation in the quadratically-damped Mathieu equation. International Journal of Non-Linear Mechanics, 39(3), 491–502. DOI 10.1016/S0020-7462(02)00218-4.
- 8. Taylor, J. H., Narendra, K. S. (1969). Stability regions for the damped Mathieu equation. SIAM Journal on Applied Mathematics, 17(2), 343–352. DOI 10.1137/0117033.
- 9. Insperger, T., Stépán, G. (2003). Stability of the damped Mathieu equation with time delay. Journal of Dynamic Systems, Measurement, and Control, 125(2), 166–171. DOI 10.1115/1.1567314.
- 10. El-Nady, A. O., Lashin, M. M. A. (2016). Approximate solution of nonlinear Duffing oscillator using Taylor expansion. Journal of Mechanical Engineering and Automation, 6(5), 110–116. DOI 10.5923/j.jmea.20160605.03.
- 11. Feng, Z., Chen, G., Hsu, S. B. (2006). A qualitative study of the damped Duffing equation and applications. Discrete and Continuous Dynamical Systems–Series B, 5(5), 1–24. DOI 10.3934/dcdsb.2006.6.1097.
- 12. Luo, A. C. J., Ma, H. (2018). Bifurcation trees of periodic motions to chaos in a parametric Duffing oscillator. International Journal of Dynamics and Control, 6(2), 425–458. DOI 10.1007/s40435-017-0314-x.
- 13. Michon, G., Manin, L., Parker, R. G., Dufour, R. (2008). Duffing oscillator with parametric excitation: Analytical and experimental investigation on a belt-pulley system. Journal of Computational and Nonlinear Dynamics, 3(3), 31001. DOI 10.1115/1.2908160.
- 14. Zivieri, R., Vergura, S., Carpentier, M. (2016). Analytical and numerical solution to the nonlinear cubic Duffng equation: An application to electrical signal analysis of distribution lines. Applied Mathematical Modelling, 40(21–22), 9152–9164. DOI 10.1016/j.apm.2016.05.043.
- 15. Moatimid, G. M. (2020). Stability analysis of a parametric duffing oscillator. Journal of Engineering Mechanics, 146(5), 5020001. DOI 10.1061/(ASCE)EM.1943-7889.0001764.
- 16. Nayfeh, A. H., Mook, D. T. (1979). Non-linear oscillations. New York, USA: John Wily & Sons.
- 17. He, J. H. (1999). Homotopy perturbation technique. Computational Methods in Applied Mechanics and Engineering, 178(3–4), 257–262. DOI 10.1016/S0045-7825(99)00018-3.
- 18. He, J. H. (2000). A coupling method of homotopy technique and a perturbation technique for non-linear problems. International Journal of Non-Linear Mechanics, 35(1), 37–43. DOI 10.1016/S0020-7462(98)00085-7.
- 19. He, J. H. (2005). Application of homotopy perturbation method to nonlinear wave equations. Chaos, Solitons and Fractals, 26(3), 295–300. DOI 10.1016/j.chaos.2005.03.006.
- 20. He, J. H. (2006). Homotopy perturbation method for solving boundary value problems. Physics Letters A, 350(1–2), 87–88. DOI 10.1016/j.physleta.2005.10.005.
- 21. He, J. H. (2012). Homotopy perturbation method with an auxiliary term. Abstract and Applied Analysis, 2012, 1–7. DOI 10.1155/2012/857612.
- 22. He, J. H. (2003). Homotopy perturbation method: A new nonlinear analytical technique. Applied Mathematics and Computation, 135(1), 73–79. DOI 10.1016/S0096-3003(01)00312-5.
- 23. He, J. H. (2003). A simple perturbation approach to Blasius equation. Applied Mathematics and Computation, 140(2–3), 217–222. DOI 10.1016/S0096-3003(02)00189-3.
- 24. Babolian, E., Saeidian, J., Azizi, A. (2009). Application of homotopy perturbation method to some nonlinear problems. Applied Mathematical Sciences, 3(45), 2215– 2226.
- 25. El-Dib, Y. O., Moatimid, G. M. (2019). 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, 44(7), 6581–6659. DOI 10.1007/s13369-018-03705-6.
- 26. He, J. H. (2014). Homotopy perturbation method with two expanding parameters. Indian Journal of Physics, 88(2), 193–196. DOI 10.1007/s12648-013-0378-1.
- 27. El-Dib, Y. O. (2018). Multi-homotopy perturbations technique for solving nonlinear partial differential equations with Laplace transforms. Nonlinear Science Letters A, 9(4), 349–359.
- 28. El-Dib, Y. O. (2017). Multiple scales homotopy perturbation method for nonlinear oscillators. Nonlinear Science Letters A, 8(4), 352–364.
- 29. Ren, Z. F., Yao, S. W., He, J. H. (2019). He’s multiple scales method for nonlinear vibrations. Journal of Low Frequency Noise, Vibration & Active Control, 38(3–4), 1708–1712. DOI 10.1177/1461348419861450.
- 30. Rabbani, M. (2013). New homotopy perturbation method to solve non-linear problems. Journal of Mathematics and Computer Science, 7(4), 272– 275. DOI 10.22436/jmcs.07.04.06.
- 31. Ayati, Z., Biazar, J., Ebrahimi, S. (2014). A new homotopy perturbation method for solving linear and nonlinear Schrödinger equations. Journal of Interpolation and Approximation in Scientific Computing, 2014(1), 1–8. DOI 10.5899/2014/jiasc-00062.
- 32. El-Dib, Y. O., Moatimid, G. M., Mady, A. A. (2019). A novelty to the nonlinear rotating Rayleigh-Taylor instability. Pramana–Journal of Physics, 93, 82. DOI 10.1007/s12043-019-1844-x.
- 33. El-Dib, Y. O. (2019). Periodic solution of the cubic nonlinear Klein-Gordon equation and the stability criteria via the He-multiple-scales method. Pramana–Journal of Physics, 92, 7. DOI 10.1007/s12043-018-1673-3.
- 34. El-Dib, Y. O., Elgazery, N. S. (2020). Effect of fractional derivative properties on the periodic solution of the nonlinear oscillations, Fractals. Fractals, 28(7), 2050095. DOI 10.1142/S0218348X20500954.
- 35. El-Dib, Y. O., Mady, A. A. (2020). He’s multiple-scale solution for the three-dimensional nonlinear KH instability of rotating magnetic fluids. International Annals of Science, 9(1), 52–69. DOI 10.21467/ias.9.1.52-69.
- 36. El-Dib, Y. O. (2020). Stability approach of a fractional-delayed duffing oscillator. Discontinuity Nonlinearity and Complexity, 9(3), 367–376. DOI 10.5890/DNC.2020.09.003.
- 37. Moatimid, G. M., El-Dib, Y. O., Zekry, M. H. (2020). The nonlinear instability of a cylindrical interface between two hydromagnetic darcian flows. Arabian Journal for Science and Engineering, 45(1), 391–409. DOI 10.1007/s13369-019-04192-z.
- 38. El-Dib, Y. O. (2020). Modified multiple scale technique for the stability of the fractional delayed nonlinear oscillator. Pramana–Journal of Physics, 94(1), 56. DOI 10.1007/s12043-020-1930-0.
- 39. Shen, Y., El-Dib, Y. O. (2020). A periodic solution of the fractional sine-Gordon equation arising in architectural engineering. Journal of Low Frequency Noise, Vibration and Active Control, 40(2), 683–691. DOI 10.1177/1461348420917565.
- 40. Elgazery, N. S. (2020). A periodic solution of the newell-whitehead-segel (NWS) wave equation via fractional calculus. Journal of Applied and Computational Mechanics, 6, 1293–1300. DOI 10.22055/jacm.2020.33778.2285.
- 41. Yu, D. N., He, J. H., Garcıa, A. G. (2019). Homotopy perturbation method with an auxiliary parameter for nonlinear oscillators. Journal of Low Frequency Noise, Vibration and Active Control, 38(3–4), 1540–1554. DOI 10.1177/1461348418811028.
- 42. He, J. H., El-Dib, Y. O. (2020). Homotopy perturbation method for Fangzhu oscillator. Journal of Mathematical Chemistry, 58(10), 2245–2253. DOI 10.1007/s10910-020-01167-6.
- 43. He, J. H., Jin, X. (2020). A short review on analytical methods for the capillary oscillator in a nanoscale deformable tube. Mathematical Methods in the Applied Sciences, 38(3), 1676. DOI 10.1002/mma.6321.
- 44. He, C. H., He, J. H., Sedighi, H. M. (2020). Fangzhu (方诸An ancient Chinese nanotechnology for water collection from air: History, mathematical insight, promises and challenges. Mathematical Methods in the Applied Sciences,2020, 1–10. DOI 10.1002/mma.6384.
- 45. He, J. H. (2020). A short review on analytical methods for to a fully fourth-order nonlinear integral boundary value problem with fractal derivatives. International Journal of Numerical Methods for Heat and Fluid Flow, 30(11), 4933–4943. DOI 10.1108/HFF-01-2020-0060.
- 46. Nash, C., Sen, S. (1983). Topology and geometry for physicists. London, UK: Academic Press.
- 47. El-Dib, Y. O. (2001). Nonlinear Mathieu equation and coupled resonance mechanism. Chaos, Solitons & Fractals, 12(4), 705–720. DOI 10.1016/S0960-0779(00)00011-4.
- 48. El-Dib, Y. O. (2016). Stability criterion for time-delay 3-dimension damped Mathieu equation. Science and Engineering Applications, 1(5), 76–88.

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. |