By analyzing the composite effects of the plastic shear strain and confining stress according to the laboratory test results, a model adopting a deteriorating elastic modulus for the intact rock is proposed. Soft rock such as coal and the strong rock such as granite with different quality are investigated. Geological Strength Index (GSI) is adopted to represent the integrity of rock, to be specific, from intact rock to the rock mass. Specifically, for the strain-softening rock mass, the elastic modulus is gained by introducing GSI into the deterioration model. The strength parameters are obtained by fitting GSI into Hoek-Brown failure criterion. Then, the elastic modulus and strength parameters of the rock mass are employed in the numerical procedure for solving out the rock deformation, the radii of the plastic softening and residual zones of the tunnels. The rationality of the numerical procedure is verified with the existing semi-analytical and analytical procedures. In the discussion, the influences of the critical plastic parameters and GSI, on the elastic modulus, rock deformation, radii of the plastic softening and residual zones are investigated. Five models for the variation of the modulus including the proposed deterioration model are defined. The rock mass deformation behaviours of the five models are compared. The results indicate that, in the plastic softening zone that is far away from the tunnel periphery, the elastic modulus is more sensitive to the plastic shear strain, whereas near the tunnel periphery, the elastic modulus is primarily affected by the confining stress. Regarding a linear decrease of the elastic modulus versus the plastic shear strain overestimates the elastic modulus to some extent, especially for the soft rock mass.