Theoretical and Computational Advances in Nonlinear Dynamical Systems
View this Special IssueResearch Article  Open Access
The Spreading Residue Harmonic Balance Method for Strongly Nonlinear Vibrations of a Restrained Cantilever Beam
Abstract
The exact solutions of the nonlinear vibration systems are extremely complicated to be received, so it is crucial to analyze their approximate solutions. This paper employs the spreading residue harmonic balance method (SRHBM) to derive analytical approximate solutions for the fifthorder nonlinear problem, which corresponds to the strongly nonlinear vibration of an elastically restrained beam with a lumped mass. When the SRHBM is used, the residual terms are added to improve the accuracy of approximate solutions. Illustrative examples are provided along with verifying the accuracy of the present method and are compared with the HAM solutions, the EBM solutions, and exact solutions in tables. At the same time, the phase diagrams and time history curves are drawn by the mathematical software. Through analysis and discussion, the results obtained here demonstrate that the SRHBM is an effective and robust technique for nonlinear dynamical systems. In addition, the SRHBM can be widely applied to a variety of nonlinear dynamic systems.
1. Introduction
A lot of problems in physical, mechanical, and aeronautical technology and even in structural applications are essentially nonlinear. Majority of the nonlinear dynamical models are mainly composed of a group of differential equations and auxiliary conditions for modeling processes [1]. In general, it is difficult to obtain the exact solution for strongly nonlinear high dimensional dynamic systems. Hence, the analytical approximate solution of the nonlinear problem has become the research object of many scholars in recent years [2–28].
Generally speaking, the fifthorder Duffing type problem with the inertial and static nonlinear terms is sophisticated all the better [29]. Recently, some scholars have tried to study this kind of nonlinear problem. For instance, Telli and Kopmaz [30] and Lai and Lim [31] used the harmonic balance method to study the linear and nonlinear springs. S.S. Chen and C.K. Chen [32] dealt with this fifthorder nonlinear problem by applying the differential transformation approach. Subsequently, Ganji et al. [33] and Mehdipour et al. [34], respectively, brought in the homotopy perturbation method, amplitudefrequency formulation, and the energy balance method. They used these methods to solve this strongly nonlinear problem, and lowerorder approximate solutions are yielded. Qian et al. [35] studied the nonlinear vibrations of cantilever beam by the HAM. Latterly, Guo et al. [36, 37] have presented the residue harmonic balance solution procedure to approximate the periodic behavior of different oscillation systems and they have obtained some more accurate results. Ju and Xue [38, 39] proposed the global residue harmonic balance method to study strongly nonlinear systems. Comparing the obtained solutions with the exact one, they discovered that the approximate results excellently agree with the exact one. Lee [40] used the multilevel residue harmonic balance method to solve a nonlinear panel coupled with extended cavity.
The principal intention of this paper is to investigate the utility of the spreading residue harmonic balance method (SRHBM) [36] for the fifthorder strongly nonlinear problem. The paper consists of the following several parts. Section 2 describes how the strongly nonlinear equation is educed from the governing equations of the cantilever beam model in a nutshell. In Section 3, the SRHBM is introduced and the solution process of different order solutions will be presented. The numerical examples of the SRHBM are rendered and compared with other solutions in Section 4. Finally, conclusion of the paper is drawn in Section 5.
2. Mathematical Formulation
An isotropic slender beam with uniform length and mass per unit length is considered, as shown in Figure 1 [25]. It is assumed that the beam thickness is much smaller than the beam length, so the effects of shear deformation and rotary inertia can be ignored. The angle of inclination is and the beam displacement is . For the boundary condition constraints, one of the conditions is hinged at the bottom of a rotational spring with stiffness , and the other condition is independent. Moreover, the intermediate lumped mass is also connected in along the beam span. By the EulerLagrange differential equation, the fifthorder Duffing type temporal problem with strongly inertial and static nonlinearities is able to be derived as follows [25]:where is the dimensionless deflection at the tip of the beam, is the maximum amplitude, the overdot indicates the derivative relative to , and , , , and are parameters. For the complete formulation of (1), readers are referred to [25] for details.
3. Solution Methodology
In the following, the spreading residue harmonic is used to solve (1). Firstly, by introducing a new variable and substituting it into (1), we can get where represents the firstorder derivative with respect to and is the unknown angular frequency of (1).
Since we discuss the existence of a periodic solution, we usually choose the base functionsand we find the expression of the steady state solutionswhere is an order parameter and are unknown.
Next, we mainly analyze the zerothorder harmonic approximation, the firstorder harmonic approximation, and the secondorder harmonic approximation.
3.1. The ZerothOrder Harmonic Approximation
To meet the initial conditions in (2), we can set the following initial guess solution of ; that is,Substituting (5) into (2), the equation is yielded:
denotes the zerothorder residual term. When , is the exact solution.
Based on the Galerkin procedure, the secular term cannot appear on the right hand side of (6). Equating the term’s coefficient to zero, we obtain a linear equation containing an unknown . Through solving that equation, we are able to work out the unknown frequency : In addition, the zerothorder approximation solution can be obtained as follows:
When the obtained zerothorder solution is substituted into (6), the terms of and generally are not zero.
3.2. The FirstOrder Harmonic Approximations
Substituting (4) into (2), the coefficient of parameter is put forward, and we can obtain
According to (3), we choose the following equation as the solution of the representation:where is the unknown.
Combining (10) with (9), we considerWe can obtainTo increase the accuracy, (11) should be added in the zerothorder residual term .
Based on the Galerkin procedure, (12) should not contain secular terms. Letting coefficients of and be zeros, hence, we obtain a linear equation set containing two unknowns and . Through solving the equation set, we can get where
Therefore, the firstorder harmonic approximation can be procured:
When the obtained firstorder solution is substituted into (12), the terms of and commonly are not zero.
3.3. The SecondOrder Harmonic Approximations
Substituting (4) to (2) and presenting the coefficient of parameter yield
Observing (16), we know it is linear with respect to and . In accordance with the form of (3), we apply where and are unknowns.
Combining (17) with (16) and calculating the equationwe can obtain
With the purpose of increasing the accuracy, (18) should be put on the firstorder residual term .
In order to prevent the right hand side of (19) from exhibiting the secular terms , , and , we make all their coefficients equal to zero. Then, we can get three linear equations containing three unknown parameters , , and . According to the three linear equations, we can solve the three unknowns. Thus, the secondorder harmonic approximation is shown:
In conclusion, we draw the thorder harmonic approximation
4. Results and Discussion
In order to make sure of the effectiveness of the current technique, we compare the results from the secondorder spreading residue harmonic balance approach with the energy balance method [34], the homotopy analysis method [35], and the exact solution [33], which are presented in Table 1, for different parameters and amplitudes of vibration , where the exact solution is computed using the numerical technique. The relative errors of vibration frequency are tabulated in Table 2.


For Mode 1 in Table 2, we observe that the relative error between and is much less than the relative error between and . The same goes for Mode 2 and Mode 4. However, the relative error between and for Mode 3 is more than the relative error between and . Hence, we conclude that the accuracy of the secondorder SRHBM solutions is improved in Mode 1, Mode 2, and Mode 4. Similarly, from Table 2, the SRHBM relative error is smaller than EBM relative error in Mode 1, Mode 2, Mode 3, and Mode 4. These results show that the approximate solutions obtained by the SRHBM are closer to the exact solutions.
To further demonstrate the accuracy of the spreading residue harmonic balance approach, the time history responses and the phase portrait are rendered for four different sets of parameters in Figures 2–5. From the phase portrait diagram, we obviously discover that the secondorder residue harmonic balance solutions are very consistent with the exact solutions. From the phase portrait, we observe that the system is a periodic motion. Moreover, in the whole range, the presented approximate solutions converge to the exact solutions. Extraordinarily, we argue for in .
5. Conclusions
In this paper, the spreading residue harmonic balance method is applied to discuss the strongly nonlinear vibration system. Particularly, we take a restrained cantilever beam as an example. The SRHBM does not need to add small parameters in the calculation process. Besides, this approach approximates the exact solution quickly and only the first or secondorder approximations. And by comparing its results with HAM and EBM for various parameters and amplitudes of vibration, it reveals that SRHBM can be used to solve a nonlinear equation with high nonlinearities. According to the figures and tables, it is effective to explain that the presented approximations are more accurate. Therefore, we can conclude that the SRHBM is more available and effective. Ultimately, we consider that the SRHBM can be used to deal with more complex strongly nonlinear vibration problems.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Authors’ Contributions
All the authors contributed equally and significantly to the writing of this paper. All the authors read and approved the final manuscript.
Acknowledgments
The authors gratefully acknowledge the support of the National Natural Science Foundation of China (NNSFC) through Grants nos. 11572288, 11302184, and 11372015.
References
 R. E. Mickens, Mathematical Methods for the Natural and Engineering Sciences, World Scientific, Singapore, 2004.
 O. Abdulaziz, N. F. M. Noor, and I. Hashim, “Homotopy analysis method for fully developed MHD micropolar fluid flow between vertical porous plates,” International Journal for Numerical Methods in Engineering, vol. 78, no. 7, pp. 817–827, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 Y. M. Chen and J. K. Liu, “Homotopy analysis method for limit cycle flutter of airfoils,” Applied Mathematics and Computation, vol. 203, no. 2, pp. 854–863, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 R. R. Pušenjak, “Extended LindstedtPoincare method for nonstationary resonances of dynamical systems with cubic nonlinearities,” Journal of Sound & Vibration, vol. 314, no. 12, pp. 194–216, 2008. View at: Google Scholar
 P. Amore and A. Aranda, “Improved LindstedtPoincaré method for the solution of nonlinear problems,” Journal of Sound and Vibration, vol. 283, no. 3–5, pp. 1115–1136, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 M. Senator and C. N. Bapat, “A perturbation technique that works even when the nonlinearity is not small,” Journal of Sound and Vibration, vol. 164, no. 1, pp. 1–27, 1993. View at: Publisher Site  Google Scholar  MathSciNet
 Y. K. Cheung, S. H. Chen, and S. L. Lau, “A modified LindstedtPoincar\'e method for certain strongly nonlinear oscillators,” International Journal of NonLinear Mechanics, vol. 26, no. 34, pp. 367–378, 1991. View at: Publisher Site  Google Scholar  MathSciNet
 Q.Q. Hu, C. W. Lim, and L.Q. Chen, “Nonlinear vibration of a cantilever with a DerjaguinMüllerToporov contact end,” International Journal of Structural Stability and Dynamics, vol. 8, no. 1, pp. 25–40, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 K. Huseyin and R. Lin, “An intrinsic multiplescale harmonic balance method for nonlinear vibration and bifurcation problems,” International Journal of NonLinear Mechanics, vol. 26, no. 5, pp. 727–740, 1991. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L.L. Ke, J. Yang, S. Kitipornchai, and Y. Xiang, “Flexural vibration and elastic buckling of a cracked timoshenko beam made of functionally graded materials,” Mechanics of Advanced Materials and Structures, vol. 16, no. 6, pp. 488–502, 2009. View at: Publisher Site  Google Scholar
 S. K. Lai, C. W. Lim, B. S. Wu, C. Wang, Q. C. Zeng, and X. F. He, “Newtonharmonic balancing approach for accurate solutions to nonlinear cubicquintic Duffing oscillators,” Applied Mathematical Modelling, vol. 33, no. 2, pp. 852–866, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 S. J. Liao, “An approximate solution technique not depending on small parameters: a special example,” International Journal of NonLinear Mechanics, vol. 30, no. 3, pp. 371–380, 1995. View at: Publisher Site  Google Scholar  MathSciNet
 S. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, vol. 2 of CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2003. View at: MathSciNet
 S. J. Liao, “Comparison between the homotopy analysis method and homotopy perturbation method,” Applied Mathematics and Computation, vol. 169, no. 2, pp. 1186–1194, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 S. J. Liao, “An optimal homotopyanalysis approach for strongly nonlinear differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 8, pp. 2003–2016, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 C. W. Lim, R. Xu, S. K. Lai, Y. M. Yu, and Q. Yang, “Nonlinear free vibration of an elasticallyrestrained beam with a point mass via the newtonharmonic balancing approach,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, no. 5, pp. 661–674, 2009. View at: Google Scholar
 T. Pirbodaghi, S. H. Hoseini, M. T. Ahmadian, and G. H. Farrahi, “Duffing equations with cubic and quintic nonlinearities,” Computers and Mathematics with Applications, vol. 57, no. 3, pp. 500–506, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. F. Seelig, “Unrestricted harmonic balance II. Application to stiff ordinary differential equations in enzyme catalysis,” Journal of Mathematical Biology, vol. 12, no. 2, pp. 187–198, 1981. View at: Publisher Site  Google Scholar  MathSciNet
 J. L. Summers and M. D. Savage, “Two timescale harmonic balance. I. Application to autonomous onedimensional nonlinear oscillators,” Philosophical Transactions of the Royal Society B: Biological Sciences, vol. 340, no. 1659, pp. 473–501, 1992. View at: Google Scholar
 R. A. Van Gorder and K. Vajravelu, “On the selection of auxiliary functions, operators, and convergence control parameters in the application of the Homotopy Analysis method to nonlinear differential equations: a general approach,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 12, pp. 4078–4089, 2009. View at: Publisher Site  Google Scholar
 H. Wagner, “Largeamplitude free vibrations of a beam,” Journal of Applied Mechanics, vol. 32, no. 4, pp. 887–892, 1965. View at: Google Scholar  MathSciNet
 C. Wang and I. Pop, “Analysis of the flow of a powerlaw fluid film on an unsteady stretching surface by means of homotopy analysis method,” Journal of NonNewtonian Fluid Mechanics, vol. 138, no. 23, pp. 161–172, 2006. View at: Publisher Site  Google Scholar
 B. Wu and P. Li, “A method for obtaining approximate analytic periods for a class of nonlinear oscillators,” Meccanica, vol. 36, no. 2, pp. 167–176, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. Yabushita, M. Yamashita, and K. Tsuboi, “An analytic solution of projectile motion with the quadratic resistance law using the homotopy analysis method,” Journal of Physics. A. Mathematical and Theoretical, vol. 40, no. 29, pp. 8403–8416, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. S. Ganji, D. D. Ganji, H. Babazadeh, and N. Sadoughi, “Application of amplitudefrequency formulation to nonlinear oscillation system of the motion of a rigid rod rocking back,” Mathematical Methods in the Applied Sciences, vol. 33, no. 2, pp. 157–166, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. D. Ganji, H. B. Rokni, M. G. Sfahani, and S. S. Ganji, “Approximate traveling wave solutions for coupled WhithamBroerKaup shallow water,” Advances in Engineering Software, vol. 41, no. 78, pp. 956–961, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. S. Ganji, A. Barari, and D. D. Ganji, “Approximate analysis of twomassspring systems and buckling of a column,” Computers and Mathematics with Applications, vol. 61, no. 4, pp. 1088–1095, 2011. View at: Publisher Site  Google Scholar
 S. S. Ganji, A. Barari, L. B. Ibsen, and G. Domairry, “Differential transform method for mathematical modeling of jamming transition problem in traffic congestion flow,” Central European Journal of Operations Research, vol. 20, no. 1, pp. 87–100, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 M. N. Hamdan and N. H. Shabaneh, “On the large amplitude free vibrations of a restrained uniform beam carrying an intermediate lumped mass,” Journal of Sound and Vibration, vol. 199, no. 5, pp. 711–736, 1997. View at: Publisher Site  Google Scholar
 S. Telli and O. Kopmaz, “Free vibrations of a mass grounded by linear and nonlinear springs in series,” Journal of Sound and Vibration, vol. 289, no. 45, pp. 689–710, 2006. View at: Publisher Site  Google Scholar
 S. K. Lai and C. W. Lim, “Accurate approximate analytical solutions for nonlinear free vibration of systems with serial linear and nonlinear stiffness,” Journal of Sound and Vibration, vol. 307, no. 3, pp. 720–736, 2007. View at: Publisher Site  Google Scholar
 S.S. Chen and C.K. Chen, “Application of the differential transformation method to the free vibrations of strongly nonlinear oscillators,” Nonlinear Analysis. Real World Applications, vol. 10, no. 2, pp. 881–888, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 S. S. Ganji, D. D. Ganji, M. G. Sfahani, and S. Karimpour, “Application of AFF and HPM to the systems of strongly nonlinear oscillation,” Current Applied Physics, vol. 10, no. 5, pp. 1317–1325, 2010. View at: Publisher Site  Google Scholar
 I. Mehdipour, D. D. Ganji, and M. Mozaffari, “Application of the energy balance method to nonlinear vibrating equations,” Current Applied Physics, vol. 10, no. 1, pp. 104–112, 2010. View at: Publisher Site  Google Scholar
 Y. H. Qian, S. K. Lai, W. Zhang, and Y. Xiang, “Study on asymptotic analytical solutions using HAM for strongly nonlinear vibrations of a restrained cantilever beam with an intermediate lumped mass,” Numerical Algorithms, vol. 58, no. 3, pp. 293–314, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Guo and W. Zhang, “The spreading residue harmonic balance study on the vibration frequencies of tapered beams,” Applied Mathematical Modelling, vol. 40, no. 15, pp. 7195–7203, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Guo and X. Ma, “Residue harmonic balance solution procedure to nonlinear delay differential systems,” Applied Mathematics and Computation, vol. 237, pp. 20–30, 2014. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. Ju and X. Xue, “Global residue harmonic balance method to periodic solutions of a class of strongly nonlinear oscillators,” Applied Mathematical Modelling, vol. 38, no. 24, pp. 6144–6152, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 P. Ju and X. Xue, “Global residue harmonic balance method for largeamplitude oscillations of a nonlinear system,” Applied Mathematical Modelling, vol. 39, no. 2, pp. 449–454, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Y. Lee, “Free vibration analysis of a nonlinear panel coupled with extended cavity using the multilevel residue harmonic balance method,” ThinWalled Structures, vol. 98, pp. 332–336, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Y. H. Qian 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.