3D DEM simulation of hard rock fracture in deep tunnel excavation induced by changes in principal stress magnitude and orientation

To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.

Numerical modeling of large deformation and nonlinear frictional contact of excavation boundary of deep soft rock tunnel

Roadways excavated in soft rocks at great depth are difficult to be maintained due to large deformation of surrounding rocks, which greatly influences the safety and efficiency of deep resources exploitation. During the excavation process of a deep soft rock tunnel, the rock wall may be compacted due to large deformation. In this paper, the technique to address this problem by a two-dimensional (2D) finite element software, large deformation engineering analyses software (LDEAS 1.0), is provided. By using the Lagrange multiplier method, the kinematic constraint of non-penetrating condition and static constraint of Coulomb friction are introduced to the governing equations in the form of incremental displacement. The numerical example demonstrates the efficiency of this technology. Deformations of a transportation tunnel in inclined soft rock strata at the depth of 1 000 m in Qishan coal mine and a tunnel excavated to three different depths are analyzed by two models, i.e. the additive decomposition model and polar decomposition model. It can be found that the deformation of the transportation tunnel is asymmetrical due to the inclination of rock strata. For extremely soft rock, large deformation can converge only for the additive decomposition model. The deformation of surrounding rocks increases with the increase in the tunnel depth for both models. At the same depth, the deformation calculated by the additive decomposition model is smaller than that by the polar decomposition model.

Zonal disintegration phenomenon in enclosing rock mass surrounding deep tunnels——Elasto-plastic analysis of stress field of enclosing rock mass

The zonal disintegration phenomenon （ZDP） is a typical phenomenon in deep block rock masses. In order to investigate the mechanism of ZDP, an improved non-linear Hock-Brown strength criterion and a bi-linear constitutive model of rock mass were used to analyze the elasto-plastic stress field of the enclosing rock mass around a deep round tunnel. The radius of the plastic region and stress of the enclosing rock mass were obtained by introducing dimensionless parameters of radial distance. The results show that tunneling in deep rock mass causes a maximum stress zone to appear in the vicinity of the boundary of the elastic and the plastic zone in the surrounding rock mass. Under the compression of a large tangential force and a small radial force, the rock mass in the maximum stress zone was in an approximate uniaxial loading state, which could lead to a split failure in the rock mass.

Zonal disintegration phenomenon in rock mass surrounding deep tunnels

Zonal disintegration is a typical static phenomenon of deep rock masses. It has been defined as alternating regions of fractured and relatively intact rock mass that appear around or in front of the working stope during excavation of a deep tunnel. Zonal disintegration phenomenon was successfully demonstrated in the laboratory with 3D tests on analogous gypsum models, two circular cracked zones were observed in the test. The linear Mohr-Coulomb yield criterion was used with a constitutive model that showed linear softening and ideal residual plastic to analyze the elasto-plastic field of the enclosing rock mass around a deep tunnel. The results show that tunneling causes a maximum stress zone to appear between an elastic and plastic zone in the surrounding rock. The zonal disintegration phenomenon is analyzed by considering the stress-strain state of the rock mass in the vicinity of the maximum stress zone. Creep instability failure of the rock due to the development of the plastic zone, and transfer of the maximum stress zone into the rock mass, are the cause of zonal disintegration. An analytical criterion for the critical depth at which zonal disintegration can occur is derived. This depth depends mainly on the character and stress concentration coefficient of the rock mass.

Analysis of mechanical behavior of soft rocks and stability control in deep tunnels

Due to the weakness in mechanical properties of chlorite schist and the high in situ stress in Jinping II hydropower station, the rock mass surrounding the diversion tunnels located in chlorite schist was observed with extremely large deformations. This may significantly increase the risk of tunnel instability during excavation. In order to assess the stability of the diversion tunnels laboratory tests were carried out in association with the petrophysical properties, mechanical behaviors and waterlweakening properties of chlorite schist. The continuous deformation of surrounding rock mass, the destruction of the support structure and a large-scale collapse induced by the weak chlorite schist and high in situ stress were analyzed. The distributions of compressive deformation in the excavation zone with large deformations were also studied. In this regard, two reinforcement schemes for the excavation of diversion tunnel bottom section were proposed accordingly. This study could offer theoretical basis for deed tunnel construction in similar geological condition~

Characterization of transient groundwater flow through a high arch dam foundation during reservoir impounding

Numerous deep underground projects have been designed and constructed in China, which are beyond the current specifications in terms of scale and construction difficulty. The severe failure problems induced by high in situ stress, such as rockburst, spalling, damage of deep surrounding rocks, and timedependent damage, were observed during construction of these projects. To address these problems, the dynamic design method for deep hard rock tunnels is proposed based on the disintegration process of surrounding rocks using associated dynamic control theories and technologies. Seven steps are basically employed：（i） determination of design objective,（ii） characteristics of site, rock mass and project, and identification of constraint conditions,（iii） selection or development of global design strategy,（iv）determination of modeling method and software,（v） preliminary design,（vi） comprehensive integrated method and dynamic feedback analysis, and（vii） final design. This dynamic method was applied to the construction of the headrace tunnels at Jinping II hydropower station. The key technical issues encountered during the construction of deep hard rock tunnels, such as in situ stress distribution along the tunnels, mechanical properties and constitutive model of deep hard rocks, determination of mechanical parameters of surrounding rocks, stability evaluation of surrounding rocks, and optimization design of rock support and lining, have been adequately addressed. The proposed method and its application can provide guidance for deep underground projects characterized with similar geological conditions.