This paper deals with a spherical cavity excavated in infinite homogeneous and isotropic strain-softening rock mass subjected to a hydrostatic initial Displacements and stresses distributions are two of the most important information for underground engineering, for that it can provide a theoretical foundation for geotechnical engineering optimism design and stability evaluation. In past few years, many scholars have developed series of (semi-)analytical and numerical solutions for circular openings and spherical cavities. . However, most of those solutions are based on the common elasto-perfectly-plastic and elasto-brittle-plastic models, because those solutions can be easily obtained using simple mathematical methods. However, most of the geotechnical materials display strain-softening behavior in the post-failure region. Therefore, the normally used elasto-perfectly plastic model and elasto-brittleplastic model can't well represent the actual failure process. Brown et al. were the earliest ones to analyze the stresses and displacements of circular openings excavated in the strain-softening rock mass. But, the assumption of constant elastic strain in the plastic region, which should obey the generalized Hooke's law, leads to a poor accuracy of the results. And Wang points out that Brown et al.'s solution is difficult to predict the plastic radius [11]. Carranza-Torres and Fairhurst's solution for H-B rock mass, which is based on the so-called self-similarity of H-B criterion, is theoretically rigorous. But, it seems rather complicated for practical use and is applicable only to the As we see, all of the solutions for strain-softening rock mass are based on semi-analytical and numerical method, no analytical solutions are available. Additionally, most of the indoor and field tests show that material properties of rock mass change during the failure process, not only strength parameters but also elastic parameters, i.e. Yong's modulus and Poisson's ratio. This is the so-called elasto-plastic coupling effect. Generally speaking, the effect of Young's modulus is analyzed in two ways. One is the pressure-dependent Young's modulus model (PDM) [26][27], the progressive failure process may be more significant. However, the two models can hardly be solved with analytical method. Considering the different tensile and compressive Young's modulus, Luo et al. obtained the closed-form solutions of spherical cavities in elastobrittle-plastic rock mass [28]. The difficulty of analytical methods for strain-softening rock mass is the coupling effect between softening index and material properties. In view of this, the Multi-step Brittle-Plastic model (MBPM) is proposed to solve the stress, strain, and displacement in the post-failure region for the spherical cavity. Moreover, the deterioration process of elastic parameters and strength parameters are also considered.