Higher Order Harmonic Balance Method Research Articles

Analysis on nonlinear stiffness isolators revealing damping thresholds

Nonlinear vibration isolators with stiffness nonlinearity show promise for a broadband isolation performance without degrading the load capacity. An analytical method is proposed to predict all the possible frequency responses of nonlinear stiffness isolators under different damping. The vibration equation is transformed into an algebraic equation through harmonic balance method. The equation is regarded as a function of the frequency instead of the amplitude, leading to a general analytical method with its complexity independent with the order of the nonlinear stiffness. A damping threshold function is proposed based on the root conditions of the frequency, and the whole damping region is divided into large damping, medium damping and small damping sub-regions accordingly. The root conditions and the frequency responses in each sub-region are analytically analyzed. All the possible responses are classified into five categories based on the damping thresholds and the existence of a jump phenomenon. Nonlinear phenomena of full-band isolation, bounded response and unbounded response are revealed according to multivalued roots of the frequency, and their sufficient and necessary conditions are presented. Simulations are performed on nonlinear isolators with cubic stiffness, high-order stiffness, and non-polynomial stiffness, respectively. The damping thresholds, frequency response categories, and the nonlinear characteristics of each category are demonstrated for the studied cases. A high-order harmonic balance method and the Runge–Kutta method are adopted to verify the effectiveness of the proposed method. This work presents a general analytical method for the prediction of frequency responses, the determination of damping thresholds, and the theoretical illustrations of unique behaviors for nonlinear stiffness isolators.

Analysis of super-harmonic resonance and periodic motion transition of fractional nonlinear vibration isolation system

The precision instruments and equipment are often utilized in low-frequency and micro-amplitude vibration systems, in which many vibration isolators of rubber materials are widely used. Ignoring the low-frequency amplitude will result in errors in the fatigue life design of the vibration isolators and predicting the dynamic response of each frequency band accurately becomes necessary. However, integer-order models cannot describe the frequency dependence of rubber materials, while the fractional-order models can describe it instead. On the other hand, the elastic restoring force is strongly nonlinear under large deformation, and the vibration of the nonlinear system contains multiple harmonic components. In order to solve those issues, the fractional nonlinear Nishimura model is used to characterize the constitutive relation of vibration isolators such as air springs, which are made of carbon black filled natural rubber. The high-order harmonic balance method is used to obtain the steady-state response of the vibration system, while the fourth-order Runge–Kutta method is applied to simulate the dynamic response of the system in the low-frequency region, and the Lyapunov exponent is used to determine the stability of the system. Furthermore, the influence of parameters on the amplitude–frequency characteristics of different frequency bands is also studied, and a method to solve the optimal damping coefficient is proposed based on the primary resonance amplitude–frequency curves. The results show that there is a diversity of periodic motions in the process of adjacent super-harmonic resonance transition. Numerical simulations also demonstrate that multi-periodic motions coexist in the system. The motion transition law between the polymorphic coexistence region and its adjacent regions is summarized.

Vibration Prediction of Bladed Disks Coupled by Friction Joints

The present review article addresses the vibration behavior of bladed disks encountered e.g. in aircraft engines as well as industrial gas and steam turbines. The utilization of the dissipative effects of dry friction in mechanical joints is a common means of the passive mitigation of structural vibrations caused by aeroelastic excitation mechanisms. The prediction of the vibration behavior is a scientific challenge due to (a) the strongly nonlinear contact interactions involving local sticking, sliding and liftoff, (b) the model order required to accurately describe the dynamic behavior of the assembly, and (c) the multi-disciplinary character of the problem associated with the need to account for structural mechanical as well as fluid dynamical effects. The purpose of this article is the overview and discussion the current state of the art of vibration prediction approaches. The modeling approaches in this work embrace the description of the rotating bladed disk, the contact modeling, the consideration of aeroelastic effects, appropriate model order reduction techniques and the exploitation of the rotationally periodic nature of the problem. The simulation approaches cover the direct computation of periodic, steady-state externally forced and self-excited vibrations using the high-order harmonic balance method, the formulation of the contact problem in the frequency domain, methods for the solution of the governing algebraic equations and advanced simulation approaches, including the concept of nonlinear modes.

Period Doubling Motions of a Nonlinear Rotating Beam at 1:1 Resonance

A rotating beam subjected to a torsional excitation is studied in this paper. Both quadratic and cubic geometric stiffening nonlinearities are retained in the equation of motion, and the reduced model is obtained via the Galerkin method. Saddle-node bifurcations and Hopf bifurcations of the period-1 motions of the model were obtained via the higher order harmonic balance method. The period-2 and period-4 solutions, which are emanated from the period-1 and period-2 motions, respectively, are obtained by the combined implementation of the harmonic balance method, Floquet theory, and Discrete Fourier transform (DFT). The analytical periodic solutions and their stabilities are verified through numerical simulation.

Dual Time Stepping Algorithms With the High Order Harmonic Balance Method for Contact Interfaces With Fretting-Wear

Contact interfaces with dry friction are frequently used in turbomachinery. Dry friction damping produced by the sliding surfaces of these interfaces reduces the amplitude of bladed-disk vibration. The relative displacements at these interfaces lead to fretting-wear which reduces the average life expectancy of the structure. Frequency response functions are calculated numerically by using the multi-harmonic balance method (mHBM). The dynamic Lagrangian frequency-time method is used to calculate contact forces in the frequency domain. A new strategy for solving nonlinear systems based on dual time stepping is applied. This method is faster than using Newton solvers. It was used successfully for solving Nonlinear CFD equations in the frequency domain. This new approach allows identifying the steady state of worn systems by integrating wear rate equations a on dual time scale. The dual time equations are integrated by an implicit scheme. Of the different orders tested, the first order scheme provided the best results.

Harmonic Balance Approach for an Airfoil with a Freeplay Control Surface

The nonlinear aeroelastic response of a two-dimensional airfoil, including a control surface with freeplay placed in an incompressible flow, is studied. The model equations are formulated as a set of first-order ordinary differential equations. First, the dynamic response is investigated by a time integration method, and the time integration results are used for the verification of the harmonic balance results. The interesting hysteresis phenomenon and the effect of initial conditions of the subcritical bifurcation are presented. A higher-order harmonic balance method is then derived to investigate the high harmonics of the airfoil motions. The harmonic balance prediction is verified by comparison to the results from a numerical time marching integration and also by comparison to results from a previous experiment. Nomenclature a = nondimensional distance of the elastic axis from the midchord, with respect to the semichord B = damping submatrix b = semichord Ch = stiffness (per unit span) of wing in deflection C(k) = generalized Theodorsen function Cα, Cβ = torsional stiffness (per unit span) of wing around a and of aileron around c c = nondimensional distance of the control surface (aileron) hinge line from the midchord, with respect to the semichord

Double Hopf Bifurcation Analysis Using Frequency Domain Methods

The dynamic behavior close to a non-resonant double Hopf bifurcation is analyzed via a frequency-domain technique. Approximate expressions of the periodic solutions are computed using the higher order harmonic balance method while their accuracy and stability have been evaluated through the calculation of the multipliers of the monodromy matrix. Furthermore, the detection of secondary Hopf or torus bifurcations (Neimark–Sacker bifurcation for maps) close to the analyzed singularity has been obtained for a coupled electrical oscillatory circuit. Then, quasi-periodic solutions are likely to exist in certain regions of the parameter space. Extending this analysis to the unfolding of the 1:1 resonant double Hopf bifurcation, cyclic fold and torus bifurcations have also been detected in a controlled oscillatory coupled electrical circuit. The comparison of the results obtained with the suggested technique, and with continuation software packages, has been included.

The Secondary Bifurcation of an Aeroelastic Airfoil Motion: Effect of High Harmonics

The nonlinear dynamical response of a two-degree-of-freedom aeroelastic airfoil motion with cubic restoring forces is investigated. A secondary bifurcation after the primary Hopf (flutter) bifurcation is detected for a cubic hard spring in the pitch degree-of-freedom. Furthermore, there is a hysteresis in the secondary bifurcation: starting from different initial conditions the motion may jump from one limit cycle to another at different fluid flow velocities. A high-order harmonic balance method is employed to investigate the possible bifurcation branches. Furthermore, a numerical time simulation procedure is used to confirm the stable and unstable bifurcation branches.