• The Journal of Engineering Geology
• /
• v.31 no.1
• /
• pp.55-66
• /
• 2021

Hydraulic-mechanical (H-M) coupled numerical modeling was used to evaluate the evolution of hydrogeological properties in response to the installation and expansion of a borehole. A domain with a discrete fracture network was adopted for discontinuum modeling to simulate changes in fracture apertures. Comparison with real hydraulic test data shows that the effects of principal stress direction and expansion of borehole diameter were reasonably simulated by H-M coupled numerical modeling. The modeling confirmed that aperture changes depended on the principal stress direction, with an increase in aperture size due to vertical displacement being the dominant effect. A concentration of shear dilation around the borehole had an additional, subsidiary, effect on the hydrogeological evolution. These results show that the permeability of fractured rock can be increased by changing the hydraulic properties of a fracture through stress redistribution caused by the installation and expansion of a borehole.

• Tunnel and Underground Space
• /
• v.20 no.1
• /
• pp.28-38
• /
• 2010

In order to investigate the tendency of general displacements and behaviors with respect to each construction process as well as the applicability of numerical analysis schemes, this research has focused on not only analyzing a variety of field observations made in a NATM tunnel, such as displacement of top and side, stress of shotcrete and axial strength of rock bolt, but also carrying out a series of numerical analyses. It was established from the investigation that the 2-dimensional continuum numerical analysis was the one which could more accurately predict displacement of crown and side in the area of one step excavation (patten, P1-P3), while the 2-dimensional discontinuum analysis was the most suitable scheme to study that of two step excavation (patten, P4-P6). In addition, the 2-dimensional continuum analysis enabled to appropriately predict the axial strength of rock bolt and stress of shotcrete in all the area of the tunnel. Finally, it has been possible to conclude from the study that the 3-dimensional continuum analysis should be applied to inspect the behavior and tendency with respect to each stage of the construction as well as in the case of joints, such as large turnouts where relaxation loads in both of horizontal and vertical direction are piled up.

• The Journal of Engineering Geology
• /
• v.20 no.1
• /
• pp.13-23
• /
• 2010

Toppling failure of a slope is defined as failure behavior accompanying the rotation of rock block which is different from other failure such as sliding along with discontinuities and so on. It generally occurs in the region that discontinuities were developed with inverse dip direction to a slope and it could play a critical role in judging stability of slope. In this study, the stability evaluation was performed about toppling failure on a jointed road cut-slope. To check the deformation behavior, numerical analysis is widely used. However common analysis programs are based on continuum model. Recently, many methods that discontinuity properties can be considered in continuum analysis are suggested. In this study, numerical analysis based on FEM(Finite Element Method) was performed using interface element applied in heterogeneous boundary to simulate effects of discontinuities.

• Tunnel and Underground Space
• /
• v.23 no.4
• /
• pp.261-270
• /
• 2013

As part of the development of the 3-D geologic modeling software, this study addresses on new development of software modules that can perform the analysis and visualization of the fracture network system in 3-D. The developed software modules, such as BOUNDARY, DISK3D, FNTWK3D, CSECT and BDM, are coded on Microsoft Visual Studio platform using the MFC and OpenGL library supported by C++ program language. Each module plays a role in construction of analysis domain, visualization of fracture geometry in 3-D, calculation of equivalent pipes, production of cross-section map and management of borehole data, respectively. The developed software modules for analysis and visualization of the 3-D fracture network system can be used to tackle the geomechanical problems related to strength, deformability and hydraulic behaviors of the fractured rock masses. All these benefits will further enhance the economic competitiveness of the domestic software industry.

• Journal of the Korean Geotechnical Society
• /
• v.17 no.4
• /
• pp.241-247
• /
• 2001

절리가 발달한 암반사면의 경우, 토사사면의 안전율 평가를 위해 널리 이용되는 한계평형법을 이용한 해석법은 파괴 활동면을 가장할 수 없기 때문에 사용이 곤란하다. 한편, 비등방 탄소성 모델(편재절리모델)을 사용하여 2조의 절리군을 가진 암반사면의 안전율을 계산할 수 있는 방법이 개발되었지만, 이 방법은 개별절리를 효과적으로 고려하지 못한다. 본 연구에서는 개별절리를 고려한 불연속체 해석에 의한 암반사면의 안정성 해석시 사면의 안전율을 평가하는 기법을 개발하였다. 이를 위해, 절리 전단강도 감소기법을 적용하였으며, UEDC의 내장 언어인 FISH를 사용하여 프로그램을 개발하였다. 또한 실제 절리 암반사면을 대상으로 절리를 측정하고, 개발된 기법에 의한 사면의 안전율을 구하였으며, 동일한 사면에 대해 등가연속체로 가정한 편재절리모델 해석을 수행하여 결과를 비교하였다.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.2
• /
• pp.119-128
• /
• 2008

In the case of a circular shaft, it is expected that asymmetric loads should apply on the surface rather than symmetric loads due to geographical factors and the non-homogeneity of the jointed rock masses. In this study, discontinuous numerical analysis was carried in order to analyze the characteristics of asymmetric load distribution on the wall of the circular shaft due to anisotropy caused by heterogeneity of rock masses affected by the discontinuities like as a Joint. And it was also analyzed that the effect of the mechanical properties varied with the rock mass rating and horizontal stress with depth had influence in the asymmetric load on the wall of the shaft. In the case of considering the effect of the joint as variable, asymmetric load ratio $(R_p)$, which was defined as the ratio of the load subtracted minimum from maximum to minimum, was below 25% in the hard rock. As regarding the variation of the rock mass rating with depth as variable, the value of $R_p$ was below than 25% in the hard rock, and the value between 30% and 40% in the soft rock. On the other hand, the $R_p$ of fractures rock was between $45{\sim}50%$ which value was much higher than that in better rock mass rating.

• Tunnel and Underground Space
• /
• v.11 no.3
• /
• pp.279-288
• /
• 2001

Jointed rock mass can be analyzed by either continuum model or discontinuum model. Finite element method or finite difference method is mainly used for continuum modelling. Although discontinuum model is very attractive in analyzing the behavior of each block in jointed blocky rock masses, it has shortcomings such that it is difficult to investigate each joint exactly with the present technology and the amount of calculation in computer becomes trio excessive. Moreover, in case of the jointed blocky rock mass which has more than 2 dominant joint sets, it is impossible to model the behavior of each block. Therefore, a model such as ubiquitous joint model theory which assumes the rock mass as a continuum, is required. In the case of tunnels, unlike slopes, it is not easy to obtain safety factor by utilizing analysis method based on limit equilibrium method because it is difficult to assume the shape of failure surface in advance. For this reason, numerical analyses for tunnels have been limited to analyzing stability rather than in calculating the safety factor. In this study, the behavior of a tunnel excavated in jointed rock mass is analyzed numerically by using ubiquitous joint model which can incorporate 2 joint sets and a method to calculate safety factor of the tunnel numerically is presented. To this end, stress reduction technique is adopted.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.1
• /
• pp.29-39
• /
• 2012

The maximum unit point resistance ($q_{max}$) of rock socketed drilled shafts subjected to axial loads was investigated by a numerical analysis. A 3D Finite Difference Method (FDM) analysis and a Distinct Element Method (DEM) analysis were performed with varying rock elastic modulus (E), discontinuity spacing ($S_j$), discontinuity dip angle ($i_j$), and pile diameter (D). Based on the results of obtained, it was found that the ultimate point resistance ($q_{max}$) increased as rock elastic modulus (E) and rock discontinuity spacing ($S_j$) increased. But, it was found that $q_{max}$ decreased as pile diameter (D) increased. As for the influence of the dip angle of rock discontinuity ($i_j$), it was shown that $q_{max}$ decreased up to 50% of maximum value within the range of $0^{\circ}$ < $i_j$ < $60^{\circ}$ due to the shear failure at rock discontinuities. Furthermore, it was found that if $20^{\circ}{\leq}i_j{\leq}40^{\circ}$, influence of $i_j$ should be taken into account because $q_{max}$ tended to approach a minimum value as $i_j$ approached a value near the friction angle of the discontinuity (${\phi}_j$).

• The Journal of Engineering Geology
• /
• v.22 no.4
• /
• pp.449-458
• /
• 2012

A sound understanding of the structural characteristics of fractured rock masses is important in designing and maintaining earth structures because their strength, deformability, and hydraulic behavior depend mainly on the characteristics of discontinuity network structures. Despite considerable progress in understanding the structural characteristics of rock masses, the complexity of discontinuity patterns has prevented satisfactory analysis based on a 3-D rock mass visualization model. This paper presents the results of studies performed to develop rock mass visualization in 3-D to analysis the mechanical and hydraulic behavior of fractured rock masses. General and particular solutions of non-linear equations of disk-shaped fractures have been derived to calculated lines of intersection and equivalent pipes. Also, program modules have been developed to perform the calculations. The procedures developed for the 3-D fractured rock mass visualization model can be used to characterize rock mass geometry and network systems effectively. The results obtained in this study will be refined and then combined for use as a tool for assessing geomechanical problems related to strength, deformability and hydraulic behaviors of the fractured rock masses.

• Explosives and Blasting
• /
• v.37 no.2
• /
• pp.14-21
• /
• 2019

Distinct element method is a powerful numerical tool for modelling the jointed rock masses. It is also a useful tool for modelling of later stage of blasting requiring large displacement. The distinct element method utilizes a rigid block idea in which the interacting force between distinct elements is calculated from contact displacement as elements penetrate slightly. The properties of joints defined as the boundaries of distinct elements are critical parameters to determine the block behavior, and affect the deformation and failure mode. However, regardless of real joint properties, joint stiffnesses have sometimes been selected without special concern just to prevent elements from penetrating too far into each other in some quasi-static problems. Depending on whether the main interest in the analysis is the prediction of the deformation with high precision, or the prediction of the block behaviour after failure, the input data such as joint stiffness may or may not have a significant effect on the results. The purpose of this study is to provide a sound understanding of the effect of the joint stiffness on the distinct element analysis results, and to help guide the selection of input data.