Study of nonlinear dynamics of a particle on a rotating parabola by using efficient mathematical techniques
A. Nikkar^{1} , H. Jodeiri^{2} , Sh. Ahari Abesh Ahmad^{3}
^{1}Young Researchers and Elite Club, Ahar Branch, Islamic Azad University, Ahar, Iran
^{2}Young Researchers and Elite Club, Tabriz Branch, Islamic Azad University, Tabriz, Iran
^{3}Faculty of Civil Engineering, University of Tabriz, Tabriz, Iran
^{2}Corresponding author
Journal of Vibroengineering, Vol. 17, Issue 7, 2015, p. 3967-3973.
Received 13 March 2015; received in revised form 8 July 2015; accepted 15 July 2015; published 15 November 2015
JVE Conferences
In this paper, nonlinear dynamics of a particle on a rotating parabola arising in the mechanical engineering is proposed. Two efficient mathematical techniques called variational iteration method and variational approach are used to solve the governing differential equation of motion. To assess the accuracy of solutions, we compare the results with the Runge-Kutta 4th order. An excellent agreement of solutions with the numerical results and published results has been demonstrated. The results reveal that the present methods are very effective and convenient in predicting the solution of such problems and find a wide application in new engineering problems.
Keywords: variational iteration method, variational approach, nonlinear problems, rotating parabola.
1. Introduction
Many phenomena exist in the environment where a particle's motion can be observed on them, such as centrifugation, centrifugal filters, industrial hopper, etc. Motion surface has different shapes especially for rotating application as it can be circular, parabolic or conical. On the other hand, the governing equations in most phenomena in applied sciences and engineering are nonlinear. In general, finding an exact or close-form solution for nonlinear problems is very difficult. So, with the rapid development of nonlinear science, there appears an ever-increasing interest of scientists in the analytical asymptotic techniques for nonlinear problems, and several analytical approximation methods have been developed; such as: homotopy perturbation method [1], variational iteration method [2], max-min approach [3], variational approach [4], Hamiltonian approach [5] and Laplace iteration method [6] .In recent years, these methods successfully applied arising in the applied sciences and engineering. For example; the vibrational behavior of quintic nonlinear in extensional beam on two-parameter elastic substrate based on the three mode assumptions is investigated by Sedighi [7]. He employed parameter expansion method to obtain the approximate expressions of nonlinear frequency-amplitude relationship for the first, second, and third modes of vibrations. A novel approach is applied to the analytic treatment of nonlinear fifth-order Equations by Saravi and Nikkar [8]. Ghaffarzadeh and Nikkar [9] applied a new analytical method called the variational iteration method-II (VIM-II) for the differential equation of the large deformation of a cantilever beam under point load at the free tip. Askari et al. [10] applied He’s energy balance method and He’s variational approach to frequency analysis of nonlinear oscillators with rational restoring force. Sedighi and Daneshmand [11] studied nonlinear transversely vibrating beams by the homotopy perturbation method with an auxiliary term. Barari et al. [12] studied non-linear vibration behavior of geometrically non-linear Euler-Bernoulli beams using variational iteration method and parameter perturbation method. Sedighi and Shirazi [13] presented a new asymptotic procedure to predict the nonlinear vibrational behaviour of classical micro-beams pre-deformed by an electric field using Parameter expansion method. Bagheri et al. [14] studied the nonlinear responses of clamped-clamped buckled beam. They used two efficient mathematical techniques called He’s variational approach and Laplace iteration method in order to obtain the responses of the beam vibrations. Salehi et al. [15] applied two efficient methods to consider large deformation of cantilever beams under point load. Younesian et al. [16] studied free oscillations of beams on nonlinear elastic foundations by VIM. Askari et al. [17] applied higher order Hamilton approach to nonlinear vibrating systems, and many other problems solved by these methods [18-21]. Recently, nonlinear dynamics of a particle on a rotating parabola are analyzed by means of the analytic and semi-analytic approaches by Mirzabeigy et al. [22].
The Variational Approach which is also called He’s variational approach, as well as the variational iteration method, was utilized here to obtain the analytical expression for the following model of nonlinear dynamics of a particle on a rotating parabola which introduced by Nayfeh and Mook [23]:
For validation of methods, we examine the advantages and disadvantages of applied methods by comparison of results with fourth-order Runge-Kutta solutions.
2. Basic concept of He’s variational approach
The variational approach to nonlinear oscillators was first proposed by Chinese mathemacian, Ji-Huan He [4]. We give a brief introduction of the method. To illustrate the method precisely, consider the following generalized nonlinear oscillations:
where $f\left(u\right)$ is a nonlinear function of $u$, $u\text{'}$ and $u\text{'}\text{'}$.
For simplicity, if function $f$ depends on $u$ only, its variational functional can be obtained as:
where $F$ is the potential, $dF/du=f$ and $T$ is the period of the nonlinear oscillator.
We assume that its approximate solution can be expressed as:
where $A$ and $\omega $ are the amplitude and frequency of the oscillator, respectively. Substituting Eq. (4) into Eq. (3), and setting $dJ/dA=0$, we can obtain an inexplicit amplitude-frequency relationship of Eq. (1).
3. Implementation of first order He’s variational approach
In this section, we solve Eq. (1) via variational approach. We can, now, easily obtain the following variational formulation:
Assume that its approximate solution can be expressed as:
where $\omega $ is the frequency to be determined and $A$ is the amplitude of oscillation. Substituting Eq. (6) into Eq. (5) results in:
Making $J$ stationary with respect to $A$, according to He’s method, we obtain:
Finally collocation at $\omega t=\pi /4$, the frequency can be approximated as:
and zero-order approximate solution from Eq. (22) is:
4. Basic concept of He’s variational iteration method
To clarify the basic ideas of He’s variational iteration method consider the following general differential equation:
where $L$ is a linear and $N$ is a nonlinear operator, $f\left(x\right)$ an inhomogeneous or forcing term. According to the variational iteration method, we can construct a correct functional as follows:
where $\lambda $ is a general Lagrange multiplier, which can be identified optimally via the variational theory, the subscript $n$ denotes the $n$th approximation, ${\stackrel{~}{u}}_{n}$ is considered as a restricted variation, i.e. $\delta {\stackrel{~}{u}}_{n}=$ 0.
For linear problems, its exact solution can be obtained by only one iteration step due to the fact that the Lagrange multiplier can be exactly identified. In this method, the problems are initially approximated with possible unknowns and it can be applied in non-linear problems without linearization or small parameters. The approximate solutions obtained by the proposed method rapidly converge to the exact solution.
5. Implementation of He’s variational iteration method
Now, for solving Eq. (1) by using VIM, first, we assume that the angular frequency of the system (1) is $\omega $, we have the following linearized equation:
So we can rewrite Eq. (1) in the form:
where $f\left(u\right)=({\epsilon}^{2}-{\omega}^{2})u+4{q}^{2}{u}^{2}\ddot{u}+4{q}^{2}u{\dot{u}}^{2}=0$.
Applying the proposed method, the following iterative formula is formed as:
Its stationary conditions can be obtained as follows:
${\left.\lambda \right|}_{\tau =t}=0,$
$\lambda \text{'}\text{'}-{\omega}^{2}{\left.\lambda \right|}_{\tau =t}=0.$
The Lagrange multiplier, therefore, can be identified as:
Substituting the identified multiplier into Eq. (15) results in the following iteration formula:
$+\frac{1}{\omega}\underset{0}{\overset{t}{\int}}\mathrm{s}\mathrm{i}\mathrm{n}\omega (\tau -t)\left({\ddot{u}}_{n}\right(\tau )+{\epsilon}^{2}{u}_{n}(\tau )+4{q}^{2}{u}_{n}(\tau {)}^{2}{\ddot{u}}_{n}\left(\tau \right)+4{q}^{2}{u}_{n}\left(\tau \right){\dot{u}}_{n}(\tau {)}^{2})d\tau .$
Assuming its initial approximate solution has the form:
And substituting Eq. (19) into Eq. (1) leads to the following residual:
$-4{q}^{2}{A}^{3}{\omega}^{2}\mathrm{c}\mathrm{o}{\mathrm{s}}^{3}\left(\omega t\right)+4{q}^{2}{A}^{2}{\omega}^{2}\mathrm{c}\mathrm{o}{\mathrm{s}}^{2}\left(\omega t\right).$
By the Eq. (18), we can obtain:
$-\left.2{\epsilon}^{2}\mathrm{s}\mathrm{i}\mathrm{n}\left(\omega t\right)t+4{A}^{2}{q}^{2}{\omega}^{2}\mathrm{s}\mathrm{i}\mathrm{n}\left(\omega t\right)t-{A}^{2}{q}^{2}\omega \mathrm{c}\mathrm{o}\mathrm{s}\left(3\omega t\right))\right).$
In the same manner, the rest of the components of the iteration formula can be obtained. In order to ensure that no secular terms appear in ${u}_{1}$, resonance must be avoided. To do so, the coefficient of $\mathrm{c}\mathrm{o}\mathrm{s}\left(\omega t\right)$ in Eq. (20) requires being zero, i.e.:
6. Results and discussion
In this stage, we show the some figures to validation of the methods. In order to verify the precision of the presented methods, current results were compared with Runge-Kutta 4th order method, presented in Figs. 1, 2. It is observable that our results are in excellent agreement with the results provided by Runge-Kutta 4th order. Also, for show the effects of different parameters on the period, we presented Figs. 3, 4. We can see that period increases when value of the $\epsilon $ decreases or value of the $q$ increases and other parameters remain constant.
Fig. 1. Comparison between the VIM, VA with the Runge-Kutta 4th order when $\epsilon =$1, $q=A=$ 1
Fig. 2. Comparison between the VIM, VA with the Runge-Kutta 4th order when $\epsilon =$1, $q=A=$ 0.5
Fig. 3. The effects of $q$ on the period of motion when $\epsilon =A=$ 1
Fig. 4. The effects of ɛ on the period of motion when $q=$ 2 and $A=$1
7. Conclusion
Study of nonlinear dynamics of a particle on a rotating parabola is main objective of this paper. We successfully applied two efficient and easy mathematical techniques called He’s variational iteration method and He’s variational approach to indicate the frequency-amplitude relationships. Comparison of results with numerical solutions are investigated. The performance and capability of applied methods are discussed by showing some figures.
References
- Omidvar M., Barari A., Momeni M., Ganji D. D. New class of solutions for water infiltration problems in unsaturated soils. Geomechanics and Geoengineering: An International Journal, Vol. 5, Issue 2, 2010, p. 127-135. [Search CrossRef]
- Salehi P., Yaghoobi H., Torabi M. Application of the differential transformation method and variational iteration method to large deformation of cantilever beams under point load. Journal of Mechanical Science and Technology, Vol. 26, Issue 9, 2012, p. 2879-2887. [Search CrossRef]
- He J. H. Some asymptotic methods for strongly nonlinear equations. International Journal of Modern Physics B, Vol. 20, Issue 10, 2006, p. 1141-1199. [Search CrossRef]
- He J. H. Variational approach for nonlinear oscillators. Chaos Solitons and Fractals, Vol. 34, 2007, p. 1430-1439. [Search CrossRef]
- He J. H. Hamiltonian approach to nonlinear oscillators. Physics Letters A, Vol. 374, 2010, p. 2312-2314. [Search CrossRef]
- Hashemi M. H., Nikkar A. An efficient iterative method for solving the governing equation of water infiltration problems in unsaturated soils. Geomechanics and Geoengineering: An International Journal, Vol. 9, Issue 4, 2014. [Search CrossRef]
- Sedighi Hamid M. Nonlinear free vibrations of quintic inextensional beams lying on Winkler elastic substrate based on three-mode assumptions. Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-body Dynamics, Vol. 228, Issue 2, 2014, p. 213-225. [Search CrossRef]
- Saravi M., Nikkar A. Promising technique for analytic treatment of nonlinear fifth-order equations. World Journal of Modelling and Simulation, Vol. 10, 2014, p. 27-33. [Search CrossRef]
- Ghaffarzadeh H., Nikkar A. Explicit solution to the large deformation of a cantilever beam under point load at the free tip using the variational iteration method-II. Journal of Mechanical Science and Technology, Vol. 27, Issue 11, 2013, p. 3433-3438. [Search CrossRef]
- Askari H., Yazdi M. K., Saadatnia Z. Frequency analysis of nonlinear oscillators with rational restoring force via He’s energy balance method and He’s variational approach. Nonlinear Science Letters A, Vol. 1, Issue 4, 2010, p. 425-430. [Search CrossRef]
- Sedighi H. M., Daneshmand F. Nonlinear transversely vibrating beams by the homotopy perturbation method with an auxiliary term. Journal of Applied and Computational Mechanics, Vol. 1, Issue 1, 2014, p. 1-9. [Search CrossRef]
- Barari A., Kaliji H. D., Ghadimi M., Domairry G. Non-linear vibration of Euler-Bernoulli beams. Latin American Journal of Solids and Structures, Vol. 8, Issue 2, 2011, p. 139-148. [Search CrossRef]
- Sedighi H. M., Shirazi K. H. Vibrations of micro-beams actuated by an electric field via parameter expansion method. Acta Astronautica, Vol. 85, 2013, p. 19-24. [Search CrossRef]
- Bagheri S., Nikkar A., Ghaffarzadeh H. Study of nonlinear vibration of Euler-Bernoulli beams by using analytical approximate techniques. Latin American Journal of Solids and Structures, Vol. 11, 2014, p. 157-168. [Search CrossRef]
- Salehi P., Yaghoobi H., Torabi M. Application of the differential transformation method and variational iteration method to large deformation of cantilever beams under point load. Journal of mechanical science and technology, Vol. 26, Issue 9, 2012, p. 2879-2887. [Search CrossRef]
- Younesian D., Saadatnia Z., Askari H. Analytical solutions for free oscillations of beams on nonlinear elastic foundations using the variational iteration method. Journal of Theoretical and Applied Mechanics, Vol. 50, Issue 2, 2012, p. 639-652. [Search CrossRef]
- Askari H., Nia Z. S., Yildirim A., Yazdi M. K., Khan Y. Application of higher order Hamilton approach to nonlinear vibrating systems. Journal of Theoretical and Applied Mechanics, Vol. 51, Issue 2, 2013, p. 287-296 [Search CrossRef]
- Sedighi H. M., Daneshmand F., Abadyan M. Modified model for instability analysis of symmetric FGM double-sided nano-bridge: corrections due to surface layer, finite conductivity and size effect. Composite Structures, Vol. 132, 2015, p. 545-557. [Search CrossRef]
- Khan Y., Taghipour R., Falahian M., Nikkar A. A new approach to modified regularized long wave equation. Neural Computing and Applications, Vol. 23, 2013, p. 1335-1341. [Search CrossRef]
- Nikkar A., Bagheri S., Saravi M. Study of nonlinear vibration of Euler-Bernoulli beams by using analytical approximate techniques. Latin American Journal of Solids and Structures, Vol. 11, 2014, p. 320-329. [Search CrossRef]
- Taghipor R., Akhlaghi T., Nikkar A. Explicit solution of the large amplitude transverse vibrations of a flexible string under constant tension. Latin American Journal of Solids and Structures, Vol. 11, 2014, p. 545-555. [Search CrossRef]
- Mirzabeigy A., Yazdi M. K., Yildirim A. Nonlinear dynamics of a particle on a rotating parabola via the analytic and semi-analytic approaches. Journal of the Association of Arab Universities for Basic and Applied Sciences, Vol. 13, 2013, p. 38-43. [Search CrossRef]
- Nayfeh A. H., Mook D. T. Nonlinear Oscillation. Willey, 1979. [Search CrossRef]