• Explosives and Blasting
• /
• v.39 no.3
• /
• pp.15-23
• /
• 2021

Ventilation shafts are pathways in mines and tunnels for the removal of dust or smoke during underground space construction and operation. In mines, blasting with long blast holes is preferred for the excavation of a ventilation shaft in the 10~20m long crown pillar section. In this case, the bottom part of the blast hole is completely drilled in order to determine the drilling error, and this causes a problem of lowering the explosive charge and blasting efficiency. It is possible to solve the problem of explosive loading and to increase the blast efficiency by covering the curb of the blasthole by using stemming material. In this study, simulations for the blasting of a ventilation shaft were performed with various stemming lengths and the blasthole diameters(45, 76mm) using AUTODYN 2D SPH(Smooth particle hydrodynamics) analysis technique. Also the optimal bottom stemming column was derived by checking the size of the boulder and burden line according to blasting. Analysis result, blasting efficiency is lessened in case of stemming length less than 30cm and the optimal length of the stemming material should be 30cm or higher to achieve high efficiency of blasting.

• Tunnel and Underground Space
• /
• v.29 no.3
• /
• pp.172-183
• /
• 2019

The generalized Hoek-Brown (GHB) failure criterion developed by Hoek et al. (2002) is a nonlinear function which defines a stress condition at failure of rock mass. The relevant strength parameter values are systematically determined using the GSI value. Since GSI index is a value quantifying the condition of in-situ rock mass, the GHB criterion is a practical failure condition which can take into the consideration of in-situ rock mass quality. Considering that most rock mechanics engineers are familiar with the linear Mohr-Coulomb criterion and that many rock engineering softwares incorporate Mohr-Coulomb criterion, the equations for the equivalent friction angle and cohesion were also proposed along with the release of the GHB criterion. The proposed equations, however, fix the lower limit of the minor principal stress range, where the linear best-fitting is performed, with the tensile strength of the rock mass. Therefore, if the tensile stress is not expected in the domain of analysis, the calculated equivalent friction angle and cohesion based on the equations in Hoek et al. (2002) could be less accurate. In order to overcome this disadvantage of the existing equations for equivalent friction angle and cohesion, this study proposes the analytical formula which can calculate optimal equivalent friction angle and cohesion in any minor principal stress interval, and verified the accuracy of the derived formula.

• Tunnel and Underground Space
• /
• v.28 no.6
• /
• pp.635-650
• /
• 2018

Many people have been recognized that the Korean Peninsula is no longer safe area from the earthquake by the recent earthquakes occurred in the country. The earthquakes that occurred at Pohang and Gyeongju appeared differently from them considered in the seismic design and researches on the seismic design method have been also conducted by many researchers. Studies on seismic loads are mainly focused on existing superstructures, and research involving them has been actively carried out in reality. However, paper regarding structural stability of reinforcement from seismic load such as soil-nails, rock-bolts, ground anchors which were constructed to ensure stability of serviced structure have been published rarely. In this study, ground anchor been effected by static load and seismic load which is settled in the weathered rock is analyzed. Results for static load are obtained from field test and seismic load is from numerical analysis. In this study, the behavioral characteristics of the ground anchor were analyzed by numerical analysis in case of seismic loading based on the result of the in-situ tensile test of the ground anchor settled weathered rock. As a result, settlement of concrete block due to application of tension force for ground anchor occurred as well as following loss of axial force for ground anchor. Also, as bond length and period of seismic load are longer, increasement of displacement is greater.

• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.509-518
• /
• 2020

This study constructed a new simulation platform, including mesh generation process, numerical simulation, and post-processing for results analysis based on exploration data to perform real-scale numerical analysis considering the actual geological structure efficiently. To build the simulation platform, we applied for open-source programs. The source code is open to be available for code modification according to the researcher's needs and compatibility with various numerical simulation programs. First, a three-dimensional model(3D) is acquired based on the exploration data obtained using a drone. Then, the domain's mesh density was adjusted to an interpretable level using Blender, the free and open-source 3D creation suite. The next step is to create a 3D numerical model by creating a tetrahedral volume mesh inside the domain using Gmsh, a finite element mesh generation program. To use the mesh information obtained through Gmsh in a numerical simulation program, a converting process to conform to the program's mesh creation protocol is required. We applied a Python code for the procedure. After we completed the stability analysis, we have created various visualization of the study using ParaView, another open-source visualization and data analysis program. We successfully performed a preliminary stability analysis on the full-scale Dokdo model based on drone-acquired data to confirm the usefulness of the proposed platform. The proposed simulation platform in this study can be of various analysis processes in future research.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.480-493
• /
• 2021

Since the generalized Hoek-Brown criterion (GHB) provides an efficient way of identifying its strength parameter values with the consideration of in-situ rock mass condition via Geological Strength Index (GSI), this criterion is recognized as one of the standard rock mass failure criteria in rock mechanics community. However, the nonlinear form of the GHB criterion makes its mathematical treatment inconvenient and limits the scope of its application. As an effort to overcome this disadvantage of the GHB criterion, the explicit approximate analytical equations for the minimum principal stress, which is associated with the maximum principal stress at failure, are formulated based on the Taylor polynomial approximation of the original GHB criterion. The accuracy of the derived approximate formula for the minimum principal stress is verified by comparing the resulting approximate minimum principal stress with the numerically calculated exact values. To provide an application example of the approximate formulation, the equivalent friction angle and cohesion for the expected plastic zone around a circular tunnel in a GHB rock mass are calculated by incorporating the formula for the approximate minimum principal stress. It is found that the simultaneous consideration of the values of mi, GSI and far-field stress is important for the accurate calculation of equivalent Mohr-Coulomb parameter values of the plastic zone.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.494-507
• /
• 2021

With the increasing use of TBM, research has recently been conducted in Korea to analyze TBM data with machine learning techniques to predict the ground in front of TBM, predict the exchange cycle of disk cutters, and predict the advance rate of TBM. In this study, classification prediction of rock characteristics of slurry shield TBM sites was made by combining traditional rock classification techniques and machine learning techniques widely used in various fields with machine data during TBM excavation. The items of rock characteristic classification criteria were set as RQD, uniaxial compression strength, and elastic wave speed, and the rock conditions for each item were classified into three classes: class 0 (good), 1 (normal), and 2 (poor), and machine learning was performed on six class algorithms. As a result, the ensemble model showed good performance, and the LigthtGBM model, which showed excellent results in learning speed as well as learning performance, was found to be optimal in the target site ground. Using the classification model for the three rock characteristics set in this study, it is believed that it will be possible to provide rock conditions for sections where ground information is not provided, which will help during excavation work.

• Tunnel and Underground Space
• /
• v.32 no.5
• /
• pp.312-326
• /
• 2022

The generalized Hoek-Brown (GHB) criterion, which is recognized as one of the standard failure conditions for rock mass, is specialized for rock engineering applications and covers a wide range of rock mass conditions. Accordingly, many research efforts have been devoted to the incorporation of this criterion into the stability analysis of rock structures. In this study, the slip-line analysis method, which is a kind of elastoplastic analysis method, is combined with the GHB failure criterion to derive analytical equations that can easily calculate the plastic radius and stress distribution in the vicinity of the circular tunnel. In the process of derivation of related formulas, it is assumed that the behavior of rock mass after failure is perfectly plastic and the in-situ stress condition is hydrostatic. In the formulation, it is revealed that the plastic radius can be calculated analytically using the two respective tangential friction angles corresponding to the stress conditions at tunnel wall and elastic-plastic boundary. It is also shown that the plastic radius and stress distribution calculated using the derived analytical equations coincide with the results of Lee & Pietruszczak's numerical method published in 2008. In the latter part of this paper, the influence of the quality of the rock mass on the size of the plastic zone, the stress distribution, and the change of the tangential friction angle was investigated using the derived analytical equations.

• Tunnel and Underground Space
• /
• v.33 no.3
• /
• pp.150-168
• /
• 2023

The anisotropic behavior of rocks is primarily attributed to the directional arrangement of rock-forming minerals and the distribution characteristics of microcracks. Notably, sedimentary and metamorphic rocks often exhibit distinct transverse isotropy in terms of their strength and deformation characteristics. Consequently, it is crucial to gain accurate insights into the deformation and failure characteristics of transversely isotropic rocks during rock mechanics design processes. The deformation of such rocks is described by five independent elastic constants, which are determined through laboratory testing. In this study, the characteristics of the directional variation of apparent elastic constants in transversely isotropic rocks were investigated using experimental data reported in the literature. To achieve this, the constitutive equation proposed by Mehrabadi & Cowin was introduced to calculate the apparent elastic constants more efficiently and systematically in a rotated Cartesian coordinate system. Four transversely isotropic rock types from the literature were selected, and the influence of changes in the orientation of the weak plane on the variations of the apparent elastic modulus, apparent shear modulus, and apparent Poisson's ratio was analyzed. Based on the investigation, a new constraint on the elastic constants has been proposed. If the proposed constraint is satisfied, the directional variation of the apparent elastic constants in transversely isotropic rocks aligns with intuitive predictions of their tendencies.

• 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.33 no.6
• /
• pp.435-444
• /
• 2023

In the licensing application for the deep geological disposal system of high-level radioactive waste in Finland, the disposal tunnel backfilling has been changed from the block/pellet (for the construction) to the granular type (for the operation). Accordingly, for establishing the design concept for backfilling, it is necessary to examine applicability to the domestic facility through analyzing problems of the existing method and improvements in the alternative design. In this paper, we first reviewed the principal studies conducted for changing the backfill method in the licensing process of the Finnish facility, and identified the expected problems in applying the block/pellet backfill method. In addition, we derived the evaluation factors to be considered in terms of technical and operational aspects for the backfilling solution, and then conducted a comparative analysis for two types of backfill methods. This analysis confirmed the overall superiority of the design change. It is expected that these results could be utilized as the technical basis for deriving the optimum design plan in development process of the Korean-specific deep disposal facility. However, applicability should be reviewed in advance based on the latest technical data for the detailed evaluation factors that must be considered for selecting the backfilling method.