Calculating the distribution of stresses and displacements around a circular tunnel excavated in infinite isotropic rock mass subjected to hydrostatic stress condition is one of the basic problems in rock engineering. While closed-form solutions for the distribution are known if rock masses are considered as elastic, perfectly plastic, or brittle-plastic media, a few numerically approximated solutions based on various simplifying assumptions have been reported for strain-softening rock mass. In this study, a simple numerical method is introduced for the analysis of strain-softening behavior of the circular tunnel in Mohr-Coulomb rock mass. The method can also applied to the analysis of the tunnel in brittle-plastic or perfectly plastic media. For the brittle-plastic case where closed-formsolution exists, the performance of the present method is verified by showing an excellent agreement between two solutions. In order to demonstrate the strain-softening behaviors predicted by the proposed method. a parameter study for a softening index is given and the construction of ground reaction curves is carried out. The importance of defining the characteristics of dilation in plastic analysis is discussed through analyzing the displacements near the surface of tunnel.