In this study two dimensional visco-plastic finite element model capable of handling the multi-step excavation was developed for investigating the effect of excavation support sequences on the behavior of underground openings in the jointed rock mass. First, the finite element model which is capable of handling the multi-step excavation is developed and verified. And then the model is combined with visco-plastic joint model. Ubiquitous joint pattern was considered in the model and joint properties in cach set were assumed to be indentical. Passive, full-grouted rockbolts were cosidered in the numerical model. The visco-plastic deformations of joints and rockbolts were assumed to be governed by Mohr-Conlomb and von Mises yield criteria, respectively. With the ability of removing elements, the model can simulate the multi-step excavation-suppport sequences. The reliability and applicability of the model to the stability analysis for the underground excavation in pratice was checked by simulating the behavior of underground crude oil storage caverns under construction.