Periodic response and stability analysis of vibro-impact systems by an enriched harmonic balance method

doi:10.1007/s10483-025-3253-8

Abstract

Abstract:

A vibro-impact system is a hot topic in the study on nonlinear dynamics due to its generality and importance in engineering. In general, the alternating frequency-time harmonic balance (AFT-HB) method can be used to solve elastic collision. However, since the system is non-smooth, the required Fourier/harmonic truncation order is high in order to achieve the theoretical convergence rate, resulting in expensive computational cost. Furthermore, for rigid body collision, the periodic response of the system cannot be solved with the AFT-HB method due to the discontinuous velocity of the system. In order to accelerate the convergence and solve highly non-smooth systems, an enriched harmonic balance (HB) method is proposed, which is derived from the AFT-HB method in the framework of event-driven Gauss quadrature. The basic idea is to augment the Fourier bases by introducing a non-smooth Bernoulli base such that the non-smooth Bernoulli base compensates for the non-smooth part of the solution and the smooth part of the solution is approximated by the Fourier bases, thus achieving accelerated convergence. Based on the enriched HB method, gear pair systems with gear backlash and oscillator systems with rigid impact are solved, and the dynamic response characteristics are analyzed in this work. Then, based on the Floquet theory, the event-driven monodromy matrix method for non-smooth systems is used to analyze the stability and bifurcation of the periodic solutions. The numerical example shows that the results obtained from the enriched HB method are consistent with those from the Runge-Kutta method, which proves that the presented method is an effective method for analyzing the dynamic response characteristic of the vibro-impact system.

Key words: non-smooth system, vibro-impact system, harmonic balance (HB) method, event-driven Gauss quadrature, non-smooth Bernoulli base

2010 MSC Number:

O322

Yu ZHOU, Li WANG, Jianliang HUANG. Periodic response and stability analysis of vibro-impact systems by an enriched harmonic balance method. Applied Mathematics and Mechanics (English Edition), 2025, 46(5): 907-926.

TrendMD

Figures/Tables 14

Fig. 1

Fig. 2

Fig. 3

Fig. 4

Fig. 5

Fig. 6

Fig. 7

Fig. 8

Fig. 9

Fig. 10

Fig. 11

Fig. 12

Fig. 13

Fig. 14

References 42

[1]	LUO, G. W., LV, X. H., and ZHU, X. F.Dynamics of vibro-impact mechanical systems with large dissipation. International Journal of Mechanical Sciences, 50, 214–232 (2008)
[2]	LEINE, R. I. and NIJMEIJER, H.Dynamics and Bifurcations of Non-Smooth Mechanical Systems, Springer, New York, 3–4 (2004)
[3]	MA, S. C., NING, X., WANG, L., JIA, W. T., and XU, W.Complex response analysis of a non-smooth oscillator under harmonic and random excitations. Applied Mathematics and Mechanics (English Edition), 42(5), 641–648 (2021) https://doi.org/10.1007/s10483-021-2731-5
[4]	SCHREYER, F. and LEINE, R. I.A mixed shooting-harmonic balance method for unilaterally constrained mechanical systems. Archive of Mechanical Engineering, 63(2), 297–314 (2016)
[5]	BUTCHER, J. C.A history of Runge-Kutta methods. Applied Numerical Mathematics, 20, 247–260 (1996)
[6]	MICKENS, R. E.Comments on the method of harmonic balance. Journal of Sound and Vibration, 94(3), 456–460 (1984)
[7]	LIN, R. Z., HOU, L., CHEN, Y., JIN, Y. H., SAEED, N. A., and CHEN, Y. S.A novel adaptive harmonic balance method with an asymptotic harmonic selection. Applied Mathematics and Mechanics (English Edition), 44(11), 1887–1910 (2023) https://doi.org/10.1007/s10483-023-3047-6
[8]	WANG, Y., ZHANG, L., and ZHOU, J.Incremental harmonic balance method for periodic forced oscillation of a dielectric elastomer balloon. Applied Mathematics and Mechanics (English Edition), 41(3), 459–470 (2020) https://doi.org/10.1007/s10483-020-2590-7
[9]	KRACK, M. and GROSS, J.Harmonic Balance for Nonlinear Vibration Problems, Springer, Cham, 62–69 (2019)
[10]	WANG, L., LU, Z. R., and LIU, J. K.Convergence rates of harmonic balance method for periodic solution of smooth and non-smooth systems. Communications in Nonlinear Science and Numerical Simulation, 99, 105826 (2021)
[11]	CAMERON, T. M. and GRIFFIN, J. H.An alternating frequency/time domain method for calculating the steady state response of nonlinear dynamic systems. Journal of Applied Mechanics, 56, 149–154 (1989)
[12]	YUAN, T. C., YANG, J., and CHEN, L. Q.A harmonic balance approach with alternating frequency/time domain progress for piezoelectric mechanical systems. Mechanical Systems and Signal Processing, 120, 274–289 (2019)
[13]	WANG, Q. T., ZHANG, Z. Y., YING, Y. H., and PANG, Z. J.Analysis of the dynamic stiffness, hysteresis resonances and complex responses for nonlinear spring systems in power-form order. Applied Sciences, 11, 7722 (2021)
[14]	TIWARI, M., GUPTA, K., and PRAKASH, O.Effect of radial internal clearance of a ball bearing on the dynamics of a balanced horizontal rotor. Journal of Sound and Vibration, 238(5), 723–756 (2000)
[15]	GUSKOV, M., SINOU, J. J., and THOUVEREZ, F.Multi-dimensional harmonic balance applied to rotor dynamics. Mechanics Research Communications, 35, 537–545 (2008)
[16]	LI, H. L., CHEN, Y. S., HOU, L., and ZHANG, Z. Y.Periodic response analysis of a misaligned rotor system by harmonic balance method with alternating frequency/time domain technique. Science China Technological Sciences, 59(11), 1717–1729 (2016)
[17]	KIM, Y. B. and NOAH, S. T.Response and bifurcation analysis of a MDOF rotor system with a strong nonlinearity. Nonlinear Dynamics, 2, 215–234 (1991)
[18]	SINOU, J. J.Non-linear dynamics and contacts of an unbalanced flexible rotor supported on ball bearings. Mechanism and Machine Theory, 44(9), 1713–1732 (2009)
[19]	ZHANG, Z. Y., CHEN, Y. S., and LI, Z. G.Influencing factors of the dynamic hysteresis in varying compliance vibrations of a ball bearing. Science China Technological Sciences, 58, 775–782 (2015)
[20]	ZHANG, Z., NIU, M. Q., YUAN, K., and ZHANG, Y. W.Research on linear/nonlinear viscous damping and hysteretic damping in nonlinear vibration isolation systems. Applied Mathematics and Mechanics (English Edition), 41(7), 983–998 (2020) https://doi.org/10.1007/s10483-020-2630-6
[21]	LI, M. and DING, H.A vertical track nonlinear energy sink. Applied Mathematics and Mechanics (English Edition), 45(6), 931–946 (2024) https://doi.org/10.1007/s10483-024-3127-6
[22]	ZHOU, Y., HUANG, J. L., and WANG, L.Event-driven Gauss quadrature and stability analysis for fast alternating frequency-time harmonic balance of non-smooth systems. Communications in Nonlinear Science and Numerical Simulation, 120, 107189 (2023)
[23]	BANERJEE, N. S. and GEER, J. F.Exponentially accurate approximations to periodic Lipschitz functions based on Fourier series partial sums. Journal of Scientific Computing, 13(4), 419–460 (1998)
[24]	NERSESSIAN, A.Acceleration of convergence of Fourier series using the phenomenon of over-convergence. Armenian Journal of Mathematics, 14(14), 1-31 (2022)
[25]	LANCZOS, C.Discourse on Fourier Series, Oliver and Boyd, Edinburgh (1966)
[26]	GOTTLIEB, D.Issues in the application of high order schemes. Algorithmic Trends in Computational Fluid Dynamics, Springer, New York, 195–218 (1991)
[27]	GOTTLIEB, D. and SHU, C. W.On the Gibbs phenomenon and its resolution. SIAM Review, 39(4), 644–668 (1997)
[28]	GOTTLIEB, D., SHU, C. W., SOLOMONOFF, A., and VANDEVEN, H.On the Gibbs phenomenon I: recovering exponential accuracy from the Fourier partial sum of a nonperiodic analytic function. Journal of Computational and Applied Mathematics, 43, 81–98 (1992)
[29]	GELB, A.A hybrid approach to spectral reconstruction of piecewise smooth functions. Journal of Scientific Computing, 15(3), 293–322 (2000)
[30]	ECKHOFF, K. S.Accurate reconstructions of functions of finite regularity from truncated Fourier series expansions. Mathematics of Computation, 64, 671–690 (1995)
[31]	ECKHOFF, K. S.Accurate and efficient reconstruction of discontinuous functions from truncated series expansions. Mathematics of Computation, 61, 745–763 (1993)
[32]	ECKHOFF, K. S.On a high order numerical method for functions with singularities. Mathematics of Computation, 67, 1063–1087 (1998)
[33]	BATENKOV, D. and YOMDIN, Y.Local and global geometry of Prony systems and Fourier reconstruction of piecewise-smooth functions. Operator-Related Function Theory and Time-Frequency Analysis, Springer, Switzerland, 57–76 (2015)
[34]	CAI, W., GOTTLIEB, D., and SHU, C. W.Essentially nonoscillatory spectral Fourier methods for shock wave calculations. Mathematics of Computation, 52, 389–410 (1989)
[35]	GOTTLIEB, D., LUSTMAN, L., and ORSZAG, S. A.Spectral calculations of one-dimensional inviscid compressible flows. SIAM Journal on Scientific and Statistical Computing, 2, 296–310 (1981)
[36]	LING, F. H. and WU, X. X.Fast Galerkin method and its application to determine periodic solutions of non-linear oscillators. International Journal of Non-Linear Mechanics, 22(2), 89–98 (1987)
[37]	KIM, W. J. and PERKINS, N. C.Harmonic balance/Galerkin method for non-smooth dynamic systems. Journal of Sound and Vibration, 261(2), 213–224 (2003)
[38]	NATALINI, P. and BERNARDINI, A.A generalization of the Bernoulli polynomials. Journal of Applied Mathematics, 2003, 155–163 (2003)
[39]	TREFETHEN, L. N.Is Gauss quadrature better than Clenshaw-Curtis?SIAM Review, 50, 67–87 (2008)
[40]	GOU, X. F., ZHU, L. Y., and CHEN, D. L.Bifurcation and chaos analysis of spur gear pair in two-parameter plane. Nonlinear Dynamics, 79, 2225–2235 (2015)
[41]	LI, H. B., HU, J. J., SHI, Y. T., and LIU, S.Dynamic behavior analysis and time delay feedback control of gear pair system with backlash non-smooth characteristic. Journal of Vibroengineering, 19, 302–313 (2017)
[42]	SEYDEL, R.Practical Bifurcation and Stability Analysis, Springer, New York, 314–316 (1994)