Numerical Simulation of Deformation and Failure Process of Tunnel Surrounding Rock Based on Different Constitutive Models

Based on five different constitutive models, the progressive failure process of deep hard rock was simulated, and the plastic zone, shear strain increment high-value zone, horizontal displacement, and maximum principal stress distribution of surrounding rock were obtained. The calculation results show that the plastic zone of surrounding rock calculated by MCSS and HBSS models is mainly concentrated in the floor and vault of the cavern, and the sidewall is not damaged. There is no shear zone in the surrounding rock, and the magnitude of horizontal displacement is 10^-4. In the L-CWFS model, the hysteresis cohesion of friction strengthening is weakened, and the initial friction angle is 0. The calculated plastic zone runs through the surrounding rock, and the whole section is destroyed. The order of magnitude of the horizontal displacement extremum is 10^-3. The MC-HB model simulates the spalling failure of the side wall, and the strain localization of the side wall is not obvious, and no shear band is formed inside the surrounding rock deeply. The N-CWFS-TSS model simulates the buckling failure process of the straight wall plate crack. A shear band is formed at the bottom of the straight wall and penetrates into the interior of the surrounding rock. The strain localization is obvious and the order of magnitude of the shear strain increment is 10^-2. The calculation results are consistent with the results of laboratory tests and field observations, and the failure mechanism of the compression-shear-tensile-shear composite in the progressive failure process of the surrounding rock is clarified.