To study the tunnel stability at various static water pressures and determine the mechanical properties and deformation behavior of surrounding rock, a modified effective stress formula was introduced into a numerical integration algorithm of elastic-plastic constitutive equation, that is, closest point projection method (CPPM). Taking the effects of water pressure and seepage into account, a CPPM-based formula was derived and a CPPM algorithm based on Drucker-Prager yield criterion considering the effect of pore water pressure was provided. On this basis, a CPPM-based elastic-plastic numerical analysis program considering pore water pressure was developed, which can be applied in the engineering of tunnels and other underground structures. The algorithm can accurately take the effects of groundwater on stability of surrounding rock mass into account and it can show the more pronounced effect of pore water pressure on stress, deformation, and the plastic zone in a tunnel. The stability of water flooding in Fusong tunnel was systematically analyzed using the developed program. The analysis results showed that the existence of groundwater seepage under tunnel construction will give rise to stress redistribution in the surrounding rock mass. Pore water pressure has a significant effect on the surrounding rock mass.
Finite element framework for closest point projection method for surface PDEs
