Stabilized forms of orthogonal residual and constant incremental work control path following methods

Stabilized versions of the orthogonal residual and constant incremental work control path following methods are introduced. These new methods overcome the difficulty observed in the original methods when traversing the equilibrium path near displacement limit points. The modification automatically relaxes the orthogonality constraint near the displacement limit points. In addition, the modified methods can ideally be used in conjunction with iterative linear solvers.