• 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.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2009.03a
• /
• pp.1-10
• /
• 2009

Based on recent studies, the side resistance of rock socketed drilled shafts was affected by unconfined compressive strength of rock, socket roughness, rock types and joints, and initial normal stress. Especially, the socket roughness was affected by rock types and joints, drilling methods, and diameters of pile. In this study, a new roughness measurement system (BKS-LRPS, Backyoung-KyungSung Laser Roughness Profiling System) usable in water was developed. Based on the laboratory model tests, an EMD (Effective Measurement Distances) according to various turbidity was proposed as EMD=$1149.2{\times}T^{-0.64}_b$.

• Tunnel and Underground Space
• /
• v.17 no.3 s.68
• /
• pp.203-215
• /
• 2007

Stability and mechanical deformation behavior of rock masses are highly dependent on the mechanical characteristics of contained discontinuities. Therefore, mechanical characteristics of the discontinuities should be considered in the design of tunnel and underground structures. In this study, direct shear tests for rough granite joints were carried out under constant normal stiffness conditions. Effects of initial normal stress, shear velocity, and surface roughness on the characteristics of shear strength and deformation behaviors were examined. Results of shear testing under constant normal stiffness conditions reveal that shear behaviors could be classified into two categories, based on the amount of decrease in shear stress at the Int peak shear stress. With initial normal stiffness increasing, it turned out that shear displacement at peak stress and the first peak shear stress increased, however friction angle and friction coefficient showed decrease. In case of shear stiffness and average friction coefficient, it turned out that they are not dependent on the initial normal stress. Minor effects of shear velocity on rough joints were observed in several shear quantities. However, the effects of shear velocity were insignificant regardless of the normal stress increase. Change of shear strength and deformation characteristics on joint roughness were examined, however, it turned out that the variations were attributed to deviation of shear test specimens.

• Geomechanics and Engineering
• /
• v.24 no.3
• /
• pp.215-226
• /
• 2021

To study the evolution mechanism of cracks in rocks with multiple defects, rock-like samples with multiple defects, such as strip-shaped through-going cracks and cavity groups, are used, and the crack propagation law and changes in AE (acoustic emission) and strain of cavity groups under different inclination angles are studied. According to the test results, an increase in the cavity group inclination angle can facilitate the initial damage degree of the rock and weaken the crack initiation stress; the initial crack initiation direction is approximately 90°, and the extension angle is approximately 75~90° from the strip-shaped through-going cracks; thus, the relationship between crack development and cavity group initiation strengthens. The specific performance is as follows: when the initiation angle is 30°, the cracks between the cavities in the cavity group develop relatively independently along the parallel direction of the external load; when the angle is 75°, the cracks between the cavities in the cavity group can interpenetrate, and slip can occur along the inclination of the cavity group under the action of the shear mechanism rupture. With the increase in the inclination angle of the cavity group, the AE energy fluctuation frequency at the peak stress increases, and the stress drop is obvious. The larger the cavity group inclination angle is, the more obvious the energy accumulation and the more severe the rock damage; when the cavity group angle is 30° or 75°, the peak strain of the local area below the strip-shaped through-going fracture plane is approximately three times that when the cavity group angle is 45° and 60°, indicating that cracks are easily generated in the local area monitored by the strain gauge at this angle, and the further development of the cracks weakens the strength of the rock, thereby increasing the probability of major engineering quality damage. The research results will have important reference value for hazard prevention in underground engineering projects through rock with natural and artificial defects, including tunnels and air-raid shelters.

• Geomechanics and Engineering
• /
• v.14 no.5
• /
• pp.459-466
• /
• 2018

By conducting uniaxial loading cycle tests on the coal rock with outburst proneness, the dilatation characteristics at different loading rates were investigated. Under uniaxial loading and unloading, the lateral deformation of coal rock increased obviously before failure, leading to coal dilatation. Moreover, the post-unloading recovery of the lateral deformation was rather small, suggesting the onset of an accelerated failure. As the loading rate increased further, the ratio of the stress at the dilatation critical point to peak-intensity increased gradually, and the pre-peak volumetric deformation decreased with more severe post-peak damage. Based on the laboratory test results, the lateral deformation of the coals at different depths in the #1302 isolated coal pillars, Yangcheng Coal Mine, was monitored using wall rock displacement meter. The field monitoring result indicates that the coal lateral displacement went through various distinct stages: the lateral displacement of the coals at the depth of 2-6 m went through an "initial increase-stabilize-step up-plateau" series. When the coal wall of the working face was 24-18 m away from the measuring point, the coals in this region entered the accelerated failure stage; as the working face continued advancing, the lateral displacement of the coals at the depth over 6 m increased steadily, i.e., the coals in this region were in the stable failure stage.

• Tunnel and Underground Space
• /
• v.23 no.3
• /
• pp.203-211
• /
• 2013

Shear behavior of joint plane has been investigated considering the magnitude of normal stresses and initial surface roughness. Shear strength of joint plane has been measured by performing the multi-stage shear test in which applied normal stress level has been increased stepwise. Multi-stage shear test within the specified normal stress range has been repeated and two types of strength parameter variation have been observed: type 1 - both cohesion and friction angle decrease, type 2 - cohesion decrease and friction angle increase. Trends of strength parameter variation for the three rock types, gneiss, granite and shale, have been investigated and the influence of initial roughness of joint plane on the sequential shear strength change for the repeated multi-stage shear tests also has been analyzed.

• Explosives and Blasting
• /
• v.23 no.4
• /
• pp.31-44
• /
• 2005

In highly stressed conditions, the excavation damage zone induced by stress redistribution and disturbance must be evaluated after tunnel excavation. Therefore, the investigation of stress-induced deformation and fracture in rock is indispensable. In this study, fracture and damage mechanisms of rock induced by the accumulation of microcracks were investigated by the moving point regression technique as well as acoustic emission measured during uniaxial compression tests. Especially, the modified procedures to determine damage thresholds more systematically were newly proposed, and successfully applied to rock. From experiments, crack initiation and track damage stress levels were estimated to be $33{\~}36\%$ and $84{\~}89\%$ of uniaxial compressive strength respectively, for both of Hwangdeung granite and Yeosan marble. However, the normalized crack closure stress level for Yeosan marble was much higher than for Hwangdeung granite. In addition, the largest proportion of total axial strain in Hwangdeung granite was attributable to elastic deformation and initial microcracking. However, the greatest part of axial deformation in Yeosan marble arose from initial crack closure and unstable cracking. Finally, it was seen that unstable cracking after the crack damage stress level played a key part in the lateral deformation in rocks under uniaxial compression.

• Tunnel and Underground Space
• /
• v.19 no.6
• /
• pp.567-582
• /
• 2009

The world's largest nucleon decay experiment facility is constructed at a depth of approximately 1,000meters, in the Kamioka mine, Japan. Because of the character as a large cavern in deep underground, in-situ stress measurements were conducted to provide basic information for design of the cavern. Three overcoring methods were used: 8-element embedding gauges developed by Japanese Central Research Institute of Electric Power Industry, hemispherical ended borehole technique with eight strain cross-gauges, and Hollow Inclusion Cell with 12 strain gauges. The principle stresses were not perfectly similar in each measurement. The average values of the 6 stress element were used to provide the direction and the magnitude of three principle stress.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.244-251
• /
• 2003

A two-dimensional BEM code, $FRACOD^{2D}$, was applied to simulate fracture initiation and propagation processes in a rock pillar during an in situ heater test of a rock pillar planned at the $\"{A}sp\"{o}$ Underground Rock laboratory of SKB, in Southern Sweden. To take the advantage of conventional BEM for simulating fracturing processes, but without efforts for domain integral transformation, a hybrid approach is developed to simulate the fracturing processes in rock pillar under coupled thermo-mechanical loading. The code FRACOD was used for simulating the fracture initiation and propagation processes with its boundary tractions reflecting the effects of the initial and redistributed thermomechanical stresses in the domain of interest at multiple excavation and heating steps were produced by a special algorithm of stress inversion, based on resultant thermo-mechanical stress fields at each excavation and heat loading step by a FEM code without considering fracturing processes. This hybrid approach can take the advantages of both types of numerical methods and avoids their shortcomings for fracturing process simulation and domain effects, respectively. In this paper, we present the hybrid approach for the stress, displacements, and fracturing processes at sequential excavation and heating steps of the in situ heater test as a predictive modelling, the formulation of the fracturing models and the predictive results. Two sections of borehole depth, 0.5 m and 1.5 m below the tunnel floor are considered. The pillar area is modelled with the FRACOD and the stress field produced by excavation and heating is transferred with corresponding boundary stresses. From the modelling results, the degree of fracturing and damage are evaluated for 120 days of heating. Dominated shear fracturing in the vicinity of the central pillar was observed from the models at both sections, but spalled area appears to be limited. Based on the modelling results, a sensitivity study for the effect of pre-existing fractures in the vicinity of the holes is also conducted, and the initiation and evolution of EDZ around the deposition holes are investigated using this particular numerical technique.

• Geomechanics and Engineering
• /
• v.12 no.3
• /
• pp.347-360
• /
• 2017

The problem of cylindrical cavity expansion incorporating deformation dependent of intermediate principal stress in rock or soil mass is investigated in the paper. Assumptions that the initial axial total strain is a non-zero constant and the axial plastic strain is not zero are defined to obtain the numerical solution of strain which incorporates deformation-dependent intermediate principal stress. The numerical solution of plastic strains are achieved by the 3-D plastic potential functions based on the M-C and generalized H-B failure criteria, respectively. The intermediate principal stress is derived with the Hook's law and plastic strains. Solution of limited expansion pressure, stress and strain during cylindrical cavity expanding are given and the corresponding calculation approaches are also presented, which the axial stress and strain are incorporated. Validation of the proposed approach is conducted by the published results.