Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 93
PROCEEDINGS OF THE TENTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY Edited by:
Paper 152
Implementation of a Direct Procedure for Critical Point Computations J. Mäkinen^{1}, R. Kouhia^{2}, A. Fedoroff^{2} and H. Marjamäki^{1}
^{1}Department of Mechanics and Design, Tampere University of Technology, Finland
, "Implementation of a Direct Procedure for Critical Point Computations", in , (Editors), "Proceedings of the Tenth International Conference on Computational Structures Technology", CivilComp Press, Stirlingshire, UK, Paper 152, 2010. doi:10.4203/ccp.93.152
Keywords: nonlinear eigenvalue problem, equilibrium equations, preconditioned iterations.
Summary
Computation of critical points on an equilibrium path requires solution of a nonlinear eigenvalue problem. When the external load is parameterized using a single parameter, the nonlinear stability eigenvalue problem consists of solving the equilibrium equations at the criticality condition. Several techniques exist for solution of such a system. Their algorithmic treatment is usually focused on direct linear solvers and thus the use of the block elimination strategy. This paper places special emphasis on a strategy which can also be used with iterative linear solvers. Due to the nonuniqueness of the critical eigenmode a normalizing condition is required. In addition, for bifurcation points, the Jacobian matrix of the augmented system is singular at the critical point and additional stabilization is required in order to maintain the quadratic convergence of the Newton's method. Depending on the normalizing condition, convergence to a critical point with a negative load parameter value can happen. The form of the normalizing equation is critically discussed and an alternative form is proposed which quarantees convergence to a positive critical load value.
If properly scaled, all the examined constraint equations yield a similar type of convergence of the extended system. When starting from a random eigenvector, one step of inverse power iteration considerably improved the robustness of the Newton's iteration. In the present study, the second derivatives of the residual vector for the geometrically exact threedimensional beam element have been computed from exact expressions. This approach guarantees the asymptotic quadratic convergence rate of the Newton's method. Numerical example computations of some spatial structures clearly demonstrates this. The proposed simple preconditioning strategy, which requires only a standard preconditioner for the stiffness matrix, seems to work well. However, a really high quality preconditioner is needed to obtain faster execution than using the block elimination strategy with a fast direct solver. This is mainly due to the small bandwidth of the shell structures discretized using the finite element method. Further studies will be focused on increasing the robustness and efficiency of the nonlinear eigenvalue procedure as well as developing a suitable preconditioning strategy. Also procedures enlarging the convergence domain will be developed. Furthermore, utilization of the information already present in the nonlinear eigenvalue algorithm to branch out to the possible secondary paths should be investigated.
purchase the fulltext of this paper (price £20)
go to the previous paper 
