Nonlinear dynamic analysis of quasi-symmetric anisotropic structures

Abstract An efficient computational strategy is presented for the nonlinear dynamic analysis of quasi-symmetric anisotropic structures. A mixed formulation is used with the fundamental unknowns consisting of stress resultants, generalized displacements and velocity components. The governing semi-discrete finite element equations consist of a mixed system of algebraic and ordinary differential equations. The temporal integration of the differential equations is performed by using an explicit half-station central difference method. The three key elements of the strategy are: 1. (a) use of mixed finite element models with independent shape functions for the stress resultants, generalized displacements and velocity components and with the stress resultants allowed to be discontinuous at interelement boundaries 2. (b) operator splitting, or additive decomposition of the different arrays in the governing equations into the contributions to a symmetrized response plus correction terms 3. (c) application of a preconditioned conjugate gradient technique to generate the unsymmetric response of the structure, at each time step, as the sum of symmetric and antisymmetric modes, each obtained using approximately half the degrees of freedom of the original model. The preconditioning matrix is taken to be the matrix associated with the symmetrized response. The effectiveness of the proposed strategy is demonstrated by means of a numerical example and the potential of the proposed strategy for solving more complex nonlinear problems is discussed.

[1]  O. Zienkiewicz,et al.  A note on mass lumping and related processes in the finite element method , 1976 .

[2]  Isaac Fried,et al.  Finite element mass matrix lumping by numerical integration with no convergence rate loss , 1975 .

[3]  Loyce M. Adams,et al.  An M-Step preconditioned Conjugate Gradient Method for Parallel Computation , 1983, ICPP.

[4]  Ahmed K. Noor,et al.  Exploiting symmetries for efficient post-buckling analysis of composite plates , 1977 .

[5]  Samuel W. Key,et al.  Transient shell response by numerical time integration , 1973 .

[6]  Ahmed K. Noor,et al.  Preconditioned conjugate gradient technique for the analysis of symmetric anisotropic structures , 1987 .

[7]  A. Greenbaum Comparison of splittings used with the conjugate gradient algorithm , 1979 .

[8]  Ahmed K. Noor,et al.  EXPLOITING SYMMETRIES FOR EFFICIENT POST-BUCKLING ANALYSIS OF COMPOSITE PLATES , 1976 .

[9]  T. Belytschko,et al.  Computational Methods for Transient Analysis , 1985 .

[10]  D. Kershaw The incomplete Cholesky—conjugate gradient method for the iterative solution of systems of linear equations , 1978 .

[11]  T. Belytschko,et al.  A Précis of Developments in Computational Methods for Transient Analysis , 1983 .

[12]  C. Chia Nonlinear analysis of plates , 1980 .

[13]  Ahmed K. Noor,et al.  Mixed models and reduced/selective integration displacement models for nonlinear shell analysis , 1982 .

[14]  Ahmed K. Noor,et al.  Symmetry considerations for anisotropic shells , 1976 .

[15]  A. Jennings,et al.  The solution of sparse linear equations by the conjugate gradient method , 1978 .

[16]  J. Whitney Analysis of anisotropic rectangular plates. , 1972 .