Application of the Homotopy Perturbation Method to Solve Nonlinear Equations Arising in Oscillatory Systems

Subject Areas : Mechanical Engineering

Mohammadjavad Mahmoodabadi ¹ , Neda Amiri ²

1 -
2 -

Received: 2024-01-08 Accepted : 2024-07-30 Published : 2025-01-19

Keywords: Duffing-Holmes Model, Homotopy Perturbation Method, Nonlinear Equations, Oscillatory Systems ,

Abstract :

In this research, the application of the homotopy perturbation method to solve nonlinear Equations arising in oscillatory systems is investigated. In this way, the performance of the Homotopy Perturbation Method (HPM) is compared with the numerical methods to find the solutions of nonlinear Equations in the vibration field. To this end, the Duffing–Holmes oscillatory model with nonlinear terms is regarded and solved by the HPM method. In order to validate the obtained solution by the HPM, the answers are compared with those of numerical methods. The results clearly depict that the homotopy perturbation method, without needing to small parameters, could present the answers near to the exact solutions and also to the numerical one.

References:

[1] Kumar, Y., The Rayleigh–Ritz Method for Linear Dynamic, Static and Buckling Behavior of Beams, Shells and Plates: A Literature Review, Journal of Vibration and Control, Vol. 24, No. 7, 2017.
[2] Pradhan, S., Chakraverty, K. K., Free Vibration of Euler and Timoshenko Functionally Graded Beams by Rayleigh–Ritz Method, Composites Part B: Engineering, Vol. 51, 2013, pp. 175-184.
[3] Yserentant, H., A Short Theory of the Rayleigh–Ritz Method, Computational Methods in Applied Mathematics, Vol. 13, No. 4, 2013, pp. 495–502.
[4] Bhat, R.B., Natural Frequencies of Rectangular Plates Using Characteristic Orthogonal Polynomials in Rayleigh-Ritz Method, Journal of Sound and Vibration, Vol. 102, No. 4, 1985, pp. 493-499.
[5] Lu, Y. Y., Belytschko, and T., Gu, L., A New Implementation of The Element Free Galerkin Method, Computer Methods in Applied Mechanics and Engineering, Vol. 113, No. (3–4), 1994, pp. 397-414.
[6] Thomas, J. R., Engel, H. G., Mats, L. M., and Larson, G., The Continuous Galerkin Method Is Locally Conservative, Journal of Computational Physics, Vol. 163, No. 2, 2000, pp. 467-488.
[7] Demkowicz, L. F., Gopalakrishnan, J., An Overview of the Discontinuous Petrov Galerkin Method, Recent Developments in Discontinuous Galerkin Finite Element Methods for Partial Differential Equations, Vol. 157, 2013, pp. 149–180.
[8] Thomas, J. R., Scovazzi Pavel, H. G., Bochev, B., and Buffa, A., A Multiscale Discontinuous Galerkin Method with The Computational Structure of a Continuous Galerkin Method, Computer Methods in Applied Mechanics and Engineering, Vol. 195, No. (19–22), 2006, pp. 2761-2787.
[9] Adomian, G., Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Academic Publishers, 1994.
[10] Wazwaz, A. M., A Reliable Modification of Adomian Decomposition Method, Applied Mathematics and Computation, Vol. 102, No. 1, 1999, pp. 77-86.
[11] Hosseini, M. M., Nasabzade, H., On the Convergence of Adomian Decomposition Method, Applied Mathematics and Computation, Vol. 182, No. 1, 2006, pp. 536-543.
[12] Allahviranloo, T., The Adomian Decomposition Method for Fuzzy System of Linear Equations, Applied Mathematics and Computation, Vol. 163, No. 2, 2005, pp. 553-563.
[13] Tatari, M., Dehghan, M., and Razzaghi, M., Application of the Adomian Decomposition Method for the Fokker–Planck Equation, Mathematical and Computer Modelling, Vol. 45, No. (5–6), 2007, pp. 639-650.
[14] Ennis, J. E., The Kantorovich and Contractive Mapping Theorems, Numerical Methods for Unconstrained Optimization and Nonlinear Equations, 1983, pp. 92–94.
[15] Donea, J., A Taylor–Galerkin Method for Convective Transport Problems, International Journal for Numerical Methods in Fluids, Vol. 20, No. 1, 1984, pp. 101-119.
[16] Timmermans, L. J. P., Van De Vosse, F. N., and Minev, P. D., Taylor-Galerkin-Based Spectral Element Methods for Convection-Diffusion Problems, Namerical Methods in Fluids, Vol. 18, No. 9, 1994, pp. 853-870.
[17] Shafiee Sarvestany, A. R., Mahmoodabadi M. J., FA-ABC: A Novel Combination of Firefly Optimization Algorithm and Artificial Bee Colony for Mathematical Test Functions and Real-World Problems, Advanced Design and Manufacturing Technology, Vol. 15, No. 2, 2022, pp. 69-82.
[18] Mahmoodabadi, M. J., Nemati, A. R., A New Optimum Numerical Method for Analysis of Nonlinear Conductive Heat Transfer Problems, Journal of the Brazilian Society of Mechanical Sciences and Engineering, Vol. 43, 2021, pp. 1-8.
[19] Mahmoodabadi, M. J., Sadeghi Googhari F., Numerical Solution of Time-Dependent Schrodinger Equation by Combination of The Finite Difference Method and Particle Swarm Optimization, Journal of Research on Many-body Systems, Vol. 11, No. 1, 2021, pp. 114-127.

Full-Text:



Int.  J.   Advanced Design and Manufacturing Technology, 2024, Vol. 17, No. 4, pp. 29-33

DOI: 10.71644/admt.2024.873801                                       ISSN: 2252-0406 https://admt.isfahan.iau.ir

 

Application of the Homotopy Perturbation Method to Solve Nonlinear Equations Arising in Oscillatory Systems

Mohammad Javad Mahmoodabadi *, Neda Amiri 

Department of Mechanical Engineering,

Sirjan University of Technology, Sirjan, Iran

E-mail: [e-mail], [e-mail]

*Corresponding author

Received: 8 January 2024, Revised: 14 April 2024, Accepted: 30 July 2024

Abstract: In this research, the application of the homotopy perturbation method to solve nonlinear Equations arising in oscillatory systems is investigated. In this way, the performance of the Homotopy Perturbation Method (HPM) is compared with the numerical methods to find the solutions of nonlinear Equations in the vibration field. To this end, the Duffing–Holmes oscillatory model with nonlinear terms is regarded and solved by the HPM method. In order to validate the obtained solution by the HPM, the answers are compared with those of numerical methods. The results clearly depict that the homotopy perturbation method, without needing to small parameters, could present the answers near to the exact solutions and also to the numerical one.

Keywords: Duffing-Holmes Model, Homotopy Perturbation Method, Nonlinear Equations, Oscillatory Systems

 

Biographical notes: Mohammad Javad Mahmoodabadi received his BSc and MSc degrees in Mechanical Engineering from Shahid Bahonar University of Kerman, Iran in 2005 and 2007, respectively. He received his PhD degree in Mechanical Engineering from the University of Guilan, Rasht, Iran in 2012. During his research, he was a scholar visitor with Robotics and Mechatronics Group, University of Twente, Enchede, Netherlands for 6 months. Now, he is an Associate Professor of Mechanical Engineering at the Sirjan University of Technology, Sirjan, Iran. His research interests include optimization algorithms, nonlinear and robust control, robotics, and computational methods. Neda Amiri received her BSc in robotic engineering from the Sirjan University of Technology, Sirjan, Iran in 2020. Her current research focuses on robotic, and homotopy perturbation methods.

 

 

 

 

 

1        Introduction

There are many engineering problems, such as those appearing in mechanical vibrations which are nonlinear, therefore most of them are solved by using numerical approaches, while the others are solved via the analytical methods. In the numerical methods, stability and convergence should be considered to avoid divergent or inappropriate results. On the other hand, many different analytical methods have been recently introduced to eliminate these issues. In this way, Kumar [1] presented a comprehensive literature review on the application of the Rayleigh-Ritz method to analyze vibration, static, and buckling characteristics of beams, shells, and plates by employing different theories. The governing Equations were obtained using the Rayleigh-Ritz method to study the effects of constituent volume fractions, slenderness ratios, and the beam theories on the natural frequencies by Pradhan and Chakraverty [2]. Yserentant [3] displayed some new error estimates for the eigenvalues and Eigen functions obtained by the Rayleigh-Ritz method, the common variational method to solve Eigen problems. Natural frequencies of rectangular plates were obtained by employing a set of beam characteristic orthogonal polynomials in the Rayleigh-Ritz method by Bhat [4]. 

A new implementation was developed by Lu et al. [5] based on a modified variational principle in which the Lagrange multipliers were replaced at the outset by their physical meaning so that the discrete Equations were banded. Thomas et al. [6] examined the conservation law structure of the continuous Galerkin method for solving the scalar, advection-diffusion Equation as a model problem. Demkowicz and Gopalakrishnan [7] discussed the principles and methodology of the discontinuous Petrov Galerkin method with optimal test functions and provided a literature review on the subject. Thomas et al. [8] developed a computational formulation that combines the advantages of discontinuous Galerkin methods with the data structure of their continuous Galerkin counterparts. Wazwaz [10] proposed a powerful modification of the Adomian decomposition method, introduced in the 1970s to the 1990s by George Adomian [9], to accelerate the rapid convergence of the series solution.  Moreover, a simple method to determine the rate of convergence of the Adomian decomposition method was introduced by Hosseini et al. [11]. In this way, the application of the Adomian method for solving fuzzy systems of linear Equations was considered by Allahviranloo [12].

The solution of an initial value problem of the parabolic type was discussed by Tatari et al. [13] to propose an alternative method of solution, one not based on finite difference or finite element or spectral methods. The Kantorovich theorem, or Newton–Kantorovich theorem, as a mathematical statement on the semi-local convergence of Newton's method, was first stated by Leonid Kantorovich to form the Banach fixed-point theorem [14]. A Taylor–Galerkin method was described to derive finite element schemes for the scalar convection Equation in one or more space dimensions based on the forward-time Taylor series expansions by Donea [15]. Several explicit Taylor-Galerkin-based time integration schemes were proposed for the solution of both linear and non-linear convection problems with the divergence-free velocity by Timmermans et al. [16]. Shafiee Sarvestany and Mahmoodabadi [17] investigated a novel combination of the firefly optimization algorithm and artificial bee colony for mathematical test functions and real-world problems. Mahmoodabadi, and Nemati [18] presented an optimum numerical method for analysis of nonlinear conductive heat transfer problems. Mahmoodabadi and Sadeghi Googhari [19] studied numerical solutions of the time-dependent Schrodinger Equation by the combination of the finite difference method and particle swarm optimization.

In this research work, the basic idea of the HPM is introduced, its application on the oscillatory Equations is studied, and a comparison with the exact solution is also made.

2        Homotopy Perturbation Method (HPM)

To illustrate the basic ideas of the HPM, the following nonlinear differential Equation is considered.

 

                                  (1)

 

With boundary conditions:

 

                                           (2)

 

Where  denotes a general differential operator,  represents a boundary operator,  signifies a known analytical function, and  is the boundary of domain . 

Operator  can be generally divided into two linear  and nonlinear  parts. Therefore, “Eq. (1)” can be rewritten as follows:

 

                                        (3)

 

Hence, Homotopy function is constructed as follows:

 

                                                                                       (4)

 

Where  denoted the homotopy parameter.

According to the homotopy perturbation idea, the approximate solution of “Eq. (4)” can be expressed as a series of the powers of , i.e:

 

                                       (5)

3        APPLICATIONS of the HPM for Duffing Holmes Oscillator

Consider the Duffing-Holmes oscillator Equation with the following state-space configuration:

 

                                                                                (6)

 

                                       (7)

 

By regarding: 

 

 

 

And

 

 

 

Then, 

 

                                                        (8)

 

Therefore, the homotopy function is defined as follows:

 

                                                                                             (9)

 

By employing  and boundary conditions  and  the following relations are obtained:

 

 

3] = 0; 

                                                     (10)

 

By arranging the terms of the above Equation with respect to the power of parameter p, the following differential Equations are reached:

 

 

                                         (11)

 

3        

                                     (12)

 

                    (13)

 

By utilizing the Laplace transformation on Equation (11), we have:

 

                                                 (14)

 

Then,

 

                                                 (15)

 

Hence,

 

                                              (16)

 

Therefore:

 

 

                                             (17)

 

By employing the initial conditions, parameters  and  would be computed as follows:

 

                                             (18)

 

By substituting  and  for  and  into Equation (12), the first term of the solution could be introduced as follows:

 

                            (19)

 

If Equation (12) is rewritten as follows:

 

3                                     (20)

 

If

 

3 =   (21)

 

Then, particular solution z is formulated as follows:

 

     (22)

 

With

 

                            (23)

 

 

 

And

 

                                (24)

 

Therefore,

 

                (25)

 

The following algebraic Equations will determine unknown parameters a, b, c, and d.

 

 

 

 

               (26)

 

Finally, the particular solution is rewritten as follows.

 

            (27)

 

By applying the general solution, the total relation for  could be introduced as follows:

 

                                             (28)

Where, constant parameters  and  would be calculated by the following initial conditions:

 

                                                                  (29)

 

                                            (30)

 

Regarding the homotopy solution, the answer to the problem would be defined as follows:

 

                                   (31)

 

Finally,

 

                                            (32)

 

In order to validate the obtained solution by the HPM, the answers related to interval time [0, 10] (s) are compared with those of the fourth-order Runge-Kutta Method (RKM) in “Fig. 1”. Although the results have a good agreement at the initial times, the differences could be obviously seen at the bigger times.

 

 

 

Fig. 1        Solutions found by the homotopy perturbation method and forth-order Runge-Kutta approach for different time intervals.

 

 

4        Conclusions

This research study implemented the homotopy perturbation method to analytically solve the nonlinear Dofing-Holmes Equation related to oscillatory dynamical systems. A closed mathematical formulation was determined to calculate the unknown parameter of the Equation at each time. The validations were performed through comparisons of the results with the numerical ones. The accuracy of the HPM was challenged by comparing the found results with those of the fourth-order Runge-Kutta technique.

References

[1]       Kumar, Y., The Rayleigh–Ritz Method for Linear Dynamic, Static and Buckling Behavior of Beams, Shells and Plates: A Literature Review, Journal of Vibration and Control, Vol. 24, No. 7, 2017.

[2]       Pradhan, S., Chakraverty, K. K., Free Vibration of Euler and Timoshenko Functionally Graded Beams by Rayleigh–Ritz Method, Composites Part B: Engineering, Vol. 51, 2013, pp. 175-184.

[3]       Yserentant, H., A Short Theory of the Rayleigh–Ritz Method, Computational Methods in Applied Mathematics, Vol. 13, No. 4, 2013, pp. 495–502.

[4]       Bhat, R.B., Natural Frequencies of Rectangular Plates Using Characteristic Orthogonal Polynomials in Rayleigh-Ritz Method, Journal of Sound and Vibration, Vol. 102, No. 4, 1985, pp. 493-499.

[5]       Lu, Y. Y., Belytschko, and T., Gu, L., A New Implementation of The Element Free Galerkin Method, Computer Methods in Applied Mechanics and Engineering, Vol. 113, No. (3–4), 1994, pp. 397-414.

[6]       Thomas, J. R., Engel, H. G., Mats, L. M., and Larson, G., The Continuous Galerkin Method Is Locally Conservative, Journal of Computational Physics, Vol. 163, No. 2, 2000, pp. 467-488.

[7]       Demkowicz, L. F., Gopalakrishnan, J., An Overview of the Discontinuous Petrov Galerkin Method, Recent Developments in Discontinuous Galerkin Finite Element Methods for Partial Differential Equations, Vol. 157, 2013, pp. 149–180.

[8]       Thomas, J. R., Scovazzi Pavel, H. G., Bochev, B., and Buffa, A., A Multiscale Discontinuous Galerkin Method with The Computational Structure of a Continuous Galerkin Method, Computer Methods in Applied Mechanics and Engineering, Vol. 195, No. (19–22), 2006, pp. 2761-2787.

[9]       Adomian, G., Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Academic Publishers, 1994.

[10]       Wazwaz, A. M., A Reliable Modification of Adomian Decomposition Method, Applied Mathematics and Computation, Vol. 102, No. 1, 1999, pp. 77-86.

[11]       Hosseini, M. M., Nasabzade, H., On the Convergence of Adomian Decomposition Method, Applied Mathematics and Computation, Vol. 182, No. 1, 2006, pp. 536-543.

[12]       Allahviranloo, T., The Adomian Decomposition Method for Fuzzy System of Linear Equations, Applied Mathematics and Computation, Vol. 163, No. 2, 2005, pp. 553-563.

[13]       Tatari, M., Dehghan, M., and Razzaghi, M., Application of the Adomian Decomposition Method for the Fokker–Planck Equation, Mathematical and Computer Modelling, Vol. 45, No. (5–6), 2007, pp. 639-650.

[14]       Ennis, J. E., The Kantorovich and Contractive Mapping Theorems, Numerical Methods for Unconstrained Optimization and Nonlinear Equations, 1983, pp. 92–94.

[15]       Donea, J., A Taylor–Galerkin Method for Convective Transport Problems, International Journal for Numerical Methods in Fluids, Vol. 20, No. 1, 1984, pp. 101-119.

[16]       Timmermans, L. J. P., Van De Vosse, F. N., and Minev, P. D., Taylor-Galerkin-Based Spectral Element Methods for Convection-Diffusion Problems, Namerical Methods in Fluids, Vol. 18, No. 9, 1994, pp. 853-870.

[17]       Shafiee Sarvestany, A. R., Mahmoodabadi M. J., FA-ABC: A Novel Combination of Firefly Optimization Algorithm and Artificial Bee Colony for Mathematical Test Functions and Real-World Problems, Advanced Design and Manufacturing Technology, Vol. 15, No. 2, 2022, pp. 69-82.

[18]       Mahmoodabadi, M. J., Nemati, A. R., A New Optimum Numerical Method for Analysis of Nonlinear Conductive Heat Transfer Problems, Journal of the Brazilian Society of Mechanical Sciences and Engineering, Vol. 43, 2021, pp. 1-8.

[19]       Mahmoodabadi, M. J., Sadeghi Googhari F., Numerical Solution of Time-Dependent Schrodinger Equation by Combination of The Finite Difference Method and Particle Swarm Optimization, Journal of Research on Many-body Systems, Vol. 11, No. 1, 2021, pp. 114-127.

 

Research paper

COPYRIGHTS

© 2024 by the authors. Licensee Islamic Azad University Isfahan Branch. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution 4.0 International (CC BY 4.0) (https://creativecommons.org/licenses/by/4.0/)

 

---