• Tunnel and Underground Space
• /
• v.15 no.4 s.57
• /
• pp.264-274
• /
• 2005

The optimum blasting pattern to excavate a quarry efficiently and economically can be determined based on the minimum production cost which is generally estimated according to rock fragmentation. Therefore it is a critical problem to predict fragment size distribution of blasted rocks over an entire quarry. By comparing various prediction models, it can be ascertained that the result obtained from Kuz-Ram model relatively coincides with that of field measurements. Kuz-Ram model uses the concept of rock factor to signify conditions of rock mass such as block size, rock jointing, strength and others. For the evaluation of total production cost, it is imperative to estimate 3-D spatial distribution of rock factor for the entire quarry. In this study, a sequential indicator simulation technique is adopted for estimation of spatial distribution of rock factor due to its higher reproducibility of spatial variability and distribution models than Kriging methods. Further, this can reduce the uncertainty of predictor using distribution information of sample data The entire quarry is classified into three types of rock mass and optimum blasting pattern is proposed for each type based on 3-D spatial distribution of rock factor. In addition, plane maps of rock factor distribution for each ground levels is provided to estimate production costs for each process and to make a plan for an optimum blasting pattern.

• Wind and Structures
• /
• v.27 no.4
• /
• pp.213-221
• /
• 2018

In the desert and Gobi regions with strong wind and large sediment discharge, sand transporting engineering is more effective than sand blocking and sand fixing measures in sand prevention. This study uses the discrete phase model of 3D numerical simulation to study the motion trail, motion state and distribution rule of sand particles with different grain diameters when the included angle between the main shaft of the feather-row lateral transportation sand barrier and the wind direction changes, and conducts a comparison in combination with the wind tunnel test and the flow field rule of common sand barrier. According to the comparison, when wind-sand incoming flow passes through a feather-row sand barrier, sand particles slow down and deposit within the deceleration area under the resistance of the feather-row sand barrier, move along the transportation area formed by the transportation force, and accumulate as a ridge at the tail of the engineering. With increasing wind speed, the eolian erosion of the sand particles to the ground and the feather-row sand barrier is enhanced, and the sand transporting quantity and throw-over quantity of the feather-row sand barrier are both increased. When sand particles with different grain diameters bypass the feather-row sand barrier, the particle size of the infiltrating sands will increase with the included angle between the main shaft of the feather-row sand barrier and the wind direction. The obtained result demonstrates that, at a constant wind speed, the flow field formed is most suitable for the lateral transportation of the wind-drift flow when the included angle between the main shaft of the feather-row sand barrier lateral transportation engineering and the wind speed is less than or equal to $30^{\circ}$.

• Tunnel and Underground Space
• /
• v.33 no.1
• /
• pp.57-69
• /
• 2023

In this study, we simulated induced seismicity in the Groningen natural gas reservoir using 2D hydro-mechanical coupled discrete element modelling (DEM). The code used is PFC2D (Particle Flow Code 2D), a commercial software developed by Itasca, and in order to apply to this study we further developed 1)initialization of inhomogeneous reservoir pressure distribution, 2)a non-linear pressure-time history boundary condition, 3)local stress field monitoring logic. We generated a 2D reservoir model with a size of 40 × 50 km² and a complex fault system, and simulated years of pressure depletion with a time range between 1960 and 2020. We simulated fault system failure induced by pressure depletion and reproduced the spatiotemporal distribution of induced seismicity and assessed its failure mechanism. Also, we estimated the ground subsidence distribution and confirmed its similarity to the field measurements in the Groningen region. Through this study, we confirm the feasibility of the presented 2D hydro-mechanical coupled DEM in simulating the deformation of a complex fault system by hydro-mechanical coupled processes.

• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.509-518
• /
• 2020

This study constructed a new simulation platform, including mesh generation process, numerical simulation, and post-processing for results analysis based on exploration data to perform real-scale numerical analysis considering the actual geological structure efficiently. To build the simulation platform, we applied for open-source programs. The source code is open to be available for code modification according to the researcher's needs and compatibility with various numerical simulation programs. First, a three-dimensional model(3D) is acquired based on the exploration data obtained using a drone. Then, the domain's mesh density was adjusted to an interpretable level using Blender, the free and open-source 3D creation suite. The next step is to create a 3D numerical model by creating a tetrahedral volume mesh inside the domain using Gmsh, a finite element mesh generation program. To use the mesh information obtained through Gmsh in a numerical simulation program, a converting process to conform to the program's mesh creation protocol is required. We applied a Python code for the procedure. After we completed the stability analysis, we have created various visualization of the study using ParaView, another open-source visualization and data analysis program. We successfully performed a preliminary stability analysis on the full-scale Dokdo model based on drone-acquired data to confirm the usefulness of the proposed platform. The proposed simulation platform in this study can be of various analysis processes in future research.

• Tunnel and Underground Space
• /
• v.26 no.6
• /
• pp.516-523
• /
• 2016

Roughness of a joint surface is one of the most important parameters that affects the mechanical and hydraulic behavior of rock mass. Therefore, various studies on making constitutive model and/or roughness quantification have been conducted in experimental and empirical manners. Advances in recent 3D printing technology can be utilized to generate a joint surface with a specific roughness. In this study, a reliable technique to generate a rough joint surface was introduced and its quantitative assessment was made. Random midpoint displacement method was applied to generate a joint surface and the distribution of $Z_2$ was investigated to assess its roughness. As a result, a certain roughness can be embodied by controlling input parameters and furthermore it was able to generate a joint surface with specific roughness anisotropy.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.320-334
• /
• 2020

This study presents an application of hydro-mechanical coupled Particle Flow Code 3D (PFC3D) to simulation of fluid injection induced fault slip experiment conducted in Mont Terri Switzerland as a part of a task in an international research project DECOVALEX-2019. We also aimed as identifying the current limitations of the modelling method and issues for further development. A fluid flow algorithm was developed and implemented in a 3D pore-pipe network model in a 3D bonded particle assembly using PFC3D v5, and was applied to Mont Terri Step 2 minor fault activation experiment. The simulated results showed that the injected fluid migrates through the permeable fault zone and induces fault deformation, demonstrating a full hydro-mechanical coupled behavior. The simulated results were, however, partially matching with the field measurement. The simulated pressure build-up at the monitoring location showed linear and progressive increase, whereas the field measurement showed an abrupt increase associated with the fault slip We conclude that such difference between the modelling and the field test is due to the structure of the fault in the model which was represented as a combination of damage zone and core fractures. The modelled fault is likely larger in size than the real fault in Mont Terri site. Therefore, the modelled fault allows several path ways of fluid flow from the injection location to the pressure monitoring location, leading to smooth pressure build-up at the monitoring location while the injection pressure increases, and an early start of pressure decay even before the injection pressure reaches the maximum. We also conclude that the clay filling in the real fault could have acted as a fluid barrier which may have resulted in formation of fluid over-pressurization locally in the fault. Unlike the pressure result, the simulated fault deformations were matching with the field measurements. A better way of modelling a heterogeneous clay-filled fault structure with a narrow zone should be studied further to improve the applicability of the modelling method to fluid injection induced fault activation.

• Tunnel and Underground Space
• /
• v.19 no.2
• /
• pp.86-94
• /
• 2009

A three dimensional rock joint element was developed considering fracture mechanics and subcritical crack growth to simulate non-linear behavior and the progressive failure of rock joints. Using this 3-D joint element, joint shear tests of rock discontinuities were simulated by a numerical method. The asperities on the joint surface began to fail at stress levels lower than the rock fracture toughness and continued progressively due to subcritical crack growth. As a result of progressive failing in each and every asperity, the joint showed non-linear stress-time behavior including stress hardening/softening and the reaching of a residual stress.

• Journal of the Korean Society of Visualization
• /
• v.7 no.2
• /
• pp.64-71
• /
• 2010

The flow and noise characteristics of wake behind wind-turbine blades have been investigated experimentally using a two-frame particle image velocimetry (PIV) technique. Experiments were carried out in a POSTECH subsonic large wind-tunnel ($1.8^W{\times}1.5^H{\times}4.3^L\;m^3$) with KBP-750D (3-blade type) wind-turbine model at a freestream velocity of $U_o\;=\;15\;m/s$ and a tip speed ratio $\lambda\;=\;6.14$ (2933 rpm). The wind-turbine blades are connected to an AC servo motor, brake, encoder and torque meter to control the rotational speed and to extract a synchronization signal for PIV measurements. The wake flow was measured at four azimuth angles ($\phi\;=\;0^{\circ}$, $30^{\circ}$, $60^{\circ}$ and $90^{\circ}$) of the wind-turbine blade. The dominant flow structure of the wake is large-scale tip vortices. The turbulent statistics such as turbulent intensity are weakened as the flow goes downstream due to turbulent dissipation. The dominant peak frequency of the noise signal is identical to the rotation frequency of blades. The noise seems to be mainly induced by the tip vortices.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.323-323
• /
• 2017

본 연구 대상인 주암댐-주암조절지댐간 도수터널은 저수지 간 연계운영을 통하여 서 남해안 일대의 용수 공급하기 위하여 건설되어 운영되었으나, 터널의 구조적 불안정성이 제기됨에 따라 향후 발생 가능성이 있는 용수공급 중단을 방지하기 위하여 신규 도수터널이 제안되었다. 계획된 신규 도수 터널의 주요 시설물은 크게 양방향 운영이 가능한 터널(D=3.3 m, L=11.23 km), 각 저수지 상황별 운영을 위한 취수문비 2개소, 도수가 이루어지는 상황에서 수문 돌발 폐쇄시 수충격을 감쇄하기 위한 배기구(air vents) 2개소가 계획되었다. 이에 따라 본 연구에서는 주암댐에서 주암조절지댐으로 최대 유량이 통수되는 상태에서, 수문폐쇄에 따른 수충격을 정량적으로 분석하고자 Joukowsky 공식에 적용하여 완폐쇄와 급폐쇄시 도수터널의 안정성을 검토하였으며, 수문 폐쇄로 인한 천이적인 흐름상태 등 수충격 모의가 가능한 1차원 ITM 모형을 적용하여, 수리분석과 수치모형과 결과 비교하고, 계획된 배기구의 유무에 따른 효과를 알아보고자 하였다. 분석 결과, 계획된 0.3 m/min으로 수문을 폐쇄할 경우, 도수터널의 안정성에는 문제가 없을 것으로 분석되었으나, 수문을 급 폐쇄 할 경우, 압력수두가 크게 증가하여 도수터널에 위험이 있을 것으로 분석되었으며 배기구 유무에 따른 수충격 검토 결과 도수터널 내의 수압상승을 적절히 조절하는 조압수조의 역할을 만족스럽게 수행할 것으로 분석되었다.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.22 no.3
• /
• pp.277-291
• /
• 2020

Shield TBM is a tunnelling method that has a wider range of applications in the poor ground condition compared to conventional tunnels (Drill and Blast). Currently, a 13.3 m large-diameter slurry shield TBM is preparing for construction to pass under the Han River. Shield TBM is divided into slurry and EPB shield TBM, and management items during construction are different depending on each characteristic. In this paper, the equipment type, origin, application case and trouble case were analyzed for slurry shield TBM, which is mainly constructed in soft ground. In addition, 2D and 3D model tests were conducted on the condition of soil depth for the possibility of slurry leakage into front of the equipment, with appropriate chamber pressure. Based on this paper, it proposed to provide basic and reference data for proper excavation surface pressure and chamber pressure during construction of slurry shield TBM under soft ground conditions, and proposed measures to minimize stability and environmental decline due to slurry ejection.