• Journal of Korean Tunnelling and Underground Space Association
• /
• v.15 no.3
• /
• pp.345-356
• /
• 2013

Numerical simulation on rock cutting by a TBM disc cutter was carried out using SPH (Smoothed Particle Hydrodynamics) code. AUTODYN3D, a commercial software program based on finite element method, was used in this study. The three-dimensional geometry of a disc cutter and a rock specimen were modeled by Lagrange and SPH code respectively. The numerical simulation was carried out for Hwangdeung granite for 10 different cutting conditions. The results of the numerical simulation, i.e. the relation between cutter force and failure behavior, had a good agreement with those from LCM test. The cutter forces measured in the numerical simulation had 10% deviation from the LCM test results. Moreover, the optimum cutter spacing was almost identical with the experimental results. These results indicate that SPH code can be successfully used had applicability for simulation on rock cutting by a TBM disc cutter. However, further study on Lagrange-SPH coupled modelling would be necessary to reduce the computation time.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.1
• /
• pp.91-104
• /
• 2008

In tunnels, safety factor concept has been suggested to estimate their stability quantitatively. It is merely limited in the framework of mechanical analysis. However safety factor concept has not been applied in hydro-mechanical coupled analyses due to their modelling complexity. Recently studies on this topic are being actively made. In this study, induced drainage modelling methods for hydro-mechanical coupled analyses are compared and analyzed to estimate safety factor of a subsea tunnel exactly. To this end, methods both controlling hydraulic characteristic of shotcrete and using a drainage well are considered. Sensitivity analysis were carried out on rock class, thickness of shotcrete, and hydraulic properties of rock mass. As the results of this study, it turned out that the induced drainage modelling using a drainage well would give more reliable results than that of controlling hydraulic characteristic of shotcrete in estimating tunnel stability in hydro-mechanical coupled analyses.

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

• 한국터널공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.213-224
• /
• 2005

The safety factor of a tunnel considering the failure of supports is important because the failure of supports might cause the collapse of the tunnel. In the previous studies, shotcrete was modelled as beam elements and the failure of the shotcrete was checked according to the allowable working stress concept. In this study, shotcrete was modelled by both beam elements and continuum (elasto-plastic) elements. Safety factors of tunnels were estimated by two dimensional numerical analysis with varying rock mass class, coefficient of lateral pressure, thickness of shotcrete, rock bolt reinforcement and excavation method. Also the study suggested not only a proper amount of supports but also modelling method.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.11a
• /
• pp.485-492
• /
• 2005

Since blast-induced vibration may cause serious problem to the rock mass as well as the nearby structures, the prediction of blast-induced vibration and the stability evaluation must be performed before blasting activities. Dynamic analysis has been Increased recently in order to analyze the effect of the blast-induced vibration. Most of the previous studies, however, were based on the continuum analysis unable to consider rock joints which significantly affect the wave propagation and attenuation characteristics. They also adopted pressure curves estimated by theoretical or empirical equations as input detonation load, thus there were very difficult to reflect the characteristics of propagating media. In this study, therefore, we suggested a discontinuum dynamic analysis technique which uses velocity waveform obtained from a test blast as an input detonation load. A distinct element program, UDEC was used to consider the effect of rock joints. In order to verify the validity of proposed method, the test blast was simulated. The predicted results from the proposed method showed a good agreement with the measured vibration data from the test blast Through the dynamic numerical modelling on the planned road tunnel and slope, we evaluated the effect of blast-induced vibration and the stability of rock slope.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.195-200
• /
• 2000

Numerical analysis using FLACSD has been conducted to estimate the stability of a large underground hall that is to be excavated in a mined area and constructed as an unit of a resort park. Numerical modelling is divided into two stages. The first stage is related to the analysis of the mechanical stability of the hall itself and the second to that of the influence of an adjacent mined cavity upon the hall. In the first stage, the stability of the hall is judged from the interpretation of numerical results in three respects: convergence of the unbalanced force of the model, occurrence of plastic zones and distribution of the displacement. In the second stage, variation of the stress state around the underground hall due to the existence of the cavity is compared to that in the case of the absence of the cavity. Through these analyses, it could be known that the large underground hall is not exposed to any mechanical problems and also not affected by the adjacent cavity. Key words : 3D numerical analysis, large underground cavern, stability analysis

• The Journal of Engineering Geology
• /
• v.1 no.1
• /
• pp.137-145
• /
• 1991

Groundwater flow in fractured rocks generates many challenging problems to scientist and engineers in the projects related to oil and geothermal reservoirs, subsurface contaminations and underground openings. To circumvent these problems, the numerical simulation of groundwater system is used as an established tool in these days. Discrete modelling approach emphasizes geometric parameters, aperture and transport properties of fracture. On the other hand, continuum modelling approach uses the parameters formulated in a way of average hydraulic property. In recent years, the results of field observations from underground opening indicate that groundwater in rock mass flows in a channel form. The channel flow is postulated as the result of the combined effects of geometric pattern and aperture variation.

• Tunnel and Underground Space
• /
• v.8 no.1
• /
• pp.8-16
• /
• 1998

The first purpose of this study is to verify the application of computer modelling to a enclosed space fire. The second one is to determine temperature distribution for the three different ventilation types in case of a enclosed space fire. The third one is to find out the ventilation direction and ventilation quantity to remove effectively heat and combustion products generated by a fire in variable air volume(VAV) system. Firstly, compared with experimental results of Lawrence Livermore National Laboratory(LLNL), numerical results show good agreements. Secondly, among three different ventilation types, the numerical analyses show the highest temperature distribution in occupied zone(up to 1.8 m from bottom) from firing moment to 100 sec. when supply ducts are placed in ceiling and extract duct is placed close to the bottom on side walls. This is due to disadvantageous position of extract duct in ventilating high temperature air which rise because of buoyancy force. Thirdly, this study finds out effective ventilation direction and ventilation quantity to remove heat and combustion products generated by a fire by using VAV system. $CO_2$ concentration is used as a fire fume removal index. As soon as a fire happens, ventilation direction is changed in order to gather and drive out fire fumes. In case of three times ventilation quantity of ordinary one, $CO_2$ concentration and temperature have begun to decrease at 120 sec. after firing, i.e.fire fumes have begun to be removed.

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

• Structural Engineering and Mechanics
• /
• v.69 no.5
• /
• pp.537-545
• /
• 2019

Stress analysis of bottom-hole rock has to be considered with much care to further understand rock fragmentation mechanism and high penetration rate. This original study establishes a fully coupled simulation model and explores the effects of overburden pressure, horizontal in-situ stresses, drilling mud pressure, pore pressure and temperature on the stress distribution in bottom-hole rock. The research finds that in air drilling, as the well depth increases, the more easily the bottom-hole rock is to be broken. Moreover, the mud pressure has a great effect on the bottom-hole rock. The bigger the mud pressure is, the more difficult to break the bottom-hole rock is. Furthermore, the maximum principal stress of the bottom-hole increases as the mud pressure, well depth and temperature difference increase. The bottom-hole rock can be divided into three main regions according to the stress state, namely a) three directions tensile area, b) two directions compression areas and c) three directions compression area, which are classified as a) easy, b) normal and c) hard, respectively, for the corresponding fragmentation degree of difficulty. The main contribution of this paper is that it presents for the first time a thorough study of the effect of related factors, including stress distribution and temperature, on the bottom-hole rock fracture rather than the well wall, using a thermo-poroelastoplasticity model.