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

作　　者：Yu ZHOU Li WANG Jianliang HUANG 

机构地区：[1]School of Aeronautics and Astronautics,Sun Yat-sen University,Shenzhen 518107,Guangdong Province,China

出　　处：《Applied Mathematics and Mechanics(English Edition)》2025年第5期907-926,共20页应用数学和力学(英文版)

基　　金：Project supported by the National Natural Science Foundation of China(No.12372028);the Guangdong Basic and Applied Basic Research Foundation(No.2022A1515011809)。

摘　　要：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.

关 键 词：non-smooth system vibro-impact system harmonic balance(HB)method event-driven Gauss quadrature non-smooth Bernoulli base 

分 类 号：O322[理学—一般力学与力学基础]