Observations on the effect of numerical dissipation on the nonlinear dynamic response of structural systems

Abstract Analysis of highly nonlinear systems subjected to dynamic loads involves continuously changing temporal solution characteristics. This is due to changes in the stiffness, and hence, the periods of vibration, with the consequence of changes in the extent of numerical dissipation (or algorithmic damping). The effect of modal contributions in the time-marching solution procedure is discussed in this paper through the implementation of an α-damping scheme within the frontal solution procedure employed in a generalized dynamic analysis program. As the theoretical background to the effectiveness of such a dissipative scheme is restricted to linear analysis, the interaction between the time integration step size and the finite element idealization errors in the nonlinear range is investigated by means of numerical examples. It is concluded that the use of a numerically dissipative scheme in conjunction with a carefully selected cutoff frequency is essential in acceleration response predictions, particularly where a highly nonlinear response is of interest. It is also shown that the displacement response is much less affected by the spurious higher modes of vibration and that the time required to achieve successful dissipation of the spurious higher mode response may involve the execution of a significant number of time-steps.

[1]  T. J.R. Hughes,et al.  ANALYSIS OF TRANSIENT ALGORITHMS WITH PARTICULAR REFERENCE TO STABILITY BEHAVIOR. , 1983 .

[2]  K. Bathe,et al.  Stability and accuracy analysis of direct integration methods , 1972 .

[3]  T. Hughes,et al.  An error analysis of truncated starting conditions in step‐by‐step time integration: Consequences for structural dynamics , 1987 .

[4]  Thomas J. R. Hughes,et al.  Improved numerical dissipation for time integration algorithms in structural dynamics , 1977 .

[5]  Edward L. Wilson,et al.  Nonlinear dynamic analysis of complex structures , 1972 .

[6]  P. J. Pahl,et al.  Development of an implicit method with numerical dissipation from a generalized ingle-step algorithm for structural dynamics , 1988 .

[7]  G. L. Goudreau,et al.  Evaluation of numerical integration methods in elastodynamics , 1973 .

[8]  P. J. Pahl,et al.  Practical performance of the θ 1 -method and comparison with other dissipative algorithms in structural dynamics , 1988 .

[9]  P. J. Pahl,et al.  Extended comparison of the Hilber-Hughes-Taylor α-method and the φ 1 -method , 1989 .

[10]  B. Izzuddin,et al.  Nonlinear dynamic analysis of framed structures , 1990 .

[11]  W. L. Wood,et al.  A unified set of single step algorithms. Part 1: General formulation and applications , 1984 .

[12]  G. Bazzi,et al.  The ρ‐family of algorithms for time‐step integration with improved numerical dissipation , 1982 .

[13]  W. L. Wood Practical Time-Stepping Schemes , 1990 .

[14]  John C. Houbolt,et al.  A Recurrence Matrix Solution for the Dynamic Response of Elastic Aircraft , 1950 .

[15]  T. Hughes,et al.  Collocation, dissipation and [overshoot] for time integration schemes in structural dynamics , 1978 .