In underground construction the multi-step excavation sequence is commonly adopted for the convenience of the underground work. A numerical simulation method which is capable of analyzing the effects of excavation sequence on the stability of the opening is greatly needed. In this study a two dimensional finite element code was developed based on the effective numerical algorithm for the multistep excavation. The practical applicability of the model was verified for the simplified excavation sequences.