Blasting induced dynamic stress concentration and failure characteristics of deep-buried rock tunnel

In this study,the dynamic stress concentration factors(DSCF)around a straight-wall arch tunnel(SWAT)were solved analytically utilizing the complex variable function methods and Duhamel’s integral.The effects of wavelength,incident angle,and blasting rising time on the DSCF distribution were analyzed.Theoretical results pointed out dynamic disturbances resulting in compressive stress concentration in the vertical direction and tensile stress in the incident direction.As the wavelength and rising time increased,there was a tendency for the amplitude of stress concentration to initially rise and then converge.Moreover,a series of 3D FEM models were established to evaluate the effect of different initial stress states on the dynamic failure of the tunnel surrounding rock.The results indicated that the failure of the surrounding rock was significantly influenced by the direction of the static maximum principal stress and the direction of the dynamic disturbance.Under the coupling of static and blasting loading,damage around the tunnel was more prone to occur in the dynamic and static stress concentration coincidence zone.Finally,the damage modes of rock tunnel under static stress and blasting disturbance from different directions were summarized and a proposed support system was presented.The results reveal the mechanisms of deep-buried rock tunnel destruction and dynamically triggered rockburst.

The role of polyurethane foam compressible layer in the mechanical behaviour of multi-layer yielding supports for deep soft rock tunnels

The polyurethane foam(PU)compressible layer is a viable solution to the problem of damage to the secondary lining in squeezing tunnels.Nevertheless,the mechanical behaviour of the multi-layer yielding supports has not been thoroughly investigated.To fill this gap,large-scale model tests were conducted in this study.The synergistic load-bearing mechanics were analyzed using the convergenceconfinement method.Two types of multi-layer yielding supports with different thicknesses(2.5 cm,3.75 cm and 5 cm)of PU compressible layers were investigated respectively.Digital image correlation(DIC)analysis and acoustic emission(AE)techniques were used for detecting the deformation fields and damage evolution of the multi-layer yielding supports in real-time.Results indicated that the loaddisplacement relationship of the multi-layer yielding supports could be divided into the crack initiation,crack propagation,strain-hardening,and failure stages.Compared with those of the stiff support,the toughness,deformability and ultimate load of the yielding supports were increased by an average of 225%,61%and 32%,respectively.Additionally,the PU compressible layer is positioned between two primary linings to allow the yielding support to have greater mechanical properties.The analysis of the synergistic bearing effect suggested that the thickness of PU compressible layer and its location significantly affect the mechanical properties of the yielding supports.The use of yielding supports with a compressible layer positioned between the primary and secondary linings is recommended to mitigate the effects of high geo-stress in squeezing tunnels.

Design Theory of Full Face Rock Tunnel Boring Machine Transition Cutter Edge Angle and Its Application

At present, the inner cutters of a full face rock tunnel boring machine (TBM) and transition cutter edge angles are designed on the basis of indentation test or linear grooving test. The inner and outer edge angles of disc cutters are characterized as symmetric to each other with respect to the cutter edge plane. This design has some practical defects, such as severe eccentric wear and tipping, etc. In this paper, the current design theory of disc cutter edge angle is analyzed, and the characteristics of the rock-breaking movement of disc cutters are studied. The researching results show that the rotational motion of disc cutters with the cutterhead gives rise to the difference between the interactions of inner rock and outer rock with the contact area of disc cutters, with shearing and extrusion on the inner rock and attrition on the outer rock. The wear of disc cutters at the contact area is unbalanced, among which the wear in the largest normal stress area is most apparent. Therefore, a three-dimensional model theory of rock breaking and an edge angle design theory of transition disc cutter are proposed to overcome the flaws of the currently used TBM cutter heads, such as short life span, camber wearing, tipping. And a corresponding equation is established. With reference to a specific construction case, the edge angle of the transition disc cutter has been designed based on the theory. The application of TBM in some practical project proves that the theory has obvious advantages in enhancing disc cutter life, decreasing replacement frequency, and making economic benefits. The proposed research provides a theoretical basis for the design of TBM three-dimensional disc cutters whose rock-breaking operation time can be effectively increased.

Wear Analysis of Disc Cutters of Full Face Rock Tunnel Boring Machine

Wear is a major factor of disc cutters’ failure. No current theory offers a standard for the prediction of disc cutter wear yet. In the field the wear prediction method commonly used is based on the excavation length of tunnel boring machine（TBM） to predict the disc cutter wear and its wear law, considering the location number of each disc cutter on the cutterhead（radius for installation）; in theory, there is a prediction method of using arc wear coefficient. However, the preceding two methods have their own errors, with their accuracy being 40% or so and largely relying on the technicians’ experience. Therefore, radial wear coefficient, axial wear coefficient and trajectory wear coefficient are defined on the basis of the operating characteristics of TBM. With reference to the installation and characteristics of disc cutters, those coefficients are modified according to penetration, which gives rise to the presentation of comprehensive axial wear coefficient, comprehensive radial wear coefficient and comprehensive trajectory wear coefficient. Calculation and determination of wear coefficients are made with consideration of data from a segment of TBM project（excavation length 173 m）. The resulting wear coefficient values, after modification, are adopted to predict the disc cutter wear in the follow-up segment of the TBM project（excavation length of 5621 m）. The prediction results show that the disc cutter wear predicted with comprehensive radial wear coefficient and comprehensive trajectory wear coefficient are not only accurate（accuracy 16.12%） but also highly congruous, whereas there is a larger deviation in the prediction with comprehensive axial wear coefficient（accuracy 41%, which is in agreement with the prediction of disc cutters’ life in the field）. This paper puts forth a new method concerning prediction of life span and wear of TBM disc cutters as well as timing for replacing disc cutters.

Development of a new type of foam concrete and its application on stability analysis of large-span soft rock tunnel

The long-term stability of large-span soft rock tunnel is influenced greatly by the creep effect of surrounding rock.The development of a new type of foam concrete which has the property of high compressibility and low ductility was introduced.And it was made as filling material of reserved deformation layer between the first lining and the second lining used in large-span soft rock tunnel.The effect of the new type of foam concrete was simulated as filling material of reserved deformation layer using numerical simulation.Through the comparison with the common large-span soft rock tunnel,the vault settlement and surrounding convergence are reduced by about 61% and 45%,respectively,after creep of 100 a.And in the second lining,the plastic zone reduces apparently and the maximum equivalent plastic strain decreases relatively.So,it can be found that the application of the new type of foam concrete as the filling material of reserved deformation layer can relieve the excessive force in second lining induced by rock creep,reduce its deformation and improve the stability of tunnel.

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.

Updating the models and uncertainty of mechanical parameters for rock tunnels using Bayesian inference

Rock mechanical parameters and their uncertainties are critical to rock stability analysis,engineering design,and safe construction in rock mechanics and engineering.The back analysis is widely adopted in rock engineering to determine the mechanical parameters of the surrounding rock mass,but this does not consider the uncertainty.This problem is addressed here by the proposed approach by developing a system of Bayesian inferences for updating mechanical parameters and their statistical properties using monitored field data,then integrating the monitored data,prior knowledge of geotechnical parameters,and a mechanical model of a rock tunnel using Markov chain Monte Carlo(MCMC)simulation.The proposed approach is illustrated by a circular tunnel with an analytical solution,which was then applied to an experimental tunnel in Goupitan Hydropower Station,China.The mechanical properties and strength parameters of the surrounding rock mass were modeled as random variables.The displacement was predicted with the aid of the parameters updated by Bayesian inferences and agreed closely with monitored displacements.It indicates that Bayesian inferences combined the monitored data into the tunnel model to update its parameters dynamically.Further study indicated that the performance of Bayesian inferences is improved greatly by regularly supplementing field monitoring data.Bayesian inference is a significant and new approach for determining the mechanical parameters of the surrounding rock mass in a tunnel model and contributes to safe construction in rock engineering.

Tunnelling performance prediction of cantilever boring machine in sedimentary hard-rock tunnel using deep belief network

Evaluating the adaptability of cantilever boring machine(CBM) through in-depth excavation and analysis of tunnel excavation data and rock mass parameters is the premise of mechanical design and efficient excavation in the field of underground space engineering.This paper presented a case study of tunnelling performance prediction method of CBM in sedimentary hard-rock tunnel of Karst landform type by using tunneling data and surrounding rock parameters.The uniaxial compressive strength(UCS),rock integrity factor(Kv),basic quality index([BQ]),rock quality index RQD,brazilian tensile strength(BTS) and brittleness index(BI) were introduced to construct a performance prediction database based on the hard-rock tunnel of Guiyang Metro Line 1 and Line 3,and then established the performance prediction model of cantilever boring machine.Then the deep belief network(DBN) was introduced into the performance prediction model,and the reliability of performance prediction model was verified by combining with engineering data.The study showed that the influence degree of surrounding rock parameters on the tunneling performance of the cantilever boring machine is UCS > [BQ] > BTS >RQD > Kv > BI.The performance prediction model shows that the instantaneous cutting rate(ICR) has a good correlation with the surrounding rock parameters,and the predicting model accuracy is related to the reliability of construction data.The prediction of limestone and dolomite sections of Line 3 based on the DBN performance prediction model shows that the measured ICR and predicted ICR is consistent and the built performance prediction model is reliable.The research results have theoretical reference significance for the applicability analysis and mechanical selection of cantilever boring machine for hard rock tunnel.

Optimization of construction scheme and supporting technology for HJS soft rock tunnel

For a soft rock tunnel under high stress in jointed and swell soft rock （HJS）, two construction schemes pilot-tunneling enlarging excavation and step-by-step excavation were optimized using FLAC20, and the deformation effects of the two construction schemes were verified by field tests. Based on engineer- ing geological investigation and mechanical analysis of large deformations, the complex deformation mechanisms of stress expansion and structural deformation of the soft rock tunnel were confirmed, and support countermeasures from the complex deformation mechanism converted to a single type were proposed, and the support parameters were optimized by field tests. These technologies were proved by engineering practice, which produced significant technical and economic benefits.

Study on the Deformation Mechanism of a Soft Rock Tunnel

The large deformation of soft rock tunnel is one of the key problems to be overcome in the tunnel construction stage.In the present study,the deformation mechanism of a representative tunnel and some related countermeasures are investigated using field tests and engineering geological analysis.Owing to the scarce performances of methods based on other criteria such as small pipe spacing,anchor bolt length and steel frame spacing,a new support scheme is implemented and optimized.Results show that shear failure and bedding sliding are produced under high horizontal stress conditions.The low strength of the surrounding rock results in the uneven convergence of both sides of the tunnel.With the aforementioned new support scheme,however,most of such problems can be mitigated leading to good stability properties and ensuing economic advantages.

Study on the new shell bolting and shotcrete support and its application in soft rock tunnels

A new kind of tunnel support was put forward on the basis of the anchor spraysupport principle. The mechanics of the new three-dimensional steel bar shotcrete liningsupport was studied and the structure's internal forces were analyzed. The model experiment was done relying on the industrial test. The conclusion of numerical calculationsproved that the ANSYS program is reasonable and creditable. It was compared to otherkinds of support that are commonly used in soft rock tunnels. The technique and economiccontrasts of the typical tunnel with support three-dimensional steel bar were completed.

Heterogeneity identification method for surrounding rock of large-section rock tunnel faces based on support vector machine

Purpose–The purpose of the study is to quickly identify significant heterogeneity of surrounding rock of tunnel face that generally occurs during the construction of large-section rock tunnels of high-speed railways.Design/methodology/approach–Relying on the support vector machine(SVM)-based classification model,the nominal classification of blastholes and nominal zoning and classification terms were used to demonstrate the heterogeneity identification method for the surrounding rock of tunnel face,and the identification calculation was carried out for the five test tunnels.Then,the suggestions for local optimization of the support structures of large-section rock tunnels were put forward.Findings–The results show that compared with the two classification models based on neural networks,the SVM-based classification model has a higher classification accuracy when the sample size is small,and the average accuracy can reach 87.9%.After the samples are replaced,the SVM-based classification model can still reach the same accuracy,whose generalization ability is stronger.Originality/value–By applying the identification method described in this paper,the significant heterogeneity characteristics of the surrounding rock in the process of two times of blasting were identified,and the identification results are basically consistent with the actual situation of the tunnel face at the end of blasting,and can provide a basis for local optimization of support parameters.

Mechanical properties and associated seismic isolation effects of foamed concrete layer in rock tunnel

Foamed concrete has a good energy absorption capability and can be used as seismic isolation material for tunnels. This study aims to investigate the mechanical properties and associated seismic isolation effects of foamed concrete layer in rock tunnel. For this, a series of uniaxial/triaxial compression tests was conducted to understand the effects of concrete density, confining stress and strain rate on the mechanical properties of foamed concrete. The direct shear tests were also performed to investigate the effects of concrete density and normal stress on the nonlinear behaviors of foamed concrete layer-lining interface. The test results showed that the mechanical properties of foamed concrete are significantly influenced by the concrete density. The foamed concrete also has high volumetric compressibility and strain-rate dependence. The peak stress. residual stress. shear stiffness and residual friction coefficient of the foamed concrete layer-lining interface are influenced by the foamed concrete density and normal stress applied. Then, a crushable foam constitutive model was constructed using ABAQUS software and a composite exponential model was also established to study the relationship between shear stress and shear displacement of the interface, in which their parameters were fitted based on the experimental results. Finally, a parametric analysis using the finite element method(FEM) was conducted to understand the influence of foamed concrete layer properties on the seismic isolation effect, including the density and thickness of the layer as well as the shear stiffness and residual friction coefficient of the interface. It was revealed that lower density and greater thickness in addition to smaller shear stiffness or residual friction coefficient of the foamed concrete layer could yield better seismic isolation effect, and the influences of the first two tend to be more significant.

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.

Tunnel failure in hard rock with multiple weak planes due to excavation unloading of in-situ stress

Natural geological structures in rock(e.g.,joints,weakness planes,defects)play a vital role in the stability of tunnels and underground operations during construction.We investigated the failure characteristics of a deep circular tunnel in a rock mass with multiple weakness planes using a 2D combined finite element method/discrete element method(FEM/DEM).Conventional triaxial compression tests were performed on typical hard rock(marble)specimens under a range of confinement stress conditions to validate the rationale and accuracy of the proposed numerical approach.Parametric analysis was subsequently conducted to investigate the influence of inclination angle,and length on the crack propagation behavior,failure mode,energy evolution,and displacement distribution of the surrounding rock.The results show that the inclination angle strongly affects tunnel stability,and the failure intensity and damage range increase with increasing inclination angle and then decrease.The dynamic disasters are more likely with increasing weak plane length.Shearing and sliding along multiple weak planes are also consistently accompanied by kinetic energy fluctuations and surges after unloading,which implies a potentially violent dynamic response around a deeply-buried tunnel.Interactions between slabbing and shearing near the excavation boundaries are also discussed.The results presented here provide important insight into deep tunnel failure in hard rock influenced by both unloading disturbance and tectonic activation.

A novel image-based approach for interactive characterization of rock fracture spacing in a tunnel face

This paper presents a novel integrated method for interactive characterization of fracture spacing in rock tunnel sections.The main procedure includes four steps:(1)Automatic extraction of fracture traces,(2)digitization of trace maps,(3)disconnection and grouping of traces,and(4)interactive measurement of fracture set spacing,total spacing,and surface rock quality designation(S-RQD)value.To evaluate the performance of the proposed method,sample images were obtained by employing a photogrammetrybased scheme in tunnel faces.Experiments were then conducted to determine the optimal parameter values(i.e.distance threshold,angle threshold,and number of fracture trace grouping)for characterizing rock fracture spacing.By applying the identified optimal parameters involved in the model,the proposed method could lead to excellent qualitative results to a new tunnel face.To perform a quantitative analysis,three methods(i.e.field,straightening,and the proposed method)were employed in the same study and comparisons were made.The proposed method agrees well with the field measurement in terms of the maximum and average values of measured spacing distribution.Overall,the proposed method has reasonably good accuracy and interactive advantage for estimating the ultimate fracture spacing and S-RQD.It can be a possible extension of existing methods for fracture spacing characterization for two-dimensional(2D)rock tunnel faces.

Face stability analysis of circular tunnels in layered rock masses using the upper bound theorem

An analysis of tunnel face stability generally assumes a single homogeneous rock mass.However,most rock tunnel projects are excavated in stratified rock masses.This paper presents a two-dimensional(2D)analytical model for estimating the face stability of a rock tunnel in the presence of rock mass stratification.The model uses the kinematical limit analysis approach combined with the block calculation technique.A virtual support force is applied to the tunnel face,and then solved using an optimization method based on the upper limit theorem of limit analysis and the nonlinear Hoek-Brown yield criterion.Several design charts are provided to analyze the effects of rock layer thickness on tunnel face stability,tunnel diameter,the arrangement sequence of weak and strong rock layers,and the variation in rock layer parameters at different positions.The results indicate that the thickness of the rock layer,tunnel diameter,and arrangement sequence of weak and strong rock layers significantly affect the tunnel face stability.Variations in the parameters of the lower layer of the tunnel face have a greater effect on tunnel stability than those of the upper layer.

Time-dependent squeezing deformation mechanism of tunnels in layered soft-rock stratum under high geo-stress

Large squeezing deformation of layered soft rock tunnel under high geo-stress has a significant time-dependent deformation behavior.In this paper,we studied the deformation mechanism during the construction period of deep-buried softrock tunnel by means of a combination of field observations and a numerical method.First,a new classification criterion for large deformations based on the power exponent variation law between the deformation and the strength-stress ratio is proposed.Then,the initial damage tensor reflecting the bedding plane(joint)distribution and an equivalent damage evolution equation derived from the viscoplastic strain are introduced based on the geometric research method,i.e.,a new rheological damage model(RDL model)of layered soft rock is established consisting of elastic,viscous,viscoelastic,viscoplastic and plastic elements.A field test was conducted on the Maoxian tunnel in Sichuan province,southwestern China,which is in broken phyllite(layered soft rock)under high geo-stress.The tunnel has experienced large deformation due to serious squeezing pressure,thus we adopted double primary support method to overcome the supporting structure failure problems.The rheological parameters of phyllite in the Maoxian tunnel were recognized by using SA-PSO optimization,and the RDL model does a good job in describing the time-dependent deformation behavior of a layered soft-rock tunnel under high geo-stress.Thus,the RDL model was used to investigate the supporting effect and bearing mechanism of the double primary support method.Compared with the single primary support method,the surrounding rock pressure,secondary lining force,surrounding rock deformation,and the depth of the damage to the rock mass was reduced by 40%-60%after the double primary support method was used.

Automatic identification of rock discontinuity and stability analysis of tunnel rock blocks using terrestrial laser scanning

Local geometric information and discontinuity features are key aspects of the analysis of the evolution and failure mechanisms of unstable rock blocks in rock tunnels.This study demonstrates the integration of terrestrial laser scanning(TLS)with distinct element method for rock mass characterization and stability analysis in tunnels.TLS records detailed geometric information of the surrounding rock mass by scanning and collecting the positions of millions of rock surface points without contact.By conducting a fuzzy K-means method,a discontinuity automatic identification algorithm was developed,and a method for obtaining the geometric parameters of discontinuities was proposed.This method permits the user to visually identify each discontinuity and acquire its spatial distribution features(e.g.occurrences,spac-ings,trace lengths)in great detail.Compared with hand mapping in conventional geotechnical surveys,the geometric information of discontinuities obtained by this approach is more accurate and the iden-tification is more efficient.Then,a discrete fracture network with the same statistical characteristics as the actual discontinuities was generated with the distinct element method,and a representative nu-merical model of the jointed surrounding rock mass was established.By means of numerical simulation,potential unstable rock blocks were assessed,and failure mechanisms were analyzed.This method was applied to detection and assessment of unstable rock blocks in the spillway and sand flushing tunnel of the Hongshiyan hydropower project after a collapse.The results show that the noncontact detection of blocks was more labor-saving with lower safety risks compared with manual surveys,and the stability assessment was more reliable since the numerical model built by this method was more consistent with the distribution characteristics of actual joints.This study can provide a reference for geological survey and unstable rock block hazard mitigation in tunnels subjected to complex geology and active rockfalls.

Failure responses of rock tunnel faces during excavation through the fault-fracture zone

It is essential to cast light on the construction risks in tunnel excavations through the fault-fracture zone(FFZ).This study adopts the material point method(MPM)to simulate the failure responses of a rock tunnel face during excavation through the FFZ.A numerical study was conducted to compare a physical model test and validate the feasibility of using the MPM in simulating tunnel face failure.One hundred ninety numerical simulation cases were constructed to represent a rock tunnel excavation project with different site con-figurations.The simulation results suggest that the cohesion and the friction angle significantly influence failure responses.The tunnel cover depth can magnify the failure responses,and the FFZ thickness significantly affects the mobilized rock mass volume when the FFZ consists of a weak rock mass.The numerical simulation results suggest three deformation patterns:face bulge,partial failure,and slide collapse.The failure responses can be characterized by stress arch,slip surface,angle of reposing,and influence range.The insights suggested by the face failure responses during excavation through the FFZ can aid field engineers in determining the scope of possible damage,and in establishing emergency measures to minimize losses if such failure occurs.