An efficient Galerkin averaging-incremental harmonic balance method for nonlinear dynamic analysis of rigid multibody systems governed by differential–algebraic equations

Author(s):  
R. Ju ◽  
W. Fan ◽  
W. D. Zhu
2021 ◽  
Author(s):  
Ren Ju ◽  
Wei Fan ◽  
Weidong Zhu

Abstract The bridge between the multibody dynamic modeling theory and nonlinear dynamic analysis theory is built for the first time in this work by introducing an efficient Galerkin averaging-incremental harmonic balance (EGA-IHB) method for steady-state nonlinear dynamic analysis of index-3 differential algebraic equations (DAEs) for general rigid multibody systems. The multibody dynamic modeling theory has made significant advances in generality and simplicity, and multibody systems are usually governed by DAEs. Since the fast Fourier transform and EGA are used, the EGA-IHB method has excellent robustness and computational efficiency. Since the Floquet theory cannot be directly used for stability analysis of periodic responses of DAEs, a new stability analysis procedure is developed, where perturbed, linearized DAEs are reduced to ordinary differential equations with use of independent generalized coordinates. A modified arc-length continuation method with a scaling strategy is used for calculating response curves and conducting parameter studies. Three examples are used to show the performance and capability of the current method. Periodic solutions of DAEs from the EGA-IHB method show excellent agreement with those from numerical integration methods. Amplitude-frequency and amplitude-parameter response curves are generated, and stability and period-doubling bifurcations are analyzed. The EGA-IHB method can be used as a universal solver and nonlinear analyzer for obtaining steady-state periodic responses of DAEs for general multibody systems.


Author(s):  
Ren Ju ◽  
Weidong Zhu

Abstract Modern computers are generally equipped with multi-core central processing units (CPUs). There has been great interest to introduce a parallelized harmonic balance method to make full use of computing resources and improve the efficiency in dealing with nonlinear dynamic analysis of high-dimensional spatially discretized models of continuous systems. In this work, an optimized efficient Galerkin averaging-incremental harmonic balance method for solving high-dimensional models of continuous systems based on parallel computing is introduced. Optimized parallel implementation based on tensor contraction is introduced in time-domain series calculations and the quasi-Newton method is used in the iteration procedure, which greatly accelerate computational speeds of both serial and parallel implementations. Especially, the parallel implementation achieves high parallel efficiency when multiple CPU cores are used. Due to its high computational efficiency and good robustness, the proposed method has the potential to be used as a powerful universal solver and analyzer for general types of continuous systems.


Author(s):  
R. Ju ◽  
W. Fan ◽  
W. D. Zhu ◽  
J. L. Huang

A modified two-timescale incremental harmonic balance (IHB) method is introduced to obtain quasi-periodic responses of nonlinear dynamic systems with combinations of two incommensurate base frequencies. Truncated Fourier coefficients of residual vectors of nonlinear algebraic equations are obtained by a frequency mapping-fast Fourier transform procedure, and complex two-dimensional (2D) integration is avoided. Jacobian matrices are approximated by Broyden's method and resulting nonlinear algebraic equations are solved. These two modifications lead to a significant reduction of calculation time. To automatically calculate amplitude–frequency response surfaces of quasi-periodic responses and avoid nonconvergent points at peaks, an incremental arc-length method for one timescale is extended for quasi-periodic responses with two timescales. Two examples, Duffing equation and van der Pol equation with quadratic and cubic nonlinear terms, both with two external excitations, are simulated. Results from the modified two-timescale IHB method are in excellent agreement with those from Runge–Kutta method. The total calculation time of the modified two-timescale IHB method can be more than two orders of magnitude less than that of the original quasi-periodic IHB method when complex nonlinearities exist and high-order harmonic terms are considered.


2019 ◽  
Vol 87 (1) ◽  
Author(s):  
X. M. Xu ◽  
J. H. Luo ◽  
Z. G. Wu

Abstract The inertia representation of a constrained system includes the formulation of the kinetic energy and its corresponding mass matrix, given the coordinates of the system. The way to find a proper inertia representation achieving better numerical performance is largely unexplored. This paper extends the modified inertia representation (MIR) to the constrained rigid multibody systems. By using the orthogonal projection, we show the possibility to derive the MIR for many types of non-minimal coordinates. We present examples of the derivation of the MIR for both planar and spatial rigid body systems. Error estimation shows that the MIR is different from the traditional inertia representation (TIR) in that its parameter γ can be used to reduce the kinetic energy error. With preconditioned γ, numerical results show that the MIR consistently presents significantly higher numerical accuracy and faster convergence speed than the TIR for the given variational integrator. The idea of using different inertia representations to improve the numerical performance may go beyond constrained rigid multibody systems to other systems governed by differential algebraic equations.


2014 ◽  
Vol 2014 ◽  
pp. 1-8
Author(s):  
Hongliang Yao ◽  
Qian Zhao ◽  
Qi Xu ◽  
Bangchun Wen

The efficiency and accuracy of common time and frequency domain methods that are used to simulate the response of a rotor system with malfunctions are compared and analyzed. The Newmark method and the incremental harmonic balance method are selected as typical representatives of time and frequency domain methods, respectively. To improve the simulation efficiency, the fixed interface component mode synthesis approach is combined with the Newmark method and the receptance approach is combined with the incremental harmonic balance method. Numerical simulations are performed for rotor systems with single and double frequency excitations. The inherent characteristic that determines the efficiency of the two methods is analyzed. The results of the analysis indicated that frequency domain methods are suitable single and double frequency excitation rotor systems, whereas time domain methods are more suitable for multifrequency excitation rotor systems.


Sign in / Sign up

Export Citation Format

Share Document