An Innovative Precise Integration Method in Solving Structural Dynamic Problems

Author(s):  
Chiun-lin Wu ◽  
Ching-Chiang Chuang

An innovative time integration method that incorporates spurious high-frequency dissipation capability into the so called “high precision direct integration algorithm” is presented, and its numerical stability and accuracy is discussed. The integration algorithm is named “high precision” to emphasize its numerical capability in reaching computer hardware precision. The proposed procedure employs the well-known state space approach to solve the simultaneous ordinary differential equations, the exact solution of which contains an exponential matrix to be efficiently computed using the truncated Taylor series expansion together with the power-of-two algorithm. The proposed method, belonging to the category of explicit methods, is found to provide better accuracy than many other existing time integration methods, and the integration scheme remains numerically stable over a wide range of frequencies of engineering interest. This paper is also devoted to study numerical accuracy of the Precise Integration Method in solving forced vibration problems, particularly near resonance conditions. The numerically obtained transfer functions are then compared with the analytical exact solution to detect spurious resonance. Finally, numerical examples are used to illustrate its high performance in numerical stability and accuracy. The proposed method carries the merit that can be directly applied to solve momentum equations of motion with exactly the same procedure.

1997 ◽  
Vol 119 (3) ◽  
pp. 334-340 ◽  
Author(s):  
Zhong Wanxie ◽  
F. W. Williams ◽  
P. N. Bennett

A precise integration algorithm has recently been proposed by Zhong (1994) for dynamic stiffness matrix computations, but he did not give a corresponding eigenvalue count method. The Wittrick-Williams algorithm gives an eigenvalue count method for pure displacement formulations, but the precise integration method uses a mixed variable formulation. Therefore the Wittrick-Williams method is extended in this paper to give the eigenvalue count needed by the precise integration method and by other methods involving mixed variable formulations. A simple Timoshenko beam example is included.


2014 ◽  
Vol 8 (1) ◽  
pp. 264-269
Author(s):  
Guangtian Shi

Through an example of rotor system which has multi-degree of freedom mounted on the nonlinear fluid film bearings, this paper analyzes the precise integration algorithm, a new numerical solution for high–dimensional nonlinear dynamics system. The precise integration method has advantages of long step, high precision and absolute stability for solving nonlinear dynamics equations. To make good use of the method, firstly, the precise integral iterative formula is given and then the mechanism of controlling high precision and efficiency is discussed. The evolution of precise integration method is an algorithm with explicit, simple form, self-start, and fast to solve nonlinear dynamics equations. High power of athwart of Hamiltonian matrix is not needed, so it is convenient in this case. The stability of period response of nonlinear rotor-bearing system is analyzed by employing the precise integration method. The bifurcation rules of the period response of the elastic rotor system with multi-degree of freedom are obtained and the chaos of the system is determined according to the fractal dimension of Poincare mapping of phase space at a certain speed.


1999 ◽  
Vol 20 (3) ◽  
pp. 269-275 ◽  
Author(s):  
Qiang Shizhong ◽  
Wang Xiaoguo ◽  
Tang Maolin ◽  
Liu Min

2013 ◽  
Vol 2013 ◽  
pp. 1-9
Author(s):  
Hai-Jun Peng ◽  
Sheng Zhang ◽  
Zhi-Gang Wu ◽  
Biao-Song Chen

The key of solving the noncooperative linear quadratic (LQ) differential game is to solve the coupled matrix Riccati differential equation. The precise integration method based on the adaptive choosing of the two parameters is expanded from the traditional symmetric Riccati differential equation to the coupled asymmetric Riccati differential equation in this paper. The proposed expanded precise integration method can overcome the difficulty of the singularity point and the ill-conditioned matrix in the solving of coupled asymmetric Riccati differential equation. The numerical examples show that the expanded precise integration method gives more stable and accurate numerical results than the “direct integration method” and the “linear transformation method”.


2018 ◽  
Vol 2018 ◽  
pp. 1-21 ◽  
Author(s):  
Dongdong He ◽  
Qiang Gao ◽  
Wanxie Zhong

Based on the finite element method (FEM), the parametric variational principle (PVP) is combined with a numerical time-domain integral method to simulate the dynamic behavior of the pantograph-catenary system. Based on PVP, formulations for the nonlinear droppers in the catenary and for the contact between the pantograph and the contact wire are proposed. The formulations can accurately determine the tension state or compression state of the nonlinear droppers and the contact state between the pantograph and the contact wire. Based on the periodicity of the catenary and the precise integration method (PIM), a numerical time-integration method is developed for the dynamic responses of the catenary. For this method, the matrix exponential of only one unit cell of the catenary is computed, which greatly improves the computational efficiency. Moreover, the validation shows that the formulations can compute the contact force accurately and represent the nonlinearity of the droppers, which demonstrates the accuracy and reliability of the proposed method. Finally, the dynamic behaviors of the pantograph-catenary system with different types of catenaries are simulated.


Author(s):  
Takeshi Fujikawa ◽  
Etsujiro Imanishi

Abstract A method of time integration algorithm is presented for solving stiff vibration and motion problems. It is absolutely stable, numerically dissipative, and much accurate than other dissipative time integration methods. It achieves high-frequency dissipation, while minimizing unwanted low-frequency dissipation. In this method change of acceleration during time step is expressed as quadratic function including some parameters, whose appropriate values are determined through numerical investigation. Two calculation examples are demonstrated to show the usefulness of this method.


Sign in / Sign up

Export Citation Format

Share Document