Search

Floor heave is a disaster that is frequently encountered in tunnel engineering. Existing explanations of floor heave mechanisms have many limitations, and the most important of which is that the stress concentration, release and transfer phenomena and the influences of the in-situ stress lateral coefficient and tensile strength are ignored. Therefore, the combined finite-discrete element method (FDEM) is used to study the failure process of floor rock masses and propose a new mechanical mechanism of floor heave that can consider the influences of the in-situ stress lateral coefficient and tensile strength. The type I fracture energy corresponding to different tensile strengths is calibrated using a direct tensile simulation test. Then, the floor heave mechanism is investigated under different lateral coefficients and tensile strengths, and five different floor heave modes are proposed under different in-situ stress lateral coefficients and tensile strengths. The floor heave modes controlled by the various in-situ stress lateral coefficients and tensile strengths are different, but they can all be explained by the maximum concentrated tangential stress. [ABSTRACT FROM AUTHOR]

A large amount of slates is distributed along Kangding No. 2 tunnel of Sichuan–Tibet Railway (China). The bedding of the Kangding slate is rich, showing apparent anisotropy of strength. For underground tunnels excavated in the anisotropic rock mass, the rock material is vulnerable to damage in tension because of its low tensile strength. Therefore, the tensile strength and fracture modes are important for the support design of layered surrounding-rock tunnel. In this study, first, a series of direct tensile tests were carried out on slate samples with five kinds of bedding dip angles. Second, the effect of bedding dip angle on the slate's tensile strength, failure modes, and acoustic emission characteristics was analyzed. Third, based on the Nova–Zaninetti criterion and the experiment results, a new phenomenological anisotropic criterion was proposed by introducing anisotropic coefficient, which could describe the nonlinear variation trend of tensile strength with foliation angles. Compared with four typical tensile strength criteria, the accuracy and universality of the new criterion were verified. Finally, the new criterion was implemented into the combined finite-discrete element method (FDEM) to conduct direct tensile simulation of slate with different foliation angles. The results indicated that the direct tensile strength followed an S-shaped increasing trend with the variation of foliation angles. The direct tensile failure modes of slate samples with different foliation angles could be divided into layer activation failure, mixed failure, and rock matrix failure. The acoustic emission process of slate with different bedding inclination under direct tension could be roughly divided into the following three stages: quiet period—stepped saltatory increase—drastic increase. The fitting effect of the new criterion was the highest among Kangding slate and other six typical layered rocks, which showed that the new criterion was suitable for predicting its Brazilian strength and direct tensile strength. At the same time, the numerical simulation results were in good agreement with the experimental results, which further verified the accuracy of the new criterion. Highlights: The acoustic emission characteristic parameters helped to reveal the direct tensile failure characteristics of layered rocks. A new phenomenological anisotropic tensile failure criterion was proposed. The direct tensile test data and the Brazilian splitting test data were used to study the accuracy and universality of the new anisotropic tensile failure criterion. FDEM was used to study the accuracy of the new anisotropic tensile failure criterion. [ABSTRACT FROM AUTHOR]

• A novel thermo-mechanical coupling scheme is developed based on the combined finite-discrete element method (FDEM) to capture the mechanical response of the rock to the isotropic/anisotropic thermal conduction process. • The mechanisms of the creation and development of the excavation damaged zone (EDZ) around the wellbore subject to excavation unloading, drilling mud support, and thermal change in the surrounding rock are analyzed. • The influences of the thermal-mechanical properties of the rock mass, such as the isotropic/anisotropic rock strength and thermal conductivity, thermal expansivity, as well as drilling mud-rock thermal interaction process, on the EDZ development, are discussed. Insights into fracture development and excavation damaged zone (EDZ) formation around wellbore excavation sites are essential for understanding wellbore stability. In this study, a novel thermo-mechanical scheme is implemented in the 2D combined finite-discrete element method (FDEM) to investigate the fracturing process around a wellbore in a high-temperature subsurface environment. The developed scheme captures the isotropic/anisotropic thermal conduction characteristics within the rock formation. The coupled scheme and stress distributions around the wellbore subject to in-situ stress, drilling mud pressure, and temperature changes are validated by closed-form solutions. Fracture development and progressive EDZ formation by various mechanisms are analyzed. The results show that for an unsupported wellbore, fractures initiate in the region of the most severe stress concentration due to excavation unloading and constitute an EDZ that resembles the logarithmic-spiral rupture zone captured by experimental observations and Mohr-Coulomb models, resulting in potential water inrush due to the significantly increased fracture transmissivity of the EDZ and wellbore collapse. The shape and extent of the EDZ for wellbore excavation in the layered rock mass are dominated by the stress redistribution and the presence of low-strength bedding planes favorably oriented for bedding slippage. The excavation unloading-induced stresses can be effectively counteracted by properly applying drilling mud pressure, and the disturbance of the stress field occurs only within a small region, greatly reducing the extension of the EDZ. Thermal contraction of the surrounding rock due to its convective interaction with the cold drilling mud induces extra fractures around the wellbore. The extension of the EDZ is greater with increased thermal expansivity of the rock formation. Overall, the results indicate that the adopted thermo-mechanical FDEM simulation can provide unique geomechanical insight into wellbore stability behavior analysis, in which explicit consideration of the fracturing and fragmentation processes is of great significance. [ABSTRACT FROM AUTHOR]

The large deformation mechanisms and predictions of soft rock tunnels have always been important but difficult to solve problems in the field of geotechnical engineering. A combined finite–discrete element numerical simulation method (FDEM) was used to study large deformation mechanism, classification and prediction. The deformations or failure mechanisms of tunnel surrounding rock were revealed, and the failure modes and displacement value prediction of surrounding rock with different strength-stress ratios were also investigated. The following conclusions were obtained: (1) Critical hysteresis damping can be adopted to simulate the progressive large deformation process of soft rock tunnels, which can obtain the final deformation of unsupported tunnels and avoid a dynamic response. (2) Under concentrated tangential stress, the surrounding rock undergoes X-shaped conjugate shear fracture, and this type of fracture network continues to propagate toward the depth of the surrounding rock until the model reaches a stable state; the reduction in tunnel cross-section is mainly caused by the macroscopic movement and volume expansion of rock fragments, and the latter is due to the generation of a large number of macroscopic voids. (3) As the strength-stress ratio decreases, the deformation or failure modes of the surrounding rock can be divided into four categories: elastic–plastic deformation, closed fracturing, shear dilation and broken expansion. (4) Finally, a prediction equation for isotropic and homogeneous soft rock unsupported tunnel deformation with general size driven by hydrostatic in situ stress is obtained, which indicates that with an increasing strength-stress ratio, the surrounding rock displacement decreases as an exponential function, with a correlation coefficient of R 2 = 0.997. In addition, the robustness and reliability of the prediction equation is verified. [ABSTRACT FROM AUTHOR]

• A new phenomenological tensile failure criterion with an anisotropic coefficient was proposed. • Direct tensile test of the Kangding slate was carried out to verify the accuracy of the new criterion. • The tensile strength data of eight typical layered rocks were collected and compared with the other six typical tensile strength criteria to verify the universality and superiority of the new criterion. • Based on the sensitivity analysis of tensile rate, the loading rate of 0.01 m/s was recommended for the samples with various foliation angles. • Layer thickness has little effect on tensile strength, but a significant effect on tensile failure modes. The tensile strength is a key parameter for the design of many geotechnical engineerings such as tunnel support, slope support, and petroleum drilling. The tensile strength in layered rocks, however, is anisotropic due to the existence of bedding planes. In this paper, to study the anisotropic tensile behaviors, a series of direct tensile tests were conducted using Kangding slate with five different foliation angles (β = 0°, 30°, 45°, 60°, and 90°). Based on the N-Z criterion, a new tensile failure criterion with the anisotropic coefficient was proposed. The tensile strength data of nine typical layered rocks were collected to evaluate the prediction ability of seven anisotropic tensile failure criteria including the new criterion. Embedding the new criterion into FDEM, the loading rate-dependent and layer thickness-dependent tensile mechanism of layered rock under direct tensile test were simulated. The experimental and numerical simulation results showed that the foliation angle significantly affected the tensile strength and failure modes of the slate. With the increase of foliation angle (β), the tensile strength of nine typical layered rocks, including Kangding slate, presented a non-linear increase trend. Compared with the other six typical anisotropic tensile failure criteria, the new criterion has the best prediction capability on the tensile strength of rocks with four different anisotropic degrees, and therefore the new criterion was recommended. The new criterion was embedded in FDEM to carry out direct tensile simulation of slate with different β. The numerical simulation results were in good agreement with the experimental results, further verifying the accuracy of the new criterion. The loading rate sensitivity analysis indicated that a tensile rate of 0.01 m/s was recommended to assure quasi-static loading of the specimen. Through the sensitivity analysis of layer thickness, it is concluded that layer thickness has very little influence on tensile strength, but has a very significant influence on the final failure modes of the specimens. [ABSTRACT FROM AUTHOR]

The combined finite-discrete element method (FDEM) is an advanced finite element and discrete element coupling method, which is commendably suitable for simulating the slope's entire evolution process from cracking, expansion, penetration, sliding, collision to deposition under seismic load. The original FDEM has apparent shortcomings in simulating the dynamic instability process of earthquake-induced landslides, and some targeted improvements are made to the original program. Numerical tests are carried out to prove the accuracy and robustness of the improved FDEM, and the influence of the structural plane on the dynamic evolution mechanism and process of the slope is studied. The results show that the improved FDEM can accurately reproduce the entire dynamic evolution processes and failure forms of different jointed rock slopes under seismic load and can quantitatively evaluate the dynamic stability of jointed rock slopes (dynamic safety factor and critical failure surface). In addition, the results also emphasize that the slope failure degree significantly affects the dynamic response of the slope and suggest that the combination of the cut-through of the failure surface and the increase of kinetic energy should be used as the criteria for slope dynamic instability simulation by FDEM. This study provides scientific and technological support for revealing the failure mechanism of jointed rock slopes and disaster reduction and prevention. [ABSTRACT FROM AUTHOR]

Rock materials are typically heterogeneous. Traditional digital image processing (DIP) technology and the grain-based model (GBM) method are difficult to apply to engineering-scale problem studies, such as tunnel excavation. Based on the combined finite-discrete element method (FDEM), a random parameter assignment method is proposed to simulate the mechanical properties and failure behavior of heterogeneous rock. The elastic moduli of triangular elements, as well as the strength parameters of quadrilateral joint elements, assigned by this method obey a Weibull distribution. The simulation results of uniaxial compression, triaxial compression and Brazilian disc show that with the increase in the heterogeneity m value, the uniaxial compressive strength, elastic modulus, triaxial compressive strength, equivalent cohesion, equivalent internal friction angle and tensile strength all increase exponentially and approach those of the homogeneous rock sample. In addition, the heterogeneous rock samples mainly shear fail along their corresponding theoretical failure angle under uniaxial and triaxial compression. The simulation results of tunnel excavation indicate that for heterogeneous rock with different m values, the overall fracture morphology of the surrounding rock remains unchanged, and X-shaped conjugate shear failures mainly occur followed by a small amount of tensile failures, which are similar to the homogeneous surrounding rock. However, with the increase in the m value, the failure degree and the maximum fracture propagation range of the surrounding rock decrease exponentially. Moreover, the simulation results of uniaxial compression and tunnel excavation with different element sizes suggest that the random parameter assignment method proposed in this paper is reliable. [ABSTRACT FROM AUTHOR]

To investigate the effects of the oval hole on the strength and the failure mechanisms of rock materials, a series of uniaxial and biaxial compression numerical simulations were carried out using the combined finite-discrete element method (FDEM). Based on the numerical simulation results, the effects of dip angle (β), major–minor ratio (ω) and position of the oval hole on rock strength characteristics crack initiation and development process and final failure modes were analysed in detail. The numerical simulation results indicate that: (1) Under uniaxial compression, the ultimate failure modes of rock containing oval holes can be mainly classified into single inclined plane shear failure, X-shaped tensile-shear composite failure and Y-shaped tensile-shear composite failure. Under biaxial compression, the ultimate failure mode is single-inclined plane shear failure dominated by shear cracks. Whether it is uniaxial compression or biaxial compression, the crack initiation position of the main fracture surface is near the end of the major axis of the oval hole. (2) Under uniaxial compression, the axial stress–strain curve of rock containing the oval hole is smooth before and after the peak, but under biaxial compression, there is obvious fluctuation in the post-peak stage due to crack initiation and propagation. In both uniaxial and biaxial compression, the pre-peak elastic stage of the axial stress–strain curve of the rock containing the oval hole is significantly shortened and the peak stress is significantly reduced compared with the intact rock, which indicates that the bearing capacity of the rock is seriously reduced by the existence of oval hole. (3) When the position and major–minor axis ratio of the oval hole remain unchanged, the uniaxial compressive strength and biaxial compressive strength increase with the increase of the dip angle of the oval hole. When the position and dip angle of the oval hole remain unchanged, the uniaxial compressive strength and biaxial compressive strength decrease with the increase of the major–minor axis ratio of the oval hole. (4) At different dip angles, the brittleness index Bi of rock containing oval holes decreases with the increase of confining pressure, and the brittleness degree shows a gradual weakening trend. At the same time, when the dip angle increases, the decreasing amplitude of brittleness index Bi decreases with the increase of confining pressure. Three types of failure modes of rock with elliptical hole are systematically summarised, and the failure mechanisms of elliptical hole rock are revealed. The influence of different positions, major–minor axis ratio, and dip angle of rock containing oval holes on uniaxial compression and biaxial compression are studied. The relationship between brittleness index and failure mode of rock containing the oval hole is investigated. [ABSTRACT FROM AUTHOR]

The combined finite-discrete element method (FDEM) is an advanced finite and discrete element coupling method, which can commendably simulate the large deformation and fracture process of materials. The explicit solution method in FDEM is naturally suitable for solving dynamic problems, and the key to dynamic analysis is to establish artificial boundary conditions that can be prepared to simulate infinite domain motion. At present, there is only the viscous boundary condition in the original FDEM, but this is only applicable for solving the internal source wave problems (e.g., blasting, mechanical vibration), and there are no artificial boundary conditions for solving the external source wave problems (e.g., earthquake). Therefore, three artificial boundary conditions have been developed in the present study to enhance the ability of FDEM to conduct dynamic response analysis. Firstly, the basic principle of FDEM and the existing viscous boundary condition are briefly introduced. Then, three newly added boundary conditions are introduced in sequence: (1) viscous-spring boundary condition, which can absorb stress wave energy at the boundary and restore the residual displacement to meet the actual engineering; (2) free-field boundary condition, which is coupled with free-field motion to absorb scattered waves at the lateral boundaries of the model, and the algorithms of mesh matching retrieval and coupling calculation are also introduced; (3) static-dynamic unified boundary condition, which can accurately transform the fixed boundary under the quasi-static analysis into the non-reflective boundary condition under the dynamic analysis. Finally, several classic models are built to verify the feasibility of the newly added boundary conditions. The comparison with the shaking table test also shows that the improved FDEM can be used for seismic response analysis under actual earthquakes. [ABSTRACT FROM AUTHOR]

The combined finite-discrete element method (FDEM) is an advanced finite and discrete element coupling method. This study further applies it to the problem of stress wave propagation in the jointed rock mass. Firstly, the fundamental theory of FDEM is briefly described and the methods commonly used in FDEM to characterize structural planes are comprehensively analyzed. Then, based on the viscous boundary condition (VBC), the viscous-spring boundary condition (VSBC) is added to absorb the reflected wave at the artificial boundary and restore the residual displacement to meet the actual engineering. In addition, the application range of the VBC and VSBC is verified, which indicates that the VBC and VSBC in the improved FDEM can be well applied not only to the continuum range but also to the case of the new fracture formed. Finally, several classic models are solved to verify the stress wave propagation across different types of joints by FDEM. The results reveal that results from FDEM agree well with analytical solutions based on the displacement discontinuity model and numerical solutions from universal distinct element code (UDEC), which means that the improved FDEM can capably and accurately simulate wave propagation in the jointed rock mass. [ABSTRACT FROM AUTHOR]

The failure mechanism and failure mode of anisotropic rock masses are very different from those of isotropic rock masses. For unsupported tunnels, the failure modes of anisotropic rock masses can be summarized as composite failure, V-shaped notch failure and spalling of bedding planes. Different failure modes are controlled by the mechanical parameters of the rock masses, the environmental occurrences and the excavation conditions. The combined finite-discrete element method (FDEM) was employed to investigate the failure mechanism and failure process of layered rock masses and the influence of different factors on the failure mode. The study results indicate that (1) composite failure is the basic failure mode of the layered rock masses and that it is caused by the combined action of shear-slip fractures F1 parallel to the bedding plane, tensile fractures F2 perpendicular to the bedding plane, and conjugate shear fractures F3; the rotation movement of slab-like rock fragments induced by fractures F1 and F3 into the tunnel causes significant deformation of the surrounding rock masses; (2) on the basis of composite failure, when fractures F1 disappear and fractures F2 obliquely intersect with the bedding plane, V-shaped notch failure occurs; furthermore, spalling of the bedding plane occurs when the propagation of fractures F3 is blocked at the bedding plane; and (3) with increasing rock mass strength, lateral pressure coefficient of in situ stress and layer thickness, or decreasing tunnel span, the surrounding rock masses change from composite failure to V-shaped notch failure and then to spalling of the bedding plane. However, the effect of the elastic modulus on the failure mode of anisotropic rock masses is weak, and the influence of the dip angle of the rock layer on the failure mode is complicated. [ABSTRACT FROM AUTHOR]

• A new constitutive model is proposed to eliminate the dependence of parameters on element size. • The influence laws of different parameters on the selection of fracture energy values are proposed, including E , ν , c , φ i , f t and h. • The reliability of the selection of fracture energy values is verified by uniaxial and triaxial compression simulations. The combined finite-discrete element method (FDEM) has been widely used in the numerical study of rock mechanics and geotechnical engineering. Selection of the correct parameter values is a prerequisite to ensure the accuracy of the simulation results. However, current calibration procedures are relatively cumbersome, and the parameters are strongly dependent on the element size. To eliminate the dependence of parameter values on the element size, this paper proposes a new constitutive model based on the stress–strain relationship. In addition, uniaxial compression and direct tension simulations are used to investigate the sensitivities of type II and type I fracture energies, i.e., G II and G I , including the influence of the G I (or G II) value on the simulation results of direct tension (or uniaxial compression) and the influence of different parameters (such as Young's modulus E , Poisson's ratio ν , cohesion c , internal friction angle φ i , tensile strength f t and element size h) on the selection of the values of G II and G I. The following results are obtained. (1) As G I or G II increases, the strength of the rock sample increases, and the plastic characteristics are more obvious, but the failure mode tends to be stable. (2) As the E value increases, the value of G I or G II decreases as a power function. As the value of f t or c increases, the value of G I or G II increases as a power function. However, other parameters, including ν , φ i and h , have no effect on the selection of the values of G II and G I. (3) Actual uniaxial compression and triaxial compression simulations prove that the calibrated G II and G I values are reasonable. In this study, only two fracture energies, G II and G I , must be calibrated; other macroparameters can be obtained through laboratory tests, and other microparameters can be obtained from previous empirical or theoretical values, which greatly improves the calibration efficiency of FDEM input parameters. [ABSTRACT FROM AUTHOR]

• A new calculation equation for normal contact stiffness is proposed. • The value of the coefficient in the new equation is proposed and the robustness of this value is verified. • The influence of different normal contact stiffness on uniaxial compression and tunnel excavation simulations are studied. The combined finite-discrete element method (FDEM) has been widely used in numerical studies in the fields of rock mechanics and geotechnical engineering. The normal contact stiffness between triangular elements is an important influencing parameter, but there is currently no effective method of measuring it. First, an equation for normal contact stiffness is proposed (P b = αP f , where P b is the basic stiffness, α is the coefficient, and P f is the joint penalty; the specific stiffness of each contact couple is determined by P b and the geometric sizes of the triangular elements together); then, the compression-shear failure of a single joint element is used to study the value range and robustness of α ; finally, uniaxial compression and tunnel excavation simulations are used to study the influence of different α values on rock crack propagation. The study results show that (1) an α value of 0.1448 is optimal and robust for a single joint element simulation and is therefore suitable for all simulation conditions; in other words, the value of α should not be too low to avoid decreasing the rock mass stiffness of existing natural fractures and deviating from the actual situation; in addition, the value of α should also not be too large to avoid the development of additional cracks, especially for hard rock simulation, in which the value of α is more sensitive; and (2) the crack topologies of the surrounding rock obtained by tunnel excavation simulation with a large α value deviate from the actual rock core (i.e., when α = 0.1448, the simulation results coincide well with the actual fractures). Through this study, the reliability of FDEM numerical simulation results is improved. [ABSTRACT FROM AUTHOR]

• The long-term degradation behavior of rock strength was systematically summarized. • The long-term degeneration laws of cohesion and internal friction angle were proposed. • FDEM was used to study the long-term crack propagation of tunnel surrounding rock. • The conclusion of this paper was verified by engineering case. The time-dependence of rock strength is an issue that cannot be ignored in the rock engineering field. This paper summarizes previous study results on the degradation of the uniaxial compressive strength of rock (including component cohesion and internal friction angle), tensile strength and Young's modulus over time. Based on previous studies, the degeneracy formulas of intact rock cohesion and internal friction angle are derived. Taking circular tunnel excavation under hydrostatic pressure as an example, the surrounding rock response is studied under different degradation degrees of rock strength and stiffness, including the radius of the excavation damage zone, damage degree of surrounding rock, lining support pressure and tunnel boundary displacement. The results show that with the increasing degree of degradation of the surrounding rock compressive strength, the above four parameters increase in the form of an exponential function and that weakening of the tensile strength had little effect on the surrounding rock response. In addition, as the degradation degree of rock stiffness increases, the radius of the excavation damage zone remains unchanged, whereas the other three parameters values are reduced. An actual engineering case simulation using the combined finite-discrete element method shows that results consistent with field observations can be obtained if the rock strength degradation law proposed in this paper is adopted. [ABSTRACT FROM AUTHOR]

Core material softening method has been widely used in two-dimensional tunnel excavation simulation by the combined finite-discrete element method (FDEM). In numerical studies, softening stress paths (generally described by critical kinetic energy (CKE), softening time step (STS) and softening curve (SC)) usually have a significant influence on simulation results and computational efficiency. There is not a direct way to determine these three parameters. In the present study, a series of numerical studies by the FDEM are performed to discuss the effects of these three parameters on the simulation results and calculation efficiency in a circular tunnel excavation project with a 4 m diameter. The results show that the CKE is closely related to the model size and rock density. The CKE can be obtained from the model kinetic energy curve, and a selection method is proposed. The STS value of 50 000 steps is optimal. To obtain a uniform and stable propagation process of the surrounding rock crack, the SC should be in an exponential function. The engineering case shows that the simulation results are consistent well with the field observations as the softening stress paths are calibrated by the proposed method. [ABSTRACT FROM AUTHOR]

Deformation and fracturing of the rock mass are heavily affected by complex thermal-hydraulic-mechanical (THM) processes in the subsurface environment. This study develops a THM coupling formulation based on the 2D combined finite-discrete element method (FDEM) to investigate the rock deformation and fracturing behavior driven by coupled THM phenomena in the rock mass. The coupled scheme consists of three subsolvers, i.e., a hydraulic solver, thermal solver, and mechanical solver. In the developed hydraulic solver, a flow network searching algorithm is proposed based on the pore network models and the special mesh discretization in the FDEM. The laminar viscous flow in the flow network is solved using the cubic law approximation, and the fluid pressure field is determined by a linear compressibility model. The developed thermal solver accounts for the complex thermal transport processes decomposed into conduction in the solid and the fluid, advection by fluid flow, and fluid-solid thermal convective transfer. The coupled THM formulation is implemented by iterating the three subsolvers using an explicit, partitioned scheme. This approach is then progressively verified by a series of benchmark problems. Finally, the deformation and branching of a single fracture embedded in the rock mass due to cold fluid injection are simulated, in which the gradual aseismic fracture slip/opening, as well as the branch fracture initiation and propagation, and their potential impact on seismic events are discussed. The results show the potential of the developed numerical tool in modeling the coupled THM process-driven deformation and fracturing of the rock mass. Improvements in the implementation are needed, such as improvements in modeling the thermo-poroelastic effects and temperature- and pressure-dependent physical properties. [ABSTRACT FROM AUTHOR]