Fast Fourier nonlinear vibration analysis

Abstract

We present an implementation of the multiharmonic balance method (MHB) where intensive use of the Fast Fourier Transform algorithm (FFT) is made at all stages of calculations. The MHB method is not modi®ed in essence, but computations are organized to obtain a very attractive method that can be applied systematically on general nonlinear vibration problems. The resulting nonlinear algebraic problem is solved by a particular implementation of a continuation method. Nonlinear vibration results are analyzed a posteriori by a Floquet method to determine their stability. The technique is applied on a series of problems of different nature, demonstrating the robustness and¯exibility of the approach. IntroductionThe multi-harmonic balance method (MHB) has been widely used to solve nonlinear vibration problems under periodic excitation. It ®nds applications in several ®elds of mechanical engineering, e.g. machine dynamics, vehicle dynamics, helicopter rotor blade analysis, structural dynamics.The method is well-known from literature. Urabe [1] investigated the convergence conditions of the method and presented numerical applications [2] using the classical Fourier transform. Lau et al. [3±5] and more recently [6] developed an incremental form of the method. Pierre et al. [7, 8] and Ferri [9] followed the same approach to solve nonlinear vibration problems involving dry-friction. Ling and Wu [10] introduced the use of the Fast Fourier Transform (FFT) algorithm in Urabe's formulation. The same path was followed by other authors [11±16]. The use of FFT permits an overwhelming gain of CPU time when computing the Fourier transform. Most authors use the FFT to switch displacements, velocities, accelerations and forces between the time and frequency domains. However, even if the original proposal by Urabe included the analytical expression of the Jacobian of the nonlinear algebraic problem, most of them avoid computing this matrix. Ling and Wu [10] used a Broyden method; Cameron and Grif®n [11] used either a Picard iteration or a Jacobian matrix evaluated by ®nite differences; Lewandowski [14] gave particular expressions for the contributions of nonlinear terms to the Jacobian for the case of geometric nonlinearities in the form of trigonometric expansions.In reference [15] we presented an approach which differs from most others in that an analytical expression for the Jacobian matrix of the nonlinear algebraic problem is developed for general applications. This fact allows to solve the nonlinear algebraic problem with utmost ef®-ciency, reaching quadratic convergence rate. We have shown how the Jacobian matrix can be computed from the Fourier transform of the time domain stiffness, damping and mass matrices of the system under analysis.The MHB leads to a nonlinear algebraic problem in which the solution should be searched in the Fourier transformed displacements versus excitation period space. The solution is found in the form of a nonlinear dynamic equilibrium path in this space for varying values of the excitation period. St...

Nonlinear Forced Response of a Stator Vane With Multiple Friction Contacts Using a Coupled Static/Dynamic Approach

“…For this reason, ad hoc numerical codes must be developed in order to compute the forced response in the frequency domain. These codes are based on the Harmonic Balance Method (HBM) [12]. The periodic variables (displacements and forces) are written as a sum of harmonic terms by Fourier analysis and then the balance of each harmonic is imposed, turning the original nonlinear differential equations in a set of nonlinear algebraic equations.…”

“…In order to reduce the calculation time of a numerical integration, the harmonic balance method (HBM) can be used to solve the equations of motion of the system [5,12,13,14]. This is possible due to the periodicity of the external excitation that involves also the periodicity of the displacements Q and the non-linear forces F nl .…”

Modelling Friction Contacts in Structural Dynamics and its Application to Turbine Bladed Disks

“…In order to reduce the calculation times typical of numerical integration of non-linear systems, the harmonic balance method (HBM) can be used to compute the steady-state response of the system (Cardona et al, 1998;Griffin, 1980;Petrov & Ewins 2003). In detail, due to the periodicity of the external excitation, also the displacements Q and the non-linear forces F NL are periodical at steady-state.…”

“…For this reason, ad hoc numerical codes must be developed in order to compute the forced response in the frequency domain. These codes are based on the Harmonic Balance Method (HBM) (Cardona et al 1998): the periodic variables (displacements and forces) are expressed as a superposition of harmonic terms by Fourier analysis and then the balance of each harmonic is imposed, turning the original nonlinear differential equations in a set of nonlinear algebraic equations. In order to couple the bodies in contact when operating in the frequency domain, several contact elements have been developed.…”

Latest investigations on underplatform damper inner mechanics

“…blade pairs). In order to reduce the calculation time typical of numerical integration of non-linear systems, the harmonic balance method (HBM) can be used to compute the steady-state response of the system [1][2][3]. In detail, due to the periodicity of the external excitation, also the displacements and the non-linear forces are periodical at steady-state, hence the displacement and friction forces can be approximated by the first terms of their Fourier series.…”