Numerical implementation of pseudo-arclength algorithm in nonlinear finite element analysis

A numerical procedure to calculate the pre-buckling and postbuckling response of general structures is presented. This procedure is based on the pseudo-arclength algorithm suggested by E. Riks et al., which has some numerical difficulties during implementation of large applied analysis programs. To overcome these difficulties, a scheme based on rank-1 modification of the matrix is proposed. Some example show this procedure behaves well in passing through the limit point and is rather efficient. recommended by Prof. Chen Yaosong and Prof. Wu Jike