• Tunnel and Underground Space
• /
• v.30 no.1
• /
• pp.39-62
• /
• 2020

In this study, Barcelona Basic Model (BBM) was implemented into TOUGH2-MP/FLAC3D for the numerical analysis of coupled thermo-hydro-mechanical (THM) behavior of unsaturated soils and the prediction of long-term behaviors. Similar to the methodology described in a previous study for the implementation of BBM into TOUGH-FLAC, the User Defined Model (UDM) of FLAC based on the Modified Cam Clay Model (MCCM) and the FISH function of FLAC3D were used to extend the existing MCCM module in FLAC3D for the implementation of BBM into TOUGH2-MP/FLAC3D. In the developed BBM module in TOUGH2-MP/FLAC3D, the plastic strains due to change in suction increase (SI) in addition to mean effective stress are calculated. In addition to loading-collapse (LC) yield surface, suction increase (SI) yield surface is changed by hardening rules in the developed BBM module. Several numerical simulations were conducted to verify and validate the implementation of BBM: using an example presented in the FLAC3D manual for the standard MCCM, simulation results using COMSOL, and experimental data presented in SKB Reports. In addition, the developed BBM analysis module was validated by simultaneously performing a series of modeling tests that were performed for the validation of the Quick tools developed for the purpose of effectively deriving BBM parameters, and by comparing the Quick tools and Code_Bright results reported in a previous study.

• Tunnel and Underground Space
• /
• v.30 no.2
• /
• pp.164-179
• /
• 2020

Applying the limit state design method to geotechnical structures, accuracy and reliability of its design are mainly affected by parameters for geotechnical site characteristics, such as unit weight, Poisson's ratio, deformation modulus, cohesion and frictional angle. When the structures are located in weathered ground, especially, cohesion and frictional angle of ground are closely related with decision of parameters for structures' load and ground's resistance. Therefore, the accurate determination of these parameters, which are commonly obtained from field measurement, such as borehole shear test, are essential for optimum design of geotechnical structures. The 38 case studies, in this study, have been analyzed for understanding the importance of these parameters in designing the ground structures. From these results, importance of field measurement was also ascertained. With these evaluations, an apparatus for determining the strength characteristics, which are fundamental in limit state design (LSD) method, have been newly developed. This apparatus has an improved function as following the ASTM suggestion. Through the field application of this apparatus, the strong point of minimizing the possibility of error occurrence during the measurement has been verified and authors summarized that the essential parameters for LSD can be qualitatively obtained by this apparatus for determination of strength characteristics of weathered ground.

• Tunnel and Underground Space
• /
• v.26 no.3
• /
• pp.235-252
• /
• 2016

This study presents a methodology to reproduce the mechanical behavior of isotropic or transversely isotropic rock using the polygonal grain-based distinct element model. A numerical technique to monitor the evolution of micro-cracks during the simulation was developed in the present study, which enabled us to examine the contribution of tensile cracking and shear cracking to the progressive process of the failure. The numerical results demonstrated good agreement with general observations from rock specimens in terms of the behavior and the evolution of micro-cracks, suggesting the capability of the model to represent the mechanical behavior of rock. We also carried out a parametric study as a fundamental work to examine the relationships between the microscopic properties of the constituents and the macroscopic behavior of the model. Depending on the micro-properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics. In addition, a numerical technique to reproduce the transversely isotropic rock was suggested and applied to Asan gneiss from Korea. The behavior of the numerical model was in good agreement with the results obtained in the laboratory-scale experiments of the rock.

• Tunnel and Underground Space
• /
• v.29 no.5
• /
• pp.292-304
• /
• 2019

Limestone joint surfaces with smooth roughness were experimented by means of both the individual direct shear tests based on the KSRM standard test method and the multi-stage direct shear test to apply the stepwise vertical stresses. Changes in the roughness of the joint surfaces before and after the shear tests were examined and the difference between the two kinds of tests mentioned above was analyzed. In both tests, the shear resistance increased as the joint roughness increased and the maximum shear stress required for shearing the joint surface increased as the vertical stress increased. The peak friction angle obtained by the multi-stage direct shear tests was only 63% of that obtained by the individual direct shear tests. In the multi-stage direct shear test, the initial engagement of the concave-convex parts changes frequently during stepwise shearing process, which deforms the original roughness of a joint surface. Accordingly, the individual direct shear test is thought to be more effective when obtaining the friction angle of the rock joint surfaces. Limestone joint surfaces with smooth roughness of JRC value 4~8 were found to have peak friction angle of $47^{\circ}$, residual friction angle of $38^{\circ}$ and cohesion of 37 kPa.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.610-622
• /
• 2021

The DECOVALEX project is one of the representative international cooperative projects to enhance the understanding of the complex Thermo-Hydro-Mechanical-Chemical(THMC) coupled behavior in the high-level radioactive waste disposal system based on the numerical simulation. DECOVALEX-2023 is the current phase consisting of 7 tasks, and Task C aims to model the THM coupled behavior in the disposal system based on the Full-scale Emplacement (FE) experiment at the Mont-Terri underground rock laboratory. This study performs the numerical simulation based on the OGS-FLAC developed for the current study. In the numerical model, we emplaced the heater with constant power horizontally based on the FE experiment and monitored the pressure development, temperature increase, and mechanical deformation at the specific monitoring points. We monitored the capillary pressure as the primary effect inducing the flow in the buffer system, and thermal stress and pressurization were dominant in the surrounding rocks' area. The results will also be compared and validated with the other participating groups and the experimental data further.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.400-414
• /
• 2021

Under high-level radioactive waste repository conditions, bentonite as an engineered barrier material undergoes thermal, hydrological, mechanical, and chemical processes. We report the applications of X-ray Computed Tomography (CT) imaging technique on the characterization and analysis of bentonite over the past decade to provide a reference of the utilization of this technique and the recent research trends. This overview of the X-ray CT technique applications includes the characterization of the bentonite either in pellets or powder form. X-ray imaging has provided a means to extract grain information at the microscale and identify crack networks responsible for the pellets' heterogeneity. Regarding samples of pellets-powder mixtures under hydration, X-ray CT allowed the identification and monitoring of heterogeneous zones throughout the test. Some results showed how zones with pellets only swell faster compared to others composed of pellets and powder. Moreover, the behavior of fissures between grains and bentonite matrix was observed to change under drying and hydrating conditions, tending to close during the former and open during the latter. The development of specializing software has allowed obtaining strain fields from a sequence of images. In more recent works, X-ray CT technique has served to estimate the dry density, water content, and particle displacement at different testing times. Also, when temperature was added to the hydration process of a sample, CT technology offered a way to observe localized and global density changes over time.

• Tunnel and Underground Space
• /
• v.33 no.2
• /
• pp.95-107
• /
• 2023

In this study, the effect of repeated traffic vibration on the long-term stability of mine openings is analyzed for re-utilization of abandoned mine galleries. The research mine in this study is an underground limestone mine which is developed by room-and-pillar mining method, and a dynamic numerical analysis is performed assuming that the research mine will be utilized as a logistics warehouse. The actual traffic vibration generated by the mining vehicles is measured directly, and its waveform is used as input data for dynamic numerical analysis, As a results of dynamic numerical analysis, after 20,000 repetitions of traffic vibration, the mine openings is analyzed to be stable, but an increase in the maximum principal stress and an additional area of plastic zone are observed in the analysis section. As shown in the changes of displacement, volumetric strain, and maximum principal stress which are measured at the mine opening walls. It is confirmed that if the repeated traffic vibration is continuously applied, the instability of the mine openings can be increased. Authors expect that the results of this study can be used as a reference for basic study on utilization of abandoned mine.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.11 no.1
• /
• pp.71-83
• /
• 2009

Shotcrete is an important primary support for tunnelling in rock. The quality control of shotcrete is a core issue in the safe construction and maintenance of tunnels. Although shotcrete may be applied well initially onto excavated rock surfaces, it is affected by blasting, rock deformation and shrinkage and can debond from the excavated surface, causing problems such as corrosion, buckling, fracturing and the creation of internal voids. This study suggests an effective non-destructive evaluation method of the tunnel shotcrete bonding state applied onto hard rocks using the impact-echo (IE) method and ground penetration radar (GPR). To verify previous numerical simulation results, experimental study carried out. Generally, the bonding state of shotcrete can be classified into void, debonded, and fully bonded. In the laboratory, three different bonding conditions were modeled. The signals obtained from the experimental IE tests were analyzed at the time domain, frequency domain, and time-frequency domain (i.e., the Short- Time Fourier transform). For all cases in the analyses, the experimental test results were in good agreement with the previous numerical simulation results, verifying this approach. Both the numerical and experimental results suggest that the bonding state of shotcrete can be evaluated through changes in the resonance frequency and geometric damping ratio in a frequency domain analysis, and through changes in the contour shape and correlation coefficient in a time-frequency analysis: as the bonding state worsens in hard rock condition, the autospectral density increases, the geometric damping ratio decreases, and the contour shape in the time-frequency domain has a long tail parallel to the time axis. The correlation coefficient can be effectively applied for a quantitative evaluation of bonding state of tunnel shotcrete. Finally, the bonding state of shotcrete can be successfully evaluated based on the process suggested in this study.

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.

• Tunnel and Underground Space
• /
• v.28 no.6
• /
• pp.670-691
• /
• 2018

This study presents the research results of the BMT(Benchmark Model Test) simulations of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to predict fault reactivation and the coupled hydro-mechanical behavior of fault. BMT scenario simulations of Task B were conducted to improve each numerical model of participating group by demonstrating the feasibility of reproducing the fault behavior induced by water injection. The BMT simulations consist of seven different conditions depending on injection pressure, fault properties and the hydro-mechanical coupling relations. TOUGH-FLAC simulator was used to reproduce the coupled hydro-mechanical process of fault slip. A coupling module to update the changes in hydrological properties and geometric features of the numerical mesh in the present study. We made modifications to the numerical model developed in Task B Step 1 to consider the changes in compressibility, Permeability and geometric features with hydraulic aperture of fault due to mechanical deformation. The effects of the storativity and transmissivity of the fault on the hydro-mechanical behavior such as the pressure distribution, injection rate, displacement and stress of the fault were examined, and the results of the previous step 1 simulation were updated using the modified numerical model. The simulation results indicate that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing interaction and collaboration with other research teams of DECOVALEX-2019 Task B and validated using the field experiment data in a further study.