Homotopy Perturbation Method for the Fractal Toda Oscillator

: The fractal Toda oscillator with an exponentially nonlinear term is extremely difficult to solve; Elias ‐ Zuniga et al. (2020) suggested the equivalent power ‐ form method. In this paper, first, the fractal variational theory is used to show the basic property of the fractal oscillator, and a new form of the Toda oscillator is obtained free of the exponential nonlinear term, which is similar to the form of the Jerk oscillator. The homotopy perturbation method is used to solve the fractal Toda oscillator, and the analytical solution is examined using the numerical solution which shows excel ‐ lent agreement. Furthermore, the effect of the order of the fractal derivative on the vibration prop ‐ erty is elucidated graphically.


Introduction
An oscillation occurs when its kinetic energy and its potential energy are changed alternatively, while the total energy remains unchanged. Its variational formulation can be expressed as [1][2][3]: where 2 1 2 u  is kinetic energy and   u f is the potential energy.
They meet the following energy balance equation [4]: where H is the Hamiltonian constant depending on the initial conditions. Generally, a nonlinear oscillator can be written in the form: with the initial conditions: , we have the well-known Duffing oscillator, which reads as follows: The potential function,   u f , will greatly affect its periodic property. In this paper, we consider a potential with an exponential form: ( ) u f u e u   and its variational formulation is: The nonlinear oscillator can be obtained from Equation (8) as a stationary condition: This is the Toda oscillator [7-10].

Fractal Toda Oscillator and Fractal Weierstrass Theorem
The Toda oscillator is well-known for its chaos and bifurcation properties [9,10], and its fractal partner was first studied in [8], where some new properties were revealed. The fractal vibration theory has attracted much attention recently and can solve what the classic vibration theory cannot solve. For example, the pull-in instability of a micro-electromechanical system (MEMS) oscillator can never be solved by differential models, and it was Tian and colleagues that made a momentous contribution to the landmark using a fractal modification [11][12][13]. Ji et al. suggested a new fractal vibration theory to avoid vibration damage [14]. A fractal oscillation revealed, for the first time, the mechanism of Fangzhu device for water collection from the air [15][16][17][18]. In a microgravity space, a fractal modification of Newton's law is a must to describe the motion [19]. Zuo established a fractal vibration model for the 3D printing process [20]. Elias-Zuniga et al. suggested a fractal model for current generation in porous electrodes [21]. Elias-Zuniga et al. extended the Duffing equation to its fractal modification [22]. He, Liu, and Gepreel found the large period property of a fractal concrete beam [23]. He et al. revealed the thermal response of a fractal concrete [24]. Various analytical methods for the frequency-amplitude relationship of a fractal oscillator have appeared in literature, among which He's frequency formulation was the simplest [25][26][27][28].
In a fractal space, the variational formulation of Equation (8) can be modified as [29]: where du dt  is the two-scale fractal derivative [30]. The fractal variational theory is helpful to establish a governing equation in a fractal space. Some useful fractal variational principles can be found in the literature, for example, the fractal variational theory for compressible fluids [31], the fractal variational principle for the telegraph equation [32], the fractal variational theory for Chaplygin-He gas [33], and the fractal variational principle for a solitary wave traveling [34][35][36][37][38].
Consider a fractal variational principle in the form: and introduce a new variable w defined as: du w dt   (12) The fractal variational formulation given in Equation (11) becomes: which is subject to the constraint of Equation (12). The fractal Weierstrass function is defined as [29]: According to the fractal Weierstrass theorem [29], Equation (11) is a strong minimal variational principle if the following conditions are satisfied: In our study for Equation (10), the fractal Weierstrass function is: It is obvious that Therefore, Equation (10) is a strong minimal variational principle. The fractal Toda oscillator can be obtained from minimizing Equation (10), which is: The initial conditions are: Using the two-scale transform [30]: The fractal Toda oscillator can be expressed in the form: where the dot is the derivative T concerning.

A Simplified Model for the Fractal Toda Oscillator
To establish a compatible form of the Toda oscillator, it is suitable to convert it to a Jerk oscillator. To accomplish this aim, an annihilation operator is required for the function u e which is found to be  , therefore, Equation (21) can be converted to the following form: This is an alternative equation free of the exponential nonlinear term. It is known as a Jerk oscillator. It lacks the linear term; therefore, it can be modified by adding the term u 3  to both sides as follows: where 3  is an unknown to be determined later; therefore, these nonlinear equations can be solved by the HPM. In [39], a reduced-order model for analysis is proposed by using the HPM. Equation (23) can be arranged in the form: Since the original equation is a second-order equation, then, it is suitable to reduce the rank of Equation (23). This can be accomplished when operating on both sides of Equation (24) by the inverse of     D as follows [40][41][42][43][44][45][46][47][48]: This is a second-order differential equation that contains quadratic nonlinear terms.

Application of the HPM
The homotopy perturbation method was first proposed to solve nonlinear differential equations [40], and it is also an effective tool for various nonlinear oscillators [41][42][43]. To solve Equation (21) by the homotopy perturbation method, we can establish the homotopy equation in the form: Accordingly, we expand the solution in the form: ..., are unknown constants to be later determined. Because there is only one unknown in the homotopy Equation (16), that is  , the first-order perturbation is enough to determine it. Substituting the expansion (27) into the homotopy Equation (26) yields: Equation (28) Removing the terms that producing secular terms gives: Accordingly, the periodic solution of Equation (30) has the form: Employing (32) and (33) The first-order approximate solution can be performed by employing (33) and (34) into the expansion (27)

Numerical Illustration
For more opportunities, the attained approximate analytical solution must validate the numerical solution. For this perseverance, the numerical solution of Equation (21) was utilized to confirm the approximate solution (35). Therefore, Figure 1 is designed to demonstrate the bounded analytical solution given from the HPM, together with the numerical solution that is given by Mathematica software. Through this figure, the usage of the sample chosen system exposed that the value of the amplitude 1 .

Conclusions
In this paper, we show the basic properties of a fractal oscillator in the frame of fractal variational theory. The fractal Tota oscillator can be obtained by minimizing the fractal variational formulation. The complexity of the Toda equation was relaxed by converting it to the jerk oscillator. The homotopy perturbation method was adopted to solve the nonlinear third-order differential equation, and the results are good enough for practical applications. The accuracy of the obtained solution was examined by comparing it with the numerical solutions which show excellent agreement even for relatively large amplitudes.  Data Availability Statement:In this study there is no data used.

Conflicts of Interest:
The author declares that there are no competing interests regarding the publication of the present paper