scholarly journals One-Parameter Controlled Non-Dissipative Unconditionally Stable Explicit Structure-Dependent Integration Methods with No Overshoot

2021 ◽  
Vol 11 (24) ◽  
pp. 12109
Author(s):  
Veerarajan Selvakumar ◽  
Shuenn-Yih Chang
Keyword(s):  
Time Integration ◽  
Unconditional Stability ◽  
Second Order ◽  
Numerical Dissipation ◽  
Implicit Methods ◽  
Time Step ◽  
Order Accuracy ◽  
Structural Dynamic ◽  
Integration Methods ◽  
New Family

Although many families of integration methods have been successfully developed with desired numerical properties, such as second order accuracy, unconditional stability and numerical dissipation, they are generally implicit methods. Thus, an iterative procedure is often involved for each time step in conducting time integration. Many computational efforts will be consumed by implicit methods when compared to explicit methods. In general, the structure-dependent integration methods (SDIMs) are very computationally efficient for solving a general structural dynamic problem. A new family of SDIM is proposed. It exhibits the desired numerical properties of second order accuracy, unconditional stability, explicit formulation and no overshoot. The numerical properties are controlled by a single free parameter. The proposed family method generally has no adverse disadvantage of unusual overshoot in high frequency transient responses that have been found in the currently available implicit integration methods, such as the WBZ-α method, HHT-α method and generalized-α method. Although this family method has unconditional stability for the linear elastic and stiffness softening systems, it becomes conditionally stable for stiffness hardening systems. This can be controlled by a stability amplification factor and its unconditional stability is successfully extended to stiffness hardening systems. The computational efficiency of the proposed method proves that engineers can do the accurate nonlinear analysis very quickly.

A New Unconditionally Stable Time Integration Method for Analysis of Nonlinear Structural Dynamics

2013 ◽  
Vol 80 (2) ◽  
Author(s):  
Ali Akbar Gholampour ◽  
Mehdi Ghassemieh ◽  
Mahdi Karimi-Rad
Keyword(s):  
Time Integration ◽  
Second Order ◽  
Numerical Dispersion ◽  
Computational Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Similar Time ◽  
Unconditionally Stable ◽  
Average Acceleration

A new time integration scheme is presented for solving the differential equation of motion with nonlinear stiffness. In this new implicit method, it is assumed that the acceleration varies quadratically within each time step. By increasing the order of acceleration, more terms of the Taylor series are used, which are expected to have responses with better accuracy than the classical methods. By considering this assumption and employing two parameters δ and α, a new family of unconditionally stable schemes is obtained. The order of accuracy, numerical dissipation, and numerical dispersion are used to measure the accuracy of the proposed method. Second order accuracy is achieved for all values of δ and α. The proposed method presents less dissipation at the lower modes in comparison with Newmark's average acceleration, Wilson-θ, and generalized-α methods. Moreover, this second order accurate method can control numerical damping in the higher modes. The numerical dispersion of the proposed method is compared with three unconditionally stable methods, namely, Newmark's average acceleration, Wilson-θ, and generalized-α methods. Furthermore, the overshooting effect of the proposed method is compared with these methods. By evaluating the computational time for analysis with similar time step duration, the proposed method is shown to be faster in comparison with the other methods.

scholarly journals On the optimization of n-sub-step composite time integration methods

2020 ◽  
Vol 102 (3) ◽  
pp. 1939-1962
Author(s):  
Huimin Zhang ◽  
Runsen Zhang ◽  
Yufeng Xing ◽  
Pierangelo Masarati
Keyword(s):  
Time Integration ◽  
Low Frequency ◽  
Unconditional Stability ◽  
Frequency Content ◽  
Order Accuracy ◽  
Integration Methods ◽  
Universal Approach ◽  
Energy Conserving ◽  
Time Integration Methods ◽  
Linear And Nonlinear Systems

AbstractA family of n-sub-step composite time integration methods, which employs the trapezoidal rule in the first $$n-1$$ n - 1 sub-steps and a general formula in the last one, is discussed in this paper. A universal approach to optimize the parameters is provided for any cases of $$n\ge 2$$ n ≥ 2 , and two optimal sub-families of the method are given for different purposes. From linear analysis, the first sub-family can achieve nth-order accuracy and unconditional stability with controllable algorithmic dissipation, so it is recommended for high-accuracy purposes. The second sub-family has second-order accuracy, unconditional stability with controllable algorithmic dissipation, and it is designed for heuristic energy-conserving purposes, by preserving as much low-frequency content as possible. Finally, some illustrative examples are solved to check the performance in linear and nonlinear systems.

Comparisons of Structure-Dependent Explicit Methods for Time Integration

2015 ◽  
Vol 15 (03) ◽  
pp. 1450055 ◽  
Author(s):  
Shuenn-Yih Chang
Keyword(s):  
Time Integration ◽  
Low Frequency ◽  
Unconditional Stability ◽  
Explicit Methods ◽  
Explicit Method ◽  
Computationally Efficient ◽  
Total Response ◽  
Time Step ◽  
Order Accuracy ◽  
Practical Applications

Chang explicit method (CEM)1,2 and CR explicit method3 (CRM) are two structure-dependent explicit methods that have been successfully developed for structural dynamics. The most important property of both integration methods is that they involve no nonlinear iterations in addition to unconditional stability and second-order accuracy. Thus, they are very computationally efficient for solving inertial problems, where the total response is dominated by low frequency modes. However, an unusual overshooting behavior for CR explicit method is identified herein and thus its practical applications might be largely limited although its velocity computing for each time step is much easier than for the CEM.

A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method

1993 ◽  
Vol 60 (2) ◽  
pp. 371-375 ◽  
Author(s):  
J. Chung ◽  
G. M. Hulbert
Keyword(s):  
Structural Dynamics ◽  
Time Integration ◽  
Low Frequency ◽  
Numerical Dissipation ◽  
Integration Algorithm ◽  
Integration Methods ◽  
New Family ◽  
Time Integration Methods ◽  
Improved Performance ◽  
Dynamics Problems

A new family of time integration algorithms is presented for solving structural dynamics problems. The new method, denoted as the generalized-α method, possesses numerical dissipation that can be controlled by the user. In particular, it is shown that the generalized-α method achieves high-frequency dissipation while minimizing unwanted low-frequency dissipation. Comparisons are given of the generalized-α method with other numerically dissipative time integration methods; these results highlight the improved performance of the new algorithm. The new algorithm can be easily implemented into programs that already include the Newmark and Hilber-Hughes-Taylor-α time integration methods.

New Explicit Integration Algorithms with Controllable Numerical Dissipation for Structural Dynamics

2018 ◽  
Vol 18 (03) ◽  
pp. 1850044 ◽  
Author(s):  
Xiaoqiong Du ◽  
Dixiong Yang ◽  
Jilei Zhou ◽  
Xiaoliang Yan ◽  
Yongliang Zhao ◽  
...  
Keyword(s):  
Nonlinear Systems ◽  
Time Integration ◽  
Dynamic Responses ◽  
Discrete Control ◽  
Numerical Dissipation ◽  
Explicit Integration ◽  
Structural Dynamic ◽  
New Family ◽  
Recursive Formulas ◽  
Integration Algorithms

This paper presents a new family of explicit time integration algorithms with controllable numerical dissipation for structural dynamic problems by utilizing the discrete control theory. Firstly, the equilibrium equation of the implicit Yu-[Formula: see text] algorithm is adopted, and the recursive formulas of velocity and displacement for the explicit CR algorithm are used in the algorithms. Then, the transfer function and characteristic equation of the algorithms with integration coefficients are obtained by the [Formula: see text] transformation. Furthermore, their integration coefficients are derived according to the poles condition. It was indicated that the proposed algorithms possess the advantages of second-order accuracy, self-starting, and unconditional stability for linear systems and nonlinear systems with softening stiffness. The numerical dissipation of the algorithms is controlled by the spectral radius at infinity [Formula: see text]. It was also shown that the proposed algorithms have the same poles as the Yu-[Formula: see text] algorithm, and thus the same numerical properties. Compared with the implicit Yu-[Formula: see text] algorithm, the proposed algorithms are explicit in terms of both the displacement and velocity formulas. Finally, the effectiveness of the proposed algorithms in reducing the undesired participation of higher modes for solving the dynamic responses of linear and nonlinear systems has been demonstrated.

An Efficient Computing Explicit Method for Structural Dynamics

2009 ◽  
Vol 131 (2) ◽  
Author(s):  
Shuenn-Yih Chang
Keyword(s):  
Structural Dynamics ◽  
Unconditional Stability ◽  
Explicit Methods ◽  
Explicit Method ◽  
Implicit Methods ◽  
Time Step ◽  
Integration Algorithm ◽  
Structural Dynamic ◽  
Linear Elastic ◽  
Linearized Stability

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.

Time Integration Method With Structure-Dependent Difference Equations

2013 ◽  
Author(s):  
Shuenn-Yih Chang ◽  
Chiu-Li Huang
Keyword(s):  
High Frequency ◽  
Integration Method ◽  
Time Integration ◽  
Low Frequency ◽  
Unconditional Stability ◽  
Total Response ◽  
Time Step ◽  
Order Accuracy ◽  
Time Integration Method ◽  
Frequency Responses

A novel family of structure-dependent integration method is proposed for time integration. This family method can have the possibility of unconditional stability, second-order accuracy and the explicitness of each time step. Since it can integrate the most important advantage of an implicit method, unconditional stability, and that of an explicit method, the explicitness of each time step, a lot of computational efforts can be saved in solving an inertial type problem, where the total response is dominated by low frequency modes and high frequency responses are of no interest.

scholarly journals Performance of Composite Implicit Time Integration Scheme for Nonlinear Dynamic Analysis

2008 ◽  
Vol 2008 ◽  
pp. 1-16 ◽  
Author(s):  
William Taylor Matias Silva ◽  
Luciano Mendes Bezerra
Keyword(s):  
Transient Response ◽  
Time Integration ◽  
Nonlinear Dynamic Analysis ◽  
Implicit Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Structural Dynamic ◽  
Implicit Time Integration ◽  
Time Integration Scheme

This paper presents a simple implicit time integration scheme for transient response solution of structures under large deformations and long-time durations. The authors focus on a practical method using implicit time integration scheme applied to structural dynamic analyses in which the widely used Newmark time integration procedure is unstable, and not energy-momentum conserving. In this integration scheme, the time step is divided in two substeps. For too large time steps, the method is stable but shows excessive numerical dissipation. The influence of different substep sizes on the numerical dissipation of the method is studied throughout three practical examples. The method shows good performance and may be considered good for nonlinear transient response of structures.

scholarly journals Investigation of Air Pocket Behavior in Pipelines Using Rigid Column Model and Contributions of Time Integration Schemes

Water ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 785
Author(s):  
Arman Rokhzadi ◽  
Musandji Fuamba
Keyword(s):  
Transient Flow ◽  
Time Integration ◽  
Driving Pressure ◽  
Implicit Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Column Model ◽  
Integration Schemes ◽  
Backward Euler

This paper studies the air pressurization problem caused by a partially pressurized transient flow in a reservoir-pipe system. The purpose of this study is to analyze the performance of the rigid column model in predicting the attenuation of the air pressure distribution. In this regard, an analytic formula for the amplitude and frequency will be derived, in which the influential parameters, particularly, the driving pressure and the air and water lengths, on the damping can be seen. The direct effect of the driving pressure and inverse effect of the product of the air and water lengths on the damping will be numerically examined. In addition, these numerical observations will be examined by solving different test cases and by comparing to available experimental data to show that the rigid column model is able to predict the damping. However, due to simplified assumptions associated with the rigid column model, the energy dissipation, as well as the damping, is underestimated. In this regard, using the backward Euler implicit time integration scheme, instead of the classical fourth order explicit Runge–Kutta scheme, will be proposed so that the numerical dissipation of the backward Euler implicit scheme represents the physical dissipation. In addition, a formula will be derived to calculate the appropriate time step size, by which the dissipation of the heat transfer can be compensated.

An Efficient Weighted Residual Time Integration Family

2021 ◽  
pp. 2150106
Author(s):  
Mohammad Rezaiee-Pajand ◽  
S. A. H. Esfehani ◽  
H. Ehsanmanesh
Keyword(s):  
Degrees Of Freedom ◽  
Time Integration ◽  
Nonlinear Damping ◽  
Implicit Methods ◽  
New Family ◽  
Error Analyses ◽  
Damping Behavior ◽  
The Family ◽  
New Strategy ◽  
Time Integration Methods

A new family of time integration methods is formulated. The recommended technique is useful and robust for the loads with large variations and the systems with nonlinear damping behavior. It is also applicable for the structures with lots of degrees of freedom, and can handle general nonlinear dynamic systems. By comparing the presented scheme with the fourth-order Runge–Kutta and the Newmark algorithms, it is concluded that the new strategy is more stable. The authors’ formulations have good results on amplitude decay and dispersion error analyses. Moreover, the family orders of accuracy are [Formula: see text] and [Formula: see text] for even and odd values of [Formula: see text], respectively. Findings demonstrate the superiority of the new family compared to explicit and implicit methods and dissipative and non-dissipative algorithms.

Sign in / Sign up

Export Citation Format

Share Document