• Journal of the Korean Geotechnical Society
• /
• v.19 no.5
• /
• pp.211-221
• /
• 2003

At Aspo hard rock laboratory in Sweden, an in-situ heater experiment called "$\"{A}"{s}"{p}"{o}$ Pillar Stability Experiment (APSE)" is prepared to assess capability to predict spatting and stability in a rock mass between deposition holes for radioactive waste. To Predict reasonably fracturing process at rock pillar under a planned configuration before testing, a boundary element code FRACOD has been applied for modelling. The code has been improved to simulate explicitly fracture evolution both at rock boundaries and in intact rocks. A new inverse stress reconstruction technique using boundary element has been also developed to transfer stress field by excavation and thermal loading into the FRACOD model. This article presents the results from predictive modelling far the planned in-situ test condition. Excavation induced stresses might cause slight fracturing in the pillar walls. Typical shear fractures have been initiated and propagated near central pillar walls during 120 days of heating, but overall rock mass remained stable under the considered configuration. The effects of pre-existing joints and properties of fractures are also discussed. It is found from the results that FRACOD can properly model essential rock spatting and propagation at deep tunnels and boreholes.at deep tunnels and boreholes.

• 한국지구물리탐사학회:학술대회논문집
• /
• 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.

• Tunnel and Underground Space
• /
• v.21 no.1
• /
• pp.66-81
• /
• 2011

In this study, it was aimed to develop a thermal-hydraulic-mechanical coupled fracture mechanics code that models a fracture initiation, propagation and failure of underground rock mass due to thermal and hydraulic loadings. The development was based on a 2D FRACOD (Shen & Stephasson, 1993), and newly developed T-M and H-M coupled analysis modules were implemented into it. T-M coupling in FRACOD employed a fictitious heat source and time-marching method, and explicit iteration method was used in H-M coupling. The validity of developed coupled modules was verified by the comparison with the analytical result, and its applicability to the fracture initiation and propagation behavior due to temperature changes and hydraulic fracturing was confirmed by test simulations.

• Geomechanics and Engineering
• /
• v.22 no.6
• /
• pp.529-540
• /
• 2020

Water inrush is a major hazard for mining and excavation in deep coal seams or rock masses. It can be attributed to the coalescence of rock fractures in rock mass due to the interaction of fractures, hydraulic flow and stress field. One of the key technical challenges is to understand the course and mechanism of fluid flows in rock joint networks and fracture propagation and hence to take measures to prevent the formation of water inrush channels caused by possible rock fracturing. Several case observations of fluid flowing in rock joint networks and coupled fracture propagation in underground coal roadways are shown in this paper. A number of numerical simulations were done using the recently developed flow coupling function in FRACOD which simulates explicitly the fracture initiation and propagation process. The study has demonstrated that the shortest path between the inlet and outlet in joint networks will become a larger fluid flow channel and those fractures nearest to the water source and the working faces become the main channel of water inrush. The fractures deeper into the rib are mostly caused by shearing, and slipping fractures coalesce with the joint, which connects the water source and eventually forming a water inrush channel.

• Tunnel and Underground Space
• /
• v.14 no.4
• /
• pp.248-260
• /
• 2004

Damage in brittle rock due to stress increase starts from initiation of microcracks, and then results in failure by forming macro failure planes due to propagation and coalescence of these discrete cracks. Conventionally, continuum approaches using macro-failure criteria or a number of elasto-plastic models have been major solution to implement rock damage and failure. However, actual brittle failure processes can be better described in phenomenological approach if initiation and propagation of discrete fractures are explicitly considered. This study presents damage and failure process of rock using a boundary element code, FRACOD, which has been developed to model fracturing process of rocks. Through a series of numerical uniaxial compressive tests, the feasibility of the developed model was verified, and realistic rock failure process was reproduced considering scale effects in rocks. In addition, the fracturing process and the corresponding rock damage in the vicinity of deep shaft in rock mass were presented as an application of this approach. This approach will be expected to contribute to finding better engineering solutions for the analysis of stability problems in brittle rock masses.

• Geomechanics and Engineering
• /
• v.16 no.1
• /
• pp.35-47
• /
• 2018

This paper investigates the mechanisms of tunnel spalling and massive tunnel failures using fracture mechanics principles. The study starts with examining the fracture propagation due to tensile and shear failure mechanisms. It was found that, fundamentally, in rock masses with high compressive stresses, tensile fracture propagation is often a stable process which leads to a gradual failure. Shear fracture propagation tends to be an unstable process. Several real case observations of spalling failures and massive shear failures in boreholes, tunnels and underground roadways are shown in the paper. A number of numerical models were used to investigate the fracture mechanisms and extents in the roof/wall of a deep tunnel and in an underground coal mine roadway. The modelling was done using a unique fracture mechanics code FRACOD which simulates explicitly the fracture initiation and propagation process. The study has demonstrated that both tensile and shear fracturing may occur in the vicinity of an underground opening. Shallow spalling in the tunnel wall is believed to be caused by tensile fracturing from extensional strain although no tensile stress exists there. Massive large scale failure however is most likely to be caused by shear fracturing under high compressive stresses. The observation that tunnel spalling often starts when the hoop stress reaches $0.4^*UCS$ has been explained in this paper by using the extension strain criterion. At this uniaxial compressive stress level, the lateral extensional strain is equivalent to the critical strain under uniaxial tension. Scale effect on UCS commonly believed by many is unlikely the dominant factor in this phenomenon.

• Tunnel and Underground Space
• /
• v.24 no.1
• /
• pp.67-80
• /
• 2014

Long-term stability and delayed failure of granite were evaluated through the laboratory test based on Wilkins method and Brazilian disc test (BDT) which yields tensile strength, mode I fracture toughness and subcritical crack growth parameters. Then, the long-term strength of granite was estimated by using analytical models and long-term stability of compressed air-energy storage (CAES) pilot cavern pressurized up to 5 ~ 6 MPa was evaluated using numerical code, FRACOD with the determined subcritical crack growth parameters. The results of test and analyses showed that the subcritical crack growth index, n was determined as 29.39 and the inner pressure of 5 ~ 6 MPa had an insignificant effect on the long-term stability of pilot cavern. It was also found that the measurement and analysis of acoustic emission events can describe the accumulation of damage due to subcritical crack growth quantitatively. That is, AE monitoring can provide the current status of rock under loading if we make an identical installation condition in the field with that of the laboratory test.

• Tunnel and Underground Space
• /
• v.30 no.2
• /
• pp.109-135
• /
• 2020

Globally, the deepening depth in the underground is a situation of the high interest for a purpose of the development of various facilities. The development of deep underground space should be based on the structural stability of rocks. Spalling is known to have an impact on the structural stability degradation in deep underground space. As an attempt to predict spalling, many researchers have proposed predicted conditions in accordance with stress states which occur around the tunnel, rock conditions, and types of rock. In addition, the analysis on spalling method has been verified by using computer modeling such as FLAC, EXAMINE, Insight 2D, UDEC and FRACOD, along with in-situ measurement results. In Canada URL (Underground Research Tunnel), CWFS model (Cohesion Weakening Frictional Strengthening) was used to precisely predict for the state of spalling, comparing spalling modeling. CWFS model has been identified as a reliable method for predicting such phenomena. This study aims to analyze several cases of spalling, and then make a comparison between the conditions for spalling occurrence and the predicted results of model CWFS. With this, it investigates the applicability of prediction of spalling, targeting pillar under deep depth condition.

• Tunnel and Underground Space
• /
• v.31 no.5
• /
• pp.309-332
• /
• 2021

Natural barrier systems surrounding the geological repository for the high-level radioactive waste should guarantee the hydraulic performance for preventing or delaying the leakage of radionuclide. In the case of the behavior of a crystalline rock, the hydraulic performance tends to be decided by the existence of discontinuities, so the coupled hydro-mechanical(HM) processes on the discontinuities should be characterized. The discontinuum modelling can describe the complicated behavior of discontinuities including creation, propagation, deformation and slip, so it is appropriate to model the behavior of a crystalline rock. This paper investigated the coupled HM processes in discontinuum modelling such as UDEC, 3DEC, PFC, DDA, FRACOD and TOUGH-UDEC. Block-based discontinuum methods tend to describe the HM processes based on the fluid flow through the discontinuities, and some methods are combined with another numerical tool specialized in hydraulic analysis. Particle-based discontinuum modelling describes the overall HM processes based on the fluid flow among the particles. The discontinuum methods that are currently available have limitations: exclusive simulations for two-dimension, low hydraulic simulation efficiency, fracture-dominated fluid flow and simplified hydraulic analysis, so it could be improper to the modelling the geological repository. Based on the concepts of various discontinuum modelling compiled in this paper, the advanced numerical tools for describing the accurate coupled HM processes of the deep geological repository should be developed.