• Two triangular shapes are used to formulate the iteration procedure required for nonlinear analysis. • The areas of these shapes are considered as two variables objective functions. • These functions are minimized with respect to each variable. • Two new constraint equations are obtained, which can be utilized for the nonlinear solver. • All of numerical experiences clearly demonstrate the merits of the authors' formulations. In this paper, two triangular shapes are used to formulate iteration procedures required for nonlinear analysis. These triangles are formed by the structural load–displacement curve. The areas of these shapes are considered as two variables objective functions. To find the optimal solutions, these functions are minimized with respect to each variable. As a result, two new constraint equations are obtained, which can be utilized as the nonlinear solver. To demonstrate the merit of authors’ scheme, some geometric nonlinear analyses of shells, frames and trusses are performed. Furthermore, present formulations are compared with the cylindrical arc-length method.