Geometric nonlinear behavior of shell structures is traced using the combination algorithm of the load incrementation and the displacement incrementation. Particularly, generation scheme of the variable load / displacement increment sizes is added to improve the efficiency of this algorithm. The element utilized in this study is the 9-node degenerated shell element with five degrees of freedom and a total Lagrangian formulation is considered for the nonlinear finite element analysis. Various numerical examples are presented in order to demonstrate the effectiveness of this improved scheme in practical analysis.