Stability of an explicit multi-time step integration algorithm for linear structural dynamics equations

1996 ◽  
Vol 18 (3) ◽  
pp. 236-244 ◽  
Author(s):  
P. Smolinski ◽  
S. Sleith ◽  
T. Belytschko
Author(s):  
Salvatore Lopez

A second-order accurate single-step time integration method for nonlinear structural dynamics is developed. The method combines algorithmic dissipation of higher modes and conservation of linear and angular momentum and is composed of two phases. In the first phase, a solution point is computed by a basic integration scheme, the generalized-[Formula: see text] method being adopted due to its higher level of high-frequency dissipation. In the second phase, a correction is hypothesized as a linear combination of the solution in the basic step and the gradient of vector components of the incremental linear and angular momentum. By solving a system composed of six linear equations, the searched for corrected solution in the time step is then provided. The novelty in the presented integration scheme lies in the way of imposing the conservation of linear and angular momentum. In fact, this imposition is carried out as a correction of the computed solution point in the time step and not through an enlarged system of equations of motion. To perform tests on plane and spatial motion of three-dimensional structural models, a small strains — finite rotations corotational formulation is also described.


2009 ◽  
Vol 131 (2) ◽  
Author(s):  
Shuenn-Yih Chang

An integration algorithm, which integrates the most important advantage of explicit methods of the explicitness of each time step and that of implicit methods of the possibility of unconditional stability, is presented herein. This algorithm is analytically shown to be unconditionally stable for any linear elastic and nonlinear systems except for the instantaneous stiffness hardening systems with the instantaneous degree of nonlinearity larger than 43 based on a linearized stability analysis. Hence, its stability property is better than the previously published algorithm (Chang, 2007, “Improved Explicit Method for Structural Dynamics,” J. Eng. Mech., 133(7), pp. 748–760), which is only conditionally stable for instantaneous stiffness hardening systems although it also possesses unconditional stability for linear elastic and any instantaneous stiffness softening systems. Due to the explicitness of each time step, the possibility of unconditional stability, and comparable accuracy, the proposed algorithm is very promising for a general structural dynamic problem, where only the low frequency responses are of interest since it consumes much less computational efforts when compared with explicit methods, such as the Newmark explicit method, and implicit methods, such as the constant average acceleration method.


Author(s):  
Clemens Bernhard Domnick ◽  
Friedrich-Karl Benra ◽  
Dieter Brillert ◽  
Hans Josef Dohmen ◽  
Christian Musch

The power output of steam turbines is controlled by steam turbine inlet valves. These valves have a large flow capacity and dissipate in throttled operation a huge amount of energy. Due to that, high dynamic forces occur in the valve which can cause undesired valve vibrations. In this paper, the structural dynamics of a valve are analysed. The dynamic steam forces obtained by previous computational fluid dynamic (CFD) calculations at different operating points are impressed on the structural dynamic finite element model (FEM) of the valve. Due to frictional forces at the piston rings and contact effects at the bushings of the valve plug and the valve stem the structural dynamic FEM is highly nonlinear and has to be solved in the time domain. Prior to the actual investigation grid and time step studies are carried out. Also the effect of the temperature distribution within the valve stem is discussed and the influence of the valve actuator on the vibrations is analysed. In the first step, the vibrations generated by the fluid forces are investigated. The effects of the piston rings on the structural dynamics are discussed. It is found, that the piston rings are able to reduce the vibration significantly by frictional damping. In the second step, the effect of the moving valve plug on the dynamic flow in the valve is analysed. The time dependent displacement of the valve is transferred to CFD calculations using deformable meshes. With this one way coupling method the response of the flow to the vibrations is analysed.


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.


2005 ◽  
Vol 18 (2) ◽  
pp. 239-244
Author(s):  
Zheng-hua Zhou ◽  
Yu-huan Wang ◽  
Quan Liu ◽  
Xiao-tao Yin ◽  
Cheng Yang

Author(s):  
Chao Zhang ◽  
Mingxiang Ling ◽  
Meng Tao

Abstract This paper puts forward a computationally-efficient parallel precise integration algorithm for solving vibration response subjected to time-variable excitation and nonlinearity, especially for non-homogenous dynamic response solution with large-scale degree of freedom. In detail, both of the nonlinear parts and time-varying inputs of the dynamic system are separated from the original dynamic equations and then simulated by employing a piecewise interpolation function within a computing time-step. A novel closed-form iteration formula is presented in conjunction with the block matrix strategy and modified increment-dimensional precise integration technique. Interestingly, the presented approach is essentially a high-accuracy and parallel algorithm, which exhibits a high prediction accuracy without the limitation of matrix inversion, higher-order derivative, periodicity requirement nor cycle oscillation and instability of high-order interpolation. At last, the feasibility and advantage of the proposed method is verified with two numerical examples.


Sign in / Sign up

Export Citation Format

Share Document