Solution of Post-Buckling & Limit Load Problems,Without Inverting the Tangent Stiffness Matrix & WithoutUsing Arc-Length Methods

Abstract

In this study, the Scalar Homotopy Methods are applied to the solution of post-buckling and limit load problems of solids and structures, as exemplified by simple plane elastic frames, considering only geometrical nonlinearities. Explicit- ly derived tangent stiffness matrices and nodal forces of large-deformation planar beam elements, with two translational and one rotational degrees of freedom at each node, are adopted following the work of (Kondoh and Atluri (1986)). By us- ing the Scalar Homotopy Methods, the displacements of the equilibrium state are iteratively solved for, without inverting the Jacobian (tangent stiffness) matrix. It is well-known that, the simple Newton's method (and the Newton-Raphson iteration method that is widely used in nonlinear structural mechanics), which necessitates the inversion of the Jacobian matrix, fails to pass the limit load as the Jacobian matrix becomes singular. Although the so called arc-length method can resolve this problem by limiting both the incremental displacements and forces, it is quite complex for implementation. Moreover, inverting the Jacobian matrix generally consumes the majority of the computational burden especially for large-scale prob- lems. On the contrary, by using the presently developed Scalar Homotopy Meth- ods, convergence near limit loads, and in the post-buckling region, can be easily achieved, without inverting the tangent stiffness matrix and without using complex arc-length methods. The present paper thus opens a promising path for conduct- ing post-buckling and limit-load analyses of nonlinear structures. While the simple Williams' toggle is considered as an illustrative example in this paper, extension