• Tunnel and Underground Space
• /
• v.25 no.5
• /
• pp.450-461
• /
• 2015

Extensive numerical experiments have been carried out to investigate effect of block size and fracture geometry on hydraulic characteristics of fractured rock masses based on connected pipe flow in DFN systems. Using two fracture sets, a total of 72 2-D fracture configurations were generated with different combinations of fracture size distribution and deterministic fracture density. The directional block conductivity including the theoretical block conductivity, principal conductivity tensor and average block conductivity for each generated fracture network system were calculated using the 2-D equivalent pipe network method. There exist significant effects of block size, orientation, density and size of fractures in a fractured rock mass on its hydraulic behavior. We have been further verified that it is more difficult to reach the REV size for the fluid flow network with decreasing intersection angle of two fracture sets, fracture plane density and fracture size distribution.

• Tunnel and Underground Space
• /
• v.31 no.1
• /
• pp.66-81
• /
• 2021

The effects of directional fracture tensor components and first invariant of fracture tensor on deformation moduli and shear moduli of fractured rock masses is analyzed based on regression analysis performed between 3-D fracture tensor parameters and deformability of DFN blocks. Using one or two deterministic joint sets, a total of 224 3-D discrete fracture network (DFN) cube blocks were generated with various configurations of deterministic density and probabilistic size distribution. The fracture tensor parameters were calculated for each generated DFN systems. Also, deformability moduli with respect to three perpendicular direction of the DFN cube blocks were estimated based on distinct element method. The larger the first invariant of fracture tensor, the smaller the values for the deformability moduli of the DFN blocks. These deformability properties present an asymptotic pattern above the certain threshold. It is found that power-law function describes the relationship between the directional deformability moduli and the corresponding fracture tensor components estimated in same direction.

• Tunnel and Underground Space
• /
• v.8 no.4
• /
• pp.342-350
• /
• 1998

In this study, fracture network in rock mass was interpreted using borehole wall images obtained by televiewer. The orientation and JRC value of major joint set were evaluated adopting image analysis techniques, of which process were written in macro-program code. As linking JRC to joint stiffness using Barton-Bandis model, fracture network map was produced for application to jointed rock modelling in numerical analysis of underground structure.

• Tunnel and Underground Space
• /
• v.19 no.2
• /
• pp.95-106
• /
• 2009

Numerical analysis of the progressive failure of a rock slope was conducted using a 3-D rock joint element considering fracture mechanics and subcritical crack growth of asperities in the rock joints. Even though the stress state in the rock slope is not changing, the elapse of time causes subcritical crack growth to break asperities in the joints. The increase of broken asperities causes failure of joints in the rock slope and the increase of failed joints results in failure of a jointed rock slope. As a result, the progressive failure of a jointed rock slope due to the gradual breaking of small asperities along joints generated by subcritical crack growth occurs at a lower stress than if rock failure occurred by exceeding the static strength or fracture toughness.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.270-288
• /
• 2021

We proposed a numerical method to simulate the hydro-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) in the paper. As a part of DECOVALEX-2023 Task G, we verified the method via benchmarks with analytical solutions. DECOVALEX-2023 Task G aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as a group of tetrahedral grains and calculated the interaction of the grains and their interfaces using 3DEC. The micro-parameters of the grains and interfaces were determined by a new methodology based on an equivalent continuum approach. In benchmark modeling, a single fracture embedded in the rock was examined for the effects of fracture inclination and roughness, the boundary stress condition and the applied pressure. The simulation results showed that the developed numerical model reasonably reproduced the fracture slip induced by boundary stress condition, the fracture opening induced by fluid injection, the stress distribution variation with fracture inclination, and the fracture roughness effect. In addition, the fracture displacements associated with the opening and slip showed good agreement with the analytical solutions. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Geomechanics and Engineering
• /
• v.17 no.6
• /
• pp.597-606
• /
• 2019

Ultrasonic Time-of-Flight Diffraction (TOFD) techniques are useful methods for non-destructive evaluation of fracture characteristics. This study focuses on the reliability and accuracy of ultrasonic diffraction methods to estimate the depth of rock fractures. The study material includes three different rock types; andesite, basalt and ignimbrite. Four different ultrasonic techniques were performed on these intact rocks. Artificial near-surface fracture depths were created in the laboratory by sawing. The reliability and accuracy of each technique was assessed by comparison of the repeated measurements at different path lengths along the rock surface. The standard error associated with the predictive equations is very small and their reliability and accuracy seem to be high enough to be utilized in estimating the depth of rock fractures. The performances of these techniques were re-evaluated after filling the artificial fractures with another material to simulate natural infills.

• Geomechanics and Engineering
• /
• v.11 no.2
• /
• pp.253-267
• /
• 2016

The stability of surrounding rock will be poor when the tunnel is excavated through nearly horizontal stratum. In this paper, the instability mechanism of local nearly horizontal stratum in super-large section and deep buried tunnel is revealed by the analysis of the macro failure and micro fracture. A structural model is proposed to explain the mechanics of surrounding rock collapse under the action of stress redistribution and shed light on the macroscopic analytical approach of the stability of surrounding rock. Then, some highly effective formulas applied in the tunnel engineering are developed according to the theory of mixed-mode micro fracture. And well-documented field case is made to demonstrate the effectiveness and accuracy of the proposed analytical methods of mixed-mode fracture. Meanwhile, in order to make the more accurate judgment about yield failure of rock mass, a series of comprehensive failure criteria are formed. In addition, the relationship between the nonlinear failure criterion and $K_I$ and $K_{II}$ of micro fracture is established to make the surrounding rock failure criterion more comprehensive and accurate. Further, the influence of the parameters related to the tension-shear mixed-mode fracture and compression-shear mixed-mode fracture on the propagation of rock crack is analyzed. Results show that ${\sigma}_3$ changes linearly with the change of ${\sigma}_1$. And the change rate is related to ${\beta}$, angle between the cracks and ${\sigma}_1$. The proposed simple analytical approach is economical and efficient, and suitable for the analysis of local nearly horizontal stratum in super-large section and deep buried tunnel.

• Tunnel and Underground Space
• /
• v.33 no.3
• /
• pp.189-207
• /
• 2023

In the present study, the thermoshearing experiment on a rough rock fracture were modeled using a three-dimensional grain-based distinct element model (GBDEM). The experiment was conducted by the Korea Institute of Construction Technology to investigate the progressive shear failure of fracture under the influence of thermal stress in a critical stress state. The numerical model employs an assembly of multiple polyhedral grains and their interfaces to represent the rock sample, and calculates the coupled thermo-mechanical behavior of the grains (blocks) and the interfaces (contacts) using 3DEC, a DEM code. The primary focus was on simulating the temperature evolution, generation of thermal stress, and shear and normal displacements of the fracture. Two fracture models, namely the mated fracture model and the unmated fracture model, were constructed based on the degree of surface matedness, and their respective behaviors were compared and analyzed. By leveraging the advantage of the DEM, the contact area between the fracture surfaces was continuously monitored during the simulation, enabling an examination of its influence on shear behavior. The numerical results demonstrated distinct differences depending on the degree of the surface matedness at the initial stage. In the mated fracture model, where the surfaces were in almost full contact, the characteristic stages of peak stress and residual stress commonly observed in shear behavior of natural rock joints were reasonably replicated, despite exhibiting discrepancies with the experimental results. The analysis of contact area variation over time confirmed that our numerical model effectively simulated the abrupt normal dilation and shear slip, stress softening phenomenon, and transition to the residual state that occur during the peak stress stage. The unmated fracture model, which closely resembled the experimental specimen, showed qualitative agreement with the experimental observations, including heat transfer characteristics, the progressive shear failure process induced by heating, and the increase in thermal stress. However, there were some mismatches between the numerical and experimental results regarding the onset of fracture slip and the magnitudes of fracture stress and displacement. This research was conducted as part of DECOVALEX-2023 Task G, and we expect the numerical model to be enhanced through continued collaboration with other research teams and validated in further studies.

• Tunnel and Underground Space
• /
• v.19 no.4
• /
• pp.366-372
• /
• 2009

Recently, in order to achieve smooth fracture plane and minimize the excavation damage zone in rock blasting, controlled blasting methods which utilize new technologies such as electronic delay detonator (EDD) and a notched charge hole have been suggested. In this study, smooth blastings utilizing three wing type notched charge holes are simulated to investigate the influence of explosive initial density on the resultant fracture plane and damage zone using dynamic fracture process analysis (DFPA) code. Finally, based on the dynamic fracture process analyses, novel smooth blasting method, ED-Notch SB (Electronic Detonator Notched Charge Hole Smooth Blasting) is suggested.

• The Journal of Engineering Geology
• /
• v.7 no.1
• /
• pp.53-62
• /
• 1997

Fracture mechanics properties of rock materials can he applied to predict the distribution of natural fractures in rock masses, and also to assess the safety of rock slopes and underground structures. In this study, rock fracture toughness and other fracture rrechanics properties of sorne lithologies showing apparently rock-property-controlled distribution of natural fractures were measured. Propagation behaviors of natural and experirnental fractures were also characterized both qualitatively and quantitatively, in terns of the propagation types and sorne statistical parameters. It was concluded that the application of fracture mechanics theories to the ge6logic materials should be based on the geological background and evidences.