We propose a new numerical analysis method which combines beam-spring model and discrete element method. This method enables stability analysis of tunnel lining with large deformation or change of analysis domain. Analysis of a box culvert is solved as an example of this method by using the program which is based on the formulation explained in this paper.
